| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||

Gustavo Hernández-Guzmán,1
José Luis Hernández-Flores,1
Jesús Murillo,3 and
Ariel Alvarez-Morales1*
Cinvestav, IPN Unidad Irapuato, Departamento de Ingeniería Genética, Irapuato, Gto., Apdo. Postal 629, CP 36500, Mexico,1 Departamento de Ciencias Agrícolas, Universidad Nacional de Colombia, Carrera 32 Chapinero, Palmira (Valle), Colombia,2 Depto. de Producción Agraria, Universidad Pública de Navarra, Pamplona, Spain3
Received 8 December 2006/ Accepted 15 January 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
(N'-sulfodiaminophosphinyl)-ornithyl-alanyl-homoarginine] (22, 25). Phaseolotoxin is a reversible inhibitor of the enzyme ornithine carbamoyltransferase, (OCTase) (EC 2.1.3.3) (7, 17) that catalyzes the formation of citrulline from ornithine and carbamoylphosphate in the sixth step of the arginine biosynthetic pathway. In planta phaseolotoxin is readily cleaved by peptidases to release N
(N'-sulfodiaminophosphinyl)-ornithine (PSOrn), the major toxic chemical species present in diseased leaf tissue (23). The production of this phaseolotoxin is temperature dependent, being optimally produced at 18 to 20°C, while no detectable amounts of the toxin are present at 30°C (9, 21, 27). Phaseolotoxin toxin is active against not only plant OCTases but also bacterial and mammalian enzymes (6, 24, 38). This has led to the development of a rapid bioassay that measures the growth inhibition of a bacterial culture exposed to this toxin (38). However, P. syringae pv. phaseolicola is insensitive to the effect of its own toxin due to the presence of a phaseolotoxin-resistant OCTase (7, 37), which is produced under conditions leading to the synthesis of phaseolotoxin (13, 14) and is encoded by the argK gene (10, 26).
The argK gene had been previously found to be part of a gene cluster contained in a cosmid clone, pRCP17, which also included genes required for the synthesis of phaseolotoxin (28, 29). In an independent study, a second clone was isolated which contained a 25-kb insert and complemented several Tox mutants. This clone, pHK120, was shown to contain eight transcriptional units that were designated phtA through phtH. Both clones pRCP17 and pHK120 were shown to be partially overlapping; however, pHK120 lacked argK but contained sequences missing in pRCP17 (43). A gene coding an amidinotransferase, amtA, and shown to be involved in the synthesis of the homoarginine residue of the tripeptide moiety of phaseolotoxin was also found to be linked to argK (12). All this evidence suggested that most or perhaps all of the genes required for the synthesis of phaseolotoxin were part of a large gene cluster.
In P. syringae pv. phaseolicola, the nucleotide sequence of argF, a gene coding for a phaseolotoxin-sensitive OCTase present under conditions not permissive for phaseolotoxin synthesis, has very low similarity with the sequence from argK; it is so low in fact that these genes do not cross-hybridize (26, 29). The sequence of argK also presents a G+C content of 49.4%, which is below that of argF (57.3%) and the hrp genes (56 to 58%). These characteristics suggest that in P. syringae pv. phaseolicola, argK did not evolve from the housekeeping gene argF, but rather seems to have been acquired from another microorganism through horizontal gene transfer. The close linkage of argK to some of the other genes involved in the synthesis of phaseolotoxin and its essential role in protection against this toxin also suggests that not only argK but also all of the genes involved in phaseolotoxin biosynthesis and immunity may have been acquired together (10).
We decided to sequence the entire region around the argK and amtA genes in P. syringae pv. phaseolicola NPS3121 in order to determine the limits of the putative phaseolotoxin gene cluster as well as to determine the transcriptional pattern of the genes comprising this cluster. The phaseolotoxin cluster (Pht cluster) is composed of 23 genes and is flanked by insertion sequences and transposases. Furthermore, we demonstrate that the mutation of 14 of the genes within the cluster results in a Tox phenotype for 11 of them, while three mutants exhibited only low levels of toxin production. The 23 genes are organized in five transcriptional units, two monocistronic and three polycistronic, with one overlapping a larger operon. The promoter regions have been determined for each transcriptional unit, and preliminary evidence for a regulatory role for the gene product of phtL is discussed.
| MATERIALS AND METHODS |
|---|
|
|
|---|
and JM103 were grown in Luria-Bertani medium at 37°C. For phaseolotoxin production, P. syringae pv. phaseolicola was grown in M9 medium at 18°C/48 h. When required, the following supplements were added: carbenicillin, 100 µg ml1; kanamycin, 50 µg ml1; rifampin, 50 µg ml1; spectinomycin, 50 µg ml1; and X-Gluc (5-bromo-4-chloro-3-indolyl-ß-D-glucuronic acid), 30 µg ml1.
|
Construction and screening of a genomic library of P. syringae pv. phaseolicola. Chromosomal DNA from P. syringae pv. phaseolicola was obtained, and fragments of 15 to 20 kb were purified from Sau3A1 partially digested DNA, ligated into BamHI-cut EMBL3 vector, and packaged in vitro according to instructions of the supplier of commercial packaging extracts (Promega). The library was amplified in E. coli LE392, plated, and hybridized according to standard methods (33). Positive recombinant phages were amplified and replated.
DNA sequencing and analysis. The DNA inserted in the recombinant phages was digested with SalI and BamHI and subcloned into pUC19 prior to sequencing. Universal and reverse primers were used for the primary reactions, and synthesized primers were used to sequence both strands. The DNA sequence was analyzed using DNAStar and Vector NTI software. Database comparisons were made using available BLAST programs (2).
Construction of P. syringae pv. phaseolicola mutants. Restriction fragments or PCR-derived amplicons containing selected open reading frames (ORFs) from the phaseolotoxin gene cluster were cloned into pUC19 or pCR4Blunt-TOPO vectors. Mutants were obtained in these constructs by the interruption of the gene of interest with different antibiotic/uidA cassettes, followed by the replacement of the wild-type allele in the P. syringae pv. phaseolicola chromosome. Mutants with polar effects on downstream ORFs were obtained by using a 2.2-kb NruI tetracycline cassette and, in some cases, either a 3.8-kb SmaI fragment containing the uidA-aph cassette or a 1.38-kb SalI fragment containing only the aph cassette from pWM6 (20) inserted in an orientation opposite to the transcription of the ORFs. Nonpolar insertions were obtained by using the 3.8-kb SmaI fragment containing the uidA-aph cassette or the 1.38-kb SalI fragment containing only the aph cassette from pWM6 inserted in line with the transcription of the ORFs. All constructs were confirmed by restriction digests and introduced by electroporation into P. syringae pv. phaseolicola. Tetracycline or kanamycin resistance was used to select for double-recombination events. The fidelity of double recombination in all cases was confirmed by Southern blot hybridization or PCR analysis.
Construction of transcriptional uidA fusions and ß-glucuronidase (GUS) assays. Plasmid pRG960sd (40), which contains a promoterless uidA gene preceded by a Shine-Dalgarno sequence (SD-uidA) downstream of a multiple cloning site, was used to identify putative promoter sequences within intergenic regions (InRes) upstream of some of the ORFs from the phaseolotoxin cluster. These fragments were obtained as restriction fragments or by PCR using primers designed to include suitable restriction sites (data available on request) and were cloned into the BamHI-SmaI sites of pRG960sd to create pKL66, pYUDF1, and pSELF1. The orientations were determined by restriction digests, and constructs were mobilized into P. syringae pv. phaseolicola via electroporation. The intergenic region between argK and phtA (InReKA) (Fig. 1C) was cloned in both orientations to assess promoter activity towards argK and phtA.
|
Phaseolotoxin bioassays. Phaseolotoxin production by P. syringae pv. phaseolicola and mutants was assayed by the E. coli growth inhibition assay (38) as previously described (12, 35). In every case, plates containing arginine were used as controls to confirm that growth inhibition was due to phaseolotoxin.
RNA extraction and Northern blot analysis.
RNA was isolated from cultures of P. syringae pv. phaseolicola grown in M9 medium at 18 or 28°C when they reached an optical density at 600 nm of 0.8. Total RNA was extracted from cells by using TRIzol reagent as recommended by the manufacturer (Invitrogen). Genomic DNA was removed by digestion with RNase-free DNase (Roche). Samples of total RNA (20 µg) were denatured by treatment with formamide and separated by electrophoresis using 1.3% denaturing agarose gels. The RNA was transferred to Hybond N+ nylon membranes (Amersham) and cross-linked by exposure to UV radiation. The hybridization of the nylon membranes was performed using Church buffer (33) and different DNA probes for each ORF. The probes were labeled with [
-32P]dCTP by using the Rediprime random primer labeling kit (Amersham). The hybridization was carried out overnight at 65°C. The membranes were washed with 2x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate)-0.1% sodium dodecyl sulfate for 3 min at room temperature, followed by a wash with 1x SSC-0.1% sodium dodecyl sulfate for 3 min at 65°C. The membranes were exposed in a Storm 860 apparatus, and the analysis was made using ImageQuant 1.1 software.
Reverse transcription-PCR analysis. DNA-free RNA was obtained from cultures grown in M9 medium for 48 h at 18 and 28°C. The RNA was checked for integrity in an agarose gel and used for reverse transcription (RT) and PCR using the SuperScript one-step kit (Invitrogen). A list of the primers used in this analysis is available on request. Controls used for each set of primers were (i) PCR without the reverse transcription step to verify the absence of DNA, (ii) RT-PCRs performed without RNA templates to detect any contaminating DNA/RNA, (iii) PCRs performed using chromosomal DNA as a template to ensure primer fidelity, and (iv) amplification of a portion of the 23S rRNA operon using suitable primers as an internal control of the reaction. A portion of the 23S rRNA operon was amplified using suitable primers. The RT reaction was performed at 50°C for 40 min, followed by PCR amplification at 94°C for 2 min for 1 cycle; 94°C for 35 s, 58°C for 30 s, and 72°C for 2.0 min for 25 cycles; and 72°C for 15 min for 1 cycle.
Mapping the site of transcription initiation.
The site of transcription initiation for each operon was obtained by sequencing the 5' end of the different mRNA transcripts using a 5' circular first-strand cDNA-mediated RACE (rapid amplification of cDNA ends) method (5' cRACE) (18) as follows. First, the synthesis of single-stranded cDNA was performed. Total DNA-free RNA was obtained from cultures grown in M9 medium for 48 h at 18°C. One microgram of this RNA was used as a template to synthesize the first-strand cDNA of the mRNA by using a cDNA synthesis kit (Invitrogen) with a gene-specific oligonucleotide primer designed to anneal within the coding region of the gene. Second, a poly(A) tail was added to the 3' end of the single-stranded cDNA. After hydrolysis of the template mRNA with 0.5 N NaOH (final concentration) at 37°C for 10 min, the single-stranded cDNA was precipitated with ethanol and redissolved in 10 µl of a reaction mixture containing 25% polyethylene glycol (molecular weight,
8,000), 10 mM Tris-HCl (pH 8.0), 1 mM hexamine cobalt chloride, 10 mM MgCl2, 0.01 mM ATP, 1 µg/ml bovine serum albumin, 10 U T4 RNA ligase (New England Biolabs), and a commercial 30mer poly(A) oligonucleotide previously phosphorylated with T4 polynucleotide kinase (New England Biolabs). This ligation mixture was incubated at 22°C for 16 h. Third, the first PCR amplification was performed. An aliquot of the previous ligation mixture was directly used as a template for the first PCR amplification. The PCR was set up in a total volume of 50 µl containing 0.5 µl of the cDNA ligation mixture, a poly(T) oligonucleotide, and a gene-specific primer. The PCR was carried out for 30 cycles with the following conditions: 94°C for 30 s denaturation, 58°C for 30 s annealing, and 72°C for a 1.5-min extension. Last, nested PCR was performed. The resulting PCR product was diluted 1,000-fold with sterile H2O, and a 2-µl aliquot was used as template for a second nested PCR amplification with a 30mer commercial poly(T) oligonucleotide and a gene-specific nested primer under the same conditions. The amplification products were cloned into a TOPO-TA cloning vector (Invitrogen) and sequenced.
| RESULTS |
|---|
|
|
|---|
clones which allowed us to sequence the regions around both of these genes. With this strategy, we were able to sequence and assemble around 25 kb before we had access to the preliminary release of the complete sequence of P. syringae pv. phaseolicola 1448A (GenBank accession no. CP000058 to CP000060). Once we had this information, the gaps in the sequence of P. syringae pv. phaseolicola NPS3121 were covered by sequencing amplicons obtained through PCR generated with primers designed from both the partial sequence of NPS3121 and the reported sequence from 1448A. A total of 30,245 bp were sequenced on both strands (GenBank DQ141263). When we aligned the sequence we obtained for NPS3121 to the corresponding region reported for strain 1448A, we found only 35 nucleotide differences, all of them occurring after position 14694 in the sequence of NPS3121. Among these changes, there were 15 bases not found in 1448A, including a stretch of 9 contiguous nucleotides, 4 bases missing in NPS3121 with respect to 1448A, and 16 single-point changes corresponding to 6 transitions and 10 transversions. Of these, the most significant change corresponds to the 9-nucleotide stretch found in NPS3121 between positions 18748 and 18756 since this introduces a stop codon in the region corresponding to locus PSPPH4306 of 1448A, and as a consequence, in strain NPS3121, we found two separate and distinct ORFs (phtM and phtN) at this locus.
The computational analysis of the 30,245-bp region of P. syringae pv. phaseolicola NPS3121 suggested the presence of 27 ORFs in this region (Fig. 1A). Of these, three ORFs at the 5' end correspond to a sir2 homolog, a "hypothetical protein," and an ISPsy25 transposase. At the 3' end, between positions 28997 and 29324, we found a 328-bp region whose translation showed 59% identity to Sir2 from P. syringae pv. tomato DC3000, 55% to that from Pseudomonas fluorescens Pf-5, and 54% to that from Pseudomonas entomophila L48 among others. This truncated region of similarity is linked through a 127-bp region to a region with homology to ISPsy5/ISPpu13 transposases. Between these ends, we defined 23 ORFs, each one preceded by a putative ribosome-binding sequence. In this work, we refer to this group of ORFs as the phaseolotoxin (Pht) cluster as previously proposed (11). A BLASTX analysis of the Pht cluster returned the same information reported for strain 1448A (16). The unambiguous genes identified at the moment are argK (coding for the phaseolotoxin-resistant ornithyl carbamoyl transferase [10, 26]), amtA (coding for an L-arginine:lysine amidinotransferase [12]), and the des gene (desI) (coding for a fatty acid desaturase [11]).
Phenotype of mutants. Mutants with polar effects on downstream genes were obtained for phtA, -B, -E, -F, -H, -K, -L, -T, amtA, and sir2, and these were unable to produce phaseolotoxin as shown by the growth inhibition assay, the only exception being the mutant on the sir2 gene located outside the Pht cluster which produces phaseolotoxin at a level comparable to that of the wild-type strain. Nonpolar mutants were obtained for desI and the pht genes phtK, -L, -M, -O, -Q, and -T. In this case, mutants in phtK, -L, -M, and -Q showed a clear Tox phenotype, whereas desI, phtO, and phtT were still capable of producing phaseolotoxin, although at a very low level (Fig. 2).
|
In order to identify the promoters responsible for the observed expression, InRes of the Pht cluster were cloned into a vector containing the promoterless uidA gene preceded by a SD sequence. These constructs were electroporated into P. syringae pv. phaseolicola to assess promoter activity in trans. Also, the nonpolar insertions in desI and the pht genes phtK, -L, -M, -O, and -Q, were obtained with a construct comprising the SD-uidA gene followed by the aph gene, allowing for the promoter of the interrupted gene/operon to drive uidA expression in cis. Results from these experiments are depicted graphically in Fig. 1C; nevertheless, we are aware that the levels of expression do not necessarily correlate to levels in vivo since the reporter in trans occurs in multiple copies (40), whereas in the fusions in cis, a direct or indirect effect on the levels of transcription due to the mutation itself could not be excluded.
Expression was unambiguously detected for the InRe between argK and phtA, both in the direction of argK (InReKAa
) and divergently into phtA (InReKAb
), suggesting the presence of two distinct promoters. Promoter activity was also detected within the InRe between phtC and phtD (InReCD), but not within the 68-bp intergenic region between phtA and phtB (InReAB), indicating that phtABC comprise an operon and that phtD is the first gene of a different transcriptional unit. The 141-bp region between desI and phtJ and the 129-bp region between phtT and phtU were also unable to drive uidA expression, indicating the lack of promoter sequences (data not shown).
With respect to the uidA fusion in cis, activity was detected within desI. Considering that the largest InRe found within the phtD-desI fragment is only 13 bp (within phtD and phtE) and the lack of promoter activity in the InRe between desI and phtJ, the result strongly suggests that phtD, -E, -F, -G, -H, desI, and phtJ are part of a single operon, which is also supported by the expression pattern (observed through Northern blot analysis) that suggests a similar pattern for both desI and phtJ. The expression of phtK as assayed with Northern blotting suggests that it could be part of the same operon as that for phtJ.
With respect to genes downstream of phtL, we observed GUS activity from the uidA fusion to phtM, phtO, and phtQ. The largest intergenic space found in the region from phtM to phtQ is 116 bp between phtP and amtA. However, our previous primer extension analyses of the region upstream of amtA consistently suggested that this gene was transcribed as part of a larger transcript (data not shown), therefore it seems reasonable to assume that phtM, -N, -O, -P, amtA, and phtQ are all part of a single operon. To analyze the operon structure downstream of phtQ, confirm the structure suggested by the expression experiments, and solve the uncertainties remaining, we conducted reverse transcription-PCR experiments.
Reverse transcription-PCR analysis. We conducted experiments in which we tried to amplify specific fragments derived from cDNA of contiguous ORFs to investigate whether these were being transcribed from the same mRNA. In every case, we included as a control the amplification within a single pht ORF besides the amplification of the 23S rRNA operon. The results are shown graphically in Fig. 3A. Primers spanning argK and phtA did not produce any amplicon, whereas amplification from argK produced the expected fragment in the same reaction, which confirms that argK and phtA belong to different transcription units.
|
|
There are 14 bp between the stop codon of phtL and the predicted ATG of phtM; however, no amplification product was obtained using primers spanning both genes (Fig. 3B). The expression of phtL was detected at temperatures permissive, 18°C, and nonpermissive, 28°C, for phaseolotoxin production. This pattern of expression was observed in Northern blot experiments (P13) (Fig. 1B), with the uidA fusion in cis (phtL) (Fig. 1C), and in the RT-PCR experiments (phtLb) (Fig. 3B). This expression pattern contrasts with that of phtM, which in similar experiments, is always negligible at 28°C (Fig. 1B and C and 3B). This result indicates that phtL and phtM cannot be transcribed from the same mRNA. Therefore, to explain the apparent lack of space to accommodate the promoter driving the expression of phtM, we assumed that the coding sequence could begin 51 or 132 bp downstream of the postulated ATG, at a second or third ATG codon in the same reading frame, leaving up to 146 bp to accommodate the promoter for phtM and having phtL as a monocistronic transcriptional unit.
In agreement with the fusion data, amplification was obtained from phtMN, phtO, amtA, phtQRST, and phtTU, indicating that most likely all of the genes from phtM to at least phtU comprise a single transcriptional unit (Fig. 3A).
Taken together, all these results suggest that there are five transcriptional units (Fig. 5), two of which are monocistronic, argK and phtL. Additionally, there is a promoter driving the expression of a large operon that includes 11 genes, from phtA to phtK, with an internal promoter immediately downstream of phtC capable of driving expression of phtD to phtK, and a further transcriptional unit that includes 10 genes, phtM to phtV.
|
70) promoter with appropriate 10 and 35 regions (41). The initiation of transcription that was determined for phtA also shows a well-conserved 10 and 35 region characteristic of a Pribnow-type promoter, and these divergent promoters show two A-rich regions on both strands that could be implicated in polymerase binding (Fig. 6A) (1).
|
factors (41) or of the HrpL
factor (8). However the promoter region of phtD and phtM share at least six obvious conserved regions within 60 bp upstream of the site of transcription initiation, suggesting that both promoters are under a common mechanism of transcriptional regulation. None of the promoters identified show, within the predicted polymerase binding region, the TRR elements previously postulated to be involved in thermoregulation of phaseolotoxin synthesis (30, 31), although these TRR sequences could possibly be located outside this region to exert their effects. Transcriptional analysis of a phtL mutant. It was interesting to note that the expression of phtL occurred at 28°C at a level significantly above that observed for the other pht genes. The observation that transcripts of phtL were present at 28°C, a temperature nonpermissive for phaseolotoxin synthesis, suggested that the gene product of phtL could have a role in regulation, such as being a temperature sensor and/or signal transducer. Consistent with this role, the protein from this gene should be present at 28°C and a mutation in this gene should result in a Tox phenotype because perhaps the temperature change would not be detected or this information would not be relayed to the cognate effector. To investigate this further, we analyzed the transcription pattern of the Pht cluster in a phtL insertion mutant through RT-PCR. Consistent with the idea that the PhtL product may be involved in regulation, all of the transcripts of the phtMV operon were clearly absent at 18°C, contrasting with the situation observed in the wild-type strain (Fig. 7).
|
| DISCUSSION |
|---|
|
|
|---|
The fact that insertion mutants in 11 different genes of the cluster, from phtA to phtT, result in Tox phenotypes clearly indicates that all genes within this region encode proteins that are required at any of the different stages of phaseolotoxin production, such as synthesis, transport, and regulation. With respect to the transcriptional organization of the Pht cluster, argK has been reported to be a monocistronic transcriptional unit (10, 26) and our results show that there are four other regions where distinct promoter activity was observed by the use of fusions to a promoterless uidA gene either in cis or in trans. The initiation of the transcription for each promoter has been determined, indicating that the 23 genes of the Pht cluster are organized in five transcriptional units, two being monocistronic, argK and phtL, with a large complex polycistronic operon comprising 11 genes, from phtA to phtK, harboring an internal promoter upstream of phtD and another large operon comprising genes from phtM to phtV.
It is interesting to note that the largest operon, comprising 11 genes from phtA to phtK, presents an internal promoter, a feature observed in operons harboring genes for the production of secondary metabolites, such as antibiotics (5, 39), or for the expression of gene regulators, such as the ntrBC two-component regulatory system involved in regulating different aspects of nitrogen metabolism in gram-negative bacteria (19).
We have shown that the promoter driving the expression of phtABC is a Pribnow-type promoter very similar to the divergent promoter of argK and that both share an interpromoter region with distinct sequence features that suggest coordinate regulation. It is interesting to note that while trying to obtain evidence for a regulatory role for the phtL gene product, we observed that a null phtL mutant completely abolishes the transcription of the phtM operon, consistent with a regulatory role for the phtL gene product, with no effect on any of the other transcriptional units (Fig. 7). However, our findings indicate that both the phtD and pthM operons are transcribed from very similar promoters, suggesting a common mechanism of regulation. The fact that transcripts from the phtD operon are present in the phtL mutant may indicate that transcription is proceeding from the phtA promoter generating a transcript from phtA to phtK and overriding the effect of the phtL mutation. This in turn suggests that PhtL does not have an effect on the regulation of the argK or phtA promoters.
Therefore, we propose the Pht cluster to be composed of five transcriptional units, which we have named after the gene immediately downstream of the promoter region and presented in the structure depicted in Fig. 5.
We are presently investigating how these promoters are being regulated and the particular sequence of the phtL promoter as well as which type of
factor is involved in recognition of the phtD, -M, and -L promoters. Also, we are currently investigating the precise role of the phtL gene product. The identity and role of most of the gene products of the Pht cluster remain to be determined as it is still uncertain which genes participate in regulation and in the biosynthesis of the inorganic moiety and the tripeptide composing phaseolotoxin.
As we have shown conclusively that all the genes within the Pht cluster (with the sole exception of phtV) are expressed and have a role in phaseolotoxin production, we propose the use of the nomenclature used in this article for further reference to these genes until more information is obtained for the proteins they encode and specific names can be assigned.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
Published ahead of print on 19 January 2007. ![]()
Present address: Facultad de Medicina Veterinaria y Zootecnia, Universidad Michoacana de San Nicolás de Hidalgo, Morelia, Mexico. ![]()
| REFERENCES |
|---|
|
|
|---|
subunit. Proc. Natl. Acad. Sci. USA 95:14652-14657.
-(N'-sulfodiaminophosphinyl)-L-ornithine, a true transition state analogue? J. Biol. Chem. 275:20012-20019.
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Appl. Environ. Microbiol. | Infect. Immun. | Eukaryot. Cell |
|---|---|---|
| Mol. Cell. Biol. | J. Virol. | Microbiol. Mol. Biol. Rev. |
| ALL ASM JOURNALS |