BioMed Research International
Volume 2016 (2016), Article ID 1327320, 7 pages
http://dx.doi.org/10.1155/2016/1327320

Research Article

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

Huixin Liu,1,2 Wei Liu,1,2 Xuexia Hou,1,2 Lin Zhang,1,2 Qin Hao,1,2 and Kanglin Wan1,2

¹State Key Laboratory for Infectious Disease Prevention and Control, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing 102206, China
²Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Hangzhou 310003, China

Received 14 June 2016; Accepted 4 October 2016

Academic Editor: Marcelo A. Soares

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.

Abstract.

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 [7–10]. 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).

Table 1: Distribution of strains in different areas of China.

Areas
Number of isolates

Jilin Province
31

Guangdong Province
2

Inner Mongolia (Neimeng)
11

Shandong Province
1

Liaoning Province
2

Guizhou Province
7

Sichuan Province
6

Heilongjiang Province
9

Xinjiang Uygur Autonomous Region
15

Beijing Municipality
3

Hebei Province
1

Hunan Province
1

Total
89

Table 2: Genotypes of 89 Chinese strains by MLSA.

MLSA genotypes
Number of strains

Borrelia burgdorferi sensu stricto
1

Borrelia garinii
67

Borrelia afzelii
16

Borrelia valaisiana
5

Total
89

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.

Table 3: Base changes (nonsynonymous mutations) in 41 kD flagellin&.

Base change
AA position
Number of changed strains
Distribution

B.b.s.s
B.g
B.a
B.v

GGC(G)-GCC(A)
17
—
—
16
—

GCA(A)-TCA(S)
105
—
—
16
5

TCT(S)-GCT(A)
142
—
—
—
5
ER

AGA(R)-AAA(K)
146
—
—
16
5
ER

GCA(A)-TCA(S)
191
—
—
16
—
ER

TCT(S)-GCT(A)
199
—
—
16
—
ER

ACT(T)-GCT(A)
205
—
66
16
—
ER

ACT(T)-TCT(S)
205
—
—
—
5
ER

GCT(A)-ACT(T)
208
—
66
—
—
ER

GAG(E)-GAC(D)
213
—
1
—
—
ER

GTT(V)-GCT(A)
215
—
66
16
5
ER

CAG(Q)-GAG(E)
216
—
—
16
—
ER

GCA(A)-ACA(T)
224
—
1
16
—
ER

TCT(S)-ACT(T)
230
—
—
16
—
ER

ATA(I)-GTG(V)
260
—
—
—
5
ER

AAT(N)-GAT(D)
279
—
50
—
—
NER

&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].

Table 4: Amino acid changes of the B-cell epitopes included in 41 kD flagellin antigen.

IEDB-ID
Epitopes
Base change

26607
IINHNTSAINASRNNGINAANLSKTQEKLSSGYRIN
GGC(G)-GCC(A)

27732
INRIADQAQY
No change

27733
INRIADQAQYNQMHMLSNKSASQNVRTAEELGMQPAKI
TCT(S)-GCT(A), AGA(R)-AAA(K)

54118
RIADQAQYNQ
No change

752
ADQAQYNQMH
No change

50350
QAQYNQMHML
No change

53012
QYNQMHMLSN
No change

45589
NQMHMLSNKS
No change

41667
MHMLSNKSAS
TCT(S)-GCT(A)

42062
MLSNKSASQN
TCT(S)-GCT(A)

59856
SNKSASQNVR
TCT(S)-GCT(A), AGA(R)-AAA(K)

7895
DEAIAVNIYAANVANLFSGEGAQTAQAAPVQEGVQQE
GCA(A)-TCA(S), TCT(S)-GCT(A),
ACT(T)-G(T)CT(AS), GCT(A)-ACT(T)
GAG(E)-GAC(D), GTT(V)-GCT(A)
CAG(Q)-GAG(E)

35977
LFSGEGAQTAQAAPVQEGVQQEGAQQPAPATAPSQGGVNSPVNVT
TCT(S)-GCT(A), ACT(T)-G(T)CT(A S)GCT(A)-ACT(T), GAG(E)-GAC(D),
GTT(V)-GCT(A), CAG(Q)-GAG(E),
GCA(A)-ACA(T), TCT(S)-ACT(T)

3341
ANLFSGEGAQTAQ
TCT(S)-GCT(A), ACT(T)-G(T)CT(A S)

58007
SGEGAQTAQAAPV
TCT(S)-GCT(A), ACT(T)-G(T)CT(A S),
GCT(A)-ACT(T)

