BioMed Research International

BioMed Research International / 2016 / Article

Research Article | Open Access

Volume 2016 |Article ID 1327320 | 7 pages |

Polymorphism of 41 kD Flagellin Gene and Its Human B-Cell Epitope in Borrelia burgdorferi Strains of China

Academic Editor: Marcelo A. Soares
Received14 Jun 2016
Accepted04 Oct 2016
Published15 Nov 2016


The 41 kD flagellin of Borrelia burgdorferi (B. burgdorferi) is a major component of periplasmic flagellar filament core and a good candidate for serodiagnosis in early stage of Lyme disease. Here, we chose 89 B. burgdorferi strains in China, amplified the gene encoding the 41 kD flagellin, and compared the sequences. The results showed that genetic diversity presented in the 41 kD flagellin genes of all 89 strains among the four genotypes of B. burgdorferi, especially in the genotype of B. garinii. Some specific mutation sites for each genotype of the 41 kD flagellin genes were found, which could be used for genotyping B. burgdorferi strains in China. Human B-cell epitope analysis showed that thirteen of 15 nonsynonymous mutations occurred in the epitope region of 41 kD flagellin and thirty of 42 B-cell epitopes were altered due to all 13 nonsynonymous mutations in the epitope region, which may affect the function of the antigen. Nonsynonymous mutations and changed human B-cell epitopes exist in 41 kD flagellin of B. burgdorferi sensu lato strains; these changes should be considered in serodiagnosis of Lyme disease.

1. Introduction

Lyme disease is a widespread, systemic disease caused by B. burgdorferi, which is transmitted to human by Ixodes ticks [1, 2]. Its incidence varies between countries, with approximately 65,500 patients annually in Europe [3, 4]. Approximately 30,000 new cases of LD occur in the United States each year [5, 6]. Many studies have confirmed that B. burgdorferi is phenotypically and genotypically heterogeneous. To date, 18 B. burgdorferi genospecies have been described; at least four of these species, B. burgdorferi sensu stricto (B.b.s.s), B. garinii, B. afzelii, and B. spielmanii, are associated with human infection [710]. The common causative LD agent in North America is Borrelia burgdorferi sensu stricto (B.b.s.s). In Europe and Asia, Borrelia garinii and Borrelia afzelii are the most abundant species [11, 12]. Clinical manifestations of Lyme disease are diverse, mainly including erythema migrans (EM) skin lesions, acrodermatitis chronica atrophicans, and neurotropic and arthritogenic symptoms. Laboratory evidence of infection, mainly serology, is essential for diagnosis, except in the case of typical EM. Immunological and molecular biological characterization of B. burgdorferi has led to the identification of several antigens that may be useful in the development of improved diagnostic methods and vaccines [13].

The 41 kD flagellin is encoded by the gene flaB and is a major component of B. burgdorferi’s periplasmic flagellar filament core. Recently, a study by Sultan et al. [14] demonstrated that flaB mutant of B. burgdorferi was nonmotile. They also found that whereas wild-type cells were motile and had a flat-wave morphology, flaB mutant cells were nonmotile and rod shaped. Hence, the 41 kD flagellin is critical for optimal survival in ticks and infection of the mammalian host by the arthropod tick vector. Moreover, the first detectable specific immunoglobulin (Ig) M and IgG responses were directed to the 41 kD flagellin in the patients with B. burgdorferi infection [15]. It makes this antigen important for serodiagnosis. The 41 kD flagellin is the most human B-cell epitope-harbored antigen containing forty-four B-cell epitopes [16], which suggest that it may play an important role in the immune reaction between B. burgdorferi and human B-cells. In order to understand the gene diversity of the 41 kD flagellin of Chinese strains, especially on the human B-cell epitopes, we sequenced and analyzed the 41 kD flagellin gene of 89 strains of four genotypes in China.

2. Materials and Methods

2.1. Strain Selection

We selected 89 B. burgdorferi strains which were isolated in Beijing Municipality and 11 provinces and autonomous regions in China (Table 1). These strains were genotyped by multilocus sequence analysis (MLSA) in previous study [17]. 89 strains belonged to four genotypes, which were B.b.s.s, B. garinii, B. afzelii, and B. valaisiana. According to the genotypic ratio of B. burgdorferi strains in China, we selected 1 B.b.s.s strain, 67 B. garinii strains, 16 B. afzelii strains, and 5 B. valaisiana strains in this study (Table 2).

AreasNumber of isolates

