Macrophage mediated recognition and clearance of Borrelia burgdorferi elicits MyD88-dependent and -independent phagosomal signals that contribute to phagocytosis and inflammation

doi:10.21203/rs.3.rs-279261/v1

PDF

Research article

Macrophage mediated recognition and clearance of Borrelia burgdorferi elicits MyD88-dependent and -independent phagosomal signals that contribute to phagocytosis and inflammation

https://doi.org/10.21203/rs.3.rs-279261/v1

This work is licensed under a CC BY 4.0 License

You are reading this latest preprint version

Background: 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).

Results: 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.

Conclusion: 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.

Allergy & Immune Disorders

Immunology

macrophage

phagocytosis

inflammation

borrelia

MyD88

Lyme disease (LD) is a highly prevalent tick-borne illness caused by the spirochetal bacterium Borrelia burgdorferi (Bb) (1–3). The disease is characterized by a wide array of clinical manifestations which vary in duration and severity between 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 the clinical findings characteristic of LD (6). 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–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 (10). 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 α₃ (13). 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 (16).

Bb is an extracellular pathogen that needs to be taken up and degraded for significant recognition by the host immune system (17). We have defined this process as “phagosomal signaling” (14). 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–19). Bb does not contain LPS and, therefore, does not engage TLR4. The cell envelope of Bb contains abundant triacylated lipoproteins (20), which are known to be recognized by TLR1/2 heterodimers (21–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 (27). 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–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 (36). 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 (38). In the context of Bb infection, lysosome maturation markers are recruited to Bb-containing phagosomes in macrophages lacking MyD88 (28). 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.

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.

Mice:

Female 6-8-week-old C57BL/6J wild type (WT) and C57BL/6J 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). Original WT breeding pairs were purchased from The Jackson Laboratory (Bar Harbor, Maine). 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 (39). 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 (40). 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.

Bacterial Strains:

Low-passage virulent wild-type Bb strain 297 (41) 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) (42) were maintained in Barbour-Stonner-Kelly (BSK)-II media supplemented with normal rabbit serum and gentamicin (50 µg/µl) (42). Cultures were grown at 23°C for at least one week prior to being shifted to 37°C as previously described (42). Spirochetes were centrifuged at 3300 x g for 20 minutes at 4°C and resuspended in either BSK-II for in vivo experiments or DMEM (Gibco, 15630-080) supplemented with sodium pyruvate (Gibco, 11360-070) and HEPES (Gibco, 15630-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 previous described (43).

BMDM Stimulation:

Bone marrow-derived macrophages (BMDMs) were isolated from 6-8-week-old WT and MyD88^−/− mice as described previously (18). Single cell macrophage suspensions were seeded into either 12-well tissue culture-treated plates at a concentration of 1 x 10⁶ cells/ml per well or 1 x 10⁵ cells/500 µL per well in 8-chamber cell microscopy slides. Slides and plates were then incubated overnight at 37 °C/5% CO₂ to allow cell adherence before experimentation. Cells were incubated for either 0.5, 1, 4 or 6 hours at 37 °C/5% CO₂ 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).

Confocal Microscopy:

After stimulation, BMDMs were fixed in 2% paraformaldehyde with 0.05% Triton-X-100 (Fisher, BP151-100) for 10 minutes. 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 was used to detect anti-MyD88, anti-TLR2, anti-ASC and anti-LAMP-1 (Life Technologies Alexa Fluor 350, A21093) (1:100). Incubations with primary and secondary antibodies were done for 1 hour 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 H₂O 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 hours) 1 hour 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 minutes. 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 minutes. 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 minutes. Protein pellets were resuspended in 2x Laemmli buffer. Lysates were boiled at 99°C for 10 minutes and run on a 12.5% SDS-PAGE gel at 140V for 1 hour (5µL per lane, 15 lanes). Proteins were then transferred to nitrocellulose membranes (Bio-Rad 162–0177) at 20V for 20 minutes. Membranes were blocked for 1 hour 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 minutes 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 hours at room temperature. Following additional washes, membranes were incubated in HyGlo spray chemilunescent substrates (Denville Scientific, E2400) for 5 minutes and imaged on a Biorad ChemiDoc MP imaging system.

