Lost in plasmids: next generation sequencing and the complex genome of the tick-borne pathogen Borrelia burgdorferi

Multi-Omics Technologies Applied to Improve Tick Research

The advancement of multi-omics technologies is crucial to deepen knowledge on tick biology. These approaches, used to study diverse phenomena, are applied to experiments that aim to understand changes in gene transcription, protein function, cellular processes, and prediction of systems at global biological levels. This review addressed the application of omics data to investigate and elucidate tick physiological processes, such as feeding, digestion, reproduction, neuronal, endocrine systems, understanding population dynamics, transmitted pathogens, control, and identifying new vaccine targets. Furthermore, new therapeutic perspectives using tick bioactive molecules, such as anti-inflammatory, analgesic, and antitumor, were summarized. Taken together, the application of omics technologies can help to understand the protein functions and biological behavior of ticks, as well as the identification of potential new antigens influencing the development of alternative control strategies and, consequently, the tick-borne disease prevention in veterinary and public health contexts. Finally, tick population dynamics have been determined through a combination of environmental factors, host availability, and genetic adaptations, and recent advances in omics technologies have improved our understanding of their ecological resilience and resistance mechanisms. Future directions point to the integration of spatial omics and artificial intelligence to further unravel tick biology and improve control strategies.

Complex exchanges among plasmids and clonal expansion of lineages shape the population structure and virulence of Borrelia burgdorferi

Background

In the United States, Borrelia burgdorferi (Bb) is the principal etiologic agent of Lyme disease. The complex structure of Bb genomes has posed challenges for genomic studies because homology among the bacterium's many plasmids, which account for ~40% of the genome by length, has made them difficult to sequence and assemble.

Results

We used long-read sequencing to generate near-complete assemblies of 62 isolates of human-derived Bb and collected public genomes with plasmid sequences. We characterized genetic diversity and population structure in the resulting set of 82 plasmid-complete Borrelia burgdorferi sensu stricto genomes. The Bb core genome is encoded by a chromosome and the conserved plasmids cp26, lp54, and lp17; the accessory genome is encoded by all other plasmids and the distal arm of the chromosome. Near-complete genomes reveal that the most granular Bb genotypes are clonal expansions of complex rearrangements among accessory genome elements. Ribosomal spacer types (RST) represent multiple collections of such genotypes, whereas OspC types are usually clonal. Structural rearrangements are non-randomly distributed throughout the genome, with cp32 plasmids undergoing dense exchanges and most linear plasmids, except lp54, sharing blocks among themselves and with the distal arm of the chromosome. OspC type A strains, known to possess greater virulence in humans, are distinguished by the presence of lp28-1 and lp56. Rearrangements among plasmids tended to preserve gene content, suggesting functional constraints among gene networks. Using k-partite graph decompositions, we identified gene sets with correlation patterns suggestive of conserved functional modules.

Conclusions

Long-read assemblies reveal that Bb population genetic structure results from clonal expansion of lineages that have undergone complex rearrangements among plasmid-encoded accessory genome elements. Genetic structure is preserved among genes even when plasmid rearrangements occur, suggesting that selection among epistatic loci maintains functional genetic networks. The analysis of near-complete genomes assembled using long-read sequencing methods advances our understanding of Bb biology and Lyme disease pathogenesis by providing the first detailed view of population variation in previously inaccessible areas of the Bb genome.

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

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

Description of genome sequences of arthropod-associated spirochetes of the genus Entomospira

Spirochetal bacteria isolated from arthropods of the genera Culex and Aedes are termed BR149, BR151 (Entomospira culicis), BR193 (Entomospira entomophila), and BR208 (Entomospira nematocerorum). Genome sizes assembled from Illumina MiSeq and Oxford Nanopore reads varied between 1.67 and 1.78 Mb containing three to six plasmids. GC content ranged from 38.5% to 45.76%.

Optimising Transformation Efficiency in Borrelia: Unravelling the Role of the Restriction-Modification System of Borrelia afzelii and Borrelia garinii

Borrelia spp. are transmitted to humans by the bite of an infected tick. In Europe, Borrelia afzelii and Borrelia garinii are the main causative agents of Lyme borreliosis, one of the most prevalent tick-borne diseases in the northern hemisphere. In bacteria such as Borrelia spp., a restriction-modification system (RMS) protects against the harmful introduction of foreign DNA. The RMS comprises two activities: methyltransferase and endonuclease. This study is aimed to characterize the RMS of B. afzelii and B. garinii. First, we identified potential RMS genes. The predicted genes were cloned into a methylase-deficient Escherichia coli strain and digested with methylation-sensitive restriction enzymes to verify methyltransferase activity. Additionally, the RMS proteins were purified to evaluate endonuclease activity. Subsequently, methylated and unmethylated plasmids were used to investigate the effect of methylation on endonuclease activity and transformation efficiency. We identified four possible RMS genes in B. afzelii and four RMS genes in B. garinii. We analyzed the presence of these genes in patient isolates and observed a high degree of heterogeneity. The restriction pattern of DNA methylated by each of the four recombinantly expressed genes provided strong evidence that all encode adenine-specific methyltransferases. After 24 h of incubation with purified RMS proteins, we observed complete digestion of unmethylated plasmid DNA, demonstrating endonuclease activity. Finally, we proved that methylation protects against endonuclease activity and increases transformation efficiency.

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

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

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

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

IMPORTANCE

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

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

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

Unprecedented genetic variability of PFam54 paralogs among Eurasian Lyme borreliosis-causing spirochetes

Lyme borreliosis (LB) is the most common vector-borne disease in the Northern Hemisphere caused by spirochetes belonging to the Borrelia burgdorferi sensu lato (Bbsl) complex. Borrelia spirochetes circulate in obligatory transmission cycles between tick vectors and different vertebrate hosts. To successfully complete this complex transmission cycle, Bbsl encodes for an arsenal of proteins including the PFam54 protein family with known, or proposed, influences to reservoir host and/or vector adaptation. Even so, only fragmentary information is available regarding the naturally occurring level of variation in the PFam54 gene array especially in relation to Eurasian-distributed species. Utilizing whole genome data from isolates (n = 141) originated from three major LB-causing Borrelia species across Eurasia (B. afzelii, B. bavariensis, and B. garinii), we aimed to characterize the diversity of the PFam54 gene array in these isolates to facilitate understanding the evolution of PFam54 paralogs on an intra- and interspecies level. We found an extraordinarily high level of variation in the PFam54 gene array with 39 PFam54 paralogs belonging to 23 orthologous groups including five novel paralogs. Even so, the gene array appears to have remained fairly stable over the evolutionary history of the studied Borrelia species. Interestingly, genes outside Clade IV, which contains genes encoding for proteins associated with Borrelia pathogenesis, more frequently displayed signatures of diversifying selection between clades that differ in hypothesized vector or host species. This could suggest that non-Clade IV paralogs play a more important role in host and/or vector adaptation than previously expected, which would require future lab-based studies to validate.

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.

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.

