Open access
Editor's Pick
Research Article
29 September 2020

A Distinct Microbiome Signature in Posttreatment Lyme Disease Patients

ABSTRACT

Lyme disease is the most common vector-borne disease in the United States, with an estimated incidence of 300,000 infections annually. Antibiotic intervention cures Lyme disease in the majority of cases; however, 10 to 20% of patients develop posttreatment Lyme disease syndrome (PTLDS), a debilitating condition characterized by chronic fatigue, pain, and cognitive difficulties. The underlying mechanism responsible for PTLDS symptoms, as well as a reliable diagnostic tool, has remained elusive. We reasoned that the gut microbiome may play an important role in PTLDS given that the symptoms overlap considerably with conditions in which a dysbiotic microbiome has been observed, including mood, cognition, and autoimmune disorders. Analysis of sequencing data from a rigorously curated cohort of patients with PTLDS revealed a gut microbiome signature distinct from that of healthy control subjects, as well as from that of intensive care unit (ICU) patients. Notably, microbiome sequencing data alone were indicative of PTLDS, which presents a potential, novel diagnostic tool for PTLDS.
IMPORTANCE Most patients with acute Lyme disease are cured with antibiotic intervention, but 10 to 20% endure debilitating symptoms such as fatigue, neurological complications, and myalgias after treatment, a condition known as posttreatment Lyme disease syndrome (PTLDS). The etiology of PTLDS is not understood, and objective diagnostic tools are lacking. PTLDS symptoms overlap several diseases in which patients exhibit alterations in their microbiome. We found that patients with PTLDS have a distinct microbiome signature, allowing for an accurate classification of over 80% of analyzed cases. The signature is characterized by an increase in Blautia, a decrease in Bacteroides, and other changes. Importantly, this signature supports the validity of PTLDS and is the first potential biological diagnostic tool for the disease.

INTRODUCTION

Lyme disease, a tick-borne infection caused by Borrelia burgdorferi, affects approximately 300,000 people annually in the United States (1, 2). The symptoms of acute Lyme disease are highly variable, and when untreated, it can progress in severity over time from malaise and flu-like symptoms to neurological disorders, cardiac complications, and, in late stages, arthritis (3). Antibiotic intervention typically cures Lyme disease; however, approximately 10 to 20% of Lyme patients develop posttreatment Lyme disease syndrome (PTLDS) with symptoms, including myalgias, chronic fatigue, and cognitive difficulties for more than 6 months after completion of antibiotic treatment (48). The etiopathology of PTLDS is unknown, but it presents with symptoms that overlap those of other diseases, including chronic fatigue syndrome (5, 9), depression, fibromyalgia (5), and multiple sclerosis (10).
Along with unknown etiopathology and a diverse range of symptoms, diagnosing PTLDS remains challenging. Although a clinical case definition proposed by the Infectious Diseases Society of America (IDSA) in 2006 has served as a specific research tool, there is no biological method to diagnose PTLDS (11). While clinical biomarkers associated with PTLDS have been observed, a suitable diagnostic method and therapy remain elusive. A positron emission tomography (PET) brain imaging study among patients with PTLDS demonstrated elevated microglial activation compared to that of controls, congruent with localized inflammation (12). Additional research has shown that a greater B. burgdorferi-specific plasmablast response prior to treatment favors a resolution of symptoms versus the development of PTLDS, which indicates that even before treatment, a patient’s immunological landscape plays an important role in the development of PTLDS (13). Compared to healthy controls, patients with PTLDS have significantly elevated expression of interferon alpha, greater antibody reactivity to brain antigens (14), increased levels of the chemokine CCL19 (15) and the cytokine interleukin 23 (IL-23) (16), and a decrease in the CD57 lymphocyte subset (17, 18). Furthermore, patients have a higher risk of developing new-onset autoimmune joint diseases after a Lyme erythema migrans rash (19). Therefore, while the etiopathology is still unknown, these markers indicate biological abnormalities among patients with PTLDS.
The microbiome has been implicated in many diseases with symptoms that overlap those of PTLDS, including autoimmune diseases. The pathogenesis of autoimmune disease is affected by environmental (20) and genetic (21) factors, as well as the gut microbiome (22, 23). The gut microbiome plays an important role in human health and has been shown to strongly influence host metabolism (24, 25) and the immune (22, 26) and nervous (27) systems, as well as provide crucial colonization resistance against a range of intestinal pathogens (28, 29). Further, microbiome compositional changes can alter immune tolerance (22, 23). For instance, members of the intestinal microbiota have been characterized as contributing to the development of the long-term sequelae of acute infection events upon disruption of tissue and immune homeostasis (30). Studies have found the microbiome to be on par with and often superior to the human genome in predicting disease states (31, 32). Indeed, many microbiome-wide association studies have established correlation, and sometimes causation, of the gut microbiome in diseases such as multiple sclerosis (33, 34), rheumatoid arthritis (35), and systemic lupus (36). Patients with PTLDS often undergo extensive antibiotic treatment (5) which likely causes adverse alterations to their microbiomes. A potential parallel exists in autoimmune disease, in which antibiotic use has been linked to an increase in disease frequency because of the dramatic impact antibiotics have on the microbiome (37).
Since PTLDS symptoms present similarly to diseases in which the microbiome is implicated, we reasoned that the same may be true for the gut microbiome of patients with PTLDS. We analyzed the gut microbiome of subjects with PTLDS from the John Hopkins Lyme Disease Research Center’s Study of Lyme Immunology and Clinical Endpoints (SLICE) cohorts, specifically drawing from a cross-sectional cohort of patients meeting the IDSA-proposed case definition for PTLDS. We performed 16s rRNA gene sequencing on stool samples from this cohort and report a gut microbiome signature associated with PTLDS. These data present a novel biomarker and potential diagnostic tool for PTLDS, as well as suggest a therapeutic avenue for PTLDS.

RESULTS

Curation of the PTLDS and control cohorts.

Fecal samples were collected from 87 patients with well-defined PTLDS in the SLICE cohort. Patients had medical record documentation of prior Lyme disease which met the CDC surveillance case definitions for definitive or probable Lyme disease and current nonspecific, patient-reported symptoms meeting an operationalized case definition for PTLDS with symptoms leading to functional impairment (8). Patients had all received appropriate antibiotic treatment at the time of their initial diagnosis of Lyme disease, and many had received subsequent antibiotics for treatment of persistent symptoms. The median time from Lyme disease symptom onset to the study visit was 1.1 years (interquartile range [IQR], 0.5 years to 3.3 years), and participants reported taking a median of 56 days (IQR, 30 days to 84 days) of antibiotics during that interval. Eight (9.2%) reported currently taking antibiotics at the time of the study visit. The mean age of this cohort sample was 48.3 years (standard deviation [SD], 14.7), and 36 (41.4%) of the subjects were female.
The healthy control cohort consisted of fecal samples collected from 17 healthy donors at Northeastern University, as well as 152 donors from a previously identified healthy subset of the American Gut (38). To control for the generally high levels of antibiotic use that could alter the microbiome in patients with PTLDS, a previously curated cohort of 123 samples of intensive care unit (ICU) patients from two time points (39) was also used as a control.
To analyze the gut microbiome composition, 16S rRNA gene sequencing was performed using the Illumina and Ion Torrent platforms; the dual-platform approach did not have a measurable effect on the data (Fig. 1). The sequences, combined with the control sequencing data, were processed using Qiita (40) and Quantitative Insights into Microbial Ecology 2 (QIIME2) (41). Closed-reference operational taxonomic units (OTUs), a common designation used instead of “species” or “genus,” were generated (97% identity) and analyzed. Closed-reference picking was performed because it allowed for increased sample size of the PTLDS cohort due to samples being processed in different platforms, but the conclusions of this study do not differ from an analysis using a subset of samples sequenced by Illumina technology and processed with Deblur to generate amplicon sequence variants (42) (see Fig. S1 in the supplemental material).
FIG 1
FIG 1 Principal-coordinate analysis of unweighted UniFrac distances of 16S rRNA genes. All samples from patients in the ICU, those with PTLDS, and the IT-Healthy and AGP Healthy cohorts are represented. The sequencing platform, Illumina, Ion Torrent, or both, used for each sample is indicated by shape.

