A bacterial genome in flux: the twelve linear and nine circular extrachromosomal DNAs in an infectious isolate of the Lyme disease spirochete Borrelia burgdorferi

Diagnostic validation of novel Borrelia antigens discovered by whole-proteome microarray: Advancing early detection and test of cure for Lyme disease

Lyme disease serodiagnosis has limited early sensitivity and cannot distinguish active from past infections. To address this, we screen a Borrelia afzelii whole-proteome microarray (1,296 proteins) using human (n = 149) and murine (n = 32) sera. We evaluate three early-stage antigens-BafPKo_A0001, BafPKo_D0016, and BafPKo_A0029. ELISA cutoffs are established using discovery cohort sera (n = 99) and validated with the validation (n = 242) and the prospective (n = 223) cohorts. A0001 demonstrates 87.8% sensitivity, outperforming C6 (69.4%) and STTT (22.5%) in the discovery cohort. In the validation cohort, A0001 reaches 90.5% sensitivity, surpassing C6 by 11.6% and STTT by 50%. In hyper-acute erythema migrans sera (from the prospective cohort), A0001 achieves 55.1% sensitivity, exceeding C6 and STTT by 14.6% and 33.3%, respectively. COMBO-3 and COMBO-2 yield the highest sensitivity of 92.9% and 66.1% in the validation and prospective cohort, respectively. A0001 and D0016 show enhanced and robust seroreversion after antibiotic treatment suggesting their potential as test of cure biomarkers in early Lyme disease.

Characterization of the family-level Borreliaceae pan-genome and development of an episomal typing protocol

The Borreliaceae family includes many obligate parasitic bacterial species etiologically associated with a myriad of zoonotic borrelioses, including Lyme disease and vector-borne relapsing fevers. Borreliaceae infections are difficult to detect by both direct and indirect methods, often leading to delayed and missed diagnoses. Efforts to improve diagnostics center around the development of molecular diagnostics (MDx), but due to deep tissue sequestration and the lack of persistent bacteremias, even MDx assays suffer from a lack of sensitivity. Additionally, the extensive genomic heterogeneity among isolates, even within the same species, contributes to the lack of assay sensitivity, as single target assays, whether nucleic acid-based or serologically based, cannot provide universal coverage. This within-species heterogeneity is partly due to differences in replicon repertoires and genomic structures that have likely arisen to support the complex Borreliaceae life cycle necessary for these parasites to survive in multiple hosts, each with unique immune responses. We constructed a Borreliaceae family-level pan-genome and characterized the phylogenetic relationships among the constituent taxa, which supports the recent, although contested, taxonomy of splitting the family into at least two genera. Gene content profiles were created for the majority of the Borreliaceae replicons, providing for the first time their unambiguous molecular typing. Our characterization of the Borreliaceae pan-genome supports the splitting of the former Borrelia genus into two genera and provides for the phylogenetic placement of several non-species designated isolates. Mining this family-level pan-genome will enable the development of precision diagnostics corresponding to gene content-driven clinical outcomes while also providing targets for interventions.

IMPORTANCE

Using whole genome sequencing, we demonstrated that the bacteria that are transmitted by ticks and other arthropod vectors that cause Lyme disease and relapsing fevers, while related, do not belong within the same genus classification. In addition, through characterization of their highly atypical genomic structure, we were able to develop a genetic typing system that will help with future studies of how they cause disease while also providing targets for medical interventions.

The peptidoglycan of Borrelia burgdorferi can persist in discrete tissues and cause systemic responses consistent with chronic illness

Persistent symptoms after an acute infection is an emerging public health concern, but the pathobiology of such conditions is not well understood. One possible scenario involves the persistence of lingering antigen. We have previously reported that patients with postinfectious Lyme arthritis often harbor the peptidoglycan (PG) cell wall of Borrelia burgdorferi, the Lyme disease agent, in the synovial fluid of their inflamed joints after treatment. However, it is not yet known how B. burgdorferi PG persists, in what form, or if it may play a role in other postinfectious complications after Lyme disease. Using a murine model, we developed a real-time in vivo system to track B. burgdorferi PG as a function of cell wall chemistry and validated our findings using both molecular and cellular approaches. Unlike typical bacterial PG, the unique chemical properties of polymeric B. burgdorferi PG drive murine liver accumulation, where the cell wall material persists for weeks. Kupffer cells and hepatocytes phagocytose and retain B. burgdorferi PG and, although liver occupancy coincides with minimal pathology, both organ-specific and secreted protein profiles produced under these conditions bear some similarities to reported proteins enriched in patients with chronic illness after acute infection. Moreover, transcriptomic profiling indicated that B. burgdorferi PG affects energy metabolism in peripheral blood mononuclear cells. Our findings provide mechanistic insights into how a pathogenic molecule can persist after agent clearance, potentially contributing to illness after infection.

Lactate dehydrogenase is the Achilles' heel of Lyme disease bacterium Borreliella burgdorferi

As a zoonotic pathogen, the Lyme disease bacterium Borreliella burgdorferi has evolved unique metabolic pathways, some of which are specific and essential for its survival and thus present as ideal targets for developing new therapeutics. B. burgdorferi dispenses with the use of thiamin as a cofactor and relies on lactate dehydrogenase (BbLDH) to convert pyruvate to lactate for balancing NADH/NAD+ ratios. This report first demonstrates that BbLDH is a canonical LDH with some unique biochemical and structural features. A loss-of-function study then reveals that BbLDH is essential for B. burgdorferi survival and infectivity, highlighting its therapeutic potential. Drug screening identifies four previously unknown LDH inhibitors with minimal cytotoxicity, two of which inhibit B. burgdorferi growth. This study provides mechanistic insights into the function of BbLDH in the pathophysiology of B. burgdorferi and lays the groundwork for developing genus-specific metabolic inhibitors against B. burgdorferi and potentially other tick-borne pathogens as well.

IMPORTANCE

Lyme disease (LD) is the most commonly reported tick-borne illness in the U.S. and Europe, and its geographic distribution is continuously expanding worldwide. Though early LD can be treated with antibiotics, chronic LD is recalcitrant to antibiotic treatments and thus requires multiple courses of antibiotic therapy. Currently, there are no human vaccines nor prophylactic antibiotics to prevent LD. As the causative agent of LD, Borreliella burgdorferi has evolved unique metabolic pathways, some of which are specific and essential for its survival and thus present as ideal targets for developing new therapeutics. By using an approach of genetics, biochemistry, structural biology, drug screening, and animal models, this report provides evidence that lactate dehydrogenase can be a potential target for developing genus-specific metabolic inhibitors against B. burgdorferi and potentially other tick-borne pathogens as well.

Borrelia burgdorferi Strain-Specific Differences in Mouse Infectivity and Pathology

Lyme disease (LD), caused by infection with the tick-borne bacteria, Borrelia burgdorferi, is associated with a wide array of symptoms in human patients. Variations in clinical manifestations are thought to be influenced by genetic differences among B. burgdorferi strains. In this study, we evaluated the infectivity, tissue bacterial load, pathology, and immunogenicity of five strains of B. burgdorferi sensu stricto (297 Ah130, Bb16-54, B31-A3, Bb16-126, JD1) in female C3H/HeN mice at three infectious doses (104, 105, 106 spirochetes). We found that strains Bb16-126 and JD1 were the most infectious, resulting in 100% infection across all the tested doses. Strain Bb16-126 caused the highest bacterial burden in the heart tissue and significant carditis, whereas JD1 exhibited the lowest spirochete load in the heart and minimal carditis. In comparison, strain B31-A3 demonstrated the highest abundance in the tibiotarsal joint. Infection with all the strains induced severe lymph node hyperplasia, with JD1 producing the greatest increase in cellularity. Using a diagnostic C6 peptide ELISA, all the strains induced significant anti-C6 IgM and IgG antibody titers at 14 days post-infection; however, strain B31-A3 elicited the highest anti-C6 IgM titers. Our findings demonstrate the importance of strain diversity in shaping B. burgdorferi pathogenesis in a mouse model and provide insights for developing strain-specific diagnostic, therapeutic, and vaccine strategies.

Divergent dimerization mechanisms and conserved DNA-binding function in PFam12 proteins of Borrelia burgdorferi

