A chemosensory-like histidine kinase is dispensable for chemotaxis in vitro but regulates the virulence of Borrelia burgdorferi through modulating the stability of RpoS.

As an enzootic pathogen, the Lyme disease bacterium Borrelia burgdorferi possesses multiple copies of chemotaxis proteins, including two chemotaxis histidine kinases (CHK), CheA1 and CheA2. Our previous study showed that CheA2 is a genuine CHK that is required for chemotaxis; however, the role of CheA1 remains mysterious. This report first compares the structural features that differentiate CheA1 and CheA2 and then provides evidence to show that CheA1 is an atypical CHK that controls the virulence of B. burgdorferi through modulating the stability of RpoS, a key transcriptional regulator of the spirochete. First, microscopic analyses using green-fluorescence-protein (GFP) tags reveal that CheA1 has a unique and dynamic cellular localization. Second, loss-of-function studies indicate that CheA1 is not required for chemotaxis in vitro despite sharing a high sequence and structural similarity to its counterparts from other bacteria. Third, mouse infection studies using needle inoculations show that a deletion mutant of CheA1 (cheA1mut) is able to establish systemic infection in immune-deficient mice but fails to do so in immune-competent mice albeit the mutant can survive at the inoculation site for up to 28 days. Tick and mouse infection studies further demonstrate that CheA1 is dispensable for tick colonization and acquisition but essential for tick transmission. Lastly, mechanistic studies combining immunoblotting, protein turnover, mutagenesis, and RNA-seq analyses reveal that depletion of CheA1 affects RpoS stability, leading to reduced expression of several RpoS-regulated virulence factors (i.e., OspC, BBK32, and DbpA), likely due to dysregulated clpX and lon protease expression. Bulk RNA-seq analysis of infected mouse skin tissues further show that cheA1mut fails to elicit mouse tnf-α, il-10, il-1ß, and ccl2 expression, four important cytokines for Lyme disease development and B. burgdorferi transmigration. Collectively, these results reveal a unique role and regulatory mechanism of CheA1 in modulating virulence factor expression and add new insights into understanding the regulatory network of B. burgdorferi.

Whole genome sequencing of human Borrelia burgdorferi isolates reveals linked blocks of accessory genome elements located on plasmids and associated with human dissemination.

Lyme disease is the most common vector-borne disease in North America and Europe. The clinical manifestations of Lyme disease vary based on the genospecies of the infecting Borrelia burgdorferi spirochete, but the microbial genetic elements underlying these associations are not known. Here, we report the whole genome sequence (WGS) and analysis of 299 B. burgdorferi (Bb) isolates derived from patients in the Eastern and Midwestern US and Central Europe. We develop a WGS-based classification of Bb isolates, confirm and extend the findings of previous single- and multi-locus typing systems, define the plasmid profiles of human-infectious Bb isolates, annotate the core and strain-variable surface lipoproteome, and identify loci associated with disseminated infection. A core genome consisting of ~900 open reading frames and a core set of plasmids consisting of lp17, lp25, lp36, lp28-3, lp28-4, lp54, and cp26 are found in nearly all isolates. Strain-variable (accessory) plasmids and genes correlate strongly with phylogeny. Using genetic association study methods, we identify an accessory genome signature associated with dissemination in humans and define the individual plasmids and genes that make up this signature. Strains within the RST1/WGS A subgroup, particularly a subset marked by the OspC type A genotype, have increased rates of dissemination in humans. OspC type A strains possess a unique set of strongly linked genetic elements including the presence of lp56 and lp28-1 plasmids and a cluster of genes that may contribute to their enhanced virulence compared to other genotypes. These features of OspC type A strains reflect a broader paradigm across Bb isolates, in which near-clonal genotypes are defined by strain-specific clusters of linked genetic elements, particularly those encoding surface-exposed lipoproteins. These clusters of genes are maintained by strain-specific patterns of plasmid occupancy and are associated with the probability of invasive infection.

Borrelia burgdorferi Outer Surface Protein C Is Not the Sole Determinant of Dissemination in Mammals.

