The pathogenic intestinal spirochaete Brachyspira pilosicoli forms a diverse recombinant species demonstrating some local clustering of related strains and potential for zoonotic spread

Introduction

The genus Brachyspira includes seven officially named and several unofficially named species of anaerobic spirochaetes that colonize the large intestine of mammals and birds [1]. The three most commonly reported pathogenic species are Brachyspira hyodysenteriae, the agent of swine dysentery, Brachyspira intermedia, a pathogen mainly of adult chickens, and Brachyspira pilosicoli, the cause of a condition that has been called ‘intestinal spirochaetosis’. B. pilosicoli has a broader host range than the other two main pathogenic species, colonizing various species of mammals and birds, as well as human beings [2,3].

Infections with B. pilosicoli are particularly common in intensively housed pigs and chickens, in which they cause depressed rates of growth and production. Colonization also commonly occurs in human beings living in crowded and unhygienic conditions in developing countries [4-7], as well as amongst homosexual males [8]. Individuals colonized with B. pilosicoli may develop focal colitis and chronic diarrhoea, with abdominal pain, failure to thrive and rectal bleeding. An in vitro study using Caco-2 cells has shown that B. pilosicoli strains initially target the cell junctions, where one cell end of the spirochaete invaginates into the Caco-2 cell membranes [9]. The whole cell surface progressively becomes colonized by attached spirochaetes, forming a “false brush border”. In this model colonized monolayers demonstrated accumulation of actin at the cell junctions, loss of tight junction integrity, condensation and fragmentation of nuclear material consistent with apoptosis, and a significant up-regulation of interleukin-1beta and interleukin-8 expression. Besides colitis, a spirochaetaemia with B. pilosicoli has been recorded in immunocompromised or debilitated human beings [10,11], and systemic spread involving the liver also has been described in experimentally infected chickens [12].

B. pilosicoli may be found in water contaminated with faeces and on foodstuffs, and hence has potential importance as a water-borne or food-borne zoonotic pathogen [5,13,14].

Earlier studies using multilocus enzyme electrophoresis (MLEE) have suggested that B. pilosicoli is a recombinant species [2], and that cross-species transmission is likely to occur [3]. A similar MLEE study with the related B. hyodysenteriae also suggested that this species is recombinant with an epidemic population structure [15]; however, more recent studies using multilocus sequence typing (MLST) have indicated that B. hyodysenteriae has a clonal population structure [16,17], as does B. intermedia [18]. These results now have left some uncertainty about the likely population structure of B. pilosicoli.

In an earlier genus-wide MLST study of Brachyspira species, 12 strains of B. pilosicoli were included in the analysis [19]; however, this was too few to deduce the population structure, and there have been no subsequent reports where MLST has been used to analyze B. pilosicoli isolates. Consequently the overall aim of the current study was to apply the previously developed but incomplete Brachyspira MLST system to a large and diverse collection of B. pilosicoli strains to improve understanding of diversity, population structure, host-specificity and geographical links between strains.

Results and discussion

