MLST of housekeeping genes captures geographic population structure and suggests a European origin of Borrelia burgdorferi

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.

Multilocus sequence typing of pathogenic treponemes isolated from contagious ovine digital dermatitis stage five lesions: Implications for disease transmission dynamics

Contagious ovine digital dermatitis (CODD) causes a severe, infectious foot disease and lameness of sheep, is common within the UK and is now also emerging in other countries. As well as causing severe animal welfare issues, huge economic losses emerge from the disease due to weight loss/lack of weight gain, and veterinary treatments. CODD lesion progress is measured, with a scoring system from 1 (early lesions) to 5 (healed). Here, using samples from an experimental flock infected by natural means, samples were taken from CODD stage 5 lesions, post treatment, and subjected to bacterial isolation and MLST using previously published methods. Sequences were compared to others from the same flock, and those from previous studies. All CODD 5 lesions produced viable Treponema spp. bacteria. High levels of variation of bacteria were seen, with 12 sequence types (STs) for T. medium phylogroup (11 new), 15 STs for T. phagedenis phylogroup (9 new) and six T. pedis STs, of which two were new. This study shows that CODD stage 5 lesions still contain viable bacteria, representing all three known pathogenic Treponema spp. phylogroups, and these may thus play a role in disease transmission and epidemiology despite appearing healed after treatment. The high level CODD treponeme variability within an infected flock where sheep were bought from different sources, as might occur in common agricultural practice, may suggest reasons as to why the bacterial disease is difficult to treat, control and eradicate, and adds further complexity to the polybacterial pathogenesis of these lesions.

Isolation and molecular characterization of a novel relapsing fever group Borrelia from the white-eared opossum Didelphis albiventris in Brazil

This study aimed to detect, isolate and to characterize by molecular methods a relapsing fever group (RFG) Borrelia in white-eared opossums (Didelphis albiventris) from Brazil. During 2015-2018, when opossums (Didelphis spp.) were captured in six municipalities of the state of São Paulo, Brazil, molecular analyses revealed the presence of a novel RFG Borrelia sp. in the blood of seven opossums (Didelphis albiventris), out of 142 sampled opossums (4.9% infection rate). All seven infected opossums were from a single location (Ribeirão Preto municipality). In a subsequent field study in Ribeirão Preto during 2021, two new opossums (D. albiventris) were captured, of which one contained borrelial DNA in its blood. Macerated tissues from this infected opossum were inoculated into laboratory animals (rodents and rabbits) and two big-eared opossums (Didelphis aurita), which had blood samples examined daily via dark-field microscopy. No spirochetes were visualized in the blood of the laboratory animals. Contrastingly, spirochetes were visualized in the blood of the two D. aurita opossums between 12 and 25 days after inoculation. Blood samples from these opossums were used for a multi-locus sequencing typing (MLST) based on six borrelial loci. Phylogenies inferred from MLST genes positioned the sequenced Borrelia genotype into the RFG borreliae clade basally to borreliae of the Asian-African group, forming a monophyletic group with another Brazilian isolate, "Candidatus B. caatinga". Based on this concatenated phylogenetic analysis, which supports that the new borrelial isolate corresponds to a putative new species, we propose the name "Candidatus Borrelia mimona".

Potential drivers of vector-borne pathogens in urban environments: European hedgehogs (Erinaceus europaeus) in the spotlight

Vector-borne diseases (VBDs) are considered as (re-)emerging, but information on the transmission cycles and wildlife reservoirs is often incomplete, particularly with regard to urban areas. The present study investigated blood samples from European hedgehogs (Erinaceus europaeus) presented at wildlife rehabilitation centres in the region of Hanover. Past exposure to B. burgdorferi sensu lato (s.l.) and tick-borne encephalitis virus (TBEV) was assessed by serological detection of antibodies, while current infections with Borrelia spp., Anaplasma phagocytophilum, Rickettsia spp., Neoehrlichia mikurensis, Bartonella spp., Babesia spp. and Spiroplasma ixodetis were investigated by (q)PCR. Of 539 hedgehogs tested for anti-Borrelia antibodies, 84.8% (457/539) were seropositive, with a higher seropositivity rate in adult than subadult animals, while anti-TBEV antibodies were detected in one animal only (0.2%; 1/526). By qPCR, 31.2% (168/539) of hedgehog blood samples were positive for Borrelia spp., 49.7% (261/525) for A. phagocytophilum, 13.0% (68/525) for Bartonella spp., 8.2% for S. ixodetis (43/525), 8.0% (42/525) for Rickettsia spp. and 1.3% (7/525) for Babesia spp., while N. mikurensis was not detected. While further differentiation of Borrelia spp. infections was not successful, 63.2% of the A. phagocytophilum infections were assigned to the zoonotic ecotype I and among Rickettsia spp. infections, 50.0% to R. helvetica by ecotype- or species-specific qPCR, respectively. Sequencing revealed the presence of a Rickettsia sp. closely related to Rickettsia felis in addition to a Bartonella sp. previously described from hedgehogs, as well as Babesia microti and Babesia venatorum. These findings show that hedgehogs from rehabilitation centres are valuable sources to identify One Health pathogens in urban areas. The hedgehogs are not only exposed to pathogens from fleas and ticks in urban areas, but they also act as potent amplifiers for these vectors and their pathogens, relevant for citizens and their pets.

BS, Raghuraman V, Rajyaguru U, Llamera KE, Andrew L, Anderson AS, Hovius JW, Liberator PA, Simon R, Hao L. Development of a sequence-based in silico OspA typing method for Borrelia burgdorferi sensu lato

Lyme disease (LD), caused by spirochete bacteria of the genus Borrelia burgdorferi sensu lato, remains the most common vector-borne disease in the northern hemisphere. Borrelia outer surface protein A (OspA) is an integral surface protein expressed during the tick cycle, and a validated vaccine target. There are at least 20 recognized Borrelia genospecies, that vary in OspA serotype. This study presents a new in silico sequence-based method for OspA typing using next-generation sequence data. Using a compiled database of over 400 Borrelia genomes encompassing the 4 most common disease-causing genospecies, we characterized OspA diversity in a manner that can accommodate existing and new OspA types and then defined boundaries for classification and assignment of OspA types based on the sequence similarity. To accommodate potential novel OspA types, we have developed a new nomenclature: OspA in silico type (IST). Beyond the ISTs that corresponded to existing OspA serotypes 1-8, we identified nine additional ISTs that cover new OspA variants in B. bavariensis (IST9-10), B. garinii (IST11-12), and other Borrelia genospecies (IST13-17). The IST typing scheme and associated OspA variants are available as part of the PubMLST Borrelia spp. database. Compared to traditional OspA serotyping methods, this new computational pipeline provides a more comprehensive and broadly applicable approach for characterization of OspA type and Borrelia genospecies to support vaccine development.

Borrelia puertoricensis in opossums (Didelphis marsupialis) from Colombia

BACKGROUND

The genus Borrelia comprises pathogenic species of bacteria that pose a significant risk to public health. Borrelia spp. are emerging or reemerging infectious agents worldwide with complex transmission cycles, and many species use rodents as vertebrate reservoir hosts. Spirochetes morphologically compatible with Borrelia have been recurrently observed in opossums; however, there is currently a lack of genetic evidence confirming infection or supporting that these marsupials are hosts of Borrelia spirochetes.

METHODS

