Plasmid diversity and phylogenetic consistency in the Lyme disease agent Borrelia burgdorferi

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.

Genotyping and Quantifying Lyme Pathogen Strains by Deep Sequencing of the Outer Surface Protein C (ospC) Locus

A mixed infection of a single tick or host by Lyme disease spirochetes is common and a unique challenge for the diagnosis, treatment, and surveillance of Lyme disease. Here, we describe a novel protocol for differentiating Lyme strains on the basis of deep sequencing of the hypervariable outer surface protein C locus (ospC). Improving upon the traditional DNA-DNA hybridization method, the next-generation sequencing-based protocol is high throughput, quantitative, and able to detect new pathogen strains. We applied the method to more than one hundred infected Ixodes scapularis ticks collected from New York State, USA, in 2015 and 2016. An analysis of strain distributions within individual ticks suggests an overabundance of multiple infections by five or more strains, inhibitory interactions among coinfecting strains, and the presence of a new strain closely related to Borreliella bissettiae A supporting bioinformatics pipeline has been developed. The newly designed pair of universal ospC primers target intergenic sequences conserved among all known Lyme pathogens. The protocol could be used for culture-free identification and quantification of Lyme pathogens in wildlife and potentially in clinical specimens.

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.

RNase III Processing of rRNA in the Lyme Disease Spirochete Borrelia burgdorferi

The rRNA genes of Borrelia (Borreliella) burgdorferi are unusually organized; the spirochete has a single 16S rRNA gene that is more than 3 kb from a tandem pair of 23S-5S rRNA operons. We generated an rnc null mutant in B. burgdorferi that exhibits a pleiotropic phenotype, including decreased growth rate and increased cell length. Here, we demonstrate that endoribonuclease III (RNase III) is, as expected, involved in processing the 23S rRNA in B. burgdorferi The 5' and 3' ends of the three rRNAs were determined in the wild type and rncBb mutants; the results suggest that RNase III in B. burgdorferi is required for the full maturation of the 23S rRNA but not for the 5S rRNA nor, curiously, for the 16S rRNA.IMPORTANCE Lyme disease, the most common tick-borne zoonosis in the Northern Hemisphere, is caused by the bacterium Borrelia (Borreliella) burgdorferi, a member of the deeply branching spirochete phylum. B. burgdorferi carries a limited suite of ribonucleases, enzymes that cleave RNA during processing and degradation. Several ribonucleases, including RNase III, are involved in the production of ribosomes, which catalyze translation and are a major target of antibiotics. This is the first study to dissect the role of an RNase in any spirochete. We demonstrate that an RNase III mutant is viable but has altered processing of rRNA.

Genomic and phenotypic characterization of Borrelia afzelii BO23 and Borrelia garinii CIP 103362

In recent years, the number of Lyme disease or borreliosis cases in Eurasia has been dramatically increasing. This tick-borne disease is caused by Borrelia burgdorferi sensu lato, which includes B. burgdorferi sensu stricto, the main species found in North America, and B. afzelii and B. garinii, which are primarily responsible for the disease in Eurasia. Currently, research on Lyme disease has focused mainly on B. burgdorferi while B. afzelii and B. garinii, which cause disease with distinctly different symptoms, are less studied. The purpose of this study is to evaluate B. afzelii BO23 and B. garinii CIP 103362 as model organisms to study Eurasian Lyme disease. To begin our analyses, we sequenced, annotated the chromosomes of both species and compared them to B. burgdorferi strain B31. We also assayed shuttle vector, pBSV2, for transformation efficacy and demonstrated that these strains can be cultured on solid media. In addition, we characterized how physicochemical parameters (e.g., oxygen, osmolarity, oxidative stress) affect both growth and motility of the bacteria. Finally, we describe each strain's antibiotic susceptibility and accessed their ability to infect mice. In conclusion, B. afzelii BO23 was more practical for in vitro and in vivo studies than B. garinii CIP 103362.

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.

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.

BpWrapper: BioPerl-based sequence and tree utilities for rapid prototyping of bioinformatics pipelines

BACKGROUND

Automated bioinformatics workflows are more robust, easier to maintain, and results more reproducible when built with command-line utilities than with custom-coded scripts. Command-line utilities further benefit by relieving bioinformatics developers to learn the use of, or to interact directly with, biological software libraries. There is however a lack of command-line utilities that leverage popular Open Source biological software toolkits such as BioPerl ( http://bioperl.org ) to make many of the well-designed, robust, and routinely used biological classes available for a wider base of end users.

RESULTS

Designed as standard utilities for UNIX-family operating systems, BpWrapper makes functionality of some of the most popular BioPerl modules readily accessible on the command line to novice as well as to experienced bioinformatics practitioners. The initial release of BpWrapper includes four utilities with concise command-line user interfaces, bioseq, bioaln, biotree, and biopop, specialized for manipulation of molecular sequences, sequence alignments, phylogenetic trees, and DNA polymorphisms, respectively. Over a hundred methods are currently available as command-line options and new methods are easily incorporated. Performance of BpWrapper utilities lags that of precompiled utilities while equivalent to that of other utilities based on BioPerl. BpWrapper has been tested on BioPerl Release 1.6, Perl versions 5.10.1 to 5.25.10, and operating systems including Apple macOS, Microsoft Windows, and GNU/Linux. Release code is available from the Comprehensive Perl Archive Network (CPAN) at https://metacpan.org/pod/Bio::BPWrapper . Source code is available on GitHub at https://github.com/bioperl/p5-bpwrapper .

CONCLUSIONS

BpWrapper improves on existing sequence utilities by following the design principles of Unix text utilities such including a concise user interface, extensive command-line options, and standard input/output for serialized operations. Further, dozens of novel methods for manipulation of sequences, alignments, and phylogenetic trees, unavailable in existing utilities (e.g., EMBOSS, Newick Utilities, and FAST), are provided. Bioinformaticians should find BpWrapper useful for rapid prototyping of workflows on the command-line without creating custom scripts for comparative genomics and other bioinformatics applications.

Whole-Genome Sequencing of Six Borrelia miyamotoi Clinical Strains Isolated in Russia

Here, we report the whole-genome sequence of six clinical Borrelia miyamotoi isolates from the Russian Federation. Using two independent next-generation sequencing platforms, we determined the complete sequence of the chromosome and several plasmids. All strains have an Asian genotype with 99.8% chromosome nucleotide similarity with B. miyamotoi strain FR64b.

Genomic insights into the ancient spread of Lyme disease across North America

Lyme disease is the most prevalent vector-borne disease in North America and continues to spread. The disease was first clinically described in the 1970s in Lyme, Connecticut, but the origins and history of spread of the Lyme disease bacteria, Borrelia burgdorferi sensu stricto (s.s.), are unknown. To explore the evolutionary history of B. burgdorferi in North America, we collected ticks from across the USA and southern Canada from 1984 to 2013 and sequenced the, to our knowledge, largest collection of 146 B. burgdorferi s.s. genomes. Here, we show that B. burgdorferi s.s. has a complex evolutionary history with previously undocumented levels of migration. Diversity is ancient and geographically widespread, well pre-dating the Lyme disease epidemic of the past ~40 years, as well as the Last Glacial Maximum ~20,000 years ago. This means the recent emergence of human Lyme disease probably reflects ecological change-climate change and land use changes over the past century-rather than evolutionary change of the bacterium.