Sample classification of PTLDS, ICU, and healthy fecal microbiomes.

We evaluated the ability of the fecal microbiome to distinguish PTLDS, ICU, and healthy cohorts using a supervised-learning random-forest classifier model to classify sample cohorts (43) (Fig. 2). Receiver operating characteristic (ROC) analysis was used to evaluate the accuracy of the model’s classifications (44). The model’s performance was quantified by reducing the two-dimensional ROC curve into a one-dimensional scalar value called the area under the ROC curve (AUROC). An AUROC is a value between 0 and 1, where 0.5 would lie along the diagonal line and indicate that the model was as effective at classifying samples as random chance. The higher the AUROC, the better the model is at differentiating the classes. Our model robustly distinguished the three cohorts with high accuracy, yielding rounded AUROC values of 1.00 (Fig. 3a), which indicates strong differences in the microbiomes of these cohorts. ICU samples versus healthy or PTLDS samples were correctly classified in 100% of samples, which is to be expected, given the heavy use of antibiotics by such patients that results in severe alteration of the microbiome (39). The model correctly classified 82.4% of PTLDS samples against ICU and healthy controls, with the remaining 17.6% of PTLDS samples being misclassified as healthy (Table 1). We report the relative abundance of the five most important features (OTUs) for sample classification (Fig. 3b). Of note, Blautia species (OTU identifiers [IDs] 4474380, 4465907, and 4327141) comprise three of the five most important features for classification and were observed at a significantly greater relative abundance in the PTLDS cohort (8.86% ± 1.26%) than in the ICU (0.070% ± 0.017%) or healthy (1.34% ± 0.18%) cohort (P value < 0.0001). Conversely, the two other top five features most important for classification were Staphylococcus aureus (OTU ID 446058), which was present at a significantly higher relative abundance in the ICU cohort (0.95% ± 0.56%) (P value < 0.0001) than in the PTLDS (0.0024% ± 0.00030%) or healthy (0.0077% ± 0.0020%) cohort, likely due to it being a widespread nosocomial pathogen (45), and a Roseburia species (OTU ID 4481427) elevated in the healthy cohort (0.29% ± 0.050%) compared to PTLDS (0.15% ± 0.045%) (not significant [NS]) or ICU (0.0024% ± 0.0013%) (P value < 0.0001) (Fig. 3b).
FIG 2
FIG 2 A simplified flow chart of the Qiime2 classifier model pipeline used to analyze the microbiome data.
FIG 3
FIG 3 (a) Receiver operating characteristic curve evaluating the ability of a random-forest classifier model to classify PTLDS, healthy, and ICU controls based on the fecal microbiome determined by 16S rRNA gene sequencing. Rounded area under the ROC curve values were 1.00 for all cohorts. Gray lines represent the null model or random chance. (b) Relative abundance plots of the 5 most important features (OTUs) for classification of PTLDS, healthy, and ICU controls based on the fecal microbiome. The first, third, and fifth ranked most important features were Blautia species (OTU IDs 4474380, 4465907, and 4327141); the relative abundance of the Blautia spp. were combined for clarity. S. aureus (OTU ID 446058) and Roseburia sp. (OTU ID 4481427) were the second and fourth most important features, respectively. Bars represent the mean relative abundance plus or minus the standard error of the mean. Statistical significance was determined using Kruskal-Wallis (nonparametric) test followed by Dunn’s multiple comparison. ****, P value < 0.0001. ns, not significant.
TABLE 1
TABLE 1 Probability to distinguish PTLDS from healthy and ICU controls based on the fecal microbiome compositiona
CohortValue for cohort
ICUPTLDSHealthy
ICU100
PTLDS00.8240.176
Healthy001
a
A random-forest classifier model was used to determine the ability to distinguish PTLDS from healthy and ICU cohorts based on the fecal microbiome. The reported number is the proportion of samples in each cohort classified into a given cohort.

Elucidation of the effect of antibiotics on patients with PTLDS.

Patients with PTLDS are treated with antibiotics, initially to eradicate B. burgdorferi and often subsequently after treatment failure, which might impact their gut microbiome composition. We therefore examined the possibility that antibiotics alone were responsible for the distinctive microbiomes observed in PTLDS. Out of the 87 patients in the PTLDS cohort, 79 (90.8%) were treated with antibiotics, for Lyme disease and other conditions, within 1 year of sample collection, 64 (73.6%) within 6 months, 36 (41.4%) within a month, and 23 (26.4%) within a week (Table 2). Doxycycline and amoxicillin, both commonly prescribed for the treatment of acute Lyme disease, were the most common antibiotics taken. Note that the patients misidentified as “healthy” in the classifier model (Table 1) received an antibiotic regimen that was no different from that of the bulk of PTLDS patients.
TABLE 2
TABLE 2 Summation of antibiotic use within the PTLDS cohorta
TimeNo. (%) for:
Antibiotic useDoxycyclineAmoxicillin
1 wk23 (26.4)6 (6.9)12 (13.8)
1 mo36 (41.4)12 (13.8)15 (17.2)
6 mo64 (73.6)35 (40.2)19 (21.8)
1 yr79 (90.8)46 (52.9)24 (27.6)
a
Time refers to the period within which antibiotics were taken prior to sample donation. Antibiotic use is the total number of patients (percent) who have used antibiotics within the time frame. Doxycycline and amoxicillin columns describe the number and percentage of PTLDS patients who have taken doxycycline and/or amoxicillin during the indicated time frame.
To investigate the role antibiotics may play in shaping the microbiome of patients with PTLDS, we used principal-coordinate analysis and identified the type of antibiotic used and the time since last antibiotic use relative to when the stool sample was collected. Importantly, patients with PTLDS did not cluster by time since antibiotic treatment (1 week, 1 month, 6 months, 1 year, or over 1 year) or by type of antibiotic (doxycycline, amoxicillin, both, other, or none) in principal-coordinates analysis (Fig. 4a). We then separated the PTLDS cohort into groups based on how recently a patient had taken antibiotics, within 1 week to 1 month or ≥6 months, and used a supervised random-forest classifier model to evaluate the ability of antibiotic history to distinguish these groups within the PTLDS cohort and healthy and ICU samples. The difference in antibiotic administration regimens did not distinguish patients within the PTLDS cohort (Fig. 4b). Therefore, while antibiotics likely affect the microbiome, our results suggest that antibiotic influence alone cannot explain the distinctive microbiomes of PTLDS patients.
FIG 4
FIG 4 (a) Principal-coordinate analysis of unweighted UniFrac distances of 16S rRNA genes. All samples from patients in the ICU (red) and IT-Healthy and AGP Healthy (blue) cohorts are represented. The time since antibiotic (Abx) treatment before sample collection is indicated by symbol color for PTLDS samples. The most recently taken type of antibiotic taken within 1 week to 1 year of sample collection is indicated by the symbol shape for PTLDS samples. (b) Area under the receiver operating characteristic curve (AUROC) evaluating the ability of a random-forest classifier model to classify the fecal microbiome in the healthy, ICU, and PTLDS cohorts, separated into two groups based on antibiotic use within the last 1 week to 1 month or equal to or greater than 6 months.