Genome diversity of Borrelia garinii in marine transmission cycles does not match host associations but reflects the strains evolutionary history

Borrelia burgdorferi sensu lato is a species complex of spirochetal bacteria that occupy different ecological niches which is reflected in their reservoir host- and vector-associations. Borrelia genomes possess numerous linear and circular plasmids. Proteins encoded by plasmid genes play a major role in host- and vector-interaction and are important for Borrelia niche adaptation. However, the plasmid composition and therewith the gene repertoire may vary even in strains of a single species. Borrelia garinii, one of the six human pathogenic species, is common in Europe (vector Ixodes ricinus), Asia (vector Ixodes persulcatus) and in marine birds (vector Ixodes uriae). For the latter, only a single culture isolate (Far04) and its genome were previously available. The genome was rather small containing only one circular and six linear plasmids with a notable absence of cp32 plasmids. To further investigate B. garinii from marine transmission cycles and to explore i) whether the small number of plasmids found in isolate Far04 is a common feature in B. garinii from marine birds and presents an adaptation to this particular niche and ii) whether there may be a correlation between genome type and host species, we initiated in vitro cultures from live I. uriae collected in 2017 and 2018 from marine avian hosts and their nests. Hosts included common guillemots, Atlantic Puffin, razorbill, and kittiwake. We obtained 17 novel isolates of which 10 were sequenced using Illumina technology, one also with Pacific Bioscience technology. The 10 genomes segregated into five different genome types defined by plasmid types (based on PFam32 loci). We show that the genomes of seabird associated B. garinii contain fewer plasmids (6-9) than B. garinii from terrestrial avian species (generally ≥10), potentially suggesting niche adaptation. However, genome type did not match an association with the diverse avian seabird hosts investigated but matched the clonal complex they originated from, perhaps reflecting the isolates evolutionary history. Questions that should be addressed in future studies are (i) how is plasmid diversity related to host- and/or vector adaptation; (ii) do the different seabird species differ in reservoir host competence, and (iii) can the genome types found in seabirds use terrestrial birds as reservoir hosts.

Evidence for Bartonella quintana in Lice Collected from the Clothes of Ethiopian Homeless Individuals

Human lice, Pediculus humanus, can transmit various pathogens, including Bartonella quintana, Borrelia recurrentis, and Rickettsia prowazekii. Xenosurveillance is an epidemiological approach to assessing human infection risks performed by screening vectors of infectious disease agents. In the proof-of-principle study reported herein, the DNA of 23 human lice was collected from the clothes of 30 homeless Ethiopian individuals. These samples were assessed using 16S rRNA gene-specific pan-eubacterial PCR for screening, followed by Bartonella genus 16S-23S internal transcribed spacer (ITS) sequence-specific PCR, Bartonella genus gltA gene-specific PCR, and 16S rRNA gene PCR with specificity for relapsing-fever-associated Borrelia spp. with subsequent sequencing of the amplicons. In one sample, the pan-eubacterial 16S rRNA gene-specific screening PCR, the Bartonella genus 16S-23S ITS sequence-specific PCR, and the Bartonella genus gltA gene-specific PCR allowed for the sequencing of B. quintana-specific amplicons. In two additional samples, Bartonella genus gltA gene-specific PCR also provided sequences showing 100% sequence identity with B. quintana. In total, 3/23 (13.0%) of the assessed lice were found to be positive for B. quintana. Correlating clinical data were not available; however, the assessment confirmed the presence of B. quintana in the local louse population and thus an associated infection pressure. Larger-sized cross-sectional studies seem advisable to more reliably quantify the infection risk of lice-infested local individuals. The need for prevention by providing opportunities to maintain standard hygiene for Ethiopian homeless individuals is stressed by the reported findings, especially in light of the ongoing migration of refugees.

Characterisation and comparative genomics of three new Varanus-associated Borrelia spp. from Indonesia and Australia

BACKGROUND

Borrelia are important disease-causing tick- and louse-borne spirochaetes than can infect a wide variety of vertebrates, including humans and reptiles. Reptile-associated (REP) Borrelia, once considered a peculiarity, are now recognised as a distinct and important evolutionary lineage, and are increasingly being discovered worldwide in association with novel hosts. Numerous novel Borrelia spp. associated with monitor lizards (Varanus spp.) have been recently identified throughout the Indo-Pacific region; however, there is a lack of genomic data on these Borrelia.

METHODS

We used metagenomic techniques to sequence almost complete genomes of novel Borrelia spp. from Varanus varius and Varanus giganteus from Australia, and used long- and short-read technologies to sequence the complete genomes of two strains of a novel Borrelia sp. previously isolated from ticks infesting Varanus salvator from Indonesia. We investigated intra- and interspecies genomic diversity, including plasmid diversity and relatedness, among Varanus-associated Borrelia and other available REP Borrelia and, based on 712 whole genome orthologues, produced the most complete phylogenetic analysis, to the best of our knowledge, of REP Borrelia to date.

RESULTS

The genomic architecture of Varanus-associated Borrelia spp. is similar to that of Borrelia spp. that cause relapsing fever (RF), and includes a highly conserved megaplasmid and numerous smaller linear and circular plasmids that lack structural consistency between species. Analysis of PF32 and PF57/62 plasmid partitioning genes indicated that REP Borrelia plasmids fall into at least six distinct plasmid families, some of which are related to previously defined Borrelia plasmid families, whereas the others appear to be unique. REP Borrelia contain immunogenic variable major proteins that are homologous to those found in Borrelia spp. that cause RF, although they are limited in copy number and variability and have low sequence identities to RF variable major proteins. Phylogenetic analyses based on single marker genes and 712 single copy orthologs also definitively demonstrated the monophyly of REP Borrelia as a unique lineage.

CONCLUSIONS

In this work we present four new genomes from three novel Borrelia, and thus double the number of REP Borrelia genomes publicly available. The genomic characterisation of these Borrelia clearly demonstrates their distinctiveness as species, and we propose the names Borrelia salvatorii, 'Candidatus Borrelia undatumii', and 'Candidatus Borrelia rubricentralis' for them.

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.

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.

Large-Scale Sequencing of Borreliaceae for the Construction of Pan-Genomic-Based Diagnostics

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.

Clostridium autoethanogenum isopropanol production via native plasmid pCA replicon

Clostridium autoethanogenum is a model gas-fermenting acetogen for commercial ethanol production. It is also a platform organism being developed for the carbon-negative production of acetone and isopropanol by gas fermentation. We have assembled a 5.5 kb pCA plasmid for type strain DSM10061 (JA1-1) using three genome sequence datasets. pCA is predicted to encode seven open-reading frames and estimated to be a low-copy number plasmid present at approximately 12 copies per chromosome. RNA-seq analyses indicate that pCA genes are transcribed at low levels and two proteins, CAETHG_05090 (putative replication protein) and CAETHG_05115 (hypothetical, a possible Mob protein), were detected at low levels during batch gas fermentations. Thiolase (thlA), CoA-transferase (ctfAB), and acetoacetate decarboxylase (adc) genes were introduced into a vector for isopropanol production in C. autoethanogenum using the native plasmid origin of replication. The availability of the pCA sequence will facilitate studies into its physiological role and could form the basis for genetic tool optimization.

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.