Lyme disease in the United States is most often caused by Borrelia burgdorferi sensu stricto. After a tick bite, the patient may develop erythema migrans at that site. If hematogenous dissemination occurs, the patient may then develop neurologic manifestations, carditis, or arthritis. Host-pathogen interactions include factors that contribute to hematogenous dissemination to other body sites. Outer surface protein C (OspC), a surface-exposed lipoprotein of B. burgdorferi, is essential during the early stages of mammalian infection. There is a high degree of genetic variation at the ospC locus, and certain ospC types are more frequently associated with hematogenous dissemination in patients, suggesting that OspC may be a major contributing factor to the clinical outcome of B. burgdorferi infection. In order to evaluate the role of OspC in B. burgdorferi dissemination, ospC was exchanged between B. burgdorferi isolates with different capacities to disseminate in laboratory mice, and these strains were then tested for their ability to disseminate in mice. The results indicated that the ability of B. burgdorferi to disseminate in mammalian hosts does not depend on OspC alone. The complete genome sequences of two closely related strains of B. burgdorferi with differing dissemination phenotypes were determined, but a specific genetic locus that could explain the differences in the phenotypes could not be definitively identified. The animal studies performed clearly demonstrated that OspC is not the sole determinant of dissemination. Future studies of the type described here with additional borrelial strains will hopefully clarify the genetic elements associated with hematogenous dissemination.

Stability of Kuramoto networks subject to large and small fluctuations from heterogeneous and spatially correlated noise.

Oscillatory networks subjected to noise are broadly used to model physical and technological systems. Due to their nonlinear coupling, such networks typically have multiple stable and unstable states that a network might visit due to noise. In this article, we focus on the assessment of fluctuations resulting from heterogeneous and spatially correlated noise inputs on Kuramoto model networks. We evaluate the typical, small fluctuations near synchronized states and connect the network variance to the overlap between stable modes of synchronization and the input noise covariance. Going beyond small to large fluctuations, we introduce the indicator mode approximation that projects the dynamics onto a single amplitude dimension. Such an approximation allows for estimating rates of fluctuations to saddle instabilities, resulting in phase slips between connected oscillators. Statistics for both regimes are quantified in terms of effective noise amplitudes that are compared and contrasted for several noise models. Bridging the gap between small and large fluctuations, we show that a larger network variance does not necessarily lead to higher rates of large fluctuations.

The chaotic milling behaviors of interacting swarms after collision.

We consider the problem of characterizing the dynamics of interacting swarms after they collide and form a stationary center of mass. Modeling efforts have shown that the collision of near head-on interacting swarms can produce a variety of post-collision dynamics including coherent milling, coherent flocking, and scattering behaviors. In particular, recent analysis of the transient dynamics of two colliding swarms has revealed the existence of a critical transition whereby the collision results in a combined milling state about a stationary center of mass. In the present work, we show that the collision dynamics of two swarms that form a milling state transitions from periodic to chaotic motion as a function of the repulsive force strength and its length scale. We used two existing methods as well as one new technique: Karhunen-Loeve decomposition to show the effective modal dimension chaos lives in, the 0-1 test to identify chaos, and then constrained correlation embedding to show how each swarm is embedded in the other when both swarms combine to form a single milling state after collision. We expect our analysis to impact new swarm experiments which examine the interaction of multiple swarms.

Lack of Convincing Evidence That Borrelia burgdorferi Infection Causes Either Alzheimer Disease or Lewy Body Dementia.

The role that microorganisms might have in the development of Alzheimer disease is a topic of considerable interest. In this article, we discuss whether there is credible evidence that Lyme disease is a cause of Alzheimer disease and critically review a recent publication that claimed that Borrelia burgdorferi sensu stricto infection, the primary cause of Lyme disease in the United States, may cause Lewy body dementia. We conclude that no convincing evidence exists that Lyme disease is a cause of either Alzheimer disease or Lewy body dementia.

Outbreak Size Distribution in Stochastic Epidemic Models.

Motivated by recent epidemic outbreaks, including those of COVID-19, we solve the canonical problem of calculating the dynamics and likelihood of extensive outbreaks in a population within a large class of stochastic epidemic models with demographic noise, including the susceptible-infected-recovered (SIR) model and its general extensions. In the limit of large populations, we compute the probability distribution for all extensive outbreaks, including those that entail unusually large or small (extreme) proportions of the population infected. Our approach reveals that, unlike other well-known examples of rare events occurring in discrete-state stochastic systems, the statistics of extreme outbreaks emanate from a full continuum of Hamiltonian paths, each satisfying unique boundary conditions with a conserved probability flux.

Multipartite Genome of Lyme Disease Borrelia: Structure, Variation and Prophages.