Grouping patients with PTLDS.

While the classifier model indicates a predictive signature of the fecal microbiome in PTLDS, the model misclassifies PTLDS for healthy samples 17.6% of the time (Table 1), suggesting that the microbiome may not be affected in all patients. Thus, we reasoned that subgrouping patients with PTLDS based on important taxonomic features, rather than averaging the entire cohort, would allow us to further identify important aspects of the microbiome in PTLDS. As Blautia was the most represented genus in the important features for classification (Fig. 3b), we used Blautia as a metric for the groupings. We observed that patients with a relative abundance of Blautia over 10% tended to have a decreased abundance of the genus Bacteroides, below 15%, compared to an average relative abundance of 23.15% in the healthy cohort (Fig. 5a). Given this, together with the known importance of Bacteroides as a common gut symbiont (46, 47) and the correlation between decreased Bacteroides in diseases with symptoms overlapping those of PTLDS, such as depression (48), we used Bacteroides as our secondary grouping metric. Plotting the relative abundance of Bacteroides versus the relative abundance of Blautia yielded three distinct subgroups in the PTLDS cohort, which we defined as follows: group 1 (G1), >10% Blautia and <15% Bacteroides; group 2 (G2), >15% Bacteroides; and group 3 (G3), <10% Blautia and <15% Bacteroides (Fig. 5b). Few samples in the healthy (1.78% of samples) and ICU (0.813%) cohorts overlapped with G1 (high Blautia and low Bacteroides), which comprised 30.92% of PTLDS samples, but greater overlap existed in G2 and G3 (Fig. 5b). In the classifier model (Fig. 3a), all samples in the test set that were defined as G1 were correctly classified as PTLDS.
FIG 5
FIG 5 (a) Relative abundances of Bacteroides in the fecal microbiomes of healthy and ICU control cohorts and patients with PTLDS. Bars indicate the mean plus or minus the standard error of the mean. Statistical significance was determined using the Kruskal-Wallis (nonparametric) test followed by Dunn’s multiple comparison (****, P value < 0.0001; ***, P value < 0.001; *, P value < 0.01). (b) Relative abundances of Blautia versus Bacteroides in the fecal microbiomes of healthy and ICU control cohorts and of patients with PTLDS separated into three groups, G1, G2, and G3. The groups were determined based on the relative abundances of Blautia and Bacteroides: G1, >10% Blautia and <15% Bacteroides; G2, >15% Bacteroides; and G3, <10% Blautia and <15% Bacteroides. (c) Relative abundances of Enterobacteriaceae in the ICU and PTLDS cohorts and the Ion Torrent subset of the healthy control cohort.
Expansion of proinflammatory Enterobacteriaceae is a common feature of disease-associated microbiomes (49), and we therefore examined the relative abundance of this group in patients with PTLDS. Although the family Enterobacteriaceae was not in the top 5 most important features for classification of the microbiome in PTLDS, healthy, and ICU cohorts, some patients with PTLDS had exceptionally high levels of Enterobacteriaceae compared to the healthy control population at Northeastern University (IT-Healthy). Of the 193 OTUs in the Enterobacteriaceae family represented in the combined data sets (PTLDS, ICU, and healthy) in this study, the mean relative abundance of Enterobacteriaceae in IT-Healthy was 1.14% (median = 0.0275%), compared to an average relative abundance of 9.20% in PTLDS subjects (median = 0.46%). Approximately one-fifth (19.5%) of patients with PTLDS presented with a relative abundance of Enterobacteriaceae over 10%. As expected, ICU patients had a higher average relative abundance of Enterobacteriaceae (31.21%) (Fig. 5c).
Microbiome-associated studies (MAS) have been found to be excellent predictors in various diseases (50), often outperforming genome-wide association studies (31), likely because the microbiome is a confluence of genetics and the environment. ROC analysis of the PTLDS cohort yielded a rounded AUROC of 1.00, correctly classifying patients with PTLDS for 82.4% of samples. These results are similar in accuracy to results for well-established microbiome-associated diseases (50) such as Clostridium difficile infection, while outperforming the predictive capabilities of other MAS, such as inflammatory bowel disease (IBD), in which abnormalities within the microbiome are strongly implicated (Fig. 6).
FIG 6
FIG 6 Ranked area under receiver operating characteristic curve (AUROC) reported by Duvallet et al. (50) for the classification of the fecal microbiome in each disease versus a healthy control cohort. ART, arthritis; ASD, autism spectrum disorder; CDI, Clostridium difficile infection; CRC, colorectal cancer; EDD, enteric diarrheal disease; HIV, human immunodeficiency virus; IBD, inflammatory bowel disease; LIV, liver disease; NASH, nonalcoholic steatohepatitis; nonCDI, non-Clostridium difficile infection; OB, obesity; PAR, Parkinson’s disease; T1D, type I diabetes.

DISCUSSION