Lyme disease, caused by the spirochete Borrelia burgdorferi, is transmitted to mammalian hosts during the feeding process of infected Ixodes ticks. Our previous studies demonstrated that the paralogous gene family 12 (PFam12) consisting of five members (BBK01, BBG01, BBH37, BBJ08, and BB0844) are non-specific DNA-binding proteins. PFam12 proteins share 31-69% sequence identity, are located either on the surface or within the periplasm and are upregulated as the tick starts its blood meal. The crystal structure of BBK01 revealed that the protein forms a homodimer, which is potentially critical for DNA binding. In this study, we determined the crystal structure of another PFam12 member, BBH37, to gain a better insight into this unique paralogous family. Although BBK01 dimerization is mediated by its C-terminal region and is thought to be critical for DNA binding, BBH37 forms dimers through an alternative mechanism where a unique disulfide bond is involved. We found that BBH37 is still able to interact with DNA with micromolar affinity. Molecular dynamics simulations and site-directed mutagenesis was conducted to characterize these unique DNA binding proteins. This study highlights the structural diversity within the PFam12, demonstrating that despite significant differences in dimerization mechanisms, these proteins retain their DNA-binding capability.

Borrelia burgdorferi lacking all cp32 prophage plasmids retains full infectivity in mice

The causative agent of Lyme disease, Borrelia burgdorferi, contains a unique, segmented genome comprising multiple linear and circular plasmids. To date, the genomes of over 63 sequenced Lyme disease Borrelia carry one or more 32 kbp circular plasmids (cp32) or cp32-like elements. The cp32 plasmids are endogenous prophages and encode, among other elements, a family of surface exposed lipoproteins termed OspEF-related proteins. These lipoproteins are synthesized during mammalian infection and are considered important components of the spirochete's adaptive response to the vertebrate host. Here, we detail the construction and infectivity of the first described B. burgdorferi strain lacking all cp32 plasmids. Despite their universal presence, our findings indicate that B. burgdorferi does not require any cp32 plasmids to complete the experimental mouse-tick-mouse infectious cycle and a total lack of cp32s does not impair spirochete infectivity.

Borrelia burgdorferi loses essential genetic elements and cell proliferative potential during stationary phase in culture but not in the tick vector

The Lyme disease agent Borrelia burgdorferi is a polyploid bacterium with a segmented genome in which both the chromosome and over 20 distinct plasmids are present in multiple copies per cell. This pathogen can survive for at least 9 months in its tick vector in an apparent dormant state between blood meals, without losing cell proliferative capability when re-exposed to nutrients. Cultivated B. burgdorferi cells grown to stationary phase or resuspended in nutrient-limited media are often used to study the effects of nutrient deprivation. However, a thorough assessment of the spirochete's ability to recover from nutrient depletion has been lacking. Our study shows that starved B. burgdorferi cultures rapidly lose cell proliferative ability. Loss of genetic elements essential for cell proliferation contributes to the observed proliferative defect in stationary phase. The gradual decline in copies of genetic elements is not perfectly synchronized between chromosomes and plasmids, generating cells that harbor one or more copies of the essential chromosome but lack all copies of one or more non-essential plasmids. This phenomenon likely contributes to the well-documented issue of plasmid loss during in vitro cultivation of B. burgdorferi. In contrast, B. burgdorferi cells from ticks starved for 14 months showed no evidence of reduced cell proliferative ability or plasmid loss. Beyond their practical implications for studying B. burgdorferi, these findings suggest that the midgut of the tick vector offers a unique environment that supports the maintenance of B. burgdorferi's segmented genome and cell proliferative potential during periods of tick fasting.IMPORTANCEBorrelia burgdorferi causes Lyme disease, a prevalent tick-borne illness. B. burgdorferi must survive long periods (months to a year) of apparent dormancy in the midgut of the tick vector between blood meals. Resilience to starvation is a common trait among bacteria. However, this study reveals that, in laboratory cultures, B. burgdorferi poorly endures starvation and rapidly loses viability. This decline is linked to a gradual loss of genetic elements required for cell proliferation. These results suggest that the persistence of B. burgdorferi in nature is likely shaped more by unique environmental conditions in the midgut of the tick vector than by an innate ability of this bacterium to endure nutrient deprivation.

Dissection of amino acid acquisition pathways in Borrelia burgdorferi uncovers unique physiological responses

Borrelia burgdorferi, the causative agent of Lyme disease, is well known for its unique morphology, physiology, and enzootic lifecycle. Building on previous work that showed peptide transport is essential for viability, we endeavored to more clearly define the impact of peptide starvation on the spirochete and directly compare peptide starvation to targeted free amino acid starvation. Herein, we confirm the ability of a putative GltP, BB0401, to facilitate transport of glutamate and aspartate as well as demonstrate its requirement for cell growth and motility. Using conditional mutants for both peptide transport and BB0401, we characterize these systems throughout the enzootic cycle, both confirming their essential role during murine infection and revealing that they are, surprisingly, dispensable during prolonged colonization of the tick midgut. We broadly define the metabolic perturbations resulting from these amino acid starvation models and show that, even under the most severe amino acid stress, B. burgdorferi is unable to modulate its physiological response via the canonical (p)ppGpp-driven stringent response.

Borrelial phosphomannose isomerase as a cell surface localized protein that retains enzymatic activity and promotes host-pathogen interaction

All organisms produce an intracellular Zn2+-dependent enzyme, phosphomannose isomerase (PMI) or mannose-6 phosphate isomerase, that catalyzes the reversible conversion of mannose-6-phosphate and fructose-6-phosphate during sugar metabolism and polysaccharide biosynthesis. Unexpectedly, we discovered an additional PMI function in Borrelia burgdorferi, the pathogen of Lyme disease, where the enzyme is localized on the cell surface and binds to collagen IV-a host extracellular matrix component predominantly found in the skin. The AlphaFold 3-based structural model of B. burgdorferi PMI (BbPMI) retains the active site with tetrahedrally-coordinated Zn2+ seen in other PMIs of known structure, residing in an elongated crevice. Ligand docking shows that the crevice can accommodate the tip trisaccharide moiety of a glycosylated asparagine residue on the collagen IV 7S domain. Low doses of a well-known PMI benzoisothiazolone inhibitor impair the growth of diverse strains of B. burgdorferi in culture, but not other tested Gram-negative or Gram-positive pathogens. Borrelia cells are even more susceptible to several other structurally related benzoisothiazolone analogs. The passive transfer of anti-BbPMI antibodies in ticks can impact spirochete transmission to mice, while the treatment of collagen IV-containing murine skin with PMI inhibitors impairs spirochete infectivity. Taken together, these results highlight a newly discovered role for BbPMI in mediating host-pathogen interactions during the spirochete infectivity process. In turn, this discovery offers an opportunity for the development of a novel therapeutic strategy to combat Lyme disease by preventing the BbPMI interaction with its host receptor, collagen IV.

IMPORTANCE

All organisms produce an intracellular enzyme, phosphomannose isomerase (PMI), that converts specific sugars during metabolism. Unexpectedly, we discovered an additional PMI function in Borrelia burgdorferi, the Lyme disease pathogen, where the enzyme is localized on the cell surface and binds to collagen IV-a host extracellular molecule mainly found in the skin. Low doses of PMI chemical inhibitors impair the growth of diverse strains of B. burgdorferi in culture, but not other tested bacterial pathogens. The passive transfer of anti-BbPMI antibodies in ticks can impact B. burgdorferi transmission to mice, while the treatment of collagen IV-containing murine skin with PMI inhibitors impairs infectivity. Taken together, these results highlight a newly discovered role for BbPMI in mediating host-pathogen interactions during infection. In turn, this discovery offers an opportunity for the development of a novel therapeutic strategy to combat Lyme disease by preventing BbPMI function and interaction with host collagen IV.

Analysis of the BadR regulon in Borrelia burgdorferi

BACKGROUND

Borrelia burgdorferi, the causative agent of Lyme disease, relies on tightly coordinated gene expression to quickly adapt and survive in the tick vector and mammalian host. BadR, an ROK (repressor, open reading frame, kinase) family transcriptional regulator, binds directly to B. burgdorferi promoter DNA, however, many questions concerning the role for BadR in gene regulation remain unanswered. In particular, there are conflicting reports concerning what genes are regulated by BadR in B. burgdorferi. Furthermore, previous studies have suggested important roles for BadR in unfed ticks, but the BadR regulon has not been defined under such conditions. Additionally, although BadR regulates rpoS expression in a growth phase-dependent manner, it remains unknown whether BadR regulates other genes during different growth phases.

RESULTS

To address these questions, we cultivated a B. burgdorferi badR mutant and wild-type strain under various conditions and analyzed the transcriptome using RNA-sequencing. When spirochetes were grown at 37 °C and collected at the mid-logarithmic and stationary phase of growth, 211 and 272 genes were differentially expressed in the badR mutant, respectively. A total of 79 genes were differentially expressed when spirochetes were grown at 23 °C. A vast majority of genes identified in this study encode proteins of unknown function.

CONCLUSIONS

