Emerg Infect DisEmerging Infect. DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention24963721407384413-093110.3201/eid2007.130931DispatchDispatchDiversity of Francisella tularensis Subsp. holarctica Lineages, ChinaFrancisella tularensis Subsp. holarctica, ChinaWangYanhuaPengYaoHaiRongXiaLianxuLiHangZhangZhikaiCaiHongLiangYingShenXiaonaYuDongzhengBirdsellDawnWagnerDavid M.KeimPaulNational Institute for Communicable Disease Control and Prevention; State Key Laboratory for Infectious Disease Prevention and Control; Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Beijing, China (Y. Wang, Y. Peng, R. Hai, L. Xia, Z. Zhang, H. Cai, Y. Liang, X. Shen, D. Yu); Beijing Genomics Institute–Shenzhen, Shenzhen, China (H. Li); Northern Arizona University, Flagstaff, Arizona, USA (D. Birdsell, D.M. Wagner, P. Keim)Address for correspondence: Rong Hai, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing 102206, China; email: hairong@icdc.cn or Paul.Keim@nau.edu7201420711911194

We analyzed 10 isolates of Francisella tularensis subspecies holarctica from China and assigned them to known clades by using canonical single-nucleotide polymorphisms. We found 4 diverse subtypes, including 3 from the most basal lineage, biovar japonica. This result indicates unprecedented levels of diversity from a single region and suggests new models for emergence.

Keywords: Francisella tularensis subsp. holarcticaphylogeographySNPcanSNPChinabacteriatularemialineage

Tularemia is a disease caused by distinct subspecies and phylogenetic groups within the bacterial species Francisella tularensis (1,2). These groups exhibit distinct phylogeographic patterns; F. tularensis subsp. tularensis (type A) is restricted to North America, whereas F. tularensis subsp. holarctica (type B) is found throughout many parts of the Northern Hemisphere (3) and has been reported recently in Tasmania (4). Both subspecies exhibit highly clonal population structures, as determined by phylogenetic analysis using data from multilocus variable number tandem repeat analysis, single-nucleotide polymorphisms (SNPs), and indels (57). The wide geographic distribution and low diversity of F. tularensis subsp. holarctica isolates have been used to argue that this clade is recently emerged and highly fit (3), but the geographic origin of its emergence has not been determined.

F. tularensis subsp. holarctica has been further subdivided by whole-genome sequencing and canonical SNP (canSNP) genotyping into multiple clades (7) (Figure 1). The most basal clade consists of strains assigned to the biovar japonica; this biovar had previously only been reported from Japan (8), but a recent report suggests that it may be found in Turkey (9). The next derived clade (B.2/3) has been described only from 2 isolates from California, USA (7). Isolates from these 2 most basal clades are rare, and apparently geographically restricted, but still provide insights into the origin of F. tularensis subsp. holarctica. The global expansion of the more derived clades is extensive, and closely related isolates are common and widely distributed. The source for emergence of the main type B has been proposed for either North America or Scandinavia, on the basis of the presence of the OSU18 clade isolates in both locations (6,7). However, a sampling bias toward both of these geographic regions has left F. tularensis subsp. holarctica diversity in much of the rest of the world poorly understood. We analyzed 10 isolates of F. tularensis subsp. holarctica from China (10) to determine their placement within the current global phylogeographic framework of this pathogen.

Evolutionary development of Francisella tularensis subsp. holarctica. Previous studies (6,7,11) defined the major lineages of this subspecies on the bases of whole-genome sequences and single-nucleotide polymorphism analysis. The other F. tularensis subspecies and closely-related Francisella species are shown at the root.

The Study

The F. tularensis subsp. holarctica isolates we analyzed were collected over a long period but have been preserved by lyophilization and have been verified every 5 years since they were isolated (Table 1, Appendix). We assigned these isolates into previously defined (6,7) phylogenetic clades and conducted a phylogeographic analysis by using a panel of 12 canSNPs specific for F. tularensis subspecies or clades within F. tularensis subsp. holarctica (Table 2); these canSNPs were obtained from previous reports (6,7). The canSNP analysis was PCR based and performed as described (7). Table 1, Appendix, lists the derived or ancestral allele status for these isolates and for 13 control isolates. These data facilitated the assignment of the 10 F. tularensis subsp. holarctica isolates to major phylogenetic subgroups previously identified within this subspecies (6,7).