During 2017, 53 serum samples of Didelphis marsupialis from the municipality of Colosó (department of Sucre, Colombia) were collected and allocated in a serum bank. DNA extracted from the serum samples was submitted to a Borrelia genus-specific real-time PCR targeting the 16S rRNA gene. Positive samples were subsequently derived from semi-nested PCR protocols to obtain large fragments of the 16S rRNA and flaB genes. Obtained amplicons were subjected to Sanger sequencing. One positive sample was randomly selected for next-generation sequencing (NGS). Obtained reads were mapped to genomes of Borrelia spp. and sequences of two genes used in a multilocus sequence typing scheme retrieved for taxonomic assignment and phylogenetic analyses.

RESULTS

Overall, 18.8% (10/53) of the samples were positive by qPCR. Of them, 80% (8/10) and 60% (6/10) were positive for the 16S rRNA and flaB genes after semi-nested PCRs, respectively. Bioinformatic analysis of one sample sequenced with NGS yielded 22 reads of genus Borrelia with different sizes. Two housekeeping genes, rplB and pyrG, were recovered. Nucleotide pairwise comparisons and phylogenetic analyses of 16S rRNA, flaB, rplB and pyrG genes showed that the Borrelia sp. found in opossums from Colosó corresponded to Borrelia puertoricensis.

CONCLUSIONS

We describe the first molecular evidence to our knowledge of B. puertoricensis in Colombia, specifically in opossums, and the first detection of this spirochete in a vertebrate host since its isolation from Ornithodoros puertoricensis in Panama. This detection is also relevant because of the epidemiological importance of opossums as reservoirs of zoonotic diseases to humans.

Molecular Typing on Human Blood Reveals the Borrelia afzelii Infection in Korea

Lyme disease is a tick-borne infection in Korea. Here, clinical samples were collected from a 72-year old patient, with sudden onset of fever on April, 2018. The patient was passed away after 3rd day of doxycycline administration. The molecular diagnostic tests, nested polymerase chain reaction targeting 5S-23S rRNA intergenic spacer region (IGS) and multilocus sequence typing (MLST), showed positive for Borrelia afzelii from blood. Further, mutations in both 5S - 23S IGS and pepX allele of MLST were determined. Herein, we report the expected first death case by B. afzelii infection in Korea.

Concurrent Infection of the Human Brain with Multiple Borrelia Species

Lyme disease (LD) spirochetes are well known to be able to disseminate into the tissues of infected hosts, including humans. The diverse strategies used by spirochetes to avoid the host immune system and persist in the host include active immune suppression, induction of immune tolerance, phase and antigenic variation, intracellular seclusion, changing of morphological and physiological state in varying environments, formation of biofilms and persistent forms, and, importantly, incursion into immune-privileged sites such as the brain. Invasion of immune-privileged sites allows the spirochetes to not only escape from the host immune system but can also reduce the efficacy of antibiotic therapy. Here we present a case of the detection of spirochetal DNA in multiple loci in a LD patient's post-mortem brain. The presence of co-infection with Borrelia burgdorferi sensu stricto and Borrelia garinii in this LD patient's brain was confirmed by PCR. Even though both spirochete species were simultaneously present in human brain tissue, the brain regions where the two species were detected were different and non-overlapping. The presence of atypical spirochete morphology was noted by immunohistochemistry of the brain samples. Atypical morphology was also found in the tissues of experimentally infected mice, which were used as a control.

Possible effect of mutations on serological detection of Borrelia burgdorferi sensu stricto ospC major groups: An in-silico study

The outer surface protein C (OspC) of the agent of Lyme disease, Borrelia burgdorferi sensu stricto, is a major lipoprotein surface-expressed during early-phase human infections. Antibodies to OspC are used in serological diagnoses. This study explored the hypothesis that serological test sensitivity decreases as genetic similarity of ospC major groups (MGs) of infecting strains, and ospC A (the MG in the strain B31 used to prepare antigen for serodiagnosis assays) decreases. We used a previously published microarray dataset to compare serological reactivity to ospC A (measured as pixel intensity) versus reactivity to 22 other ospC MGs, within a population of 55 patients diagnosed by two-tier serological testing using B. burgdorferi s.s. strain B31 as antigen, in which the ospC MG is OspC A. The difference in reactivity of sera to ospC A and reactivity to each of the other 22 ospC MGs (termed 'reactivity difference') was the outcome variable in regression analysis in which genetic distance of the ospC MGs from ospC A was the explanatory variable. Genetic distance was computed for the whole ospC sequence, and 9 subsections, from Neighbour Joining phylogenetic trees of the 23 ospC MGs. Regression analysis was conducted using genetic distance for the full ospC sequence, and the subsections individually. There was a significant association between the reactivity difference and genetic distance of ospC MGs from ospC A: increased genetic distance reduced reactivity to OspC A. No single ospC subsection sequence fully explained the relationship between genetic distance and reactivity difference. An analysis of single nucleotide polymorphisms supported a biological explanation via specific amino acid modifications likely to change protein binding affinity. This adds support to the hypothesis that genetic diversity of B. burgdorferi s.s. (here specifically OspC) may impact serological diagnostic test performance. Further prospective studies are necessary to explore the clinical implications of these findings.

Seroprevalence and risk factors of Borrelia burgdorferi sensu lato and Rickettsia species infection in humans in Mongolia, 2016-2020

Borrelia burgdorferi sensu lato and Rickettsia spp. are worldwide causes of tick-borne infections. We aimed to estimate the seroprevalence of immunoglobulin G (IgG) antibodies against different tick-borne diseases (TBDs) and determine risk factors among Mongolians from 2016 to 2020. Blood samples were obtained from voluntary participants with a history of suspected tick bite who visited our hospital, and IgG antibodies against Rickettsia and Borrelia were detected using enzyme-linked immunosorbent assay (ELISA). The IgG antibody seropositivity rate against Rickettsia was 21.8% (1032/4724), while 3.4% (162/4724) of participants tested positive for serum IgG antibodies against Borrelia by ELISA.Binary logistic regression analysis was performed to evaluate risk factors for tick-borne rickettsiosis (TBR) and tick-borne borreliosis (TBB) using IgG serum sample. Age, occupation, and residence were significantly associated with these diseases; however, sex did not show any significant association. Seroprevalence was significantly higher among herders (40.6%, 95% confidence interval [CI]: 35.5-45.8; odds ratio [OR] 0.61; P < 0.001) and students (32.8%, 95% CI: 30.2-35.4; OR 0.75; P < 0.001) than among individuals with other occupations. The 25-29 age group had a slightly higher seroprevalence (35.1%, 95% CI: 28.1-42.6; OR 0.61; P < 0.006) than those in other age groups. Province was a stronger predictor of TBR than occupation and age group. In univariate subgroup analysis by age group, occupation, and residence were significantly associated with TBR seroprevalence, whereas age and province were associated with TBB seroprevalence. Thus, risk factors for TBD include residence, occupation, and age group. This study was conducted using samples from all Mongolian provinces and the capital city, and the risk factors and prevalence of Rickettsia and Borreliaare highlighted.

Lyme borreliosis diagnosis: state of the art of improvements and innovations