Jilin Province31
Guangdong Province2
Inner Mongolia (Neimeng)11
Shandong Province1
Liaoning Province2
Guizhou Province7
Sichuan Province6
Heilongjiang Province9
Xinjiang Uygur Autonomous Region15
Beijing Municipality3
Hebei Province1
Hunan Province1

MLSA genotypesNumber of strains

Borrelia burgdorferi sensu stricto1
Borrelia garinii67
Borrelia afzelii16
Borrelia valaisiana5

2.2. Culture and PCR

These strains were cultured in Barbour-Stoenner-Kelly (BSK) medium, collected by centrifuging at 13000 rpm/min, and then heat-inactivated at 100°C. DNA obtained by this method was used as a template for amplifying the gene of 41 kDa flagellin. The nucleotide sequences of the primers used in this study were designed with Primer 5 software according to B31 genome sequence and were as follows: 5′-TTATCATGGAGGAATGATAT-3′ and 5′-ACCCTACTCAAAGCAAACTC-3′. PCR was performed in a total volume of 50 μL. The PCR mix contained 25 μL PCR buffer, 20 pM of each primer, 2.5 mM each of four dNTPs, and 1 U DNA Taq Polymerase (Takara). Amplification was performed for 10 min of initial denaturation at 94°C; 35 cycles under the following conditions: 1 min of denaturation at 94°C, 1 min of annealing at 52°C, and 1 min of extension at 72°C; 10 min of final extension at 72°C. Negative control (reagent only, no DNA) was included when the PCR was performed. The positive control was 300 ng DNA from the B31 of B.b.s.s genotype, which was the standard strain in the United States. The presence and size of PCR products were determined by electrophoresis on 1.5% agarose gel in Tris-boric acid-EDTA buffer followed by staining with goldview. We performed all of the PCRs at least twice to validate the reproducibility.

2.3. Analytical Methods

The sequences of all PCR products were determined with an ABI 3730xl DNA Analysis. Distances were calculated using the neighbor-joining method. The sequences which contained 46–1011 bp of 41 kD flagellin were compared by MEGA5.10 software [18].

3. Results

We amplified the gene encoding the 41 kDa flagellin, obtained PCR products of all 89 strains, and then compared the sequences based on the 41 kDa flagellin gene sequence of B31. As a result, the nucleotide and amino acids sequences of 41 kD flagellin in B.b.s.s strain CS4 were exactly identical to B31, whereas there were 133 single nucleotide polymorphisms, consisting of 15 nonsynonymous mutations and 118 synonymous mutations in 41 kDa flagellin genes of the remaining 88 strains (Table 3). As shown in Table 3, except AA position 105 and 279 mutations, 13 nonsynonymous mutations located in epitope region. AA positions 205 and 215 displayed higher polymorphism among 15 nonsynonymous mutations. Strains with other three genotypes, B. garinii, B. afzelii, and B. valaisiana, were almost all changed at AA positions 205 and 215. In addition, AA position 205 had two variants, T(Thr)-A(Ala) and T(Thr)-S(Ser). Moreover, the changes of AA positions 208, 213, and 279 were unique to B. garinii strains; AA positions 17, 191, 199, 216, and 230 were unique to B. afzelii strains and AA positions 142 and 260 were unique to B. valaisiana. The specific mutation sites for each genotype could be used for genotyping Chinese B. burgdorferi strains.

Base changeAA positionNumber of changed strainsDistribution


&The CDS of 41 kD flagellin of Borrelia burgdorferi sensu stricto B31 has been used as the reference sequence.
: epitope region; : nonepitope region.

The 41 kD flagellin harbored forty-four human B-cell epitopes [16]. Because two of forty-four human B-cell epitopes were not mapped to B31, we analyzed the changes of forty-two epitopes (Table 4). In our study, thirteen of 15 nonsynonymous mutations occurred in the epitope region of 41 kD flagellin and thirty of 42 human B-cell epitopes were altered due to all 13 nonsynonymous mutations in the epitope region. The changed epitopes mainly concentrated in central region of 41 kD flagellin, which was about AA positions 131–266 [19].

IEDB-IDEpitopesBase change

27732INRIADQAQYNo change
54118RIADQAQYNQNo change
752ADQAQYNQMHNo change
50350QAQYNQMHMLNo change
53012QYNQMHMLSNNo change
45589NQMHMLSNKSNo change
57265SDQRANLGAFNo change
52186QRANLGAFQNNo change
3342ANLGAFQNRLNo change
36055LGAFQNRLESNo change
1379AFQNRLESIKNo change

