Trans-Atlantic exchanges have shaped the population structure of the Lyme disease agent Borrelia burgdorferi sensu stricto

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.

Case report of disseminated borrelial lymphocytoma with isolation of Borrelia burgdorferi sensu stricto in chronic lymphatic leukemia stage Binet A-an 11 year follow up

We report a rare manifestation of cutaneous borreliosis in a patient with pre-existing malignant lymphoproliferative disease, in particular chronic lymphocytic B cell leukemia (B-CLL). The patient's cutaneous lesions were initially diagnosed histologically as leukemia cutis. Distribution pattern of the skin lesions were in typical localizations for borrelial lymphocytoma. Borrelia burgdorferi sensu stricto was isolated and cultured from two sites (ear, mammilla). Antibiotic therapy improved the cutaneous lesions and the general condition of the patient. However, a second round of antibiotic therapy was required to resolve the lesions. At eleven years of follow-up the patient's skin was clear and she still had a stable condition of B-CLL without chemotherapy. In conclusion, the patient suffered from Lyme borreliosis (Borrelia lymphocytoma) and the cutaneous symptoms were aggravated by the underlying condition of chronic B-CLL condition.

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.

On the road to genomically defining bacterial intra-species units

Over almost three decades, average nucleotide identity (ANI) analysis has been instrumental in operationally defining species in bacteria. However, barely any attention has been paid to soundly defining intra-species units employing ANI analyses until recently. Notably, some very recent publications are good steps forward in that direction. The level of granularity provided by these intra-species units will be relevant to understanding the eco-evolutionary dynamics and transmission of bacterial lineages and mobile genetic elements, antibiotic resistance, and virulence genes. These intra-species units will undoubtedly advance the genomic epidemiology of many bacterial pathogens. In the coming years, we anticipate that many studies will implement ANI-based definitions of different intra-species units, such as strains or sequence types, for many different bacterial species.

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.

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.

Phylogenomic Diversity Elucidates Mechanistic Insights into Lyme Borreliae-Host Association

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.

Accessory Genomic Epidemiology of Cocirculating Acinetobacter baumannii Clones

Acinetobacter baumannii has become one of the most important multidrug-resistant nosocomial pathogens all over the world. Nonetheless, very little is known about the diversity of A. baumannii lineages coexisting in hospital settings. Here, using whole-genome sequencing, epidemiological data, and antimicrobial susceptibility tests, we uncover the transmission dynamics of extensive and multidrug-resistant A. baumannii in a tertiary hospital over a decade. Our core genome phylogeny of almost 300 genomes suggests that there were several introductions of lineages from international clone 2 into the hospital. The molecular dating analysis shows that these introductions happened in 2006, 2007, and 2013. Furthermore, using the accessory genome, we show that these lineages were extensively disseminated across many wards in the hospital. Our results demonstrate that accessory genome variation can be a very powerful tool for conducting genomic epidemiology. We anticipate future studies employing the accessory genome along with the core genome as a powerful phylogenomic strategy to track bacterial transmissions over very short microevolutionary scales. IMPORTANCE Whole-genome sequencing for epidemiological investigations (genomic epidemiology) has been of paramount importance to understand the transmission dynamics of many bacterial (and nonbacterial) pathogens. Commonly, variation in the core genome, single nucleotide polymorphisms (SNPs), is employed to carry out genomic epidemiology. However, at very short periods of time, the core genome might not have accumulated enough variation (sufficient SNPs) to tell apart isolates. In this scenario, gene content variation in the accessory genome can be an option to conduct genomic epidemiology. Here, we used the accessory genome, as well as the core genome, to uncover the transmission dynamics of extensive and multidrug-resistant A. baumannii in a tertiary hospital for a decade. Our study shows that accessory genome variation can be a very powerful tool for conducting genomic epidemiology.

The Population Structure of Borrelia lusitaniae Is Reflected by a Population Division of Its Ixodes Vector