All members of the Borrelia genus that have been examined harbour a linear chromosome that is about 900 kbp in length, as well as a plethora of both linear and circular plasmids in the 5-220 kbp size range. Genome sequences for 27 Lyme disease Borrelia isolates have been determined since the elucidation of the B. burgdorferi B31 genome sequence in 1997. The chromosomes, which carry the vast majority of the housekeeping genes, appear to be very constant in gene content and organization across all Lyme disease Borrelia species. The content of the plasmids, which carry most of the genes that encode the differentially expressed surface proteins that interact with the spirochete's arthropod and vertebrate hosts, is much more variable. Lyme disease Borrelia isolates carry between 7-21 different plasmids, ranging in size from 5-84 kbp. All strains analyzed to date harbor three plasmids, cp26, lp54 and lp17. The plasmids are unusual, as compared to most bacterial plasmids, in that they contain many paralogous sequences, a large number of pseudogenes, and, in some cases, essential genes. In addition, a number of the plasmids have features indicating that they are prophages. Numerous methods have been developed for Lyme disease Borrelia strain typing. These have proven valuable for clinical and epidemiological studies, as well as phylogenomic and population genetic analyses. Increasingly, these approaches have been displaced by whole genome sequencing techniques. Some correlations between genome content and pathogenicity have been deduced, and comparative whole genome analyses promise future progress in this arena.

The Lyme Disease Biobank: Characterization of 550 Patient and Control Samples from the East Coast and Upper Midwest of the United States.

Lyme disease (LD) is an increasing public health problem. Current laboratory testing is insensitive in early infection, the stage at which appropriate treatment is most effective in preventing disease sequelae. The Lyme Disease Biobank (LDB) collects samples from individuals with symptoms consistent with early LD presenting with or without erythema migrans (EM) or an annular, expanding skin lesion and uninfected individuals from areas of endemicity. Samples were collected from 550 participants (298 cases and 252 controls) according to institutional review board-approved protocols and shipped to a centralized biorepository. Testing was performed to confirm the presence of tick-borne pathogens by real-time PCR, and a subset of samples was tested for Borrelia burgdorferi by culture. Serology was performed on all samples using the CDC's standard two-tiered testing algorithm (STTTA) for LD. LD diagnosis was supported by laboratory testing in 82 cases, including positive results by use of the STTTA, PCR, or culture or positive results by two enzyme-linked immunosorbent assays for cases presenting with EM lesion sizes of >5 cm. The remaining 216 cases had negative laboratory testing results. For the controls, 43 were positive by at least one of the tiers and 6 were positive by use of the STTTA. The results obtained with this collection highlight and reinforce the known limitations of serologic testing in early LD, with only 29% of individuals presenting with EM lesion sizes of >5 cm yielding a positive result using the STTTA. Aliquots of whole blood, serum, and urine from clinically characterized patients with and without LD are available to investigators in academia and industry for evaluation or development of novel diagnostic assays for LD, to continue to improve upon currently available methods.

Rejection of the name Borreliella and all proposed species comb. nov. placed therein.

Rejection (nomen rejiciendum) of the name Borreliella and all new combinations therein is being requested on grounds of risk to human health and patient safety (Principle 1, subprinciple 2 and Rule 56a) and violation to aim for stability of names, to avoid useless creation of names (Principle 1, subprinciple 1 and 3) and that names should not be changed without sufficient reason (Principle 9 of the International Code of Nomenclature of Prokaryotes).

Delay induced swarm pattern bifurcations in mixed reality experiments.

Swarms of coupled mobile agents subject to inter-agent wireless communication delays are known to exhibit multiple dynamic patterns in space that depend on the strength of the interactions and the magnitude of the communication delays. We experimentally demonstrate communication delay-induced bifurcations in the spatiotemporal patterns of robot swarms using two distinct hardware platforms in a mixed reality framework. Additionally, we make steps toward experimentally validating theoretically predicted parameter regions where transitions between swarm patterns occur. We show that multiple rotation patterns persist even when collision avoidance strategies are incorporated, and we show the existence of multi-stable, co-existing rotational patterns not predicted by usual mean field dynamics. Our experiments are the first significant steps toward validating existing theory and the existence and robustness of the delay-induced patterns in real robotic swarms.

Torus bifurcations of large-scale swarms having range dependent communication delay.