52474
QTAQAAPVQEGVQ
ACT(T)-G(T)CT(A S), GCT(A)-ACT(T)
GAG(E)-GAC(D), GTT(V)-GCT(A)

62967
TAQAAPVQEG
ACT(T)-G(T)CT(A S), GCT(A)-ACT(T)
GAG(E)-GAC(D)

62968
TAQAAPVQEGVQQEGAQQPAPA
ACT(T)-G(T)CT(A S), GCT(A)-ACT(T)
GAG(E)-GAC(D), GTT(V)-GCT(A)
CAG(Q)-GAG(E), GCA(A)-ACA(T)

50213
QAAPVQEGVQ
GCT(A)-ACT(T), GAG(E)-GAC(D)
GTT(V)-GCT(A)

379
AAPVQEGVQQEGA
GCT(A)-ACT(T), GAG(E)-GAC(D)
GTT(V)-GCT(A), CAG(Q)-GAG(E)

3893
APVQEGVQQE
GAG(E)-GAC(D), GTT(V)-GCT(A)
CAG(Q)-GAG(E)

70549
VQEGVQQEGA
GAG(E)-GAC(D), GTT(V)-GCT(A)
CAG(Q)-GAG(E)

70550
VQEGVQQEGAQQP
GAG(E)-GAC(D), GTT(V)-GCT(A)
CAG(Q)-GAG(E)

50610
QEGVQQEGAQQPAP
GAG(E)-GAC(D), GTT(V)-GCT(A)
CAG(Q)-GAG(E), GCA(A)-ACA(T)

12306
EGVQQEGAQQ
GAG(E)-GAC(D), GTT(V)-GCT(A),
CAG(Q)-GAG(E)

12307
EGVQQEGAQQPA
GAG(E)-GAC(D), GTT(V)-GCT(A)
CAG(Q)-GAG(E), GCA(A)-ACA(T)

70628
VQQEGAQQPA
GTT(V)-GCT(A), CAG(Q)-GAG(E),
GCA(A)-ACA(T)

52010
QQEGAQQPAPATA
CAG(Q)-GAG(E), GCA(A)-ACA(T)

50597
QEGAQQPAPA
CAG(Q)-GAG(E), GCA(A)-ACA(T)

51693
QPAPATAPSQGGV
GCA(A)-ACA(T), TCT(S)-ACT(T)

3820
APSQGGVNSPVNV
TCT(S)-ACT(T)

3821
APSQGGVNSPVNVTTTVDANTSLAKIENAIRMISDQRANLGAFQNR
TCT(S)-ACT(T), ATA(I)-GTG(V)

25890
IENAIRMISD
ATA(I)-GTG(V)

43169
NAIRMISDQR
ATA(I)-GTG(V)

28340
IRMISDQRAN
ATA(I)-GTG(V)

41672
MISDQRANLG
ATA(I)-GTG(V)

57265
SDQRANLGAF
No change

52186
QRANLGAFQN
No change

3342
ANLGAFQNRL
No change

36055
LGAFQNRLES
No change

1379
AFQNRLESIK
No change

65265
TMTDEVVAATTNSILTQSAMAMIAQANQVPQYVLSLLR
No 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.

Figure 1: Phylogenetic analysis of 89 strains based on amino acid sequence of 41 kD flagellin.

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 [25–28]. A few studies have demonstrated that VlsE and OspC were critical for immune evasion of B. burgdorferi [1, 29–32]. 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.

Acknowledgments

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).

References

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.

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.

Z. Hubálek, “Epidemiology of lyme borreliosis,” Current Problems in Dermatology, vol. 37, pp. 31–50, 2009.

R. A. Sykes and P. Makiello, “An estimate of Lyme borreliosis incidence in Western Europe,” Journal of Public Health, 2016.

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.

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.

B. Wilske, “Epidemiology and diagnosis of Lyme borreliosis,” Annals of Medicine, vol. 37, no. 8, pp. 568–579, 2005.

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.

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.

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.

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.

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.

G. Stanek, G. P. Wormser, J. Gray, and F. Strle, “Lyme borreliosis,” The Lancet, vol. 379, no. 9814, pp. 461–473, 2012.

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.

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.

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.

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.

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.

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.

R. J. Limberger, “The periplasmic flagellum of spirochetes,” Journal of Molecular Microbiology and Biotechnology, vol. 7, no. 1-2, pp. 30–40, 2004.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.