With almost 700 000 estimated cases each year in the United States and Europe, Lyme borreliosis (LB), also called Lyme disease, is the most common tick-borne illness in the world. Transmitted by ticks of the genus Ixodes and caused by bacteria Borrelia burgdorferi sensu lato, LB occurs with various symptoms, such as erythema migrans, which is characteristic, whereas others involve blurred clinical features such as fatigue, headaches, arthralgia, and myalgia. The diagnosis of Lyme borreliosis, based on a standard two-tiered serology, is the subject of many debates and controversies, since it relies on an indirect approach which suffers from a low sensitivity depending on the stage of the disease. Above all, early detection of the disease raises some issues. Inappropriate diagnosis of Lyme borreliosis leads to therapeutic wandering, inducing potential chronic infection with a strong antibody response that fails to clear the infection. Early and proper detection of Lyme disease is essential to propose an adequate treatment to patients and avoid the persistence of the pathogen. This review presents the available tests, with an emphasis on the improvements of the current diagnosis, the innovative methods and ideas which, ultimately, will allow more precise detection of LB.

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.

Structural evolution of an immune evasion determinant shapes pathogen host tropism

Modern infectious disease outbreaks often involve changes in host tropism, the preferential adaptation of pathogens to specific hosts. The Lyme disease-causing bacterium Borrelia burgdorferi (Bb) is an ideal model to investigate the molecular mechanisms of host tropism, because different variants of these tick-transmitted bacteria are distinctly maintained in rodents or bird reservoir hosts. To survive in hosts and escape complement-mediated immune clearance, Bb produces the outer surface protein CspZ that binds the complement inhibitor factor H (FH) to facilitate bacterial dissemination in vertebrates. Despite high sequence conservation, CspZ variants differ in human FH-binding ability. Together with the FH polymorphisms between vertebrate hosts, these findings suggest that minor sequence variation in this bacterial outer surface protein may confer dramatic differences in host-specific, FH-binding-mediated infectivity. We tested this hypothesis by determining the crystal structure of the CspZ-human FH complex, and identifying minor variation localized in the FH-binding interface yielding bird and rodent FH-specific binding activity that impacts infectivity. Swapping the divergent region in the FH-binding interface between rodent- and bird-associated CspZ variants alters the ability to promote rodent- and bird-specific early-onset dissemination. We further linked these loops and respective host-specific, complement-dependent phenotypes with distinct CspZ phylogenetic lineages, elucidating evolutionary mechanisms driving host tropism emergence. Our multidisciplinary work provides a novel molecular basis for how a single, short protein motif could greatly modulate pathogen host tropism.

Borrelia miyamotoi Infection in Immunocompromised Man, California, USA, 2021

Infection with Borrelia miyamotoi in California, USA, has been suggested by serologic studies. We diagnosed B. miyamotoi infection in an immunocompromised man in California. Diagnosis was aided by plasma microbial cell-free DNA sequencing. We conclude that the infection was acquired in California.

Borrelia burgdorferi Outer Surface Protein C Is Not the Sole Determinant of Dissemination in Mammals

Lyme disease in the United States is most often caused by Borrelia burgdorferi sensu stricto. After a tick bite, the patient may develop erythema migrans at that site. If hematogenous dissemination occurs, the patient may then develop neurologic manifestations, carditis, or arthritis. Host-pathogen interactions include factors that contribute to hematogenous dissemination to other body sites. Outer surface protein C (OspC), a surface-exposed lipoprotein of B. burgdorferi, is essential during the early stages of mammalian infection. There is a high degree of genetic variation at the ospC locus, and certain ospC types are more frequently associated with hematogenous dissemination in patients, suggesting that OspC may be a major contributing factor to the clinical outcome of B. burgdorferi infection. In order to evaluate the role of OspC in B. burgdorferi dissemination, ospC was exchanged between B. burgdorferi isolates with different capacities to disseminate in laboratory mice, and these strains were then tested for their ability to disseminate in mice. The results indicated that the ability of B. burgdorferi to disseminate in mammalian hosts does not depend on OspC alone. The complete genome sequences of two closely related strains of B. burgdorferi with differing dissemination phenotypes were determined, but a specific genetic locus that could explain the differences in the phenotypes could not be definitively identified. The animal studies performed clearly demonstrated that OspC is not the sole determinant of dissemination. Future studies of the type described here with additional borrelial strains will hopefully clarify the genetic elements associated with hematogenous dissemination.

Isolation and multilocus sequence typing of Borrelia burgdorferi from Ixodes scapularis collected from dogs in Ontario, Canada

OBJECTIVE

To identify the multilocus sequence typing (MLST) sequence types of Borrelia burgdorferi from Ixodes scapularis in Ontario, Canada.

RESULTS

One hundred and eighty-five I. scapularis ticks were submitted from 134 dogs via participating clinics from April 1, 2019, to March 31, 2020. Seventeen MLST sequence types of B. burgdorferi were detected from fifty-eight cultured isolates from 21 ticks. The most common MLST sequence types were 12 and 16. Mixed infections of two MLST sequence types were detected in four ticks. Three sequence types (48, 317, 639) were new detections in Ontario.

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.

A Novel Relapsing Fever Group Borrelia Isolated from Ornithodoros Ticks of the Brazilian Caatinga

Tick-borne relapsing fever group (RFG) borreliosis remains neglected as a human disease and little is known on its maintenance in ticks and vertebrates, especially in South America. Therefore, this study investigated borrelial infection in Ornithodoros ticks collected in rodent-inhabited rock formations in the Brazilian semiarid region, within the Caatinga biome. Collected ticks (Ornithodoros rietcorreai and Ornithodoros cf. tabajara) were allowed to feed under laboratory conditions on guinea pigs, which had blood samples examined daily by dark-field microscopy. No spirochetes were visualized in the blood of any of four O. rietcorreai-infested guinea pigs. Contrastingly, spirochetes were visualized between 9 and 39 days after tick feeding in the blood of three guinea pigs, each infested with O. cf. tabajara ticks from a different locality. Guinea pig infection was confirmed by passages into experimental animals and by generating DNA sequences of Borrelia spp. from the blood of spirochetemic guinea pigs. Three O. cf. tabajara populations were infected by the same borrelial organism, which was characterized as a novel RFG agent (named as 'Candidatus Borrelia caatinga') based on 10 Borrelia loci (rrs, flaB, glpQ, gyrB, clpX, pepX, pyrG, recG, rplB and uvrA). We demonstrated that O. cf. tabajara is a competent vector of the novel Borrelia sp. isolates, although none of the infected rodents developed clinical illness.

Out of Asia? Expansion of Eurasian Lyme borreliosis causing genospecies display unique evolutionary trajectories

Vector-borne pathogens exist in obligate transmission cycles between vector and reservoir host species. Host and vector shifts can lead to geographic expansion of infectious agents and the emergence of new diseases in susceptible individuals. Three bacterial genospecies (Borrelia afzelii, Borrelia bavariensis, and Borrelia garinii) predominantly utilize two distinct tick species as vectors in Asia (Ixodes persulcatus) and Europe (Ixodes ricinus). Through these vectors, the bacteria can infect various vertebrate groups (e.g., rodents, birds) including humans where they cause Lyme borreliosis, the most common vector-borne disease in the Northern hemisphere. Yet, how and in which order the three Borrelia genospecies colonized each continent remains unclear including the evolutionary consequences of this geographic expansion. Here, by reconstructing the evolutionary history of 142 Eurasian isolates, we found evidence that the ancestors of each of the three genospecies probably have an Asian origin. Even so, each genospecies studied displayed a unique substructuring and evolutionary response to the colonization of Europe. The pattern of allele sharing between continents is consistent with the dispersal rate of the respective vertebrate hosts, supporting the concept that adaptation of Borrelia genospecies to the host is important for pathogen dispersal. Our results highlight that Eurasian Lyme borreliosis agents are all capable of geographic expansion with host association influencing their dispersal; further displaying the importance of host and vector association to the geographic expansion of vector-borne pathogens and potentially conditioning their capacity as emergent pathogens.