Complex transcriptional regulation mechanisms coordinate the expression of genes required for the survival of B. burgdorferi throughout its tick-mammal enzootic lifecycle. As part of this process, BadR functions as a global regulatory protein and regulates B. burgdorferi virulence gene expression. Combined, this work supports a role for BadR in global B. burgdorferi gene regulation by modulating expression of different sets of genes at different stages of the enzootic lifecycle. We anticipate that investigating the function of genes in the BadR regulon will lead to the identification of novel virulence factors for therapeutic and vaccine development.

Lactate dehydrogenase is the Achilles' heel of Lyme disease bacterium Borreliella burgdorferi

As a zoonotic pathogen, the Lyme disease bacterium Borreliella burgdorferi has evolved unique metabolic pathways, some of which are specific and essential for its survival and thus present as ideal targets for developing new therapeutics. B. burgdorferi dispenses with the use of thiamin as a cofactor and relies on lactate dehydrogenase (BbLDH) to convert pyruvate to lactate for balancing NADH/NAD + ratios. This report first demonstrates that BbLDH is a canonical LDH with some unique biochemical and structural features. A loss-of-function study then reveals that BbLDH is essential for B. burgdorferi survival and infectivity, highlighting its therapeutic potential. Drug screening identifies four previously unknown LDH inhibitors with minimal cytotoxicity, two of which inhibit B. burgdorferi growth. This study provides mechanistic insights into the function of BbLDH in the pathophysiology of B. burgdorferi and lays the groundwork for developing genus-specific metabolic inhibitors against B. burgdorferi and potentially other tick-borne pathogens as well.

Evolutionary rescue of bacterial populations by heterozygosity on multicopy plasmids

Bacterial plasmids and other extrachromosomal DNA elements frequently carry genes with important fitness effects for their hosts. Multicopy plasmids can additionally carry distinct alleles of host-fitness-relevant genes on different plasmid copies, allowing for heterozygosity not possible for loci on haploid chromosomes. Plasmid-mediated heterozygosity may increase the fitness of bacterial cells in circumstances where there is an advantage to having multiple distinct alleles (heterozyogote advantage); however, plasmid-mediated heterozygosity is also subject to constant loss due to random segregation of plasmid copies on cell division. We analyze a multitype branching process model to study the evolution and maintenance of plasmid-mediated heterozygosity under a heterozygote advantage. We focus on an evolutionary rescue scenario in which a novel mutant allele on a plasmid must be maintained together with the wild-type allele to allow population persistance (although our results apply more generally to the maintenance of heterozygosity due to heterozygote advantage). We determine the probability of rescue and derive an analytical expression for the threshold on the fitness of heterozygotes required to overcome segregation and make rescue possible; this threshold decreases with increasing plasmids copy number. We further show that the formation of cointegrates from the fusion of plasmid copies increases the probability of rescue. Overall, our results provide a rigorous quantitative assessment of the conditions under which bacterial populations can adapt to multiple stressors through plasmid-mediated heterozygosity. Many of the results are furthermore applicable to the related problem of the maintenance of incompatible plasmids in the same cell under selection for both.

A core genome MLST scheme for Borrelia burgdorferi sensu lato improves insights into the evolutionary history of the species complex

Multi-locus sequence typing (MLST) based on eight genes has become the method of choice for Borrelia typing and is extensively used for population studies. Whole-genome sequencing enables studies to scale up to genomic levels but necessitates extended schemes. We have developed a 639-loci core genome MLST (cgMLST) scheme for Borrelia burgdorferi sensu lato (s.l.) that enables unambiguous genotyping and improves the robustness of phylogenies and lineage resolution within species. Notably, all inner nodes of the cgMLST phylogenies had consistently high statistical support. Analyses of the genetically homogeneous European B. bavariensis population support the notion that cgMLST provides high discriminatory power even for closely related isolates. While isolates differed maximally in one MLST locus, there were up to 179 cgMLST loci differences. Thus, the developed cgMLST scheme for B. burgdorferi s.l. resolves lineages at a finer resolution than MLST and improves insights into the evolutionary history of the species complex.

Borrelia PeptideAtlas: A proteome resource of common Borrelia burgdorferi isolates for Lyme research

Lyme disease is caused by an infection with the spirochete Borrelia burgdorferi, and is the most common vector-borne disease in North America. B. burgdorferi isolates harbor extensive genomic and proteomic variability and further comparison of isolates is key to understanding the infectivity of the spirochetes and biological impacts of identified sequence variants. Here, we applied both transcriptome analysis and mass spectrometry-based proteomics to assemble peptide datasets of B. burgdorferi laboratory isolates B31, MM1, and the infective isolate B31-5A4, to provide a publicly available Borrelia PeptideAtlas. Included are total proteome, secretome, and membrane proteome identifications of the individual isolates. Proteomic data collected from 35 different experiment datasets, totaling 386 mass spectrometry runs, have identified 81,967 distinct peptides, which map to 1,113 proteins. The Borrelia PeptideAtlas covers 86% of the total B31 proteome of 1,291 protein sequences. The Borrelia PeptideAtlas is an extensible comprehensive peptide repository with proteomic information from B. burgdorferi isolates useful for Lyme disease research.

Analysis of bb0556 Expression and Its Role During Borrelia burgdorferi Mammalian Infection

In Borrelia burgdorferi, BB0556 was annotated as a conserved hypothetical protein. We herein investigated gene expression and the importance of this protein during infection. Our data support that bb0556 forms an operon with five other genes. A transcriptional start site and the associated σ70-type promoter were identified in the sequences upstream of bb0554, and luciferase reporter assays indicated that this promoter is functional in B. burgdorferi. Furthermore, the sequences upstream of bb0556 contain an internal promoter to drive gene expression. bb0556 expression was affected by various environmental factors such as changes in temperature, pH, and cell density when B. burgdorferi was grown in vitro. Surprisingly, significant differences were observed for bb0556 expression between B. burgdorferi strains B31-A3 and CE162, likely due to the different cis- and trans-acting factors in these strains. Moreover, bb0556 was found to be highly expressed by B. burgdorferi in infected mice tissues, suggesting that this gene plays an important role during animal infection. To test this hypothesis, we generated a bb0556 deletion mutant in a virulent bioluminescent B. burgdorferi strain. The mutant grew normally in the medium and displayed no defect in the resistance to environmental stresses such as reactive oxygen species, reactive nitrogen species, and osmotic stress. However, when the infectivity was compared between the mutant and its parental strain using in vivo bioluminescence imaging as well as analyses of spirochete recovery and bacterial burdens in animal tissues, our data showed that, contrary to the parental strain, the mutant was unable to infect mice. Complementation of bb0556 in cis fully restored the infectious phenotype to wild-type levels. Taken together, our study demonstrates that the hypothetical protein BB0556 is a novel virulence factor essential for B. burgdorferi mammalian infection.

MCP5, a methyl-accepting chemotaxis protein regulated by both the Hk1-Rrp1 and Rrp2-RpoN-RpoS pathways, is required for the immune evasion of Borrelia burgdorferi

Borrelia (or Borreliella) burgdorferi, the causative agent of Lyme disease, is a motile and invasive zoonotic pathogen adept at navigating between its arthropod vector and mammalian host. While motility and chemotaxis are well known to be essential for its enzootic cycle, the role of each methyl-accepting chemotaxis proteins (MCPs) in the infectious cycle of B. burgdorferi remains unclear. In this study, we show that mcp5, a gene encoding one of the most abundant MCPs in B. burgdorferi, is differentially expressed in response to environmental signals and at distinct stages of the pathogen's enzootic cycle. Notably, mcp5 expression is regulated by the Hk1-Rrp1 and Rrp2-RpoN-RpoS pathways, two key regulatory pathways that are critical for the spirochete's colonization of the tick vector and mammalian host, respectively. Infection experiments with an mcp5 mutant revealed that spirochetes lacking MCP5 were unable to establish infections in either C3H/HeN mice or Severe Combined Immunodeficiency (SCID) mice, which are deficient in adaptive immunity, underscoring MCP5's critical role in mammalian infection. However, the mcp5 mutant was able to establish infection and disseminate in NOD SCID Gamma (NSG) mice, which are deficient in both adaptive and most innate immune responses, suggesting that MCP5 plays an important role in evading host innate immunity. Moreover, NK cell depletion in C3H and SCID mice restored the infectivity of the mcp5 mutant, further highlighting MCP5's role in evading NK cell-associated immunity. Co-culture assays with NK cells and macrophages revealed that the mcp5 mutant enhanced interferon-gamma production by NK cells. In the tick vector, the mcp5 mutants survived feeding but failed to transmit to mice. These findings reveal that MCP5, regulated by both the Rrp1 and Rrp2 pathways, is critical for establishing infection in mammalian hosts by evading NK cell-mediated host innate immunity and is important for the transmission of spirochetes from ticks to mammalian hosts. This work provides a foundation for further elucidation of chemotactic signals sensed by MCP5 that facilitate B. burgdorferi in evading host defenses.