In this study 131 B. pilosicoli strains and isolates from various countries and animal species that had been collected over three decades were used in an MLST scheme for B. pilosicoli. Between 50-112 alleles were identified at the seven MLST loci tested, and a total of 127 sequence type (ST) profiles were obtained (ST01 to ST127). These results demonstrated that high rates of genetic variation occur within the population. The data are summarized in Table 1, with allelic profiles for individual strains shown in Additional file 1: Table S1. The raw sequence data were deposited at the PubMLST site (http://pubmlst.org/ Brachyspira/). After the translation of nucleotides into amino acids, 16-72 alleles were identified at the various loci and 123 amino acid type (AAT) profiles were present (Table 1).

The mean genetic diversity (h value) was 0.977, with diversity at the individual loci varying from 0.913 to 0.989 (Table 1). The extensive diversity that was identified in the population agreed with the results of the earlier MLEE study on B. pilosicoli [2]. Multilocus variable number tandem repeat analysis of B. pilosicoli also has shown considerable diversity, but the frequent occurrence of null alleles limits the use of the technique for detailed analysis of relationships between isolates [20].

The results of the Shimodaira-Hasegawa (SH) test for the seven loci are recorded in Additional file 2: Table S2, and they indicate that each tree had the best topology to explain the genetic relationship of the loci tested. The 35 concatenated trees constructed by using different combinations of three alleles were distinctively different from each another. Results for the four trees that showed the greatest difference are presented in Additional file 3: Table S3. These SH tests indicated that there is substantial recombination in the evolutionary history of B. pilosicoli, and that each gene analyzed was independently evolving. Thus, for each gene there was a significant difference in the Δ - ln L values of each tree and, furthermore, for each of the seven genes the maximum likelihood (ML) trees were no more similar in likelihood than the 200 random trees for each data set. Hence significant phylogenetic incongruence was revealed, implying that frequent recombination has obscured phylogenetic signals expected from direct inheritance of genes in the population.

The standardized index of association (I^S_A) value was calculated as 0.1568 (P < 0.0005), with a small but significant linkage disequilibrium being present in the population. The values are listed in Table 2. Despite the evidence for the population being recombinant, the value suggested that there was a limited degree of clonality within the species that was not masked by high rates of recombination. Consistent with this, all trees that were constructed showed deep branching but with a few small clusters of related isolates (see Additional file 4: Figure S1 as an example). Clustering could be most easily seen in a minimum spanning (MS) tree, which also is marked in colour to show the species of origin (Figure 1) and geographical origin (Figure 2).

Table 1 Number of alleles, genetic diversity, GC content, and variable sites at the seven loci tested

The pathogenic intestinal spirochaete Brachyspira pilosicoli forms a diverse recombinant species demonstrating some local clustering of related strains and potential for zoonotic spread - Image 1

The MS tree was divided into two main but linked parts, with more strains and STs located on the right hand side than on the left. Generally isolates from the same species and geographical origins were distributed throughout the tree, although all the isolates from chickens were located on the right hand side. This distribution might have been influenced by the relatively restricted range of isolates available for analysis, and more robust results will be obtained when more B. pilosicoli strains from different hosts and geographic areas are added to the PubMLST database. There was a limited degree of clustering of isolates (nine clusters: defined as isolates with allelic differences at only one or two loci), and in all but one case where this occurred the clustered isolates were from the same species and from the same geographical origin or farm/community. Hence they are likely to have represented a clonal group that has been transmitted locally. The exception was the largest cluster around ST68 that consisted of isolates from dispersed geographical locations in Australia, and from different host species. This cluster consisted of a isolate from a human child in the Kimberley region in the north of Western Australia, an isolate from a pig in Victoria, two isolates from pigs in the same piggery in the southwest of Western Australia, and two isolates from chickens in Queensland. The occurrence of isolates from one cluster in different species does suggest the possibility of recent cross-species transmission, although it is unlikely to have occurred recently in this case due to the wide geographical distances between the sites where the isolates originated. Possible mechanisms would be transmission through migratory bird species, or mechanical transmission associated with human activities.

By contrast to B. pilosicoli, when using MLST the species B. hyodysenteriae and B. intermedia both have been deduced to be essentially clonal [16,18]. Hence these three important pathogenic species in the same genus have different population structures. One interpretation could be that the latter two species have evolved relatively recently from single stable strains or clones that were derived from a highly recombinant ancestral species such as B. pilosicoli, and which have been successful in finding suitable specialized niches in specific host species. Another possibility could be that the recombinant B. pilosicoli developed from a more stable clonal ancestor following development or acquisition of more effective means for gene transfer and recombination.

The source of the variation amongst Brachyspira species and strains is of considerable interest. Based on the high degree of conservation in the 16S rDNA sequences of the Brachyspira species it has been suggested that they have evolved relatively recently [1]. The location of the seven loci used in the MLST scheme mapped on seven complete Brachyspira genomes is shown as Figure 3. Not only do the relative positions of the loci vary greatly between species, but there are also remarkable differences between the locations in the four sequenced B. pilosicoli strains. These extensive genomic rearrangements within and across species demonstrate the plasticity of Brachyspira genomes.

Table 2 Index of association values generated in the START2 program

Figure 1 Minimum spanning tree showing the MLST profiles of 131 Brachyspira pilosicoli strains with the host species of origin marked. Each node corresponds to a sequence type (ST). The lines between STs show inferred phylogenetic relationships and are represented by bold, continuous, continuous thin, dashed and dotted lines according to the number of allelic mismatches between profiles (1, 2, 3, 4 and 5 or more, respectively). Host species of origin are indicated with coloured text (human (red circle), pig (black circle), chicken (blue circle), dog (violet circle), horse (green circle)).

Click here to enlarge the image

One potential means for recombination may be the activity of bacteriophage-like gene transfer agents (GTA) that have been detected in various Brachyspira species, and which have the potential to facilitate gene transduction within or possibly even across species [21-23]. In addition, in the case of B. pilosicoli, recent analysis of the genomes of three sequenced strains identified genome rearrangements that largely correlated with the positions of mobile genetic elements [24]. Novel bacteriophages also were detected in the newly sequenced genomes, and clearly such genetic elements may have the potential to transduce genetic information and contribute to the recombination that has been recorded here. Interestingly the sizes of the genomes of three sequenced B. pilosicoli strains (B2904, WesB and 95/1000) were ~2,765, 2.890 and 2.596 Mb, respectively [24], while the genome of strain P43/6/78^T has been recorded as 2.56 Mb [25]. This variation in genome size with accompanying loss or gain of genes provides clear evidence for the genomic plasticity of B. pilosicoli.

Conclusions

This study has confirmed that B. pilosicoli has a strongly recombinant population structure that contrasts to the more highly clonal population structures of the related pathogenic species B. hyodysenteriae and B. intermedia. Brachyspira species showed evidence of extensive rearrangement of MLST loci on their genomes, including across four previously sequenced B. pilosicoli strains. The greater genomic plasticity of the recombinant B. pilosicoli may help to explain why it can colonize the large intestines of a wider range of hosts compared to other Brachyspira species. The MLST system that was used was sufficiently sensitive to be able to detect a number of instances where closely related strains (clones) of B. pilosicoli were present in individual animals or people from the same farms or communities, as well as providing evidence for the potential for cross-species and zoonotic transmission by related B. pilosicoli strains.

Figure 2 Minimum spanning tree showing the origin of the B. pilosicoli strains. The country of isolation is shown in coloured text (Australia (red circle), Sweden (black circle), Papua New Guinea (blue circle), USA (violet circle), Canada (green circle), Italy (light blue circle), France (yellow circle), UK (olive green circle), New Zealand (orange circle)).

Click here to enlarge the image

Figure 3 Genome maps of B. hyodysenteriae WA1, B. intermedia PWS/A, B. murdochii DSM 12563 and the four publically available complete Brachyspira pilosicoli genomes (95/1000; B2904; P43/6/78T ; WesB) showing the relative positions of the seven genes targeted for MLST.

Click here to enlarge the image

Methods

Brachyspira pilosicoli strains and isolates

A total of 119 well-characterized strains and isolates of B. pilosicoli were obtained as frozen stock from the culture collection at the Reference Centre for Intestinal Spirochaetes at Murdoch University. They originated from different States of Australia (n = 66), Papua New Guinea (n = 29), the United States of America (n = 8), Canada (n = 5), Italy (n = 5), the United Kingdom (n = 3), France (n = 2) and New Zealand (n = 1). Sequence data for 12 Scandinavian and European strains (AN4170/01, AN991/02, AN76/92, AN497/93, C62, AN984/03, AN1085/02, AN652/ 02, AN2248/02, AN738/02, AN953/02 and C162) that were previously used in a Brachyspira genus-wide MLST study [19] were obtained from the PubMLST website (http://pubmlst.org/) and were included in this analysis.

The full collection, representing 131 isolates, came from a range of species, and consisted of 58 from pigs, 44 from human beings, 24 from chickens, five from dogs and two from horses. The names and origins of the isolates are listed in Additional file 1: Table S1. The identity of the isolates was confirmed using a species-specific PCR for B. pilosicoli incorporating 16S rDNA primers [26].

Spirochaete culture and DNA extraction

The spirochaetes were propagated at 37°C for 5 days in Kunkle’s pre-reduced anaerobic broth containing 2% foetal bovine serum and a 1% ethanolic cholesterol solution [27]. Cells were harvested from culture by centrifuging at 10,000 g, and counted in a haemocytometer chamber under a phase contrast microscope at 40 times magnification.

For each strain, 10 ml of Trypticase Soy broth containing ~10⁸ cells/ml of B. pilosicoli was centrifuged at 5000 g. The supernatant was discarded and the pellet resuspended in an equal volume of phosphate buffered saline (pH 7.4) and heated at 95°C for 15 min to release the DNA, before storing at -20°C. The solution containing the extracted DNA was used as the template for the PCR reactions.

Multilocus sequence typing (MLST)

The seven loci used in MLST were the same as those previously described for use with members of the genus Brachyspira [19]. These were the genes encoding for the conserved “housekeeping” genes alcohol dehydrogenase (adh), alkaline phosphatase (alp), esterase (est), glutamate dehydrogenase (gdh), glycerol kinase (glpK), phosphoglucomutase (pgm), and acetyl-CoA acetyltransferase (thi). The PCR primers used were the same as those used previously [16]. To confirm the conservation of these loci their positions were plotted on genomes of the available single strains of B. hyodysenteriae (WA1), B. intermedia (PWS/A^T ), B. murdochii (DSM12563) and four strains of B. pilosicoli (95/1000; B2904; WesB; P43/6/78^T) [24,28-30].

PCR was performed on DNA from all the B. pilosicoli isolates, using 0.2 μl Taq DNA polymerase, 5 μl of 10× PCR buffers, 3 μl of 25 mM MgCl₂, 5 μl of 8 mM dNTP, 5 μl of the forward and reverse primers, 12 μl of cresol red solution and 2 μl of template, with the reaction mix topped up with PCR-grade water to a final 50 μl volume. Each PCR reaction set included DNA from B. pilosicoli strain 95/1000 as a positive control and distilled water as a negative control. The PCR conditions were 95°C for 2 min, followed by 33 cycles at 95°C for 30 sec, 50°C for 15 sec, 72°C for 1 min for every 1 kbp of product, and a final extension period of 5 min at 72°C before holding at 14°C.

The PCR products were subjected to electrophoresis in a 1% agarose gel in a Bio-Rad Sub-Cell® GT Agarose Gel electrophoresis unit at 120 V for 30 min. A 1 Kbp ladder marker was added to the first and last well of each row to allow estimates of the molecular masses of the samples. The gel was stained by immersion in an ethidium bromide solution at a concentration of 0.5 μg/ml for 30 mins, and the DNA was visualized over a UV illuminator (Biorad Chem Doc XRS Universal Hood).

For sequencing, the PCR products were purified with the Wizard® SV Gel and PCR Clean-Up System Kit (Promega) following the manufacturer’s instructions, then PCR was performed on the purified products using the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, USA) in a 96 well plate, according to the manufacturer’s instructions, using 10 μl of a single primer instead of 5 μl of both forward and reverse primer in each reaction. The amplified products were purified using ethanol precipitation and the pellet was held at 4°C before being sequenced with the ABI 373A sequencing system (Applied Biosystems).

Analysis

The sequences were analyzed and assembled using the Bioedit Sequence Alignment Editor [31]. The sequences for each locus were aligned using the ClustalW program (EMBL-EBI, European Bioinformatics Institute [http:// www.ebi.ac.uk/Tools/msa/clustalw2/]) and B. pilosicoli strain 95/1000 sequence as the standard for the process.

Two methods were used to generate phylogenetic trees. In the first the aligned sequences for each of the seven loci were analyzed using the non-redundant databases (NRDB) program (http://pubmlst.org/analysis/) to identify strain sequences that were identical to each other. Each unique nucleotide sequence was then assigned with a different allele number. The allelic profile for each isolate was determined and consisted of a line listing the allele number for each locus in turn. Isolates were assigned a sequence type (ST) according to their allelic profiles. Isolates were considered genetically identical and belonging to the same ST if their sequences were identical at all seven loci.

The allelic profile was then entered into the dataset of the START2 program and rooted phylogenetic trees (“consensus trees”) with 1000 bootstrap replicates were generated from the data matrix using the ‘Unweighted Pair Group Method with Arithmetic Mean’ (UPGMA) and ‘Neighbour-Joining’ (NJ) method with the ‘Maximum likelihood’ (ML) models [32]. A minimum spanning tree also was generated using the Bionumerics Software (version 7.1, Applied Maths), and colour coded according to the species of origin and geographical origin of the isolates.

The allelic profile was used to calculate genetic diversity, as previously described [33]. To help determine whether recombination had occurred within the B. pilosicoli population, the START2 program was used to estimate the degree of linkage disequilibrium in the population by calculating the index of association (I_A) and the standardized index of association (I^S_A) [34].

The second method of generating phylogenetic trees was by concatenating the nucleotide sequences for the seven genes of each isolates in the order adh, pgm, est, glp, gdh, thi and alp (the same order previously used for other Brachyspira species). All sequences were placed in a single FASTA formatted file and aligned with ClustalW before being converted to the MEGA format (http://ccg.murdoch.edu.au/tools/ clustalw2mega/). UPGMA and NJ trees were constructed from the aligned DNA sequences using the MEGA v4.0.2 program [35].

To verify the topology of the phylogenetic trees, the Shimodaira-Hasegawa (SH) test was carried out using the Phylip v3.69 program [36] to detect significant differences amongst the trees extrapolated for each gene. This analysis was carried out by estimating the maximum likelihood (ML) trees for each of the seven genes, and then comparing, in turn, the difference in log likelihood (Δ - ln L) between each of the seven topologies on each of the seven genes. Randomization tests were used to assess the extent of congruence amongst the seven ML gene trees, and the Δ - ln L values for each of the seven ML trees fitted to each of the seven genes were compared to the equivalent values computed for 200 random trees created from each gene.

To test whether the genetic variation at different loci were independent of one another, 35 ML trees were constructed from concatenated sequences of sets of three random loci for the 127 STs. It was expected that if there were associations between the loci, there would be little variation between the different trees.

Additional files

Additional file 1: Table S1. The names of the 131 isolates, species from which they were isolated, the country of isolation, the sequence type (ST) to which they were assigned in the study and the allelic number assigned to the seven loci. The shaded boxes represent nine sets of isolates in different adjacent STs that differ at only one or two loci, and were each defined as a cluster.

Additional file 2: Table S2. Results of the Shimodaira-Hasegawa test for the seven loci.

Additional file 3: Table S3. Results of the Shimodaira-Hasegawa test on the four concatenated trees that showed the greatest difference with combinations of three loci.

Additional file 4: Figure S1. Neighbour joining tree using the consensus sequences of the 131 B. pilosicoli isolates. A few localized clusters of isolates can be seen, with the largest being ST68 – ST73.

Competing interests

All authors declare that they have no conflicts of interests.

Authors’ contributions

EN, TL and DJH designed the experiments and they were performed by EN, NDP, TL and MYA. EN, TL and DJH analysed the data and DJH and EN wrote the manuscript. All authors approved the manuscript for publication.

Acknowledgements

The research was supported by Murdoch University.

Author details

1 School of Veterinary and Life Sciences, Murdoch University, Murdoch, 6150 Western Australia, Australia. 2 Faculty of Medicine, Hamadan University of Medical Sciences, Hamadan, Iran.

Received: 11 July 2013

Accepted: 14 August 2013 Published: 16 August 2013

This article was originally published in Gut Pathogens 2013, 5:24. DOI: 10.1186/1757-4749-5-24. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0).

References

1. Stanton TB: The genus Brachyspira. In The Prokaryotes. Volume 7. Edited by Falkow S, Rosenberg SE, Schleifer KH, Stackebrant E. New York: Springer; 2006:330–356.

2. Trott DJ, Mikosza ASJ, Combs BG, Oxberry SL, Hampson DJ: Population genetic analysis of Serpulina pilosicoli and its molecular epidemiology in villages in the Eastern Highlands of Papua New Guinea. Int J Syst Bacteriol 1998, 48:659–668.

3. Hampson DJ, Oxberry SL, La T: Potential for zoonotic transmission of Brachyspira pilosicoli. Emerg Infect Dis 2006, 12:869–870.

4. Trott DJ, Combs BG, Oxberry SL, Mikosza ASJ, Robertson ID, Passey M, Taime J, Sehuko R, Hampson DJ: The prevalence of Serpulina pilosicoli in humans and domestic animals in the Eastern Highlands of Papua New Guinea. Epidemiol Infect 1997, 119:369–379.

5. Margawani KR, Robertson ID, Brooke CJ, Hampson DJ: Prevalence, risk factors and molecular epidemiology of Brachyspira pilosicoli in humans on the island of Bali, Indonesia. J Med Microbiol 2004, 53:325–332.

6. Munshi MA, Traub RJ, Robertson ID, Mikosza ASJ, Hampson DJ: Colonization and risk factors for Brachyspira aalborgi and Brachyspira pilosicoli in humans and dogs on tea-estates in Assam, India. Epidemiol Infect 2004, 132:137–144.

7. Nelson EJ, Tanudra A, Chowdhury A, Kane AV, Qadri F, Calderwood SB, Coburn J, Camilli A: High prevalence of spirochetosis in cholera patients, Bangladesh. Emerg Infect Dis 2009, 15:571–573.

8. Trivett-Moore NL, Gilbert GL, Law CLH, Trott DJ, Hampson DJ: Isolation of Serpulina pilosicoli from rectal biopsy specimens showing evidence of intestinal spirochetosis. J Clin Microbiol 1998, 36:261–265.

9. Naresh R, Song Y, Hampson DJ: The intestinal spirochete Brachyspira pilosicoli attaches to cultured Caco-2 cells and induces pathological changes. PLoS ONE 2009, 4(12):e8352.

10. Trott DJ, Jensen NS, Saint Girons I, Oxberry SL, Stanton TB, Lindquist D, Hampson DJ: Identification and characterization of Serpulina pilosicoli isolates recovered from the blood of critically ill patients. J Clin Microbiol 1997, 35:482–485.

11. Bait-Merabet L, Thille A, Legrand P, Brun-Buisson C, Cattoir V: Brachyspira pilosicoli bloodstream infections: case report and review of the literature. Ann Clin Microbiol Antimicrob 2008, 7:19.

12. Mappley LJ, Tchórzewska MA, Nunez A, Woodward MJ, La Ragione RM: Evidence for systemic spread of the potentially zoonotic intestinal spirochaete Brachyspira pilosicoli in experimentally challenged laying chickens. J Med Microbiol 2013, 62:297–302.

13. Oxberry SL, Trott DJ, Hampson DJ: Serpulina pilosicoli, water birds and water: potential sources of infection for humans and other animals. Epidemiol Infect 1998, 121:219–225.

14. Verlinden M, Pasmans F, Garmyn A, De Zutter L, Haesebrouck F, Martel A: Occurrence of viable Brachyspira spp. on carcasses of spent laying hens from supermarkets. Food Microbiol 2012, 32:321–324.

15. Trott DJ, Oxberry SL, Hampson DJ: Evidence for Serpulina hyodysenteriae being recombinant, with an epidemic population structure. Microbiology 1997, 143:3357–3365.

16. La T, Phillips ND, Harland BL, Wanchanthuek P, Bellgard MI, Hampson DJ: Multilocus sequence typing as a tool for studying the molecular epidemiology and population structure of Brachyspira hyodysenteriae. Vet Microbiol 2009, 138:330–338.

17. Osorio J, Carvajal A, Naharro G, La T, Phillips ND, Rubio P, Hampson DJ: Dissemination of clonal groups of Brachyspira hyodysenteriae amongst pig farms in Spain, and their relationships to isolates from other countries. PLoS ONE 2012, 7(6):e39082.

18. Phillips ND, La T, Amin MM, Hampson DJ: Brachyspira intermedia strain diversity and relationships to the other indole-positive Brachyspira species. Vet Microbiol 2010, 143:246–254.

19. Råsbäck T, Johansson K-E, Jansson DS, Fellström C, Alikhani Y, La T, Dunn DS, Hampson DJ: Development of a multilocus sequence typing scheme for intestinal spirochaetes of the genus Brachyspira. Microbiology 2007, 153:4074–4087.

20. Neo E, La T, Phillips ND, Hampson DJ: Multiple locus variable number tandem repeat analysis (MLVA) of the pathogenic intestinal spirochaete Brachyspira pilosicoli. Vet Microbiol 2013, 163:299–344.

21. Humphrey SB, Stanton TB, Jensen NS, Zuerner RL: Purification and characterization of VSH-1, a gereralized tansducing bacteriophage of Serpulina hyodysenteriae. J Bacteriol 1997, 179:323–329.

22. Matson EG, Thompson MG, Humphrey SB, Zuerner RL, Stanton TB: Identification of genes of VSH-1, a prophage-like gene transfer agent of Brachyspira hyodysenteriae. J Bacteriol 2005, 187:5885–5892.

23. Motro Y, La T, Bellgard MI, Dunn DS, Phillips ND, Hampson DJ: Identification of genes associated with prophage-like gene transfer agents in the pathogenic intestinal spirochaetes Brachyspira hyodysenteriae, Brachyspira pilosicoli and Brachyspira intermedia. Vet Microbiol 2009, 134:340–345.

24. Mappley LJ, Black ML, Abuoun M, Darby AC, Woodward MJ, Darby A, Turner AK, Parkhill J, Bellgard MI, La T, La Ragione RM, Phillips ND, Hampson DJ: Comparative genomics of Brachyspira pilosicoli strains: genome rearrangements, reductions and correlation of genetic compliment with phenotypic diversity. BMC Genomics 2012, 13(1):454.

25. Lin C, den Bakker HC, Suzuki H, Lefébure T, Ponnala L, Sun Q, Stanhope MJ, Wiedmann M, Duhamel GE: Complete genome sequence of the porcine strain Brachyspira pilosicoli P43/6/78(T.). Genome Announc 2013, 1:1.

26. La T, Phillips ND, Hampson DJ: Development of a duplex PCR assay for the detection of Brachyspira hyodysenteriae and Brachyspira pilosicoli in pig feces. J Clin Microbiol 2003, 41:3372–3375.

27. Kunkle RA, Harris DL, Kinyon JM: Autoclaved liquid medium for propagation of Treponema hyodysenteriae. J Clin Microbiol 1986, 24:669–671.

28. Bellgard MI, Wanchanthuek P, La T, Ryan K, Moolhuijzen P, Albertyn Z, Shaban B, Motro Y, Dunn DS, Schibeci D, Hunter A, Barrero R, Phillips ND, Hampson DJ: Genome sequence of the pathogenic intestinal spirochete Brachyspira hyodysenteriae reveals adaptations to its lifestyle in the porcine large intestine. PLoS ONE 2009, 4(3):e4641.

29. Wanchanthuek P, Bellgard MI, La T, Ryan K, Moolhuijzen P, Chapman B, Black M, Schibeci D, Hunter A, Barrero R, Phillips ND, Hampson DJ: The complete genome sequence of the pathogenic intestinal spirochete Brachyspira pilosicoli and comparison with other Brachyspira genomes. PLoS ONE 2010, 5(7):e11455.

30. Håfström T, Jansson DS, Segerman B: Complete genome sequence of Brachyspira intermedia reveals unique genomic features in Brachyspira species and phage-mediated horizontal gene transfer. BMC Genomics 2011, 12:395.

31. Hall TA: Bioedit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser 1999, 41:95–98.

32. Jolley KA, Feil EJ, Chan MS, Maiden MC: Sequence type analysis and recombinational tests (START). Bioinformatics 2001, 17:1230–1231.

33. Nei M: F-statistics and analysis of gene diversity in subdivided populations. Annals Human Genet 1977, 41:225–233.

34. Haubold B, Hudson RR: LIAN 3.0: detecting linkage disequilibrium in multilocus data. Bioinfo Applic Notes 2000, 16:847–848.

35. Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 2007, 24:1596–1599.

36. Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.69. Seattle: Department of Genetics, University of Washington; 1993.