Lipoproteome screening of the Lyme disease agent identifies inhibitors of antibody-mediated complement killing

Spirochetal pathogens, such as the causative agent of Lyme disease, Borrelia burgdorferi sensu lato, encode an abundance of lipoproteins; however, due in part to their evolutionary distance from more well-studied bacteria, such as Proteobacteria and Firmicutes, few spirochetal lipoproteins have assigned functions. Indeed, B. burgdorferi devotes almost 8% of its genome to lipoprotein genes and interacts with its environment primarily through the production of at least 80 surface-exposed lipoproteins throughout its tick vector–vertebrate host lifecycle. Several B. burgdorferi lipoproteins have been shown to serve roles in cellular adherence or immune evasion, but the functions for most B. burgdorferi surface lipoproteins remain unknown. In this study, we developed a B. burgdorferi lipoproteome screening platform utilizing intact spirochetes that enables the identification of previously unrecognized host interactions. As spirochetal survival in the bloodstream is essential for dissemination, we targeted our screen to C1, the first component of the classical (antibody-initiated) complement pathway. We identified two high-affinity C1 interactions by the paralogous lipoproteins, ElpB and ElpQ (also termed ErpB and ErpQ, respectively). Using biochemical, microbiological, and biophysical approaches, we demonstrate that ElpB and ElpQ bind the activated forms of the C1 proteases, C1r and C1s, and represent a distinct mechanistic class of C1 inhibitors that protect the spirochete from antibody-mediated complement killing. In addition to identifying a mode of complement inhibition, our study establishes a lipoproteome screening methodology as a discovery platform for identifying direct host–pathogen interactions that are central to the pathogenesis of spirochetes, such as the Lyme disease agent.

Utilizing Two Borrelia bavariensis Isolates Naturally Lacking the PFam54 Gene Array To Elucidate the Roles of PFam54-Encoded Proteins

Lyme borreliosis is the most common vector-borne disease in the Northern Hemisphere, caused by spirochetes belonging to the Borrelia burgdorferi sensu lato species complex, which are transmitted by ixodid ticks. B. burgdorferi sensu lato species produce a family of proteins on the linear plasmid 54 (PFam54), some of which confer the functions of cell adhesion and inactivation of complement, the first line of host defense. However, the impact of PFam54 in promoting B. burgdorferi sensu lato pathogenesis remains unclear because of the hurdles to simultaneously knock out all PFam54 proteins in a spirochete. Here, we describe two Borrelia bavariensis strains, PBN and PNi, isolated from patients naturally lacking PFam54 but maintaining the rest of the genome with greater than 95% identity to the reference B. bavariensis strain, PBi. We found that PBN and PNi less efficiently survive in human serum than PBi. Such defects were restored by introducing two B. bavariensis PFam54 recombinant proteins, BGA66 and BGA71, confirming the role of these proteins in providing complement evasion of B. bavariensis. Further, we found that all three strains remain detectable in various murine tissues 21 days post-subcutaneous infection, supporting the nonessential role of B. bavariensis PFam54 in promoting spirochete persistence. This study identified and utilized isolates deficient in PFam54 to associate the defects with the absence of these proteins, building the foundation to further study the role of each PFam54 protein in contributing to B. burgdorferi sensu lato pathogenesis. IMPORTANCE To establish infections, Lyme borreliae utilize various means to overcome the host's immune system. Proteins encoded by the PFam54 gene array play a role in spirochete survival in vitro and in vivo. Moreover, this gene array has been described in all currently available Lyme borreliae genomes. By investigating the first two Borrelia bavariensis isolates naturally lacking the entire PFam54 gene array, we showed that both patient isolates display an increased susceptibility to human serum, which can be rescued in the presence of two PFam54 recombinant proteins. However, both isolates remain infectious to mice after intradermal inoculation, suggesting the nonessential role of PFam54 during the long-term, but may differ slightly in the colonization of specific tissues. Furthermore, these isolates show high genomic similarity to type strain PBi (>95%) and could be used in future studies investigating the role of each PFam54 protein in Lyme borreliosis pathogenesis.

Comparative Genomics Analyses Support the Reclassification of Bisgaard Taxon 40 as Mergibacter gen. nov., With Mergibacter septicus sp. nov. as Type Species: Novel Insights Into the Phylogeny and Virulence Factors of a Pasteurellaceae Family Member Associated With Mortality Events in Seabirds

The Pasteurellaceae family has been associated with fatal diseases in numerous avian species. Several new taxa within this family, including Bisgaard taxon 40, have been recently described in wild birds, but their genomic characteristics and pathogenicity are not well understood. We isolated Bisgaard taxon 40 from four species of seabirds, including one sampled during a mass, multi-species mortality event in Florida, United States. Here, we present a comprehensive phenotypic and genetic characterization of Bisgaard taxon 40 and comparative genomic analysis with reference strains from the Pasteurellaceae family, aiming at determining its phylogenetic position, antimicrobial susceptibility profile, and identifying putative virulence factors. In silico multilocus sequence-based and whole-genome-based phylogenetic analysis clustered all Bisgaard taxon 40 strains together on a distinct branch separated from the other members of the Pasteurellaceae family, indicating that Bisgaard taxon 40 could represent a new genus. These findings were further supported by protein similarity analyses using the concatenation of 31 conserved proteins and other taxonomic approaches such as the percentage of conserved protein test. Additionally, several putative virulence factors were identified, including those associated with adhesion (capsule, ompA, ompH) and colonization (exbD, fur, galU, galE, lpxA, lpxC, and kdsA) of the host and a cytolethal distending toxin (cdt), which may have played a role in disease development leading to the mortality event. Considerably low minimum inhibitory concentrations (MICs) were found for all the drugs tested, in concordance with the absence of antimicrobial resistance genes in these genomes. The novel findings of this study highlight genomic and phenotypic characteristics of this bacterium, providing insights into genome evolution and pathogenicity. We propose a reclassification of these organisms within the Pasteurellaceae family, designated as Mergibacter gen. nov., with Mergibacter septicus sp. nov. as the type species. The type strain is Mergibacter septicus A25201^T (=DSM 112696).

Recent Progress in Lyme Disease and Remaining Challenges

Lyme disease (also known as Lyme borreliosis) is the most common vector-borne disease in the United States with an estimated 476,000 cases per year. While historically, the long-term impact of Lyme disease on patients has been controversial, mounting evidence supports the idea that a substantial number of patients experience persistent symptoms following treatment. The research community has largely lacked the necessary funding to properly advance the scientific and clinical understanding of the disease, or to develop and evaluate innovative approaches for prevention, diagnosis, and treatment. Given the many outstanding questions raised into the diagnosis, clinical presentation and treatment of Lyme disease, and the underlying molecular mechanisms that trigger persistent disease, there is an urgent need for more support. This review article summarizes progress over the past 5 years in our understanding of Lyme and tick-borne diseases in the United States and highlights remaining challenges.