Populations of vector-borne pathogens are shaped by the distribution and movement of vector and reservoir hosts. To study what impact host and vector association have on tick-borne pathogens, we investigated the population structure of Borrelia lusitaniae using multilocus sequence typing (MLST). Novel sequences were acquired from questing ticks collected in multiple North African and European locations and were supplemented by publicly available sequences at the Borrelia Pubmlst database (accessed on 11 February 2020). Population structure of B. lusitaniae was inferred using clustering and network analyses. Maximum likelihood phylogenies for two molecular tick markers (the mitochondrial 16S rRNA locus and a nuclear locus, Tick-receptor of outer surface protein A, trospA) were used to confirm the morphological species identification of collected ticks. Our results confirmed that B. lusitaniae does indeed form two distinguishable populations: one containing mostly European samples and the other mostly Portuguese and North African samples. Of interest, Portuguese samples clustered largely based on being from north (European) or south (North African) of the river Targus. As two different Ixodes species (i.e., I. ricinus and I. inopinatus) may vector Borrelia in these regions, reference samples were included for I. inopinatus but did not form monophyletic clades in either tree, suggesting some misidentification. Even so, the trospA phylogeny showed a monophyletic clade containing tick samples from Northern Africa and Portugal south of the river Tagus suggesting a population division in Ixodes on this locus. The pattern mirrored the clustering of B. lusitaniae samples, suggesting a potential co-evolution between tick and Borrelia populations that deserve further investigation.

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.

Origin of OXA-23 Variant OXA-239 from a Recently Emerged Lineage of Acinetobacter baumannii International Clone V

A. baumannii is a major cause of nosocomial infections all over the world. Although many isolates from developed countries have been studied in terms of their genome sequence, isolates from Latin America have been much less studied. In this study, using a population genomics approach considering the whole genomes of 148 isolates, we describe the recent emergence of the lineage ST758 endemic to Latin America and the inception of the OXA-239 carbapenemase. Our study highlights the urgent need to investigate recently emerged lineages of this species in Latin America and elsewhere, as these might harbor novel antibiotic resistance genes.

Over the last few decades, carbapenemase-producing Acinetobacter baumannii has become a major cause of nosocomial infections all over the world. However, the genome identity of lineages of this species in Latin America has not been studied as much as in developed countries. Here, through a population genomics approach considering the whole genomes of 148 isolates (almost 40 from Mexico and Honduras), we describe the recent emergence of the lineage sequence type 758 (ST758), which belongs to the international clone V and has spread out to Canada, Mexico, Honduras, and Colombia. Notably, this lineage was found to coexist with other A. baumannii lineages in hospitals in Mexico and Honduras. Isolates from this lineage show considerable variation in antibiotic resistance profiles, but most of them are resistant to carbapenems. Moreover, we found a variety of acquired oxacillinase (OXA) families within this lineage and tracked the very recent inception, and subsequent horizontal transmission, of the OXA-239 carbapenemase. This work highlights the urgent need to investigate recently emerged lineages of this species in Latin America and elsewhere, as these might harbor novel antibiotic resistance genes.

IMPORTANCEA. baumannii is a major cause of nosocomial infections all over the world. Although many isolates from developed countries have been studied in terms of their genome sequence, isolates from Latin America have been much less studied. In this study, using a population genomics approach considering the whole genomes of 148 isolates, we describe the recent emergence of the lineage ST758 endemic to Latin America and the inception of the OXA-239 carbapenemase. Our study highlights the urgent need to investigate recently emerged lineages of this species in Latin America and elsewhere, as these might harbor novel antibiotic resistance genes.

Identification of Borrelia bissettii in Ixodes scapularis ticks from New Brunswick, Canada

Lyme disease is a tick-borne disease that is emerging in Canada. The disease is caused by spirochetes of the Lyme borreliosis group, which is expanding as new species are discovered. In Canada, Lyme disease risk has so far been assessed primarily by detection of Borrelia burgdorferi sensu stricto. Of Ixodes scapularis ticks collected between 2014 and 2016 in New Brunswick, Canada, 7 were shown to be infected with Borrelia bissettii by nested PCR and sequencing of 5 B. bissettii genes. Since different Borrelia species are associated with different clinical manifestations and are not detected with the same diagnostic tests, the identification of a previously undocumented or underreported pathogenic Borrelia species has important implications for public and veterinary medicine.

Whole genome sequencing and phylogenetic analysis of strains of the agent of Lyme disease Borrelia burgdorferi from Canadian emergence zones