Molecular Surveillance for Vector-Borne Bacteria in Rodents and Tree Shrews of Peninsular Malaysia Oil Palm Plantations

Many human clinical cases attributed to vector-borne pathogens are underreported in Malaysia, especially in rural localities where healthcare infrastructures are lacking. Here, 217 small mammals, consisting of rodents and tree shrews, were trapped in oil palm plantations in the Peninsular Malaysia states of Johor and Perak. Species identification was performed using morphological and DNA barcoding analyses, and 203 small mammals were included in the detection of selected vector-borne bacteria. The DNA extracted from the spleens was examined for Orientia tsutsugamushi, Borrelia spp., Bartonella spp. and Rickettsia spp. using established PCR assays. The small mammals collected in this study included Rattus tanezumi R3 mitotype (n = 113), Rattus argentiventer (n = 24), Rattus tiomanicus (n = 22), Rattus exulans (n = 17), Rattus tanezumi sensu stricto (n = 1) and Tupaia glis (n = 40). Orientia tsutsugamushi, Borrelia spp. and Bartonella phoceensis were detected in the small mammals with the respective detection rates of 12.3%, 5.9% and 4.9%. Rickettsia spp., however, was not detected. This study encountered the presence of both Lyme disease and relapsing fever-related borreliae in small mammals collected from the oil palm plantation study sites. All three microorganisms (Orientia tsutsugamushi, Borrelia spp. and Bartonella phoceensis) were detected in the R. tanezumi R3 mitotype, suggesting that the species is a competent host for multiple microorganisms. Further investigations are warranted to elucidate the relationships between the ectoparasites, the small mammals and the respective pathogens.

Detection and genetic characterization of a putative novel Borrelia genospecies in Ixodes apronophorus / Ixodes persulcatus / Ixodes trianguliceps sympatric areas in Western Siberia

Four genospecies from the Borrelia burgdorferi sensu lato complex were detected in Ixodes persulcatus and Ixodes pavlovskyi ticks from Siberia and genetically characterized. The presence of Borrelia spp. in Ixodes apronophorus and Ixodes trianguliceps ticks found in Asia has never been studied. In this study, genetic diversity of B. burgdorferi s.l. was investigated in three I. persulcatus / I. trianguliceps / I. apronophorus sympatric habitats in Western Siberia. Three groups of samples were examined: (i) ticks that were taken from rodents and molted in a laboratory; (ii) non-molted ticks collected from rodents; (iii) specimens from small mammals. Expectedly, Borrelia afzelii and Borrelia bavariensis were detected in I. persulcatus and in small mammals from the studied locations. Borrelia bavariensis was first found in molted I. apronophorus and I. trianguliceps. Identical genovariants of B. bavariensis were found in I. apronophorus, I. trianguliceps, and I. persulcatus. In addition, a new Borrelia genovariant was discovered in non-molted and molted I. apronophorus and non-molted I. persulcatus and I. trianguliceps, as well as in small mammals. This new genovariant was genetically characterized using MLST and single locus sequence analysis, which indicated that the new Borrelia genovariant significantly differs from all known Borrelia species. We propose the name "Candidatus Borrelia sibirica" for this putative new species.

Outer surface protein E (OspE) mediates Borrelia burgdorferi sensu stricto strain-specific complement evasion in the eastern fence lizard, Sceloporus undulatus

In North America, Lyme disease is primarily caused by the spirochetal bacterium Borrelia burgdorferi sensu stricto (Bb), which is transmitted between multiple vertebrate hosts and ixodid ticks, and is a model commonly used to study host-pathogen interactions. While Bb is consistently observed in its mammalian and avian reservoirs, the bacterium is rarely isolated from North American reptiles. Two closely related lizard species, the eastern fence lizard (Sceloporus undulatus) and the western fence lizard (Sceloporus occidentalis), are examples of reptiles parasitized by Ixodes ticks. Vertebrates are known to generate complement as an innate defense mechanism, which can be activated before Bb disseminate to distal tissues. Complement from western fence lizards has proven lethal against one Bb strain, implying the role of complement in making those lizards unable to serve as hosts to Bb. However, Bb DNA is occasionally identified in distal tissues of field-collected eastern fence lizards, suggesting some Bb strains may overcome complement-mediated clearance in these lizards. These findings raise questions regarding the role of complement and its impact on Bb interactions with North American lizards. In this study, we found Bb seropositivity in a small population of wild-caught eastern fence lizards and observed Bb strain-specific survivability in lizard sera. We also found that a Bb outer surface protein, OspE, from Bb strains viable in sera, promotes lizard serum survivability and binds to a complement inhibitor, factor H, from eastern fence lizards. Our data thus identify bacterial and host determinants of eastern fence lizard complement evasion, providing insights into the role of complement influencing Bb interactions with North American lizards.

Genomic Confirmation of Borrelia garinii, United States

Lyme disease is a multisystem disorder primarily caused by Borrelia burgdorferi sensu lato. However, B. garinii, which has been identified on islands off the coast of Newfoundland and Labrador, Canada, is a cause of Lyme disease in Eurasia. We report isolation and whole-genome nucleotide sequencing of a B. garinii isolate from a cotton mouse (Peromyscus gossypinus) in South Carolina, USA. We identified a second B. garinii isolate from the same repository. Phylogenetic analysis does not associate these isolates with the previously described isolates of B. garinii from Canada.

Optimized multilocus sequence analysis for laboratory identification of pathogens of ixodid tick-borne borreliosis

Introduction. The most common etiological agents of ixodid tick-borne borreliosis (ITBB) in Russia are Borrelia garinii, B. afzelii, B. bavariensis. Multilocus sequence typing and multilocus sequence analysis (MLSA) have been used in recent studies for Borrelia species identification. The results of using the MLSA scheme for identification of pathogens causing erythemic forms of ITBB have been presented earlier. The purpose of the study was to explore the possibility of MLSA optimization for laboratory identification of ITBB pathogens. Objectives: comparative analysis of nucleotide sequences of 6 conserved genes (rrs, hbb, fla, groEL, recA, ospA) and the rrfA-rrlB intergenic spacer, which are recommended by the MLSA protocol; identification of the minimum set of genes, the concatenated sequences of which are essential for species identification of Borrelia isolates. Materials and methods. The sequences of the above loci of 23 reference isolates collected from patients with ITBB and assigned, using MLSA, to B. bavariensis were compared with the sequences of similar genes of other Borrelia species available in international databases. The UPGMA method was used to build and analyze dendrograms based on the obtained data. Results. The sequences of ospA gene loci of reference species demonstrated the greatest difference (not less than 8.5%) from the sequences of the above gene in other analyzed species of Borrelia; approximately similar species-related differences (not less than 6.7%) were demonstrated by the comparison of recA gene sequences. The sequences of the identified variants of these two genes in B. bavariensis differed from the sequences of the similar genes in the most closely related species B. garinii. The dendrogram of the concatenated nucleotide sequences of recA and ospA genes demonstrated that it was totally consistent with the results of identification of the isolates based on the MLSA protocol. Conclusion. The optimized approach to MLSA of the B. burgdorferi sensu lato group suggests that species identification should be based on the concatenated analysis of loci of only two genes (recA and ospA) out of 7 loci recommended by the MLSA protocol.

A search for piroplasmids and spirochetes in threatened pudu (Pudu puda) and associated ticks from Southern Chile unveils a novel Babesia sp. and a variant of Borrelia chilensis

