Previous Article | Next Article ![]()
Journal of Virology, April 2009, p. 3127-3137, Vol. 83, No. 7
0022-538X/09/$08.00+0 doi:10.1128/JVI.02474-08
Copyright © 2009, American Society for Microbiology. All Rights Reserved.

Alain Vanderplasschen,1
Michel Georges,2 and
Etienne Thiry1*
Department of Infectious and Parasitic Diseases, Virology and Immunology, Faculty of Veterinary Medicine, University of Liège, B-4000 Liège, Belgium,1 Department of Animal Production, Animal Genomics, Faculty of Veterinary Medicine, University of Liège, B-4000 Liège, Belgium,2 Institut National de la Recherche Agronomique, UR1282, Infectiologie Animale et Santé Publique, F-37380 Nouzilly (Tours), France3
Received 2 December 2008/ Accepted 8 January 2009
|
|
|---|
|
|
|---|
Two types of recombination are described for herpesviruses: illegitimate and homologous recombinations (48). However, the precise mechanism of herpesvirus recombination is poorly understood. It is most likely coupled with viral DNA replication and may require both viral and cellular factors (4, 11, 43). Homologous recombination can occur only between closely related genomes (16, 27). Nonetheless, the creation of new genomes requires some sequence heterogeneity in the genomes involved in recombination. These sequence variations are used as recombination markers. As described above, an intrinsic characteristic of herpesviruses is their high level of intraspecies gene conservation. In the current study, we detected rare single-nucleotide polymorphisms (SNPs) differentiating two wild-type strains of an alphaherpesvirus, bovine herpesvirus 1 (BoHV-1). BoHV-1 is a major pathogen in cattle (31). Primary infection is associated with various clinical manifestations such as infectious bovine rhinotracheitis, infectious pustular vulvovaginitis, abortion, and systemic infection (21, 46). According to antigenic and genomic characteristics, BoHV-1 has been subdivided into two distinct but closely related subtypes: subtypes 1 (BoHV-1.1) and 2 (BoHV-1.2) (31). Although most BoHV-1.1 strains have been isolated from respiratory tract diseases or abortion cases and most BoHV-1.2 strains have been isolated from genital organ lesions, the only reliable distinctive criterion has emerged through the results of viral DNA analysis by restriction endonuclease fingerprinting (13, 42). To date, natural BoHV-1 recombinants generated from the two BoHV-1 subtypes have never been observed.
In this study, we aimed to investigate the impact of recombination on genetic diversity in the BoHV-1 progeny issued from coinfection with two wild-type BoHV-1 strains distinguished by SNP. SNP for virus populations was previously defined as the nucleotide changes shared by a minimum of two isolates in a virus population (14). TaqMan genotyping assays were set up, allowing the allelic discrimination of seven SNPs that spanned the BoHV-1 genome and differentiating two wild-type strains of each of the BoHV-1 subtypes. These SNPs were further used to characterize the progeny virus genotype of around 300 virions randomly sorted from an in vitro coinfection experiment involving BoHV-1.1 and BoHV-1.2. The tool developed in this study overcomes several biases induced by previously used recombination markers such as deletion markers (15, 27, 28, 40), marker rescue techniques (3, 20, 32), and restriction endonuclease (RE) sites (4, 47). The main limits of the latter markers are possible recognition failures of some of the RE sites and a star effect with some of the REs. Moreover, the number of RE sites that differentiates wild-type genomes is not exhaustive. When they are artificially introduced, RE sites do not allow the investigation of recombination between natural variants. In contrast, TaqMan genotyping assays developed in our study showed high sensitivity, specificity, and robustness to investigate recombination between natural BoHV-1 isolates.
This original recombination assay contributed to the discovery of unknown characteristics of recombination in herpesvirus. We found that simultaneous infection with two wild-type strains of BoHV-1 results in a highly diversified recombination mosaic. The high frequencies of multiple recombinants observed in the progeny provided the first evidence of a negative genetic interference in an animal DNA virus. In addition, our findings support the idea that the localization of a gene in the BoHV-1 genome influences its recombination rate.
|
|
|---|
Coinfection experiment. Monolayers of MDBK cells prepared in 25-cm2 flasks were coinfected with the parental viruses BoHV-1.1 Ciney and BoHV-1.2 K22 at a total MOI of 10 using a ratio of 50:50 of each the parental strains. Virus attachment was allowed for 1 h at 4°C. Cells were then further incubated at 37°C. Two hours after the temperature shift, the inoculum was removed, and the cells were washed with phosphate-buffered saline and treated for 1 min using low-pH inactivation of extracellular virions (with a solution at pH 3 containing 40 mM citric acid, 10 mM KCl, and 135 mM NaCl). Cells were then washed twice with Earle's minimal essential medium (MEM) and further incubated at 37°C in 6 ml of MEM supplemented with 5% heat-inactivated fetal calf serum. Twenty-four hours after the temperature shift (when the monolayers showed extensive cytopathic effects), the culture medium was removed, briefly sonicated, clarified twice by centrifugation (1,000 x g), aliquoted, and stored at –80°C.
Isolation and screening of progeny viruses resulting from coinfections. An aliquot of the supernatant from the situation of coinfection was diluted serially in MEM in order to identify the appropriate dilution for individual plaque isolation. Each dilution was used to infect MDBK monolayers cultured in six-well plates (Multiwell; Becton Dickinson). After a 1-h incubation at 37°C, the supernatant was removed, and cell monolayers were overlaid with MEM containing penicillin-streptomycin, 1% (wt/vol) agarose (agar bacteriological; Oxoid), and 10% bovine serum containing BoHV-1-neutralizing antibodies. After 72 h of incubation, individual plaques were picked and propagated individually by the inoculation of MDBK cells grown in 24-well plates with MEM containing penicillin-streptomycin and 2% horse serum.
Virus growth analysis. MDBK cells were infected with the respective viruses at an MOI of 10. After 1 h at 4°C, prewarmed medium was added, and cells were further incubated for 2 h at 37°C to allow virus penetration. The inoculum was then removed, and the cells were washed with phosphate-buffered saline and treated for 1 min using low-pH inactivation of extracellular virions. Cells were then washed twice with MEM and further incubated at 37°C in MEM supplemented with 5% heat-inactivated fetal calf serum. Immediately thereafter and after 5, 10, 15, 20, 25, 30, and 35 h of incubation at 37°C, an aliquot of the culture medium was removed and clarified twice by centrifugation (1,000 x g). Titrations of these supernatants determined by plaque assay on MDBK cells determine extracellular virus production. Infected cell monolayers were frozen with the rest of the medium after the same incubation times. Titrations of the unfrozen cell lysates gave the total virus production. Intracellular virus production was obtained by subtracting the extracellular production from the total production.
Viral DNA preparation. Viral DNA of the parental strains was prepared from clarified supernatants of infected MDBK cell cultures as previously described (22). To characterize the genotype of the progeny viruses obtained in the coinfection experiment, DNA was extracted together with the cellular DNA from approximately 200,000 infected cells by using a small-scale extraction-purification procedure (Qiamp DNA minikit; Qiagen). DNA extracts were stored at –20°C until use.
PCR and sequencing. Seven genes were selected along the BoHV-1 genome (Table 1): UL44 (gC), UL27 (gB), UL11 (tegument protein), LRORF2 (latency-related product), US2 (US2 protein), US6 (gD), and US8 (gE). In order to maximize the rate of detection of single point mutations differentiating BoHV-1.1 and BoHV-1.2, the PCR primers used were designed to target the portions of the genes that showed the greatest divergence between the BoHV-1 and BoHV-5 genomes, with BoHV-5 being the virus most closely related to BoHV-1 (GenBank accession numbers AJ004801 and AY261359, respectively) (Table 1) (7). PCRs were carried out using 5 ng of purified viral DNA as a template in 50 µl containing 15 pmol of both forward and reverse primers, 200 µM of each deoxynucleoside triphosphate, 3% dimethyl sulfoxide, 10 mM KCl, 10 mM (NH4)2SO4, 20 mM Tris HCl, 2 mM MgSO4, 0.1% Triton X-100, and 1 U Taq DNA polymerase (New England Biolabs). PCR started with denaturation for 4 min at 95°C, followed by 35 cycles as follows: denaturation for 30 s at 95°C, annealing for 30 s, and elongation for 1 min at 72°C. Annealing temperatures were set up for each primer pair (Table 1). PCR products were purified by using Wizard PCR preps (Promega). Sequencing reactions for purified PCR products were performed with a BigDye Terminator v3.0 apparatus (Applied Biosystems). An ABI Prism 3730 DNA analyzer (Applied Biosystems) was used for analysis of the samples. All PCR product sequences were verified on both strands.
|
View this table: [in a new window] |
TABLE 1. PCR amplification and sequencing primers used to detect mutations in seven BoHV-1 genes
|
SNP genotyping assay. According to the localization of the SNPs enabling the differentiation between BoHV-1.1 and BoHV-1.2, primers and minor groove binder probes for TaqMan assay were designed using Primer Express v 2.0 (Applied Biosystems, Applera Corporation, CA). Sequences for the primers and probes used are summarized in Table 3. Probes specific for BoHV-1.1 SNPs were labeled with VIC dye, and those specific for BoHV-1.2 SNPs were labeled with 6-carboxyfluorescein (FAM). PCRs were carried out using 10 ng of dried cellular/viral DNA in a final volume of 5 µl composed of 2.5 µl of TaqMan PCR master mix (Applied Biosystems), 0.125 µl of one SNP genotyping assay mixture containing the two primers and the two probes designed for one target SNP, and 2.375 µl of nuclease-free water. The PCR profile was 4 min at 50°C, 15 min at 95°C (activation of the polymerase), and 45 cycles of 15 s at 95°C and 1 min at 60°C. Endpoint reading of the fluorescence generated during PCR amplification was made using the ABI Prism 7900HT apparatus, and genotype assignments were obtained with Sequence Detection System (SDS) software. For each SNP locus, results were plotted on a two-dimensional scatter plot of the BoHV-1.1 "allele" versus the BoHV-1.2 "allele."
|
View this table: [in a new window] |
TABLE 3. Primers and probes used in TaqMan assays to discriminate BoHV-1.1 and BoHV-1.2 at seven SNP positions
|
= 0.05) of the individual progeny configuration were calculated by establishing the lower and upper limits of each of the intervals in a binomial sampling distribution. To test the existence of interference in the three-factor crosses, Fisher's exact test was used. The number of observed double recombinants was compared to the number of expected recombinants under a null hypothesis of independence of two neighboring recombination events. The P value represented the sum of the P values calculated from the hypergeometric distribution in the contingency tables of double-recombination events equal to or higher than the observed frequency. |
|
|---|
![]() View larger version (17K): [in a new window] |
FIG. 1. Growth kinetic curves of two wild-type BoHV-1 strains selected as parental strains in coinfection experiments. Data for intra- and extracellular growth kinetics of BoHV-1.1 (strain Ciney) and BoHV-1.2 (strain K22) were obtained after infection of MDBK cells at an MOI of 10. Titers, expressed as PFU, were recorded after 2, 5, 10, 15, 20, 25, 30, and 35 h postinfection (pi). The data presented are mean values ± 1.96 standard deviations of triplicates ( = 0.05).
|
|
View this table: [in a new window] |
TABLE 2. Mutations differentiating BoHV-1.1 Ciney and BoHV-1.2 K22 obtained by sequence alignment of PCR products in seven less conserved regions of BoHV-1 genome
|
![]() View larger version (21K): [in a new window] |
FIG. 2. TaqMan genotyping assays discriminate seven SNPs differentiating BoHV-1.1 and BoHV-1.2 genomes. (A) Organization of the BoHV-1 genome including two unique sequences, a long one (UL) and a short one (US). The latter is flanked by two repeated and inverted sequences (IR, internal repeat; TR, terminal repeat). Shown are the localizations of the seven point mutations detected in seven BoHV-1 genes that are targeted by TaqMan genotyping assays. The nucleotides indicate the allelic version discriminating BoHV-1.1 (in boldface type) and BoHV-1.2 at each position. (B) Endpoint reading of the fluorescence generated during PCR amplification of seven target genes discriminating BoHV-1.1 ( ) and BoHV-1.2 ( ). BoHV-1.1 DNA hybridized with the VIC-labeled probes, giving rise to an intense fluorescence along the x axis during amplification. BoHV-1.2 alleles hybridized with the FAM-labeled probes, resulting in a high fluorescence signal along the y axis. Different ratios of BoHV-1.1 and -1.2 DNA (1/10, 1/1, and 10/1) were mixed together and submitted to the analysis; these samples ( ) were hybridized with both the FAM- and VIC-labeled probes, resulting in increased fluorescence signals along both axes. Dot plot data were generated from endpoint readings of the fluorescence recorded after TaqMan PCR assays were performed on total DNA quantities ranging from 1 to 20 ng.
|
A total of 304 randomly selected progeny viruses were plaque purified and propagated individually from the coinfection supernatant. DNAs were individually extracted and submitted to the seven SNP genotyping assays capable of discriminating the allelic versions specific for BoHV-1.1 and BoHV-1.2 in the US8, US6, US2, UL0, UL11, UL27, and UL44 genes. With seven SNP sites that could be counted as originating from BoHV-1.1 or BoHV-1.2, there were 27, i.e., 128, possible progeny genomes (2 parental genomes and 126 recombinants). With regard to the signal obtained by allelic discrimination, each of the 304 progeny was counted as either P1, P2, or belonging to 1 of the 126 recombinant classes (Fig. 3). Of the progeny, six BoHV-1 virions were removed from the analysis because both genotypes were assigned to at least one SNP position. These viruses were considered as mixes of either the two parental (n = 1) or the two recombinant (n = 5) BoHV-1 populations. When both genotypes were detected at the seven SNP positions, the population was considered to be a mix of the parental virus. When both genotypes were assigned to six of the SNP positions at maximum, the population was considered to be a mix of two recombinants. Of the 298 monoclonal progeny viruses, 119 were characterized as recombinants, giving an overall frequency of 39.9% of recombinants, and 179 were defined as parental (60.1% of the progeny) (Fig. 3). From the binomial sampling distribution, any configuration that was not represented in the sample (n = 0 in 298 events) had a confidence interval ranging from 0 to 1.2%. Among the parental structures, BoHV-1.2 was five times more frequent than BoHV-1.1 (Fig. 3). Among the recombinant configurations, a high level of diversity was observed, since 35 of the 126 possible configurations were represented. The class of single recombinants was the most abundant, representing 24.8% of the progeny, where recombinants generated from double, triple, quadruple, and quintuple recombination events represented 11.4, 3.0, 0.3, and 0.3% of the progeny, respectively. Taken together, these results show that several mosaic structures combining up to six different segments were generated following the simultaneous in vitro infection with two wild-type BoHV-1 strains. It is also obvious that recombination events occurred along the entire genome, since at least one configuration combining two neighboring SNP markers was detected within both US and UL segments.
![]() View larger version (34K): [in a new window] |
FIG. 3. Characterization and division of the progeny virus genotypes issued from in vitro coinfection with BoHV-1.1 and BoHV-1.2. A total of 298 randomly plaque-purified virions issued from a BoHV-1.1 and BoHV-1.2 coinfection experiment were submitted to TaqMan genotyping assays at seven selected SNP sites in BoHV-1, the US8 (1), US6 (2), US2 (3), UL0 (4), UL11 (5), UL27 (6), and UL44 (7) genes. The allelic discrimination data are represented by small squares at each SNP position. Filled squares ( ) indicate the alleles inherited from BoHV-1.1; empty squares ( ) indicate the alleles inherited from BoHV-1.2. Progeny virus genotypes were distributed in 37 different SNP combinations: 2 parental (P1 and P2) and 35 recombinants (R1 to R35). The configurations of recombinants were further classified as single, double, triple, quadruple, and quintuple recombinants by counting the recombination events present in the virus genotype. Quantitative data are presented as the absolute number of genotyped viruses classified in each configuration and as the relative frequency in the total progeny. The confidence intervals ( = 0.05) were calculated by determining the lower (pi) and upper (ps) limits of frequency in a binomial sampling distribution.
|
|
View this table: [in a new window] |
TABLE 4. Interference in the three-factor crossesc
|
|
View this table: [in a new window] |
TABLE 5. rF and RF values between any SNP site calculated from the progeny virus genotypes generated after BoHV-1.1 and BoHV-1.2 coinfectionb
|
|
|
|---|
The high frequency of multiple recombinants observed in the progeny provided the first evidence of negative genetic interference in herpesviruses. In the two BoHV-1 regions where significant negative interferences were measured, the occurrence of a first recombination event positively influenced the occurrence of a second recombination event in the neighboring region. This fundamental observation would need further investigations both to determine its influence on the generation of natural recombinants in herpesvirus populations and to characterize the mechanism involved in this interference. Recombination analyses suggest that varicella-zoster virus strains belonging to the M1 and M2 genotypes are mosaic recombinant strains originating from ancestral isolates belonging to the E and J genotypes through recombination on multiple occasions (35). Similar conclusions were drawn from recombination analyses performed using other natural recombinant isolates of varicella-zoster virus (37), human simplex virus 1 (HSV-1), and HSV-2 (1, 33, 34). From the negative interference obtained in our recombination assay, it would be interesting to investigate whether these recombinants were generated from a single coinfection event producing multiple recombinant genomes or through successive occasions of single recombination events.
Two main hypothesizes could explain the rise of an excess of multiple recombinants in our analyses. This could result either from double-CO events, whose cooccurrence would have been favored by an unidentified factor, or from gene conversion, a process occurring in the heteroduplex region formed during homologous recombination (5). The relatively short distance separating the recombination markers used in our analysis argues in favor of the latter hypothesis. Indeed, recombination analysis of closely linked markers in both prokaryotic and eukaryotic organisms has generally revealed a value of c much larger than 1 and produced a large negative value of I. This type of result has been named high negative interference (50). In these situations, double recombinants have been shown to appear much more frequently than expected. The high frequency of apparent double recombinants is likely to be the result of a gene conversion occurring at mismatched positions in the region of the heteroduplex expanding from the initial recombination process (19, 39) (Fig. 4). Following a single homologous recombination, repairs convert the nucleotide at the neighboring SNP included in the heteroduplex of the initial CO, giving an "artifactual" result of double recombinant (5).
![]() View larger version (16K): [in a new window] |
FIG. 4. Schematic representation of gene conversion occurring in the heteroduplex region formed during homologous recombination. (A) The two parental DNA molecules (P1 and P2) are aligned before recombination. The initial CO takes place by an exchange of single strands between the parental DNAs. (B) Once the cross-bridge has formed, its position can diffuse along the paired molecules. This process can create long regions of heteroduplex DNA. Mismatched positions in the heteroduplex are indicated by asterisks. (C) The Holliday structure can be redrawn, as indicated in B, after a rotation of 180° around the cross-bridge. The Holliday structure might be resolved either by West-East or by North-South cuts. (D) West-East cut yields linear molecules that are not recombinant for external genetic markers. However, repairs of mismatches within heteroduplex-combining alleles of the two parental DNAs might create apparent double recombinants through gene conversion in favor of one of the parental templates. (E) North-South cut yields linear molecules that are recombinant for parental genetic markers on either side of heteroduplex. Repairs of mismatches might also create multiple recombinants (R4) through gene conversion.
|
Our analysis of 300 progeny BoHV-1 strains issued from coinfection between two wild-type BoHV-1 strains differentiated by seven SNP markers revealed an overall rF of 40%. Compared to previous results obtained with nonselected recombination markers, this is twice as high as the overall rF observed in 1,019 HSV-2 progeny (21.8%) by Brown and collaborators (4), who used two viruses originating from the same strain background but differing by five RE sites as unselected markers. However, our result is lower than the rF observed for 93 HSV-1 progeny (54.8%) by Umene (47), who used two strains differing in eight RE sites. More interesting comparisons can be made at the level of RF expressed as a percentage of recombination per kbp. In the present study, we found an overall frequency of 0.58% per kbp in BoHV-1, with significant differences between the L and S segments. This frequency is lower than the 0.7% per kbp previously observed in HSV-1 by Umene (47), who did not obtain any evidence to support the existence of enhanced recombination events in any region of HSV-1. Brown and collaborators (4) previously observed lower overall frequencies of around 0.43% per kbp in HSV-2 and 0.45% per kbp in HSV-1, although they used parental strains that were assumed to be highly homologous. The discrepancies among these three observations could be due to (i) virus species variations (HSV-1 and -2 belong to the genus Simplexvirus and have a class E genome structure, and BoHV-1 belongs to the genus Varicellovirus and has a class D genome structure), (ii) technical limits in the analyses based on RE sites, (iii) the effect of the strains involved in the respective studies, and (iv) an effect of the cell line on which the coinfections have been performed. The latter factor was previously identified as having a deep impact on the recombination process efficiency by Brown et al. (4). Nonetheless, the RF levels found in all the above-mentioned studies are in relatively good agreement and correlate with the RF of 0.5% observed in another DNA virus, human adenovirus type 5 (51).
The present study describes the use of isolated nucleotide mutations as markers in investigating recombination dynamics in herpesviruses. The use of SNPs as recombination markers is advantageous, as large numbers of markers can be detected in any kind of gene, essential or nonessential. Moreover, this is a robust method, as the presence of one allele version at one SNP locus is an all-or-nothing condition. Reversion is also expected to be a rare event, and it was not observed under our conditions. Another benefit of the tool developed in this study is its capacity to detect viruses showing "heterozygosis." Six plaque-purified viruses were withdrawn from our analysis, as they were considered to be mixes of either the two parental (n = 1) or the two recombinant BoHV-1 (n = 5) populations. However it can be reasonably postulated that some of these configurations were originally produced from a monoclonal recombinant virus harboring a "nonrepaired" mismatched position(s) in the heteroduplex region involved in recombination. Under this hypothesis, the packaging and the release of viral genomes harboring heterozygosis would participate in the genetic diversification following virus replication.
SNP genotyping is also a useful method that allows the investigation of recombination at multiple positions. Although we restricted our investigation to seven loci, we were able to detect several recombination events and recombinant structures that were invisible by the use of two deletion markers (27, 28, 30, 40, 41). In addition, the use of seven SNP markers enabled us to partially answer this question of "invisible recombinants," but an increased number of target SNPs along the BoHV-1 genome would be needed to resolve this definitively. The better accuracy of this method for investigating recombination is demonstrated by the differences observed between rF and RF in the quantitative assessment of recombination between widely separated markers (such as US8 and UL44). It is obvious that the underestimation provided by the rF results from the invisible recombinants that cannot be detected when distantly separated markers are used.
Regarding the numerous genetic structures (35 out of the 126 possible recombinant combinations) obtained in vitro, the question of the in vivo rise of such recombinants needs to be addressed. A previous study showed that recombinant viruses represent up to 30% of the progeny population isolated from cattle infected with two BoHV-1 mutants (40). Therefore, the rise and the evolution of a recombination mosaic could be monitored following the coinoculation of cattle with the two parental BoHV-1 strains used in this study. The in vivo replication context is much more complex. First, multiple growth cycles are observed during primary BoHV-1 infection. If recombinants are present in the first progeny, they will have the opportunity to coinfect new target cells, giving rise to new gene combinations. Second, several selection factors (the presence of immune effectors and the adaptability of the strains to cattle) may influence the fitness of progeny viruses arising from these coinfections. In vitro, BoHV-1.2 seemed to be favored when cells were coinfected by both strains. The situation could be totally different in vivo. A BoHV-1.1 epidemiological overdominance was observed in the 1970s when BoHV-1.1 strains were introduced into European cattle herds infected hitherto by BoHV-1.2 strains (12). Finally, the establishment of a latent infection with both strains increases the likelihood of cellular coinfection following a reactivation stimulus acting as a virus synchronizer. It was previously demonstrated that cattle can support the latent infection of two distinguishable BoHV-1 strains, namely, BoHV-1.2 K22 and BoHV-1.1 Cooper (49).
In conclusion, the use of SNPs as recombination markers detected by TaqMan assays provides an accurate tool for the further investigation of recombination in herpesvirus both in vitro and in vivo. The data obtained in this study demonstrate that in vitro, (i) BoHV-1 recombines at a high rate, (ii) multiple recombination events occur within one BoHV-1 growth cycle, (iii) BoHV-1 displays a negative interference, and (iv) levels of RF are different within the US and UL sequences.
This study was supported by grants from the Fonds National de la Recherche Scientifique (grants FRFC 2.4532.98, FRFC 2.4508.02, and FRFC 1.5.105.03). B. Muylkens is a postdoctoral researcher for the FNRS.
Published ahead of print on 19 January 2009. ![]()
Present address: Division of Virology, GlaxoSmithKline Biologicals, Rixensart, Belgium. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»