Host association of Borrelia burgdorferi sensu lato: A review

Borrelia burgdorferi sensu lato (Bbsl) is a bacterial species complex that includes the etiological agents of the most frequently reported vector-borne disease in the Northern hemisphere, Lyme borreliosis. It currently comprises > 20 named and proposed genospecies that use vertebrate hosts and tick vectors for transmission in the Americas and Eurasia. Host (and vector) associations influence geographic distribution and speciation in Bbsl, which is of particular relevance to human health. To target gaps in knowledge for future efforts to understand broad patterns of the Bbsl-tick-host system and how they relate to human health, the present review aims to give a comprehensive summary of the literature on host association in Bbsl. Of 465 papers consulted (404 after exclusion criteria were applied), 96 sought to experimentally establish reservoir competence of 143 vertebrate host species for Bbsl. We recognize xenodiagnosis as the strongest method used, however it is infrequent (20% of studies) probably due to difficulties in maintaining tick vectors and/or wild host species in the lab. Some well-established associations were not experimentally confirmed according to our definition (ex: Borrelia garinii, Ixodes uriae and sea birds). We conclude that our current knowledge on host association in Bbsl is mostly derived from a subset of host, vector and bacterial species involved, providing an incomplete knowledge of the physiology, ecology and evolutionary history of these interactions. More studies are needed on all host, vector and bacterial species globally involved with a focus on non-rodent hosts and Asian Bbsl complex species, especially with experimental research that uses xenodiagnosis and genomics to analyze existing host associations in different ecosystems.

Quantitative multiplexed strategies for human Lyme disease serological testing

Lyme disease, which is primarily caused by infection with the bacterium Borrelia burgdorferi in the United States or other Borrelia species internationally, presents an ongoing challenge for diagnostics. Serological testing is the primary means of diagnosis but testing approaches differ widely, with varying degrees of sensitivity and specificity. Moreover, there is currently no reliable test to determine disease resolution following treatment. A distinct challenge in Lyme disease diagnostics is the variable patterns of human immune response to a plurality of antigens presented by Borrelia spp. during the infection. Thus, multiplexed testing approaches that capture these patterns and detect serological response against multiple antigens may be the key to prompt, accurate Lyme disease diagnosis. In this review, current state-of-the-art multiplexed diagnostic approaches are presented and compared with respect to their diagnostic accuracy and their potential for monitoring response to treatment.

Antibiotic resistance plasmid composition and architecture in Escherichia coli isolates from meat

Resistance plasmids play a crucial role in the transfer of antimicrobial resistance from the veterinary sector to human healthcare. In this study plasmids from foodborne Escherichia coli isolates with a known (ES)BL or tetracycline resistance were sequenced entirely with short- and long-read technologies to obtain insight into their composition and to identify driving factors for spreading. Resistant foodborne E. coli isolates often contained several plasmids coding for resistance to various antimicrobials. Most plasmids were large and contained multiple resistance genes in addition to the selected resistance gene. The majority of plasmids belonged to the IncI, IncF and IncX incompatibility groups. Conserved and variable regions could be distinguished in each of the plasmid groups. Clusters containing resistance genes were located in the variable regions. Tetracycline and (extended spectrum) beta-lactamase resistance genes were each situated in separate clusters, but sulphonamide, macrolide and aminoglycoside formed one cluster and lincosamide and aminoglycoside another. In most plasmids, addiction systems were found to maintain presence in the cell.

Lyme Disease in Humans

Lyme disease (Lyme borreliosis) is a tick-borne, zoonosis of adults and children caused by genospecies of the Borrelia burgdorferi sensu lato complex. The ailment, widespread throughout the Northern Hemisphere, continues to increase globally due to multiple environmental factors, coupled with increased incursion of humans into habitats that harbor the spirochete. B. burgdorferi sensu lato is transmitted by ticks from the Ixodes ricinus complex. In North America, B. burgdorferi causes nearly all infections; in Europe, B. afzelii and B. garinii are most associated with human disease. The spirochete's unusual fragmented genome encodes a plethora of differentially expressed outer surface lipoproteins that play a seminal role in the bacterium's ability to sustain itself within its enzootic cycle and cause disease when transmitted to its incidental human host. Tissue damage and symptomatology (i.e., clinical manifestations) result from the inflammatory response elicited by the bacterium and its constituents. The deposition of spirochetes into human dermal tissue generates a local inflammatory response that manifests as erythema migrans (EM), the hallmark skin lesion. If treated appropriately and early, the prognosis is excellent. However, in untreated patients, the disease may present with a wide range of clinical manifestations, most commonly involving the central nervous system, joints, or heart. A small percentage (~10%) of patients may go on to develop a poorly defined fibromyalgia-like illness, post-treatment Lyme disease (PTLD) unresponsive to prolonged antimicrobial therapy. Below we integrate current knowledge regarding the ecologic, epidemiologic, microbiologic, and immunologic facets of Lyme disease into a conceptual framework that sheds light on the disorder that healthcare providers encounter.

Genetic Manipulation of Borrelia

Genetic studies in Borrelia require special consideration of the highly segmented genome, complex growth requirements and evolutionary distance of spirochetes from other genetically tractable bacteria. Despite these challenges, a robust molecular genetic toolbox has been constructed to investigate the biology and pathogenic potential of these important human pathogens. In this review we summarize the tools and techniques that are currently available for the genetic manipulation of Borrelia, including the relapsing fever spirochetes, viewing them in the context of their utility and shortcomings. Our primary objective is to help researchers discern what is feasible and what is not practical when thinking about potential genetic experiments in Borrelia. We have summarized published methods and highlighted their critical elements, but we are not providing detailed protocols. Although many advances have been made since B. burgdorferi was first transformed over 25 years ago, some standard genetic tools remain elusive for Borrelia. We mention these limitations and why they persist, if known. We hope to encourage investigators to explore what might be possible, in addition to optimizing what currently can be achieved, through genetic manipulation of Borrelia.

High conservation combined with high plasticity: genomics and evolution of Borrelia bavariensis

Background

Borrelia bavariensis is one of the agents of Lyme Borreliosis (or Lyme disease) in Eurasia. The genome of the Borrelia burgdorferi sensu lato species complex, that includes B. bavariensis, is known to be very complex and fragmented making the assembly of whole genomes with next-generation sequencing data a challenge.

Results