Cervids are important hosts for ticks and although they are refractory to some tick-borne agents such as Borrelia, they do act as reservoirs for others such as Babesia. Babesia and Borrelia are commonly transmitted by Ixodes spp. associated with deer, and most of the knowledge on their biological cycles comes from northern latitudes of the globe. In this study, we performed genetic screenings to detect tick-borne agents in blood and Ixodes stilesi ticks collected from an insular population of threatened pudu (Pudu puda), a pygmy deer species that inhabits temperate rainforests of southern South America. Inferred by phylogenetic analyses for 18S rRNA, COI and cytb genes, our results unveiled a novel genospecies of Babesia (Babesia sp. pudui) genetically related to Babesia odocoilei, a species that infects Odocoileus virginianus deer in North America. Although blood of the deer was negative for Borrelia infection, multilocus sequencing typing performed in one I. stilesi tick revealed the occurrence of a novel genetic variant of Borrelia chilensis, differing 0.93% and 0.18% in flaB and pepX genes with the type of strain for the species, respectively. Such a genetic divergence could be the result of thousands of years of isolation because of recent glaciation events that separated pudus and their tick populations at Chiloé Island approximately 437,000 years ago. The finding of a Babesia sp. has no precedents for wild and domestic ungulates in Chile and shows a novel piroplasmid that must be considered now on in rehabilitation centres and zoos that attend pudu deer. Further research is now necessary to confirm pathogenic roles.

Borrelia burgdorferi strain and host sex influence pathogen prevalence and abundance in the tissues of a laboratory rodent host

Experimental infections with different pathogen strains give insight into pathogen life history traits. The purpose of the present study was to compare variation in tissue infection prevalence and spirochete abundance among strains of Borrelia burgdorferi in a rodent host (Mus musculus, C3H/HeJ). Male and female mice were experimentally infected via tick bite with one of 12 strains. Ear tissue biopsies were taken at days 29, 59 and 89 postinfection, and seven tissues were collected at necropsy. The presence and abundance of spirochetes in the mouse tissues were measured by quantitative polymerase chain reaction. To determine the frequencies of our strains in nature, their multilocus sequence types were matched to published data sets. For the infected mice, 56.6% of the tissues were infected with B. burgdorferi. The mean spirochete load in the mouse necropsy tissues varied 4.8-fold between the strains. The mean spirochete load in the ear tissue biopsies decreased rapidly over time for some strains. The percentage of infected tissues in male mice (65.4%) was significantly higher compared to female mice (50.5%). The mean spirochete load in the seven tissues was 1.5× higher in male mice compared to female mice; this male bias was 15.3× higher in the ventral skin. Across the 11 strains, the mean spirochete loads in the infected mouse tissues were positively correlated with the strain-specific frequencies in their tick vector populations. The study suggests that laboratory-based estimates of pathogen abundance in host tissues can predict the strain composition of this important tick-borne pathogen in nature.

Use of mammalian museum specimens to test hypotheses about the geographic expansion of Lyme disease in the southeastern United States

Lyme disease, caused primarily in North America by the bacterium Borrelia burgdorferi sensu stricto, is the most frequently reported vector-borne disease in North America and its geographic extent is increasing in all directions from foci in the northeastern and north central United States. Several southeastern states, including Virginia and North Carolina, have experienced large increases in Lyme disease incidence in the past two decades, with the biggest changes in incidence occurring in the western portion of each state. We tested the hypothesis that B. burgdorferi s.s. was present in western Virginia and North Carolina Peromyscus leucopus populations prior to the recent emergence of Lyme disease. Specifically, we examined archived P. leucopus museum specimens, sampled between 1900 and 2000, for B. burgdorferi s.s. DNA. After confirming viability of DNA extracted from ear punch biopsies from P. leucopus study skins collected between 1945 and 2000 in 19 Virginia counties and 17 North Carolina counties, we used qPCR of two species-specific loci to test for the presence of B. burgdorferi s.s. DNA. Ten mice, all collected from the Eastern Shore of Virginia in 1989, tested positive for presence of B. burgdorferi; all of the remaining 344 specimens were B. burgdorferi-negative. Our results suggest that B. burgdorferi s.s was not common in western Virginia or North Carolina prior to the emergence of Lyme disease cases in the past two decades. Rather, the emergence of Lyme disease in this region has likely been driven by the relatively recent expansion of B. burgdorferi s.s. in southward-moving ticks and reservoir hosts in the mountainous counties of these two states.

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.

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.

Universal Mitochondrial Multi-Locus Sequence Analysis (mtMLSA) to Characterise Populations of Unanticipated Plant Pest Biosecurity Detections

Simple Summary

Agricultural and environmental sustainability requires effective biosecurity responses that prevent the establishment or spread of exotic insect pests. Understanding where new detections may have come from or if recurrent detections are connected contributes to this. Suitable population genetic markers use relatively rapidly evolving gene regions which render the PCR method species-specific at best. Because resource limitations mean these are pre-emptively developed for the highest risk species, populations of other exotic pests are unable to be characterised at the time. Here we have developed a generic method that is useful across species within the same taxonomic Order, including where there is little or no prior knowledge of their gene sequences. Markers are formed by concomitant sequencing of four gene regions. Sequence concatenation was shown to retrieve higher resolution signatures than standard DNA barcoding. The method is encouragingly universal, as illustrated across species in ten fly and 11 moth superfamilies. Although as-yet untested in a biosecurity situation, this relatively low-tech, off-the-shelf method makes a proactive contribution to the toolbox of quarantine agencies at the time of detection without the need for impromptu species-specific research and development.

Abstract

Biosecurity responses to post-border exotic pest detections are more effective with knowledge of where the species may have originated from or if recurrent detections are connected. Population genetic markers for this are typically species-specific and not available in advance for any but the highest risk species, leaving other less anticipated species difficult to assess at the time. Here, new degenerate PCR primer sets are designed for within the Lepidoptera and Diptera for the 3′ COI, ND3, ND6, and 3′ plus 5′ 16S gene regions. These are shown to be universal at the ordinal level amongst species of 14 and 15 families across 10 and 11 dipteran and lepidopteran superfamilies, respectively. Sequencing the ND3 amplicons as an example of all the loci confirmed detection of population-level variation. This supported finding multiple population haplotypes from the publicly available sequences. Concatenation of the sequences also confirmed that higher population resolution is achieved than for the individual genes. Although as-yet untested in a biosecurity situation, this method is a relatively simple, off-the-shelf means to characterise populations. This makes a proactive contribution to the toolbox of quarantine agencies at the time of detection without the need for unprepared species-specific research and development.

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.

A Novel Genospecies of Borrelia burgdorferi Sensu Lato Associated with Cricetid Rodents in Brazil

Borrelia burgdorferi sensu lato (Bbsl) spirochetes thrive in sylvatic transmission cycles infecting vertebrates and their ticks. Rodents and ticks of the genus Ixodes are important hosts of these spirochetes globally. Although evidence suggests that Borrelia burgdorferi sensu stricto does not exist in South America, genospecies of the group (Bbsl) can be found in this region but have been poorly characterized from a genetic viewpoint, and data on their ecoepidemiology are still incipient. Aiming to detect the natural foci of Borrelia in Brazil, we targeted small mammals inhabiting seven forests fragments during a period of three years (2015–2018). Organs (lung) from two Oligoryzomys rodents over a total of 382 sampled mammals were positive, and we performed a molecular characterization of 10 borrelial genes to achieve a robust analysis. Phylogenetic trees inferred from 16S rRNA, flaB, ospC, and seven MLST loci (clpA, nifS, pepX, pyrG, recG, rlpB, and uvrA) support the characterization of a novel genospecies of Bbsl that we herein name “Candidatus Borrelia paulista” Rp42. Remarkably, “Ca. B. paulista” is phylogenetically related to Borrelia carolinensis, a genospecies that infects Ixodes ticks and cricetid rodents in North America. A previous study performed in the same area identified Ixodes schulzei feeding on Oligoryzomys rodents. Although this tick species could be considered a probable host for this novel Borrelia sp., further research is needed to confirm this hypothesis.

