| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Previous Article | Next Article ![]()
Journal of Bacteriology, April 2008, p. 2597-2606, Vol. 190, No. 7
0021-9193/08/$08.00+0 doi:10.1128/JB.01695-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
Richard A. Gibbs,4
Bálint Csörg
,5
György Pósfai,5
George M. Weinstock,4 and
Frederick R. Blattner1,2
DNAStar, Inc., Madison, Wisconsin 53705,1 Department of Genetics,2 Biotechnology Center, University of Wisconsin, Madison, Wisconsin 53706,3 Human Genome Sequencing Center, Baylor College of Medicine, Houston, Texas 77030,4 Institute of Biochemistry, Biological Research Center, H-6726 Szeged, Hungary5
Received 22 October 2007/ Accepted 22 January 2008
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
-complementation. However, analysis of sequenced bacterial artificial chromosome (BAC) clones derived from DH10B shows a high incidence of insertion sequence (IS) transposition from the chromosome into the cloned fragment (25).
The genome of DH10B was constructed before the modern era of molecular biology, through a series of genetic manipulations (Fig. 1). The progenitors were all K-12 strains, with the exception of D7091F, in which a region surrounding the
(araA-leu)7697 deletion had been derived from E. coli B SB3118 by P1 transduction (John Wertz, personal communication). Ultimately, MC1061 (9) served as a starting point for Hanahan and coworkers to replace alleles by using a series of P1 transductions that resulted in DH10B (17). Among the engineered gene replacements were recA1 to improve clone stability by inhibiting the homologous recombination system; endA1, which inactivates the encoded periplasmic DNA-specific endonuclease, thereby enhancing DNA stability during transformation; and a
80 derivative containing the lacZ
M15 mutation for screening by
-complementation. Recombination functions for the latter two steps were provided by expressing RecA from a plasmid that was subsequently cured by treatment with coumermycin (Frederic Bloom, personal communication). The resultant strain, DH10, was also reported to contain an unspecified mutation in deoR which was hypothesized to increase transformation efficiency, though no explanation was offered (19). A DH10 derivative containing the allele mdoB::Tn10 (zjj202::Tn10 in reference 22) was subjected to fusaric acid treatment to counterselect against the tetracycline resistance gene in the transposon, again in the presence of the RecA-expressing plasmid. DH10B was one isolate from this selection with a deletion spanning the marker and the flanking region, including the MDRS loci (mrr, mcrA, mcrB, and mcrC). Deletion of the MDRS loci was shown to improve the cloning efficiency of mammalian DNA, in which cytosine is commonly methylated (17). Later, a strain with a spontaneous mutation in tonA that confers resistance to the bacteriophages T1, T5, and
80 was isolated, and that strain, DH10B tonA, was also commercialized (Invitrogen Corp.).
|
| MATERIALS AND METHODS |
|---|
|
|
|---|
80 (unpublished data) sequences, in order to obtain regions present in DH10B but absent in MG1655. During the finishing phase, gaps between contigs were closed by the direct sequencing of PCR products generated from primers flanking the gaps, and for larger gaps, small insert (SMIL) libraries were generated from the PCR products spanning the gaps. Appropriate primer pairs were generated from the assembly using SeqMan software (DNAStar, Inc.).
Sequence assembly. Genome assembly was done primarily using the SMGA (DNAStar, Inc.). Sequences were preprocessed within SMGA to trim vector and low-quality sequences and to check for more than 900 known repeats. Once the sequences are processed, assembly consists of three steps. First, each sequence read is parsed into overlapping 25-base oligomers and a running count for each oligomer is tabulated. Oligomers occurring more than once in a sequence or more than 1.5 times the expected coverage are marked as repeat oligomers. Sequences are then scanned against the oligomer table, and one nonrepeat oligomer for each 20-base window is chosen (oligomer tag). Oligomer tags are placed in a bucket-and-chain table that includes the position and orientation for the oligomer in each sequence (1). Second, oligomer tags are reorganized into a sequence-sequence overlap table. Each sequence read is scanned against the overlap table to generate a list of sequences that have compatible oligomer tags. Third, contigs are built by sequential addition of overlapping reads. Initial alignments are started using pairs of nonrepeated, overlapping sequences which have mate pair reads that also overlap. Additional constrained sequence pairs are then added, followed by nonrepeat sequences lacking mate pair information. Next, sequences with one nonrepeated end and one repeated end are added. Fully repeated sequences are added to an existing contig only if they have a nonrepeated mate pair within that contig. Finally, repeat sequences that cannot be unambiguously placed into existing contigs are combined into repeat contigs and incorporated manually.
Following the initial assembly, contigs were ordered within SeqMan Pro (DNAStar, Inc.), first using available mate pair data and then alignment information with MG1655 from the Mauve genome aligner (10). When possible, adjacent contigs were merged, resulting in a 36-contig scaffold spanning the genome. Gaps were filled by sequencing of directed PCR products and, where necessary, the corresponding MG1655 sequence. The draft sequence was then used as a template for a complete reassembly by SMGA to maximize coverage and quality. Repeat regions were manually incorporated to ensure accuracy. A second round of PCR-based directed sequencing was used to complete the sequence.
Annotation. The consensus sequence of the first complete (single-contig) assembly was imported into the Web-based ASAP database (15, 16) to aid analysis and annotation by multiple users in different locations. Annotated Feature coordinates were updated when refinements to the assembly produced new genome versions. Genes identical to those in MG1655 were identified using genome alignment data generated by Mauve (10); single-nucleotide polymorphisms (SNPs) and other sequence differences from MG1655 were also identified from Mauve alignments. Pseudogenes were identified by manual inspection and carefully verified. IS elements were identified using RepeatMasker (A. F. A. Smit, R. Hubley, and P. Green, unpublished data; http://repeatmasker.org) searching against the ISFinder database (40). The genome and annotations are freely accessible for viewing in ASAP (https://asap.ahabs.wisc.edu), and annotations may be added or updated by members of the community upon registration.
Growth measurements. Cultures were grown in MOPS (morpholinethanesulfonic acid) minimal medium (32) supplemented with 0.1% glucose and amino acids (40 µg/ml each) as described previously. For each growth curve, a single colony was used to inoculate 50 ml fresh medium in a 250-ml baffled flask and grown at 37°C with shaking in an orbital water bath. Optical density measurements at 600 nm (OD600) were taken every minute using an automated system.
Analysis of mutation rates.
D-Cycloserine resistance assays were performed as previously described (12, 33). Briefly, in a fluctuation assay, 20 tubes of 1 ml mineral salts medium (17a) supplemented with glucose and thiamine were inoculated with 104 cells each, and cultures were grown to early stationary phase. Fifty-microliter aliquots from each tube were then spread on minimal plates containing D-cycloserine (0.04 mM). The estimated number of mutations per tube (m) was calculated from the number of colonies by using the Ma-Sandri-Sarkar maximum-likelihood method (39). Equation 41 from the report of Stewart et al. (41) was used to extrapolate the obtained m value, valid for 50 µl, to 1 ml. Statistical comparisons of m values were made only when the difference in total cell number was negligible (<3%, P
0.6, with a two-tailed, unpaired t test). The total number of cells in a tube was calculated by spreading dilutions from three tubes onto nonselective plates. Dividing the number of mutations per tube by the average total number of cells in a tube gives the mutational rate (mutation/cell/generation).
In a second protocol, appropriate for detection of base substitutions, cells resistant to rifampin were selected. Resistant cells carry base substitutions in rpoB (24). Twenty tubes of 1 ml LB were inoculated with 104 cells each, and cultures were grown to early stationary phase. Appropriate dilutions were seeded onto LB agar plates and LB agar plates containing rifampin (100 µg/ml), and colony counts were performed after 24 or 48 h, respectively. Mutation frequencies are reported as a proportion of the number of rifampin-resistant colonies relative to the total viable count. The results correspond to the mean value obtained in two independent experiments.
Analysis of mutation types. Analysis of the mutational spectrum of the cycA gene has been described previously (12). Briefly, a 1,877-bp genomic segment encompassing the entire gene was amplified from mutant cells using the primer pair cycA1/cycA2. A representative sample was obtained by analyzing 5 colonies from each parallel plate, yielding a total of 96 samples per experiment. The amplified fragments were resolved on an agarose gel and compared to a fragment generated from the wild-type template. Identical sizes indicated a mutation affecting only one or a few nucleotides, a decrease in size or failure of amplification indicated a deletion, and a detectable increase in size indicated an IS insertion. Where further analysis of the insertion mutants was desired, the identity of the ISs was determined by PCR using combinations of oppositely oriented IS-specific primers and primers flanking cycA.
Primers for PCR analysis of mutations. The primers used for PCR analysis of mutations (33) included the following: cycA1, 5'-CTGATGCCGGTAGGTTCT-3'; cycA2, 5'-GCGCCATCCAGCATGATA-3'; for IS1, IS1A1 (5'-TCGCTGTCGTTCTCA-3') and IS1A2 (5'-AAGCCACTGGAGCAC-3'); for IS2, UK1R (5'-TCGCAGGCATACCATCAA-3') and UK2R (5'-CAGACGGGTTAACGGCA-3'); for IS5, IS5ki3 (5'-ATAGGCTGATTCAAGGCA-3') and IS5ki2 (5'-GCTCGATGACTTCCACCA-3'); and for IS150, IS150ki1 (5'-ACGTGCCGAGATGATCCT-3') and IS150ki2 (5'-GCGCCATCCAGCATGATA-3').
| RESULTS |
|---|
|
|
|---|
150 kb) ligated into a BAC vector. During the construction of sublibraries from BAC clones, small amounts of DH10B genomic DNA are present that are occasionally cloned and sequenced. Sequence reads from such "contaminating" clones are then identified in silico and filtered out by comparison against a database containing the MG1655 reference genome (7), known repeated bacterial sequence elements, phages, and other elements. In this way, almost 200,000 sequence reads were collected either entirely from DH10B or from chimeras of bacterial and bovine DNA. All reads were generated by Sanger sequencing technology, and mate pair data were available for most reads. Based on an average read length of 700 bp and a DH10B genome size of 4.6 million bases, these data sets represent an average depth of coverage of approximately 25. The complete genome sequence was then assembled de novo using a new desktop sequence assembly package, the SMGA. The core of SMGA consists of an efficient sequence analyzer that breaks down sequences into a bucket-and-chain hash table of overlapping subsequences, or oligomers (1), and a segmented banded aligner that can rapidly align sequences of any length with nearly constant memory usage (see Materials and Methods). The assembler makes full use of mate pair information, when available, and handles repeat elements in a conservative manner. SMGA is also capable of assembling data sets, either individually or in combination, from traditional Sanger sequencing or massively parallel, next-generation technologies, such as pyrosequencing (29).
The Sanger sequence data set for DH10B was assembled into contigs using SMGA. Scaffolds spanning the genome were organized first using mate pair information and, second, by alignment with the MG1655 genome (7, 35) using the Mauve genome aligner (10). Following gap filling, the resulting draft consensus sequence was then used as a template to reassemble the entire data set using SMGA. Remaining low-coverage areas and consensus ambiguities were resolved using directed sequencing of PCR products corresponding to the areas in question. Finally, the sequence was independently verified with SOLiD (Sequencing by Oligonucleotide Ligation and Detection) technology using the final consensus as a template for assembly (K. McKernan, J. Malek, H. Peckham, F. R. Blattner, and G. M. Weinstock, unpublished results).
Overview of the DH10B genome. The circular genome of DH10B is 4,686,137 base pairs in length (Fig. 2) and can be readily aligned end to end with the wild-type MG1655 sequence with a genome aligner such as Mauve (Fig. 3). As expected, the extensive colinearity between DH10B and MG1655 is also reflected in the gene content. Among the 4,305 protein-encoding genes present in MG1655, 4,058 have identical counterparts in DH10B and another 30 genes contain one or more synonymous SNPs. The remaining 217 genes each differ from their MG1655 counterpart at the protein level, ranging from single-amino-acid substitutions to complete deletions. Other than IS and phage elements, we have not detected any protein-encoding genes that are not also present in MG1655.
|
|
Large-scale alterations of the DH10B genome.
DH10B contains several large structural differences relative to MG1655 (Fig. 3; Table 1). Four known large-scale deletions,
(ara leu)7697,
lacX74,
(mrr-hsdRMS-mcrBC), and precise excision of the e14 prophage, together remove 135,044 bp and 121 genes either completely or partially (Table 1). Excision of e14 also restores a wild-type icd (b1136), although the gene produces a functional protein in the presence of the phage as well (21). These deletions remove multiple operons encoding diverse cellular functions in addition to those targeted for a desired phenotype (see Table S1 in the supplemental material).
|
80dlacZ
M15 insertion is a mosaic element described in detail below. Second, a 113-kb region of genome (genomic coordinates 514341 to 627601) is precisely duplicated in tandem. IS5 elements immediately flanking this segment presumably provided the homology necessary for creating the duplication. The 106 duplicated genes encode functions ranging from membrane components and transporters to transcriptional regulators. Finally, an 11-kb segment at coordinates 3087540 to 3098670 is inverted relative to MG1655. This region is flanked by IS10 elements and can function as a transposon, as evidenced by its identification in the mouse BAC clone AC125523 (GenBank GI:28626891). We confirmed the authenticity of the inversion by PCR.
Polymorphisms. In addition to the 30 genes carrying synonymous SNPs in DH10B, 66 genes have missense mutations relative to MG1655 and 5 genes have nonsense mutations (Fig. 4A; see also Table S2 in the supplemental material). Eleven genes with missense mutations also contain a synonymous SNP. In addition, 42 SNPs are in intergenic regions. While some clustering of SNPs is apparent, particularly in regions of directed strain construction, polymorphisms are found across the genome.
|
The five nonsense mutations (in chiA, gatZ, tonA, yigA, and ygcG), all result in substantially truncated protein products. Approximately half the tonA sequence reads are wild type and half are mutant, indicating that the bovine BAC library was made with competent cells from both the original DH10B strain and DH10B tonA. There is also a sequence ambiguity in the rpsT gene (causing a synonymous CAG-to-CAA change at genomic coordinate 20897) that may represent a mixed population as well.
Four genes are affected by frameshift mutations relative to MG1655: flhC, mglA, fruB, and rph. The 1-bp insertion in rph restores the wild-type reading frame from the rph-1 frameshift found in both MG1655 and W3110, thus alleviating the partial requirement for pyrimidine that is characteristic of both wild-type strains (23).
IS elements and phage elements. The DH10B genome contains significantly more IS elements than MG1655 (63 elements compared to 43) (see Table S3 in the supplemental material). Among the 63 elements in DH10B, 26 are located within coding regions, including all 11 (disrupting 10 genes) found in MG1655 (Fig. 4B). Similarly, all 32 IS elements in MG1655 intergenic regions are the same elements in the same intervals in DH10B, except that the IS5 element located near the flhD promoter in DH10B is replaced by IS2 in MG1655. The DH10B duplication results in an extra copy of both IS5 and IS186B. Three IS elements in DH10B intergenic regions are not present at those sites in MG1655.
Six of the DH10B disruptions, including the previously known galK16 allele, affect carbon source uptake and metabolic pathways. The relA1 allele, which was present in the early progenitor strains (Fig. 1), inactivates the major GDP/GTP pyrophosphokinase responsible for producing (p)ppGpp during a stringent response (30). DH10B also harbors the spoT1 allele (Table 1) which increases basal (p)ppGpp levels by reducing the hydrolase activity of the bifunctional enzyme (14).
MG1655 contains 10 defective prophages, of which 8 are identical in DH10B with the exception of a synonymous SNP in the intR gene of Rac. As noted above, e14 is excised from DH10B; prophage Qin has two alterations, a missense mutant-encoding SNP in ydfU and an IS2 disruption of intQ. The
80lacZ
M15 prophage was added during construction and is not present in MG1655.
In the final sequence assembly, we observed a contig that corresponds to a precisely excised, circular episome of the Rac defective phage element. The ability of Rac to exist as the oriJ plasmid has been documented previously (11). As oriJ plasmids do not replicate in the presence of a Rac prophage (11), the integrated and plasmid forms of the element must be derived from different DH10B cell subpopulations. The oriJ plasmid copy number has not been reported, but the depth of coverage of the episomal contig is twice that of the prophage, arguing that Rac excision and maintenance as a plasmid occur with significant frequency.
80lacZ
M15.
80 is a temperate lambdoid phage and shares a common genetic organization with
(38), although there is limited homology at the nucleotide level (13). Integration of the phage occurs at the attP site in the 5' end of yciI (27). Defective
80 phages have also been used as transfer vectors to move desired alleles between strains (5, 6, 8). Using the complete
80 sequence (G. Plunkett III, unpublished data) as a reference, we evaluated the structure of the
80lacZ
M15 locus in DH10B (Fig. 5). The defective prophage consists of the three following segments: (i) 28,648 bp of the
80 prophage ending within gene 5, (ii) a 12,744-bp piece of the E. coli chromosome extending from the middle of cynS to mhpD which also contains the lac region including the lacZ
M15 allele, and (iii) the 5,983-bp Tn1000 segment of the F plasmid. The 5' end of the element is properly recombined at the attP site, whereas the 3' end consists of the Tn1000 3'-end terminal repeat res, integrated into a partial copy of kch. Downstream of the disrupted kch gene is an intact copy of yciI including an attP site. This organization is consistent with the aberrant excision and recombination events necessary both to generate the defective phage element and integrate it into the recipient genome (5).
|
Comparison of the DH10B sequence with the two wild-type strains argues that DH10B is more closely related to MG1655. DH10B does not contain the W3110 inversion and is identical to MG1655 at 16 out of 21 sites of divergence between the two wild-type strains (Table 2). The exception is the IS5 insertion in the flhD promoter in both DH10B and W3110 rather than the IS2 found in MG1655. Interestingly, DH10B lacks the unique 23S rRNA SNPs of the two other strains [rrlE(G2256A) in W3110, and rrlD(A2547G T2548A G2549T) in MG1655].
|
(ara leu)7697], and indeed, when leucine is added to MOPS minimal medium with 0.1% glucose, DH10B grows aerobically with a doubling time of 76 min while MG1655 has a doubling time of 69 min in the same medium. The disparity in doubling times between DH10B and MG1655 was also observed with additional amino acid supplements. These results demonstrate that DH10B is strictly auxotrophic only for leucine and suggest that the lower growth rates are not due to partial requirements for other amino acids. DH10B mutation rate. The observations of Kovarik and coworkers regarding the high incidence of IS transposition (25), and the number and distribution of SNPs described here, raised the possibility that mutation frequencies are elevated in DH10B compared to those in the wild type. To test this hypothesis, the spontaneous mutation rates of both DH10B and MG1655 were determined using two different assays, namely, (i) a D-cycloserine resistance assay, detecting all types of mutations in the cycA gene (37, 43), and (ii) a rifampin resistance assay, detecting point mutations in the essential rpoB gene (24).
Mutation rate measurements based on D-cycloserine resistance were significantly different between DH10B and MG1655. Data from four independent experiments showed that the total spontaneous mutation rate of cycA in growing DH10B cultures was 13.5 times higher than that of MG1655 (1.07 x 10–6 compared to 7.90 x 10–8, respectively). The mutation spectra obtained from the two strains are also dramatically different (Fig. 6). In MG1655, 74% of the mutations obtained were point mutations, 24% were IS insertions, and 2% were deletions. About half of the IS-mediated disruptions were caused by IS150. In contrast, in DH10B, 95% of cycA mutations were IS insertions, 72% of which were IS150 mediated. Although the proportions of point mutations were vastly different (74% in MG1655 versus 5% in DH10B), the actual rates of point mutations were similar in the two strains (5.86 x 10–8 for MG1655 versus 5.02 x 10–8 in DH10B). No deletions were found among the cycA alleles in DH10B.
|
Together, the results are consistent with the high level of IS transposition observed previously for DH10B (25) and imply that the numerous SNPs found in DH10B are a result of the extensive genetic manipulation that went into the construction of the strain rather than an increased point mutation rate.
| DISCUSSION |
|---|
|
|
|---|
Although the alleles targeted during the various construction steps are basically as expected, the deoR gene is a significant exception. The mutation of deoR was thought to be responsible for the enhanced transformation efficiency of DH10B (17, 19), but its sequence is unambiguously wild type. Mutations of deoR were originally isolated by selecting for mutants that grew rapidly on inosine but not uridine, due to the constitutive activation of the deoCABD operon (31). Using a similar selection scheme, Hanahan observed that the fast-growing strains also had higher transformation efficiencies of large plasmids and assumed that this was also caused by mutation of deoR (19a). The wild-type deoR locus indicates either that the two phenotypes (fast growth on inosine and high transformation efficiency) are completely separable or that another undefined locus (or loci) is responsible. In favor of the latter possibility, the same selection scheme was used to independently isolate the highly transformable DH5 strain (19a) which has now also been shown to contain a wild-type deoR gene (Invitrogen Corp., unpublished results). Interestingly, the multiple-deletion strain, MDS42, has transformation properties similar to those of DH10B (33), raising the possibility that they may share a common subset of mutated genes accounting for the phenotype. Even when pseudogenes and phage elements are excluded, there are still 52 mutated genes in common to be investigated. A systematic investigation of this set and its effect on transformation efficiency is now possible.
The 13.5-fold-higher mutation rate in DH10B than in MG1655 is entirely due to increased IS transposition. This is consistent with previous findings showing a high incidence of IS transposition into eukaryotic BAC library clones (25). In those studies, IS10 was the most frequently observed element in the BAC clones, while IS150 transposition predominated in our study. It is important to note that cycA does not contain a preferred IS10 target site (18, 25), so that IS10 transpositions are expected to be rare. No target site specificity has yet been reported for IS150 or other IS3 family members, but this is under investigation in a separate project. IS150 transposase levels and/or activity could be elevated in DH10B. IS150 transposase production is regulated by a highly efficient programmed translational frameshifting mechanism (42), although precise details of the frameshifting mechanism are only now emerging. No obvious connections with the DH10B genotype are evident.
Conserved IS5 elements were likely involved in creating the large tandem duplication that doubles the gene dosage of 106 genes. Such duplications are quite common in E. coli and Salmonella enterica serovar Typhimurium, probably due to RecA-dependent unequal-sister-strand exchanges between repeated sequences (36), although they are lost at high frequency unless they confer some selective advantage under the given growth conditions. While recA1 in DH10B would allow fixation of the duplication even in the absence of a selective advantage, the three construction steps following introduction of recA1 (Fig. 1) employed a wild-type RecA-expressing plasmid that was subsequently cured from the strain. This implies either that the duplication arose very late in the construction and was fixed by curing of the recA plasmid or that it confers a selective advantage for growth on complex media (e.g., Luria-Bertani or tryptone broth) that are generally used for culturing DH10B. One candidate operon for positive selection is gltLKJI, which encodes the glutamate-aspartate ABC family transporter (28). Cells growing in complex media consume available amino acids in a sequential fashion, with serine and aspartate being used during exponential growth and others such as glutamate used in the transition to stationary phase (4, 34). Doubling the expression of gltLKJI could enhance the uptake of these amino acids, providing a growth advantage.
The range of nutrients that DH10B can utilize is limited by the deletion of numerous metabolic pathways. Nevertheless, DH10B requires only leucine as a supplement for growth on minimal media with a suitable carbon source. The consistently lower growth rates observed with DH10B cultures compared to MG1655 are likely a consequence of elevated basal (p)ppGpp levels caused by the spoT1 allele as seen in different backgrounds (14, 26). Consistent with the "relaxed" phenotype imparted by relA1, however, DH10B does exhibit extensive growth lags during nutrient downshifts (data not shown).
Classical genetic strain construction has been an invaluable tool in elucidating much basic molecular biology. DH10B can be considered an extreme case by the extent of manipulations and the resulting changes to the genome. While most of the changes do not impact the strain itself, IS transposition into cloned fragments propagated in this host sends a strong cautionary signal regarding uncharacterized genomes. Recent advances in sequencing technology and strain construction are now allowing such issues to be eliminated (33).
| ACKNOWLEDGMENTS |
|---|
This work was supported by NIH grant U54 HG003273 at the Baylor College of Medicine HGSC and by an OTKA grant to G.P.
| FOOTNOTES |
|---|
Published ahead of print on 1 February 2008. ![]()
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
Present address: Dupont Agriculture and Nutrition, Johnston, IA 50131. ![]()
| REFERENCES |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Appl. Environ. Microbiol. | Infect. Immun. | Eukaryot. Cell |
|---|---|---|
| Mol. Cell. Biol. | J. Virol. | Microbiol. Mol. Biol. Rev. |
| ALL ASM JOURNALS |