We present a genome reconstruction for 33 B. bavariensis isolates from Eurasia based on long-read (Pacific Bioscience, for three isolates) and short-read (Illumina) data. We show that the combination of both sequencing techniques allows proper genome reconstruction of all plasmids in most cases but use of a very close reference is necessary when only short-read sequencing data is available. B. bavariensis genomes combine a high degree of genetic conservation with high plasticity: all isolates share the main chromosome and five plasmids, but the repertoire of other plasmids is highly variable. In addition to plasmid losses and gains through horizontal transfer, we also observe several fusions between plasmids. Although European isolates of B. bavariensis have little diversity in genome content, there is some geographic structure to this variation. In contrast, each Asian isolate has a unique plasmid repertoire and we observe no geographically based differences between Japanese and Russian isolates. Comparing the genomes of Asian and European populations of B. bavariensis suggests that some genes which are markedly different between the two populations may be good candidates for adaptation to the tick vector, (Ixodes ricinus in Europe and I. persulcatus in Asia).

Conclusions

We present the characterization of genomes of a large sample of B. bavariensis isolates and show that their plasmid content is highly variable. This study opens the way for genomic studies seeking to understand host and vector adaptation as well as human pathogenicity in Eurasian Lyme Borreliosis agents.

Solution NMR structure of Borrelia burgdorferi outer surface lipoprotein BBP28, a member of the mlp protein family

Lyme disease is the most widespread vector-transmitted disease in North America and Europe, caused by infection with Borrelia burgdorferi sensu lato complex spirochetes. We report the solution NMR structure of the B. burgdorferi outer surface lipoprotein BBP28, a member of the multicopy lipoprotein (mlp) family. The structure comprises a tether peptide, five α-helices and an extended C-terminal loop. The fold is similar to that of Borrelia turicatae outer surface protein BTA121, which is known to bind lipids. These results contribute to the understanding of Lyme disease pathogenesis by revealing the molecular structure of a protein from the widely found mlp family.

Comparative Genomic Analysis of Food-Originated Coagulase-Negative Staphylococcus: Analysis of Conserved Core Genes and Diversity of the Pan-Genome

To shed light on the genetic differences among food-originated coagulase-negative Staphylococcus (CNS), we performed pan-genome analysis of five species: Staphylococcus carnosus (two strains), Staphylococcus equorum (two strains), Staphylococcus succinus (three strains), Staphylococcus xylosus (two strains), and Staphylococcus saprophyticus (one strain). The pan-genome size increases with each new strain and currently holds about 4,500 genes from 10 genomes. Specific genes were shown to be strain dependent but not species dependent. Most specific genes were of unknown function or encoded restriction-modification enzymes, transposases, or prophages. Our results indicate that unique genes have been acquired or lost by convergent evolution within individual strains.

Analysis of COMPASS, a New Comprehensive Plasmid Database Revealed Prevalence of Multireplicon and Extensive Diversity of IncF Plasmids

Plasmids are genetic elements that enable rapid adaptation and evolution by transferring genes conferring selective advantages to their hosts. Conjugative plasmids are predominantly responsible for the global dissemination of antimicrobial resistance, representing an important threat to global health. As the number of plasmid sequences grows exponentially, it becomes critical to depict the global diversity and decipher the distribution of circulating plasmids in the bacterial community. To this end, we created COMPASS, a novel and comprehensive database compiling 12,084 complete plasmids with associated metadata from 1571 distinct species isolated worldwide over more than 100 years. The curation of the database allowed us to identify identical plasmids across different bacteria revealing mainly intraspecies dissemination and rare cases of horizontal transmission. We outlined and analyzed all relevant features, plasmid properties, host range and characterized their replication and mobilization systems. After an exhaustive comparison of PlasmidFinder and MOB-typer, the MOB-typer-based analysis revealed that the current knowledge embedded in the current typing schemes fails to classify all the plasmid sequences collected in COMPASS. We were able to categorize 6828 and 5229 plasmids by replicon and MOB typing, respectively, mostly associated with Proteobacteria and Firmicutes. We then searched for the presence of multiple core genes involved in replication and propagation. Our results showed that 2403 plasmids carried multiple replicons that were distributed in 206 bacterial species. The co-integration of replicon types from different incompatibility (Inc) groups is an adaptive mechanism, which plays an important role in plasmid survival and dissemination by extending their host range. Our results highlight the crucial role of IncF alleles (present in 56% of all multireplicons) and revealed that IncH, IncR, and IncU replicons were also frequently carried in multireplicons. Here, we provided a comprehensive picture of the different IncF subtypes by identifying 20 different profiles in 849 IncF multireplicons, which were mostly associated with Enterobacteriaceae. These results could provide the basis for a novel IncF plasmid nomenclature based on different allelic profiles.

Borrelia maritima sp. nov., a novel species of the Borrelia burgdorferi sensu lato complex, occupying a basal position to North American species

Borrelia species are vector-borne parasitic bacteria with unusual, highly fragmented genomes that include a linear chromosome and linear as well as circular plasmids that differ numerically between and within various species. Strain CA690^T, which was cultivated from a questing Ixodes spinipalpis nymph in the San Francisco Bay area, CA, was determined to be genetically distinct from all other described species belonging to the Borrelia burgdorferi sensu lato complex. The genome, including plasmids, was assembled using a hybrid assembly of short Illumina reads and long reads obtained via Oxford Nanopore Technology. We found that strain CA690^T has a main linear chromosome containing 902176 bp with a blast identity ≤91 % compared with other Borrelia species chromosomes and five linear and two circular plasmids. A phylogeny based on 37 single-copy genes of the main linear chromosome and rooted with the relapsing fever species Borrelia duttonii strain Ly revealed that strain CA690^T had a sister-group relationship with, and occupied a basal position to, species occurring in North America. We propose to name this species Borrelia maritima sp. nov. The type strain, CA690^T, has been deposited in two national culture collections, DSMZ (=107169) and ATCC (=TSD-160).

Changing of the guard: How the Lyme disease spirochete subverts the host immune response

Lyme disease, also known as Lyme borreliosis, is the most common tick-transmitted disease in the Northern Hemisphere. The disease is caused by the bacterial spirochete Borrelia burgdorferi and other related Borrelia species. One of the many fascinating features of this unique pathogen is an elaborate system for antigenic variation, whereby the sequence of the surface-bound lipoprotein VlsE is continually modified through segmental gene conversion events. This perpetual changing of the guard allows the pathogen to remain one step ahead of the acquired immune response, enabling persistent infection. Accordingly, the vls locus is the most evolutionarily diverse genetic element in Lyme disease-causing borreliae. Small stretches of information are transferred from a series of silent cassettes in the vls locus to generate an expressed mosaic vlsE gene version that contains genetic information from several different silent cassettes, resulting in ∼1040 possible vlsE sequences. Yet, despite its extreme evolutionary flexibility, the locus has rigidly conserved structural features. These include a telomeric location of the vlsE gene, an inverse orientation of vlsE and the silent cassettes, the presence of nearly perfect inverted repeats of ∼100 bp near the 5' end of vlsE, and an exceedingly high concentration of G runs in vlsE and the silent cassettes. We discuss the possible roles of these evolutionarily conserved features, highlight recent findings from several studies that have used next-generation DNA sequencing to unravel the switching process, and review advances in the development of a mini-vls system for genetic manipulation of the locus.