Role of Dual Specificity Phosphatase 1 (DUSP1) in influencing inflammatory pathways in macrophages modulated by Borrelia burgdorferi lipoproteins

Borrelia burgdorferi (Bb), the spirochetal agent of Lyme disease, has a large array of lipoproteins that play a significant role in mediating host-pathogen interactions within ticks and vertebrates. Although there is substantial information on the effects of B. burgdorferi lipoproteins (BbLP) on immune modulatory pathways, the application of multi-omics methodologies to decode the transcriptional and proteomic patterns associated with host cell responses induced by lipoproteins in murine bone marrow-derived macrophages (BMDMs) has identified additional effectors and pathways. Single-cell RNA-Seq (scRNA-Seq) performed on BMDMs treated with various concentrations of borrelial lipoproteins revealed macrophage subsets within the BMDMs. Differential expression analysis showed that genes encoding various receptors, type I IFN-stimulated genes, signaling chemokines, and mitochondrial genes are altered in BMDMs in response to lipoproteins. Unbiased proteomics analysis of lysates of BMDMs treated with lipoproteins corroborated several of these findings. Notably, dual specificity phosphatase 1 (Dusp1) gene was upregulated during the early stages of BMDM exposure to BbLP. Pre-treatment with benzylidene-3-cyclohexylamino-1-indanone hydrochloride (BCI), an inhibitor of both DUSP1 and 6 prior to exposure to BbLP, demonstrated that DUSP1 negatively regulates NLRP3-mediated pro-inflammatory signaling and positively regulates the expression of interferon-stimulated genes and those encoding Ccl5, Il1b, and Cd274. Moreover, DUSP1, IkB kinase complex and MyD88 also modulate mitochondrial changes in BMDMs treated with borrelial lipoproteins. These findings advance the potential for exploiting DUSP1 as a therapeutic target to regulate host responses in reservoir hosts to limit survival of B. burgdorferi during its infectious cycle between ticks and mammalian hosts.

Borrelia burgdorferi loses essential genetic elements and cell proliferative potential during stationary phase in culture but not in the tick vector

The Lyme disease agent Borrelia burgdorferi is a polyploid bacterium with a segmented genome in which both the chromosome and over 20 distinct plasmids are present in multiple copies per cell. This pathogen can survive at least nine months in its tick vector in an apparent dormant state between blood meals, without losing cell proliferative capability when re-exposed to nutrients. Cultivated B. burgdorferi cells grown to stationary phase or resuspended in nutrient-limited media are often used to study the effects of nutrient deprivation. However, a thorough assessment of the spirochete's ability to recover from nutrient depletion has been lacking. Our study shows that starved B. burgdorferi cultures rapidly lose cell proliferative. Loss of genetic elements essential for cell proliferation contributes to the observed proliferative defect in stationary phase. The gradual decline in copies of genetic elements is not perfectly synchronized between chromosomes and plasmids, generating cells that harbor one or more copies of the essential chromosome but lack all copies of one or more non-essential plasmids. This phenomenon likely contributes to the well-documented issue of plasmid loss during in vitro cultivation of B. burgdorferi. In contrast, B. burgdorferi cells from ticks starved for 14 months showed no evidence of reduced cell proliferative ability or plasmid loss. Beyond their practical implications for studying B. burgdorferi, these findings suggest that the midgut of the tick vector offers a unique environment that supports the maintenance of B. burgdorferi's segmented genome and cell proliferative potential during periods of tick fasting.

The plasmid-encoded members of paralogous gene family 52 are dispensable to the enzootic cycle of Borrelia burgdorferi

Lyme disease, the leading vector-borne disease in the United States and Europe, develops after infection with Borrelia burgdorferi sensu lato bacteria. Transmission of the spirochete from the tick vector to a vertebrate host requires global changes in gene expression that are controlled, in part, by the Rrp2/RpoN/RpoS alternative sigma factor cascade. Transcriptional studies defining the B. burgdorferi RpoS regulon have suggested that RpoS activates the transcription of paralogous family 52 (PFam52) genes. In strain B31, PFam52 genes (bbi42, bbk53, and bbq03) encode a set of conserved hypothetical proteins with >89% amino acid identity that are predicted to be surface-localized. Extensive homology among members of paralogous families complicates studies of protein contributions to pathogenicity as the potential for functional redundancy will obfuscate findings. Using a sequential mutagenesis approach, we generated clones expressing a single PFam52 paralog, as well as a strain deficient in all three. The single paralog expressing strains were used to confirm BBI42, BBK53, and BBQ03 surface localization and RpoS regulation. Surprisingly, the PFam52-deficient strain was able to infect mice and complete the enzootic cycle similar to the wild-type parental strain. Indeed, the presence of numerous pseudogenes that contain frameshifts or internal stop codons among the PFam52 genes suggests that they may be subjected to gene loss in B. burgdorferi's reduced genome. Alternatively, the lack of phenotype might reflect the limitations of the experimental mouse infection model.

Meta-analysis of the Vmp-like sequences of Lyme disease Borrelia: evidence for the evolution of an elaborate antigenic variation system

VMP-like sequence (vls) antigenic variation systems are present in every Lyme disease Borrelia strain with complete genome sequences. The linear plasmid-encoded vls system consists of a single expression site (vlsE) and contiguous array(s) of silent cassettes that have ~90% identity with the central cassette region of the cognate vlsE gene; antigenic variation occurs through random, segmental, and unidirectional recombination of vls silent cassette sequences into the vlsE expression site. Automated annotation programs do not accurately recognize vls silent cassette sequences, so these regions are not correctly annotated in most genomic sequences. In this study, the vls sequences were re-analyzed in the genomic sequences of 31 available Lyme disease Borrelia and one relapsing fever Borrelia organisms, and this information was utilized to systematically compare the vls systems in different species and strains. In general, the results confirm the conservation of the overall architecture of the vls system, such as the head-to-head arrangement of vlsE and a contiguous series of vlsS silent cassette sequences and presence of inverted repeat sequences between the two regions. However, the data also provide evidence for the divergence of the vls silent cassette arrays through point mutations, short indels, duplication events, and rearrangements. The probable occurrence of convergent evolution toward a vls system-like locus is exemplified by Borrelia turcica, a variable large protein (Vlp) expressing organism that is a member of the relapsing fever Borrelia group.

Combining short- and long-read sequencing unveils geographically structured diversity in Borrelia miyamotoi

Borrelia miyamotoi is an emerging Ixodes tick-borne human pathogen in the Northern hemisphere. The aim of the current study was to compare whole genome sequences of B. miyamotoi isolates from different continents. Using a combination of Illumina and PacBio platforms and a novel genome assembly and plasmid typing pipeline, we reveal that the 21 sequenced B. miyamotoi isolates and publically available B. miyamotoi genomes from North America, Asia, and Europe form genetically distinct populations and cluster according to their geographical origin, where distinct Ixodes species are endemic. We identified 20 linear and 17 circular plasmid types and the presence of specific plasmids for isolates originating from different continents. Linear plasmids lp12, lp23, lp41, and lp72 were core plasmids found in all isolates, with lp41 consistently containing the vmp expression site. Our data provide insights into the genetic basis of vector competence, virulence, and pathogenesis of B. miyamotoi.

Hitchhiker's Guide to Borrelia burgdorferi

Don't Panic. In the nearly 50 years since the discovery of Lyme disease, Borrelia burgdorferi has emerged as an unlikely workhorse of microbiology. Interest in studying host-pathogen interactions fueled significant progress in making the fastidious microbe approachable in laboratory settings, including the development of culture methods, animal models, and genetic tools. By developing these systems, insight has been gained into how the microbe is able to survive its enzootic cycle and cause human disease. Here, we discuss the discovery of B. burgdorferi and its development as a model organism before diving into the critical lessons we have learned about B. burgdorferi biology at pivotal stages of its lifecycle: gene expression changes during the tick blood meal, colonization of a new vertebrate host, and developing a long-lasting infection in that vertebrate until a new tick feeds. Our goal is to highlight the advancements that have facilitated B. burgdorferi research and identify gaps in our current understanding of the microbe.

Natural selection and recombination at host-interacting lipoprotein loci drive genome diversification of Lyme disease and related bacteria

