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

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