<italic>Francisella tularensis</italic> isolates used in study of lineages in China*
Strain IDYear of isolationRegion and countrySubspecies or speciesSourceF.3 910179†B.2 5162‡B.3 470841‡B.4 823672‡B.5 1853655‡B.6 713647‡B.12 109781‡B.16 608245§B.17 1743207§B.19 1373999§B.20 1396082, 1789417§B.23 253120§Major subgroup¶Subgroup#
4101081962Tibet, Chinasubsp. holarcticaHumanGCGACGATCCTCBVB.16
4101091962Tibet, Chinasubsp. holarcticaHumanGCGACGATCCTCBVB.16
4101111964Tibet, Chinasubsp. holarcticaTick (D. everestianus) GCGACGATCCTCBVB.16
FSC022**1950Japansubsp. holarcticaHumanGCGACGATCCTCBVB.16
410107Pre-1992Heilongjiang, Chinasubsp. holarcticaHareGAAATGTGCACCBIB.20
4101161983Xinjiang, Chinasubsp. holarcticaTickGAAATGTGCACCBIB.20
410117Pre-1992Xinjiang, Chinasubsp. holarcticaTickGAAATGTGCACCBIB.20
FSC200**1998Swedensubsp. holarcticaHumanGAAATGTGCACCBIB.20
RC503/ FSC257**1949Russiasubsp. holarcticaTick (D. pictus)GAAATGTGCATABIIIB.23
LVS**UNKRussiasubsp. holarcticaUNKGAAATGTGCATABIIIB.23
4101121971Tibet, Chinasubsp. holarcticaTick (D. everestianus) GAATCGAGACTCBIIB.4
4101131964Tibet, Chinasubsp. holarcticaTick (D. everestianus) GAATCGAGACTCBIIB.4
920607Pre-1992Tibet, Chinasubsp. holarcticaUNKGAATCGAGACTCBIIB.4
OSU18**1978Oklahoma, United Statessubsp. holarcticaBeaverGAATCGAGACTCBIIB.4
MI00-1730**2000Michigan, United Statessubsp. holarcticaHumanGAATCGAGACTCBIIB.4
4101051962Tibet, Chinasubsp. holarcticaHumanGAAATAAGCCTCBIVB.6
FTNF002-00**UNKFrancesubsp. holarcticaUNKGAAATAAGCCTCBIVB.6
URFT1**1997Francesubsp. holarcticaHumanGAAATAAGCCTCBIVB.6
F92**2004Germanysubsp. holarcticaMarmosetGAAATAAGCCTCBIVB.6
OR96-0246**1996Oregon, United Statessubsp. holarcticaUNKGAAATAAGCCTCBIVB.6
Schu S4**1941Ohio, United Statessubsp. tularensisHumanACGACGAGCCTCNANA
FSC147**1965Kazakhstansubsp. mediasiaticaGerbilACGACGAGCCTCNANA
U112**1950Utah, United StatesF. novicidaWaterACGACGAGCCTCNANA

*Clade designations are underlined. Gray shading indicates derived SNPs; other SNPs are ancestral. ID, Identification; UNK, unknown; NA, not applicable; SNP, single-nucleotide polymorphism.
†SNP from (10). SNP position based on the reference Schu S4 (NC_006570)  
‡Canonical SNP from (7). SNP position based on the reference Schu S4 (NC_006570).
§Canonical SNP from (6), SNP position based on the reference Schu S4 (NC_006570).
¶Major subgroup based on (6).
#Subgroup based on (6,7).
**Reference strain genome published in the National Center for Biotechnology Information (www.ncbi.nlm.nih.gov/genome/genomes/511).