Lyme disease, caused by spirochetes in the Borrelia burgdorferi sensu lato clade within the Borrelia genus, is transmitted by Ixodes ticks and is currently the most prevalent and rapidly expanding tick-borne disease in Europe and North America. We report complete genome sequences of 47 isolates that encompass all established species in this clade while highlighting the diversity of the widespread human pathogenic species B. burgdorferi. A similar set of plasmids has been maintained throughout Borrelia divergence, indicating that they are a key adaptive feature of this genus. Phylogenetic reconstruction of all sequenced Borrelia genomes revealed the original divergence of Eurasian and North American lineages and subsequent dispersals that introduced B. garinii, B. bavariensis, B. lusitaniae, B. valaisiana, and B. afzelii from East Asia to Europe and B. burgdorferi and B. finlandensis from North America to Europe. Molecular phylogenies of the universally present core replicons (chromosome and cp26 and lp54 plasmids) are highly consistent, revealing a strong clonal structure. Nonetheless, numerous inconsistencies between the genome and gene phylogenies indicate species dispersal, genetic exchanges, and rapid sequence evolution at plasmid-borne loci, including key host-interacting lipoprotein genes. While localized recombination occurs uniformly on the main chromosome at a rate comparable to mutation, lipoprotein-encoding loci are recombination hotspots on the plasmids, suggesting adaptive maintenance of recombinant alleles at loci directly interacting with the host. We conclude that within- and between-species recombination facilitates adaptive sequence evolution of host-interacting lipoprotein loci and contributes to human virulence despite a genome-wide clonal structure of its natural populations.

IMPORTANCE

Lyme disease (also called Lyme borreliosis in Europe), a condition caused by spirochete bacteria of the genus Borrelia, transmitted by hard-bodied Ixodes ticks, is currently the most prevalent and rapidly expanding tick-borne disease in the United States and Europe. Borrelia interspecies and intraspecies genome comparisons of Lyme disease-related bacteria are essential to reconstruct their evolutionary origins, track epidemiological spread, identify molecular mechanisms of human pathogenicity, and design molecular and ecological approaches to disease prevention, diagnosis, and treatment. These Lyme disease-associated bacteria harbor complex genomes that encode many genes that do not have homologs in other organisms and are distributed across multiple linear and circular plasmids. The functional significance of most of the plasmid-borne genes and the multipartite genome organization itself remains unknown. Here we sequenced, assembled, and analyzed whole genomes of 47 Borrelia isolates from around the world, including multiple isolates of the human pathogenic species. Our analysis elucidates the evolutionary origins, historical migration, and sources of genomic variability of these clinically important pathogens. We have developed web-based software tools (BorreliaBase.org) to facilitate dissemination and continued comparative analysis of Borrelia genomes to identify determinants of human pathogenicity.

Analysis of the Borreliaceae Pangenome Reveals a Distinct Genomic Architecture Conserved Across Phylogenetic Scales

The family Borreliaceae contains arthropod-borne spirochetes that cause two widespread human diseases, Lyme disease and relapsing fever. Lyme disease is a subacute, progressive illness with variable stage and tissue manifestations. Relapsing fever is an acute febrile illness with prominent bacteremia that may recur and disseminate, particularly to the nervous system. Clinical heterogeneity is a hallmark of both diseases. While human clinical manifestations are influenced by a wide variety of factors, including immune status and host genetic susceptibility, there is evidence that Borreliaceae microbial factors influence the clinical manifestations of human disease caused by this family of spirochetes. Despite these associations, the spirochete genes that influence the severity and manifestations of human disease are, for the most part, unknown. Recent work has identified lineage-specific expansions of lipoproteome-rich accessory genome elements in virulent clones of Borrelia burgdorferi. Using publicly available genome assemblies, it is shown that all Borreliaceae lineages for which sufficient sequence data are available harbor a similar pattern of strongly structured, lineage-specific expansions in their accessory genomes, particularly among lipoproteins, and that this pattern holds across phylogenetic scales including genera, species, and genotypes. The relationships among pangenome elements suggest that infrequent episodes of marked genomic change followed by clonal expansion in geographically and enzootically structured populations may account for the unique lineage structure of Borreliaceae. This analysis informs future genotype-phenotype studies among Borreliaceae and lays a foundation for studies of individual gene function guided by phylogenetic patterns of conservation, diversification, gain, and/or loss.

Characterization of the family-level Borreliaceae pan-genome and development of an episomal typing protocol

Background

The Borreliaceae family includes many obligate parasitic bacterial species which are etiologically associated with a myriad of zoonotic borrelioses including Lyme disease and vector-borne relapsing fevers. Infections by the Borreliaceae are difficult to detect by both direct and indirect methods, often leading to delayed and missed diagnoses. Efforts to improve diagnoses center around the development of molecular diagnostics (MDx), but due to deep tissue sequestration of the causative spirochaetes and the lack of persistent bacteremias, even MDx assays suffer from a lack of sensitivity. Additionally, the highly extensive genomic heterogeneity among isolates, even within the same species, contributes to the lack of assay sensitivity as single target assays cannot provide universal coverage. This within-species heterogeneity is partly due to differences in replicon repertoires and genomic structures that have likely arisen to support the complex Borreliaceae lifecycle in which these parasites have to survive in multiple hosts each with unique immune responses.

Results

We constructed a Borreliaceae family-level pangenome and characterized the phylogenetic relationships among the constituent taxa which supports the recent taxonomy of splitting the family into at least two genera. Gene content pro les were created for the majority of the Borreliaceae replicons, providing for the first time their unambiguous molecular typing.

Conclusion

Our characterization of the Borreliaceae pan-genome supports the splitting of the former Borrelia genus into two genera and provides for the phylogenetic placement of several non-species designated isolates. Mining this family-level pangenome will enable precision diagnostics corresponding to gene content-driven clinical outcomes while also providing targets for interventions.

MCP5, a methyl-accepting chemotaxis protein regulated by both the Hk1-Rrp1 and Rrp2-RpoN-RpoS pathways, is required for the immune evasion of Borrelia burgdorferi

Borrelia (or Borreliella) burgdorferi, the causative agent of Lyme disease, is a motile and invasive zoonotic pathogen, adept at navigating between its arthropod vector and mammalian host. While motility and chemotaxis are well established as essential for its enzootic cycle, the function of methyl-accepting chemotaxis proteins (MCPs) in the infectious cycle of B. burgdorferi remains unclear. In this study, we demonstrate that MCP5, one of the most abundant MCPs in B. burgdorferi, is differentially expressed in response to environmental signals as well as at different stages of the pathogen's enzootic cycle. Specifically, the expression of mcp5 is regulated by the Hk1-Rrp1 and Rrp2-RpoN-RpoS pathways, which are critical for the spirochete's colonization of the tick vector and mammalian host, respectively. Infection experiments with an mcp5 mutant revealed that spirochetes lacking MCP5 could not establish infections in either C3H/HeN mice or Severe Combined Immunodeficiency (SCID) mice, which are defective in adaptive immunity, indicating the essential role of MCP5 in mammalian infection. However, the mcp5 mutant could establish infection and disseminate in NOD SCID Gamma (NSG) mice, which are deficient in both adaptive and most innate immune responses, suggesting a crucial role of MCP5 in evading host innate immunity. In the tick vector, the mcp5 mutants survived feeding but failed to transmit to mice, highlighting the importance of MCP5 in transmission. Our findings reveal that MCP5, regulated by the Rrp1 and Rrp2 pathways, is critical for the establishment of infection in mammalian hosts by evading host innate immunity and is important for the transmission of spirochetes from ticks to mammalian hosts, underscoring its potential as a target for intervention strategies.

Unguarded liabilities: Borrelia burgdorferi's complex amino acid dependence exposes unique avenues of inhibition

Recent reports from the Centers for Disease Control and Prevention approximate 500,000 cases of Lyme disease in the United States yearly, a significant economic burden on the healthcare system. The standard treatment for Lyme disease includes broad-spectrum antibiotics, which may be administered for extensive periods of time and result in significant impacts to the patient. Recently, we demonstrated that Borrelia burgdorferi, the causative agent of Lyme disease, is uniquely dependent upon peptide acquisition via an oligopeptide transport (Opp) system. This dependence appears unique to the spirochete; thus, the Opp system may constitute a novel and specific target for the inhibition of B. burgdorferi. For proof of concept, we conducted a pilot screen to determine if the Opp system constitutes a viable inhibitor target. OppA2 was utilized as our target protein as it is the most prolific peptide-binding protein throughout the enzootic cycle. We validated a thermal shift assay (TSA) to detect ligand binding against OppA2 and performed a high-throughput screen of 2,240 molecules from a diversity set library. The TSA results identified eight compounds (C1-8) demonstrating potential binding to OppA2, and growth assays identified C2 and C7 as inhibitors of B. burgdorferi growth. We confirmed by TSA that these two compounds interact with additional B. burgdorferi OppAs, potentially resulting in a cumulative inhibitory effect. Additionally, we showed that these compounds have no effect on Escherichia coli, a bacterium that encodes a dispensable Opp system which serves only as an ancillary nutrient transporter. These data demonstrate that the Opp system of B. burgdorferi acts as a viable drug target, with the potential for targeting multiple OppAs with a single compound. Moreover, the lack of inhibition against E. coli suggests that selective targeting of B. burgdorferi via the Opp system may be possible.