Dynamical emergent patterns of swarms are now fairly well established in nature and include flocking and rotational states. Recently, there has been great interest in engineering and physics to create artificial self-propelled agents that communicate over a network and operate with simple rules, with the goal of creating emergent self-organizing swarm patterns. In this paper, we show that when communicating networks have range dependent delays, rotational states, which are typically periodic, undergo a bifurcation and create swarm dynamics on a torus. The observed bifurcation yields additional frequencies into the dynamics, which may lead to quasi-periodic behavior of the swarm.

Bacterial Content of Fortified and Unfortified Holder Pasteurized Donor Human Milk During Prolonged Refrigerated Storage.

OBJECTIVE: To assess the microbial growth in unfortified and fortified Holder pasteurized donor human milk (HPDHM) during 96 hours of refrigerated storage in a clinical setting. METHODS: Thirty-six unfortified samples and 77 fortified samples of HPDHM were prepared in a neonatal intensive care milk preparation room and stored in the NICU refrigerator at 4°C to simulate a real-life feeding environment. One milliliter aliquots were removed at 24-hour intervals and cultured in duplicate for bacterial growth on solid blood agar medium. Viable bacterial colonies were characterized using standard microbiological methods. RESULTS: 96.5% of milk samples manipulated in a vertical laminar flow hood were negative for bacterial growth. In the remainder 3.5% of the samples, the maximum growth was 1 colony forming unit/0.1 ml plated. Higher colony counts were observed when the laminar hood was not used. In all cases, the colonies represented common skin bacteria and demonstrated an inconsistent and unsustained growth. Fortifier status and storage time were not significantly associated with increased bacterial growth (P > 0.05). CONCLUSIONS: Unfortified and fortified HPDHM remain largely free of bacterial growth for up to 96 hours of refrigerated storage in NICU settings. Sample handling techniques are important for preventing microbial contamination.

Infection with the Lyme disease pathogen suppresses innate immunity in mice with diet-induced obesity.

Obesity is a major global public health concern. Immune responses implicated in obesity also control certain infections. We investigated the effects of high-fat diet-induced obesity (DIO) on infection with the Lyme disease bacterium Borrelia burgdorferi in mice. DIO was associated with systemic suppression of neutrophil- and macrophage-based innate immune responses. These included bacterial uptake and cytokine production, and systemic, progressive impairment of bacterial clearance, and increased carditis severity. B. burgdorferi-infected mice fed normal diet also gained weight at the same rate as uninfected mice fed high-fat diet, toll-like receptor 4 deficiency rescued bacterial clearance defects, which greater in female than male mice, and killing of an unrelated bacterium (Escherichia coli) by bone marrow-derived macrophages from obese, B. burgdorferi-infected mice was also affected. Importantly, innate immune suppression increased with infection duration and depended on cooperative and synergistic interactions between DIO and B. burgdorferi infection. Thus, obesity and B. burgdorferi infection cooperatively and progressively suppressed innate immunity in mice.

Rare slips in fluctuating synchronized oscillator networks.

We study rare phase slips due to noise in synchronized Kuramoto oscillator networks. In the small-noise limit, we demonstrate that slips occur via large fluctuations to saddle phase-locked states. For tree topologies, slips appear between subgraphs that become disconnected at a saddle-node bifurcation, where phase-locked states lose stability generically. This pattern is demonstrated for sparse networks with several examples. Scaling laws are derived and compared for different tree topologies. On the other hand, for dense networks slips occur between oscillators on the edges of the frequency distribution. If the distribution is discrete, the probability-exponent for large fluctuations to occur scales linearly with the system size. However, if the distribution is continuous, the probability is a constant in the large network limit, as individual oscillators fluctuate to saddles while all others remain fixed. In the latter case, the network's coherence is approximately preserved.

Transcriptome Assessment of Erythema Migrans Skin Lesions in Patients With Early Lyme Disease Reveals Predominant Interferon Signaling.

Background: The most common clinical manifestation of early Lyme disease is the erythema migrans (EM) skin lesion that develops at the tick bite site typically between 7 and 14 days after infection with Borreliella burgdorferi. The host-pathogen interactions that occur in the skin may have a critical role in determining outcome of infection. Methods: Gene arrays were used to characterize the global transcriptional alterations in skin biopsy samples of EM lesions from untreated adult patients with Lyme disease in comparison to controls. Results: The transcriptional pattern in EM biopsies consisted of 254 differentially regulated genes (180 induced and 74 repressed) characterized by the induction of chemokines, cytokines, Toll-like receptors, antimicrobial peptides, monocytoid cell activation markers, and numerous genes annotated as interferon (IFN)-inducible. The IFN-inducible genes included 3 transcripts involved in tryptophan catabolism (IDO1, KMO, KYNU) that play a pivotal role in immune evasion by certain other microbial pathogens by driving the differentiation of regulatory T cells. Conclusions: This is the first study to globally assess the human skin transcriptional response during early Lyme disease. Borreliella burgdorferi elicits a predominant IFN signature in the EM lesion, suggesting a potential mechanism for spirochetal dissemination via IDO1-mediated localized immunosuppression.