Primers used in sequencing to obtain canonical SNP loci for <italic>Francisella tularensis</italic> subsp. <italic>holarctica</italic> isolates*
SNPSCHU S4 SNP position†SNP state (D/A)‡Primer sequence, 5′ → 3′Annealing temperature, °C§
F.3910179G/AF: GCTGTATCATCATTTAATAAACTGCTG55
R: TTGGGAAGCTTGTATCATGGCACT
B.25162A/CF: TTAGTCTATGAGCAGCCAG50
R: TAATATCACCAAGGTAGCC
B.3470841A/GF: ACGCTAGGTGTCTTGGT50
R: CTATATCCGCTCAACAT
B.4823672T/AF: TAGACGCACTGGATTTAGGT53.5
R: AACCATCACGCCACCATAAG
B.51853655T/CF: TGGATCAAACAACCGT50
R: TCTCAAGAGCTGGTGC
B.6713647A/GF: AGTAGTGGTAGCGAGGC53.5
R: TACCGTTAGCCCAACAG
B.12109781T/AF: TACTGCCCAACATAGAG55
R: ATCGTGATAAGGCTGGA
B.16608245T/GF: ATGCTAGCAAATTACCATCAAAAG57
R: AACTCTTCTCGCCATCAACTTCTAT
B.171743207A/CF: CCAAGAGCTAAATTAGCTTCAA53.5
R: TGACCAAGAAGGTAGAGGTATTGGTT
B.191373999A/CF: TTGCTACTGATGGTTTAACT57
R: CAATACGTCACTTATGCAGTGAT
B.201396082, 1789417C/TF: ATGGGTCGGACTATCACATC56
R: ATTATTGTTAAACGGCATCG
B.23253120A/CF: GGCAACAGCAGATTCGTGAG56
R: TGAAAGCAGGTTTAGAAGGACAG

*SNP, single-nucleotide polymorphism; D, derived; A, ancestral; F, forward; R, reverse.
†SNP position based on the reference isolate Schu S4 (NC_006570).
‡Top strand orientation of SCHU S4.
§Sequencing conditions as described in (10).

The isolates were of wide phylogenetic diversity for isolates from a single country. The 10 isolates we analyzed were assigned to 4 distinct phylogenetic clades: 3 were assigned to the basal japonica clade (B.16), 3 to the OSU18 clade (B.4), 3 to the FSC200 clade (B.20), and 1 to clade B.6 (Figure 2; Table 1, Appendix). Two of these clades are very basal (B.16 and B.4; Figure 1), whereas the other 2 are relatively derived (B.6 and B.20). Regardless, these results demonstrate the presence of multiple distinct F. tularensis subsp. holarctica lineages in China. Within China, isolates from the Tibetan plateau in the areas bordering Nepal, Bhutan, India, and central Asia were particularly diverse; all 7 strains assigned to clades B.4, B.6, and B.16 were from this region.

Phylogeography of Francisella tularensis (F.t.) subsp. holarctica. The global distribution of different clades (indicated by colored stars, circles, and circle sections) and their phylogenetic relationships (tree) are shown as described (6,7,11). Stars indicate sequenced reference strains. The phylogenetic positions of the 10 isolates from China (boxes on tree) and their sites of isolation (circles within China) are indicated. The exact lineage of strain 410105 (black circle) was not determined.

The substantial diversity of F. tularensis subsp. holarctica from the Tibetan region provides evidence for an Asian ancestral focus of this subspecies. With the exception of the rare B.2/3 California group, all major basal lineages were represented in this small sample from this region. The center of diversity rationale would suggest that F. tularensis subsp. holarctica diversified in Asia and then spread outward to the rest of the world. The presence of representatives of the basal japonica and OSU18 lineages further argues for ancestral populations in this region. In contrast, the derived and terminal position of the B.20 (FSC200) isolates in northern China suggests that this lineage was introduced to this region from other regions—perhaps Europe, in which B.20 is found (6)—after ancestral strains dispersed to other regions from Asia and diversified in these new locations. The analyses used in this current study show that the B.4 (OSU18) isolates from China are indistinguishable from B.4 isolates from North America or Europe and could represent an ancestral population or a reintroduction after global dissemination.

Although strong evidence shows that F. tularensis subsp. holarctica is a highly fit and recently emerged clone (3), we know little about the basis for its great fitness. It is possible that certain, as yet unidentified, adaptive features developed that led to an increase in its fitness. Alternatively, a stochastic event may have led to the emergence and subsequently circumpolar expansion of this subspecies. However, our understanding of the ecology of F. tularensis subsp. holarctica is severely limited, so the dispersal mechanisms that led to its wide geographic distribution have yet to be identified.