Cytokine Analysis:

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 hours. Following stimulation, RNA was isolated using the total RNA isolation kit (Macherey-Nagel) and was used as input for the Ovation RNA-seq V1 kit (NuGen, San Carlos, CA). cDNA output was analyzed for correct size distribution with an Experion Standard Sensitivity RNA chip and quantified using a Qubit Fluorometer (Invitrogen, Carlsbad, CA). 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) (44) 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 (45) 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) (46). 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 (47), 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 (48) 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 (49). To extract the desired subnetworks, we used OCSANA (50) within the BiNOM plugin (51) 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 (7). Macrophages have also been shown as part of the inflammatory cell infiltrate in heart and joint tissue of mice experimentally infected with Bb (10, 52), and the importance of MyD88 in Bb clearance from mouse tissues has been previously reported (40, 53, 54). It has also been well established that MyD88 enhances phagocytosis of multiple bacterial species by macrophages (28, 32, 34, 35, 55). 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 hours. 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 1C). 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 (56–59). Importantly, we have demonstrated that in human monocytes TLR2 and TLR8 co-localize to endosomes containing Bb (19). In addition, other investigators have shown a prominent role for TLR7 in the Bb inflammatory response (60). Murine TLR8, unlike murine TLR7 and human TLR8, does not seem to utilize ssRNA as its ligand (61). 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 2B, 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 (17). This process, known as phagosome maturation, requires reduction of phagosome pH and fusion with lysosomes (62). 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-hour stimulation at MOI 10:1 showed that both WT and MyD88^−/− BMDMs contained degraded GFP + Bb within the cell actin matrix (Fig. 3A and 3B). 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 (63). Both WT and MyD88^−/− BMDMs showed comparable LAMP-1 and Bb colocalization in phagosomes (Fig. 3A and 3B, 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 hours 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 (64). However, in discordance with previous studies (65), 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 minutes or 6 hours (Figure S4). As previously reported (66) (Figure S4A and S4C), ASC activation was observed with Sa, but no ASC was observed in BMDMs stimulated with Bb at 30 minutes or 6 hours (Figure S4B and S4D). 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 hours. 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 4E), 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 (67). 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 (67). 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 (68) and the absence of ApoE increases Bb burdens in experimentally infected mice (68). We then used OCSANA, a specialized package available in Cytoscape (50) 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–19, 29). The findings from these studies, as well as others (28), 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. Of particular novelty, phagosomal signaling cascades induced by Bb ligands in macrophages trigger a number of inflammatory and chemotactic pathways. 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.

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. 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 (38). 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 minutes (55). 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–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 (36), 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 (36), 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. Rhoq, activated by Cdc42 codes for TC10, a protein involved in generating long filopodia protrusions (69). The gene Cyfip1, which encodes a part of the WAVE complex that regulates actin polymerization (70), was also up-regulated according to the network through Rac1 protein interactions. The WAVE complex has higher involvement with lamellipodia formations (71). 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, 72). More recently, the leukotriene LTB₄ has been shown to promote phagocytosis of Bb by macrophages (73), 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 (74). 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 (68), 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 (75). 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. (72). 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.

Lyme disease (LD), Borrelia burgdorferi (Bb), Toll-like receptor (TLR), pathogen associated molecular patterns (PAMPs), wild type (WT), UConn Health (UCH), green fluorescent protein (GFP), Barbour Stoner Kelly (BSK), Staphylococcus aureus (Sa), bone marrow derived macrophages (BMDM), bovine serum albumin (BSA), differentially expressed genes (DEG), Gene Ontology (GE), multiplicities of infection (MOI), Apoptosis-associated speck-like protein containing a CARD (ASC),

APPENDIX

Not Applicable

ACKNOWLEDGEMENTS

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.

COMPETING INTERESTS

The authors declared that they have no conflict of interest.

FUNDING

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.

AUTHORS CONTRIBUTIONS

All authors have read and approved the manuscript.

SB: Conceptualization, Formal Analysis, Investigation, Methodology, Validation, Visualization, Original Draft Preparation, Review and Editing

KLH: Methodology, Project Administration, Supervision, Review and Editing

PVL: Conceptualization, Formal Analysis, Project Administration, Supervision, Validation, Review and Editing

CJL: Investigation, Methodology, Review and Editing

JLC: Investigation

YR: Resources

JDR: Conceptualization, Project Administration, Review and Editing

JCS: Conceptualization, Funding Acquisition, Project Administration, Resources, Supervision, Original Draft Preparation, Review and Editing.

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.