Lyme disease is emerging in southern Canada due to range expansion of the tick vector, followed by invasion of the agent of Lyme disease Borrelia burgdorferi sensu stricto. Strain diversity, as determined by Multi Locus Sequence Typing, occurs in this zone of emergence, and this may have its origins in adaptation to ecological niches, and have phenotypic consequences for pathogenicity and serological test performance. Sixty-four unique strains were cultured from ticks collected in southern Canada and the genomes sequenced using the Illumina MiSeq platform. A maximum likelihood phylogenetic tree of the chromosome revealed two large clades with multiple subclades. Consistent with previous studies on this species, the clades were not geographically defined, and some Canadian strains were highly divergent from previously sequenced US strains. There was evidence for recombination in the chromosome but this did not affect the phylogeny. Analysis of chromosomal genes indicated that these are under intense purifying selection. Phylogenies of the accessory genome and chromosome were congruent. Therefore strain differences identified in the phylogeny of chromosomal genes likely act as a proxy for genetic determinants of phenotypic differences amongst strains that are harboured in the accessory genome. Further studies on health implications of strain diversity are needed.

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.

Association of Borrelia and Rickettsia spp. and bacterial loads in Ixodes ricinus ticks

In recent years, awareness of coinfections has increased as synergistic or antagonistic effects on interacting bacteria have been observed. To date, several reports on coinfections of ticks with Rickettsia and Borrelia spp. are available. However, associations are rarely described and studies are based on rather low sample sizes. In the present study, coinfections of Ixodes ricinus with these pathogens were investigated by determining their association in a meta-analysis. A total of 5079 tick samples examined for Rickettsia and Borrelia spp. via probe-based quantitative real-time PCR in previous prevalence studies or as submitted diagnostic material were included. In Borrelia-positive ticks, genospecies were determined by Reverse Line Blot. Determination of bacterial loads resulted in an increase between developmental tick stages with highest mean bacterial loads in female ticks (7.96×10⁴ in Borrelia single-infected, 4.87×10⁵ in Rickettsia single-infected and 3.22×10⁵ in Borrelia-Rickettsia coinfected females). The determined Borrelia-Rickettsia tick coinfection rate was 12.3% (626/5079) with a significant difference to the expected coinfection rate of 9.0% (457/5079). A significant slight association as well as correlation between Borrelia and Rickettsia were determined. In addition, a significant interrelation of the bacterial load in coinfected ticks was shown. At the level of Borrelia genospecies, significant weak associations with Rickettsia spp. were detected for B. afzelii, B. garinii/bavariensis, B. valaisiana and B. lusitaniae. The positive association provides evidence for interactions between Borrelia and Rickettsia spp. in the tick vector, presumably resulting in higher bacterial replication rates in the tick vector and possibly the reservoir host. However, coinfection may impact the vector negatively as indicated by an absent increase in coinfection rates from nymphs to adults. Future studies are needed to investigate the underlying mechanisms of the positive association in ticks and possible associations in the vertebrate host as well as the potential influence of environmental factors.

Rapid Gene Turnover as a Significant Source of Genetic Variation in a Recently Seeded Population of a Healthcare-Associated Pathogen

Genome sequencing has been useful to gain an understanding of bacterial evolution. It has been used for studying the phylogeography and/or the impact of mutation and recombination on bacterial populations. However, it has rarely been used to study gene turnover at microevolutionary scales. Here, we sequenced Mexican strains of the human pathogen Acinetobacter baumannii sampled from the same locale over a 3 year period to obtain insights into the microevolutionary dynamics of gene content variability. We found that the Mexican A. baumannii population was recently founded and has been emerging due to a rapid clonal expansion. Furthermore, we noticed that on average the Mexican strains differed from each other by over 300 genes and, notably, this gene content variation has accrued more frequently and faster than the accumulation of mutations. Moreover, due to its rapid pace, gene content variation reflects the phylogeny only at very short periods of time. Additionally, we found that the external branches of the phylogeny had almost 100 more genes than the internal branches. All in all, these results show that rapid gene turnover has been of paramount importance in producing genetic variation within this population and demonstrate the utility of genome sequencing to study alternative forms of genetic variation.

Completed Genome Sequences of Borrelia burgdorferi Sensu Stricto B31(NRZ) and Closely Related Patient Isolates from Europe

Borrelia burgdorferi sensu stricto is a causative agent of human Lyme borreliosis in the United States and Europe. We report here the completed genome sequences of strain B31 isolated from a tick in the United States and two closely related strains from Europe, PAli and PAbe, which were isolated from patients with erythema migrans and neuroborreliosis, respectively.

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