A conservative estimate of the prevalence of PTLDS in the United States in 2020 is 800,000 (51). The etiopathology of PTLDS is unknown, and the disease lacks suitable laboratory diagnostics and therapeutic methods (5). We aimed to address this ambiguity by investigating the microbiome, which has been shown to play an important role in diseases with symptom overlap with PTLDS. Through 16S rRNA gene sequencing, we identified alterations in the gut microbiome in a cohort of well-characterized patients with PTLDS compared to two healthy control cohorts and an ICU control cohort. The majority of the ICU patients were on antibiotic treatment at the time of collection; this served as a control for antibiotic use among patients with PTLDS. Using a receiver operating characteristic analysis based on a random-forest classifier model, we found that the microbiome of PTLDS patients is distinct from those of ICU patients and healthy controls. Blautia species represented three of the five most important features for this classification, and the relative abundance was elevated in the PTLDS cohort compared to those in the healthy and ICU cohorts. Interestingly, an increased relative abundance of Blautia has been observed in several other diseases. In patients with type 1 diabetes, elevated Blautia abundance was observed and was correlated with increased IA-2 tyrosine phosphatase autoantibodies, important markers of autoimmunity (52). Increased Blautia has also been seen in obesity (53), Alzheimer’s disease (54), nonalcoholic fatty liver disease (55), and multiple sclerosis (56).
In addition, approximately one-fifth of patients with PTLDS had a relative abundance of Enterobacteriaceae over 10%, while the average relative abundance of Enterobacteriaceae in the healthy cohort collected at Northeastern University was 1.14%, consistent with reports of Enterobacteriaceae in the healthy gut microbiome (57, 58). Members of the Enterobacteriaceae family have a proinflammatory lipopolysaccharide (LPS) in the outer membrane and can exacerbate inflammation (59). A high relative abundance of Enterobacteriaceae is reported in inflammatory bowel diseases (60), in metabolic disorders like type 2 diabetes, and in immune diseases like cancers (49). Concomitantly, we report a depletion of Bacteroides in the G1 and G3 subsets of the PTLDS cohort which together comprised 64.4% of the PTLDS cohort. Bacteroides is a common member of the gut microbiome and plays important, symbiotic roles, such as modulation and regulation of the immune system, maintenance of intestinal integrity, and carbohydrate digestion (46, 47). We have previously reported that γ-aminobutyric acid (GABA) production by human-derived Bacteroides is widespread, and there is a correlation between brain signatures of depression and fecal Bacteroides levels in patients with major depressive disorder (48). Moreover, Bacteroides organisms are major producers of short-chain fatty acids, which have been shown to support host immune homeostasis both locally and systemically (61, 62).
While it is possible given the nature of these aberrations that the microbiome is causal or contributory to PTLDS, establishing this relationship is difficult, as no animal model of PTLDS exists. Our results suggest an intriguing opportunity to test causality by using fecal microbiota transplant (FMT) or defined symbiotic consortia to treat patients with PTLDS. FMT has been successfully used to treat Clostridium difficile infection in patients (6365). Furthermore, FMT or the administration of symbiotic bacteria has also been shown to be efficacious in treating multiple sclerosis (34), Parkinson’s disease (66), Alzheimer’s disease (67), and rheumatoid arthritis (68) in animal models of disease.
As well as suggesting a potential diagnostic tool through the interrogation of the fecal microbiome, the robustness of these results reinforces the validity of PTLDS by showing strong distinctions between the fecal microbiomes of a rigorously curated cohort of patients with PTLDS, patients in the ICU, and healthy controls. Previously reported biomarkers further validate PTLDS and provide an opportunity for the field to progress. These biomarkers include quantitative immune alterations; patients with PTLDS present with elevated levels of the T-cell chemokine CCL19 compared to patients with acute Lyme disease who have returned to health (15), an increase in the cytokine IL-23 (16), a decrease in the CD57 lymphocyte subset (17, 18), and a decreased plasmablast response prior to treatment (13). Furthermore, a pilot study used [11C]DPA-713 PET imaging to study cerebral glial activation and found that several brain regions had higher [11C]DPA-713 binding in patients with PTLDS than in healthy controls (12). In addition to these biomarkers, Fallon et al. (69) developed a survey, the General Symptom Questionnaire-30 (GSQ-30), to assess symptom burden and changes; patients with PTLDS reported higher GSQ-30 scores before treatment and maintained these scores until 6 months posttreatment. The GSQ-30 could be a powerful tool to accompany biomarkers like the gut microbiome in PTLDS. The existence of these biomarkers, along with the microbiome signature that we report, contributes to the evidence for a biological basis for PTLDS. It also supports clinical and accumulating research evidence, first published for treatment trials and population-based studies (7072), that persistent symptoms after treatment of Lyme disease are common and can significantly impact quality of life. The lack of sensitivity of PTLDS symptoms such as fatigue, pain, and cognitive dysfunction can lead to the conclusion that they are not different than the background noise levels in the general population (73). However, studies operationalizing the proposed case definition for PTLDS (11) which utilize standardized symptom and quality-of-life measures have shown that the prevalence and magnitude of these symptoms are often more severe (4, 74).
Antibiotic use likely affected the microbiome of patients with PTLDS; 9.2% of patients with PTLDS were on antibiotic treatment during the time of collection, and many had extensive antibiotic treatment in their recent health history. While antibiotics likely alter the microbiome composition in patients with PTLDS, our data show that the PTLDS microbiome does not cluster by antibiotic history, and having PTLDS is a better classifier than antibiotic history. Therefore, it is unlikely that antibiotic use alone explains the distinct PTLDS microbiome. Importantly, the microbiome of PTLDS patients was distinct from the microbiome of patients in the ICU. Regardless of the cause of the disruption observed in the microbiome in PTLDS, our data suggest that therapeutic intervention targeting the microbiome may ameliorate PTLDS symptoms. In conclusion, we report that a cohort of patients with PTLDS have microbiomes distinct from those of healthy and ICU controls. Furthermore, we show that through machine learning we can use the microbiome as a high-fidelity indicator of PTLDS. We reinforce the validity of this disease by showing strong distinctions between a rigorously curated cohort of patients with PTLDS and controls.

MATERIALS AND METHODS

Ethics statement.

The Northeastern University Institutional Review Board (IRB no. 16-02-22 and 08–11-16), Johns Hopkins Medicine Institutional Review Board (IRB no. 00035457), and the University of California San Diego Institutional Review Board (IRB no. 141853X) approved the collection of feces from human subjects. Informed written consent was obtained from all participants.

Study cohorts and sample collection.

(i) Posttreatment Lyme disease syndrome cohort. The posttreatment Lyme disease syndrome (PTLDS) cohort is part of the Study of Lyme disease Immunology and Clinical Events (SLICE) curated at the Johns Hopkins Lyme Disease Research Center. Detailed enrollment and eligibility criteria for this cohort have been previously described (4). Briefly, patients with PTLDS had medical record documentation of prior Lyme disease meeting the CDC surveillance case definitions with appropriate treatment and had current patient-reported symptoms of fatigue, cognitive dysfunction, and/or musculoskeletal pain resulting in functional impairment. Many of those enrolled had received subsequent antibiotics for treatment of their persistent symptoms, and participants were permitted to be actively taking antibiotics for their condition at the time of enrollment. Patients with PTLDS were also excluded for a range of preexisting or comorbid conditions with significant PTLDS symptom overlap and/or immunosuppressive effects. Information on appropriate antibiotic treatment for Lyme disease was abstracted from the medical record; subsequent antibiotic use was recorded as part of the research study visit.
Subjects were provided with stool collection containers containing 9 ml of 20% glycerol and BBL culture swabs (Becton, Dickinson and Company, Sparks, MD). From a single stool sample produced at any time of day, stool was self-collected into the collection container to reach ∼10 ml and swabs were taken; samples were returned to the Johns Hopkins Lyme Disease Research Center (MD) and stored at −80°C. Samples in stool collection containers were sequenced using Ion Torrent technology as described below, and swabs were sequenced using Illumina technology at the Knight Lab at University of California San Diego (UCSD). Additional metadata on prior treatment and current antibiotic use from participants with PTLDS were gathered as part of the larger clinical case series study.
(ii) Healthy control cohort. The healthy control cohort consisted of two healthy populations: a healthy cohort at Northeastern University (IT-Healthy; Boston, MA) and 152 donors from a healthy subset of the American Gut Project (38) (AGP Healthy). Sample processing for these cohorts was performed according to Earth Microbiome Project protocols (75). Using stool collection vessels (Medline Industries), one fresh stool sample was self-collected from 17 healthy adult donors. Donors were excluded if they were currently taking antibiotics or if they had taken antibiotics for at least 2 weeks at the time of collection. A sample of the stool was immediately placed in 9 ml of oxygen-prereduced phosphate-buffered saline (PBS) to a total of ∼10 ml of slurry in a 50-ml collection tube (Fisher Scientific). The stool slurry was quickly homogenized in a Coy anaerobic vinyl chamber (Coy Laboratory Products, Inc.) in 5% hydrogen, 10% CO2, and 85% nitrogen at 37°C. Samples were stored at −80°C and sequenced using Ion Torrent technology as described below. A healthy subset of the American Gut Project was identified as previously described (38); 152 samples were randomly selected from the healthy subset. Samples were collected and sequenced as previously reported (38).
(iii) ICU cohort. The ICU cohort consists of 123 samples from two time points (within 48 h of ICU admission and at ICU discharge or on ICU day 10) from critically ill patients in the intensive care unit in four centers across the United States and Canada as reported previously (39). Sample collection and processing for this cohort were performed according to Earth Microbiome Project protocols (75). This cohort served as a control for the effect of antibiotics on the microbiome, as the ICU patients had omnipresent antibiotic use. All ICU patients were treated with differing antibiotic regimens.

