Previous Article | Next Article ![]()
Journal of Virology, January 2005, p. 67-78, Vol. 79, No. 1
0022-538X/05/$08.00+0 doi:10.1128/JVI.79.1.67-78.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Molecular Biology,1 Department of Medical Microbiology and Immunology, University of Aarhus, Aarhus, Denmark,3 Department of Comparative Medicine, GSF-National Research Center for Environment and Health, Neuherberg, Germany2
Received 31 May 2004/ Accepted 26 August 2004
|
|
|---|
|
|
|---|
In contrast to the retrovirus-like Ty retrotransposons in yeast, which are very selective in the choice of integration sites (as reviewed in reference 9), retroviruses integrate throughout the chromosomes (10, 22, 33, 37, 42, 43, 46, 49, 50, 63, 65, 67, 70, 75, 76). Although studied in vivo and by use of simplified in vitro models during the last decades (for a review, see references 8 and 36), integration site selection still remains poorly understood. Based on these reports, factors such as nucleosomal structure, DNase I-hypersensitive sites, and methylation seem to affect integration (44, 55-58, 60, 73). Moreover, genes appear to be favored targets for both murine leukemia viruses (MLVs) and human immunodeficiency virus type 1 (HIV-1) as examined in cell cultures (63, 76). In contrast to MLV, which prefers integration near the start of transcriptional units, the entire transcriptional unit except upstream of the transcriptional start is favored by HIV-1 (76).
Simple non-acutely transforming retroviruses induce hematopoietic malignancies by a complex process including insertional mutagenesis of host genes (as reviewed in references 41 and 51). Extensive analyses of proviral integration sites in mice, cats, rats, and birds have revealed that c-myc is one of the most frequently targeted genes (for a review, see references 20, 41, and references therein). In chickens, 3' promoter insertion is the predominant form of activation, while c-myc expression is deregulated primarily by enhancer activation in mammals. In mice, both MLVs of the gammaretrovirus genus (e.g., Moloney MLV [Mo-MLV], SL3-3, and MCF 69L1) and the thymotropic betaretroviral leukemia virus (TBLV), which is closely related to mouse mammary tumor virus (5), target this proto-oncogene, giving rise to hematopoietic malignancies such as T-cell lymphomas and erythroleukemias (3, 6, 19, 33, 41, 43, 49, 54, 59, 67, 70). In SL3-3-induced T-cell lymphomas, the proviral insertions are predominantly located upstream of the first exon and the majority of insertions are integrated in the transcriptionally opposite orientation relative to that of c-myc (3, 43, 67).
One of the major leukemogenic determinants of MLVs is the transcriptional enhancer containing densely packed binding sites for a variety of transcription factors. MLVs are able to infect a broad spectrum of cell types; however, tumor development is favored in target cells expressing the array of transcription factors that matches the profile of the enhancer framework of the infecting retrovirus (41). As we have previously reported, the introduction of an extra wild-type repeat in combination with the deletion of two NF1 (nuclear factor 1) binding sites in the SL3-3 enhancer region generates a potent inducer of T-cell tumors, the SL3-3(turbo) virus (also known as SL3-3[2
18-3]) (26). In contrast to wild-type SL3-3 [SL3-3(wt)]-induced tumors, of which 20 to 25% display clonal rearrangements in the c-myc locus due to proviral insertions (35, 48, 53, 67), initial data from lymphomas induced by this enhancer variant demonstrated no such rearrangements (23). However, subsequent PCR analyses of proviral integration sites in 12 tumors induced by SL3-3(turbo) revealed that 92% of these harbored nonclonal insertions in c-myc (A. A. Nielsen, A. B. Sørensen, and F. S. Pedersen, unpublished data). In order to make a thorough examination of this observation, about 250 lymphomas induced by SL3-3(wt) and a variety of enhancer variants of SL3-3 were tested for insertions in the promoter region of this proto-oncogene. Results from this study point to (i) integration site hot spots in the c-myc promoter region, (ii) frequent cases of various atypical provirus-host junction structures, (iii) a virus-dependent distribution of integration sites, and (iv) a correlation between latency and the number of proviral insertions detected per tumor.
|
|
|---|
DNA. Genomic DNA was extracted from frozen tumor tissues with a DNeasy tissue kit (QIAGEN) according to the manufacturer's instructions.
PCR amplification. c-myc-specific PCRs were performed with a 50-µl volume containing 5 µl of 10x Taq buffer (Invitrogen), a 0.2 mM concentration of each deoxynucleoside triphosphate (Invitrogen), 1.25 U of Taq DNA polymerase (5 U/µl; Invitrogen), 1.5 mM MgCl2 (Invitrogen), and 10 pmol of each primer (see below). One hundred to 1,000 ng of genomic tumor DNA was used as the template in the PCRs. The primers were as follows: v1, 5'-XGAATTCGATATC GATCCCCGGTCATCTGGG-3'; v2, 5'-XTGCGGCCGCGATTCCCAGATGACCGGGGATC-3' (the underlined sequences in v1 and v2 anneal to the viral sequence, and the remainder of the primer sequences consists of linker sequences added for other purposes, with X being biotin); myc1, 5'-TGTGTATGTATACGTTTGGGGATTGTAC-3'; and myc2, 5'-CACTCCAGCACCTCCGGTTCGGACT-3'. The proviral primers v1 and v2 correspond to positions 8204 to 8187 and 8187 to 8204, respectively, of GenBank and EMBL database accession number AF169256. The two gene-specific primers myc1 and myc2 correspond to positions 77 to 104 and 710 to 686, respectively, of GenBank and EMBL database accession number M12345. Oligonucleotides were synthesized at DNA Technology ApS, Aarhus, Denmark. The fragments were amplified in a TouchDown thermal cycler (Hybaid) with the following program: 1 cycle of denaturation at 94°C for 3 min and then 40 cycles of denaturation at 94°C for 1 min, annealing at 62°C for 1 min, and extension at 72°C for 3 min, and finally 1 cycle of extension at 72°C for 10 min. Subsequently, the amplification products were visualized on ethidium bromide-stained 1.5% agarose (Invitrogen Life Technologies) gels in 0.5x Tris-borate-EDTA (Invitrogen Life Technologies) buffer. The GeneRuler 1-kb DNA ladder and the GeneRuler 100-bp DNA ladder were purchased from Fermentas.
Purification of PCR products. Amplified PCR products were purified by using streptavidin-coated magnetic beads (Dynabead M280-streptavidin; Dynal AS, Oslo, Norway), the Wizard DNA clean-up System (Promega), or the GFX PCR DNA and GelBand purification kit (Amersham Biosciences) according to the manufacturer's instructions.
Sequencing procedure and sequence comparison. The amplified PCR products were sequenced with the DYEnamic ET terminator cycle sequencing kit (Amersham Pharmacia Biotech) by following the manufacturer's recommendations, and reaction products were analyzed on an automated DNA sequencer (genetic analyzer model 3100; Applied Biosystems Inc.). Amplification products harboring the 5' long terminal repeat (LTR) or 3'-LTR regions were sequenced with primer v3 (5'-CTCTGGTATTTTCCCATG-3') and primer v4 (5'-TCCGAATCGTGGTCTCGCTGATCCTTGG-3'), respectively. Sequencing primers v3 and v4 correspond to positions 7904 to 7886 and 69 to 96, respectively, of GenBank and EMBL database accession number AF169256. Both oligonucleotides were purchased from DNA Technology ApS, Aarhus, Denmark. Sequencing of enhancer regions was performed by using the viral primer v1. Sequences were edited by use of the program Sequencer, version 3.0 (Gene Codes Corporation), and Vector NTI (InforMax, Inc.). Edited sequences were compared with available sequences in databases by using the nucleotide-nucleotide BLAST (basic local alignment search tool) search tool with comparison to sequences in the nonredundant database (http://www.ncbi.nlm.nih.gov/BLAST/), the mouse genome BLAST search tool (http://www.ncbi.nlm.nih.gov/genome/seq/MmBlast.html), and the Ensembl (http://www.ensembl.org/Mus_musculus/) and University of California at Santa Cruz (assembly date, October 2003) (http://genome.UCSC.edu/) assemblies.
Statistical analysis. The difference in the distributions of target sites in c-myc selected for during tumor development induced by SL3-3(wt) and SL3-3(turbo) was tested by use of Fisher's exact test for count data (1).
The proviral insertions in c-myc per tumors were analyzed in groups with respect to the latency periods, assuming a Poisson distribution. This assumption was tested by calculating a conditional chi-square statistic. This test statistic was shown to work well even for small Poisson
parameters (7). The Poisson parameter is simply estimated by the average number of proviral insertions (mi) in each group. The confidence intervals for the Poisson parameters are based on the following equations:
1 =
T,2:
/2/(2n) and
u =
T + 1,2;1
/2/(2n), where 
,ß;
is the deviate associated with the lower tail probability
of the gamma distribution with a shape parameter
and scale parameter ß (27).
To test whether the parameters of two independent Poisson random variables (
1 and
2) are equal, a binomial test is used (40). This test is based on the number of successes (n1 · m1), the number of trials (n1 · m1 + n2 · m2), and the hypothesized probability of success [n1/(n1 + n2)], where mi and ni are the average number of insertions and the sample size in group i (i = 1 or 2).
All the statistical calculations were done with Ra language and environment for statistical computing and graphics which is available as free software (www.r-project.org) (38).
|
|
|---|
![]() View larger version (29K): [in a new window] |
FIG. 1. Enhancer structures of wild-type and mutant SL3-3 MLVs used in the study (24, 25, 26). (A) The SL3-3(wt) enhancer located in U3 consists of a 72-bp direct repeat followed by a third repetition of 34 bp containing binding sites for a variety of transcription factors. Both SL3-3(turbo) and SL3-3(2 18-2) harbor two identical 18-bp deletions encompassing the NF1 site. In addition, the SL3-3(turbo) enhancer contains an extra 72-bp wt repeat. SL3-3(GTT) and SL3-3(TUMdm) are both Runx binding site mutants. The SL3-3(GTT) enhancer contains the GTT-to-TGG mutation in both Runx site I binding sites. The SL3-3(TUMdm) enhancer harbors several alterations compared to the SL3-3(wt) enhancer structure. In addition to an extra 72-bp repeat, the enhancer harbors two identical 28-bp deletions encompassing the NF1 site. Furthermore, GTT-to-TGG and GAC-to-TCA mutations are present in Runx site I and site II, respectively. A T-to-G base substitution is present in the last Runx site I. (B) The enhancer structures are illustrated schematically with boxes; deletions are marked by gaps, and Runx binding site mutations are indicated by letters. WT, wild type.
|
![]() View larger version (59K): [in a new window] |
FIG. 2. Detection of proviral integrations in the c-myc promoter region. (A) Provirus integrations in the c-myc promoter region in the SL3-3(wt) and variant-induced tumors were detected by using a PCR method employing two virus-specific (v1 and v2) and two gene-specific (myc1 and myc2) primers. PCRs with the four primer combinations v1-myc1, v2-myc1, v1-myc2, and v2-myc2 were performed in order to detect viruses integrated in either orientation. The four base pairs b1-b4' to b4-b1' duplicated during proviral integrations of MLVs are shown. We have defined the first base pair (in the c-myc sense orientation) in this duplicated stretch of bases (b1-b4') to be the site of integration (boxed). (B) The amplified PCR products were visualized on ethidium bromide-stained agarose gels. As examples, amplification products obtained from 10 SL3-3(turbo)-induced tumors in male mice are shown. M1, 1-kb ladder (2, 1.5, 1, 0.75, 0.5, and 0.25 kb); M2, 100-bp ladder (1,031, 900, 800, 700, 600, 500, 400, 300, and 200 bp, respectively); N, negative control; P, positive control. Tumor numbers and primer combinations are listed above the lanes.
|
|
View this table: [in a new window] |
TABLE 1. Frequency of tumors harboring integrations into c-myc in different experimental series
|
![]() View larger version (18K): [in a new window] |
FIG. 3. Integration sites in the c-myc promoter region. Integration sites in the c-myc promoter obtained by gene-specific PCRs are shown relative to the promoter sequence with the GenBank and EMBL accession number NT_078782. The region from approximately 800 to 1,800 bp upstream of exon 1 is enlarged in order to illustrate the densely packed integration sites. Each proviral integration is indicated by a box, and the orientations are contrasted, with white indicating the sense direction and grey indicating the antisense direction.
|
|
View this table: [in a new window] |
TABLE 2. A typical virus-chromosome junction structures
|
Generally, the single point mutations introduced a V-to-T (V denotes A, C, or G) nucleotide alteration at the position juxtaposed to the virus or at the second nearby position; i.e., for viruses integrated in the antisense orientation, the nucleotide of the target sequence at position b1 or b2 (the nomenclature is shown in Fig. 2A) is altered from V to T when the mutation is detected at the 5'-LTR-flanking junction, and the nucleotide at position b1' or b2' is altered from V to T at the 3'-LTR-flanking junction. The opposite is seen for viruses integrated in the sense orientation. At all three junctions harboring a dinucleotide alteration, two identical nucleotides, CC or GG, were mutated. The mutations were observed equally in 5'-LTR-flanking and 3'-LTR-flanking repeat structures (45 and 55%, respectively). All together, the observed nucleotide changes suggest that defects in one of the steps in the integration process, perhaps due to an imperfect removal of the viral 5' dinucleotide left over, introduce mutations in the stretch of 4 bp duplicated at the MLV integration.
A [5'-(A/C/G)TA(C/G/T)-3'] consensus sequence for SL3-3 MLV target sites. In total, the 164 insertions were distributed on 83 different integration sites. To clarify whether the selection of target sites correlates with a consensus sequence, the repeats at the integration sites displaying a typical MLV-generated duplication of 4 bp (Table 3) were compared. At all four positions, the percentage of each nucleotide was calculated and the results are presented in Tables 4 to 6. The calculations were based on nucleotides located at the first to fourth positions from the 5' LTR of the virus and into the flanking DNA; i.e., for viruses in the sense orientation, the four nucleotides are 5'-b1'b2'b3'b4'-3', and for viruses in the antisense orientation, the nucleotides are 5'-b1b2b3b4-3' (the nomenclature is shown in Fig. 2A). As seen, a T nucleotide at the second position and an A or T nucleotide at the third position is favored for viruses in either orientation (Tables 4 and 5). Several integration sites are repeatedly targeted, and when we include this fact in the calculations (Table 6), the data suggest a weak consensus sequence [5'-(A/C/G)TA(C/G/T)-3'] for SL3-3 MLV target sequences. The apparent palindromic nature of the consensus sequence indicates that IN does not distinguish between 5'-LTR or 3'-LTR att elements at the strand transfer step.
|
View this table: [in a new window] |
TABLE 3. Sequences of the base pair repeats identified at the 83 different integration sites
|
|
View this table: [in a new window] |
TABLE 4. Distribution of each nucleotide at the four positions in 21 target sequences from viruses in the sense orientationa
|
|
View this table: [in a new window] |
TABLE 5. Distribution of each nucleotide at the four positions in 72 target sequences from viruses in the antisense orientationa
|
|
View this table: [in a new window] |
TABLE 6. Distribution of each nucleotide as calculated for all of the 160 proviral insertions with 4-bp duplicationsa
|
![]() View larger version (25K): [in a new window] |
FIG. 4. Virus-dependent integration pattern in the c-myc promoter region. (A) The promoter region is divided into sections of 30 bp, and within each section the percentages of all virus insertions are calculated. (B) Similarly, the percentages of SL3-3(wt) and SL3-3(turbo) insertions within 30-bp windows from positions 827 to 1756 are shown (this region is indicated by lines in panel A; i.e., a total of 10 insertions outside this region are not included). (C) The numbers of SL3-3(wt) and SL3-3(turbo) insertions at individual positions from positions 1396 to 1427 are illustrated.
|
![]() View larger version (32K): [in a new window] |
FIG. 5. Correlation between tumor latency and number of proviral integrations in the c-myc promoter region. (A) The number of provirus insertions in the individual tumors is shown as a function of the latency period. Identical data points (latency period and number of integrations) are shown above each other. The numbers of different integrations in c-myc in the individual tumors are obtained by counting the numbers of amplification products achieved in the repeated rounds of PCRs. (B) The average number of proviral integrations in the c-myc gene as a function of latency period to disease development is shown along with 95% confidence intervals. The total observation period of 300 days has been divided into subgroups of 10 days each; however, due to the small amount of data samples for each latency period above 80 days, these latency periods have been grouped. N, number of mice within a given subgroup.
|
i of the Poisson distribution in subgroup i. Figure 5B shows these estimators in common with 95% confidence intervals based on the gamma distribution (see Materials and Methods). There is a strong correlation between the average number of insertions and the grouped latency periods. Comparing two adjacent latency periods showed significant differences (binomial test P values were 0.031, 0.006, 0.246, and 0.043). In view of these stringent conditions, there was no significant difference in Poisson parameters between the latency period from 60 to 69 days and the latency period from 70 to 79 days. However, this may be due to the relatively small number of samples in these groups. |
|
|---|
The integration sites are not uniformly distributed along the promoter sequence (Fig. 3), and the insertion pattern both reflects target site selection by IN and represents sites selected for during cancer development. Several of the targeted positions are separated by a few base pairs only, and we believe that this is due to intrinsic features of IN regarding target site selection rather than divergences in the potential of the integrated proviruses to deregulate the expression of c-myc.
Single point mutations in the direct repeat element were observed at either the 5'-LTR- or 3'-LTR-flanking junction for 17 of the integration sites (Table 2). In thirteen of the mutated direct repeats, the nucleotide immediately adjacent to the virus was altered, and in four integration events, the second nearby nucleotide was changed. In all cases except for one, a V-to-T mutation was introduced. Based on this observation, we propose the model for error-prone gap repair illustrated in Fig. 6. By this model, mutations are introduced when the repeated target sequence harbors a T nucleotide at either the first or second position from the virus that is able to base pair with one of the two A nucleotides in the viral 5' overhang, thereby generating a region of microhomology. The DNA polymerase involved in gap repair fills in the dinucleotide gap, and subsequent editing of the mismatching base introduces the mutation.
![]() View larger version (38K): [in a new window] |
FIG. 6. Model for error-prone gap repair at proviral integration. The figure illustrates the introduction of mutations at either virus-chromosome junction. (A) Schematic representation of the SL3-3 proviral DNA copy of the genome. The 3'-LTR and 5'-LTR att sites are shown. (B) The integration step is initiated by 3' processing of two T nucleotides. (C) During strand transfer, the two DNA ends are joined in concert to the chromosomal target DNA 4 bp apart. The target sequence VTAB (V denotes A, C, or G; B denotes C, G, or T) is shown, and the T nucleotide located at the second position is underlined. Due to base pairing, the 5' dinucleotide overhangs are not removed during the gap repair process. (D) Subsequent editing of the mismatching base pairs (marked with an asterisk) generates the V-to-T mutation at either junction. The figure illustrates the introduction of mutations at both junctions for the same proviral integration; however, it is noted that this is never observed in our data set.
|
The gene-specific PCRs were performed with a regular Taq polymerase without proofreading capacity (see Materials and Methods). However, as the mutations are position specific (never observed at the third or fourth position from the virus) and nucleotide specific (V to T) and only a very few differences were detected during comparison with sequences in publicly available databases, we can rule out the possibility that the observed base changes resulted from errors generated during the elongation step.
The gap repair system involved at proviral integration still remains to be determined. Based on studies with cell culture, both DNA-dependent protein kinase participating in the repair of double-stranded DNA breaks and the poly(ADP-ribose) polymerase-1 (PARP-1) implicated in the repair of single-stranded as well as double-stranded breaks have been proposed to be involved (16-18, 29, 31). However, in both cases, investigators have obtained deviating results (4, 13, 21, 28, 47). In addition, both viral and host polymerases have been shown to repair gapped DNA substrates in concert with flap endonuclease and host ligases in vitro (77). The point mutations introduced at error-prone gap repair may be generated due to the mismatch repair system responsible for removal of base mismatches introduced during replication or by other processes (reviewed in the work of Christmann et al. [12]). The mechanism of strand discrimination in eukaryotes is still not clear, and mismatch repair may correct the base present at either strand of the mismatching base pair, leading to introduction of point mutations in some cases only. Alternatively, DNA replication may have taken place prior to repair of the mismatching base pair, thereby generating daughter cells harboring different junction structures.
The 17 single-nucleotide mutations generated by error-prone gap repair were detected at 10 different integration sites dispersed regularly along the c-myc promoter region, and 9 of the cases were located at three target site hot spots (at positions 1306, 1199, and 1112). At position 1306, error-prone gap repair was detected at seven out of nine insertion events. The target sequence at this hot spot is 5'-ATAC-3' (Table 3), which provides the opportunity of base pairing between target sequence and the viral 5' overhang at either junction (Fig. 6), a fact that may partly explain this high frequency. In contrast, none of the nine proviral integrations located at the integration site hot spot at position 1397 also harboring an optimal target sequence for error-prone gap repair (5'-GTAC-3') (Table 3) were correlated with point mutations in the junction structures. This finding raises the possibility that the location of the insertion site may play a role in the frequency of error-prone gap repair.
Atypical lengths of the duplicated repeats at virus-host DNA junctions have been reported previously (14, 30, 39, 52, 71, 72). Likewise, one case of point mutations similar to those described in this report was identified at proviral integration sites by a Mo-MLV-based vector in fibroblasts by Jin et al. (39). Of the 38 integration sites analyzed by Jin et al. (39), 6 were correlated with a duplication of a 5-bp repeat, 1 displayed single-nucleotide point mutations, and 1 was correlated with another type of aberrant junction structure. Thus, an atypical length of repeats is the foremost alteration in the study by Jin et al. (39), while point mutations predominate in our data set. Differences in the cell types, chromosomal regions affected, and integrating mutagens, i.e., virus versus vector, may explain this divergence. Differences in the IN enzyme among MLVs might also play a role. Our data stem from T cells, in contrast to the fibroblasts used with the Mo-MLV-based vector. Nonhomologous end joining plays a central role in T cells due to the involvement in rearrangements of the T-cell receptor, a fact that may relate to the higher frequency of error-prone gap repair in our data set. Moreover, these cells are tumor cells in which gap repair systems may be impaired. In addition, the importance of the organization of chromatin in different cell types cannot be excluded.
The different target sites conformed to a weak consensus sequence [5'-(A/C/G)TA(C/G/T)-3'], in which a T and an A nucleotide at the second and third positions, respectively, are favored (Table 6). Previously, the weak consensus sequence [5'-GT(A/T)AC-3'] was identified for HIV-1 integration (10, 69). As with our observation, a bias for a T nucleotide and an A nucleotide at the second and fourth positions, respectively, was observed (10). Likewise, biases for A's and T's at the central positions of the target sequence have been reported for Mo-MLV integration in vitro and in cell culture (39, 57) and for human T-lymphotropic virus type 1 insertions isolated from patient samples (11, 45). Based on these observations, both SL3-3 and HIV-1 INs target rotationally symmetric sequences, and both enzymes prefer a combination of T and A at the central positions of the targeted stretch of host DNA.
As seen in Fig. 5, a clear correlation between tumor latency and the number of insertions in the c-myc promoter was observed. In several tumors, more than two insertions in c-myc were identified (Fig. 5A). The clonality status of the insertions remains to be analyzed by Southern blotting analysis. Yet, the vast majority of the insertions may be polyclonally integrated, as demonstrated by Southern blotting analysis for the 20 insertions detected in 12 SL3-3(turbo)-induced lymphomas in female mice (Nielsen et al., unpublished). Similarly, tumors displaying single or multiple polyclonal TBLV insertions in c-myc were observed in the study by Broussard et al. (6). The lymphomas induced by wild-type and enhancer variants of SL3-3 may be generated from numerous clones in which the c-myc gene is targeted with similar frequencies, or the tumors may be made up of a smaller number of clones, in which the frequency of insertions in c-myc varies. To address this question, gene-specific PCRs have to be performed on DNA isolated from single cells.
In conclusion, the presented study demonstrates (i) a [5'-(A/C/G)TA(C/G/T)-3'] integration target site consensus sequence; (ii) nonuniform dispersion of proviral integration sites in the c-myc promoter, including repeated targeting of specific nucleotide positions; (iii) dependence of tumorigenic target site selection upon the integrating virus; (iv) correlation of faster latency periods with higher average numbers of c-myc integrations per tumor; and (v) frequent mutation of the provirus-flanking repeat structures upon error-prone gap repair.
This project was supported by the Karen Elise Jensen Foundation, the Danish Cancer Society, the Novo Nordic Foundation, and the Danish Health Sciences Research Council.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»