Changes in the epidemiological characteristics of human brucellosis in Shaanxi Province from 2008 to 2020

In the present study, surveys of case numbers, constituent ratios, conventional biotyping, and multilocus sequence typing (MLST) were applied to characterize the incidence rate and epidemiological characteristics of human brucellosis in Shaanxi Province, China. A total of 12,215 human brucellosis cases were reported during 2008-2020, for an annual average incidence rate of 2.48/100,000. The most significant change was that the county numbers of reported cases increased from 36 in 2008 to 84 in 2020, with a geographic expansion trend from northern Shaanxi to Guanzhong, and southern Shaanxi regions; the incidence rate declined in previous epidemic northern Shaanxi regions while increasing each year in Guanzhong and southern Shaanxi regions such as Hancheng and Xianyang. The increased incidence was closely related to the development of large-scale small ruminants (goats and sheep) farms in Guanzhong and some southern Shaanxi regions. Another significant feature was that student cases (n = 261) were ranked second among all occupations, accounting for 2.14% of the total number of cases, with the majority due to drinking unsterilized goat milk. Three Brucella species were detected (B. melitensis (bv. 1, 2, 3 and variant), B. abortus bv. 3/6, and B. suis bv. 1) and were mainly distributed in the northern Shaanxi and Guanzhong regions. Three known STs (ST8, ST2, and ST14) were identified based on MLST analysis. The characteristics that had not changed were that B. melitensis strains belonging to the ST8 population were the dominant species and were observed in all nine regions during the examined periods. Strengthened human and animal brucellosis surveillance and restriction of the transfer of infected sheep (goats) as well as students avoiding drinking raw milk are suggested as optimal control strategies.

Presence of Human Pathogens of the Borrelia burgdorferi sensu lato Complex Shifts the Sequence Read Abundances of Tick Microbiomes in Two German Locations

The distribution of human Lyme borreliosis (LB) is assumed random in Germany, indicating that the human pathogenic species of the Borrelia burgdorferi sensu lato complex (Bb) are similarly distributed as part of the tick microbiome. The aim of this study was to differentiate if the presence of Bb occurs with a defined tick microbiome composition. Furthermore, the effect of location on tick microbiome composition was addressed for two German locations. Therefore, nucleic acid extracts from 82 Borrelia-positive and 118 Borrelia-negative Ixodes ricinus ticks sampled from human hosts in both districts were selected. Nucleic acid extracts were used for human pathogenic Bb species diagnostics based on qPCR and multilocus sequence typing (MLST) and bacterial 16S rRNA gene amplicon sequencing followed by network analyses. As a result, the presence of Bb shifted the sequence read abundances of Candidatus Midichloria, Rickettsia, Pseudomonas, Staphylococcus, and Candidatus Neoehrlichia and their topological roles in the tick microbiome. Moreover, the location was less important in the tick microbiome composition but shifted significantly sequence read abundances of Pseudomonas and Wolbachia as well as the topological role of microbial members. Since the presence of human pathogenic Bb species with other tick-associated pathogens varies regionally, we suggest that a bacterial 16S rRNA gene-based microbiome survey should be implemented in the routine diagnostics for both tick and host if human pathogenic species of Bb were detected. This diagnostic extension will help to optimize therapeutic approaches against Bb infection and co-occurring pathogens.

Surveillance of Borrelia miyamotoi-carrying ticks and genomic analysis of isolates in Inner Mongolia, China

Background

Borrelia miyamotoi is a newly described relapsing fever spirochete transmitted by ixodid tick species. Little is known about the prevalence of B. miyamotoi infections in humans and ticks in Inner Mongolia, China. Therefore, we investigated the prevalence of B. miyamotoi in Ixodes persulcatus ticks, and we aimed to isolateB. miyamotoi from I. persulcatus from four regions of Greater Khingan, Inner Mongolia, China.

Methods

From May to June each year during the period 2016–2019, host-seeking adult I. persulcatus ticks were collected from vegetation. Genomic DNA was prepared from half of each tick body for PCR template, and the remaining half was used to cultivate B. miyamotoi in BSK-M medium. We employed quantitative real-time PCR (qPCR) to detect Borrelia DNA in the ticks and to calculate the prevalence of B. miyamotoi and infections with other borreliae. For characterization of the isolated B. miyamotoi, we performed draft genome sequencing and multilocus sequencing analysis (MLSA).

Results

A total of 2656 adult I. persulcatus ticks were collected. The overall prevalence of relapsing fever (RF) borreliae in ticks was 5.0% (134/2656) and that of Lyme disease (LD) borreliae was 43.8% (1164/2656). Co-infection with RF and LD borreliae was observed in 63 ticks (2.4%). Ticks that were positive for RF borreliae by qPCR were subjected to glycerophosphodiester diester phosphodiesterase gene (glpQ) PCR amplification and sequencing, through which we identified the RF borrelia specimens as B. miyamotoi. Furthermore, the B. miyamotoi strain Hetao-1 was isolated from I. persulcatus, and a draft genome sequence was obtained from the isolate. Sequencing determined the strain Hetao-1 genome to be approximately 906.1 kbp in length (28.9% average GC content), and MLSA identified the strain as ST633, which has previously been reported in Japan and Mongolia.

Conclusion

We detected B. miyamotoi from I. persulcatus ticks collected in Inner Mongolia, and successfully isolated a B. miyamotoi strain. To our knowledge, this is the first study to culture a B. miyamotoi isolate from China. The data on the prevalence of B. miyamotoi and other borreliae in I. persulcatus ticks will be fundamental for future epidemiological studies of B. miyamotoi disease in Inner Mongolia.

Simultaneous Detection and Differentiation of Clinically Relevant Relapsing Fever Borrelia with Semimultiplex Real-Time PCR

Bacterial vector-borne diseases, including Borrelia species, present a significant diagnostic, clinical, and public health challenge due to their overlapping symptoms and the breadth of causative agents and arthropod vectors. The relapsing fever (RF) borreliae encompass both established and emerging pathogens and are transmitted to humans by soft ticks, hard ticks, or lice. We developed a real-time semimultiplex PCR assay that detects multiple RF borreliae causing human illness and classifies them into one of three groups. The groups are based on genetic similarity and include agents of soft-tick relapsing fever (Borrelia hermsii and others), the emerging hard-tick-transmitted pathogen B. miyamotoi, and the agent of louse-borne relapsing fever (B. recurrentis). The real-time PCR assay uses a single primer pair designed to amplify all known pathogenic RF borreliae and multiple TaqMan probes to allow the detection of and differentiation among the three groups. The assay detects all RF borreliae tested, with an analytical limit of detection below 15 genome equivalents per reaction. Thirty isolates of RF borreliae encompassing six species were accurately identified. Thirty-nine of 41 residual specimens (EDTA whole blood, serum, or plasma) from patients with RF were detected and correctly classified. None of 42 clinical samples from patients with other infections and 46 culture specimens from non-RF bacteria were detected. The development of a single-assay real-time PCR approach will help to improve the diagnosis of RF by simplifying the selection of tests to aid in the clinical management of acutely ill RF patients.

