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

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.

SUMMARY

Lyme disease is the most commonly reported arthropod-borne illness in the US, Europe, and Asia. The causative agent of Lyme disease, Borrelia burgdorferi , is maintained in an enzootic cycle involving arthropod vectors ( Ixodes ticks) and rodent mammalian hosts. Understanding how B. burgdorferi moves within this natural cycle is crucial for developing new strategies to combat Lyme disease. The complex nature of the enzootic cycle necessitates sensory-guided movement in response to environmental stimuli. B. burgdorferi possesses a unique and intricate chemotaxis signaling system, with methyl-accepting chemotaxis proteins (MCPs) at its core. These proteins are responsible for sensing environmental signals and guiding bacterial movement toward or away from stimuli. This study found that one of the MCPs, MCP5, is highly expressed and differentially regulated during the enzootic cycle by the Hk1-Rrp1 and Rrp2-RpoN-RpoS pathways. MCP5 is crucial for mammalian infection, aiding in immune evasion and transmission from ticks to mammals, providing a foundation for further research into B. burgdorferi 's navigation within its hosts.

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.

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.

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.

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.

Organization and replicon interactions within the highly segmented genome of Borrelia burgdorferi

Borrelia burgdorferi, a causative agent of Lyme disease, contains the most segmented bacterial genome known to date, with one linear chromosome and over twenty plasmids. How this unusually complex genome is organized, and whether and how the different replicons interact are unclear. We recently demonstrated that B. burgdorferi is polyploid and that the copies of the chromosome and plasmids are regularly spaced in each cell, which is critical for faithful segregation of the genome to daughter cells. Regular spacing of the chromosome is controlled by two separate partitioning systems that involve the protein pairs ParA/ParZ and ParB/Smc. Here, using chromosome conformation capture (Hi-C), we characterized the organization of the B. burgdorferi genome and the interactions between the replicons. We uncovered that although the linear chromosome lacks contacts between the two replication arms, the two telomeres are in frequent contact. Moreover, several plasmids specifically interact with the chromosome oriC region, and a subset of plasmids interact with each other more than with others. We found that Smc and the Smc-like MksB protein mediate long-range interactions on the chromosome, but they minimally affect plasmid-chromosome or plasmid-plasmid interactions. Finally, we found that disruption of the two partition systems leads to chromosome restructuring, correlating with the mis-positioning of chromosome oriC. Altogether, this study revealed the conformation of a complex genome and analyzed the contribution of the partition systems and SMC family proteins to this organization. This work expands the understanding of the organization and maintenance of multipartite bacterial genomes.

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

Lyme disease, caused by an infection with the spirochete Borrelia burgdorferi, is the most common vector-borne disease in North America. B. burgdorferi strains harbor extensive genomic and proteomic variability and further comparison is key to understanding the spirochetes infectivity and biological impacts of identified sequence variants. To achieve this goal, both transcript and mass spectrometry (MS)-based proteomics was applied to assemble peptide datasets of laboratory strains B31, MM1, B31-ML23, infective isolates B31-5A4, B31-A3, and 297, and other public datasets, to provide a publicly available Borrelia PeptideAtlas http://www.peptideatlas.org/builds/borrelia/. Included is information on total proteome, secretome, and membrane proteome of these B. burgdorferi strains. Proteomic data collected from 35 different experiment datasets, with a total of 855 mass spectrometry runs, identified 76,936 distinct peptides at a 0.1% peptide false-discovery-rate, which map to 1,221 canonical proteins (924 core canonical and 297 noncore canonical) and covers 86% of the total base B31 proteome. The diverse proteomic information from multiple isolates with credible data presented by the Borrelia PeptideAtlas can be useful to pinpoint potential protein targets which are common to infective isolates and may be key in the infection process.

Strict Conservation yet Non-Essential Nature of Plasmid Gene bba40 in the Lyme Disease Spirochete Borrelia burgdorferi

The highly segmented genome of Borrelia burgdorferi, the tick-borne bacterium that causes Lyme disease, is composed of a linear chromosome and more than 20 co-existing endogenous plasmids. Many plasmid-borne genes are unique to B. burgdorferi and some have been shown to provide essential functions at discrete points of the infectious cycle between a tick vector and rodent host. In this study, we investigated the role of bba40, a highly conserved and differentially expressed gene on a ubiquitous linear plasmid of B. burgdorferi. In a prior genome-wide analysis, inactivation of bba40 by transposon insertion was linked with a noninfectious phenotype in mice, suggesting that conservation of the gene in the Lyme disease spirochete reflected a critical function of the encoded protein. To address this hypothesis, we moved the bba40::Tn allele into a similar wild-type background and compared the phenotypes of isogenic wild-type, mutant and complemented strains in vitro and throughout the in vivo mouse/tick infectious cycle. In contrast to the previous study, we identified no defect in the ability of the bba40 mutant to colonize the tick vector or murine host, or to be efficiently transmitted between them. We conclude that bba40 joins a growing list of unique, highly conserved, yet fully dispensable plasmid-borne genes of the Lyme disease spirochete. We infer that the experimental infectious cycle, while including the tick vector and murine host, lacks key selective forces imposed during the natural enzootic cycle. IMPORTANCE The key finding of this study contradicts our premise that the ubiquitous presence and strict sequence conservation of a unique gene in the Lyme disease spirochete, Borrelia burgdorferi, reflect a critical role in either the murine host or tick vector in which these bacteria are maintained in nature. Instead, the outcome of this investigation illustrates the inadequate nature of the experimental infectious cycle currently employed in the laboratory to fully model the enzootic cycle of the Lyme disease spirochete. This study also highlights the importance of complementation for accurate interpretation of mutant phenotypes in genetic studies of Borrelia burgdorferi.

