ABSTRACT
Lyme disease, caused by infection with members of the Lyme borreliosis group of Borrelia spirochete bacteria, is increasing in frequency and distribution worldwide. Epigenetic interactions between the mammalian host, tick, and bacterial pathogen are poorly understood. In this study, high-throughput next-generation sequencing (NGS) allowed for the in vitro study of the transcriptome, non-coding RNAs, and methylome in human host cells in response to Borrelia burgdorferi infection. We tested the effect of the Borrelia burgdorferi strain B31 on a human primary cell line (HUVEC) and an immortalized cell line (HEK-293) for 72 h, a long-duration time that might allow for epigenetic responses in the exposed human host cells. Differential gene expression was detected in both cell models in response to B. burgdorferi. More differentially expressed genes were found in HUVECs compared to HEK-293 cells. Borrelia burgdorferi exposure significantly induced genes in the interferon, in addition to cytokine and other immune response signaling in HUVECs. In HEK-293 cells, pre-NOTCH processing in Golgi was significantly downregulated in Borrelia-exposed cells. Other significantly altered gene expressions were found in genes involved in the extracellular matrix. No significant global methylation changes were detected in HUVECs or HEK-293 cells exposed to B. burgdorferi; however, two long non-coding RNAs and a pseudogene were deregulated in response to B. burgdorferi in HUVECs, suggesting that other epigenetic mechanisms may be initiated by infection.
Subject(s)
Borrelia burgdorferi , Lyme Disease , RNA, Long Noncoding , Transcriptome , Humans , Borrelia burgdorferi/genetics , RNA, Long Noncoding/genetics , Lyme Disease/microbiology , Lyme Disease/genetics , HEK293 Cells , Human Umbilical Vein Endothelial Cells/metabolism , Human Umbilical Vein Endothelial Cells/microbiology , Epigenome , DNA Methylation , High-Throughput Nucleotide Sequencing , Host-Pathogen Interactions/genetics , Epigenesis, GeneticABSTRACT
BACKGROUND: Postinfectious Lyme arthritis (LA) is associated with dysregulated immunity and autoreactive T- and B-cell responses in joints. Here we explored the role of host genetic variation in this outcome. METHODS: The frequency of 253 702 single-nucleotide polymorphisms (SNPs) was determined in 147 patients with LA (87 with postinfectious LA and 60 with antibiotic-responsive LA), and for comparison in 90 patients with erythema migrans or the general population (n = 2504). Functional outcome of candidate SNPs was assessed by evaluating their impact on clinical outcome and on immune responses in blood and synovial fluid in patients with LA. RESULTS: Six SNPs associated with late cornified envelope (LCE3) genes were present at greater frequency in patients with postinfectious LA compared to those with antibiotic-responsive LA (70% vs 30%; odds ratio, 2; P < .01). These SNPs were associated with heightened levels of inflammatory Th17 cytokines in serum but lower levels of interleukin 27, a regulatory cytokine, implying that they may contribute to dysregulated Th17 immunity in blood. Moreover, in patients with postinfectious LA, the levels of these Th17 mediators correlated directly with autoantibody responses in synovial fluid, providing a possible link between LCE3 SNPs, maladaptive systemic Th17 immunity, and autoreactive responses in joints. CONCLUSIONS: Variation in the LCE3 locus, a known genetic risk factor in psoriasis and psoriatic arthritis, is associated with dysregulated systemic Th17 immunity and heightened autoantibody responses in joints. These findings underscore the importance of host genetic predisposition and systemic Th17 immunity in the pathogenesis of postinfectious (antibiotic-refractory) Lyme arthritis.
Subject(s)
Lyme Disease , Polymorphism, Single Nucleotide , Th17 Cells , Humans , Lyme Disease/genetics , Lyme Disease/immunology , Th17 Cells/immunology , Male , Female , Adult , Middle Aged , Synovial Fluid/immunology , Aged , Cytokines/genetics , Cytokines/metabolism , Arthritis, Infectious/genetics , Arthritis, Infectious/immunology , Young AdultABSTRACT
BACKGROUNDAntibiotic-Refractory Lyme Arthritis (ARLA) involves a complex interplay of T cell responses targeting Borrelia burgdorferi antigens progressing toward autoantigens by epitope spreading. However, the precise molecular mechanisms driving the pathogenic T cell response in ARLA remain unclear. Our aim was to elucidate the molecular program of disease-specific Th cells.METHODSUsing flow cytometry, high-throughput T cell receptor (TCR) sequencing, and scRNA-Seq of CD4+ Th cells isolated from the joints of patients with ARLA living in Europe, we aimed to infer antigen specificity through unbiased analysis of TCR repertoire patterns, identifying surrogate markers for disease-specific TCRs, and connecting TCR specificity to transcriptional patterns.RESULTSPD-1hiHLA-DR+CD4+ effector T cells were clonally expanded within the inflamed joints and persisted throughout disease course. Among these cells, we identified a distinct TCR-ß motif restricted to HLA-DRB1*11 or *13 alleles. These alleles, being underrepresented in patients with ARLA living in North America, were unexpectedly prevalent in our European cohort. The identified TCR-ß motif served as surrogate marker for a convergent TCR response specific to ARLA, distinguishing it from other rheumatic diseases. In the scRNA-Seq data set, the TCR-ß motif particularly mapped to peripheral T helper (TPH) cells displaying signs of sustained proliferation, continuous TCR signaling, and expressing CXCL13 and IFN-γ.CONCLUSIONBy inferring disease-specific TCRs from synovial T cells we identified a convergent TCR response in the joints of patients with ARLA that continuously fueled the expansion of TPH cells expressing a pathogenic cytokine effector program. The identified TCRs will aid in uncovering the major antigen targets of the maladaptive immune response.FUNDINGSupported by the German Research Foundation (DFG) MO 2160/4-1; the Federal Ministry of Education and Research (BMBF; Advanced Clinician Scientist-Program INTERACT; 01EO2108) embedded in the Interdisciplinary Center for Clinical Research (IZKF) of the University Hospital Würzburg; the German Center for Infection Research (DZIF; Clinical Leave Program; TI07.001_007) and the Interdisciplinary Center for Clinical Research (IZKF) Würzburg (Clinician Scientist Program, Z-2/CSP-30).
Subject(s)
HLA-DRB1 Chains , Lyme Disease , T-Lymphocytes, Helper-Inducer , Humans , Lyme Disease/immunology , Lyme Disease/pathology , Lyme Disease/genetics , HLA-DRB1 Chains/genetics , HLA-DRB1 Chains/immunology , Female , Male , T-Lymphocytes, Helper-Inducer/immunology , Borrelia burgdorferi/immunology , Middle Aged , Receptors, Antigen, T-Cell/immunology , Receptors, Antigen, T-Cell/genetics , Adult , Receptors, Antigen, T-Cell, alpha-beta/genetics , Receptors, Antigen, T-Cell, alpha-beta/immunologyABSTRACT
The incidence of Lyme borreliosis has risen, accompanied by persistent symptoms. The innate immune system and related cytokines are crucial in the host response and symptom development. We characterized cytokine production capacity before and after antibiotic treatment in 1,060 Lyme borreliosis patients. We observed a negative correlation between antibody production and IL-10 responses, as well as increased IL-1Ra responses in patients with disseminated disease. Genome-wide mapping the cytokine production allowed us to identify 34 cytokine quantitative trait loci (cQTLs), with 31 novel ones. We pinpointed the causal variant at the TLR1-6-10 locus and validated the regulation of IL-1Ra responses at transcritpome level using an independent cohort. We found that cQTLs contribute to Lyme borreliosis susceptibility and are relevant to other immune-mediated diseases. Our findings improve the understanding of cytokine responses in Lyme borreliosis and provide a genetic map of immune function as an expanded resource.
Subject(s)
Cytokines , Lyme Disease , Quantitative Trait Loci , Lyme Disease/immunology , Lyme Disease/genetics , Lyme Disease/microbiology , Humans , Cytokines/genetics , Cytokines/metabolism , Male , Female , Interleukin-10/genetics , Adult , Genome-Wide Association Study , Middle Aged , Interleukin 1 Receptor Antagonist Protein/genetics , Borrelia burgdorferi/immunology , Borrelia burgdorferi/genetics , Anti-Bacterial Agents , Polymorphism, Single Nucleotide , Genetic Predisposition to Disease , AgedABSTRACT
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.
Subject(s)
Bacteriophages , Borrelia burgdorferi , Lyme Disease , Humans , Borrelia burgdorferi/genetics , Bacteriophages/genetics , Plasmids/genetics , Lyme Disease/genetics , Genomics , DNAABSTRACT
The complement system serves as the first line of defense against invading pathogens by promoting opsonophagocytosis and bacteriolysis. Antibody-dependent activation of complement occurs through the classical pathway and relies on the activity of initiating complement proteases of the C1 complex, C1r and C1s. The causative agent of Lyme disease, Borrelia burgdorferi, expresses two paralogous outer surface lipoproteins of the OspEF-related protein family, ElpB and ElpQ, that act as specific inhibitors of classical pathway activation. We have previously shown that ElpB and ElpQ bind directly to C1r and C1s with high affinity and specifically inhibit C2 and C4 cleavage by C1s. To further understand how these novel protease inhibitors function, we carried out a series of hydrogen-deuterium exchange mass spectrometry (HDX-MS) experiments using ElpQ and full-length activated C1s as a model of Elp-protease interaction. Comparison of HDX-MS profiles between unbound ElpQ and the ElpQ/C1s complex revealed a putative C1s-binding site on ElpQ. HDX-MS-guided, site-directed ElpQ mutants were generated and tested for direct binding to C1r and C1s using surface plasmon resonance. Several residues within the C-terminal region of ElpQ were identified as important for protease binding, including a single conserved tyrosine residue that was required for ElpQ- and ElpB-mediated complement inhibition. Collectively, our study identifies key molecular determinants for classical pathway protease recognition by Elp proteins. This investigation improves our understanding of the unique complement inhibitory mechanism employed by Elp proteins which serve as part of a sophisticated complement evasion system present in Lyme disease spirochetes.
Subject(s)
Bacterial Outer Membrane Proteins , Borrelia burgdorferi , Complement Pathway, Classical , Humans , Bacterial Outer Membrane Proteins/chemistry , Bacterial Outer Membrane Proteins/genetics , Bacterial Outer Membrane Proteins/metabolism , Borrelia burgdorferi/immunology , Borrelia burgdorferi/metabolism , Borrelia burgdorferi/genetics , Complement C1r/metabolism , Complement C1r/genetics , Complement C1s/metabolism , Complement C1s/genetics , Complement C1s/chemistry , Complement Pathway, Classical/immunology , Lipoproteins/metabolism , Lipoproteins/genetics , Lipoproteins/chemistry , Lipoproteins/immunology , Lyme Disease/genetics , Lyme Disease/immunology , Lyme Disease/microbiology , Protein BindingABSTRACT
BACKGROUND: Genetic variation underly inter-individual variation in host immune responses to infectious diseases, and may affect susceptibility or the course of signs and symptoms. METHODS: We performed genome-wide association studies in a prospective cohort of 1138 patients with physician-confirmed Lyme borreliosis (LB), the most common tick-borne disease in the Northern hemisphere caused by the bacterium Borrelia burgdorferi sensu lato. Genome-wide variants in LB patients-divided into a discovery and validation cohort-were compared to two healthy cohorts. Additionally, ex vivo monocyte-derived cytokine responses of peripheral blood mononuclear cells to several stimuli including Borrelia burgdorferi were performed in both LB patient and healthy control samples, as were stimulation experiments using mechanistic/mammalian target of rapamycin (mTOR) inhibitors. In addition, for LB patients, anti-Borrelia antibody responses were measured. Finally, in a subset of LB patients, gene expression was analysed using RNA-sequencing data from the ex vivo stimulation experiments. RESULTS: We identified a previously unknown genetic variant, rs1061632, that was associated with enhanced LB susceptibility. This polymorphism was an eQTL for KCTD20 and ETV7 genes, and its major risk allele was associated with upregulation of the mTOR pathway and cytokine responses, and lower anti-Borrelia antibody production. In addition, we replicated the recently reported SCGB1D2 locus that was suggested to have a protective effect on B. burgdorferi infection, and associated this locus with higher Borrelia burgdorferi antibody indexes and lower IL-10 responses. CONCLUSIONS: Susceptibility for LB was associated with higher anti-inflammatory responses and reduced anti-Borrelia antibody production, which in turn may negatively impact bacterial clearance. These findings provide important insights into the immunogenetic susceptibility for LB and may guide future studies on development of preventive or therapeutic measures. TRIAL REGISTRATION: The LymeProspect study was registered with the International Clinical Trials Registry Platform (NTR4998, registration date 2015-02-13).
Subject(s)
Borrelia burgdorferi Group , Borrelia burgdorferi , Borrelia , Lyme Disease , Humans , Genome-Wide Association Study , Prospective Studies , Leukocytes, Mononuclear , Disease Susceptibility , Lyme Disease/genetics , Lyme Disease/diagnosis , Borrelia burgdorferi/genetics , Cytokines/genetics , TOR Serine-Threonine Kinases/genetics , TOR Serine-Threonine Kinases/therapeutic use , Borrelia burgdorferi Group/genetics , Secretoglobins/geneticsABSTRACT
The σ54-σS sigma factor cascade plays a central role in regulating differential gene expression during the enzootic cycle of Borreliella burgdorferi, the Lyme disease pathogen. In this pathway, the primary transcription of rpoS (which encodes σS) is under the control of σ54 which is activated by a bacterial enhancer-binding protein (EBP), Rrp2. The σ54-dependent activation in B. burgdorferi has long been thought to be unique, requiring an additional factor, BosR, a homologue of classical Fur/PerR repressor/activator. However, how BosR is involved in this σ54-dependent activation remains unclear and perplexing. In this study, we demonstrate that BosR does not function as a regulator for rpoS transcriptional activation. Instead, it functions as a novel RNA-binding protein that governs the turnover rate of rpoS mRNA. We further show that BosR directly binds to the 5' untranslated region (UTR) of rpoS mRNA, and the binding region overlaps with a region required for rpoS mRNA degradation. Mutations within this 5'UTR region result in BosR-independent RpoS production. Collectively, these results uncover a novel role of Fur/PerR family regulators as RNA-binding proteins and redefine the paradigm of the σ54-σS pathway in B. burgdorferi.
Subject(s)
Bacterial Proteins , Borrelia burgdorferi , Gene Expression Regulation, Bacterial , RNA Stability , RNA-Binding Proteins , Sigma Factor , Sigma Factor/metabolism , Sigma Factor/genetics , Borrelia burgdorferi/genetics , Borrelia burgdorferi/metabolism , Bacterial Proteins/metabolism , Bacterial Proteins/genetics , RNA Stability/genetics , RNA-Binding Proteins/metabolism , RNA-Binding Proteins/genetics , 5' Untranslated Regions , Lyme Disease/microbiology , Lyme Disease/genetics , Repressor Proteins/metabolism , Repressor Proteins/genetics , RNA, Messenger/metabolism , RNA, Messenger/genetics , RNA Polymerase Sigma 54/metabolism , RNA Polymerase Sigma 54/geneticsABSTRACT
Introduction: Lyme disease, the most common tick-borne infectious disease in the US, is caused by a spirochetal pathogen Borrelia burgdorferi (Bb). Distinct host responses are observed in susceptible and resistant strains of inbred of mice following infection with Bb reflecting a subset of inflammatory responses observed in human Lyme disease. The advent of post-genomic methodologies and genomic data sets enables dissecting the host responses to advance therapeutic options for limiting the pathogen transmission and/or treatment of Lyme disease. Methods: In this study, we used single-cell RNA-Seq analysis in conjunction with mouse genomics exploiting GFP-expressing Bb to sort GFP+ splenocytes and GFP- bystander cells to uncover novel molecular and cellular signatures that contribute to early stages of immune responses against Bb. Results: These data decoded the heterogeneity of splenic neutrophils, macrophages, NK cells, B cells, and T cells in C3H/HeN mice in response to Bb infection. Increased mRNA abundance of apoptosis-related genes was observed in neutrophils and macrophages clustered from GFP+ splenocytes. Moreover, complement-mediated phagocytosis-related genes such as C1q and Ficolin were elevated in an inflammatory macrophage subset, suggesting upregulation of these genes during the interaction of macrophages with Bb-infected neutrophils. In addition, the role of DUSP1 in regulating the expression of Casp3 and pro-inflammatory cytokines Cxcl1, Cxcl2, Il1b, and Ccl5 in Bb-infected neutrophils were identified. Discussion: These findings serve as a growing catalog of cell phenotypes/biomarkers among murine splenocytes that can be exploited for limiting spirochetal burden to limit the transmission of the agent of Lyme disease to humans via reservoir hosts.
Subject(s)
Borrelia burgdorferi , Lyme Disease , Mice , Humans , Animals , Borrelia burgdorferi/genetics , Transcriptome , Spleen , Single-Cell Gene Expression Analysis , Mice, Inbred C3H , Lyme Disease/geneticsABSTRACT
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.
Subject(s)
Borrelia burgdorferi , Lyme Disease , Ticks , Animals , Mice , Histidine Kinase/genetics , Histidine Kinase/metabolism , Bacterial Proteins/genetics , Bacterial Proteins/metabolism , Virulence , Chemotaxis , Lyme Disease/genetics , Lyme Disease/microbiology , Ticks/microbiology , Virulence Factors/genetics , Virulence Factors/metabolism , Gene Expression Regulation, Bacterial , Sigma Factor/genetics , Sigma Factor/metabolismABSTRACT
Ixodes scapularis, the black-legged tick, is the principal vector of the Lyme disease spirochete, Borrelia burgdorferi, and is responsible for most of the â¼470,000 estimated Lyme disease cases annually in the USA. Ixodes scapularis can transmit six additional pathogens of human health significance. Because of its medical importance, I. scapularis was the first tick genome to be sequenced and annotated. However, the first assembly, I. scapularis Wikel (IscaW), was highly fragmented because of the technical challenges posed by the long, repetitive genome sequences characteristic of arthropod genomes and the lack of long-read sequencing techniques. Although I. scapularis has emerged as a model for tick research because of the availability of new tools such as embryo injection and CRISPR-Cas9-mediated gene editing yet the lack of chromosome-scale scaffolds has slowed progress in tick biology and the development of tools for their control. Here we combine diverse technologies to produce the I. scapularis Gulia-Nuss (IscGN) genome assembly and gene set. We used DNA from eggs and male and female adult ticks and took advantage of Hi-C, PacBio HiFi sequencing, and Illumina short-read sequencing technologies to produce a chromosome-level assembly. In this work, we present the predicted pseudochromosomes consisting of 13 autosomes and the sex pseudochromosomes: X and Y, and a markedly improved genome annotation compared with the existing assemblies and annotations.
Subject(s)
Borrelia burgdorferi , Ixodes , Lyme Disease , Animals , Male , Female , Humans , Ixodes/genetics , Lyme Disease/genetics , Borrelia burgdorferi/genetics , Genome/genetics , High-Throughput Nucleotide SequencingABSTRACT
Transcription termination is an essential and dynamic process that can tune gene expression in response to diverse molecular signals. Yet, the genomic positions, molecular mechanisms, and regulatory consequences of termination have only been studied thoroughly in model bacteria. Here, we use several RNA-seq approaches to map RNA ends for the transcriptome of the spirochete Borrelia burgdorferi - the etiological agent of Lyme disease. We identify complex gene arrangements and operons, untranslated regions and small RNAs. We predict intrinsic terminators and experimentally test examples of Rho-dependent transcription termination. Remarkably, 63% of RNA 3' ends map upstream of or internal to open reading frames (ORFs), including genes involved in the unique infectious cycle of B. burgdorferi. We suggest these RNAs result from premature termination, processing and regulatory events such as cis-acting regulation. Furthermore, the polyamine spermidine globally influences the generation of truncated mRNAs. Collectively, our findings provide insights into transcription termination and uncover an abundance of potential RNA regulators in B. burgdorferi.
Subject(s)
Borrelia burgdorferi , Lyme Disease , Humans , Borrelia burgdorferi/genetics , Transcriptome , RNA , Lyme Disease/genetics , Lyme Disease/microbiology , Transcription, Genetic , Gene Expression Regulation, Bacterial , Bacterial Proteins/genetics , Bacterial Proteins/metabolismABSTRACT
PlzA is a c-di-GMP-binding protein crucial for adaptation of the Lyme disease spirochete Borrelia (Borreliella) burgdorferi during its enzootic life cycle. Unliganded apo-PlzA is important for vertebrate infection, while liganded holo-PlzA is important for survival in the tick; however, the biological function of PlzA has remained enigmatic. Here, we report that PlzA has RNA chaperone activity that is inhibited by c-di-GMP binding. Holo- and apo-PlzA bind RNA and accelerate RNA annealing, while only apo-PlzA can strand displace and unwind double-stranded RNA. Guided by the crystal structure of PlzA, we identified several key aromatic amino acids protruding from the N- and C-terminal domains that are required for RNA-binding and unwinding activity. Our findings illuminate c-di-GMP as a switch controlling the RNA chaperone activity of PlzA, and we propose that complex RNA-mediated modulatory mechanisms allow PlzA to regulate gene expression during both the vector and host phases of the B. burgdorferi life cycle.
Subject(s)
Borrelia burgdorferi Group , Borrelia burgdorferi , Ixodes , Lyme Disease , Bacterial Proteins/metabolism , Borrelia burgdorferi/metabolism , Borrelia burgdorferi Group/genetics , Lyme Disease/genetics , RNA/metabolismABSTRACT
The blacklegged tick (Ixodes scapularis (Journal of the Academy of Natural Sciences of Philadelphia, 1821, 2, 59)) is a vector of Borrelia burgdorferi sensu stricto (s.s.) (International Journal of Systematic Bacteriology, 1984, 34, 496), the causative bacterial agent of Lyme disease, part of a slow-moving epidemic of Lyme borreliosis spreading across the northern hemisphere. Well-known geographical differences in the vectorial capacity of these ticks are associated with genetic variation. Despite the need for detailed genetic information in this disease system, previous phylogeographical studies of these ticks have been restricted to relatively few populations or few genetic loci. Here we present the most comprehensive phylogeographical study of genome-wide markers in I. scapularis, conducted by using 3RAD (triple-enzyme restriction-site associated sequencing) and surveying 353 ticks from 33 counties throughout the species' range. We found limited genetic variation among populations from the Northeast and Upper Midwest, where Lyme disease is most common, and higher genetic variation among populations from the South. We identify five spatially associated genetic clusters of I. scapularis. In regions where Lyme disease is increasing in frequency, the I. scapularis populations genetically group with ticks from historically highly Lyme-endemic regions. Finally, we identify 10 variable DNA sites that contribute the most to population differentiation. These variable sites cluster on one of the chromosome-scale scaffolds for I. scapularis and are within identified genes. Our findings illuminate the need for additional research to identify loci causing variation in the vectorial capacity of I. scapularis and where additional tick sampling would be most valuable to further understand disease trends caused by pathogens transmitted by I. scapularis.
Subject(s)
Borrelia burgdorferi , Ixodes , Lyme Disease , Animals , Ixodes/genetics , Ixodes/microbiology , Phylogeography , Lyme Disease/genetics , Lyme Disease/microbiology , Borrelia burgdorferi/genetics , BacteriaABSTRACT
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.
Subject(s)
Borrelia burgdorferi , Lyme Disease , Animals , Humans , Bacterial Proteins/metabolism , Borrelia burgdorferi/genetics , Borrelia burgdorferi/metabolism , Lyme Disease/genetics , Mammals/metabolism , TranscriptomeABSTRACT
The RNA polymerase alternative σ factor RpoS in Borrelia burgdorferi (Bb), the Lyme disease pathogen, is responsible for programmatic-positive and -negative gene regulation essential for the spirochete's dual-host enzootic cycle. RpoS is expressed during tick-to-mammal transmission and throughout mammalian infection. Although the mammalian-phase RpoS regulon is well described, its counterpart during the transmission blood meal is unknown. Here, we used Bb-specific transcript enrichment by tick-borne disease capture sequencing (TBDCapSeq) to compare the transcriptomes of WT and ΔrpoS Bb in engorged nymphs and following mammalian host-adaptation within dialysis membrane chambers. TBDCapSeq revealed dramatic changes in the contours of the RpoS regulon within ticks and mammals and further confirmed that RpoS-mediated repression is specific to the mammalian-phase of Bb's enzootic cycle. We also provide evidence that RpoS-dependent gene regulation, including repression of tick-phase genes, is required for persistence in mice. Comparative transcriptomics of engineered Bb strains revealed that the Borrelia oxidative stress response regulator (BosR), a noncanonical Fur family member, and the cyclic diguanosine monophosphate (c-di-GMP) effector PlzA reciprocally regulate the function of RNA polymerase complexed with RpoS. BosR is required for RpoS-mediated transcription activation and repression in addition to its well-defined role promoting transcription of rpoS by the RNA polymerase alternative σ factor RpoN. During transmission, ligand-bound PlzA antagonizes RpoS-mediated repression, presumably acting through BosR.
Subject(s)
Borrelia burgdorferi , Borrelia , Lyme Disease , Ticks , Mice , Animals , Borrelia burgdorferi/genetics , Borrelia/metabolism , Bacterial Proteins/genetics , Bacterial Proteins/metabolism , Ticks/genetics , Sigma Factor/genetics , Sigma Factor/metabolism , Lyme Disease/genetics , Mammals/metabolism , Gene Expression Regulation, BacterialABSTRACT
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.
Subject(s)
Borrelia burgdorferi , Lyme Disease , Humans , Borrelia burgdorferi/genetics , Polyploidy , DNA , Lyme Disease/genetics , Chromosome SegregationABSTRACT
Intracellular compartmentalization of ligands, receptors and signaling molecules has been recognized as an important regulator of inflammation. The toll-like receptor (TLR) 2 pathway utilizes the trafficking molecule adaptor protein 3 (AP-3) to activate interleukin (IL)-6 signaling from within phagosomal compartments. To better understand the vesicular pathways that may contribute to intracellular signaling and cooperate with AP-3, we performed a vesicular siRNA screen. We identified Rab8 and Rab11 GTPases as important in IL-6 induction upon stimulation with the TLR2 ligand Pam3 CSK4 or the pathogen, Borrelia burgdorferi (Bb), the causative agent of Lyme disease. These Rabs were recruited to late and lysosomal stage phagosomes and co-transported with TLR2 signaling adaptors and effectors, such as MyD88, TRAM and TAK1, in an AP-3-dependent manner. Our data support a model where AP-3 mediates the recruitment of recycling and secretory vesicles and the assembly of signaling complexes at the phagosome.
Subject(s)
Borrelia burgdorferi , Lyme Disease , Adaptor Proteins, Signal Transducing/metabolism , Borrelia burgdorferi/metabolism , Ligands , Lyme Disease/genetics , Lyme Disease/metabolism , Myeloid Differentiation Factor 88/genetics , Myeloid Differentiation Factor 88/metabolism , Phagosomes/metabolism , Toll-Like Receptor 2/genetics , Toll-Like Receptor 2/metabolism , rab GTP-Binding Proteins , Animals , MiceABSTRACT
The acceleration of climate change has been associated with an alarming increase in the prevalence and geographic range of tick-borne diseases (TBD), many of which have severe and long-lasting effects-particularly when treatment is delayed principally due to inadequate diagnostics and lack of physician suspicion. Moreover, there is a paucity of treatment options for many TBDs that are complicated by diagnostic limitations for correctly identifying the offending pathogens. This review will focus on the biology, disease pathology, and detection methodologies used for the Borreliaceae family which includes the Lyme disease agent Borreliella burgdorferi. Previous work revealed that Borreliaceae genomes differ from most bacteria in that they are composed of large numbers of replicons, both linear and circular, with the main chromosome being the linear with telomeric-like termini. While these findings are novel, additional gene-specific analyses of each class of these multiple replicons are needed to better understand their respective roles in metabolism and pathogenesis of these enigmatic spirochetes. Historically, such studies were challenging due to a dearth of both analytic tools and a sufficient number of high-fidelity genomes among the various taxa within this family as a whole to provide for discriminative and functional genomic studies. Recent advances in long-read whole-genome sequencing, comparative genomics, and machine-learning have provided the tools to better understand the fundamental biology and phylogeny of these genomically-complex pathogens while also providing the data for the development of improved diagnostics and therapeutics.
Subject(s)
Borrelia burgdorferi , Lyme Disease , Borrelia burgdorferi/genetics , Genome, Bacterial , Genomics/methods , Humans , Lyme Disease/genetics , Lyme Disease/microbiology , PhylogenyABSTRACT
Host association-the selective adaptation of pathogens to specific host species-evolves through constant interactions between host and pathogens, leaving a lot yet to be discovered on immunological mechanisms and genomic determinants. The causative agents of Lyme disease (LD) are spirochete bacteria composed of multiple species of the Borrelia burgdorferi sensu lato complex, including B. burgdorferi (Bb), the main LD pathogen in North America-a useful model for the study of mechanisms underlying host-pathogen association. Host adaptation requires pathogens' ability to evade host immune responses, such as complement, the first-line innate immune defense mechanism. We tested the hypothesis that different host-adapted phenotypes among Bb strains are linked to polymorphic loci that confer complement evasion traits in a host-specific manner. We first examined the survivability of 20 Bb strains in sera in vitro and/or bloodstream and tissues in vivo from rodent and avian LD models. Three groups of complement-dependent host-association phenotypes emerged. We analyzed complement-evasion genes, identified a priori among all strains and sequenced and compared genomes for individual strains representing each phenotype. The evolutionary history of ospC loci is correlated with host-specific complement-evasion phenotypes, while comparative genomics suggests that several gene families and loci are potentially involved in host association. This multidisciplinary work provides novel insights into the functional evolution of host-adapted phenotypes, building a foundation for further investigation of the immunological and genomic determinants of host association. IMPORTANCE Host association is the phenotype that is commonly found in many pathogens that preferential survive in particular hosts. The Lyme disease (LD)-causing agent, B. burgdorferi (Bb), is an ideal model to study host association, as Bb is mainly maintained in nature through rodent and avian hosts. A widespread yet untested concept posits that host association in Bb strains is linked to Bb functional genetic variation conferring evasion to complement, an innate defense mechanism in vertebrate sera. Here, we tested this concept by grouping 20 Bb strains into three complement-dependent host-association phenotypes based on their survivability in sera and/or bloodstream and distal tissues in rodent and avian LD models. Phylogenomic analysis of these strains further correlated several gene families and loci, including ospC, with host-specific complement-evasion phenotypes. Such multifaceted studies thus pave the road to further identify the determinants of host association, providing mechanistic insights into host-pathogen interaction.