Examining Prevalence and Diversity of Tick-Borne Pathogens in Questing Ixodes pacificus Ticks in California

Tick-borne diseases in California include Lyme disease (caused by Borrelia burgdorferi), infections with Borrelia miyamotoi, and human granulocytic anaplasmosis (caused by Anaplasma phagocytophilum). We surveyed multiple sites and habitats (woodland, grassland, and coastal chaparral) in California to describe spatial patterns of tick-borne pathogen prevalence in western black-legged ticks (Ixodes pacificus). We found that several species of Borrelia-B. burgdorferi, Borrelia americana, and Borrelia bissettiae-were observed in habitats, such as coastal chaparral, that do not harbor obvious reservoir host candidates. Describing tick-borne pathogen prevalence is strongly influenced by the scale of surveillance: aggregating data from individual sites to match jurisdictional boundaries (e.g., county or state) can lower the reported infection prevalence. Considering multiple pathogen species in the same habitat allows a more cohesive interpretation of local pathogen occurrence. IMPORTANCE Understanding the local host ecology and prevalence of zoonotic diseases is vital for public health. Using tick-borne diseases in California, we show that there is often a bias to our understanding and that studies tend to focus on particular habitats, e.g., Lyme disease in oak woodlands. Other habitats may harbor a surprising diversity of tick-borne pathogens but have been neglected, e.g., coastal chaparral. Explaining pathogen prevalence requires descriptions of data on a local scale; otherwise, aggregating the data can misrepresent the local dynamics of tick-borne diseases.

Detecting and phasing minor single-nucleotide variants from long-read sequencing data

Cellular genetic heterogeneity is common in many biological conditions including cancer, microbiome, and co-infection of multiple pathogens. Detecting and phasing minor variants play an instrumental role in deciphering cellular genetic heterogeneity, but they are still difficult tasks because of technological limitations. Recently, long-read sequencing technologies, including those by Pacific Biosciences and Oxford Nanopore, provide an opportunity to tackle these challenges. However, high error rates make it difficult to take full advantage of these technologies. To fill this gap, we introduce iGDA, an open-source tool that can accurately detect and phase minor single-nucleotide variants (SNVs), whose frequencies are as low as 0.2%, from raw long-read sequencing data. We also demonstrate that iGDA can accurately reconstruct haplotypes in closely related strains of the same species (divergence ≥0.011%) from long-read metagenomic data.

Genetic diversity of Borrelia burgdorferi sensu stricto: Novel strains from Mexican wild rodents

Borrelia burgdorferi s.s. is a Gram-negative spirochaete, the aetiological agent of Lyme disease, the most common vector-borne disease in the Northern hemisphere. Reports on the presence of B. burgdorferi in central Mexico have been strongly criticized, since these were based only on unspecific serological methods. Furthermore, the worldwide genetic diversity of B. burgdorferi s.s. has not been evaluated. For this reason, the aim of the present study was to confirm the presence of B. burgdorferi in the central area of Mexico and to evaluate its relationship with regard to the global genetic diversity of B. burgdorferi s.s. To achieve this, fragments of the flagellin and the outer surface protein A genes were amplified from ear biopsies of the arboreal wild endemic mice Habromys schmidlyi. With these sequences, a concatenated Bayesian analysis was performed to confirm the identity of B. burgdorferi s.s. Afterwards, the global genetic diversity of this bacterial species was evaluated using our sequences and those available in GenBank. A prevalence of 10.4% (5/48) of H. schmidlyi infected with Borrelia sp. was detected, and the phylogenetic analyses confirmed the identity of B. burgdorferi s.s. Using both genes, the genetic diversity was low. However, genetic structuring analyses revealed that populations of western United States and those from Mexico formed slightly different genetic groups, separated from the populations of the rest of the world. Our study not only confirms the presence of this bacterium in central Mexico, but also shows the most southern record of this bacterium so far. It also highlights the importance of H. schmidlyi as a new potential host of this bacterial species. Our study also provides first genetic data on an incipient process of divergence in B. burgdorferi s.s. populations of eastern United States and central Mexico.

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.

Monitoring of pathogenic Borrelia burgdorferi sensu lato in the Almaty oblast, Kazakhstan

Lyme borreliosis (LB) is one of the most common vector-borne diseases transmitted by ticks. It is caused by the Borrelia burgdorferi sensu lato (s.l.) genospecies. The Almaty oblast of Kazakhstan is considered endemic for LB. Nevertheless, there are significant gaps in the tick surveillance for LB agents in the region. We evaluated B. burgdorferi s.l. prevalence in 413 ixodid ticks collected from bitten people and the seroprevalence of antibodies to B. burgdorferi s.l. in 589 residents of the Almaty oblast, Kazakhstan. All samples were collected between 2018 and 2020. Borrelia burgdorferi s.l. DNA was detected in 24 % (9/38) of I. persulcatus ticks that attacked humans in the city of Almaty and the Talgar and Karasay districts. Multilocus sequence typing identified two novel B. afzelii sequence types in I. persulcatus. The seroprevalence of IgG antibodies against B. burgdorferi s.l. in the analyzed population was detected to be 5.8 %. The obtained results confirm active circulation of B. burgdorferi s.l. in the Almaty oblast of Kazakhstan and raise concern regarding LB preventive measures in the region.

Bats are key hosts in the radiation of mammal-associated Bartonella bacteria

Bats are notorious reservoirs of several zoonotic diseases and may be uniquely tolerant of infection among mammals. Broad sampling has revealed the importance of bats in the diversification and spread of viruses and eukaryotes to other animal hosts. Vector-borne bacteria of the genus Bartonella are prevalent and diverse in mammals globally and recent surveys have revealed numerous Bartonella lineages in bats. We assembled a sequence database of Bartonella strains, consisting of nine genetic loci from 209 previously characterized Bartonella lineages and 121 new cultured isolates from bats, and used these data to perform a comprehensive phylogenetic analysis of the Bartonella genus. This analysis included estimation of divergence dates using a molecular clock and ancestral reconstruction of host associations and geography. We estimate that Bartonella began infecting mammals 62 million years ago near the Cretaceous-Paleogene boundary. Additionally, the radiation of particular Bartonella clades correlate strongly to the timing of diversification and biogeography of mammalian hosts. Bats were inferred to be the ancestral hosts of all mammal-associated Bartonella and appear to be responsible for the early geographic expansion of the genus. We conclude that bats have had a deep influence on the evolutionary radiation of Bartonella bacteria and their spread to other mammalian orders. These results support a 'bat seeding' hypothesis that could explain similar evolutionary patterns in other mammalian parasite taxa. Application of such phylogenetic tools as we have used to other taxa may reveal the general importance of bats in the ancient diversification of mammalian parasites.

Laboratory Diagnosis of Lyme Borreliosis

Lyme borreliosis is caused by a growing list of related, yet distinct, spirochetes with complex biology and sophisticated immune evasion mechanisms. It may result in a range of clinical manifestations involving different organ systems, and can lead to persistent sequelae in a subset of cases. The pathogenesis of Lyme borreliosis is incompletely understood, and laboratory diagnosis, the focus of this review, requires considerable understanding to interpret the results correctly. Direct detection of the infectious agent is usually not possible or practical, necessitating a continued reliance on serologic testing. Still, some important advances have been made in the area of diagnostics, and there are many promising ideas for future assay development. This review summarizes the state of the art in laboratory diagnostics for Lyme borreliosis, provides guidance in test selection and interpretation, and highlights future directions.