Background

Borrelia (B.) burgdorferi sensu lato, including the tick-transmitted agents of human Lyme borreliosis, have particularly complex genomes, consisting of a linear main chromosome and numerous linear and circular plasmids. The number and structure of plasmids is variable even in strains within a single genospecies. Genes on these plasmids are known to play essential roles in virulence and pathogenicity as well as host and vector associations. For this reason, it is essential to explore methods for rapid and reliable characterisation of molecular level changes on plasmids.

In this study we used three strains: a low passage isolate of B. burgdorferi sensu stricto strain B31(−NRZ) and two closely related strains (PAli and PAbe) that were isolated from human patients. Sequences of these strains were compared to the previously sequenced reference strain B31 (available in GenBank) to obtain proof-of-principle information on the suitability of next generation sequencing (NGS) library construction and sequencing methods on the assembly of bacterial plasmids. We tested the effectiveness of different short read assemblers on Illumina sequences, and of long read generation methods on sequence data from Pacific Bioscience single-molecule real-time (SMRT) and nanopore (Oxford Nanopore Technologies) sequencing technology.

Results

Inclusion of mate pair library reads improved the assembly in some plasmids as did prior enrichment of plasmids. While cp32 plasmids remained refractory to assembly using only short reads they were effectively assembled by long read sequencing methods. The long read SMRT and nanopore sequences came, however, at the cost of indels (insertions or deletions) appearing in an unpredictable manner. Using long and short read technologies together allowed us to show that the three B. burgdorferi s.s. strains investigated here, whilst having similar plasmid structures to each other (apart from fusion of cp32 plasmids), differed significantly from the reference strain B31-GB, especially in the case of cp32 plasmids.

Conclusion

Short read methods are sufficient to assemble the main chromosome and many of the plasmids in B. burgdorferi. However, a combination of short and long read sequencing methods is essential for proper assembly of all plasmids including cp32 and thus, for gaining an understanding of host- or vector adaptations. An important conclusion from our work is that the evolution of Borrelia plasmids appears to be dynamic. This has important implications for the development of useful research strategies to monitor the risk of Lyme disease occurrence and how to medically manage it.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3804-5) contains supplementary material, which is available to authorized users.

Nested coevolutionary networks shape the ecological relationships of ticks, hosts, and the Lyme disease bacteria of the Borrelia burgdorferi (s.l.) complex

Background

The bacteria of the Borrelia burgdorferi (s.l.) (BBG) complex constitute a group of tick-transmitted pathogens that are linked to many vertebrate and tick species. The ecological relationships between the pathogens, the ticks and the vertebrate carriers have not been analysed. The aim of this study was to quantitatively analyse these interactions by creating a network based on a large dataset of associations. Specifically, we examined the relative positions of partners in the network, the phylogenetic diversity of the tick’s hosts and its impact on BBG circulation. The secondary aim was to evaluate the segregation of BBG strains in different vectors and reservoirs.

Results

BBG circulates through a nested recursive network of ticks and vertebrates that delineate closed clusters. Each cluster contains generalist ticks with high values of centrality as well as specialist ticks that originate nested sub-networks and that link secondary vertebrates to the cluster. These results highlighted the importance of host phylogenetic diversity for ticks in the circulation of BBG, as this diversity was correlated with high centrality values for the ticks. The ticks and BBG species in each cluster were not significantly associated with specific branches of the phylogeny of host genera (R² = 0.156, P = 0.784 for BBG; R² = 0.299, P = 0.699 for ticks). A few host genera had higher centrality values and thus higher importance for BBG circulation. However, the combined contribution of hosts with low centrality values could maintain active BBG foci. The results suggested that ticks do not share strains of BBG, which were highly segregated among sympatric species of ticks.

Conclusions

We conclude that BBG circulation is supported by a highly redundant network. This network includes ticks with high centrality values and high host phylogenetic diversity as well as ticks with low centrality values. This promotes ecological sub-networks and reflects the high resilience of BBG circulation. The functional redundancy in BBG circulation reduces disturbances due to the removal of vertebrates as it allows ticks to fill other biotic niches.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-016-1803-z) contains supplementary material, which is available to authorized users.