- About this Journal ·
- Abstracting and Indexing ·
- Advance Access ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
BioMed Research International
Volume 2013 (2013), Article ID 390354, 9 pages
Dynamic of Mutational Events in Variable Number Tandem Repeats of Escherichia coli O157:H7
1Laboratorio de Inmunoquímica y Biotecnología, Centro de Investigación Veterinaria de Tandil (CIVETAN), Facultad de Ciencias Veterinarias, UNCPBA, Campus Universitario, Paraje Arroyo Seco S/N, 7000 Tandil, Argentina
Received 18 February 2013; Accepted 22 July 2013
Academic Editor: Dietmar Quandt
Copyright © 2013 A. V. Bustamante 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.
VNTRs regions have been successfully used for bacterial subtyping; however, the hypervariability in VNTR loci is problematic when trying to predict the relationships among isolates. Since few studies have examined the mutation rate of these markers, our aim was to estimate mutation rates of VNTRs specific for verotoxigenic E. coli O157:H7. The knowledge of VNTR mutational rates and the factors affecting them would make MLVA more effective for epidemiological or microbial forensic investigations. For this purpose, we analyzed nine loci performing parallel, serial passage experiments (PSPEs) on 9 O157:H7 strains. The combined 9 PSPE population rates for the 8 mutating loci ranged from 4.4 × 10−05 to 1.8 × 10−03 mutations/generation, and the combined 8-loci mutation rate was of 2.5 × 10−03 mutations/generation. Mutations involved complete repeat units, with only one point mutation detected. A similar proportion between single and multiple repeat changes was detected. Of the 56 repeat mutations, 59% were insertions and 41% were deletions, and 72% of the mutation events corresponded to O157-10 locus. For alleles with up to 13 UR, a constant and low mutation rate was observed; meanwhile longer alleles were associated with higher and variable mutation rates. Our results are useful to interpret data from microevolution and population epidemiology studies and particularly point out that the inclusion or not of O157-10 locus or, alternatively, a differential weighting data according to the mutation rates of loci must be evaluated in relation with the objectives of the proposed study.
Repetitive DNA, which occurs in large quantities in eukaryotic cells, has been increasingly identified in prokaryotes . Repeats organized in tandem, representing a single locus and showing interindividual unit number variability, are designated variable-number tandem repeat (VNTR) loci . Repeat copy number variation at these loci is the consequence of mutations resulting in the gain or loss of a certain number of tandem repeats (TRs) and can lead to a very large number of alleles . These mutations are thought to occur via an intramolecular slipped-strand mispairing (SSM), which may occur in combination with inadequate DNA mismatch repair pathways . Recombination may also play a role, especially in mutations involving large numbers of repeat units , although SSM is thought to be the predominant mutational mechanism [5, 6].
The VNTR regions have been successfully used for subtyping purposes [1, 7], and multiple-locus VNTR analysis (MLVA) involves determination of the number of repeat copy units in multiple loci . Molecular epidemiology, the integration of molecular typing and conventional epidemiological studies, likely adds significant value to analyses of infections caused by pathogenic bacteria . MLVA has led to the development of a highly effective typing system for use in molecular epidemiology and forensic analyses [10, 11]. In recent emergent pathogens, TRs are a source of very informative markers for strain genotyping [12, 13]. However, any change in the rate of hypervariability in these TRs is problematic when trying to predict the phylogenetic relationships among isolates . Noller et al.  observed single locus variants during some outbreaks, and, therefore, one potential concern is that some VNTRs evolve so rapidly that multiple MLVA profiles would emerge during an outbreak initially caused by a single clone. In relation to it, Hyytiä-Trees et al.  proposed that during an outbreak two isolates differing in one repeat unit at one or two loci could be considered as related isolates. This is supported by the observation that isolates with five or more repeat differences at one or more loci were epidemiologically unrelated. Therefore, care has to be taken to select those repeat regions that show adequate stability over the timeframe in which the epidemiological investigations take place .
The necessity of more information on mutations and their rates in different loci to develop appropriate models to interpret MLVA data has been pointed out in several studies [8, 17–19]. Understanding of VNTR mutational rates and the factors affecting them would make MLVA more effective for epidemiological or microbial forensic investigations. Mutation rate data are valuable resources for probabilistic modeling of genetic relatedness, assessing the significance of finding genotype matches, near-matches, and mismatches [18, 20]. This information can be useful to determine the source of an outbreak, to asses true versus fortuitous disease clusters and to identify sites of exposure for human infection, among other purposes [20, 21].
In the present work our aim was to estimate mutation rates of Escherichia coli O157:H7 specific VNTR loci. For this purpose, we analyzed nine loci performing parallel, serial passage experiments (PSPEs) on nine verotoxigenic E. coli (VTEC) O157:H7 strains.
2. Materials and Methods
VTEC strains used in this study are listed in Table 1. The strains belong to O157:H7 serotype and were selected from those previously studied by Bustamante et al.  because each one of them represents a unique MLVA profile. Strains were stored at −80°C.
2.2. PSPE and Detection of Mutations
Parallel, serial passage experiments (PSPEs) have been previously described for Yersinia pestis  and for Escherichia coli O157:H7 . We performed similar PSPEs on 9 E. coli O157:H7 strains (Table 1), generating a set of in vitro populations representing ~211.500 total generations (23.5 generations/colony × 100 lineages × 10 passages × 9 PSPEs). In our work, each PSPE consisted of ~100 independent clonal lineages that were each serially passaged 10 times. For each PSPE, a single isolated colony of each strain () was used to start ~100 independent clonal lineages by streaking for single colonies on LB agar plates. All cultures were grown at 37°C for 24 h before the next passage. Each lineage was then serially passaged 10 times by streaking from a single colony from the previous passage. DNA was extracted from all ~100 lineages at (day 10) in the 9 PSPEs by lysis of a single colony. Mutational events for each strain were visualized by MLVA, as previously described by Bustamante et al. . Each lineage was tested by MLVA before and after all the passages. In order to know if the mutations were due to changes in repeat copy number or to point mutation events, the mutational products and parental alleles were sequenced. Mutational events and products for the nine strains were determined from an analysis of the populations.
2.3. Mutation Rate Calculations and Genetic Diversity Analysis after the PSPE
Single-locus and combined 9-loci mutation rates (μ) were estimated for each strain by dividing the observed number of mutations by the number of total generations (GT) . GT was calculated as the average number of generations per colony × the number of lineages with usable data () × the number of transfers. The number of generations per colony was determined to be 23.5 using an average of viable plate counts. In the combined single locus mutation rates, was based on an average of across all nine strains for each locus. To calculate the mutation rate in the combined 9 loci, was based on an average of across all 9 loci for each individual strain calculation and an average of across all 9 loci and all nine strains for the combined strain calculations.
In order to analyze the relationship between mutation rate and repeat number, progressive linear regression analyses (PROC REG, Statistical Analysis System, Version 9.2) were performed to detect the range of repeat numbers where the slope of the regression line was not significantly different from zero.
A dendrogram showing the genetic variability cumulated during the PSPE, taking into account each strain at with its derived mutated lineages, was constructed by UPGMA clustering using START version 1.0.5 .
3. Results and Discussion
Mutation is the major mechanism for generating diversity in clonal organisms, and VNTR loci are ones of the fastest mutation sequences that allow detecting genetic differences. In this study, we report the mutation dynamics of nine VNTR loci used for VTEC O157 subtyping. Eight of the nine VNTRs of this MLVA protocol have repeat units (RU) of six or seven bp and the remaining one (Vhec2), 18 bp. The genomic regions susceptible to intramolecular slipped-strand mispairing (SSM) are those that contain a short, contiguous homogenous or heterogeneous repetitive DNA sequence of 6 bp or less . Repeat motifs longer than eight base pairs show typically less unit number variation but are more likely to have point mutations . In this study, from the nine VNTRs analyzed, eight mutated across all PSPEs; meanwhile the Vhec2 locus did not exhibit changes in repeat number, in accordance with data from other authors [14, 17, 20]. Noller et al. , using a different experimental model, analyzed seven of the 9 VNTRs studied here but they only observed mutation in three VNTRs (Vhec4, O157-3, and O157-10).
In order to detect rare occurrence of mutation before , in some PSPEs (5 out 9) DNA was also extracted at (day 5); however, we did not detect any mutations at this earlier time point.
In the nine in vitro PSPE populations, we observed 57 mutational events corresponding to a combined 8-loci mutation rate of 2.5 × 10−03 mutations/generation, and noticeably, 41 (72%) of the mutation events corresponded to locus O157-10 (Table 2). Comparison between mutated and parental allele sequences confirmed that mutations were changes involving complete repeat units, with only one exception at a single locus in one PSPE (Vhec4, FC O157) in which the change was due to a point mutation (Table 2). Of the 56 repeat mutations, 33 (59%) were insertions and 23 (41%) were deletions (Table 2). This proportion of insertion/deletion of repeats is similar to that observed by Vogler et al. , 62% insertions and 38% deletions, respectively. They explained the insertion bias as due to the influence of a single locus, O157-10, in a single PSPE, but in our study, on the contrary, when O157-10 was excluded from the analysis the proportion of insertions increased to 67%.
Several single and multiple repeat events could be observed (Figures 1 and 2), 29 (52%) mutation events were changes due to single repeats (either insertion or deletion), and the remaining 27 (48%) involved multiple repeat changes (Table 2). This similar proportion between single and multiple repeat changes differs with one of the general properties of the mutation model for bacterial VNTRs proposed by Vogler et al. , who indicated that the majority of mutations involve single repeats, although a certain percentage of events consist of multiple repeats. In agreement with the results of Vogler et al. , the frequency of mutations followed a decreasing pattern from 1 to 3 RU (or to 5 RU if only the deletions are taken into account), and the frequency of deletions of more than 5 RU did not seem to follow any pattern (Figure 2). In the multiple repeat events the changes involved up to 23 RU, and 14 of the 27 multiple repeat mutations were “large repeat copy number” events. This last term was defined by Vogler et al.  in order to name mutational events involving greater than four repeat units. The occurrence of “large repeat copy number” events observed in this study (14/27 = 52%) was markedly higher than the one obtained by Vogler et al.  (12/47 = 26%). Interestingly, all “large repeat copy number” events detected by us were at locus O157-10, with insertions and deletions occurring at nearly the same frequency. Noticeably, two of the deletions comprised the complete VNTR; meanwhile in all the insertion events, 5 RU were inserted.
Long microsatellites are likely to mutate to shorter ones, and short microsatellites are likely to mutate to longer ones . Considering the MLVA profile of the different strains analyzed, the largest repeat copy numbers were present at locus O157-10, with 7 alleles comprising more than 22 RU, while in the remaining analyzed loci the largest allele was of 15 RU. VNTR alleles of very large size could be detrimental or just highly unstable, leading to the observation of a greater number of deletion events , but in our study, although we observed that the VNTR O157-10 was the most unstable, we did not observe a trend towards the deletions (23 insertions versus 18 deletions).
The combined 9 PSPE population rates for the 8 mutating loci ranged from 4.4 × 10−05 to 1.8 × 10−03 mutations/generation (Table 2). Our experimental design allowed us to detect mutation rates equal or higher than 4.4 × 10−05 mutations/generation. No mutation was observed at Vhec2 suggesting that its mutation rate is near to or less than that value, which is the detection limit given the number of generations in the PSPE population.
The mutation rates at locus O157-10 ranged from 1.3 × 10−04 to 4.8 × 10−04 mutations/generation among the nine PSPE strains (Table 2), and this locus also presented the highest mutation rate when the results of all the strains were taken into account (1.8 × 10−03 mutations/generation) (Table 2). Forty-one of the 57 (72%) observed mutations occurred at this locus, and when mutations associated to this locus were removed from the combined data set, the mutation rate (excluding O157-10) decreased noticeably (6.9 × 10−04 mutations/generation). Thus, our results, as those from Vogler et al. , showed that allele sizes at locus O157-10 were the most variable and therefore impacted strongly on the combined locus mutation rates.
It has been described in different bacterial species that mutation rates correlate with the number of repeats present in the VNTRs [18, 20, 28]. Taking into account the present results, the distribution of the mutation rates versus the number of RU present in the corresponding strain at is shown in Figure 3. It can be observed that for alleles with up to 13 UR a more or less constant and low mutation rate corresponds to them; meanwhile longer alleles are associated with higher and variable mutation rates (values higher than 1 × 10−04). Furthermore, multiple repeat changes predominated among this second group of alleles, which correspond mainly to the locus O157-10. When we compared our results with those obtained by other authors who analysed VNTRs from VTEC O157:H7, Yersinia pestis, and Vibrio parahaemolyticus [18, 20, 28], we noticed that the distribution was similar, although they did not discuss this and instead focused on the correlation between mutation rates and repeat copy numbers. However, as O157-10 was the only locus which presented large alleles, its influence cannot be ruled out, and therefore we cannot make generalized conclusions about allele size and mutation rate.
In the present study, when results were compared among strains, strain EDL 933 noticeably showed mutations in 6 VNTRs, followed by Mat167/6 and FC O157 which presented variation in only 3 loci. On the other hand, two strains showed repeat number variation in only one locus (O157-10). Mutation rates were not so different among strains, but when they were recalculated without the inclusion of O157-10 locus, EDL 933 strain showed the highest rate followed by FC O157. On the contrary, Vogler et al.  did not observe a pure “strain” effect on mutation rate and indicated that instead allele size at individual loci did affect the mutation rate. Furthermore, the MLVA profile reported for the EDL 933 strain differed among several studies [8, 22, 23, 29] and also with an in silico analysis with the GenBank sequence (accession no. AE005174). These differences are probably due to repeated propagation of the strain and could also reflect a high mutation degree of the EDL 933 strain.
MLVA is applied to study the relationships among different isolates and as an aid to determine the presence of an outbreak. The analysis of our in vitro PSPE population showed that each parental strain and its derived lineages became clustered evidencing their clonal relationship. The dendrogram showing the original MLVA profiles (at ) of the strains and those of the mutated lineages is shown in Figure 4.
Cooley et al.  studied mutations in stressed bacteria and concluded that conditions such as elevated temperature and starvation are associated with mutations. They observed that tandem repeats were stable after limited replication, but in our study the strains were not subjected to selection pressure, and, however, we observed mutations in several VNTRs. Interestingly, their results showed that one locus (O157-10) mutated in all the conditions tested. Several studies recognized O157-10 as the most polymorphic locus [15, 19, 24, 30–33]. Furthermore, when the Pulse Net USA subtyping network recently established a standardized protocol in order to characterize verotoxigenic Escherichia coli O157 by MLVA, this locus was excluded because it showed too much variability during outbreaks, generating data that confounded epidemiological investigations . However, Keys et al.  previously suggested that markers with high diversities are crucial in discriminating between closely related isolates, and they highlighted the importance of including O157-10, as an example of this kind of crucial markers. In relation with it, Urdahl et al.  observed that removing this locus from the analysis resulted in an MLVA that did not discriminate between closely related strains.
In summary, taking into account different points of view of several authors and our own results, care must be taken to select loci in order to be included in MLVA and to interpret data obtained from highly discriminating loci. Particularly, some aspects such as the inclusion or not of O157-10 locus in the MLVA or alternatively, a differential weighting data according to the mutation rates, have to be considered in relation with the objectives of a certain study.
Our set of in vitro populations showed that mutation rates detected for nine VNTR loci used for VTEC O157 subtyping were according to those reported previously for other bacterial VNTRs, with a similar proportion between single and multiple repeat changes. A markedly high proportion of “large repeat copy number” events was observed in this study.
Long alleles were associated with high and variable mutation rates and corresponded to locus O157-10, which presented an outstanding high mutation rate.
Our results are useful to interpret data from microevolution and population epidemiology studies and suggest, together with previously published ones, that a differential weighting of VNTR data should be applied according to the available mutation rates of loci in order to get more accurate phylogenetic relationships.
Conflict of Interests
The authors declare that there is no conflict of interests.
The authors thank María Rosa Ortiz for technical assistance and Juan Passucci and Edgardo Rodríguez for statistical analysis. This work was supported by grants from the CONICET, FONCYT, and SECAT-UNICEN.
- A. Van Belkum, S. Scherer, L. Van Alphen, and H. Verbrugh, “Short-sequence DNA repeats in prokaryotic genomes,” Microbiology and Molecular Biology Reviews, vol. 62, no. 2, pp. 275–293, 1998.
- Y. Nakamura, M. Leppert, and P. O'Connell, “Variable number of tandem repeat (VNTR) markers for human gene mapping,” Science, vol. 235, no. 4796, pp. 1616–1622, 1987.
- S. Strand, T. A. Prolla, R. M. Liskay, and T. D. Petes, “Destabilization of tracts of simple repetitive DNA in yeast by mutations affecting DNA mismatch repair,” Nature, vol. 365, pp. 274–276, 1993.
- G. Levinson and G. A. Gutman, “Slipped-strand mispairing: a major mechanism for DNA sequence evolution,” Molecular Biology and Evolution, vol. 4, no. 3, pp. 203–221, 1987.
- H. Ellegren, “Microsatellite mutations in the germline: implications for evolutionary inference,” Trends in Genetics, vol. 16, no. 12, pp. 551–558, 2000.
- C. Schlötterer, “Evolutionary dynamics of microsatellite DNA,” Chromosoma, vol. 109, pp. 365–371, 2000.
- Y. Kashi, D. King, and M. Soller, “Simple sequence repeats as a source of quantitative genetic variation,” Trends in Genetics, vol. 13, no. 2, pp. 74–78, 1997.
- E. Hyytiä-Trees, S. C. Smole, P. A. Fields, B. Swaminathan, and E. M. Ribot, “Second generation subtyping: a proposed PulseNet protocol for multiple-locus variable-number tandem repeat analysis of Shiga toxin-producing Escherichia coli O157 (STEC O157),” Foodborne Pathogens and Disease, vol. 3, no. 1, pp. 118–131, 2006.
- F. Denœud and G. Vergnaud, “Identification of polymorphic tandem repeats by direct comparison of genome sequence from different bacterial strains: a web-based resource,” BMC Bioinformatics, vol. 5, article 4, 2004.
- J. M. Doll, P. S. Zeitz, P. Ettestad, A. L. Bucholtz, T. Davis, and K. Gage, “Cat-transmitted fatal pneumonic plague in a person who traveled from Colorado to Arizona,” American Journal of Tropical Medicine and Hygiene, vol. 51, no. 1, pp. 109–114, 1994.
- A. M. Klevytska, L. B. Price, J. M. Schupp, P. L. Worsham, J. Wong, and P. Keim, “Identification and characterization of variable-number tandem repeats in the Yersinia pestis genome,” Journal of Clinical Microbiology, vol. 39, no. 9, pp. 3179–3185, 2001.
- A. Van Belkum, S. Scherer, W. Van Leeuwen, D. Willemse, L. Van Alphen, and H. Verbrugh, “Variable number of tandem repeats in clinical strains of Haemophilus influenzae,” Infection and Immunity, vol. 65, no. 12, pp. 5017–5027, 1997.
- P. Le Flèche, M. Fabre, F. Denoeud, J.-L. Koeck, and G. Vergnaud, “High resolution, on-line identification of strains from the Mycobacterium tuberculosis complex based on tandem repeat typing,” BMC Microbiology, vol. 2, no. 1, article 37, 2002.
- M. B. Cooley, D. Carychao, K. Nguyen, L. Whitehand, and R. Mandrell, “Effects of environmental stress on stability of tandem repeats in Escherichia coli O157:H7,” Applied and Environmental Microbiology, vol. 76, no. 10, pp. 3398–3400, 2010.
- A. C. Noller, M. C. McEllistrem, A. G. F. Pacheco, D. J. Boxrud, and L. H. Harrison, “Multilocus variable-number tandem repeat analysis distinguishes outbreak and sporadic Escherichia coli O157:H7 isolates,” Journal of Clinical Microbiology, vol. 41, pp. 5389–5397, 2003.
- A. Van Belkum, “Tracing isolates of bacterial species by multilocus variable number of tandem repeat analysis (MLVA),” FEMS Immunology and Medical Microbiology, vol. 49, no. 1, pp. 22–27, 2007.
- A. C. Noller, M. C. McEllistrem, K. A. Shutt, and L. H. Harrison, “Locus-specific mutational events in a multilocus variable-number tandem repeat analysis of Escherichia coli O157:H7,” Journal of Clinical Microbiology, vol. 44, no. 2, pp. 374–377, 2006.
- A. J. Vogler, C. E. Keys, C. Allender et al., “Mutations, mutation rates, and evolution at the hypervariable VNTR loci of Yersinia pestis,” Mutation Research, vol. 616, no. 1-2, pp. 145–158, 2007.
- A. M. Urdahl, N. J. C. Strachan, Y. Wasteson, M. MacRae, and I. D. Ogden, “Diversity of Escherichia coli O157 in a longitudinal farm study using multiple-locus variable-number tandem-repeat analysis,” Journal of Applied Microbiology, vol. 105, no. 5, pp. 1344–1353, 2008.
- A. J. Vogler, C. Keys, Y. Nemoto, R. E. Colman, Z. Jay, and P. Keim, “Effect of repeat copy number on variable-number tandem repeat mutations in Escherichia coli O157:H7,” Journal of Bacteriology, vol. 188, no. 12, pp. 4253–4263, 2006.
- J. L. Lowell, D. M. Wagner, B. Atshabar et al., “Identifying sources of human exposure to plague,” Journal of Clinical Microbiology, vol. 43, pp. 650–656, 2005.
- B.-A. Lindstedt, E. Heir, E. Gjernes, T. Vardund, and G. Kapperud, “DNA fingerprinting of shiga-toxin producing Escherichia coli O157 based on multiple-locus variable-number tandem-repeats analysis (MLVA),” Annals of Clinical Microbiology and Antimicrobials, vol. 2, article 12, 2003.
- C. Keys, S. Kemper, and P. Keim, “Highly diverse variable number tandem repeat loci in the E. coli O157:H7 and O55:H7 genomes for high-resolution molecular typing,” Journal of Applied Microbiology, vol. 98, no. 4, pp. 928–940, 2005.
- A. V. Bustamante, P. M. A. Lucchesi, and A. E. Parma, “Molecular characterization of verocytotoxigenic Escherichia coli O157:H7 isolates from Argentina by multiple-loci VNTR analysis (MLVA),” Brazilian Journal of Microbiology, vol. 40, no. 4, pp. 927–932, 2009.
- J. M. Girard, D. M. Wagner, A. J. Vogler et al., “Differential plague-transmission dynamics determine Yersinia pestis population genetic structure on local, regional, and global scales,” Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 22, pp. 8408–8413, 2004.
- K. A. Jolley, E. J. Feil, M.-S. Chan, and M. C. J. Maiden, “Sequence type analysis and recombinational tests (START),” Bioinformatics, vol. 17, no. 12, pp. 1230–1231, 2002.
- Y. Lai and F. Sun, “The relationship between microsatellite slippage mutation rate and the number of repeat units,” Molecular Biology and Evolution, vol. 20, no. 12, pp. 2123–2131, 2003.
- K. García, R. G. Gavilán, M. G. Höfle, J. Martínez-Urtaza, and R. T. Espejo, “Microevolution of pandemic Vibrio parahaemolyticus assessed by the number of repeat units in short sequence tandem repeat regions,” PLoS One, vol. 7, no. 1, Article ID e30823, 2012.
- L. Yun, Y. Gu, L. Zha, F. Zhu, and Y. Hou, “Utility of multilocus variable number tandem repeat analysis as a microbial forensic tool for subtyping Chinese Escherichia coli O157: H7 strains,” Forensic Science International, vol. 3, no. 1, supplement 3, pp. e293–e294, 2011.
- B.-A. Lindstedt, T. Vardund, and G. Kapperud, “Multiple-locus variable- number tandem-repeats analysis of Escherichia coli O157 using PCR multiplexing and multi-colored capillary electrophoresis,” Journal of Microbiological Methods, vol. 58, no. 2, pp. 213–222, 2004.
- F. Kawamori, M. Hiroi, T. Harada et al., “Molecular typing of Japanese Escherichia coli O157:H7 isolates from clinical specimens by multilocus variable-number tandem repeat analysis and PFGE,” Journal of Medical Microbiology, vol. 57, no. 1, pp. 58–63, 2008.
- M. Murphy, D. Minihan, J. F. Buckley, M. O'mahony, P. Whyte, and S. Fanning, “Multiple-Locus Variable number of tandem repeat Analysis (MLVA) of Irish verocytotoxigenic Escherichia coli O157 from feedlot cattle: uncovering strain dissemination routes,” BMC Veterinary Research, vol. 4, article 2, 2008.
- E. Hyytiä-Trees, P. Lafon, P. Vauterin, and E. M. Ribot, “Multilaboratory validation study of standardized multiple-locus variable-number tandem repeat analysis protocol for Shiga toxin-producing Escherichia coli O157: a novel approach to normalize fragment size data between capillary electrophoresis platforms,” Foodborne Pathogens and Disease, vol. 7, no. 2, pp. 129–136, 2010.