Borrelia in neotropical bats: Detection of two new phylogenetic lineages

The genus Borrelia encompasses 50 spirochetal species, several of which are pathogenic and have been detected in a wide range of mammals, especially rodents and cervids. Although the order Chiroptera is the second most diverse mammalian order, and borreliosis represents a human and veterinary health problem in endemic countries, few studies have previously reported infections of Borrelia in these flying mammals. For this reason, the aim of the present study was to detect the presence of, and to analyze the diversity of Borrelia species in several bat species from Mexico. A total of 69 bats belonging to 11 species were collected and molecular detection of Borrelia was performed by amplifying three genes using specific primers. Only five individuals of four bat species (Saccopteryxbilineata, Choeroniscus godmani, Sturnira parvidens and Lasiurus cinereus) tested positive for Borrelia DNA. We now show the first Borrelia record in Mexican bats from two different ecosystems, where previously several potential vector species of the genus Ixodes and Ornithodoros had been reported. The Borrelia sequences obtained from the bats revealed two new putative lineages, one from the relapsing fever group and the second one belonging to the Borrelia burgdorferi s.l. complex, both of which are related to zoonotic species. These results highlight the importance of bats as potential hosts of Borrelia, and the imperative need of active surveillance in flying mammals in order to understand their potential role in the life cycle of this bacteria genus.

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.

Novel Borrelia genotypes in bats from the Macaregua Cave, Colombia

Bats have been implicated as reservoirs of relapsing fever group spirochaetes since the beginning of the last century. Recently, bat-associated spirochaetes have been reported as human pathogens. In 1968, a spirochaete was detected in blood of the bat Natalus tumidirostris captured inside the Macaregua cave, Colombia. Data on this microorganism were never published again. The aim of this study was to evaluate the presence of Borrelia DNA in blood from bats of Macaregua cave. We performed molecular analyses using a genus-specific real-time PCR targeting the 16S rRNA to detect DNA of Borrelia in blood samples from 46 bats captured in the Macaregua cave. Positive samples were submitted to a battery of PCRs aiming to amply Borrelia 16S rRNA, flaB, glpQ, p66, ospC, clpA, clpX, nifS, pepX, pyrG, recG, rplB and uvrA genes. Seventeen samples were positive for Borrelia after real-time PCR. With the exception of flaB gene, attempts to amplify further loci were unsuccessful. Nucleotide and amino acid divergences of four flaB haplotypes characterized from blood of Carollia perspicillata showed Borrelia burgdorferi sensu lato (Bbsl) as the most closely related group. A phylogenetic tree including 74 sequences of the genus confirmed this trend, since Borrelia genotypes detected in bats from Macaregua formed a monophyletic group basally positioned to Bbsl. Our results suggest that Borrelia genotypes characterized from bats roosting in the Macaregua cave might constitute a new taxon within the genus. This is the first molecular characterization of a Borrelia sp. in Colombia.

Detection of Borrelia burgdorferi Sensu Lato and Relapsing Fever Borrelia in Feeding Ixodes Ticks and Rodents in Sarawak, Malaysia: New Geographical Records of Borrelia yangtzensis and Borrelia miyamotoi

Members of the Borrelia burgdorferi sensu lato (Bbsl) complex are etiological agents of Lyme disease (LD), and Borrelia miyamotoi is one of the relapsing fever Borrelia (RFB). Despite the serological evidence of LD in Malaysia, there has been no report from Sarawak, Malaysian Borneo. Thus, this study aimed to detect and characterize Borrelia in rodents and Ixodes ticks from primary forests and an oil palm (OP) plantation in Sarawak. Borrelia yangtzensis (a member of the Bbsl complex) was detected in 43.8% (14/32) of Ixodes granulatus; most of the positive ticks were from the OP plantation (13/14). Out of 56 rodents, B. yangtzensis was detected in four Rattus spp. from the OP plantation and B. miyamotoi was detected in one rodent, Sundamys muelleri, from the primary forest. Further, the positive samples of B. yangtzensis were randomly selected for multilocus sequence analysis (MLSA). The MLSA results of successfully amplified tick samples revealed a clustering with the sequences isolated from Japan and China. This study is the first evidence of B. miyamotoi, a known human pathogen in Malaysia, and B. yangtzensis, which is circulating in ticks and rodents in Sarawak, Malaysian Borneo, and presenting a new geographical record of the Borrelia spp.

Spatial variability in prevalence and genospecies distributions of Borrelia burgdorferi sensu lato from ixodid ticks collected in southern Germany

Lyme borreliosis (LB) is the most common arthropod-borne disease in Europe and North America and is caused by members of the Borrelia burgdorferi sensu lato (Bbsl) species complex. These bacteria are transmitted by ixodid tick vectors and therefore human LB risk is influenced by the prevalence and distribution of Bbsl genospecies within tick vectors throughout the wild. These distributions can easily change over spatiotemporal scales and, to understand LB risk fully, up to date information on prevalence and distribution of Bbsl is required. The last survey of Bbsl in southern Germany, including parts of the Munich metropolitan area, was completed in 2006 and new data is needed. Ixodid ticks were collected in seven plots located in and around Munich, Germany, from March to July 2019 and were screened for Bbsl. Borrelia burgdorferi s. l. positive ticks (52 adults, 158 nymphs) were found in all plots and adults (0-61.5 % Bbsl positive/plot) and nymphs (17.4-59.5 % Bbsl positive/plot) did not differ significantly in their overall Bbsl prevalence. The number of Bbsl positive nymphs did vary significantly between plots but the number of positive adults did not. In total, six Bbsl genospecies were located with B. afzelii and B. garinii dominating. Additionally, the relapsing-fever species B. miyamotoi was found in two sampling plots. Our results highlight the variability in Bbsl prevalence and genospecies distribution over short geographic distances and aid in understanding LB risk in and around the Munich metropolitan area.

Ticks (Acari: Ixodidae) on birds migrating to the island of Ponza, Italy, and the tick-borne pathogens they carry

Seasonal migration of birds between breeding and wintering areas can facilitate the spread of tick species and tick-borne diseases. In this study, 151 birds representing 10 different bird species were captured on Ponza Island, an important migratory stopover off the western coast of Italy and screened for tick infestation. Ticks were collected and identified morphologically. Morphological identification was supported through sequencing a fragment of the 16S mitochondrial gene. In total, 16 captured birds carried ticks from four tick species: Hyalomma rufipes (n = 14), Amblyomma variegatum (n = 1), Amblyomma sp. (n = 1), and Ixodes ventalloi (n = 2). All specimens were either larvae (n = 2) or nymphs (n = 16). All ticks were investigated for tick-borne pathogens using published molecular methods. Rickettsia aeschlimannii was detected in six of the 14 collected H. rufipes ticks. Additionally, the singular A. variegatum nymph tested positive for R. africae. In all 14 H. rufipes specimens (2 larvae and 12 nymphs), Francisella-like endosymbionts were detected. Four H. rufipes ticks tested positive for Borrelia burgdorferi sensu lato in a screening PCR but did not produce sufficient amplicon amounts for species identification. All ticks tested negative for tick-borne encephalitis virus, Crimean-Congo hemorrhagic fever virus, Coxiella burnetii, Coxiella-like organisms, Babesia spp., and Theileria spp. This study confirms the role of migratory birds in the spread and establishment of both exotic tick species and tick-borne pathogens outside their endemic range.