Note: letters in bold font in the table represent the mutated amino acids in epitopes.

From Table 3 and Figure 1, we can see that the amino acid sequence of strains within the genotype (16 B. afzelii and 5 B. valaisiana) showed almost 100% identity, respectively. Nevertheless, 67 B. garinii strains showed more sequence variations and could be further divided into two major subgroups. Compared with the same genotype strains, three of 67 B. garinii strains, JL13, JC1-13, and JC1-15, showed obvious diversity. The only one strain with mutation at AA position 213 was JC1-13. For the mutation site of AA position 224, JC1-15 was the only one of 67 B. garinii strains.

4. Discussion

In this study, we selected 89 B. burgdorferi strains in China that derived from very large geographical areas and had four MLSA genotyping patterns. Hence, the data provided by them was representative of the genetic diversity that might be presented in China, at least to some extent.

The 41 kD flagellin protein is the predominant component of B. burgdorferi’s flagellar filament core [20, 21]. As reported by Namba et al. [22], X-ray fiber diffraction analysis of the secondary flagellin structure indicated that the terminal portions of the protein influence the structure of the filament and are located toward the center of the filament, whereas the central region of the protein is not involved in conformation and is located on the outside of the filament structure. Some studies have demonstrated that this protein and its central region could elicit the immune response to B. burgdorferi [15, 23, 24]. In addition, the central region of 41 kD flagellin was specific to B. burgdorferi and had no cross-reactive with other Borrelia species, such as Borrelia recurrentis [24]. However, it has been reported that flagellin genes cloned from several B. burgdorferi strains had been shown to be highly conserved [15]. In our study, fifteen nonsynonymous mutations were discovered among 89 strains including four B. burgdorferi genotypes; thirteen of 15 nonsynonymous mutations were found in the epitope region, and some of them were unique to each genotype. It is suggested that the sensitivity including 41 kD flagellin of four genotypes may be improved in serodiagnosis assay.

Studies in human pathogenic viruses, bacteria, and protozoa have revealed that the genes encoding antigens tend to be highly variable as a consequence of diversifying selection to evade host immunity [2528]. A few studies have demonstrated that VlsE and OspC were critical for immune evasion of B. burgdorferi [1, 2932]. The 41 kD flagellin protein containing forty-four human B-cell epitopes may play an important role in humoral immunity between B. burgdorferi and human B-cells. We found that the epitope regions of this protein were diverse; 71.43% (30/42) of the B-cell epitopes presented AA changes, which may reflect ongoing immune evasion. The distribution of these nonsynonymous mutations was located mainly in amino acid positions 142–279, consistent with the central region of 41 kD flagellin antigen.

The protein function may be affected by the changed AA’s polarity. The changes of AA positions 17, 142, 199, and 205 (T-A) caused the decreasing of hydrophilicity. On the contrary, AA positions 105, 191, 208, and 224 led to the increase of hydrophilicity. The changes of AA’s polarity may have a slight impact on the capacity of flagellin penetration, resulting in differences in protein function of each genotype strain.

Based on the diversity of amino acid in 41 kD flagellin, 89 B. burgdorferi strains were clustered into four groups (Figure 1). These four groups by clustering were consistent with MLSA genotyping pattern. As shown in Figure 1, the amino acid sequences of strains were conserved within each genotype of B.b.s.s, B. afzelii, and B. valaisiana, whereas the strains of B. garinii and all strains among four genotypes showed obvious diversity. In particular, the strain of JC1-15 belonging to B. garinii was clustered in Borrelia burgdorferi sensu stricto. However, further investigations will be needed to determine whether the observed changes are due to immune pressure, other selection pressure, or mere random genetic mutation.

5. Conclusion

In China, nonsynonymous mutations and changed human B-cell epitopes exist in 41 kD flagellin of B. burgdorferi sensu lato strains that might affect the function of the antigen and reflect ongoing immune evasion. On the other hand, these changes should be considered in serodiagnosis of Lyme disease.

Competing Interests

The authors declare that they have no competing interests.