Preparation of DNA and 16S rRNA sequencing protocols.

(i) Ion Torrent sequencing. DNA extraction and sequencing were performed by MR DNA (Shallowater, TX) on an Ion Torrent PGM. The V4 variable region was amplified using PCR primers 515/806 (515F, GTGCCAGCMGCCGCGGTAA, and 806R, GGACTACVSGGGTATCTAAT) in a single-step 30 cycle PCR with the HotStarTaq Plus master mix kit (Qiagen, USA). The following conditions were used: 94°C for 3 min and 30 cycles of 94°C for 30 s, 53°C for 40 s, and 72°C for 1 min, followed by a final elongation step at 72°C for 5 min.
(ii) Illumina sequencing. Using the primers 515f/806rB, the V4 region was amplified and was sequenced as previously described (76) using an Illumina MiSeq (Knight Lab, UCSD). Sequencing data for the ICU cohort and the American Gut project were obtained in Qiita (40) (study IDs 2136 and 10317).

16S rRNA data analysis.

Raw sequencing data were uploaded and processed in Qiita (40) (study ID 11673); the sequences were demultiplexed and trimmed to 150 bp, and closed-reference OTUs were picked with Greengenes 13-8 (77) on an OTU similarity level of 97%. The OTU table was rarefied to 10,000 reads. Data were subsequently analyzed using the software package QIIME2 (41). Since Ion Torrent and Illumina sequencing followed the Earth Microbiome Project protocol, the sequencing platform did not have a measurable effect on the data (Fig. 1) and we were able to combine the sequencing platforms for analysis. To assess the ability of the PTLDS microbiome to be distinguished from healthy and ICU controls, the sample classifier tool in QIIME2 was used (43). A random-forest classifier was trained and evaluated. ROC curves were generated to summarize the true- versus false-positive rates; the area under the curve was calculated and reflects the ability of the classifier to distinguish between cohorts. The top five most important features for distinguishing the microbiomes were reported.

Data availability.