Whole genome sequencing of Borrelia miyamotoi isolate Izh-4: reference for a complex bacterial genome

BACKGROUND

The genus Borrelia comprises spirochaetal bacteria maintained in natural transmission cycles by tick vectors and vertebrate reservoir hosts. The main groups are represented by a species complex including the causative agents of Lyme borreliosis and relapsing fever group Borrelia. Borrelia miyamotoi belongs to the relapsing fever group of spirochetes and forms distinct populations in North America, Asia, and Europe. As all Borrelia species B. miyamotoi possess an unusual and complex genome consisting of a linear chromosome and a number of linear and circular plasmids. The species is considered an emerging human pathogen and an increasing number of human cases are being described in the Northern hemisphere. The aim of this study was to produce a high quality reference genome that will facilitate future studies into genetic differences between different populations and the genome plasticity of B. miyamotoi.

RESULTS

We used multiple available sequencing methods, including Pacific Bioscience single-molecule real-time technology (SMRT) and Oxford Nanopore technology (ONT) supplemented with highly accurate Illumina sequences, to explore the suitability for whole genome assembly of the Russian B. miyamotoi isolate, Izh-4. Plasmids were typed according to their potential plasmid partitioning genes (PF32, 49, 50, 57/62). Comparing and combining results of both long-read (SMRT and ONT) and short-read methods (Illumina), we determined that the genome of the isolate Izh-4 consisted of one linear chromosome, 12 linear and two circular plasmids. Whilst the majority of plasmids had corresponding contigs in the Asian B. miyamotoi isolate FR64b, there were only four that matched plasmids of the North American isolate CT13-2396, indicating differences between B. miyamotoi populations. Several plasmids, e.g. lp41, lp29, lp23, and lp24, were found to carry variable major proteins. Amongst those were variable large proteins (Vlp) subtype Vlp-α, Vlp-γ, Vlp-δ and also Vlp-β. Phylogenetic analysis of common plasmids types showed the uniqueness in Russian/Asian isolates of B. miyamotoi compared to other isolates.

CONCLUSIONS

We here describe the genome of a Russian B. miyamotoi clinical isolate, providing a solid basis for future comparative genomics of B. miyamotoi isolates. This will be a great impetus for further basic, molecular and epidemiological research on this emerging tick-borne pathogen.

Lyme Disease Frontiers: Reconciling Borrelia Biology and Clinical Conundrums

Lyme disease is a complex tick-borne zoonosis that poses an escalating public health threat in several parts of the world, despite sophisticated healthcare infrastructure and decades of effort to address the problem. Concepts like the true burden of the illness, from incidence rates to longstanding consequences of infection, and optimal case management, also remain shrouded in controversy. At the heart of this multidisciplinary issue are the causative spirochetal pathogens belonging to the Borrelia Lyme complex. Their unusual physiology and versatile lifestyle have challenged microbiologists, and may also hold the key to unlocking mysteries of the disease. The goal of this review is therefore to integrate established and emerging concepts of Borrelia biology and pathogenesis, and position them in the broader context of biomedical research and clinical practice. We begin by considering the conventions around diagnosing and characterizing Lyme disease that have served as a conceptual framework for the discipline. We then explore virulence from the perspective of both host (genetic and environmental predispositions) and pathogen (serotypes, dissemination, and immune modulation), as well as considering antimicrobial strategies (lab methodology, resistance, persistence, and clinical application), and borrelial adaptations of hypothesized medical significance (phenotypic plasticity or pleomorphy).

Rapid, multiplexed, whole genome and plasmid sequencing of foodborne pathogens using long-read nanopore technology

U.S. public health agencies have employed next-generation sequencing (NGS) as a tool to quickly identify foodborne pathogens during outbreaks. Although established short-read NGS technologies are known to provide highly accurate data, long-read sequencing is still needed to resolve highly-repetitive genomic regions and genomic arrangement, and to close the sequences of bacterial chromosomes and plasmids. Here, we report the use of long-read nanopore sequencing to simultaneously sequence the entire chromosome and plasmid of Salmonella enterica subsp. enterica serovar Bareilly and Escherichia coli O157:H7. We developed a rapid and random sequencing approach coupled with de novo genome assembly within a customized data analysis workflow that uses publicly-available tools. In sequencing runs as short as four hours, using the MinION instrument, we obtained full-length genomes with an average identity of 99.87% for Salmonella Bareilly and 99.89% for E. coli in comparison to the respective MiSeq references. These nanopore-only assemblies provided readily available information on serotype, virulence factors, and antimicrobial resistance genes. We also demonstrate the potential of nanopore sequencing assemblies for rapid preliminary phylogenetic inference. Nanopore sequencing provides additional advantages as very low capital investment and footprint, and shorter (10 hours library preparation and sequencing) turnaround time compared to other NGS technologies.

Genomic analysis on broiler-associated Clostridium perfringens strains and exploratory caecal microbiome investigation reveals key factors linked to poultry necrotic enteritis

Background

Clostridium perfringens is a key pathogen in poultry-associated necrotic enteritis (NE). To date there are limited Whole Genome Sequencing based studies describing broiler-associated C. perfringens in healthy and diseased birds. Moreover, changes in the caecal microbiome during NE is currently not well characterised. Thus, the aim of this present study was to investigate C. perfringens virulence factors linked to health and diseased chickens, including identifying putative caecal microbiota signatures associated with NE.

Results

We analysed 88 broiler chicken C. perfringens genomes (representing 66 publicly available genomes and 22 newly sequenced genomes) using different phylogenomics approaches and identified a potential hypervirulent and globally-distributed clone spanning 20-year time-frame (1993-2013). These isolates harbored a greater number of virulence genes (including toxin and collagen adhesin genes) when compared to other isolates. Further genomic analysis indicated exclusive and overabundant presence of important NE-linked toxin genes including netB and tpeL in NE-associated broiler isolates. Secondary virulence genes including pfoA, cpb2, and collagen adhesin genes cna, cnaA and cnaD were also enriched in the NE-linked C. perfringens genomes. Moreover, an environmental isolate obtained from farm animal feeds was found to encode netB, suggesting potential reservoirs of NetB-positive C. perfringens strains (toxinotype G). We also analysed caecal samples from a small sub-set of 11 diseased and healthy broilers for exploratory microbiome investigation using 16S rRNA amplicon sequencing, which indicated a significant and positive correlation in genus Clostridium within the wider microbiota of those broilers diagnosed with NE, alongside reductions in beneficial microbiota members.

Conclusions

These data indicate a positive association of virulence genes including netB, pfoA, cpb2, tpeL and cna variants linked to NE-linked isolates. Potential global dissemination of specific hypervirulent lineage, coupled with distinctive microbiome profiles, highlights the need for further investigations, which will require a large worldwide sample collection from healthy and NE-associated birds.

Complement Evasion in Borrelia spirochetes: Mechanisms and Opportunities for Intervention