Conclusions

Wide diversity in F. tularensis subsp. holarctica strains, including basal lineages, has been observed in China and underscores a lack of phylogeographic knowledge of this subspecies. Previous arguments (1) about the emergence of this highly fit subspecies have been based on highly biased sampling of strains in North America, Europe, and Japan. Our data suggest a broader distribution in Asia of the japonica clade (B.16) in particular. The OSU18 clade (B.4) also appears to have a broader distribution in Asia than has been observed from both North America and Europe. These clades are thought to be basal to the highly fit clonal expansion on these continents. Sampling of additional regions in Asia and characterization of those isolates would greatly advance the literature on the phylogeography of F. tularensis subsp. holarctica.

Suggested citation for this article: Wang Y, Peng Y, Hai R, Xia L, Li H, Zhang Z, et al. Diversity of Francisella tularensis subsp. holarctica lineages, China. Emerg Infect Dis [Internet]. 2014 Jul [date cited]. http://dx.doi.org/10.3201/eid2007.130931

This work was supported by funding from the China Mega-Project for Infectious Disease (no. 2012ZX10004219 and no. 2011ZX10004-001) and the US Department of Homeland Security Science and Technology Directorate (award HSHQDC-10-C-US).

Dr Wang is a professor and a principal investigator of F. tularensis at the National Institute for Communicable Disease Control and Prevention, Chinese Centers for Disease Control and Prevention. Her research interests include the evolution, epidemiology, and control of bacterial zoonoses.

ReferencesKeim P, Johansson A, Wagner DM. Molecular epidemiology, evolution, and ecology of Francisella. Ann N Y Acad Sci. 2007;1105:3066 10.1196/annals.1409.01117435120Petersen JM, Schriefer ME. Tularemia: emergence/re-emergence. Vet Res. 2005;36:455 10.1051/vetres:200500615845234Keim PS, Wagner DM. Humans and evolutionary and ecological forces shaped the phylogeography of recently emerged diseases. Nat Rev Microbiol. 2009;7:81321 10.1038/nrmicro221919820723Jackson J, McGregor A, Cooley L, Ng J, Brown M, Ong CW, Francisella tularensis subspecies holarctica, Tasmania, Australia, 2011. Emerg Infect Dis. 2012;18:14846 10.3201/eid1809.11185622931809Johansson A, Farlow J, Larsson P, Dukerich M, Chambers E, Bystrom M, Worldwide genetic relationships among Francisella tularensis isolates determined by multiple-locus variable-number tandem repeat analysis. J Bacteriol. 2004;186:580818 10.1128/JB.186.17.5808-5818.200415317786Svensson K, Granberg M, Karlsson L, Neubauerova V, Forsman M, Johansson A. A real-time PCR array for hierarchical identification of Francisella isolates. PLoS ONE. 2009;4:e8360 10.1371/journal.pone.000836020027310Vogler AJ, Birdsell D, Price LB, Bowers JR, Beckstrom-Sternberg SM, Auerbach RK, Phylogeography of Francisella tularensis: global expansion of a highly fit clone. J Bacteriol. 2009;191:247484 10.1128/JB.01786-0819251856Olsufjev NG, Meshcheryakova IS. Subspecific taxonomy of Francisella tularensis McCoy and Chapin 1912. Int J Syst Bacteriol. 1983;33:8724 10.1099/00207713-33-4-872Kiliç S, Çelebi B, Acar B, Ataş M. In vitro susceptibility of isolates of Francisella tularensis from Turkey. Scand J Infect Dis. 2013;45:33741 10.3109/00365548.2012.75112523249114Wang Y, Hai R, Zhang Z, Xia L, Cai H, Liang Y, Genetic relationship between Francisella tularensis strains from China and from other countries. Biomed Environ Sci. 2011;24:310421784318Gyuranecz M, Birdsell DN, Splettstoesser W, Seibold E, Beckstrom-Sternberg SM, Makrai L, Phylogeography of Francisella tularensis subsp. holarctica, Europe. Emerg Infect Dis. 2012;18:2903 10.3201/eid1802.11130522305204