Borrelia burgdorferi PlzA is a cyclic-di-GMP dependent DNA and RNA binding protein

The PilZ domain-containing protein, PlzA, is the only known cyclic di-GMP binding protein encoded by all Lyme disease spirochetes. PlzA has been implicated in the regulation of many borrelial processes, but the effector mechanism of PlzA was not previously known. Here, we report that PlzA can bind DNA and RNA and that nucleic acid binding requires c-di-GMP, with the affinity of PlzA for nucleic acids increasing as concentrations of c-di-GMP were increased. A mutant PlzA that is incapable of binding c-di-GMP did not bind to any tested nucleic acids. We also determined that PlzA interacts predominantly with the major groove of DNA and that sequence length and G-C content play a role in DNA binding affinity. PlzA is a dual-domain protein with a PilZ-like N-terminal domain linked to a canonical C-terminal PilZ domain. Dissection of the domains demonstrated that the separated N-terminal domain bound nucleic acids independently of c-di-GMP. The C-terminal domain, which includes the c-di-GMP binding motifs, did not bind nucleic acids under any tested conditions. Our data are supported by computational docking, which predicts that c-di-GMP binding at the C-terminal domain stabilizes the overall protein structure and facilitates PlzA-DNA interactions via residues in the N-terminal domain. Based on our data, we propose that levels of c-di-GMP during the various stages of the enzootic life cycle direct PlzA binding to regulatory targets.

Members of the paralogous gene family 12 from the Lyme disease agent Borrelia burgdorferi are non-specific DNA-binding proteins

Lyme disease is the most prevalent vector-borne infectious disease in Europe and the USA. Borrelia burgdorferi, as the causative agent of Lyme disease, is transmitted to the mammalian host during the tick blood meal. To adapt to the different encountered environments, Borrelia has adjusted the expression pattern of various, mostly outer surface proteins. The function of most B. burgdorferi outer surface proteins remains unknown. We determined the crystal structure of a previously uncharacterized B. burgdorferi outer surface protein BBK01, known to belong to the paralogous gene family 12 (PFam12) as one of its five members. PFam12 members are shown to be upregulated as the tick starts its blood meal. Structural analysis of BBK01 revealed similarity to the coiled coil domain of structural maintenance of chromosomes (SMC) protein family members, while functional studies indicated that all PFam12 members are non-specific DNA-binding proteins. The residues involved in DNA binding were identified and probed by site-directed mutagenesis. The combination of SMC-like proteins being attached to the outer membrane and exposed to the environment or located in the periplasm, as observed in the case of PFam12 members, and displaying the ability to bind DNA, represents a unique feature previously not observed in bacteria.

Characterization and genomic analysis of the Lyme disease spirochete bacteriophage ϕBB-1

Lyme disease is a tick-borne infection caused by the spirochete Borrelia (Borreliella) burgdorferi. Borrelia species have highly fragmented genomes composed of a linear chromosome and a constellation of linear and circular plasmids some of which are required throughout the enzootic cycle. Included in this plasmid repertoire by almost all Lyme disease spirochetes are the 32-kb circular plasmid cp32 prophages that are capable of lytic replication to produce infectious virions called ϕBB-1. While the B. burgdorferi genome contains evidence of horizontal transfer, the mechanisms of gene transfer between strains remain unclear. While we know that ϕBB-1 transduces cp32 and shuttle vector DNA during in vitro cultivation, the extent of ϕBB-1 DNA transfer is not clear. Herein, we use proteomics and long-read sequencing to further characterize ϕBB-1 virions. Our studies identified the cp32 pac region and revealed that ϕBB-1 packages linear cp32s via a headful mechanism with preferential packaging of plasmids containing the cp32 pac region. Additionally, we find ϕBB-1 packages fragments of the linear chromosome and full-length plasmids including lp54, cp26, and others. Furthermore, sequencing of ϕBB-1 packaged DNA allowed us to resolve the covalently closed hairpin telomeres for the linear B. burgdorferi chromosome and most linear plasmids in strain CA-11.2A. Collectively, our results shed light on the biology of the ubiquitous ϕBB-1 phage and further implicates ϕBB-1 in the generalized transduction of diverse genes and the maintenance of genetic diversity in Lyme disease spirochetes.

BadR directly represses the expression of the glycerol utilization operon in the Lyme disease pathogen

Glycerol utilization as a carbohydrate source by Borreliella burgdorferi, the Lyme disease spirochete, is critical for its successful colonization and persistence in the tick vector. The expression of the glpFKD (glp) operon, which encodes proteins for glycerol uptake/utilization, must be tightly regulated during the enzootic cycle of B. burgdorferi. Previous studies have established that the second messenger cyclic di-GMP (c-di-GMP) is required for the activation of glp expression, while an alternative sigma factor RpoS acts as a negative regulator for glp expression. In the present study, we report identification of a cis element within the 5´ untranslated region of glp that exerts negative regulation of glp expression. Further genetic screen of known and predicted DNA-binding proteins encoded in the genome of B. burgdorferi uncovered that overexpressing Borrelia host adaptation regulator (BadR), a known global regulator, dramatically reduced glp expression. Similarly, the badR mutant significantly increased glp expression. Subsequent electrophoretic mobility shift assay analyses demonstrated that BadR directly binds to this cis element, thereby repressing glp independent of RpoS-mediated repression. The efficiency of BadR binding was further assessed in the presence of c-di-GMP and various carbohydrates. This finding highlights multi-layered positive and negative regulatory mechanisms employed by B. burgdorferi to synchronize glp expression throughout its enzootic cycle.IMPORTANCEBorreliella burgdorferi, the Lyme disease pathogen, must modulate its gene expression differentially to adapt successfully to its two disparate hosts. Previous studies have demonstrated that the glycerol uptake and utilization operon, glpFKD, plays a crucial role in spirochetal survival within ticks. However, the glpFKD expression must be repressed when B. burgdorferi transitions to the mammalian host. In this study, we identified a specific cis element responsible for the repression of glpFKD. We further pinpointed Borrelia host adaptation regulator as the direct binding protein to this cis element, thereby repressing glpFKD expression. This discovery paves the way for a deeper exploration of how zoonotic pathogens sense distinct hosts and switch their carbon source utilization during transmission.

Human Tick-Borne Diseases and Advances in Anti-Tick Vaccine Approaches: A Comprehensive Review

This comprehensive review explores the field of anti-tick vaccines, addressing their significance in combating tick-borne diseases of public health concern. The main objectives are to provide a brief epidemiology of diseases affecting humans and a thorough understanding of tick biology, traditional tick control methods, the development and mechanisms of anti-tick vaccines, their efficacy in field applications, associated challenges, and future prospects. Tick-borne diseases (TBDs) pose a significant and escalating threat to global health and the livestock industries due to the widespread distribution of ticks and the multitude of pathogens they transmit. Traditional tick control methods, such as acaricides and repellents, have limitations, including environmental concerns and the emergence of tick resistance. Anti-tick vaccines offer a promising alternative by targeting specific tick proteins crucial for feeding and pathogen transmission. Developing vaccines with antigens based on these essential proteins is likely to disrupt these processes. Indeed, anti-tick vaccines have shown efficacy in laboratory and field trials successfully implemented in livestock, reducing the prevalence of TBDs. However, some challenges still remain, including vaccine efficacy on different hosts, polymorphisms in ticks of the same species, and the economic considerations of adopting large-scale vaccine strategies. Emerging technologies and approaches hold promise for improving anti-tick vaccine development and expanding their impact on public health and agriculture.

Borrelia burgdorferi PlzA is a cyclic-di-GMP dependent DNA and RNA binding protein