Lyme disease (LD) is an increasingly prevalent, climate change-accelerated, vector-borne infectious disease with significant morbidity and cost in a proportion of patients who experience ongoing symptoms after antibiotic treatment, a condition known as post-treatment Lyme disease syndrome (PTLDS). Spirochetal bacteria of Borrelia species are the causative agents of LD. These obligate parasites have evolved sophisticated immune evasion mechanisms, including the ability to defeat the innate immune system’s complement cascade. Research on complement function and Borrelia evasion mechanisms, focusing on human disease, is reviewed, highlighting opportunities to build on existing knowledge. Implications for the development of new antibiotic therapies having the potential to prevent or cure PTLDS are discussed. It is noted that a therapy enabling the complement system to effectively counter Borrelia might have lower cost and fewer side-effects and risks than broad-spectrum antibiotic use and could avert the need to develop and administer a vaccine.

DksA Controls the Response of the Lyme Disease Spirochete Borrelia burgdorferi to Starvation

The pathogenic spirochete Borrelia burgdorferi senses and responds to changes in the environment, including changes in nutrient availability, throughout its enzootic cycle in Ixodes ticks and vertebrate hosts. This study examined the role of DnaK suppressor protein (DksA) in the transcriptional response of B. burgdorferi to starvation. Wild-type and dksA mutant B. burgdorferi strains were subjected to starvation by shifting cultures grown in rich complete medium, Barbour-Stoenner-Kelly II (BSK II) medium, to a defined mammalian tissue culture medium, RPMI 1640, for 6 h under microaerobic conditions (5% CO₂, 3% O₂). Microarray analyses of wild-type B. burgdorferi revealed that genes encoding flagellar components, ribosomal proteins, and DNA replication machinery were downregulated in response to starvation. DksA mediated transcriptomic responses to starvation in B. burgdorferi, as the dksA-deficient strain differentially expressed only 47 genes in response to starvation compared to the 500 genes differentially expressed in wild-type strains. Consistent with a role for DksA in the starvation response of B. burgdorferi, fewer CFU of dksA mutants were observed after prolonged starvation in RPMI 1640 medium than CFU of wild-type B. burgdorferi spirochetes. Transcriptomic analyses revealed a partial overlap between the DksA regulon and the regulon of Rel_Bbu, the guanosine tetraphosphate and guanosine pentaphosphate [(p)ppGpp] synthetase that controls the stringent response; the DksA regulon also included many plasmid-borne genes. Additionally, the dksA mutant exhibited constitutively elevated (p)ppGpp levels compared to those of the wild-type strain, implying a regulatory relationship between DksA and (p)ppGpp. Together, these data indicate that DksA, along with (p)ppGpp, directs the stringent response to effect B. burgdorferi adaptation to its environment.IMPORTANCE The Lyme disease bacterium Borrelia burgdorferi survives diverse environmental challenges as it cycles between its tick vectors and various vertebrate hosts. B. burgdorferi must withstand prolonged periods of starvation while it resides in unfed Ixodes ticks. In this study, the regulatory protein DksA is shown to play a pivotal role controlling the transcriptional responses of B. burgdorferi to starvation. The results suggest that DksA gene regulatory activity impacts B. burgdorferi metabolism, virulence gene expression, and the ability of this bacterium to complete its natural life cycle.

Core genome phylogenetic analysis of the avian associated Borrelia turdi indicates a close relationship to Borrelia garinii

Borrelia burgdorferi sensu lato comprises a species complex of tick-transmitted bacteria that includes the agents of human Lyme borreliosis. Borrelia turdi is a genospecies of this complex that exists in cryptic transmission cycles mainly between ornithophilic tick vectors and their avian hosts. The species has been originally discovered in avian transmission cycles in Asia but has increasingly been found in Europe. Next generation sequencing was used to sequence the genome of B. turdi isolates obtained from ticks feeding on birds in Portugal to better understand the evolution and phylogenetic relationship of this avian and ornithophilic tick-associated genospecies. Here we use draft genomes of these B. turdi isolates for comparative analysis and to determine the taxonomic position within the B. burgdorferi s.l. species complex. The main chromosomes showed a maximum similarity of 93% to other Borrelia species whilst most plasmids had lower similarities. All three isolates had nine or 10 plasmids and, interestingly, one plasmid with a novel partitioning protein; this plasmid was termed lp30. Phylogenetic analysis of multilocus sequence typing housekeeping genes and 113 single copy orthologous genes revealed that the isolates clustered according to their classification as B. turdi. In phylogenies generated from these 113 genes the isolates cluster together with other Eurasian genospecies and form a sister clade to the avian associated B. garinii and the rodent associated B. bavariensis. These findings show that Borrelia species maintained in cryptic ecological cycles need to be included to fully understand the complex ecology and evolutionary history of this bacterial species complex.

Genome-wide analysis of Borrelia turcica and 'Candidatus Borrelia tachyglossi' shows relapsing fever-like genomes with unique genomic links to Lyme disease Borrelia

Borrelia are tick-borne bacteria that in humans are the aetiological agents of Lyme disease and relapsing fever. Here we present the first genomes of B. turcica and B. tachyglossi, members of a recently described and rapidly expanding Borrelia clade associated with reptile (B. turcica) or echidna (B. tachyglossi) hosts, transmitted by hard ticks, and of unknown pathogenicity. Borrelia tachyglossi and B. turcica genomes are similar to those of relapsing fever Borrelia species, containing a linear ~ 900 kb chromosome, a single long (> 70 kb) linear plasmid, and numerous short (< 40 kb) linear and circular plasmids, as well as a suite of housekeeping and macronutrient biosynthesis genes which are not found in Lyme disease Borrelia. Additionally, both B. tachyglossi and B. turcica contain paralogous vsp and vlp proteins homologous to those used in the multiphasic antigen-switching system used by relapsing fever Borrelia to evade vertebrate immune responses, although their number was greatly reduced compared to human-infectious species. However, B. tachyglossi and B. turcica chromosomes also contain numerous genes orthologous to Lyme disease Borrelia-specific genes, demonstrating a unique evolutionary, and potentially phenotypic link between these groups. Borrelia tachyglossi and B. turcica genomes also have unique genetic features, including degraded and deleted tRNA modification genes, and an expanded range of macronutrient salvage and biosynthesis genes compared to relapsing fever and Lyme disease Borrelia. These genomes and genomic comparisons provide an insight into the biology and evolutionary origin of these Borrelia, and provide a valuable resource for future work.

High Throughput Sequencing and Network Analysis Disentangle the Microbial Communities of Ticks and Hosts Within and Between Ecosystems