Stage-specific global alterations in the transcriptomes of Lyme disease spirochetes during tick feeding and following mammalian host adaptation.

Borrelia burgdorferi, the agent of Lyme disease, is maintained in nature within an enzootic cycle involving a mammalian reservoir and an Ixodes sp. tick vector. The transmission, survival and pathogenic potential of B. burgdorferi depend on the bacterium's ability to modulate its transcriptome as it transits between vector and reservoir host. Herein, we employed an amplification-microarray approach to define the B. burgdorferi transcriptomes in fed larvae, fed nymphs and in mammalian host-adapted organisms cultivated in dialysis membrane chambers. The results show clearly that spirochetes exhibit unique expression profiles during each tick stage and during cultivation within the mammal; importantly, none of these profiles resembles that exhibited by in vitro grown organisms. Profound shifts in transcript levels were observed for genes encoding known or predicted lipoproteins as well as proteins involved in nutrient uptake, carbon utilization and lipid synthesis. Stage-specific expression patterns of chemotaxis-associated genes also were noted, suggesting that the composition and interactivities of the chemotaxis machinery components vary considerably in the feeding tick and mammal. The results as a whole make clear that environmental sensing by B. burgdorferi directly or indirectly drives an extensive and tightly integrated modulation of cell envelope constituents, chemotaxis/motility machinery, intermediary metabolism and cellular physiology. These findings provide the necessary transcriptional framework for delineating B. burgdorferi regulatory pathways throughout the enzootic cycle as well as defining the contribution(s) of individual genes to spirochete survival in nature and virulence in humans.

Borrelia burgdorferi induces a type I interferon response during early stages of disseminated infection in mice.

BACKGROUND: Lyme borrelia genotypes differ in their capacity to cause disseminated disease. Gene array analysis was employed to profile the host transcriptome induced by Borrelia burgdorferi strains with different capacities for causing disseminated disease in the blood of C3H/HeJ mice during early infection. RESULTS: B. burgdorferi B515, a clinical isolate that causes disseminated infection in mice, differentially regulated 236 transcripts (P < 0.05 by ANOVA, with fold change of at least 2). The 216 significantly induced transcripts included interferon (IFN)-responsive genes and genes involved in immunity and inflammation. In contrast, B. burgdorferi B331, a clinical isolate that causes transient skin infection but does not disseminate in C3H/HeJ mice, stimulated changes in only a few genes (1 induced, 4 repressed). Transcriptional regulation of type I IFN and IFN-related genes was measured by quantitative RT-PCR in mouse skin biopsies collected from the site of infection 24 h after inoculation with B. burgdorferi. The mean values for transcripts of Ifnb, Cxcl10, Gbp1, Ifit1, Ifit3, Irf7, Mx1, and Stat2 were found to be significantly increased in B. burgdorferi strain B515-infected mice relative to the control group. In contrast, transcription of these genes was not significantly changed in response to B. burgdorferi strain B331 or B31-4, a mutant that is unable to disseminate. CONCLUSIONS: These results establish a positive association between the disseminating capacity of B. burgdorferi and early type I IFN induction in a murine model of Lyme disease.

Epidemic Extinction and Control in Heterogeneous Networks.

We consider epidemic extinction in finite networks with a broad variation in local connectivity. Generalizing the theory of large fluctuations to random networks with a given degree distribution, we are able to predict the most probable, or optimal, paths to extinction in various configurations, including truncated power laws. We find that paths for heterogeneous networks follow a limiting form in which infection first decreases in low-degree nodes, which triggers a rapid extinction in high-degree nodes, and finishes with a residual low-degree extinction. The usefulness of our approach is further demonstrated through optimal control strategies that leverage the dependence of finite-size fluctuations on network topology. Interestingly, we find that the optimal control is a mix of treating both high- and low-degree nodes based on theoretical predictions, in contrast to methods that ignore dynamical fluctuations.