The authors thank the staffs of the respective institutes in Beijing Municipality and 11 provinces and autonomous regions in China for their excellent contribution to this study. This work was funded by grants from the National Science and Technology Major Project for Infectious Diseases (2016ZX10004001-004 and 2012ZX10004-215).


  1. D. Brisson, D. Drecktrah, C. H. Eggers, and D. S. Samuels, “Genetics of Borrelia burgdorferi,” Annual Review of Genetics, vol. 46, pp. 515–536, 2012. View at: Publisher Site | Google Scholar
  2. A. C. Steere, J. Coburn, and L. Glickstein, “The emergence of Lyme disease,” Journal of Clinical Investigation, vol. 113, no. 8, pp. 1093–1101, 2004. View at: Publisher Site | Google Scholar
  3. Z. Hubálek, “Epidemiology of lyme borreliosis,” Current Problems in Dermatology, vol. 37, pp. 31–50, 2009. View at: Publisher Site | Google Scholar
  4. R. A. Sykes and P. Makiello, “An estimate of Lyme borreliosis incidence in Western Europe,” Journal of Public Health, 2016. View at: Publisher Site | Google Scholar
  5. B. M. Kuehn, “CDC estimates 300000 US cases of lyme disease annually,” Journal of the American Medical Association, vol. 310, no. 11, p. 1110, 2013. View at: Publisher Site | Google Scholar
  6. A. F. Hinckley, N. P. Connally, J. I. Meek et al., “Lyme disease testing by large commercial laboratories in the United States,” Clinical Infectious Diseases, vol. 59, no. 5, pp. 676–681, 2014. View at: Publisher Site | Google Scholar
  7. B. Wilske, “Epidemiology and diagnosis of Lyme borreliosis,” Annals of Medicine, vol. 37, no. 8, pp. 568–579, 2005. View at: Publisher Site | Google Scholar
  8. X.-B. Ni, N. Jia, B.-G. Jiang et al., “Lyme borreliosis caused by diverse genospecies of Borrelia burgdorferi sensu lato in northeastern China,” Clinical Microbiology and Infection, vol. 20, no. 8, pp. 808–814, 2014. View at: Publisher Site | Google Scholar
  9. D. Richter, D. B. Schlee, R. Allgöwer, and F.-R. Matuschka, “Relationships of a novel Lyme disease spirochete, Borrelia spielmani sp. nov., with its hosts in Central Europe,” Applied and Environmental Microbiology, vol. 70, no. 11, pp. 6414–6419, 2004. View at: Publisher Site | Google Scholar
  10. G. Margos, S. A. Vollmer, N. H. Ogden, and D. Fish, “Population genetics, taxonomy, phylogeny and evolution of Borrelia burgdorferi sensu lato,” Infection, Genetics and Evolution, vol. 11, no. 7, pp. 1545–1563, 2011. View at: Publisher Site | Google Scholar
  11. S. Etti, R. Hails, S. M. Schäfer et al., “Habitat-specific diversity of Borrelia burgdorferi sensu lato in Europe, exemplified by data from Latvia,” Applied and Environmental Microbiology, vol. 69, no. 5, pp. 3008–3010, 2003. View at: Publisher Site | Google Scholar
  12. E. Ozdenerol, “GIS and remote sensing use in the exploration of lyme disease epidemiology,” International Journal of Environmental Research and Public Health, vol. 12, no. 12, pp. 15182–15203, 2015. View at: Publisher Site | Google Scholar
  13. G. Stanek, G. P. Wormser, J. Gray, and F. Strle, “Lyme borreliosis,” The Lancet, vol. 379, no. 9814, pp. 461–473, 2012. View at: Publisher Site | Google Scholar
  14. S. Z. Sultan, A. Manne, P. E. Stewart et al., “Motility is crucial for the infectious life cycle of borrelia burgdorferi,” Infection and Immunity, vol. 81, no. 6, pp. 2012–2021, 2013. View at: Publisher Site | Google Scholar
  15. W. Jiang, B. J. Luft, W. Schubach, R. J. Dattwyler, and P. D. Gorevic, “Mapping the major antigenic domains of the native flagellar antigen of Borrelia burgdorferi,” Journal of Clinical Microbiology, vol. 30, no. 6, pp. 1535–1540, 1992. View at: Google Scholar
  16. R. Vita, L. Zarebski, J. A. Greenbaum et al., “The immune epitope database 2.0,” Nucleic Acids Research, vol. 38, supplement 1, pp. D854–D862, 2009. View at: Publisher Site | Google Scholar
  17. Q. Hao, X. Hou, Z. Geng, and K. Wan, “Distribution of Borrelia burgdorferi sensu lato in China,” Journal of Clinical Microbiology, vol. 49, no. 2, pp. 647–650, 2011. View at: Publisher Site | Google Scholar
  18. K. Tamura, D. Peterson, N. Peterson, G. Stecher, M. Nei, and S. Kumar, “MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods,” Molecular Biology and Evolution, vol. 28, no. 10, pp. 2731–2739, 2011. View at: Publisher Site | Google Scholar
  19. C. Rasiah, E. Schiltz, J. Reichert, and A. Vogt, “Purification and characterization of a tryptic peptide of Borrelia burgdorferi flagellin, which reduces cross-reactivity in immunoblots and ELISA,” Journal of General Microbiology, vol. 138, no. 1, pp. 147–154, 1992. View at: Publisher Site | Google Scholar
  20. R. J. Limberger, “The periplasmic flagellum of spirochetes,” Journal of Molecular Microbiology and Biotechnology, vol. 7, no. 1-2, pp. 30–40, 2004. View at: Publisher Site | Google Scholar
  21. C. W. Wolgemuth, N. W. Charon, S. F. Goldstein, and R. E. Goldstein, “The flagellar cytoskeleton of the spirochetes,” Journal of Molecular Microbiology and Biotechnology, vol. 11, no. 3–5, pp. 221–227, 2006. View at: Publisher Site | Google Scholar
  22. K. Namba, I. Yamashita, and F. Vonderviszt, “Structure of the core and central channel of bacterial flagella,” Nature, vol. 342, no. 6250, pp. 648–654, 1989. View at: Publisher Site | Google Scholar
  23. R. Wallich, S. E. Moter, M. M. Simon, K. Ebnet, A. Heiberger, and M. D. Kramer, “The Borrelia burgdorferi flagellum-associated 41-kilodalton antigen (flagellin): molecular cloning, expression, and amplification of the gene,” Infection and Immunity, vol. 58, no. 6, pp. 1711–1719, 1990. View at: Google Scholar
  24. J. M. Robinson, T. J. Pilot-Matias, S. D. Pratt, C. B. Patel, T. S. Bevirt, and J. C. Hunt, “Analysis of the humoral response to the flagellin protein of Borrelia burgdorferi: cloning of regions capable of differentiating Lyme disease from syphilis,” Journal of Clinical Microbiology, vol. 31, no. 3, pp. 629–635, 1993. View at: Google Scholar
  25. P. Farci, A. Shimoda, A. Coiana et al., “The outcome of acute hepatitis C predicted by the evolution of the viral quasispecies,” Science, vol. 288, no. 5464, pp. 339–344, 2000. View at: Publisher Site | Google Scholar
  26. D. C. Jeffares, A. Pain, A. Berry et al., “Genome variation and evolution of the malaria parasite Plasmodium falciparum,” Nature Genetics, vol. 39, no. 1, pp. 120–125, 2007. View at: Publisher Site | Google Scholar
  27. R. Urwin, J. E. Russell, E. A. L. Thompson, E. C. Holmes, I. M. Feavers, and M. C. J. Maiden, “Distribution of surface protein variants among hyperinvasive meningococci: implications for vaccine design,” Infection and Immunity, vol. 72, no. 10, pp. 5955–5962, 2004. View at: Publisher Site | Google Scholar
  28. Y. Jiang, X. Dou, W. Zhang et al., “Genetic diversity of antigens Rv2945c and Rv0309 in Mycobacterium tuberculosis strains may reflect ongoing immune evasion,” FEMS Microbiology Letters, vol. 347, no. 1, pp. 77–82, 2013. View at: Publisher Site | Google Scholar
  29. S. J. Norris, “vls Antigenic variation systems of lyme disease borrelia: eluding host immunity through both random, segmental gene conversion and framework heterogeneity,” Microbiology Spectrum, vol. 2, no. 6, 2014. View at: Publisher Site | Google Scholar
  30. A. S. Rogovskyy and T. Bankhead, “Variable VlsE is critical for host reinfection by the Lyme disease spirochete,” PLoS ONE, vol. 8, no. 4, Article ID e61226, 2013. View at: Publisher Site | Google Scholar
  31. S. V. Seemanapalli, Q. Xu, K. McShan, and F. T. Liang, “Outer surface protein c is a dissemination-facilitating factor of borrelia burgdorferi during mammalian infection,” PLoS ONE, vol. 5, no. 12, Article ID e15830, 2010. View at: Publisher Site | Google Scholar
  32. Q. Xu, K. McShan, and F. T. Liang, “Identification of an ospC operator critical for immune evasion of Borrelia burgdorferi,” Molecular Microbiology, vol. 64, no. 1, pp. 220–231, 2007. View at: Publisher Site | Google Scholar

Copyright © 2016 Huixin Liu et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

More related articles

717 Views | 267 Downloads | 0 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19. Sign up here as a reviewer to help fast-track new submissions.