We aimed to develop a framework, based on graph theory, to capture the ecological meaning behind pure pair comparisons of microbiome-derived data. As a proof of concept, we applied the framework to analyze the co-occurrence of bacteria in either Ixodes ricinus ticks or the spleen of one of their main hosts, the vole Myodes glareolus. As a secondary lymphoid organ, the spleen acts as a filter of blood and represents well the exposure to microorganisms circulating in the blood; including those acquired and transmitted by ticks during feeding. The microbiome of 301 and 269 individual tick and vole samples, respectively, were analyzed using next generation sequencing (NGS) of 16S rRNA. To assess the effect of habitat on ecological communities of bacteria associated to ticks and voles, two different biotopes were included in the study, forest, and ecotone. An innovative approach of NGS data analysis combining network analysis and phylogenies of co-occuring of bacteria was used to study associations between bacteria in individual samples. Of the 126 bacterial genera found in ticks and voles, 62% were shared by both species. Communities of co-occurring bacteria were always more phylogenetically diverse in ticks than in voles. Interestingly, ~80% of bacterial phylogenetic diversity was found in ~20% of ticks. This pattern was not observed in vole-associated bacteria. Results revealed that the microbiome of I. ricinus is only slightly related to that of M. glareolus and that the biotope plays the most important role in shaping the bacterial communities of either ticks or voles. The analysis of the phylogenetic signal of the network indexes across the 16S rRNA-derived tree of bacteria suggests that the microbiome of both ticks and voles has high phylogenetic diversity and that closest bacterial genera do not co-occur. This study shows that network analysis is a promising tool to unravel complex microbial communities associated to arthropod vectors and vertebrate hosts.

High Throughput Sequencing and Network Analysis Disentangle the Microbial Communities of Ticks and Hosts Within and Between Ecosystems

We aimed to develop a framework, based on graph theory, to capture the ecological meaning behind pure pair comparisons of microbiome-derived data. As a proof of concept, we applied the framework to analyze the co-occurrence of bacteria in either Ixodes ricinus ticks or the spleen of one of their main hosts, the vole Myodes glareolus. As a secondary lymphoid organ, the spleen acts as a filter of blood and represents well the exposure to microorganisms circulating in the blood; including those acquired and transmitted by ticks during feeding. The microbiome of 301 and 269 individual tick and vole samples, respectively, were analyzed using next generation sequencing (NGS) of 16S rRNA. To assess the effect of habitat on ecological communities of bacteria associated to ticks and voles, two different biotopes were included in the study, forest, and ecotone. An innovative approach of NGS data analysis combining network analysis and phylogenies of co-occuring of bacteria was used to study associations between bacteria in individual samples. Of the 126 bacterial genera found in ticks and voles, 62% were shared by both species. Communities of co-occurring bacteria were always more phylogenetically diverse in ticks than in voles. Interestingly, ~80% of bacterial phylogenetic diversity was found in ~20% of ticks. This pattern was not observed in vole-associated bacteria. Results revealed that the microbiome of I. ricinus is only slightly related to that of M. glareolus and that the biotope plays the most important role in shaping the bacterial communities of either ticks or voles. The analysis of the phylogenetic signal of the network indexes across the 16S rRNA-derived tree of bacteria suggests that the microbiome of both ticks and voles has high phylogenetic diversity and that closest bacterial genera do not co-occur. This study shows that network analysis is a promising tool to unravel complex microbial communities associated to arthropod vectors and vertebrate hosts.

Whole genome sequence and comparative analysis of Borrelia burgdorferi MM1

Lyme disease is caused by spirochaetes of the Borrelia burgdorferi sensu lato genospecies. Complete genome assemblies are available for fewer than ten strains of Borrelia burgdorferi sensu stricto, the primary cause of Lyme disease in North America. MM1 is a sensu stricto strain originally isolated in the midwestern United States. Aside from a small number of genes, the complete genome sequence of this strain has not been reported. Here we present the complete genome sequence of MM1 in relation to other sensu stricto strains and in terms of its Multi Locus Sequence Typing. Our results indicate that MM1 is a new sequence type which contains a conserved main chromosome and 15 plasmids. Our results include the first contiguous 28.5 kb assembly of lp28-8, a linear plasmid carrying the vls antigenic variation system, from a Borrelia burgdorferi sensu stricto strain.

Laboratory diagnosis of Lyme borreliosis: Current state of the art and future perspectives

This review is directed at physicians and laboratory personnel in private practice and clinics who treat and diagnose Lyme borreliosis (LB) in patients as part of their daily work. A major objective of this paper is to bring together background information on Borrelia (B.) burgdorferi sensu lato (s.l.) and basic clinical knowledge of LB, which is one of the most frequently reported vector-borne diseases in the Northern Hemisphere. The goal is to provide practical guidance for clinicians and for laboratory physicians, and scientists for a better understanding of current achievements and ongoing obstacles in the laboratory diagnosis of LB, an infectious disease that still remains one of the diagnostic chameleons of modern clinical medicine. Moreover, in bringing together current scientific information from guidelines, reviews, and original papers, this review provides recommendations for selecting the appropriate tests in relation to the patient's stage of disease to achieve effective, stage-related application of current direct and indirect laboratory methods for the detection of B. burgdorferi s.l. Additionally, the review aims to discuss the current state of the art concerning the diagnostic potential and limitations of the assays and test methods currently in use to optimize LB patient management and provide insight into the possible future prospects of this rapidly changing area of laboratory medicine.

Primordial origin and diversification of plasmids in Lyme disease agent bacteria

BACKGROUND

With approximately one-third of their genomes consisting of linear and circular plasmids, the Lyme disease agent cluster of species has the most complex genomes among known bacteria. We report here a comparative analysis of plasmids in eleven Borreliella (also known as Borrelia burgdorferi sensu lato) species.

RESULTS

We sequenced the complete genomes of two B. afzelii, two B. garinii, and individual B. spielmanii, B. bissettiae, B. valaisiana and B. finlandensis isolates. These individual isolates carry between seven and sixteen plasmids, and together harbor 99 plasmids. We report here a comparative analysis of these plasmids, along with 70 additional Borreliella plasmids available in the public sequence databases. We identify only one new putative plasmid compatibility type (the 30th) among these 169 plasmid sequences, suggesting that all or nearly all such types have now been discovered. We find that the linear plasmids in the non-B. burgdorferi species have undergone the same kinds of apparently random, chaotic rearrangements mediated by non-homologous recombination that we previously discovered in B. burgdorferi. These rearrangements occurred independently in the different species lineages, and they, along with an expanded chromosomal phylogeny reported here, allow the identification of several whole plasmid transfer events among these species. Phylogenetic analyses of the plasmid partition genes show that a majority of the plasmid compatibility types arose early, most likely before separation of the Lyme agent Borreliella and relapsing fever Borrelia clades, and this, with occasional cross species plasmid transfers, has resulted in few if any species-specific or geographic region-specific Borreliella plasmid types.

CONCLUSIONS

The primordial origin and persistent maintenance of the Borreliella plasmid types support their functional indispensability as well as evolutionary roles in facilitating genome diversity. The improved resolution of Borreliella plasmid phylogeny based on conserved partition-gene clusters will lead to better determination of gene orthology which is essential for prediction of biological function, and it will provide a basis for inferring detailed evolutionary mechanisms of Borreliella genomic variability including homologous gene and plasmid exchanges as well as non-homologous rearrangements.