Availability of Data and Material

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.

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.
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.
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.
Bockenstedt JJWaLK. Host Response. In: Radolf DSSaJD, editor. Borrelia: Molecular Biology, Host Interaction and Pathogenesis. Norfolk UK: 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.
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 research & therapy. 2010;12(5):R168.
Lasky CE, Olson RM, Brown CR. Macrophage Polarization during Murine Lyme Borreliosis. Infect Immun. 2015;83(7):2627-35.
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.
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. Proceedings of the National Academy of Sciences of the United States of America. 2012;109(4):1228-32.
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.
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.
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.
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.
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.
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.
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.
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. Proceedings of the National Academy of Sciences of the United States of America. 2011;108(9):3683-8.
Brandt ME, Riley BS, Radolf JD, Norgard MV. Immunogenic integral membrane proteins of Borrelia burgdorferi are lipoproteins. Infect Immun. 1990;58(4):983-91.
Schenk M, Belisle JT, Modlin RL. TLR2 looks at lipoproteins. Immunity. 2009;31(6):847-9.
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.
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.
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.
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. Nature medicine. 2002;8(8):878-84.
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.
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.
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.
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. Journal of leukocyte biology. 2013;94(6):1231-41.
Akira S, Takeda K, Kaisho T. Toll-like receptors: critical proteins linking innate and acquired immunity. Nat Immunol. 2001;2(8):675-80.
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.
Blander JM. Coupling Toll-like receptor signaling with phagocytosis: potentiation of antigen presentation. Trends Immunol. 2007;28(1):19-25.
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.
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.
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.
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.
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.
Blander JM, Medzhitov R. Regulation of phagosome maturation by signals from toll-like receptors. Science. 2004;304(5673):1014-8.
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.
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.
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.
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. The Journal of clinical investigation. 2009;119(12):3652-65.
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.
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.
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.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106.
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.
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.
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.
Vera-Licona P, Bonnet E, Barillot E, Zinovyev A. OCSANA: optimal combinations of interventions from network analysis. Bioinformatics. 2013;29(12):1571-3.
Zinovyev A, Viara E, Calzone L, Barillot E. BiNoM: a Cytoscape plugin for manipulating and analyzing biological networks. Bioinformatics. 2008;24(6):876-7.
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. The Journal of infectious diseases. 2012;206(2):283-91.
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.
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.
Yates RM, Russell DG. Phagosome maturation proceeds independently of stimulation of toll-like receptors 2 and 4. Immunity. 2005;23(4):409-17.
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.
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.
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.
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. Journal of leukocyte biology. 2008;84(1):280-91.
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.
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.
Blander JM, Medzhitov R. On regulation of phagosome maturation and antigen presentation. Nat Immunol. 2006;7(10):1029-35.
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.
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.
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. The Journal of infectious diseases. 2010;201(12):1849-58.
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.
Naj X, Linder S. Actin-Dependent Regulation of Borrelia burgdorferi Phagocytosis by Macrophages. Curr Top Microbiol Immunol. 2017;399:133-54.
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. Proceedings of the National Academy of Sciences of the United States of America. 2015;112(17):5491-6.
Murphy GA, Solski PA, Jillian SA, Perez de la Ossa P, 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.
Levin R, Grinstein S, Canton J. The life cycle of phagosomes: formation, maturation, and resolution. Immunol Rev. 2016;273(1):156-79.
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.
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.
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.
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).
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.

AuthorChecklistE10only.pdf
BetaactinCelllysates.tif
Caspase1Celllystates.tif
Caspase1Supernatant.tif
IL1betaCelllystates.tif
IL1betaSupernatant.tif
SupplementalFile1.xlsx
SupplementalFile2.xlsx
SupplementalFile3.xlsx
SupplementalFile4.xlsx
SupplementalFile5.xlsx
SupplementalFile6.xlsx
FigS1.png
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 hours (A). MyD88-/- BMDMs were stimulated at a MOI 100:1 for 6 hours (B). To calculate differential expression, reads from the stimulated samples were normalized to unstimulated BMDMs of the same genotype; WT (C) or MyD88-/- (D).
FigS2.png
Secondary Control for TLR2. Confocal 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.
FigS3.png
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 hours. N= 3 mouse BMDM per genotype *p-value<0.05, **p-value<0.01, ***p-value<0.001, NS=not significant
FigS4.png
Bb does not induce ASC formation in BMDMs. (A-D) Confocal images (40x) of BMDMs stimulated with either Sa (A and C) or Bb (B and D) for 30 minutes (A-B) or 6 hours (C-D). Blue is nucleus, red is ASC and green is the bacteria species.
FigS5.png
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.

PDF

Review #2 received at journal
15 Mar, 2021
Editorial decision: Minor revision
15 Mar, 2021
Review #1 received at journal
08 Mar, 2021
Reviewer #2 agreed at journal
28 Feb, 2021
Editor assigned by journal
23 Feb, 2021
Reviewers invited by journal
23 Feb, 2021
Reviewer #1 agreed at journal
23 Feb, 2021
Submission checks completed at journal
23 Feb, 2021
Editor invited by journal
23 Feb, 2021

You are reading this latest preprint version

Macrophage mediated recognition and clearance of Borrelia burgdorferi elicits MyD88-dependent and -independent phagosomal signals that contribute to phagocytosis and inflammation

Status:

Version 1

Abstract

Figures

Background

Methods

Mice:

Bacterial Strains:

BMDM Stimulation:

Confocal Microscopy:

Western Blotting of BMDM Supernatants and Lysates:

Cytokine Analysis:

Identification of Differentially Expressed Genes by RNA-Seq:

Identification of Enriched Transcription Binding Sites and Master Regulator Analysis:

Gene Ontology (GO) Enrichment Analysis:

Network Reconstruction and Network Analysis:

Results

Discussion

Conclusions

Abbreviations

Declarations

References

Supplementary Files

Status:

Version 1