The PilZ domain-containing protein, PlzA, is the only known cyclic di-GMP binding protein encoded by all Lyme disease spirochetes. PlzA has been implicated in the regulation of many borrelial processes, but the effector mechanism of PlzA was not previously known. Here we report that PlzA can bind DNA and RNA and that nucleic acid binding requires c-di-GMP, with the affinity of PlzA for nucleic acids increasing as concentrations of c-di-GMP were increased. A mutant PlzA that is incapable of binding c-di-GMP did not bind to any tested nucleic acids. We also determined that PlzA interacts predominantly with the major groove of DNA and that sequence length plays a role in DNA binding affinity. PlzA is a dual-domain protein with a PilZ-like N-terminal domain linked to a canonical C-terminal PilZ domain. Dissection of the domains demonstrated that the separated N-terminal domain bound nucleic acids independently of c-di-GMP. The C-terminal domain, which includes the c-di-GMP binding motifs, did not bind nucleic acids under any tested conditions. Our data are supported by computational docking, which predicts that c-di-GMP binding at the C-terminal domain stabilizes the overall protein structure and facilitates PlzA-DNA interactions via residues in the N-terminal domain. Based on our data, we propose that levels of c-di-GMP during the various stages of the enzootic life cycle direct PlzA binding to regulatory targets.

Characterization and genomic analysis of the Lyme disease spirochete bacteriophage ϕBB-1

Lyme disease is a tick-borne infection caused by the spirochete Borrelia (Borreliella) burgdorferi. Borrelia species have highly fragmented genomes composed of a linear chromosome and a constellation of linear and circular plasmids some of which are required throughout the enzootic cycle. Included in this plasmid repertoire by almost all Lyme disease spirochetes are the 32-kb circular plasmid cp32 prophages that are capable of lytic replication to produce infectious virions called ϕBB-1. While the B. burgdorferi genome contains evidence of horizontal transfer, the mechanisms of gene transfer between strains remain unclear. While we know that ϕBB-1 transduces cp32 and shuttle vector DNA during in vitro cultivation, the extent of ϕBB-1 DNA transfer is not clear. Herein, we use proteomics and long-read sequencing to further characterize ϕBB-1 virions. Our studies identified the cp32 pac region and revealed that ϕBB-1 packages linear cp32s via a headful mechanism with preferentially packaging of plasmids containing the cp32 pac region. Additionally, we find ϕBB-1 packages fragments of the linear chromosome and full-length plasmids including lp54, cp26, and others. Furthermore, sequencing of ϕBB-1 packaged DNA allowed us to resolve the covalently closed hairpin telomeres for the linear B. burgdorferi chromosome and most linear plasmids in strain CA-11.2A. Collectively, our results shed light on the biology of the ubiquitous ϕBB-1 phage and further implicates ϕBB-1 in the generalized transduction of diverse genes and the maintenance of genetic diversity in Lyme disease spirochetes.

Developing a Prospective Gestational Lyme Disease Study

Lyme disease in pregnancy is understudied. The few available reports of Borrelia infection during pregnancy collecting clinical outcomes, with or without confirmed fetal infection both in utero and neonatal, are limited to case reports and small series. Population-based studies are not available. We propose a prospective study of Borrelia infection during pregnancy based in obstetrical practices in both endemic and nonendemic areas, with long term follow-up of pregnancy outcomes and development assessment of offspring infected or exposed to Borrelia in utero using current serological, microscopic, culture, and molecular techniques. In addition to detection of Borrelia burgdorferi sensu stricto, additional Borrelia species and other pathogens known to be transmitted by ticks will be tested. Serial biospecimens including maternal and cord blood, maternal peripheral blood mononuclear cells and urine, and, when clinically indicated, amniotic fluid, chorionic villi, intrauterine cord blood, will be collected with clinical data, imaging, and for infections treatment medications. Offspring will be followed until age 5 years with annual developmental assessments to assess pregnancy outcomes. The study will require parallel development of a biorepository with strategies for management, data security and data sharing. A public-private partnership will be required to support the study.

Genome Rearrangement Analysis : Cut and Join Genome Rearrangements and Gene Cluster Preserving Approaches

Genome rearrangements are mutations that change the gene content of a genome or the arrangement of the genes on a genome. Several years of research on genome rearrangements have established different algorithmic approaches for solving some fundamental problems in comparative genomics based on gene order information. This review summarizes the literature on genome rearrangement analysis along two lines of research. The first line considers rearrangement models that are particularly well suited for a theoretical analysis. These models use rearrangement operations that cut chromosomes into fragments and then join the fragments into new chromosomes. The second line works with rearrangement models that reflect several biologically motivated constraints, e.g., the constraint that gene clusters have to be preserved. In this chapter, the border between algorithmically "easy" and "hard" rearrangement problems is sketched and a brief review is given on the available software tools for genome rearrangement analysis.

Structural studies of chromosomally encoded outer surface lipoprotein BB0158 from Borrelia burgdorferi sensu stricto

Lyme disease, or also known as Lyme borreliosis, is caused by the spirochetes belonging to the Borrelia burgdorferi sensu lato complex, which can enter the human body following the bite of an infected tick. Many membrane lipid-bound proteins, also known as lipoproteins, are located on the surface of B. burgdorferi sensu lato and play a crucial role in the spirochete to interact with its environment, whether in ticks or mammals. Since the spirochete needs to perform various tasks, such as resisting the host's immune system or spreading throughout the organism, it is not surprising that numerous surface proteins have been found to be essential for B. burgdorferi sensu lato complex bacteria in causing Lyme disease. In this study, we have determined (at 2.4 Å resolution) and characterized the 3D structure of BB0158, one of the few chromosomally encoded outer surface proteins from B. burgdorferi sensu stricto. BB0158 belongs to the paralogous gene family 44 (PFam44), consisting of four other members (BB0159, BBA04, BBE09 and BBK52). The characterization of BB0158, which appears to form a domain-swapped dimer, in conjunction with the characterization of the corresponding PFam44 members, certainly contribute to our understanding of B. burgdorferi sensu stricto proteins.

Role of bacteriophages in shaping gut microbial community

Phages are the most diversified and dominant members of the gut virobiota. They play a crucial role in shaping the structure and function of the gut microbial community and consequently the health of humans and animals. Phages are found mainly in the mucus, from where they can translocate to the intestinal organs and act as a modulator of gut microbiota. Understanding the vital role of phages in regulating the composition of intestinal microbiota and influencing human and animal health is an emerging area of research. The relevance of phages in the gut ecosystem is supported by substantial evidence, but the importance of phages in shaping the gut microbiota remains unclear. Although information regarding general phage ecology and development has accumulated, detailed knowledge on phage-gut microbe and phage-human interactions is lacking, and the information on the effects of phage therapy in humans remains ambiguous. In this review, we systematically assess the existing data on the structure and ecology of phages in the human and animal gut environments, their development, possible interaction, and subsequent impact on the gut ecosystem dynamics. We discuss the potential mechanisms of prophage activation and the subsequent modulation of gut bacteria. We also review the link between phages and the immune system to collect evidence on the effect of phages on shaping the gut microbial composition. Our review will improve understanding on the influence of phages in regulating the gut microbiota and the immune system and facilitate the development of phage-based therapies for maintaining a healthy and balanced gut microbiota.

Metabolic modeling predicts unique drug targets in Borrelia burgdorferi

IMPORTANCE

Lyme disease is often treated using long courses of antibiotics, which can cause side effects for patients and risks the evolution of antimicrobial resistance. Narrow-spectrum antimicrobials would reduce these risks, but their development has been slow because the Lyme disease bacterium, Borrelia burgdorferi, is difficult to work with in the laboratory. To accelerate the drug discovery pipeline, we developed a computational model of B. burgdorferi's metabolism and used it to predict essential enzymatic reactions whose inhibition prevented growth in silico. These predictions were validated using small-molecule enzyme inhibitors, several of which were shown to have specific activity against B. burgdorferi. Although the specific compounds used are not suitable for clinical use, we aim to use them as lead compounds to develop optimized drugs targeting the pathways discovered here.

Structure of the Borrelia Bacteriophage φBB1 Procapsid

Bacteriophages of Borrelia burgdorferi are a biologically important but under-investigated feature of the Lyme disease-causing spirochete. No virulent borrelial viruses have been identified, but all B. burgdorferi isolates carry a prophage φBB1 as resident circular plasmids. Like its host, the φBB1 phage is quite distinctive and shares little sequence similarity with other known bacteriophages. We expressed φBB1 head morphogenesis proteins in Escherichia coli which resulted in assembly of homogeneous prolate procapsid structures and used cryo-electron microscopy to determine the three-dimensional structure of these particles. The φBB1 procapsids consist of 415 copies of the major capsid protein and an equal combined number of three homologous capsid decoration proteins that form trimeric knobs on the outside of the particle. One of the end vertices of the particle is occupied by a portal assembled from twelve copies of the portal protein. The φBB1 scaffolding protein is entirely α-helical and has an elongated shape with a small globular domain in the middle. Within the tubular section of the procapsid, the internal scaffold is built of stacked rings, each composed of 32 scaffolding protein molecules, which run in opposite directions from both caps with a heterogeneous part in the middle. Inside the portal-containing cap, the scaffold is organized asymmetrically with ten scaffolding protein molecules bound to the portal. The φBB1 procapsid structure provides better insight into the vast structural diversity of bacteriophages and presents clues of how elongated bacteriophage particles might be assembled.

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.

