- Research article
- Open Access
Macrophage mediated recognition and clearance of Borrelia burgdorferi elicits MyD88-dependent and -independent phagosomal signals that contribute to phagocytosis and inflammation
BMC Immunology volume 22, Article number: 32 (2021)
Macrophages play prominent roles in bacteria recognition and clearance, including Borrelia burgdorferi (Bb), the Lyme disease spirochete. To elucidate mechanisms by which MyD88/TLR signaling enhances clearance of Bb by macrophages, we studied wildtype (WT) and MyD88−/− Bb-stimulated bone marrow-derived macrophages (BMDMs).
MyD88−/− BMDMs exhibit impaired uptake of spirochetes but comparable maturation of phagosomes following internalization of spirochetes. RNA-sequencing of infected WT and MyD88−/− BMDMs identified a large cohort of differentially expressed MyD88-dependent genes associated with re-organization of actin and cytoskeleton during phagocytosis along with several MyD88-independent chemokines involved in inflammatory cell recruitment. We computationally generated networks which identified several MyD88-dependent intermediate proteins (Rhoq and Cyfip1) that are known to mediate inflammation and phagocytosis respectively.
Our findings show that MyD88 signaling enhances, but is not required, for bacterial uptake or phagosomal maturation and provide mechanistic insights into how MyD88-mediated phagosomal signaling enhances Bb uptake and clearance.
Lyme disease (LD) is a highly prevalent tick-borne illness caused by the spirochetal bacterium Borrelia burgdorferi (Bb) [1,2,3]. The disease is characterized by a wide array of clinical manifestations which vary in duration and severity among patients. Early clinical manifestations of LD include the characteristic “bullseye” rash known as erythema migrans and flu-like symptoms, while late manifestations include arthritis, carditis and neurological compromise [4, 5]. The invading spirochete induces both innate and adaptive immune responses, and it is believed that the innate immune response to Bb contributes to the development of clinical findings in LD . The macrophage is a principal cellular element of the innate immune response to the bacterium at sites of infection in both humans and mice [7,8,9]. Macrophages also play a prominent role in the pathogenesis of murine Lyme carditis, and their recruitment to heart tissue is important in spirochetal clearance . Macrophages have the phagocytic and signaling machinery necessary to bind, engulf, and degrade Bb. Binding of Bb to macrophages is mediated by surface integrins, such as Complement Receptor 3 (CR3) [11, 12] and α3 . Once attached, phagocytosis of Bb is complex and can occur by either a sinking or coiling mechanism [14, 15]. Both cases require rearrangements of the actin cytoskeleton to internalize Bb into the endosome, where degradation takes place .
Bb is an extracellular pathogen that needs to be taken up and degraded for significant recognition by the host immune system . We have defined this process as “phagosomal signaling” . Spirochete degradation exposes borrelial pathogen-associated molecular patterns (PAMPs), such as lipoproteins and nucleic acids, to endosomal toll-like receptors (TLRs) for recognition, resulting in signaling cascades which induce pro-inflammatory cytokine production [17,18,19]. Bb does not contain LPS and, therefore, does not engage TLR4. The cell envelope of Bb contains abundant triacylated lipoproteins , which are known to be recognized by TLR1/2 heterodimers [21,22,23,24,25,26]. However, the three fatty acid chains in the N-terminus of Bb lipoproteins, which serve as the TLR2/1 PAMP, are tethered in the outer membrane . We have shown that this results in minimal recognition of lipoproteins in intact spirochetes at the cell surface [17, 18, 21, 28]. Instead, principal recognition of Bb TLR2 ligands occurs within macrophage endosomal structures after the spirochete is phagocytosed and degraded [17, 28]. Bacterial degradation results in exposure of both lipoprotein ligands and nucleic acids, which are recognized by endosomal TLR2 and TLRs 7, 8 and 9 respectively [18, 19, 29]. Signaling cascades initiated by engagement of these TLRs utilize the adaptor protein MyD88 [14, 30], indicating that this adaptor protein is a crucial element in mediating the inflammatory response to Bb.
A role for MyD88 has been implicated in each of the four general steps associated with phagocytic clearance of bacterial pathogens: uptake, phagosome maturation, degradation and cytokine production. Murine macrophages lacking MyD88 show markedly diminished uptake of several bacterial species, including Bb [28, 31,32,33,34,35]. In WT macrophages, prior studies have shown that Bb-induced MyD88 signaling results in increased PI3K activation and when PI3K is inhibited Bb uptake is decreased . In addition, formin proteins (FMNL1, mDia1, and Daam1) have been shown to play a critical role in mediating phagocytosis of Bb [15, 37]. Whether MyD88 increases activation of these formins, and the role of PI3K signaling in this process, has not been established. Degradation of bacteria is impaired in the absence of MyD88 due to inefficient acidification of phagosomes . In the context of Bb infection, lysosome maturation markers are recruited to Bb-containing phagosomes in macrophages lacking MyD88 . However, the degree of phagosome maturation and acidification required to expose Bb ligands from the bacteria cell envelope for recognition has not been studied. Murine macrophages lacking MyD88 also show markedly diminished production of NFκB-triggered pro-inflammatory cytokines, such as TNFα and IL-6, when stimulated with different bacterial species, including Bb [28, 31]. Nevertheless, the key host components involved downstream of these MyD88-mediated phagosome signals and their effects have not been well studied in the context of Bb infection.
The objective of these studies is to examine which downstream effects of MyD88 phagosomal signaling potentially enhance clearance of Bb. Using an ex vivo murine macrophage system, we show that MyD88 signaling enhances, but is not required, for bacterial uptake or phagosomal maturation. Through RNA-sequencing analysis, we provide evidence that MyD88 signaling drives transcription of multiple genes involved in phagocytosis and identify potential intermediate proteins that facilitate the association between MyD88 and bacterial uptake. We also demonstrate that internalization of Bb by macrophages induces robust MyD88-independent inflammatory responses via production of chemokines. Our findings highlight the importance of MyD88 in efficient uptake of the Lyme disease spirochete by macrophages and provide potential mechanistic insight into how MyD88 mediates this process.
Female 6–8-week-old C57BL/6 J wild type (WT) and C57BL/6 J MyD88−/− (MyD88−/−) mice used in these studies were obtained from breeding colonies maintained in the UConn Health (UCH) Center for Comparative Medicine facility according to guidelines set by the UCH Institutional Animal Care and Use Committee (IACUC). Female mice were used exclusively in these studies due to differences in expression of TLR7 between males and females . Original WT breeding pairs were purchased from The Jackson Laboratory. Original MyD88−/− breeding pairs were kindly provided by Dr. Egil Lien at the University of Massachusetts with permission from Dr. S. Akira in Osaka, Japan. Disruption of the murine MyD88 gene was confirmed through PCR . Both WT and MyD88−/− breeding colonies are maintained on the antibiotic Sulfatrim (sulfomethoxazole [40 mg/mL] + trimethoprim [8 mg/mL]) diluted in water 1:50, which has been previously shown to not impact the degree of Bb infection . In preparation for euthanasia, individual mice were exposed to isoflurane to calm the animal prior to sedation by injection with an overdose of an anesthetic cocktail [Ketamine 50–75 mg/kg at 15 mg/ml, Xylazine 5–7.5 mg/kg at 2.5 mg/ml, and Acepromazine 0.5–1.25 mg/kg]. Following euthanasia by cervical dislocation, each animal was confirmed to lack a heartbeat in accordance with the approved IACUC protocol.
Low-passage virulent wild-type Bb strain 297  or a strain 297 isolate containing a stably-inserted copy of green fluorescent protein (GFP) under the control of the constitutively-expressed flaB promoter (Bb914)  were maintained in Barbour-Stonner-Kelly (BSK)-II media supplemented with normal rabbit serum and gentamicin (50 μg/μl) . Cultures were grown at 23 °C for at least 1 week prior to being shifted to 37 °C as previously described . Spirochetes were centrifuged at 3300 x g for 20 min at 4 °C and resuspended in either BSK-II for in vivo experiments or DMEM (Gibco, 15,630–080) supplemented with sodium pyruvate (Gibco, 11,360–070) and HEPES (Gibco, 15,630–080) for ex vivo experiments. After resuspension, the cultures were counted by dark-field microscopy using a Petroff-Hausser counting chamber (Hausser Scientific) and diluted accordingly. Staphylococcus aureus (Sa) was cultured and fluorescently labeled with FITC as previously described .
Bone marrow-derived macrophages (BMDMs) were isolated from 6 to 8-week-old WT and MyD88−/− mice as described previously . Single cell macrophage suspensions were seeded into either 12-well tissue culture-treated plates at a concentration of 1 × 106 cells/ml per well or 1 × 105 cells/500 μL per well in 8-chamber cell microscopy slides. Cells were then incubated overnight at 37 °C/5% CO2 to allow cell adherence before experimentation. Cells were incubated for either 0.5, 1, 4 or 6 h at 37 °C/5% CO2 with live GFP-Bb or labeled Sa at multiplicities of infection (MOIs) of either 10 or 100. Stimulation media was DMEM supplemented with 1% sodium pyruvate and 1% HEPES. At the end of the incubation period, culture supernatants were collected and stored at − 80 °C until cytokine analysis. Cells stimulated in chamber slides were processed for confocal microscopy. Cells stimulated in 12-well plates were processed for RNA extraction. All culture media and reagents were confirmed free of LPS contamination (< 10 pg/ml) by Limulus amoebocyte lysate assay quantification (Cambrex, MA).
After stimulation, BMDMs were fixed in 2% paraformaldehyde with 0.05% Triton-X-100 (Fisher, BP151–100) for 10 min. Slide wells were then incubated with 5% bovine serum albumin (BSA) solution in PBS overnight at 4 °C to block non-specific antibody binding. The next day, cells were stained with different combinations of anti-GFP (Thermo Scientific A-21311, 1:100), phalloidin conjugated with Alexa Fluor 647 (Biolegend 424205, 1:20), anti-MyD88 (Santa Cruz 11356, 1:100), anti-TLR2 (eBioscience 14–9021-82, 1:100), anti-TLR7 (R&D MAB7156, 1:100), anti-ASC (Santa Cruz 22514-R, 1:100) and anti-LAMP-1 (eBioscience 14–1071-82, 1:100). A secondary antibody, Alexa Fluor 350, was used to detect anti-MyD88, anti-TLR2, anti-ASC and anti-LAMP-1 (Life Technologies A21093, 1:100). Incubations with primary and secondary antibodies were done for 1 h each at room temperature; slide wells were washed following each incubation three times with PBS supplemented with 0.5% Tween-20, with a final wash in distilled H2O before mounting. After antibody staining, slides were mounted using Vectashield (Vector H-1000) and imaged using a Zeiss 880 confocal microscope. Image processing and analysis were performed using ImageJ (NIH, v1.41b). Colocalization values were determined by first analyzing profile plots in ImageJ (Plug-in: “Plot Profile”) across ten different phagosomes for each cell genotype and then calculating the average difference between the fluorescence intensity curves of the markers of interest (i.e., LAMP-1 and Bb). Binding percentages were calculated by imaging 100–200 cells using a confocal microscope and then measuring the ratio of cells containing at least one surface-bound or internalized spirochete to the total number of cells imaged for each condition, represented as %BMDMs interacting w/Bb. Uptake percentages were calculated by imaging 100–200 cells using a confocal microscope and then measuring the ratio of cells containing at least one internalized spirochete to the total number of cells imaged for each condition, represented as %BMDMs w/internalized Bb.
Western blotting of BMDM supernatants and lysates
Protein lysates were generated from BMDM cell culture lysates and supernatants after Bb stimulation. In these experiments, adenosine triphosphate (ATP) (Sigma, 3A6419-1G) was added to WT BMDMs (already stimulated with Bb for 5 h) 1 h prior to harvest for generation of lysates. Supernatants were treated with an equal volume of methanol and ¼ volume of chloroform, vortexed and spun at 16000 x g for 10 min. After removal of the upper phase, 500 μL of methanol was added to the intermediate phase, which was then vortexed and spun at 16000 x g for 10 min. The pellets were then dried at room temperature, resuspended in 30 μL of 2x Laemmli buffer and incubated in a 37 °C water bath until proteins became soluble. BMDMs were lysed using RIPA buffer at − 80 °C and spun at maximum speed for 10 min. Protein pellets were resuspended in 2x Laemmli buffer. Lysates were boiled at 99 °C for 10 min and run on a 12.5% SDS-PAGE gel at 140 V for 1 h (5 μL per lane, 15 lanes). Proteins were then transferred to nitrocellulose membranes (Bio-Rad 162–0177) at 20 V for 20 min. Membranes were blocked for 1 h in milk block solution and then incubated overnight at 4 °C with primary antibodies for either β-actin (Sigma A5441, 1:2000), IL-1β (R&D AF401NA, 1:800) or caspase-1 (Adipogen AG-20B-0042, 1:1000) diluted in milk block solution. Membranes were then washed 5 times for 5 min each in wash buffer (PBS supplemented with 0.5% Tween-20) and incubated with goat anti-mouse HRP-conjugated IgG (GE NA931) diluted 1:5000 (β-actin and Caspase-1) or 1:1000 (IL-1β) in milk block for 2 h at room temperature. Following additional washes, membranes were incubated in HyGlo spray chemilunescent substrates (Denville Scientific, E2400) for 5 min and imaged on a Biorad ChemiDoc MP imaging system.
The Cytokine Bead Array Mouse Inflammation kit (BD Biosciences 552364) was used according to manufacturer’s instructions for simultaneous measurement of IL-6, IL-10, CCL2, IFNγ, TNFα, and IL-12p70 in supernatants from stimulated BMDMs. General statistical analysis was performed using GraphPad Prism 4.0 (GraphPad Software, San Diego, CA), using an unpaired Student t test. For each experiment, both the standard deviation and the standard error of the mean were calculated. P-values of < 0.05 were considered significant.
Identification of differentially expressed genes by RNA-Seq
Total RNA was extracted from three biological replicates of WT and MyD88−/− BMDMs, either unstimulated or stimulated with Bb at MOI 10:1 or MOI 100:1 for 6 h. Following stimulation, RNA was isolated using the Macherey-Nagel total RNA isolation kit (Takara, 740955) and was used as input for the NuGen Ovation RNA-seq V1 kit. cDNA output was analyzed for correct size distribution with an Experion Standard Sensitivity RNA chip and quantified using a Qubit Fluorometer. Sequencing libraries were produced using the NuGen Encore NGS Library I kit. Libraries were multiplexed and sequenced at The Jackson Laboratory for Genomic Medicine Sequencing Core with an Illumina HiSeq 2500 as 2X50bp pair end reads. RNA-Seq reads from each individual library were mapped with Tophat2 RNA-Seq spliced reads mapper (version 2.0.5)  to mouse genome build mm9 with parameter settings adjusted to suit strand-specific pair-end RNA-Seq reads. The mapping result bam files were used as input to the HTSeq high-throughput sequencing data analysis package  to quantify the read counts mapped to all genes in UCSC mm9 mouse gene annotation set. The expression levels of genes represented as mapped read counts were normalized using the DESeq2 RNA-Seq analysis package (function: estimateSizeFactor) . Genes were considered expressed if the number of reads was above the 25th percentile for the normalized data set. For quality control, only replicates with Pearson correlation coefficient above 0.9 on their FPKM values were considered (Figure S1). Expressed genes were then further analyzed for differential gene expression using the DEseq2 package with FDR cutoff: 0.1. Differential gene expression was calculated in WT BMDMs stimulated 10:1 with Bb relative to unstimulated WT BMDMs and MyD88−/− BMDMs stimulated 100:1 with Bb relative to unstimulated MyD88−/− BMDMs. Differentially expressed genes (DEGs) were classified as either up-regulated or down-regulated based on the log2 of the fold change compared to the unstimulated control, which was calculated in R statistical software using package “DESeq2”. Determined DEGs were then separated into five groups based on their expression profiles; WT (all DEGs in WT BMDMs), MyD88−/− (all DEGs in MyD88−/− BMDMs), MyD88-dependent (all DEGs in WT but not MyD88−/− BMDMs), MyD88-independent (all DEGs in both WT and MyD88−/− BMDMs), and MyD88-privative (all DEGs in MyD88−/− but not in WT BMDMs).
Identification of enriched transcription binding sites and master regulator analysis
Transcription factor binding sites in promoters of differentially-expressed genes were analyzed using known DNA-binding motifs described in the TRANSFAC library , release 2017.2, available in the GeneXplain software (http://genexplain.com). Binding site enrichment analysis for each one of our sets of DEGs was carried out as part of a GeneXplain dedicated workflow. The background consisted of 300 mouse house-keeping genes and the TRANSFAC mouse Positional Weight Matrices PWM (motifs) for binding site prediction with p-value < 0.001 score cutoff. Promoters were extracted by the workflow with a length of 600 bp (− 500 to + 100) and an enrichment fold of 1.0.
Master regulatory molecules were searched for in signal transduction pathways upstream of the identified transcription factors. The GeneXplain workflow available for this analysis was used in conjunction with the GeneWays database. Parameters set included a maximum radius of 10 steps upstream of the transcription factor nodes, the DEG lists from the respective group as context genes and a z-score cutoff of 1.0. All transcription factors and master regulators used in the network analysis had confirmed expression in respective conditions using the total gene expression lists from the RNA-sequencing data set.
Gene ontology (GO) enrichment analysis
A Gene Ontology (GO) enrichment analysis was performed for the different sets of DEGs, transcription factors, and master regulators using the TRANSPATH  database through GeneXplain software. Input sets were the DEGs, transcription factors, or master regulators from either the MyD88-dependent, MyD88-independent, or MyD88-privative groups. Focus was directed to the GO biological processes output. GO biological processes related to Bb uptake, inflammation, and chemotaxis were identified by first reviewing previous studies for any genes involved in response to Bb relating to these phenotypes. Enrichment analysis was performed on these genes to identify GO biological processes that hit at least 60% of the genes on the list, generating a list of relevant GO biological processes. Then an intersection was performed between the list of GO biological processes identified using our DEG, transcription factor, or master regulator lists, and the GO biological processes identified from the relevant genes. Heat maps of expressed genes hits in each biological process were done in R statistical software using package “ggplots”.
Network reconstruction and network analysis
Networks were constructed joining the three identified layers on the networks: DEGs, transcription factors, and master regulators. The subnetworks were extracted from identified master regulators of interest. From the MyD88-dependent master regulator group, effort was directed on linking MyD88 with transcription factors that had binding sites in the promoter regions of the MyD88-dependent DEGs enriched in uptake biological processes. These transcription factors were identified using the TRANSPATH database with the enriched DEGs of interest as input. The output list of transcription factors was intersected with the list of transcription factors that were only expressed in WT BMDMs. Networks were assembled and analyzed using Cytoscape software . To extract the desired subnetworks, we used OCSANA  within the BiNOM plugin  in Cytoscape 2.8.3. MyD88 was considered as a source node and transcription factors from the intersected list as target nodes. For MyD88-privative chemotaxis subnetwork construction the same analysis pipeline was applied. MyD88-privative master regulators significantly enriched in chemotaxis were used as source nodes and MyD88-privative transcription factors enriched in chemotaxis were used as targets.
MyD88-deficient macrophages show comparable binding but reduced uptake of Bb
The macrophage is an essential cellular element of the human inflammatory response to the LD spirochete . Macrophages have also been shown as part of the inflammatory cell infiltrate in heart and joint tissue of mice experimentally infected with Bb [10, 53], and the importance of MyD88 in Bb clearance from mouse tissues has been previously reported [41, 54, 55]. It has also been well established that MyD88 enhances phagocytosis of multiple bacterial species by macrophages [28, 32, 34, 35, 56]. To better understand the contribution of MyD88 to spirochete binding, uptake and degradation by macrophages, we utilized an ex vivo macrophage model using WT and MyD88−/− BMDMs co-incubated with Bb at MOIs of either 10:1 or 100:1 for 1, 4 or 6 h. To quantify binding percentages, we imaged macrophages by confocal microscopy and determined the number of cells with spirochetes either attached to the surface or internalized because internalized spirochetes had to bind to macrophages before being taken up (Fig. 1a, yellow and white arrows respectively). We used the same confocal images and total cell numbers to quantify uptake percentages based on the number of cells with internalized spirochetes. The percentages of cells with spirochetes either bound or internalized were comparable between WT and MyD88−/− BMDMs at all three time points irrespective of MOI (Fig. 1b and c). While macrophages of both genotypes were able to phagocytose Bb, MyD88−/− BMDMs showed significantly reduced spirochete uptake compared to WT BMDMs at MOI 10:1 (Fig. 1d). Increasing the MOI to 100:1 significantly enhanced uptake in both cell genotypes, but MyD88−/− BMDMs never reached the phagocytic potential of their WT counterparts (Fig. 1e). These results further support the necessity of MyD88 signaling for efficient phagocytosis of Bb, irrespective of contact time with the spirochete.
TLR2, TLR7 and MyD88 are recruited to Bb-containing phagosomes in macrophages
Once spirochetes are phagocytosed by macrophages, recruitment of TLR and MyD88 proteins to the phagosome is essential to trigger MyD88-dependent signaling cascades [57,58,59,60]. Importantly, we have demonstrated that in human monocytes TLR2 and TLR8 co-localize to endosomes containing Bb . In addition, other investigators have shown a prominent role for TLR7 in the Bb inflammatory response . Murine TLR8, unlike murine TLR7 and human TLR8, does not seem to utilize ssRNA as its ligand . We therefore next characterized co-localization of TLR2, TLR7 and MyD88 with phagosomes containing Bb in BMDMs. By confocal microscopy, we observed that in WT BMDMs there is colocalization of MyD88 (Fig. 2a), TLR2 (Fig. 2b) and TLR7 (Fig. 2c) with Bb-containing phagosomes. Signals from MyD88 and TLR2 distinctly overlap with Bb GFP signals from phagosomes showing evidence of coiled or degraded spirochetes (Fig. 2a and b, graphs), but the intensity of MyD88 or TLR2 signal observed was higher with phagosomes containing degraded spirochetes. We also noted that TLR2 was expressed on the cell membrane and showed colocalization with surface-bound spirochetes (Fig. 2b). The absence of fluorescence in controls with secondary antibody only confirmed that this colocalization was not due to spectral overlap between color channels (Figure S2). TLR7 only showed strong signal with phagosomes containing partially degraded Bb but did not colocalize with surface-bound or recently internalized spirochetes (Fig. 2c). Taken together, these data confirm that endosomal TLR2, TLR7 and MyD88 colocalize to Bb-containing phagosomes to facilitate recognition of bacterial ligands and early response to infection.
Lack of MyD88 does not affect degradation of Bb in the phagosome
Degradation of the spirochete in the phagosome is crucial to expose bacterial ligands for recognition by endosomal TLRs . This process, known as phagosome maturation, requires reduction of phagosome pH and fusion with lysosomes . Given that both WT and MyD88−/− BMDMs bind and internalize Bb, we next sought to determine if spirochetes are similarly degraded in phagosomes with and without MyD88. Confocal images taken after a 6-h stimulation at MOI 10:1 showed that both WT and MyD88−/− BMDMs contained degraded GFP+ Bb within the cell actin matrix (Fig. 3a and b). To assess phagosome maturation, we quantitated recruitment of LAMP-1 to Bb-containing phagosomes by looking at colocalization of LAMP-1 and GFP fluorescence intensity . Both WT and MyD88−/− BMDMs showed comparable LAMP-1 and Bb colocalization in phagosomes (Fig. 3a and b, graphs). Colocalization between Bb and LAMP-1 was measured in multiple phagosomes in BMDMs from both genotypes and no significant differences were found (Fig. 3c). To confirm MyD88 signaling in response to Bb we also measured cytokine secretion after 1, 4 and 6 h of incubation with spirochetes. WT BMDMs showed significant increase in IL-6, TNFα and IL-10 secretion in the presence of spirochetes, whereas MyD88−/− BMDMs did not (Figure S3A-C). Consistent with prior studies by Behera et al. (2006), both WT and MyD88−/− BMDMs secrete the macrophage chemokine CCL2 (Figure S3D).
Bb ligand recognition appears to occur solely from within the phagosome
To test for the presence of bacterial products in the cytosol, we measured cleaved caspase-1, which is indicative of inflammasome activation. Western blot analysis of WT BMDM cell lysates and supernatants showed no activation of caspase-1 by stimulation of Bb alone (Fig. 3d), which is consistent with previously published studies . However, in discordance with previous studies , we did not see cleavage of IL-1β (Fig. 3d) unless exogenous ATP was added to the stimulation. To further confirm lack of NLRP3 inflammasome activation, we assessed Apoptosis-associated speck-like protein containing a CARD (ASC) in BMDMs stimulated with either Bb or Staphylococcus aureus (Sa) for 30 min or 6 h (Figure S4). As previously reported  (Figure S4A and C), ASC activation was observed with Sa, but no ASC was observed in BMDMs stimulated with Bb at 30 min or 6 h (Figure S4B and D). Thus, recognition of Bb ligands appears to occur solely within the phagosome.
MyD88-dependent signaling causes differential expression of genes in macrophages that promote the inflammatory response
Our results above show that MyD88 expression in macrophages enhances their capacity to phagocytose spirochetes (Fig. 1). To gain a better understanding of events that occur downstream of signaling by MyD88 which result in this phenotype presentation, we performed RNA-sequencing on WT and MyD88−/− BMDMs stimulated with Bb for 6 h. This time point was selected based on our data in Fig. 3 showing comparable maturation in both WT and MyD88−/− BMDM phagosomes. We sequenced RNA from WT BMDMs at a MOI of 10:1 and MyD88−/− BMDMs at a MOI of 100:1 for a comparative analysis because the uptake percentages were not significantly different between the two cell phenotypes under these conditions (Fig. 4a). Both WT and MyD88−/− BMDMs showed differentially expressed genes (DEGs) when compared to their respective unstimulated controls. We noted that the number of DEGs in WT BMDMs was much higher than in MyD88−/− BMDMs (2818 genes vs 141 genes respectively) (Fig. 4b). We saw similar numbers of up- and down-regulated DEGs in WT BMDMs (52 and 48%) (Fig. 4b). In the MyD88−/− BMDMs, approximately 83% of the DEGs were up-regulated (Fig. 4b). We classified the DEGs into three categories for further analysis: genes differentially expressed only in WT BMDMs (MyD88-dependent); genes differentially expressed in both WT and MyD88−/− BMDMs (MyD88-independent); and genes that were differentially expressed only in MyD88−/− BMDMs (MyD88-privative) (Fig. 4c).
Similar inflammatory and chemotactic processes are enriched regardless of MyD88-mediated signaling but utilize different regulatory proteins
MyD88-dependent mechanisms of inflammation have been well characterized, but little work has been done to understand the drivers of Bb-induced inflammation in the absence of MyD88. To address this issue, we next completed a comprehensive bioinformatics analysis to gain insight into how the DEGs are regulated within Bb-infected macrophages, both in the presence or absence of MyD88. We first identified transcription factors with potential binding sites in the promoter regions of the DEGs for each of the three subsets (66 for MyD88-dependent, 201 for MyD88-independent, and 39 for MyD88-privative). We then identified master regulator proteins upstream of these transcription factors and performed a Gene Ontology (GO) enrichment analysis of each group. Because data shown in Fig. 1 and Figure S3 indicate that in macrophages MyD88 affects both the inflammatory response and uptake of spirochetes, we focused our analysis on identifying whether any master regulators enriched to inflammatory and/or phagocytic biological processes in the MyD88-dependent and -privative conditions. Interestingly, similar inflammatory biological processes enriched to both the MyD88-dependent (including MyD88, Irak2 and Ly96) and MyD88-privative (including Vcam1 and Cxcl2) master regulators (Fig. 4d and e), but the individual master regulators involved were different for each subset (Fig. 4e). Importantly, over three times as many master regulators were identified for the MyD88-dependent DEGs than the MyD88-privative DEGs (Fig. 4e), suggesting that MyD88 signaling controls activation of more master regulators in the cell to control expression DEGs and enables the cell to perform unique processes in response to bacterial pathogens such as Bb.
MyD88-privative master regulators are involved in multiple chemotactic biological processes not enriched in WT BMDMs
We also observed significant overlap between the chemotactic biological processes enriched in MyD88-dependent and MyD88-privative master regulators. However, MyD88-privative master regulators significantly enriched to multiple biological processes involved in chemotaxis that were not enriched in MyD88-dependent master regulators (Figure S5A), suggesting that the lack of MyD88 signaling allows for increased up-regulation of processes to facilitate cell migration into the tissues. The MyD88-privative master regulators involved in these chemotactic processes also enriched to inflammatory processes (Figure S5B and Fig. 4e), suggesting that Bb may trigger other signaling cascades which induce inflammation more skewed to cell recruitment and localization.
MyD88 is a master regulator for transcription factors that control the MyD88-dependent DEGs enriched in uptake processes
Based on our observation that the presence of MyD88 enhances phagocytosis (Fig. 1), we also analyzed whether any of the MyD88-dependent DEGs enriched to biological processes related to uptake. We identified 164 MyD88-dependent DEGs that enriched to five different biological processes relating to phagocytosis (Actin Filament Polymerization, Regulation of Cell Shape, Actin Cytoskeleton Organization, Cytoskeleton Organization, and Actin Filament Organization). Of particular interest, Daam1 and Fmnl1, encoding two proteins known to play a role in phagocytosis of Bb [15, 37], were differentially expressed in an MyD88-dependent manner. Daam1, which was up-regulated, is a formin protein that bundles actin fibers together to increase stability of coiling pseudopods, which are more adept at capturing the highly motile spirochetes . In contrast to Daam1, Fmnl1 was down-regulated in response to Bb. Fmnl1 is also a formin protein that severs actin branches to promote polymerization and increase filopodia protrusion . To determine whether MyD88 is a master regulator in any of these processes, we first identified transcription factors that map to promoter regions of the enriched DEGs. Analysis of these transcription factors revealed that Zic1 and Zeb1 have the capacity to bind to the promoter regions of several of the MyD88-dependent DEGs that significantly enriched to processes associated with bacterial uptake (Fig. 5). Zic1 is controlled by the intermediate protein ApoE, which is known to play a role in cholesterol metabolism in macrophages  and the absence of ApoE increases Bb burdens in experimentally infected mice . We then used OCSANA, a specialized package available in Cytoscape  to link MyD88, as a master regulator, with transcription factors that map to DEGs in this specific subset. Based on this information we constructed a network illustrating potential links between MyD88-mediated signaling and up-regulated DEGs that may contribute to enhanced phagocytic capability seen in WT cells. The network (Fig. 5) shows Rhoa, Akt1, Rac1 and Cdc42 as genes that code for proteins which appear as intermediates on the network, meaning that their genes weren’t differentially expressed in our analysis. Daam1 regulates Rhoa activity, which controls Cdc42, Rac1 and Akt1. Cdc42 activates a Rho GTPase, Rhoq, which is up-regulated in response to Bb. Rac1 and Akt1, when translated, both activate multiple proteins whose corresponding genes are also up-regulated, indicating that while the genes for these intermediate proteins aren’t differentially expressed, they are still active in macrophages that have been stimulated with Bb. Taken together, these data suggest that MyD88 signaling upregulates multiple gene products involved in regulating macrophage membrane protrusions. Upregulation of these genes likely contributes to the reorganization of cell machinery that enhances the capability of the WT macrophage to take up spirochetes.
Previous studies by our group have emphasized that uptake and degradation of Bb by phagocytic cells, including monocytes and macrophages, are critical in eliciting the inflammatory response to the bacterium [14, 17,18,19, 29]. The findings from these studies, as well as others , show that the adaptor protein MyD88 plays a critical role in bacterial uptake and phagosomal signaling in macrophages. In the current study, we provide further evidence that the macrophage is a key driver of inflammation, even in the absence of MyD88. We also show that while MyD88 has a significant impact on spirochetal uptake, phagosome maturation and bacterial degradation are not affected (Fig. 1). Of particular novelty, phagosomal signaling cascades induced by Bb ligands in macrophages trigger a number of inflammatory and chemotactic pathways (Figs. 4 and S5). Moreover, the inflammatory processes are mediated by different regulatory proteins depending on whether MyD88 is present or absent, while induction of several chemotactic processes occurs independently of MyD88. In-depth analysis of these signaling cascades allowed us to identify previously underappreciated MyD88-dependent transcription factors which could lead to enhanced spirochetal uptake and clearance (Fig. 5).
To better understand the contribution of MyD88 to spirochete binding, uptake, degradation and signaling by macrophages, we used an ex vivo model. Murine macrophages lacking MyD88 show a phagocytic defect when stimulated with Bb ex vivo compared to WT macrophages. This defect in uptake has been previously demonstrated in macrophage stimulation experiments with other bacteria strains [31, 34, 35, 38]. Our results reveal that binding of Bb is not affected and that this phagocytic defect is not dependent on length of stimulation (i.e. time dependent) and is only slightly rescued by increasing the MOI. Thus, in the absence of MyD88, macrophages are still capable of binding and taking up the LD spirochete, but MyD88 signaling enhances the efficiency of Bb phagocytosis by macrophages. We also show here that when stimulated ex vivo, the macrophage response to Bb is driven by the signaling cascades induced by MyD88 as a result of bacterial ligands engaging TLR2 and TLR7 receptors in the phagosome (Figs. 2 and 4). Recognition in the phagosome is driven by degradation of bacteria, since more TLR2, TLR7 and MyD88 marker intensity were observed colocalizing with degraded spirochetes. Our results in Fig. 3 indicate that in the context of Bb infection, MyD88 is not required for phagosome maturation, evidenced by the recruitment of LAMP-1 to Bb-containing phagosomes in MyD88−/− BMDMs. This is in contrast to Blander et al., who published that MyD88−/− BMDMs infected with Sa or E. coli did not colocalize with either Lysotracker or LAMP-1 to the same degree as WT BMDMs . One possible explanation for the different findings with our study is that the recruitment of LAMP-1 is delayed in MyD88−/− BMDMs, given that Blander et al. measured phagosome maturation at an earlier time point than in our studies. This explanation is supported by Yates et al. who showed slightly delayed acidification in MyD88−/− BMDMs stimulated with TLR2 or TLR4 ligands for 40 min . The fragility of the Bb membranes also suggests that perhaps less acidification of the phagosome is needed to expose Bb PAMPs.
Reduced uptake of bacteria in macrophages lacking MyD88 is a phenotypic trait that has been extensively detailed [28, 31,32,33,34,35], but not well understood. To better understand the relationship between MyD88 signaling and phagocytosis, we used a computational systems biology approach. In prior studies, addition of TLR3 ligands to Bb stimulation of MyD88−/− BMDMs significantly rescues uptake , suggesting that in the absence of MyD88 TRIF signaling can activate pathways that result in similar actin rearrangement in the cell. This signaling was shown to be mediated through PI3K , but interestingly PI3K was not differentially expressed in our macrophage stimulation. However, our network analysis from RNA-sequencing data identified DEGs that are up-regulated downstream of common phagocytosis effector proteins (Fig. 5). Rhoq, activated by Cdc42 codes for TC10, a protein involved in generating long filopodia protrusions . The gene Cyfip1, which encodes a part of the WAVE complex that regulates actin polymerization , was also up-regulated according to the network through Rac1 protein interactions. The WAVE complex has higher involvement with lamellipodia formations . It is likely that MyD88 controls transcription factors that upregulate these genes to promote phagocytosis through formation of coiling pseudopods, which are more similar to lamellipodia, rather than through straight filopodia protrusions. In addition to MyD88, studies indicating that TLR2 can utilize TRIF have also been completed, but this interaction only appears to contribute to the inflammatory response rather than spirochete uptake [18, 73]. More recently, the leukotriene LTB4 has been shown to promote phagocytosis of Bb by macrophages , but in our BMDM sequencing data we did not find differential expression of Ltb4 or its receptor Ltb4r1. This could possibly be due to the later time point we selected for sequencing. It has also been shown that spleen tyrosine kinase (Syk) has an important role in phagocytosis of Bb via integrin binding . The Syk gene (Syk) is significantly up-regulated in an MyD88-dependent manner, suggesting that MyD88 drives over-expression of Syk to increase phosphorylation and activation of proteins involved in generating actin branches. However, in our GO analysis Syk was not one of the 164 MyD88-dependent genes that enriched to uptake biological processes, and the transcription factor Zic1, which has binding sites in the promoter regions of a significant number of these genes, is not predicted to bind in the promoter region of Syk. Zic1 was of particular interest to us because it appeared downstream of MyD88 in our network analysis (Fig. 5) and is controlled by the intermediate protein ApoE. Mice lacking ApoE have increased bacterial burdens when infected with Bb , suggesting that ApoE signaling plays a role in cell remodeling processes necessary to enhance uptake. In addition, a link between Bb phagocytosis and cholesterol has been postulated by Hawley et al. who showed that CR3, a known phagocytic receptor for Bb, is recruited to lipid rafts with the co-receptor CD14 . Thus, it is possible that MyD88 upregulates ApoE to enhance lipid rafts on the macrophage membrane, which can potentiate signaling to enhance uptake and provide scaffolding for proteins involved in actin remodeling.
Our computational analysis also supported that there are non-canonical sources of inflammation in MyD88−/− mice. Our results suggest that there is possibly another receptor recruited to the phagosome that initiates chemokine production upon recognition of a Bb ligand. Another mechanism for triggering chemokine production may be that the TLR receptors are utilizing another adaptor protein to transmit signals out of the phagosome, as postulated by Petnicki-Ocwieja et al. . Network analysis of DEGs from macrophages identified multiple master regulators that could be controlling production of these chemokines, but further investigation is needed to determine if these master regulators are in fact active in macrophages containing Bb. Additional studies to test whether acidification of the phagosome is required for Bb-induced chemokine production will also give insight into which ligand-receptor interaction induces this response.
In summary, our results emphasize that the macrophage has a very important role in both recognition and clearance of Bb and is at the epicenter of the immunologic response to spirochete infection. The findings from these studies have also advanced our understanding of how phagosomal signaling drives spirochete uptake, recognition and inflammation. The adaptor protein MyD88 plays a critical role in these processes. Initial phagocytosis of Bb by macrophages does not require MyD88, but once taken up, recognition of Bb ligands exposed upon spirochete degradation occurs through endosomal TLRs which trigger MyD88-mediated signaling cascades. This signaling results in cell remodeling to enhance phagocytosis, as indicated by our ex vivo data, and allows macrophages to more efficiently internalize and clear the highly motile spirochetes by using more dynamic membrane protrusions. Further studies using the targets identified in these experiments may also provide insight into understanding the importance of phagocytosis in other bacterial infections. We can use similar techniques to look at the role of the macrophage response and MyD88 signaling in human macrophages, with the goal of increasing our understanding of the clinical spectrum associated with Lyme disease pathogenesis.
Availability of data and materials
Catalog numbers and suppliers for commercially available materials used in these studies are provided. Non-commercially available materials can be provided upon request. ImageJ software can be downloaded at no cost from nih.gov. Cytoscape can be downloaded at no cost from cytoscape.org. GeneXplain is a subscription software that can be purchased from genexplain.com.
- Bb :
Pathogen associated molecular patterns
Green fluorescent protein
Barbour Stoner Kelly
- Sa :
Bone marrow derived macrophages
Bovine serum albumin
Differentially expressed genes
Multiplicities of infection
Apoptosis-associated speck-like protein containing a CARD
Benach JL, Bosler EM, Hanrahan JP, Coleman JL, Habicht GS, Bast TF, et al. Spirochetes isolated from the blood of two patients with Lyme disease. N Engl J Med. 1983;308(13):740–2. PubMed PMID: 6828119. https://doi.org/10.1056/NEJM198303313081302.
CDC. Lyme Disease Incidence Rates by State 2004-2015 2016. Available from: http://www.cdc.gov/lyme/stats/chartstables/incidencebystate.html.
Steere AC. Lyme disease. N Engl J Med. 1989;321(9):586–96. PubMed PMID: 2668764. Epub 1989/08/31. https://doi.org/10.1056/NEJM198908313210906.
CDC. Signs and Symptoms of Untreated Lyme Disease 2016. Available from: http://www.cdc.gov/lyme/signs_symptoms/index.html.
Steere AC, Strle F, Wormser GP, Hu LT, Branda JA, Hovius JW, et al. Lyme borreliosis. Nat Rev Dis Primers. 2016;2:16090 PubMed PMID: 27976670. PMCID: PMC5539539. Epub 2016/12/16.
Bockenstedt JJWaLK. Host Response. In: Radolf DSSaJD, editor. Borrelia: molecular biology, host interaction and pathogenesis. Norfolk: Caister Academic Press; 2010. p. 413–41.
Salazar JC, Pope CD, Sellati TJ, Feder HM Jr, Kiely TG, Dardick KR, et al. Coevolution of markers of innate and adaptive immunity in skin and peripheral blood of patients with erythema migrans. J Immunol. 2003;171(5):2660–70. PubMed PMID: 12928420. https://doi.org/10.4049/jimmunol.171.5.2660.
Shin JJ, Strle K, Glickstein LJ, Luster AD, Steere AC. Borrelia burgdorferi stimulation of chemokine secretion by cells of monocyte lineage in patients with Lyme arthritis. Arthritis Res Ther. 2010;12(5):R168 PubMed PMID: 20828409. PMCID: 2990995.
Lasky CE, Olson RM, Brown CR. Macrophage Polarization during Murine Lyme Borreliosis. Infect Immun. 2015;83(7):2627–35 PubMed PMID: 25870230. PMCID: PMC4468556. Epub 2015/04/15.
Montgomery RR, Booth CJ, Wang X, Blaho VA, Malawista SE, Brown CR. Recruitment of macrophages and polymorphonuclear leukocytes in Lyme carditis. Infect Immun. 2007;75(2):613–20 PubMed PMID: 17101663. PMCID: PMC1828503. Epub 2006/11/15.
Hawley KL, Olson CM Jr, Iglesias-Pedraz JM, Navasa N, Cervantes JL, Caimano MJ, et al. CD14 cooperates with complement receptor 3 to mediate MyD88-independent phagocytosis of Borrelia burgdorferi. Proc Natl Acad Sci U S A. 2012;109(4):1228–32 PubMed PMID: 22232682. PMCID: 3268315.
Cinco M, Murgia R, Presani G, Perticarari S. Integrin CR3 mediates the binding of nonspecifically opsonized Borrelia burgdorferi to human phagocytes and mammalian cells. Infect Immun. 1997;65(11):4784–9 PubMed PMID: 9353065. PMCID: PMC175686. Epub 1997/11/14.
Behera AK, Hildebrand E, Uematsu S, Akira S, Coburn J, Hu LT. Identification of a TLR-independent pathway for Borrelia burgdorferi-induced expression of matrix metalloproteinases and inflammatory mediators through binding to integrin alpha 3 beta 1. J Immunol. 2006;177(1):657–64. PubMed PMID: 16785564. https://doi.org/10.4049/jimmunol.177.1.657.
Cervantes JL, Hawley KL, Benjamin SJ, Weinerman B, Luu SM, Salazar JC. Phagosomal TLR signaling upon Borrelia burgdorferi infection. Front Cell Infect Microbiol. 2014;4:55 PubMed PMID: 24904837. PMCID: PMC4033037.
Naj X, Hoffmann AK, Himmel M, Linder S. The formins FMNL1 and mDia1 regulate coiling phagocytosis of Borrelia burgdorferi by primary human macrophages. Infect Immun. 2013;81(5):1683–95 PubMed PMID: 23460512. PMCID: PMC3647995. Epub 2013/03/06.
Naj X, Linder S. ER-coordinated activities of Rab22a and Rab5a drive Phagosomal compaction and intracellular processing of Borrelia burgdorferi by macrophages. Cell Rep. 2015;12(11):1816–30. PubMed PMID: 26344766. Epub 2015/09/08. https://doi.org/10.1016/j.celrep.2015.08.027.
Moore MW, Cruz AR, LaVake CJ, Marzo AL, Eggers CH, Salazar JC, et al. Phagocytosis of Borrelia burgdorferi and Treponema pallidum potentiates innate immune activation and induces gamma interferon production. Infect Immun. 2007;75(4):2046–62 PubMed PMID: 17220323. PMCID: 1865718.
Salazar JC, Duhnam-Ems S, La Vake C, Cruz AR, Moore MW, Caimano MJ, et al. Activation of human monocytes by live Borrelia burgdorferi generates TLR2-dependent and -independent responses which include induction of IFN-beta. PLoS Pathog. 2009;5(5):e1000444 PubMed PMID: 19461888. PMCID: 2679197.
Cervantes JL, Dunham-Ems SM, La Vake CJ, Petzke MM, Sahay B, Sellati TJ, et al. Phagosomal signaling by Borrelia burgdorferi in human monocytes involves Toll-like receptor (TLR) 2 and TLR8 cooperativity and TLR8-mediated induction of IFN-beta. Proc Natl Acad Sci U S A. 2011;108(9):3683–8 PubMed PMID: 21321205. PMCID: 3048123.
Brandt ME, Riley BS, Radolf JD, Norgard MV. Immunogenic integral membrane proteins of Borrelia burgdorferi are lipoproteins. Infect Immun. 1990;58(4):983–91 PubMed PMID: 2318538. PMCID: PMC258571.
Schenk M, Belisle JT, Modlin RL. TLR2 looks at lipoproteins. Immunity. 2009;31(6):847–9. PubMed PMID: 20064441. Epub 2010/01/13. https://doi.org/10.1016/j.immuni.2009.11.008.
Hirschfeld M, Kirschning CJ, Schwandner R, Wesche H, Weis JH, Wooten RM, et al. Cutting edge: inflammatory signaling by Borrelia burgdorferi lipoproteins is mediated by toll-like receptor 2. J Immunol. 1999;163(5):2382–6 PubMed PMID: 10452971.
Wooten RM, Ma Y, Yoder RA, Brown JP, Weis JH, Zachary JF, et al. Toll-like receptor 2 is required for innate, but not acquired, host defense to Borrelia burgdorferi. J Immunol. 2002;168(1):348–55. PubMed PMID: 11751980. https://doi.org/10.4049/jimmunol.168.1.348.
Radolf JD, Norgard MV, Brandt ME, Isaacs RD, Thompson PA, Beutler B. Lipoproteins of Borrelia burgdorferi and Treponema pallidum activate cachectin/tumor necrosis factor synthesis. Analysis using a CAT reporter construct. J Immunol. 1991;147(6):1968–74 PubMed PMID: 1890308.
Alexopoulou L, Thomas V, Schnare M, Lobet Y, Anguita J, Schoen RT, et al. Hyporesponsiveness to vaccination with Borrelia burgdorferi OspA in humans and in TLR1- and TLR2-deficient mice. Nat Med. 2002;8(8):878–84. PubMed PMID: 12091878. Epub 2002/07/02. https://doi.org/10.1038/nm732.
Radolf JD, Arndt LL, Akins DR, Curetty LL, Levi ME, Shen Y, et al. Treponema pallidum and Borrelia burgdorferi lipoproteins and synthetic lipopeptides activate monocytes/macrophages. J Immunol. 1995;154(6):2866–77 PubMed PMID: 7876555.
Radolf JD, Goldberg MS, Bourell K, Baker SI, Jones JD, Norgard MV. Characterization of outer membranes isolated from Borrelia burgdorferi, the Lyme disease spirochete. Infect Immun. 1995;63(6):2154–63 PubMed PMID: 7768594. PMCID: PMC173280.
Shin OS, Isberg RR, Akira S, Uematsu S, Behera AK, Hu LT. Distinct roles for MyD88 and Toll-like receptors 2, 5, and 9 in phagocytosis of Borrelia burgdorferi and cytokine induction. Infect Immun. 2008;76(6):2341–51 PubMed PMID: 18378636. PMCID: 2423091.
Cervantes JL, La Vake CJ, Weinerman B, Luu S, O'Connell C, Verardi PH, et al. Human TLR8 is activated upon recognition of Borrelia burgdorferi RNA in the phagosome of human monocytes. J Leukoc Biol. 2013;94(6):1231–41 PubMed PMID: 23906644. PMCID: 3828603.
Akira S, Takeda K, Kaisho T. Toll-like receptors: critical proteins linking innate and acquired immunity. Nat Immunol. 2001;2(8):675–80. PubMed PMID: 11477402. Epub 2001/07/31. https://doi.org/10.1038/90609.
Edelson BT, Unanue ER. MyD88-dependent but toll-like receptor 2-independent innate immunity to Listeria: no role for either in macrophage listericidal activity. J Immunol. 2002;169(7):3869–75. PubMed PMID: 12244184. https://doi.org/10.4049/jimmunol.169.7.3869.
Blander JM. Coupling toll-like receptor signaling with phagocytosis: potentiation of antigen presentation. Trends Immunol. 2007;28(1):19–25. PubMed PMID: 17126600. https://doi.org/10.1016/j.it.2006.11.001.
Takeuchi O, Hoshino K, Akira S. Cutting edge: TLR2-deficient and MyD88-deficient mice are highly susceptible to Staphylococcus aureus infection. J Immunol. 2000;165(10):5392–6. PubMed PMID: 11067888. Epub 2000/11/09. https://doi.org/10.4049/jimmunol.165.10.5392.
Ip WK, Sokolovska A, Charriere GM, Boyer L, Dejardin S, Cappillino MP, et al. Phagocytosis and phagosome acidification are required for pathogen processing and MyD88-dependent responses to Staphylococcus aureus. J Immunol. 2010;184(12):7071–81 PubMed PMID: 20483752. PMCID: PMC2935932.
Shen Y, Kawamura I, Nomura T, Tsuchiya K, Hara H, Dewamitta SR, et al. Toll-like receptor 2- and MyD88-dependent phosphatidylinositol 3-kinase and Rac1 activation facilitates the phagocytosis of Listeria monocytogenes by murine macrophages. Infect Immun. 2010;78(6):2857–67 PubMed PMID: 20368346. PMCID: PMC2876581.
Shin OS, Miller LS, Modlin RL, Akira S, Uematsu S, Hu LT. Downstream signals for MyD88-mediated phagocytosis of Borrelia burgdorferi can be initiated by TRIF and are dependent on PI3K. J Immunol. 2009;183(1):491–8 PubMed PMID: 19542460. PMCID: PMC2772066. Epub 2009/06/23.
Hoffmann AK, Naj X, Linder S. Daam1 is a regulator of filopodia formation and phagocytic uptake of Borrelia burgdorferi by primary human macrophages. FASEB J. 2014;28(7):3075–89. PubMed PMID: 24696301. Epub 2014/04/04. https://doi.org/10.1096/fj.13-247049.
Blander JM, Medzhitov R. Regulation of phagosome maturation by signals from toll-like receptors. Science. 2004;304(5673):1014–8. PubMed PMID: 15143282. https://doi.org/10.1126/science.1096158.
Shen N, Fu Q, Deng Y, Qian X, Zhao J, Kaufman KM, et al. Sex-specific association of X-linked Toll-like receptor 7 (TLR7) with male systemic lupus erythematosus. Proc Natl Acad Sci U S A. 2010;107(36):15838–43 PubMed PMID: 20733074. PMCID: PMC2936646. Epub 2010/08/25.
Adachi O, Kawai T, Takeda K, Matsumoto M, Tsutsui H, Sakagami M, et al. Targeted disruption of the MyD88 gene results in loss of IL-1- and IL-18-mediated function. Immunity. 1998;9(1):143–50. PubMed PMID: 9697844. Epub 1998/08/11. https://doi.org/10.1016/S1074-7613(00)80596-8.
Liu N, Montgomery RR, Barthold SW, Bockenstedt LK. Myeloid differentiation antigen 88 deficiency impairs pathogen clearance but does not alter inflammation in Borrelia burgdorferi-infected mice. Infect Immun. 2004;72(6):3195–203 PubMed PMID: 15155621. PMCID: 415708.
Steere AC, Grodzicki RL, Craft JE, Shrestha M, Kornblatt AN, Malawista SE. Recovery of Lyme disease spirochetes from patients. Yale J Biol Med. 1984;57(4):557–60 PubMed PMID: 6393606. PMCID: PMC2589990. Epub 1984/07/01.
Dunham-Ems SM, Caimano MJ, Pal U, Wolgemuth CW, Eggers CH, Balic A, et al. Live imaging reveals a biphasic mode of dissemination of Borrelia burgdorferi within ticks. J Clin Invest. 2009;119(12):3652–65 PubMed PMID: 19920352. PMCID: 2786795.
Sokolovska A, Becker CE, Stuart LM. Measurement of phagocytosis, phagosome acidification, and intracellular killing of Staphylococcus aureus. Curr Protoc Immunol. 2012;Chapter 14(Unit14):30 PubMed PMID: 23129153. Epub 2012/11/07.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36 PubMed PMID: 23618408. PMCID: PMC4053844. Epub 2013/04/27.
Lister R, Gregory BD, Ecker JR. Next is now: new technologies for sequencing of genomes, transcriptomes, and beyond. Curr Opin Plant Biol. 2009;12(2):107–18 PubMed PMID: 19157957. PMCID: PMC2723731. Epub 2009/01/23.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106 PubMed PMID: 20979621. PMCID: PMC3218662. Epub 2010/10/29.
Matys V, Kel-Margoulis OV, Fricke E, Liebich I, Land S, Barre-Dirrie A, et al. TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 2006;34(Database issue):D108–10 PubMed PMID: 16381825. PMCID: PMC1347505. Epub 2005/12/31.
Krull M, Pistor S, Voss N, Kel A, Reuter I, Kronenberg D, et al. TRANSPATH: an information resource for storing and visualizing signaling pathways and their pathological aberrations. Nucleic Acids Res. 2006;34(Database issue):D546–51 PubMed PMID: 16381929. PMCID: PMC1347469. Epub 2005/12/31.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504 PubMed PMID: 14597658. PMCID: PMC403769. Epub 2003/11/05.
Vera-Licona P, Bonnet E, Barillot E, Zinovyev A. OCSANA: optimal combinations of interventions from network analysis. Bioinformatics. 2013;29(12):1571–3. PubMed PMID: 23626000. Epub 2013/04/30. https://doi.org/10.1093/bioinformatics/btt195.
Zinovyev A, Viara E, Calzone L, Barillot E. BiNoM: a Cytoscape plugin for manipulating and analyzing biological networks. Bioinformatics. 2008;24(6):876–7. PubMed PMID: 18024474. Epub 2007/11/21. https://doi.org/10.1093/bioinformatics/btm553.
Hawley K, Navasa N, Olson CM Jr, Bates TC, Garg R, Hedrick MN, et al. Macrophage p38 mitogen-activated protein kinase activity regulates invariant natural killer T-cell responses during Borrelia burgdorferi infection. J Infect Dis. 2012;206(2):283–91 PubMed PMID: 22551807. PMCID: PMC3490691. Epub 2012/05/04.
Bolz DD, Sundsbak RS, Ma Y, Akira S, Kirschning CJ, Zachary JF, et al. MyD88 plays a unique role in host defense but not arthritis development in Lyme disease. J Immunol. 2004;173(3):2003–10. PubMed PMID: 15265935. https://doi.org/10.4049/jimmunol.173.3.2003.
Behera AK, Hildebrand E, Bronson RT, Perides G, Uematsu S, Akira S, et al. MyD88 deficiency results in tissue-specific changes in cytokine induction and inflammation in interleukin-18-independent mice infected with Borrelia burgdorferi. Infect Immun. 2006;74(3):1462–70 PubMed PMID: 16495516. PMCID: 1418660.
Yates RM, Russell DG. Phagosome maturation proceeds independently of stimulation of toll-like receptors 2 and 4. Immunity. 2005;23(4):409–17. PubMed PMID: 16226506. https://doi.org/10.1016/j.immuni.2005.09.007.
Underhill DM, Ozinsky A, Hajjar AM, Stevens A, Wilson CB, Bassetti M, et al. The toll-like receptor 2 is recruited to macrophage phagosomes and discriminates between pathogens. Nature. 1999;401(6755):811–5. PubMed PMID: 10548109. Epub 1999/11/05. https://doi.org/10.1038/44605.
Petnicki-Ocwieja T, Kern A, Killpack TL, Bunnell SC, Hu LT. Adaptor Protein-3-Mediated Trafficking of TLR2 Ligands Controls Specificity of Inflammatory Responses but Not Adaptor Complex Assembly. J Immunol. 2015;195(9):4331–40 PubMed PMID: 26423153. PMCID: PMC4610863. Epub 2015/10/02.
Bergstrom B, Aune MH, Awuh JA, Kojen JF, Blix KJ, Ryan L, et al. TLR8 senses Staphylococcus aureus RNA in human primary monocytes and macrophages and induces IFN-beta production via a TAK1-IKKbeta-IRF5 signaling pathway. J Immunol. 2015;195(3):1100–11. PubMed PMID: 26085680. Epub 2015/06/19. https://doi.org/10.4049/jimmunol.1403176.
Nilsen NJ, Deininger S, Nonstad U, Skjeldal F, Husebye H, Rodionov D, et al. Cellular trafficking of lipoteichoic acid and Toll-like receptor 2 in relation to signaling: role of CD14 and CD36. J Leukoc Biol. 2008;84(1):280–91 PubMed PMID: 18458151. PMCID: PMC3178504.
Petzke MM, Brooks A, Krupna MA, Mordue D, Schwartz I. Recognition of Borrelia burgdorferi, the Lyme disease spirochete, by TLR7 and TLR9 induces a type I IFN response by human immune cells. J Immunol. 2009;183(8):5279–92. PubMed PMID: 19794067. https://doi.org/10.4049/jimmunol.0901390.
Heil F, Hemmi H, Hochrein H, Ampenberger F, Kirschning C, Akira S, et al. Species-specific recognition of single-stranded RNA via toll-like receptor 7 and 8. Science. 2004;303(5663):1526–9. PubMed PMID: 14976262. Epub 2004/02/21. https://doi.org/10.1126/science.1093620.
Blander JM, Medzhitov R. On regulation of phagosome maturation and antigen presentation. Nat Immunol. 2006;7(10):1029–35. PubMed PMID: 16985500. https://doi.org/10.1038/ni1006-1029.
Westphal A, Cheng W, Yu J, Grassl G, Krautkramer M, Holst O, et al. Lysosomal trafficking regulator Lyst links membrane trafficking to toll-like receptor-mediated inflammatory responses. J Exp Med. 2017;214(1):227–44 PubMed PMID: 27881733. PMCID: PMC5206490. Epub 2016/11/25.
Oosting M, van de Veerdonk FL, Kanneganti TD, Sturm P, Verschueren I, Berende A, et al. Borrelia species induce inflammasome activation and IL-17 production through a caspase-1-dependent mechanism. Eur J Immunol. 2011;41(1):172–81 PubMed PMID: 21182088. PMCID: 3939779.
Oosting M, Berende A, Sturm P, Ter Hofstede HJ, de Jong DJ, Kanneganti TD, et al. Recognition of Borrelia burgdorferi by NOD2 is central for the induction of an inflammatory reaction. J Infect Dis. 2010;201(12):1849–58. PubMed PMID: 20441518. https://doi.org/10.1086/652871.
Sokolovska A, Becker CE, Ip WK, Rathinam VA, Brudner M, Paquette N, et al. Activation of caspase-1 by the NLRP3 inflammasome regulates the NADPH oxidase NOX2 to control phagosome function. Nat Immunol. 2013;14(6):543–53 PubMed PMID: 23644505. PMCID: PMC3708594.
Naj X, Linder S. Actin-dependent regulation of Borrelia burgdorferi phagocytosis by macrophages. Curr Top Microbiol Immunol. 2017;399:133–54. PubMed PMID: 27744511. Epub 2016/10/17. https://doi.org/10.1007/82_2016_26.
Toledo A, Monzon JD, Coleman JL, Garcia-Monco JC, Benach JL. Hypercholesterolemia and ApoE deficiency result in severe infection with Lyme disease and relapsing-fever Borrelia. Proc Natl Acad Sci U S A. 2015;112(17):5491–6 PubMed PMID: 25870274. PMCID: PMC4418910. Epub 2015/04/15.
Murphy GA, Solski PA, Jillian SA, de la Ossa PP, D'Eustachio P, Der CJ, et al. Cellular functions of TC10, a rho family GTPase: regulation of morphology, signal transduction and cell growth. Oncogene. 1999;18(26):3831–45. PubMed PMID: 10445846. https://doi.org/10.1038/sj.onc.1202758.
Levin R, Grinstein S, Canton J. The life cycle of phagosomes: formation, maturation, and resolution. Immunol Rev. 2016;273(1):156–79. PubMed PMID: 27558334. Epub 2016/08/26. https://doi.org/10.1111/imr.12439.
Takenawa T, Suetsugu S. The WASP-WAVE protein network: connecting the membrane to the cytoskeleton. Nat Rev Mol Cell Biol. 2007;8(1):37–48. PubMed PMID: 17183359. Epub 2006/12/22. https://doi.org/10.1038/nrm2069.
Petnicki-Ocwieja T, Chung E, Acosta DI, Ramos LT, Shin OS, Ghosh S, et al. TRIF mediates Toll-like receptor 2-dependent inflammatory responses to Borrelia burgdorferi. Infect Immun. 2013;81(2):402–10 PubMed PMID: 23166161. PMCID: PMC3553797. Epub 2012/11/21.
Zhang Y, Olson RM, Brown CR. Macrophage LTB4 drives efficient phagocytosis of Borrelia burgdorferi via BLT1 or BLT2. J Lipid Res. 2017;58(3):494–503 PubMed PMID: 28053185. PMCID: PMC5335579. Epub 2017/01/06.
Killpack TL, Ballesteros M, Bunnell SC, Bedugnis A, Kobzik L, Hu LT, et al. Phagocytic Receptors Activate Syk and Src Signaling during Borrelia burgdorferi Phagocytosis. Infect Immun. 2017;85(10) PubMed PMID: 28717031. PMCID: PMC5607427. Epub 2017/07/19.
Hawley KL, Martin-Ruiz I, Iglesias-Pedraz JM, Berwin B, Anguita J. CD14 targets complement receptor 3 to lipid rafts during phagocytosis of Borrelia burgdorferi. Int J Biol Sci. 2013;9(8):803–10 PubMed PMID: 23983613. PMCID: PMC3753444. Epub 2013/08/29.
The authors would like to thank Dr. Melissa Caimano, Dr. Ashley Groshong, Ms. Anna Allard, and Ms. Morgan Ledoyt for their contribution and support to execution of the experiments. We would also like to thank Oscar Luo for his input in the analysis of the RNA-sequencing data.
This work was funded by both the National Institutes of Health and Connecticut Children’s Medical Center. This work was supported by NIAID grants R01 AI090166 (JCS), the Robert E. Leet and Clara Guthrie Patterson Trust (KLH), and research funds provided by Connecticut Children’s (KLH, MJC, JDR and JCS). The funders had no role in study design, collection and analysis of data, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
All experiments in this work involving the use of animals or isolation of primary cells from animals were monitored and approved by the UConn Health Institutional Animal Care and Use Committee under protocol #101388–0819. The UConn Health Institutional Animal Care and Use Committee is accredited by The Association for Assessment and Accreditation of Laboratory Animal Care (AAALAC). Consent to participate not applicable.
Consent for publication
All authors have read and approved the manuscript.
The authors declared that they have no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Supplemental File 1.
Gene ontology analysis of differentially expressed genes.
Additional file 2: Supplemental File 2.
Gene ontology analysis of identified transcription factors.
Additional file 3: Supplemental File 3.
Gene ontology analysis of identified master regulators.
Additional file 4: Supplemental File 4.
Identified differentially expressed genes.
Additional file 5: Supplemental File 5.
Identified transcription factors.
Additional file 6: Supplemental File 6.
Identified master regulators.
Additional file 7: Figure S1.
Correlation of reads of RNA sequencing data. (A-D) Correlation of reads of RNA isolated from Bb-stimulated BMDMs. WT BMDMs were stimulated at an MOI 10:1 for 6 h (A). MyD88−/− BMDMs were stimulated at a MOI 100:1 for 6 h (B). To calculate differential expression, reads from the stimulated samples were normalized to unstimulated BMDMs of the same genotype; WT (C) or MyD88−/− (D).
Additional file 8: Figure S2.
Secondary Control for TLR2. Confocal 40x images of internalized Bb with AF350 Secondary Antibody in WT BMDMs after stimulation at MOI 10:1. Green is Bb, blue is AF350 Secondary Antibody, and red is actin.
Additional file 9: Figure S3.
MyD88−/− BMDMs stimulated with Bb show abrogated cytokine production. (A-D) Quantification of IL-6 (A), TNFα (B), IL-10 (C) and CCL2 (D) proteins in supernatant from WT and MyD88−/− (MyD) BMDMs stimulated with Bb at MOI 10:1 for 1, 4 or 6 h. N = 3 mouse BMDM per genotype *p-value< 0.05, **p-value< 0.01, ***p-value< 0.001, NS = not significant.
Additional file 10: Figure S4.
Bb does not induce ASC formation in BMDMs. (A-D) Confocal 40x images (40x) of BMDMs stimulated with either Sa (A and C) or Bb (B and D) for 30 min (A-B) or 6 h (C-D). Blue is nucleus, red is ASC and green is the bacteria species.
Additional file 11: Figure S5.
MyD88-privative MRs are significantly enriched in biological processes related to chemotaxis. (A) Venn diagram comparing biological processes (BP) relating to chemotaxis significantly enriched between MyD88-dependent (light gray) and MyD88-privative (dark gray) master regulators. (B) Heat map showing fold change of master regulators enriched in chemotaxis biological processes in MyD88−/− (yellow) BMDMs. GO numbers for significantly enriched BP are indicated on the x-axis.
Additional file 12.
Additional file 13.
Additional file 14.
Additional file 15.
Additional file 16.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Benjamin, S.J., Hawley, K.L., Vera-Licona, P. et al. Macrophage mediated recognition and clearance of Borrelia burgdorferi elicits MyD88-dependent and -independent phagosomal signals that contribute to phagocytosis and inflammation. BMC Immunol 22, 32 (2021). https://doi.org/10.1186/s12865-021-00418-8