Combined direct/indirect detection allows identification of DNA termini in diverse sequencing datasets and supports a multiple-initiation-site model for HIV plus-strand synthesis

Replication of genetic material involves the creation of characteristic termini. Determining these termini is important to refine our understanding of the mechanisms involved in maintaining the genomes of cellular organisms and viruses. Here we describe a computational approach combining direct and indirect readouts to detect termini from next-generation short-read sequencing. While a direct inference of termini can come from mapping the most prominent start positions of captured DNA fragments, this approach is insufficient in cases where the DNA termini are not captured, whether for biological or technical reasons. Thus, a complementary (indirect) approach to terminus detection can be applied, taking advantage of the imbalance in coverage between forward and reverse sequence reads near termini. A resulting metric ("strand bias") can be used to detect termini even where termini are naturally blocked from capture or ends are not captured during library preparation (e.g., in tagmentation-based protocols). Applying this analysis to datasets where known DNA termini are present, such as from linear double-stranded viral genomes, yielded distinct strand bias signals corresponding to these termini. To evaluate the potential to analyze a more complex situation, we applied the analysis to examine DNA termini present early after HIV infection in a cell culture model. We observed both the known termini expected based on standard models of HIV reverse transcription (the U5-right-end and U3-left-end termini) as well as a signal corresponding to a previously described additional initiation site for plus-strand synthesis (cPPT [central polypurine tract]). Interestingly, we also detected putative terminus signals at additional sites. The strongest of these are a set that share several characteristics with the previously characterized plus-strand initiation sites (the cPPT and 3' PPT [polypurine tract] sites): (i) an observed spike in directly captured cDNA ends, an indirect terminus signal evident in localized strand bias, (iii) a preference for location on the plus-strand, (iv) an upstream purine-rich motif, and (v) a decrease in terminus signal at late time points after infection. These characteristics are consistent in duplicate samples in two different genotypes (wild type and integrase-lacking HIV). The observation of distinct internal termini associated with multiple purine-rich regions raises a possibility that multiple internal initiations of plus-strand synthesis might contribute to HIV replication.

Gac Is a Transcriptional Repressor of the Lyme Disease Spirochete's OspC Virulence-Associated Surface Protein

The OspC outer-surface lipoprotein is essential for the Lyme disease spirochete's initial phase of vertebrate infection. Bacteria within the midguts of unfed ticks do not express OspC but produce high levels when ticks begin to ingest blood. Lyme disease spirochetes cease production of OspC within 1 to 2 weeks of vertebrate infection, and bacteria that fail to downregulate OspC are cleared by host antibodies. Thus, tight regulation of OspC levels is critical for survival of Lyme borreliae and, therefore, an attractive target for development of novel treatment strategies. Previous studies determined that a DNA region 5' of the ospC promoter, the ospC operator, is required for control of OspC production. Hypothesizing that the ospC operator may bind a regulatory factor, DNA affinity pulldown was performed and identified binding by the Gac protein. Gac is encoded by the C-terminal domain of the gyrA open reading frame from an internal promoter, ribosome-binding site, and initiation codon. Our analyses determined that Gac exhibits a greater affinity for ospC operator and promoter DNAs than for other tested borrelial sequences. In vitro and in vivo analyses demonstrated that Gac is a transcriptional repressor of ospC. These results constitute a substantial advance to our understanding of the mechanisms by which the Lyme disease spirochete controls production of OspC. IMPORTANCE Borrelia burgdorferi sensu lato requires its surface-exposed OspC protein in order to establish infection in humans and other vertebrate hosts. Bacteria that either do not produce OspC during transmission or fail to repress OspC after infection is established are rapidly cleared by the host. Herein, we identified a borrelial protein, Gac, that exhibits preferential affinity to the ospC promoter and 5' adjacent DNA. A combination of biochemical analyses and investigations of genetically manipulated bacteria demonstrated that Gac is a transcriptional repressor of ospC. This is a substantial advance toward understanding how the Lyme disease spirochete controls production of the essential OspC virulence factor and identifies a novel target for preventative and curative therapies.

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.

Organization and replicon interactions within the highly segmented genome of Borrelia burgdorferi

Borrelia burgdorferi , a causative agent of Lyme disease, contains the most segmented bacterial genome known to date, with one linear chromosome and over twenty plasmids. How this unusually complex genome is organized, and whether and how the different replicons interact are unclear. We recently demonstrated that B. burgdorferi is polyploid and that the copies of the chromosome and plasmids are regularly spaced in each cell, which is critical for faithful segregation of the genome to daughter cells. Regular spacing of the chromosome is controlled by two separate partitioning systems that involve the protein pairs ParA/ParZ and ParB/SMC. Here, using chromosome conformation capture (Hi-C), we characterized the organization of the B. burgdorferi genome and the interactions between the replicons. We uncovered that although the linear chromosome lacks contacts between the two replication arms, the two telomeres are in frequent contact. Moreover, several plasmids specifically interact with the chromosome oriC region, and a subset of plasmids interact with each other more than with others. We found that SMC and the SMC-like MksB protein mediate long-range interactions on the chromosome, but they minimally affect plasmid-chromosome or plasmid-plasmid interactions. Finally, we found that disruption of the two partition systems leads to chromosome restructuring, correlating with the mis-positioning of chromosome oriC . Altogether, this study revealed the conformation of a complex genome and analyzed the contribution of the partition systems and SMC family proteins to this organization. This work expands the understanding of the organization and maintenance of multipartite bacterial genomes.

Erp and Rev Adhesins of the Lyme Disease Spirochete's Ubiquitous cp32 Prophages Assist the Bacterium during Vertebrate Infection

Almost all spirochetes in the genus Borrelia (sensu lato) naturally contain multiple variants of closely related prophages. In the Lyme disease borreliae, these prophages are maintained as circular episomes that are called circular plasmid 32 kb (cp32s). The cp32s of Lyme agents are particularly unique in that they encode two distinct families of lipoproteins, namely, Erp and Rev, that are expressed on the bacterial outer surface during infection of vertebrate hosts. All identified functions of those outer surface proteins involve interactions between the spirochetes and host molecules, as follows: Erp proteins bind plasmin(ogen), laminin, glycosaminoglycans, and/or components of complement and Rev proteins bind fibronectin. Thus, cp32 prophages provide their bacterial hosts with surface proteins that can enhance infection processes, thereby facilitating their own survival. Horizontal transfer via bacteriophage particles increases the spread of beneficial alleles and creates diversity among Erp and Rev proteins.

Endogenous Linear Plasmids lp28-4 and lp25 Are Required for Infectivity and Restriction Protection in the Lyme Disease Spirochete Borrelia mayonii

Borrelia mayonii is a newly recognized causative agent of Lyme disease in the Upper Midwestern United States, with distinct clinical presentations compared to classical Lyme disease caused by other Lyme Borrelia species. However, little is known about the B. mayonii genetic determinants required for establishing infection or perpetuating disease in mammals. Extrachromosomal plasmids in Borrelia species often encode proteins necessary for infection and pathogenesis, and spontaneous loss of these plasmids can lead to the identification of virulence determinant genes. Here, we describe infection of Lyme disease-susceptible C3H mice with B. mayonii, and show bacterial dissemination and persistence in peripheral tissues. Loss of endogenous plasmids, including lp28-4, lp25, and lp36 correlated with reduced infectivity in mice. The apparent requirement for lp28-4 during murine infection suggests the presence of a novel virulence determinant, as this plasmid does not encode homologs of any known virulence determinant. We also describe transformation and stable maintenance of a self-replicating shuttle vector in B. mayonii, and show that loss of either lp25 or lp28-4 correlated with increased transformation competency. Finally, we demonstrate that linear plasmids lp25 and lp28-4 each encode functional restriction modification systems with distinct but partially overlapping target modification sequences, which likely accounts for the observed decrease in transformation efficiency when those plasmids are present. Taken together, this study describes a role for endogenous plasmids in mammalian infection and restriction protection in the Lyme disease spirochete Borrelia mayonii.

Pleomorphic Variants of Borreliella (syn. Borrelia) burgdorferi Express Evolutionary Distinct Transcriptomes

Borreliella (syn. Borrelia) burgdorferi is a spirochete bacterium that causes tick-borne Lyme disease. Along its lifecycle B. burgdorferi develops several pleomorphic forms with unclear biological and medical relevance. Surprisingly, these morphotypes have never been compared at the global transcriptome level. To fill this void, we grew B. burgdorferi spirochete, round body, bleb, and biofilm-dominated cultures and recovered their transcriptomes by RNAseq profiling. We found that round bodies share similar expression profiles with spirochetes, despite their morphological differences. This sharply contrasts to blebs and biofilms that showed unique transcriptomes, profoundly distinct from spirochetes and round bodies. To better characterize differentially expressed genes in non-spirochete morphotypes, we performed functional, positional, and evolutionary enrichment analyses. Our results suggest that spirochete to round body transition relies on the delicate regulation of a relatively small number of highly conserved genes, which are located on the main chromosome and involved in translation. In contrast, spirochete to bleb or biofilm transition includes substantial reshaping of transcription profiles towards plasmids-residing and evolutionary young genes, which originated in the ancestor of Borreliaceae. Despite their abundance the function of these Borreliaceae-specific genes is largely unknown. However, many known Lyme disease virulence genes implicated in immune evasion and tissue adhesion originated in this evolutionary period. Taken together, these regularities point to the possibility that bleb and biofilm morphotypes might be important in the dissemination and persistence of B. burgdorferi inside the mammalian host. On the other hand, they prioritize the large pool of unstudied Borreliaceae-specific genes for functional characterization because this subset likely contains undiscovered Lyme disease pathogenesis genes.

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 a uniquely complex genome consisting of a linear chromosome and circular and linear plasmids. The plasmids harbor genes important for the vector-host life cycle of these tick-borne bacteria. The role of Lyme disease causing Borrelia plasmids is more refined compared to RF spirochetes because of limited plasmid-resolved genomes for RF spirochetes. We recently addressed this limitation and found that three linear plasmid families (F6, F27, and F28) were syntenic across all species. 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 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 and 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. By utilizing Oxford Nanopore Technologies, we recently generated plasmid-resolved genomes for seven Borrelia spp. found in the Western Hemisphere. This current study is a more in-depth investigation into the linear plasmids that were conserved and syntenic across species. This analysis determined differences in genome structure and, importantly, in antigenic variation systems between species. This work is an important step in identifying crucial plasmid-borne genetic elements essential for the life cycle of RF spirochetes.

Whole genome sequencing of Borrelia burgdorferi isolates reveals linked clusters of plasmid-borne accessory genome elements associated with virulence

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 patient-derived B. burgdorferi sensu stricto ( Bbss ) isolates from patients in the Eastern and Midwestern US and Central Europe. We develop a WGS-based classification of Bbss isolates, confirm and extend the findings of previous single- and multi-locus typing systems, define the plasmid profiles of human-infectious Bbss isolates, annotate the core and strain-variable surface lipoproteome, and identify loci associated with disseminated infection. A core genome consisting of ∼800 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 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, are associated with increased rates of dissemination. OspC type A strains possess a unique constellation of strongly linked genetic changes 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. The patterns of OspC type A strains typify a broader paradigm across Bbss isolates, in which genetic structure is defined by correlated groups of strain-variable genes located predominantly on plasmids, particularly for expression of surface-exposed lipoproteins. These clusters of genes are inherited in blocks through strain-specific patterns of plasmid occupancy and are associated with the probability of invasive infection.

Inducible CRISPRi-Based Operon Silencing and Selective in Trans Gene Complementation in Borrelia burgdorferi

To accelerate genetic studies on the Lyme disease pathogen Borrelia burgdorferi, we developed an enhanced CRISPR interference (CRISPRi) approach for isopropyl-β-d-thiogalactopyranoside (IPTG)-inducible repression of specific B. burgdorferi genes. The entire system is encoded on a compact 11-kb shuttle vector plasmid that allows for inducible expression of both the sgRNA module and a nontoxic codon-optimized dCas9 protein. We validated this CRISPRi system by targeting the genes encoding OspA and OspB, abundant surface lipoproteins coexpressed by a single operon, and FlaB, the major subunit forming the periplasmic flagella. As in other systems, single guide RNAs (sgRNAs) complementary to the nontemplate strand were consistently effective in gene repression, with 4- to 994-fold reductions in targeted transcript levels and concomitant reductions of protein levels. Furthermore, we showed that ospAB knockdowns could be selectively complemented in trans for OspA expression via the insertion of CRISPRi-resistant, synonymously or nonsynonymously mutated protospacer adjacent motif (PAM*) ospA alleles into a unique site within the CRISPRi plasmid. Together, this establishes CRISPRi PAM* as a robust new genetic tool to simplify the study of B. burgdorferi genes, bypassing the need for gene disruptions by allelic exchange and avoiding rare codon toxicity from the heterologous expression of dCas9. IMPORTANCE Borrelia burgdorferi, the spirochetal bacterium causing Lyme disease, is a tick-borne pathogen of global importance. Here, we expand the genetic toolbox for studying B. burgdorferi physiology and pathogenesis by establishing a single plasmid-based, fully inducible, and nontoxic CRISPR interference (CRISPRi) system for transcriptional silencing of B. burgdorferi genes and operons. We also show that alleles of CRISPRi-targeted genes with mutated protospacer-adjacent motif (PAM*) sites are CRISPRi resistant and can be used for simultaneous in trans gene complementation. The CRISPRi PAM* system will streamline the study of essential Borrelia proteins and accelerate investigations into their structure-function relationships.

Borrelia miyamotoi: A Comprehensive Review

Borrelia miyamotoi is an emerging tick-borne pathogen in the Northern Hemisphere and is the causative agent of Borrelia miyamotoi disease (BMD). Borrelia miyamotoi is vectored by the same hard-bodied ticks as Lyme disease Borrelia, yet phylogenetically groups with relapsing fever Borrelia, and thus, has been uniquely labeled a hard tick-borne relapsing fever Borrelia. Burgeoning research has uncovered new aspects of B. miyamotoi in human patients, nature, and the lab. Of particular interest are novel findings on disease pathology, prevalence, diagnostic methods, ecological maintenance, transmission, and genetic characteristics. Herein, we review recent literature on B. miyamotoi, discuss how findings adapt to current Borrelia doctrines, and briefly consider what remains unknown about B. miyamotoi.

Spatial correlated incidence modeling with zero inflation

Disease mapping models have been popularly used to model disease incidence with spatial correlation. In disease mapping models, zero inflation is an important issue, which often occurs in disease incidence datasets with high proportions of zero disease count. It is originated from limited survey coverage or unadvanced testing equipment, which makes some regions have no observed patients. Then excessive zeros recorded in the disease incidence dataset would mess up the true distributions of disease incidence and lead to inaccurate estimates. To address this issue, a zero-inflated disease mapping model is developed in this work. In this model, a zero-inflated process using Bernoulli indicators is assumed to characterize whether the zero inflation occurs for each region. For regions without zero inflation, a coherent and generative disease mapping model is applied for mapping the spatially correlated disease incidence. Independent spatial random effects are incorporated in both processes to account for the spatial patterns of zero inflation and disease incidence. External covariates are also considered in both processes to better explain the disease count data. To estimate the model, a Markov chain Monte Carlo algorithm is proposed. We evaluate model performance via a variety of simulation experiments. Finally, a Lyme disease dataset of Virginia is analyzed to illustrate the application of the proposed model.

Coupled induction of prophage and virulence factors during tick transmission of the Lyme disease spirochete

The alternative sigma factor RpoS plays a central role in the critical host-adaptive response of the Lyme disease spirochete, Borrelia burgdorferi. We previously identified bbd18 as a negative regulator of RpoS but could not inactivate bbd18 in wild-type spirochetes. In the current study we employed an inducible bbd18 gene to demonstrate the essential nature of BBD18 for viability of wild-type spirochetes in vitro and at a unique point in vivo. Transcriptomic analyses of BBD18-depleted cells demonstrated global induction of RpoS-dependent genes prior to lysis, with the absolute requirement for BBD18, both in vitro and in vivo, circumvented by deletion of rpoS. The increased expression of plasmid prophage genes and the presence of phage particles in the supernatants of lysing cultures indicate that RpoS regulates phage lysis-lysogeny decisions. Through this work we identify a mechanistic link between endogenous prophages and the RpoS-dependent adaptive response of the Lyme disease spirochete.

A Brief History of Plasmids

In the late 1950s, a number of laboratories took up the study of plasmids once the discovery was made that extrachromosomal antibiotic resistance (R) factors are the responsible agents for the transmissibility of multiple antibiotic resistance among the enterobacteria. The use of incompatibility for the classification of plasmids is now widespread. It seems clear now on the basis of the limited studies to date that the number of incompatibility groups of plasmids will likely be extremely large when one includes plasmids obtained from bacteria that are normal inhabitants of poorly studied natural environments. The presence of both linear chromosomes and linear plasmids is now established for several Streptomyces species. One of the more fascinating developments in plasmid biology was the discovery of linear plasmids in the 1980s. A remarkable feature of the Ti plasmids of Agrobacterium tumefaciens is the presence of two DNA transfer systems. A definitive demonstration that plasmids consisted of duplex DNA came from interspecies conjugal transfer of plasmids followed by separation of plasmid DNA from chromosomal DNA by equilibrium buoyant density centrifugation. The formation of channels for DNA movement and the actual steps involved in DNA transport offer many opportunities for the discovery of proteins with novel activities and for establishing fundamentally new concepts of macromolecular interactions between DNA and specific proteins, membranes, and the peptidoglycan matrix.

Cas9-mediated endogenous plasmid loss in Borrelia burgdorferi

The spirochete Borrelia burgdorferi, which causes Lyme disease, has the most segmented genome among known bacteria. In addition to a linear chromosome, the B. burgdorferi genome contains over 20 linear and circular endogenous plasmids. While many of these plasmids are dispensable under in vitro culture conditions, they are maintained during the natural life cycle of the pathogen. Plasmid-encoded functions are required for colonization of the tick vector, transmission to the vertebrate host, and evasion of host immune defenses. Different Borrelia strains can vary substantially in the type of plasmids they carry. The gene composition within the same type of plasmid can also differ from strain to strain, impeding the inference of plasmid function from one strain to another. To facilitate the investigation of the role of specific B. burgdorferi plasmids, we developed a Cas9-based approach that targets a plasmid for removal. As a proof-of-principle, we showed that targeting wild-type Cas9 to several loci on the endogenous plasmids lp25 or lp28-1 of the B. burgdorferi type strain B31 results in sgRNA-specific plasmid loss even when homologous sequences (i.e., potential sequence donors for DNA recombination) are present nearby. Cas9 nickase versions, Cas9D10A or Cas9H840A, also cause plasmid loss, though not as robustly. Thus, sgRNA-directed Cas9 DNA cleavage provides a highly efficient way to eliminate B. burgdorferi endogenous plasmids that are non-essential in axenic culture.

Polyploidy, regular patterning of genome copies, and unusual control of DNA partitioning in the Lyme disease spirochete

Borrelia burgdorferi, the tick-transmitted spirochete agent of Lyme disease, has a highly segmented genome with a linear chromosome and various linear or circular plasmids. Here, by imaging several chromosomal loci and 16 distinct plasmids, we show that B. burgdorferi is polyploid during growth in culture and that the number of genome copies decreases during stationary phase. B. burgdorferi is also polyploid inside fed ticks and chromosome copies are regularly spaced along the spirochete's length in both growing cultures and ticks. This patterning involves the conserved DNA partitioning protein ParA whose localization is controlled by a potentially phage-derived protein, ParZ, instead of its usual partner ParB. ParZ binds its own coding region and acts as a centromere-binding protein. While ParA works with ParZ, ParB controls the localization of the condensin, SMC. Together, the ParA/ParZ and ParB/SMC pairs ensure faithful chromosome inheritance. Our findings underscore the plasticity of cellular functions, even those as fundamental as chromosome segregation.

Riboflavin salvage by Borrelia burgdorferi supports carbon metabolism and is essential for survival in the tick vector

The Lyme disease agent, Borrelia burgdorferi, harbors a significantly reduced genome and relies on the scavenging of critical nutrients from its tick and mammalian hosts for survival. Riboflavin salvage has been shown to be important for B. burgdorferi infection of mice, yet the contributions of riboflavin to B. burgdorferi metabolism and survival in the tick remain unknown. Using a targeted mass spectrometry approach, we confirmed the importance of bb0318, the putative ATPase component of an ABC-type riboflavin transporter, for riboflavin salvage and the production of FMN and FAD. This analysis further revealed that Δbb0318 B. burgdorferi displayed increased levels of glycerol 3-phosphate compared to the wild-type. The glycerol 3-phosphate dehydrogenase activity of GlpD was found to be FAD-dependent and the transcription and translation of glpD were significantly decreased in Δbb0318 B. burgdorferi. Finally, gene bb0318 was found to be important for maximal spirochete burden in unfed larvae and essential for survival in feeding ticks. Together, these data demonstrate the importance of riboflavin salvage for B. burgdorferi carbon metabolism and survival in ticks.

Borrelia burgdorferi, the Lyme disease spirochete, possesses genetically-encoded responses to doxycycline, but not to amoxicillin

Some species of bacteria respond to antibiotic stresses by altering their transcription profiles, in order to produce proteins that provide protection against the antibiotic. Understanding these compensatory mechanisms allows for informed treatment strategies, and could lead to the development of improved therapeutics. To this end, studies were performed to determine whether Borrelia burgdorferi, the spirochetal agent of Lyme disease, also exhibits genetically-encoded responses to the commonly prescribed antibiotics doxycycline and amoxicillin. After culturing for 24 h in a sublethal concentration of doxycycline, there were significant increases in a substantial number of transcripts for proteins that are involved with translation. In contrast, incubation with a sublethal concentration of amoxicillin did not lead to significant changes in levels of any bacterial transcript. We conclude that B. burgdorferi has a mechanism(s) that detects translational inhibition by doxycycline, and increases production of mRNAs for proteins involved with translation machinery in an attempt to compensate for that stress.

FtlA and FtlB Are Candidates for Inclusion in a Next-Generation Multiantigen Subunit Vaccine for Lyme Disease

Lyme disease (LD) is a tick-transmitted bacterial infection caused by Borreliella burgdorferi and other closely related species collectively referred to as the LD spirochetes. The LD spirochetes encode an uncharacterized family of proteins originally designated protein family twelve (PF12). In B. burgdorferi strain B31, PF12 consists of four plasmid-carried genes, encoding BBK01, BBG01, BBH37, and BBJ08. Henceforth, we designate the PF12 proteins family twelve lipoprotein (Ftl) A (FtlA) (BBK01), FtlB (BBG01), FtlC (BBH37), and FtlD (BBJ08). The goal of this study was to assess the potential utility of the Ftl proteins in subunit vaccine development. Immunoblot analyses of LD spirochete cell lysates demonstrated that one or more of the Ftl proteins are produced by most LD isolates during cultivation. The Ftl proteins were verified to be membrane associated, and nondenaturing PAGE revealed that FtlA, FtlB, and FtlD formed dimers, while FtlC formed hexamers. Analysis of serum samples from B. burgdorferi antibody (Ab)-positive client-owned dogs (n = 50) and horses (n = 90) revealed that a majority were anti-Ftl Ab positive. Abs to the Ftl proteins were detected in serum samples from laboratory-infected dogs out to 497 days postinfection. Anti-FtlA and FtlB antisera displayed potent complement-dependent Ab-mediated killing activity, and epitope localization revealed that the bactericidal epitopes reside within the N-terminal domain of the Ftl proteins. This study suggests that FtlA and FtlB are potential candidates for inclusion in a multivalent vaccine for LD.

Comparative genome analysis of Vagococcus fluvialis reveals abundance of mobile genetic elements in sponge-isolated strains

BACKGROUND

Vagococcus fluvialis is a species of lactic acid bacteria found both free-living in river and seawater and associated to hosts, such as marine sponges. This species has been greatly understudied, with no complete genome assembly available to date, which is essential for the characterisation of the mobilome.

RESULTS

We sequenced and assembled de novo the complete genome sequences of five V. fluvialis isolates recovered from marine sponges. Pangenome analysis of the V. fluvialis species (total of 17 genomes) showed a high intraspecific diversity, with 45.5% of orthologous genes found to be strain specific. Despite this diversity, analyses of gene functions clustered all V. fluvialis species together and separated them from other sequenced Vagococcus species. V. fluvialis strains from different habitats were highly similar in terms of functional diversity but the sponge-isolated strains were enriched in several functions related to the marine environment. Furthermore, sponge-isolated strains carried a significantly higher number of mobile genetic elements (MGEs) compared to previously sequenced V. fluvialis strains from other environments. Sponge-isolated strains carried up to 4 circular plasmids each, including a 48-kb conjugative plasmid. Three of the five strains carried an additional circular extrachromosomal sequence, assumed to be an excised prophage as it contained mainly viral genes and lacked plasmid replication genes. Insertion sequences (ISs) were up to five times more abundant in the genomes of sponge-isolated strains compared to the others, including several IS families found exclusively in these genomes.

CONCLUSIONS

Our findings highlight the dynamics and plasticity of the V. fluvialis genome. The abundance of mobile genetic elements in the genomes of sponge-isolated V. fluvialis strains suggests that the mobilome might be key to understanding the genomic signatures of symbiosis in bacteria.

Identifying platelet-derived factors as amplifiers of B. burgdorferi-induced cytokine production

Previous studies have shown that monocytes can be 'trained' or tolerized by certain stimuli to respond stronger or weaker to a secondary stimulation. Rewiring of glucose metabolism was found to be important in inducing this phenotype. As we previously found that Borrelia burgdorferi (B. burgdorferi), the causative agent of Lyme borreliosis (LB), alters glucose metabolism in monocytes, we hypothesized that this may also induce long-term changes in innate immune responses. We found that exposure to B. burgdorferi decreased cytokine production in response to the TLR4-ligand lipopolysaccharide (LPS). In addition, B. burgdorferi exposure decreased baseline levels of glycolysis, as assessed by lactate production. Using GWAS analysis, we identified a gene, microfibril-associated protein 3-like (MFAP3L) as a factor influencing lactate production after B. burgdorferi exposure. Validation experiments proved that MFAP3L affects lactate- and cytokine production following B. burgdorferi stimulation. This is mediated by functions of MFAP3L, which includes activating ERK2 and through activation of platelet degranulation. Moreover, we showed that platelets and platelet-derived factors play important roles in B. burgdorferi-induced cytokine production. Certain platelet-derived factors, such chemokine C-X-C motif ligand 7 (CXCL7) and (C-C motif) ligand 5 (CCL5), were elevated in the circulation of LB patients in comparison to healthy individuals.

Comparative genomics of the Western Hemisphere soft tick-borne relapsing fever borreliae highlights extensive plasmid diversity

BACKGROUND

Tick-borne relapsing fever (TBRF) is a globally prevalent, yet under-studied vector-borne disease transmitted by soft and hard bodied ticks. While soft TBRF (sTBRF) spirochetes have been described for over a century, our understanding of the molecular mechanisms facilitating vector and host adaptation is poorly understood. This is due to the complexity of their small (~ 1.5 Mb) but fragmented genomes that typically consist of a linear chromosome and both linear and circular plasmids. A majority of sTBRF spirochete genomes' plasmid sequences are either missing or are deposited as unassembled sequences. Consequently, our goal was to generate complete, plasmid-resolved genomes for a comparative analysis of sTBRF species of the Western Hemisphere.

RESULTS

Utilizing a Borrelia specific pipeline, genomes of sTBRF spirochetes from the Western Hemisphere were sequenced and assembled using a combination of short- and long-read sequencing technologies. Included in the analysis were the two recently isolated species from Central and South America, Borrelia puertoricensis n. sp. and Borrelia venezuelensis, respectively. Plasmid analyses identified diverse sequences that clustered plasmids into 30 families; however, only three families were conserved and syntenic across all species. We also compared two species, B. venezuelensis and Borrelia turicatae, which were isolated ~ 6,800 km apart and from different tick vector species but were previously reported to be genetically similar.

CONCLUSIONS

To truly understand the biological differences observed between species of TBRF spirochetes, complete chromosome and plasmid sequences are needed. This comparative genomic analysis highlights high chromosomal synteny across the species yet diverse plasmid composition. This was particularly true for B. turicatae and B. venezuelensis, which had high average nucleotide identity yet extensive plasmid diversity. These findings are foundational for future endeavors to evaluate the role of plasmids in vector and host adaptation.

The Consistent Tick-Vertebrate Infectious Cycle of the Lyme Disease Spirochete Enables Borrelia burgdorferi To Control Protein Expression by Monitoring Its Physiological Status

The Lyme disease spirochete, Borrelia burgdorferi, persists in nature by alternatingly cycling between ticks and vertebrates. During each stage of the infectious cycle, B. burgdorferi produces surface proteins that are necessary for interactions with the tick or vertebrate tissues it encounters while also repressing the synthesis of unnecessary proteins. Among these are the Erp surface proteins, which are produced during vertebrate infection for interactions with host plasmin, laminin, glycosaminoglycans, and components of the complement system. Erp proteins are not expressed during tick colonization but are induced when the tick begins to ingest blood from a vertebrate host, a time when the bacteria undergo rapid growth and division. Using the erp genes as a model of borrelial gene regulation, our research group has identified three novel DNA-binding proteins that interact with DNA to control erp transcription. At least two of those regulators are, in turn, affected by DnaA, the master regulator of chromosome replication. Our data indicate that B. burgdorferi has evolved to detect the change from slow to rapid replication during tick feeding as a signal to begin expression of Erp and other vertebrate-specific proteins. The majority of other known regulatory factors of B. burgdorferi also respond to metabolic cues. These observations lead to a model in which the Lyme spirochete recognizes unique environmental conditions encountered during the infectious cycle to "know" where they are and adapt accordingly.

Controlled Proteolysis of an Essential Virulence Determinant Dictates Infectivity of Lyme Disease Pathogens

The Borrelia burgdorferi BB0323 protein undergoes a complex yet poorly defined proteolytic maturation event that generates N-terminal and C-terminal proteins with essential functions in cell growth and infection. Here, we report that a borrelial protease, B. burgdorferi high temperature requirement A protease (BbHtrA), cleaves BB0323 between asparagine (N) and leucine (L) at positions 236 and 237, while the replacement of these residues with alanine in the mutant protein prevents its cleavage, despite preserving its normal secondary structure. The N-terminal BB0323 protein binds BbHtrA, but its cleavage site mutant displays deficiency in such interaction. An isogenic borrelial mutant with NL-to-AA substitution in BB0323 (referred to as Bb^bb0323NL) maintains normal growth yet is impaired for infection of mice or transmission from infected ticks. Notably, the BB0323 protein is still processed in Bb^bb0323NL, albeit with lower levels of mature N-terminal BB0323 protein and multiple aberrantly processed polypeptides, which could result from nonspecific cleavages at other asparagine and leucine residues in the protein. The lack of infectivity of Bb^bb0323NL is likely due to the impaired abundance or stoichiometry of a protein complex involving BB0238, another spirochete protein. Together, these studies highlight that a precise proteolytic event and a particular protein-protein interaction, involving multiple borrelial virulence determinants, are mutually inclusive and interconnected, playing essential roles in the infectivity of Lyme disease pathogens.

Deletion of a Genetic Region of lp17 Affects Plasmid Copy Number in Borrelia burgdorferi

Borrelia burgdorferi, the Lyme disease pathogen, is maintained in its enzootic life cycle through complex gene regulatory pathways encoded on its uniquely fragmented genome. This genome consists of over 20 plasmids, and the regulatory mechanisms of plasmid maintenance and replication are largely unknown. The bbd21 gene, encoded on lp17 and a member of the paralogous family 32 proteins, was originally proposed to be a putative parA orthologue involved with plasmid partitioning; however, this function has not been confirmed to date. To determine the role of bbd21 in B. burgdorferi, we utilized targeted gene deletion and discovered bbd21 and bbd22 are co-transcribed. The effects of bbd21 and bbd22 deletion on plasmid copy number and mammalian infectivity were assessed. By qPCR, lp17 copy number did not differ amongst strains during mid-exponential and stationary growth phases. However, after in vitro passaging, the mutant strain demonstrated an 8-fold increase in lp17 copies, suggesting a cumulative defect in plasmid copy number regulation. Additionally, we compared lp17 copy number between in vitro and mammalian host-adapted conditions. Our findings showed 1) lp17 copy number was significantly different between these growth conditions for both the wild type and bbd21-bbd22 deletion mutant and 2) under mammalian host-adapted cultivation, the absence of bbd21-bbd22 resulted in significantly decreased copies of lp17. Murine infection studies using culture and qPCR demonstrated bbd21-bbd22 deletion resulted in a tissue colonization defect, particularly in the heart. Lastly, we showed bbd21 transcription appears to be independent of direct rpoS regulation based on similar expression levels in wild type and ΔrpoS. Altogether, our findings indicate the bbd21-bbd22 genetic region is involved with regulation of lp17 plasmid copy number. Furthermore, we propose the possibility that lp17 plasmid copy number is important for microbial pathogenesis by the Lyme disease spirochete.