Comparative genomics analysis of three conserved plasmid families in the Western Hemisphere soft tick-borne relapsing fever borreliae provides insight into variation in genome structure and antigenic variation systems

Borrelia spirochetes, causative agents of Lyme disease and relapsing fever (RF), have uniquely complex genomes, consisting of a linear chromosome and both circular and linear plasmids. The plasmids harbor genes important for the vector-host life cycle of these tick-borne bacteria. The role of plasmids from Lyme disease causing spirochetes is more refined compared to RF Borrelia because of limited plasmid-resolved genome assemblies for the latter. We recently addressed this limitation and found that three linear plasmid families (F6, F27, and F28) were syntenic across all the RF Borrelia species that we examined. Given this conservation, we further investigated the three plasmid families. The F6 family, also known as the megaplasmid, contained regions of repetitive DNA. The F27 was the smallest, encoding genes with unknown function. The F28 family encoded the putative expression locus for antigenic variation in all species except Borrelia hermsii and Borrelia anserina. Taken together, this work provides a foundation for future investigations to identify essential plasmid-localized genes that drive the vector-host life cycle of RF Borrelia. IMPORTANCE Borrelia spp. spirochetes are arthropod-borne bacteria found globally that infect humans and other vertebrates. RF borreliae are understudied and misdiagnosed pathogens because of the vague clinical presentation of disease and the elusive feeding behavior of argasid ticks. Consequently, genomics resources for RF spirochetes have been limited. Analyses of Borrelia plasmids have been challenging because they are often highly fragmented and unassembled in most available genome assemblies. By utilizing Oxford Nanopore Technologies, we recently generated plasmid-resolved genome assemblies for seven Borrelia spp. found in the Western Hemisphere. This current study is an in-depth investigation into the linear plasmids that were conserved and syntenic across species. We identified differences in genome structure and, importantly, in antigenic variation systems between species. This work is an important step in identifying crucial plasmid-localized genetic elements essential for the life cycle of RF spirochetes.

Mobile Genetic Element Flexibility as an Underlying Principle to Bacterial Evolution

Mobile genetic elements are key to the evolution of bacteria and traits that affect host and ecosystem health. Here, we use a framework of a hierarchical and modular system that scales from genes to populations to synthesize recent findings on mobile genetic elements (MGEs) of bacteria. Doing so highlights the role that emergent properties of flexibility, robustness, and genetic capacitance of MGEs have on the evolution of bacteria. Some of their traits can be stored, shared, and diversified across different MGEs, taxa of bacteria, and time. Collectively, these properties contribute to maintaining functionality against perturbations while allowing changes to accumulate in order to diversify and give rise to new traits. These properties of MGEs have long challenged our abilities to study them. Implementation of new technologies and strategies allows for MGEs to be analyzed in new and powerful ways.

Lyme borreliosis diagnosis: state of the art of improvements and innovations

With almost 700 000 estimated cases each year in the United States and Europe, Lyme borreliosis (LB), also called Lyme disease, is the most common tick-borne illness in the world. Transmitted by ticks of the genus Ixodes and caused by bacteria Borrelia burgdorferi sensu lato, LB occurs with various symptoms, such as erythema migrans, which is characteristic, whereas others involve blurred clinical features such as fatigue, headaches, arthralgia, and myalgia. The diagnosis of Lyme borreliosis, based on a standard two-tiered serology, is the subject of many debates and controversies, since it relies on an indirect approach which suffers from a low sensitivity depending on the stage of the disease. Above all, early detection of the disease raises some issues. Inappropriate diagnosis of Lyme borreliosis leads to therapeutic wandering, inducing potential chronic infection with a strong antibody response that fails to clear the infection. Early and proper detection of Lyme disease is essential to propose an adequate treatment to patients and avoid the persistence of the pathogen. This review presents the available tests, with an emphasis on the improvements of the current diagnosis, the innovative methods and ideas which, ultimately, will allow more precise detection of LB.

The Lyme disease spirochete, Borrelia burgdorferi, as a model vector-borne pathogen: insights on regulation of gene and protein expression

The Lyme disease spirochete persists in nature through cycles between ticks and vertebrates. Although the spirochete interacts with numerous, distinct tissues and environmental conditions during its infectious cycle, Borrelia burgdorferi appears to possess a limited ability to sense its external environment. This apparent paradox is being resolved through detailed investigations of the molecular mechanisms through which B. burgdorferi controls production of virulence-associated factors such as the Erp outer surface proteins. The results have led to development of a model for how B. burgdorferi controls expression of its diverse proteins, wherein physiological and metabolic states that are unique to specific points in the infectious cycle trigger changes in gene and protein expression levels.

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.

A high fidelity approach to assembling the complex Borrelia genome

BACKGROUND

Bacteria of the Borrelia burgdorferi sensu lato (s.l.) complex can cause Lyme borreliosis. Different B. burgdorferi s.l. genospecies vary in their host and vector associations and human pathogenicity but the genetic basis for these adaptations is unresolved and requires completed and reliable genomes for comparative analyses. The de novo assembly of a complete Borrelia genome is challenging due to the high levels of complexity, represented by a high number of circular and linear plasmids that are dynamic, showing mosaic structure and sequence homology. Previous work demonstrated that even advanced approaches, such as a combination of short-read and long-read data, might lead to incomplete plasmid reconstruction. Here, using recently developed high-fidelity (HiFi) PacBio sequencing, we explored strategies to obtain gap-free, complete and high quality Borrelia genome assemblies. Optimizing genome assembly, quality control and refinement steps, we critically appraised existing techniques to create a workflow that lead to improved genome reconstruction.

RESULTS

Despite the latest available technologies, stand-alone sequencing and assembly methods are insufficient for the generation of complete and high quality Borrelia genome assemblies. We developed a workflow pipeline for the de novo genome assembly for Borrelia using several types of sequence data and incorporating multiple assemblers to recover the complete genome including both circular and linear plasmid sequences.

CONCLUSION

Our study demonstrates that, with HiFi data and an ensemble reconstruction pipeline with refinement steps, chromosomal and plasmid sequences can be fully resolved, even for complex genomes such as Borrelia. The presented pipeline may be of interest for the assembly of further complex microbial genomes.

Longitudinal map of transcriptome changes in the Lyme pathogen Borrelia burgdorferi during tick-borne transmission

Borrelia burgdorferi (Bb), the causative agent of Lyme disease, adapts to vastly different environments as it cycles between tick vector and vertebrate host. During a tick bloodmeal, Bb alters its gene expression to prepare for vertebrate infection; however, the full range of transcriptional changes that occur over several days inside of the tick are technically challenging to capture. We developed an experimental approach to enrich Bb cells to longitudinally define their global transcriptomic landscape inside nymphal Ixodes scapularis ticks during a transmitting bloodmeal. We identified 192 Bb genes that substantially change expression over the course of the bloodmeal from 1 to 4 days after host attachment. The majority of upregulated genes encode proteins found at the cell envelope or proteins of unknown function, including 45 outer surface lipoproteins embedded in the unusual protein-rich coat of Bb. As these proteins may facilitate Bb interactions with the host, we utilized mass spectrometry to identify candidate tick proteins that physically associate with Bb. The Bb enrichment methodology along with the ex vivo Bb transcriptomes and candidate tick interacting proteins presented here provide a resource to facilitate investigations into key determinants of Bb priming and transmission during the tick stage of its unique transmission cycle.

A mathematician's guide to plasmids: an introduction to plasmid biology for modellers

Plasmids, extrachromosomal DNA molecules commonly found in bacterial and archaeal cells, play an important role in bacterial genetics and evolution. Our understanding of plasmid biology has been furthered greatly by the development of mathematical models, and there are many questions about plasmids that models would be useful in answering. In this review, we present an introductory, yet comprehensive, overview of the biology of plasmids suitable for modellers unfamiliar with plasmids who want to get up to speed and to begin working on plasmid-related models. In addition to reviewing the diversity of plasmids and the genes they carry, their key physiological functions, and interactions between plasmid and host, we also highlight selected plasmid topics that may be of particular interest to modellers and areas where there is a particular need for theoretical development. The world of plasmids holds a great variety of subjects that will interest mathematical biologists, and introducing new modellers to the subject will help to expand the existing body of plasmid theory.