Data generated in this study are available in Qiita (https://qiita.ucsd.edu/study/description/11673) and the European Bioinformatics Institute (https://www.ebi.ac.uk/ena/browser/view/ERP122507).

ACKNOWLEDGMENTS

We thank Daniel McDonald, scientific director of The Microsetta Initiative at UC San Diego School of Medicine, for his helpful comments on the AGP and ICU data sets and reviewing the manuscript and Gail Ackermann for her help with sample metadata. We also thank Marcia Daniela Villega de Flores, Thelio Sewell, Erica Mihm, and Susan Joseph for their efforts in patient recruitment, consent, sample processing, and archiving.
This work was supported by the Steven and Alexandra Cohen Foundation.

Footnote

This article is a direct contribution from Kim Lewis, a Fellow of the American Academy of Microbiology, who arranged for and secured reviews by Janis Weis, University of Utah, and Olga Lomovskaya, Qpex Biopharma.

Supplemental Material

File (mbio.02310-20-sf001.tif)
ASM does not own the copyrights to Supplemental Material that may be linked to, or accessed through, an article. The authors have granted ASM a non-exclusive, world-wide license to publish the Supplemental Material files. Please contact the corresponding author directly for reuse.

REFERENCES

1.
Borchers AT, Keen CL, Huntley AC, Gershwin ME. 2015. Lyme disease: a rigorous review of diagnostic criteria and treatment. J Autoimmun 57:82–115.
2.
Hinckley AF, Connally NP, Meek JI, Johnson BJ, Kemperman MM, Feldman KA, White JL, Mead PS. 2014. Lyme disease testing by large commercial laboratories in the United States. Clin Infect Dis 59:676–681.
3.
Steere AC, Strle F, Wormser GP, Hu LT, Branda JA, Hovius JW, Li X, Mead PS. 2016. Lyme borreliosis. Nat Rev Dis Primers 2:16090.
4.
Rebman AW, Bechtold KT, Yang T, Mihm EA, Soloski MJ, Novak CB, Aucott JN. 2017. The clinical, symptom, and quality-of-life characterization of a well-defined group of patients with posttreatment Lyme disease syndrome. Front Med (Lausanne) 4:224.
5.
Rebman AW, Aucott JN. 2020. Post-treatment Lyme disease as a model for persistent symptoms in Lyme disease. Front Med (Lausanne) 7:57.
6.
Aucott JN. 2015. Posttreatment Lyme disease syndrome. Infect Dis Clin North Am 29:309–323.
7.
Marques A. 2008. Chronic Lyme disease: a review. Infect Dis Clin North Am 22:341–360. vii–viii.
8.
Aucott JN, Crowder LA, Kortte KB. 2013. Development of a foundation for a case definition of post-treatment Lyme disease syndrome. Int J Infect Dis 17:e443–e449.
9.
Gaudino EA, Coyle PK, Krupp LB. 1997. Post-Lyme syndrome and chronic fatigue syndrome. Neuropsychiatric similarities and differences. Arch Neurol 54:1372–1376.
10.
Fallon BA, Levin ES, Schweitzer PJ, Hardesty D. 2010. Inflammation and central nervous system Lyme disease. Neurobiol Dis 37:534–541.
11.
Wormser GP, Dattwyler RJ, Shapiro ED, Halperin JJ, Steere AC, Klempner MS, Krause PJ, Bakken JS, Strle F, Stanek G, Bockenstedt L, Fish D, Dumler JS, Nadelman RB. 2006. The clinical assessment, treatment, and prevention of Lyme disease, human granulocytic anaplasmosis, and babesiosis: clinical practice guidelines by the Infectious Diseases Society of America. Clin Infect Dis 43:1089–1134.
12.
Coughlin JM, Yang T, Rebman AW, Bechtold KT, Du Y, Mathews WB, Lesniak WG, Mihm EA, Frey SM, Marshall ES, Rosenthal HB, Reekie TA, Kassiou M, Dannals RF, Soloski MJ, Aucott JN, Pomper MG. 2018. Imaging glial activation in patients with post-treatment Lyme disease symptoms: a pilot study using [(11)C]DPA-713 PET. J Neuroinflammation 15:346.
13.
Blum LK, Adamska JZ, Martin DS, Rebman AW, Elliott SE, Cao RRL, Embers ME, Aucott JN, Soloski MJ, Robinson WH. 2018. Robust B cell responses predict rapid resolution of Lyme disease. Front Immunol 9:1634.
14.
Jacek E, Fallon BA, Chandra A, Crow MK, Wormser GP, Alaedini A. 2013. Increased IFNalpha activity and differential antibody response in patients with a history of Lyme disease and persistent cognitive deficits. J Neuroimmunol 255:85–91.
15.
Aucott JN, Soloski MJ, Rebman AW, Crowder LA, Lahey LJ, Wagner CA, Robinson WH, Bechtold KT. 2016. CCL19 as a chemokine risk factor for posttreatment Lyme disease syndrome: a prospective clinical cohort study. Clin Vaccine Immunol 23:757–766.
16.
Strle K, Stupica D, Drouin EE, Steere AC, Strle F. 2014. Elevated levels of IL-23 in a subset of patients with post-Lyme disease symptoms following erythema migrans. Clin Infect Dis 58:372–380.
17.
Stricker RB, Winger EE. 2001. Decreased CD57 lymphocyte subset in patients with chronic Lyme disease. Immunol Lett 76:43–48.
18.
Stricker RB, Burrascano J, Winger E. 2002. Longterm decrease in the CD57 lymphocyte subset in a patient with chronic Lyme disease. Ann Agric Environ Med 9:111–113.
19.
Arvikar SL, Crowley JT, Sulka KB, Steere AC. 2017. Autoimmune arthritides, rheumatoid arthritis, psoriatic arthritis, or peripheral spondyloarthritis following Lyme disease. Arthritis Rheumatol 69:194–202.
20.
Miller FW, Pollard KM, Parks CG, Germolec DR, Leung PS, Selmi C, Humble MC, Rose NR. 2012. Criteria for environmentally associated autoimmune diseases. J Autoimmun 39:253–258.
21.
Goodnow CC, Sprent J, Fazekas de St Groth B, Vinuesa CG. 2005. Cellular and genetic mechanisms of self tolerance and autoimmunity. Nature 435:590–597.
22.
Belkaid Y, Hand TW. 2014. Role of the microbiota in immunity and inflammation. Cell 157:121–141.
23.
Shamriz O, Mizrahi H, Werbner M, Shoenfeld Y, Avni O, Koren O. 2016. Microbiota at the crossroads of autoimmunity. Autoimmun Rev 15:859–869.
24.
Sanz Y, Olivares M, Moya-Perez A, Agostoni C. 2015. Understanding the role of gut microbiome in metabolic disease risk. Pediatr Res 77:236–244.
25.
Rowland I, Gibson G, Heinken A, Scott K, Swann J, Thiele I, Tuohy K. 2018. Gut microbiota functions: metabolism of nutrients and other food components. Eur J Nutr 57:1–24.
26.
Brestoff JR, Artis D. 2013. Commensal bacteria at the interface of host metabolism and the immune system. Nat Immunol 14:676–684.
27.
Fung TC, Olson CA, Hsiao EY. 2017. Interactions between the microbiota, immune and nervous systems in health and disease. Nat Neurosci 20:145–155.
28.
Buffie CG, Jarchum I, Equinda M, Lipuma L, Gobourne A, Viale A, Ubeda C, Xavier J, Pamer EG. 2012. Profound alterations of intestinal microbiota following a single dose of clindamycin results in sustained susceptibility to Clostridium difficile-induced colitis. Infect Immun 80:62–73.
29.
Abt MC, Osborne LC, Monticelli LA, Doering TA, Alenghat T, Sonnenberg GF, Paley MA, Antenus M, Williams KL, Erikson J, Wherry EJ, Artis D. 2012. Commensal bacteria calibrate the activation threshold of innate antiviral immunity. Immunity 37:158–170.
30.
Fonseca DM, Hand TW, Han SJ, Gerner MY, Glatman Zaretsky A, Byrd AL, Harrison OJ, Ortiz AM, Quinones M, Trinchieri G, Brenchley JM, Brodsky IE, Germain RN, Randolph GJ, Belkaid Y. 2015. Microbiota-dependent sequelae of acute infection compromise tissue-specific immunity. Cell 163:354–366.
31.
Ross EM, Moate PJ, Marett LC, Cocks BG, Hayes BJ. 2013. Metagenomic predictions: from microbiome to complex health and environmental phenotypes in humans and cattle. PLoS One 8:e73056.
32.
Gilbert JA, Quinn RA, Debelius J, Xu ZZ, Morton J, Garg N, Jansson JK, Dorrestein PC, Knight R. 2016. Microbiome-wide association studies link dynamic microbial consortia to disease. Nature 535:94–103.
33.
Probstel AK, Baranzini SE. 2018. The role of the gut microbiome in multiple sclerosis risk and progression: towards characterization of the “MS microbiome.” Neurotherapeutics 15:126–134.
34.
Berer K, Gerdes LA, Cekanaviciute E, Jia X, Xiao L, Xia Z, Liu C, Klotz L, Stauffer U, Baranzini SE, Kumpfel T, Hohlfeld R, Krishnamoorthy G, Wekerle H. 2017. Gut microbiota from multiple sclerosis patients enables spontaneous autoimmune encephalomyelitis in mice. Proc Natl Acad Sci U S A 114:10719–10724.
35.
Zhang X, Zhang D, Jia H, Feng Q, Wang D, Liang D, Wu X, Li J, Tang L, Li Y, Lan Z, Chen B, Li Y, Zhong H, Xie H, Jie Z, Chen W, Tang S, Xu X, Wang X, Cai X, Liu S, Xia Y, Li J, Qiao X, Al-Aama JY, Chen H, Wang L, Wu QJ, Zhang F, Zheng W, Li Y, Zhang M, Luo G, Xue W, Xiao L, Li J, Chen W, Xu X, Yin Y, Yang H, Wang J, Kristiansen K, Liu L, Li T, Huang Q, Li Y, Wang J. 2015. The oral and gut microbiomes are perturbed in rheumatoid arthritis and partly normalized after treatment. Nat Med 21:895–905.
36.
Lopez P, de Paz B, Rodriguez-Carrio J, Hevia A, Sanchez B, Margolles A, Suarez A. 2016. Th17 responses and natural IgM antibodies are related to gut microbiota composition in systemic lupus erythematosus patients. Sci Rep 6:24072.
37.
Rosser EC, Mauri C. 2016. A clinical update on the significance of the gut microbiota in systemic autoimmunity. J Autoimmun 74:85–93.
38.
McDonald D, Hyde E, Debelius JW, Morton JT, Gonzalez A, Ackermann G, Aksenov AA, Behsaz B, Brennan C, Chen Y, DeRight Goldasich L, Dorrestein PC, Dunn RR, Fahimipour AK, Gaffney J, Gilbert JA, Gogul G, Green JL, Hugenholtz P, Humphrey G, Huttenhower C, Jackson MA, Janssen S, Jeste DV, Jiang L, Kelley ST, Knights D, Kosciolek T, Ladau J, Leach J, Marotz C, Meleshko D, Melnik AV, Metcalf JL, Mohimani H, Montassier E, Navas-Molina J, Nguyen TT, Peddada S, Pevzner P, Pollard KS, Rahnavard G, Robbins-Pianka A, Sangwan N, Shorenstein J, Smarr L, Song SJ, Spector T, Swafford AD, Thackray VG, Thompson LR, Tripathi A, Vázquez-Baeza Y, Vrbanac A, Wischmeyer P, Wolfe E, Zhu Q, The American Gut Consortium, Knight R. 2018. American Gut: an open platform for citizen science microbiome research. mSystems 3:e00031-18.
39.
McDonald D, Ackermann G, Khailova L, Baird C, Heyland D, Kozar R, Lemieux M, Derenski K, King J, Vis-Kampen C, Knight R, Wischmeyer PE. 2016. Extreme dysbiosis of the microbiome in critical illness. mSphere 1:e00199-16.
40.
Gonzalez A, Navas-Molina JA, Kosciolek T, McDonald D, Vazquez-Baeza Y, Ackermann G, DeReus J, Janssen S, Swafford AD, Orchanian SB, Sanders JG, Shorenstein J, Holste H, Petrus S, Robbins-Pianka A, Brislawn CJ, Wang M, Rideout JR, Bolyen E, Dillon M, Caporaso JG, Dorrestein PC, Knight R. 2018. Qiita: rapid, web-enabled microbiome meta-analysis. Nat Methods 15:796–798.
41.
Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, Alexander H, Alm EJ, Arumugam M, Asnicar F, Bai Y, Bisanz JE, Bittinger K, Brejnrod A, Brislawn CJ, Brown CT, Callahan BJ, Caraballo-Rodriguez AM, Chase J, Cope EK, Da Silva R, Diener C, Dorrestein PC, Douglas GM, Durall DM, Duvallet C, Edwardson CF, Ernst M, Estaki M, Fouquier J, Gauglitz JM, Gibbons SM, Gibson DL, Gonzalez A, Gorlick K, Guo J, Hillmann B, Holmes S, Holste H, Huttenhower C, Huttley GA, Janssen S, Jarmusch AK, Jiang L, Kaehler BD, Kang KB, Keefe CR, Keim P, Kelley ST, Knights D, et al. 2019. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol 37:852–857.
42.
Amir A, McDonald D, Navas-Molina JA, Kopylova E, Morton JT, Zech Xu Z, Kightley EP, Thompson LR, Hyde ER, Gonzalez A, Knight R. 2017. Deblur rapidly resolves single-nucleotide community sequence patterns. mSystems 2:e00191-16.
43.
Bokulich NA, Dillon MR, Bolyen E, Kaehler BD, Huttley GA, Caporaso JG. 2018. q2-sample-classifier: machine-learning tools for microbiome classification and regression. J Open Res Softw 3:934.
44.
Fawcett T. 2006. An introduction to ROC analysis. Pattern Recognition Lett 27:861–874.
45.
Dancer SJ, Adams CE, Smith J, Pichon B, Kearns A, Morrison D. 2019. Tracking Staphylococcus aureus in the intensive care unit using whole-genome sequencing. J Hosp Infect 103:13–20.
46.
Tan H, Zhai Q, Chen W. 2019. Investigations of Bacteroides spp. towards next-generation probiotics. Food Res Int 116:637–644.
47.
Rios-Covian D, Salazar N, Gueimonde M, de Los Reyes-Gavilan CG. 2017. Shaping the metabolism of intestinal Bacteroides population through diet to improve human health. Front Microbiol 8:376.
48.
Strandwitz P, Kim KH, Terekhova D, Liu JK, Sharma A, Levering J, McDonald D, Dietrich D, Ramadhar TR, Lekbua A, Mroue N, Liston C, Stewart EJ, Dubin MJ, Zengler K, Knight R, Gilbert JA, Clardy J, Lewis K. 2019. GABA-modulating bacteria of the human gut microbiota. Nat Microbiol 4:396–403.
49.
Shin NR, Whon TW, Bae JW. 2015. Proteobacteria: microbial signature of dysbiosis in gut microbiota. Trends Biotechnol 33:496–503.
50.
Duvallet C, Gibbons SM, Gurry T, Irizarry RA, Alm EJ. 2017. Meta-analysis of gut microbiome studies identifies disease-specific and shared responses. Nat Commun 8:1784.
51.
DeLong A, Hsu M, Kotsoris H. 2019. Estimation of cumulative number of post-treatment Lyme disease cases in the US, 2016 and 2020. BMC Public Health 19:352.
52.
Qi CJ, Zhang Q, Yu M, Xu JP, Zheng J, Wang T, Xiao XH. 2016. Imbalance of fecal microbiota at newly diagnosed type 1 diabetes in Chinese children. Chin Med J (Engl) 129:1298–1304.
53.
Hall AB, Tolonen AC, Xavier RJ. 2017. Human genetic variation and the gut microbiome in disease. Nat Rev Genet 18:690–699.
54.
Vogt NM, Kerby RL, Dill-McFarland KA, Harding SJ, Merluzzi AP, Johnson SC, Carlsson CM, Asthana S, Zetterberg H, Blennow K, Bendlin BB, Rey FE. 2017. Gut microbiome alterations in Alzheimer’s disease. Sci Rep 7:13537.
55.
Shen F, Zheng RD, Sun XQ, Ding WJ, Wang XY, Fan JG. 2017. Gut microbiota dysbiosis in patients with non-alcoholic fatty liver disease. Hepatobiliary Pancreat Dis Int 16:375–381.
56.
Chen J, Chia N, Kalari KR, Yao JZ, Novotna M, Paz Soldan MM, Luckey DH, Marietta EV, Jeraldo PR, Chen X, Weinshenker BG, Rodriguez M, Kantarci OH, Nelson H, Murray JA, Mangalam AK. 2016. Multiple sclerosis patients have a distinct gut microbiota compared to healthy controls. Sci Rep 6:28484.
57.
Human Microbiome Project Consortium. 2012. Structure, function and diversity of the healthy human microbiome. Nature 486:207–214.
58.
Winter SE, Lopez CA, Baumler AJ. 2013. The dynamics of gut-associated microbial communities during inflammation. EMBO Rep 14:319–327.
59.
Park BS, Song DH, Kim HM, Choi BS, Lee H, Lee JO. 2009. The structural basis of lipopolysaccharide recognition by the TLR4-MD-2 complex. Nature 458:1191–1195.
60.
Zeng MY, Inohara N, Nunez G. 2017. Mechanisms of inflammation-driven bacterial dysbiosis in the gut. Mucosal Immunol 10:18–26.
61.
Arpaia N, Campbell C, Fan X, Dikiy S, van der Veeken J, deRoos P, Liu H, Cross JR, Pfeffer K, Coffer PJ, Rudensky AY. 2013. Metabolites produced by commensal bacteria promote peripheral regulatory T-cell generation. Nature 504:451–455.
62.
Furusawa Y, Obata Y, Fukuda S, Endo TA, Nakato G, Takahashi D, Nakanishi Y, Uetake C, Kato K, Kato T, Takahashi M, Fukuda NN, Murakami S, Miyauchi E, Hino S, Atarashi K, Onawa S, Fujimura Y, Lockett T, Clarke JM, Topping DL, Tomita M, Hori S, Ohara O, Morita T, Koseki H, Kikuchi J, Honda K, Hase K, Ohno H. 2013. Commensal microbe-derived butyrate induces the differentiation of colonic regulatory T cells. Nature 504:446–450.
63.
Kelly CR, Khoruts A, Staley C, Sadowsky MJ, Abd M, Alani M, Bakow B, Curran P, McKenney J, Tisch A, Reinert SE, Machan JT, Brandt LJ. 2016. Effect of fecal microbiota transplantation on recurrence in multiply recurrent Clostridium difficile infection: a randomized trial. Ann Intern Med 165:609–616.
64.
Brumbaugh DE, De Zoeten EF, Pyo-Twist A, Fidanza S, Hughes S, Dolan SA, Child J, Dominguez SR. 2018. An intragastric fecal microbiota transplantation program for treatment of recurrent Clostridium difficile in children is efficacious, safe, and inexpensive. J Pediatr 194:123–127.e1.
65.
Nicholson MR, Mitchell PD, Alexander E, Ballal S, Bartlett M, Becker P, Davidovics Z, Docktor M, Dole M, Felix G, Gisser J, Hourigan SK, Jensen MK, Kaplan JL, Kelsen J, Kennedy M, Khanna S, Knackstedt E, Leier M, Lewis J, Lodarek A, Michail S, Oliva-Hemker M, Patton T, Queliza K, Russell GH, Singh N, Solomon A, Suskind DL, Werlin S, Kellermayer R, Kahn SA. 2020. Efficacy of fecal microbiota transplantation for Clostridium difficile infection in children. Clin Gastroenterol Hepatol 18:612–619.e1.
66.
Sun M-F, Zhu Y-L, Zhou Z-L, Jia X-B, Xu Y-D, Yang Q, Cui C, Shen Y-Q. 2018. Neuroprotective effects of fecal microbiota transplantation on MPTP-induced Parkinson’s disease mice: gut microbiota, glial reaction and TLR4/TNF-α signaling pathway. Brain Behav Immun 70:48–60.
67.
Sun J, Xu J, Ling Y, Wang F, Gong T, Yang C, Ye S, Ye K, Wei D, Song Z, Chen D, Liu J. 2019. Fecal microbiota transplantation alleviated Alzheimer’s disease-like pathogenesis in APP/PS1 transgenic mice. Transl Psychiatry 9:189.
68.
Amdekar S, Singh V, Singh R, Sharma P, Keshav P, Kumar A. 2011. Lactobacillus casei reduces the inflammatory joint damage associated with collagen-induced arthritis (CIA) by reducing the pro-inflammatory cytokines: Lactobacillus casei: COX-2 inhibitor. J Clin Immunol 31:147–154.
69.
Fallon BA, Zubcevik N, Bennett C, Doshi S, Rebman AW, Kishon R, Moeller JR, Octavien NR, Aucott JN. 2019. The General Symptom Questionnaire-30 (GSQ-30): a brief measure of multi-system symptom burden in Lyme disease. Front Med (Lausanne) 6:283.
70.
Shadick NA, Phillips CB, Logigian EL, Steere AC, Kaplan RF, Berardi VP, Duray PH, Larson MG, Wright EA, Ginsburg KS, Katz JN, Liang MH. 1994. The long-term clinical outcomes of Lyme disease. A population-based retrospective cohort study. Ann Intern Med 121:560–567.
71.
Asch ES, Bujak DI, Weiss M, Peterson MG, Weinstein A. 1994. Lyme disease: an infectious and postinfectious syndrome. J Rheumatol 21:454–461.
72.
Steere AC, Hutchinson GJ, Rahn DW, Sigal LH, Craft JE, DeSanna ET, Malawista SE. 1983. Treatment of the early manifestations of Lyme disease. Ann Intern Med 99:22–26.
73.
Seltzer EG, Gerber MA, Cartter ML, Freudigman K, Shapiro ED. 2000. Long-term outcomes of persons with Lyme disease. JAMA 283:609–616.
74.
Aucott JN, Rebman AW, Crowder LA, Kortte KB. 2013. Post-treatment Lyme disease syndrome symptomatology and the impact on life functioning: is there something here? Qual Life Res 22:75–84.
75.
Gilbert JA, Jansson JK, Knight R. 2014. The Earth Microbiome Project: successes and aspirations. BMC Biol 12:69.
76.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, Owens SM, Betley J, Fraser L, Bauer M, Gormley N, Gilbert JA, Smith G, Knight R. 2012. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J 6:1621–1624.
77.
McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, Andersen GL, Knight R, Hugenholtz P. 2012. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J 6:610–618.

Information & Contributors

Information

Published In

cover image mBio
mBio
Volume 11Number 527 October 2020
eLocator: 10.1128/mbio.02310-20
Editor: Michael S. Gilmore, Harvard Medical School
PubMed: 32994327

History

Received: 17 August 2020
Accepted: 20 August 2020
Published online: 29 September 2020

Keywords

  1. Lyme disease
  2. diagnostics
  3. microbial communities
  4. microflora
  5. tick-borne pathogens

Contributors

Authors

Madeleine Morrissette
Antimicrobial Discovery Center, Department of Biology, Northeastern University, Boston, Massachusetts, USA
Norman Pitt
Antimicrobial Discovery Center, Department of Biology, Northeastern University, Boston, Massachusetts, USA
Antonio González
Department of Pediatrics, University of California San Diego, La Jolla, California, USA
Philip Strandwitz
Antimicrobial Discovery Center, Department of Biology, Northeastern University, Boston, Massachusetts, USA
Present address: Philip Strandwitz, Holobiome, Cambridge, Massachusetts, USA; Mariaelena Caboni, Holobiome, Cambridge, Massachusetts, USA.
Mariaelena Caboni
Antimicrobial Discovery Center, Department of Biology, Northeastern University, Boston, Massachusetts, USA
Present address: Philip Strandwitz, Holobiome, Cambridge, Massachusetts, USA; Mariaelena Caboni, Holobiome, Cambridge, Massachusetts, USA.
Alison W. Rebman
Lyme Disease Research Center, Division of Rheumatology, Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
Rob Knight
Department of Pediatrics, University of California San Diego, La Jolla, California, USA
Department of Computer Science and Engineering, University of California, San Diego, La Jolla, California, USA
Center for Microbiome Innovation, University of California, San Diego, La Jolla, California, USA
Anthony D’Onofrio
Antimicrobial Discovery Center, Department of Biology, Northeastern University, Boston, Massachusetts, USA
Lyme Disease Research Center, Division of Rheumatology, Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
Mark J. Soloski
Lyme Disease Research Center, Division of Rheumatology, Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
Antimicrobial Discovery Center, Department of Biology, Northeastern University, Boston, Massachusetts, USA

Editor

Michael S. Gilmore
Editor
Harvard Medical School

Notes

Address correspondence to Kim Lewis, [email protected].
Madeleine Morrissette and Norman Pitt contributed equally to this work. Author order was determined by seniority.

Metrics & Citations

Metrics

Note:

  • For recently published articles, the TOTAL download count will appear as zero until a new month starts.
  • There is a 3- to 4-day delay in article usage, so article usage will not appear immediately after publication.
  • Citation counts come from the Crossref Cited by service.

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. For an editable text file, please select Medlars format which will download as a .txt file. Simply select your manager software from the list below and click Download.

View Options

Figures

Tables

Media

Share

Share

Share the article link

Share with email

Email a colleague

Share on social media

American Society for Microbiology ("ASM") is committed to maintaining your confidence and trust with respect to the information we collect from you on websites owned and operated by ASM ("ASM Web Sites") and other sources. This Privacy Policy sets forth the information we collect about you, how we use this information and the choices you have about how we use such information.
FIND OUT MORE about the privacy policy