|
|
||||||||
1 Otto-von-Guericke-University, Institute of Medical Microbiology, Magdeburg, Germany
2 Robert Koch Institute, Wernigerode, Germany
3 West Midlands Public Health Laboratory, Health Protection Agency, Birmingham, UK
Correspondence
B. Ghebremedhin
beniam.ghebremedhin{at}medizin.uni-magdeburg.de
Received 10 August 2006
Accepted 20 October 2006
Abbreviations: CIP, ciprofloxacin; CLI, clindamycin; EMRSA, epidemic meticillin-resistant Staphylococcus aureus; ERY, erythromycin; fAFLP, fluorescent amplified fragment length polymorphism; FAM, 6-carboxyfluorescein; FUS, fusidic acid; GEN, gentamicin; LIN, linezolid; MLST, multilocus sequence typing; MRSA, meticillin-resistant Staphylococcus aureus; OXA, oxacillin; PEN, penicillin; PVL, PantonValentine leukocidin; RCCF, rehabilitation and chronic care facility; RIF, rifampicin; SCCmec, staphylococcal cassette chromosome mec; SIRU, staphylococcal interspersed repeat unit; SXT, trimethoprim/sulfamethoxazole; TEI, teicoplanin; TET, tetracycline; VAN, vancomycin; VNTR, variable-number tandem repeat.
| INTRODUCTION |
|---|
|
|
|---|
In a previous study, we analysed 199 MRSA isolates from patients at University Hospital in Magdeburg, Germany, and the rehabilitation and chronic care facility (RCCF) in Flechtingen, Germany (Ghebremedhin et al., 2005) and observed the emergence of the Barnim EMRSA ST22 strain. An understanding of the genetic microevolution and clonal spread of strains is necessary to understand the transmission of MRSA and to control its spread. Phenotypic typing methods are mostly restricted to the resistance pattern of the respective MRSA strains (Berger-Bächi, 2002). At present, there are numerous typing methods available, but there is no consensus regarding which method is the best to apply.
For epidemiological typing and for studying the dynamics of MRSA evolution in a hospital environment, knowledge of both strain identity and variability within the strain allow us to make judgements on whether direct cross-infection or independent acquisition has taken place. Thus, we need slowly [e.g. multilocus sequence typing (MLST), agr, SCCmec typing] as well as more rapidly [e.g. PFGE, fluorescent amplified fragment length polymorphisms (fAFLP), staphylococcal interspersed repeat units (SIRUs), spa typing] evolving genetic markers to study hospital-acquired MRSA. Methods that index neutral variation, which accumulates slowly in the genome, are better suited to answering global epidemiological questions.
PFGE is a highly discriminatory typing method that is well suited to the investigation of microevolution and recent transmission within a hospital and within a narrow time schedule. However, it has several limitations, including ease of reproducibility and interpretation of band patterns, even when using standardized methods (Murchan et al., 2003). Guidelines and interpretative criteria for the analysis of discrete sets of isolates obtained during outbreaks have been proposed, but these are not appropriate for the analysis of large populations of organisms collected over extended periods. Somehow, PFGE is too discriminatory for studies assessing long-term epidemiology, as factors that may dramatically alter macrorestriction banding profiles (such as insertions, deletions and inversions) may obscure similarity among epidemiologically related isolates (Tenover et al., 1995).
MLST is carried out by sequencing an internal fragment of seven unlinked housekeeping genes. Allelic profiles can be determined, and a clone is defined as a group of isolates having a strictly identical sequence for all seven genes (Enright et al., 2000). MLST examines the slowly evolving genomic core and is a useful tool for defining chromosomal ancestry (Enright et al., 2000). Although MLST is useful for defining evolutionary events, it does not provide a good marker for microevolutionary events within a hospital (Urwin & Maiden, 2003).
Variable-number tandem repeats (VNTRs) have been used as markers for the strain typing of various bacteria (Francois et al., 2005; Hawkey et al., 2003; Onteniente et al., 2003; van Belkum, 1999). Determination of the number of repeats at each locus produces a digital profile, providing a highly portable typing method allowing comparisons between laboratories. Hardy et al. (2006) previously described the presence of VNTRs in staphylococci. Seven VNTRs, termed SIRUs, distributed around the genome are described, occurring in both unique and multiple sites and varying in length from 48 to 159 bp.
fAFLP is a DNA fingerprinting method based on restriction cutting of DNA and stringent PCR amplification of the resulting fragments. fAFLP analysis represents a highly discriminatory and reproducible tool for subtyping genetically homogeneous genomes and identifying outbreak genotypes within bacterial genera (Desai et al., 1998; Hookey et al., 1999; Nair et al., 2000; Thyssen et al., 2000).
The agr locus regulates the production of most staphylococcal exoproteins, including exoenzymes, toxins, surface proteins and other virulence factors, by means of a density-dependent autoinducible signal transduction system driven by a short, post-translationally processed peptide (Ji et al., 1995, 1997).
The mecA gene encoding the alternative penicillin-binding protein (PBP2a) is found on the chromosome in a genomic island designated staphylococcal cassette chromosome mec (SCCmec) (Katayama et al., 2000). At least five major varieties of SCCmec exist and are defined by the specific mec locus present (composed of mecA and its two regulatory genes, mecI and mecRI), as well as the variety of site-specific recombinase genes present (either ccrAB or ccrC) (Katayama et al., 2000; Ma et al., 2002). SCCmec types IIII are larger elements (3467 kb) that tend to contain resistance determinants in addition to mecA and are found more frequently in hospital-acquired MRSA strains (Katayama et al., 2000; Ito et al., 2004). In contrast, our ST22-MRSA strains harboured SCCmec type IV, which can be subtyped again into IVad (Okuma et al., 2002; Oliveira & de Lencastre, 2002), now increasingly recognized as causative agents of community-associated disease in patients that are not affiliated with a hospital or healthcare environment (Zetola et al., 2005).
In the present study, we conducted an analysis of ST22-MRSA-IV isolates to evaluate and compare widely used molecular techniques, namely, MLST, SCCmec typing and subtyping, fAFLP and PFGE, for their suitability and usefulness in the typing and fingerprinting of bacterial isolates and in tracking clonal evolution.
| METHODS |
|---|
|
|
|---|
DNA extraction. Strains were grown on brain heart infusion agar or in brain heart infusion broth at 37 °C overnight. Genomic DNA used as target for all molecular methods was extracted using a Qiagen DNA extraction kit according to the manufacturer's suggestions with the modification that 30 µg lysostaphin ml1 and 1 mg lysozyme ml1 were added at the cell-lysis step. The concentration of DNA was estimated spectrophotometrically (Ghebremedhin et al., 2005).
MLST. MLST was carried out as described by Enright et al. (2000). Analysis of the allelic profile of S. aureus isolates was performed using the database of the MLST home web page (http://www.mlst.net), where seven numbers depicting the allelic profile were assigned that defined the MLST type.
PFGE. All isolates were subjected to PFGE analysis of SmaI-digested chromosomal DNA, which was performed as described by Witte et al. (2004). Fragments of sizes in the range of the ladder (501000 kb) were included in the analysis. The criteria for PFGE pattern analysis state that up to three PFGE fragment differences compared with the prototype pattern B1 indicate a close relationship and thus that the isolate is probably a subclone of the prototype (Tenover et al., 1995).
PCR for SCCmec IV subtyping. Analysis of SCCmec type IV was performed as described by Oliveira & de Lencastre (2002). For differentiation between the SCCmec types IVa and IVb, we used the protocol of Okuma et al. (2002).
agr group-specific multiplex PCR, and PantonValentine leukocidin (PVL) gene and toxin gene detection. Extracted genomic DNA was used as a template to amplify specific agr alleles (GenBank accession nos X52543, AF001782, AF001783, AF288215, Z49220, AF346724 and AF346725). For multiplex PCR, one primer set was prepared to amplify the four specific S. aureus agr alleles using the primers described by Lina et al. (2003). Details are given in Ghebremedhin et al. (2005). The primers agrB and agrC were used to amplify the variable part of the agr operon (1060 bp) with primer annealing conducted at 50 °C (Strommenger et al., 2004).
The genes seasee, segsej, tsst-1, eta, etb, hlgA, hlgC/B, lukE/D and lukS-lukF-PV (lukPV) were detected by PCR, as described previously (Ghebremedhin et al., 2005).
spa gene typing. The procedures described by Harmsen et al. (2003) were used for sequencing and analysis of the spa gene of all of the MRSA isolates.
fAFLP analysis. fAFLP analysis was performed as described by van den Braak et al. (2004). Briefly, 200 ng chromosomal DNA was digested with EcoRI and MseI (both from New England Biolabs) at 37 °C. The DNA fragments were ligated using double-stranded EcoRI and MseI adapters for 2 h at 37 °C. A pre-selective PCR with adapter-specific primers EcoRI (5'-GACTGCGTACCAGCTT-3') and MseI (5'-GATGAGTCCTGATCGC-3') was then performed. An aliquot was subjected to a selective PCR in which one of the primers contained additional nucleotides (primer MseI+CC or primer MseI+AC) at the 3' end, whilst primer EcoRI was labelled with 6-carboxyfluorescein (FAM) at the 5' end. PCR amplifications were performed with a Perkin-Elmer model 2400 thermocycler using the following protocol: an initial step of 72 °C for 2 min, 20 cycles of 94 °C for 1 s, 56 °C for 30 s and 72 °C for 2 min, and a final extension step at 72 °C for 7 min. PCR products were stored at 20 °C until analysis of the fAFLP using an automated gel electrophoresis and detection system.
fAFLP fragments were analysed on an ABI Prism 3100 DNA sequencer by adding 2 µl of the PCR product to 24 µl deionized formamide. To assess the size of the obtained fragments accurately, 1 µl GeneScan-500 marker (Applied Biosystems) labelled with the red fluorescent dye 6-carboxy-x-rhodamine was used as internal size standards.
The raw data from the GeneScan were converted to .ab1 format for analysis using BioNumerics software (version 3.0; Applied Maths). Cluster analysis of the fAFLP patterns was performed with the Dice similarity coefficient using optimization at 1 % and position tolerance and change towards the end of the fingerprint at 1 %. Dendrograms were calculated using the unweighted pair-group method with arithmetic mean (UPGMA).
Strains were analysed individually by calculating the number of isolates that clustered together with the number of unique isolates to give a crude measure of the heterogeneity of strain populations.
Tandem repeat locus identification (SIRU typing). SIRU typing of the isolates was conducted according to the method of Hardy et al. (2004) and interpretation of the results was performed according to Hardy et al. (2004, 2006).
Statistical analysis. The Hunter (1990) discriminatory index was used to calculate the level of discrimination of each typing method.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
SCCmec typing and subtyping
To study additional differences among the isolates according to SCCmec IV, subtype analysis was performed by PCR. SCCmec subtypes could be determined for all 42 ST22-MRSA strains by PCR and were found to be type IVa.
spa typing
The spa gene of S. aureus encodes protein A and was used for typing of the 42 ST22-MRSA-IV isolates using Ridom StaphType software. The majority of the ST22 isolates were typed as spa t032 (n=39; 92.9 %) with the repeat succession 26-23-23-13-23-31-29-17-31-29-17-25-17-25-16-28, two isolates were typed as spa t474 (with the repeat succession 26-23-13-23-31-05-17-25-16-28) and one isolate from the RCCF in Flechtingen was typed as spa t022 with the repeat succession 26-23-13-23-31-29-17-31-29-17-25-17-25-16-28 (Table 1
). The spa type t022 had one deletion in the repeat succession compared with the prototype t032, whereas t474 differed from t032 by one insertion and seven deletions. spa type t474 was reflected by the PFGE profile B2, whilst spa type t022 belonged to PFGE prototype B1.
spa typing has been shown to be discriminatory enough for outbreak investigation but also to be efficient for global epidemiological studies. The isolates that belonged to spa type t032 revealed a heterogeneous resistance profile, whilst the two isolates with spa type t474 shared the same resistance phenotype (antibiotype III) (see Table 1
). From this result, we concluded that the ST22-MRSA clone may indeed be heterogeneous. Our data support the previous findings that sequence information applied as MLST or spa types is useful for both evolutionary studies and global epidemiological analysis, and to a lesser extent for short-term or local epidemiological analyses (Tang et al., 2000). However, spa typing is more discriminatory compared with MLST data.
PFGE pattern analysis
PGFE analysis of the 42 ST22-MRSA strains revealed two different restriction patterns, A and B (distinguishable by four to five fragment bands), with most of the isolates exhibiting pattern B (n=37; 88.1 %) (Fig. 1
). Within pattern A (n=4), we detected no variations. Within the major pattern B, minor variations (up to three bands) allowed their differentiation into the prototype pattern B1 and the variant patterns B2, B3, B5, B6 and B7 according to O'Neill et al. (2001). One strain showed a pattern variation not described by O'Neill et al. (2001), but fitted PFGE pattern P10 as described by Grady et al. (1999). This new pattern type differed by four bands from pattern B1. The prototype B1 was isolated over the time of the study from different departments at University Hospital, whilst A1 occurred exclusively in the Department of Dermatology between February and March 2003. The profile B2 (n=2) was isolated from two different patients within 2 days in the Anaesthesiology ICU. Strains exhibiting the profile B3 (n=2) were detected at different time points in two specimens that originated from the RCCF in Flechtingen. The additional profiles B5 (n=1), B6 (n=1) and B7 (n=1) were isolated from different departments at University Hospital only once during the study period.
|
SIRU typing and correlation with PFGE patterns
PCR was performed on the ST22-MRSA isolates to analyse the discriminatory power of the different loci, as described in Methods. Among the 42 isolates, 6 different SIRU profiles (1303h-2, 1401h-2, 1403b-2, 1403g-2, 1403h-2 and 1603h-2) were detected (Tables 1 and 2![]()
). The SIRU profile 1403h-2 was the predominant type (n=32; 76.2 %). Of these 32 strains, 27 belonged to PFGE profile B1, 2 strains belonged to PFGE profile B3 and 1 strain to PFGE profile A1. The B7 strain and the unique (P10) strain also belonged to SIRU type 1403h-2. The SIRU type 1603h-2 (n=1) correlated with PFGE profile B5 (n=1), SIRU type 1403b-2 (n=2) was associated with both strains of PFGE profile B2 exhibiting the spa type t474, and SIRU type 1401h-2 (n=1) correlated with PFGE profile B6 (n=1). SIRU type 1403g-2 (n=1) was found to belong to the pulsotype B1 strain but with spa type t022. Surprisingly, SIRU profile 1303h-2 (n=5) was found within the two major PFGE restriction patterns, B1 (n=2) and A1 (n=3) (Table 2
).
|
Cluster analysis of fAFLP data
fAFLP was used to subtype isolates of the epidemic ST22-MRSA. The optimum combination of enzymes and primers for analysis of ST22-MRSA had been predicted by modelling fAFLP using different primer combinations. The primer combinations that were tested initially (EcoRI+A and MseI, EcoRI+C and MseI, EcoRI+AT and MseI, EcoRI and MseI+C, EcoRI and MseI+AC, and EcoRI and MseI+CC) generated fragments ranging in size from 80 to 400 bp. We decided to conduct fAFLP using the selective primers MseI+AC or MseI+CC and FAM-labelled EcoRI due to the discriminatory power revealed by these primer combinations. The presence and absence of fragments in each profile were scored. For reasons of reproducibility, the DNA was reextracted and the samples were again subjected to fAFLP. The fragment profiles from the different DNA extracts of the same isolate were shown to be nearly identical (data not shown).
According to the fAFLP MseI+AC and FAM-labelled EcoRI fingerprints, analysis of the banding patterns with 1 % optimization and 1 % position tolerance revealed three major clusters and one outlier cluster (Fig. 2
). Cluster I (n=12) contained the major relatively homogeneous group of isolates and consisted exclusively of proto-pulsotypes. At 95 % similarity, the cluster could be further subdivided into three subclusters, Ia (n=10), Ib (n=1) and Ic (n=1). Cluster II (n=12) contained the type strain of the prototype B1 (n=8), and A1 (n=1), P10 (n=1) and B2 (n=2) strains, which were obtained from patients from different departments and in different isolation periods. The B2 strains that originated from the Anaesthesiology ICU also grouped in this cluster.
|
Cluster IV (n=1) included only one B1 strain. Cluster IV joined the other major clusters at 75.7 % similarity. Thus, A1 strains belonged to two main clusters. In this regard, three isolates of A1 were grouped into cluster III and the remaining one into cluster II.
Using the primer combination MseI+CC and FAM-labelled EcoRI in fAFLP fingerprints, and using 1 % optimization and 1 % position tolerance, three major clusters and one outlier were detected (Fig. 3
). Cluster I (n=15) consisted of members of the proto-pulsotype B1 (n=9), A1 (n=3), B2 (n=1), B7 (n=1) and unique P10 (n=1), whereas the other isolates of B2 grouped into cluster III (n=6), which also included the proto-pulsotype B1 isolates (n=5). Cluster II (n=20) also included the proto-pulsotype B1 isolates (n=16), B3 (n=2), B5 (n=1) and B6 (n=1). Thus, the two isolates that were typed as spa type t474 with SIRU profile 1403b-2 did not cluster in one group. Within the different clusters, further subdivisions could be performed for the isolates of cluster II with a similarity of 88.2 %. The outlier isolate joined the other clusters with a similarity of 64.7 %.
|
Thus, fAFLP analysis can be expected to yield insights into the evolution of MRSA as a species. The fAFLP technique samples the whole genome sequence in a predictable and rigorous fashion by monitoring base substitutions accumulating throughout the genome, rather than being based on the fingerprints of a mobile genetic element. This could provide a valuable measure of microevolutionary change, e.g. within one clone such as ST22-MRSA-IV. The discriminatory power of fAFLP is determined by the choice of restriction enzymes and by the selectivity of the primers used. In case of the two isolates of pulsotype B2, the fAFLP MseI+CC and EcoRI fingerprints distinguished between these isolates, whereas the fAFLP MseI+AC and EcoRI pattern banding analysis put both isolates of pulsotype B2 in the same cluster II, which was supported by the SIRU profile pattern and spa typing. In summary, genome-sequence-derived fAFLP was broadly congruent with PFGE for typing MRSA. Its resolving power generally appeared to be superior to the other techniques used in this study (Table 3
). This is particularly the case where strains contain few genetic discriminatory elements and where fAFLP might be capable of revealing chains of transmission not demonstrable by PFGE.
|
Classical EMRSA-15 is characterized by the production of enterotoxin C (sec), although enterotoxin C-negative variants of EMRSA-15 have been reported (O'Neill et al., 2001). In our study, all ST22-MRSA isolates carried the gene for enterotoxin C and showed a homogeneous profile of toxin genes. We observed acquisition of the sea gene in four MRSA strains and acquisition of the PVL gene in five cases. It is known that enterotoxin genes are located on mobile elements. In this regard, the staphylococcal enterotoxin A gene is encoded by phage (Betley & Mekalanos, 1985). The additional presence of the sea or PVL gene could not be correlated with a specific PFGE pattern profile. All PVL- and sea-positive strains belonged to PFGE type B1 with the exception of one sea-positive strain, which showed PFGE profile B7. The acquisition of toxin genes on mobile elements can go in parallel with extended antibiotic resistance. In our study, the presence of the sea or PVL gene did not change the resistance pattern compared with sea- or PVL-negative strains. As expected, the presence of sea or the PVL gene did not change the SIRU profile or the respective spa type. The fAFLP pattern also was not specifically affected by the acquired toxin genes. In summary, neither of the typing techniques indicated a linkage between genotype and virulence (Hacker et al., 1997).
agr group specificity and agrB-D-C sequencing
On the basis of the agrC 5' end sequence, all isolates in this study could be grouped into agr group I. Sequencing the variable part (agrB-D-C) of the agr locus revealed no nucleotide changes within this region for any of the ST22-MRSA isolates (data not shown).
Antimicrobial resistance phenotypes of the ST22-MRSA isolates
Initially, we examined the antibiotic resistance phenotypes of all of the ST22-MRSA-IV isolates. All were uniformly sensitive to VAN, TEI and FUS. We found four different resistance phenotypes using a panel of antibiotics (PEN, OXA, SXT, TET, ERY, CLI, GEN, VAN, TEI, LIN, FUS, CIP and RIF). The most common phenotype, PEN-OXA-CIP-ERY-CLI, named antibiotype I, was present in 78.6 % (n=33) of the ST22 isolates under study. The resistance pattern PEN-OXA-CIP-ERY, named antibiotype II, was detected in 7.1 % (n=3) of the isolates, whereas the resistance pattern PEN-OXA-CIP (antibiotype III) was involved in 14.3 % (n=6) of the isolates. As is apparent from Table 1
, the most prominent resistance phenotype, PEN-OXA-CIP-ERY-CLI, was predominant in the Department of Dermatology, University Hospital.
Antibiotype I was observed throughout the study period, whereas antibiotype II was only found in the ST22-MRSA isolates in February 2003. Antibiotype III was characteristic of the ST22-MRSA isolates that were obtained between mid-June and November 2003. However, there was no correlation with PFGE pattern, SIRU typing or fAFLP profile, or with toxin profile (Table 2
).
With regard to the strains under study, we detected three different antibiotypes, PEN-OXA-CIP, PEN-OXA-CIP-ERY, PEN-OXA-CIP-ERY-CLI. In 1996, the Barnim epidemic MRSA first appeared in German hospitals with the resistance phenotype PEN-OXA-ERY-CIP-CLI. The majority of the studied strains with PFGE profile B1 (80 %) and all with profile A1, the major restriction patterns of ST22-MRSA-IV, showed the original resistance pattern. Thus, 20 % of PFGE prototype B1 (n=6) and some variants of prototype B1, namely B2, B3 and B7, showed a narrowed resistance pattern compared with the main Barnim MRSA strains. Thus, we did not see a strict concordance between antibiotic susceptibility pattern and PFGE pattern. Changes in DNA content due to deleted resistance genes may be too marginal to change the pulsotype (Lipsitch, 2001). This also applies to SIRU typing, spa typing and fAFLP analysis.
Discriminatory power of the typing methods
Reproducibility was satisfactory with all phenotypic and genotypic typing methods. In accordance with Hunter (1990), a calculator (http://insilico.ehu.es/) was used to determine the discriminatory power of each typing technique used in the study. The most discriminatory methods seemed to be PFGE and fAFLP, followed by SIRU typing and antibiotic resistance phenotyping. The results are given in Table 3
. Although there was concordance in general between the different typing methods, it seemed that the more methods used, the more confusion there was in terms of typing. Due to the fact that each method of typing looked at different genetic markers, no single typing technique was really able to decide which was an outbreak and which was an endemic MRSA strain. Nonetheless, to identify microevolutionary events, we propose a combination of two partially contradictory approaches, classification using the most stable markers by MLST and spa and SIRU typing, and subtyping of the bacterial organism by PFGE and fAFLP, in order to yield maximal information about bacterial dissemination.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
B. Ghebremedhin, M. O. Olugbosi, A. M. Raji, F. Layer, R. A. Bakare, B. Konig, and W. Konig Emergence of a Community-Associated Methicillin-Resistant Staphylococcus aureus Strain with a Unique Resistance Profile in Southwest Nigeria J. Clin. Microbiol., September 1, 2009; 47(9): 2975 - 2980. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. S. Soliman, G. Phillips, P. Whitty, and D. H. Edwards Distribution of meticillin-resistant Staphylococcus aureus spa types isolated from health-care workers and patients in a Scottish university teaching hospital J. Med. Microbiol., September 1, 2009; 58(9): 1190 - 1195. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Conceicao, M. Aires de Sousa, and H. de Lencastre Staphylococcal Interspersed Repeat Unit Typing of Staphylococcus aureus: Evaluation of a New Multilocus Variable-Number Tandem-Repeat Analysis Typing Method J. Clin. Microbiol., May 1, 2009; 47(5): 1300 - 1308. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. J. Wolter, A. Chatterjee, M. Varman, and R. V. Goering Isolation and Characterization of an Epidemic Methicillin-Resistant Staphylococcus aureus 15 Variant in the Central United States J. Clin. Microbiol., October 1, 2008; 46(10): 3548 - 3549. [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | J MED MICROBIOL | MICROBIOLOGY | J GEN VIROL | ALL SGM JOURNALS |