Previous Article | Next Article ![]()
Molecular and Cellular Biology, July 2003, p. 4439-4448, Vol. 23, No. 13
0270-7306/03/$08.00+0 DOI: 10.1128/MCB.23.13.4439-4448.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
and Jun Ma*
Graduate Program in Molecular and Developmental Biology, Division of Developmental Biology, Children's Hospital Research Foundation, University of Cincinnati College of Medicine, Cincinnati, Ohio 45229
Received 6 February 2003/ Returned for modification 25 March 2003/ Accepted 1 April 2003
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Previous studies have shown that Bcd binds DNA in a highly cooperative manner (19). Such cooperativity is likely to play an important role in facilitating the sharp on/off switches of target genes, such as hb, in response to the Bcd gradient in embryos. Cooperative DNA binding by Bcd is achieved through interactions between Bcd molecules (19), primarily relying on protein sequences outside the homeodomain (35), including the amino-terminal domain (residues 1 to 91). In particular, the homeodomain alone, which fails to self associate efficiently in vitro, has little cooperative DNA binding activity (35, 37). Experiments that used an altered-specificity mutant of another homeodomain protein, Ftz(Q50K), further suggest that protein-protein interaction mediated by Bcd sequences outside its homeodomain represents an important mechanism for selecting specific target genes for activation (36, 37).
In addition to protein-protein interaction and cooperative DNA binding, the evolutionarily conserved amino-terminal domain of Bcd also provides a self-inhibitory function (38). In particular, Bcd derivatives lacking the amino-terminal domain exhibit a greatly increased ability to activate a hb-CAT reporter gene containing the Bcd-responsive hb enhancer element in Drosophila S2 cells. A mutant protein, Bcd(A52-56), which has the critical residues 52 to 56 mutated to alanines, exhibits a dominant effect in wild-type embryos, causing a posterior shift of the fatemap and patterning defects (38). The expression of Bcd target genes and other segmentation genes is also significantly altered in embryos containing Bcd(A52-56) (38; unpublished data). These results suggest that the self-inhibitory function provided by the amino-terminal domain of Bcd is essential to embryogenesis.
In this report we demonstrate that, depending on enhancer sequences, the amino-terminal domain of Bcd is preferentially utilized for one of two distinct functions: cooperativity or self inhibition. Specifically, while cooperative DNA binding to the kni enhancer element is highly dependent on the amino-terminal domain of Bcd, this domain is largely dispensable for cooperative binding to the hb enhancer element. In contrast, while activation from the hb enhancer element is highly sensitive to the self-inhibitory function, mutations affecting this function only cause limited effects on the kni enhancer element. Our further analysis of kni enhancer derivatives suggests that enhancer structure plays a critical role in determining the characteristics in responding to Bcd. We propose that a coordinated execution of the two functions provided by the amino-terminal domaincooperativity and self inhibitionhelps define distinctive responses to the Bcd gradient.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
-32P]dCTP. In our experiments the kni probe was isolated as an XbaI-SacI fragment from pCZ72, and the hb probe was isolated as an XhoI-XbaI fragment from pMAX1 (37). The kni enhancer derivatives kni(3R), 3HH, and 3TT were isolated as EcoRI-NotI fragments from pDF501, pDF510, and pDF511, respectively. Probes for single pairs of head-to-head (HH) and tail-to-tail (TT) sites were isolated as NotI-BamHI fragments from pDF510 and pDF511, respectively. The probes were diluted to a final concentration of 1.6 x 10-11 M. The experimental procedures and conditions for gel shift assays were described previously (37). Wild-type Bcd and its derivatives used in this assay were expressed in vitro by using TnT quick coupled transcription/translation systems (Promega). The active protein concentrations were estimated by using high concentrations (4 x 10-8 M) of a 32P-labeled consensus Bcd binding site, A1 (7). The same amounts of active Bcd proteins were used in gel shift experiments. The Molecular Dynamics PhosphorImager ImageQuant program was used for quantitative analysis. Yeast strain and ß-galactosidase liquid assays. The yeast strain used in this study is CY26 containing the integrated hb-lacZ or kni-lacZ reporter gene (37). The Bcd-VP16 effector plasmids were introduced into yeast cells by using the lithium acetate method, and three independent transformants were assayed for ß-galactosidase activities as described previously (37).
Transient transfection experiments. Drosophila S2 cells (Invitrogen) were grown at 25°C in Schneider's Drosophila medium (GIBCO) supplemented with 10% fetal bovine serum (GIBCO). The cells were seeded in 60-mm-diameter tissue culture plates at 4 x 106/plate 24 h prior to transfection. Transfection was performed by the calcium phosphate transfection system following the protocol from GIBCO. Each plate was transfected with 1 µg of reporter plasmid, 1 µg of Copia-lacZ internal control plasmid, and the indicated amount of effector plasmids. Salmon sperm carrier DNA (GIBCO) was used to bring the total amount of DNA to 10 µg for each transfection. The cells were harvested 48 h later, and whole-cell lysates were prepared by using the lysis buffer (50 mM HEPES [pH 7.5], 300 mM NaCl, 0.5% Triton X-100, 1 mM dithiothreitol, 1 mM phenylmethylsulfonyl fluoride) following the freeze-thaw procedure (3). The transfection efficiency was determined by monitoring the ß-galactosidase activity, and the amounts of lysates used in chloramphenicol acetyltransferase (CAT) and Western blotting assays were normalized accordingly. CAT assays were performed as described previously by using three independently transfected samples for each experiment (3). For Western blotting, cell lysates were separated on 10% polyacrylamide sodium dodecyl sulfate gels and were transferred to nitrocellulose membrane (Bio-Rad), which was blotted with anti-HA monoclonal antibody (HA-11; 1:500 final dilution; Babco).
| RESULTS |
|---|
|
|
|---|
|
The amino-terminal domain of Bcd is required for cooperative binding to the kni but not hb enhancer. The positive effect of the amino-terminal domain of Bcd on kni-CAT expression was observed predominantly at low Bcd concentrations (Fig. 1D), suggesting that this domain may play an important role in facilitating the protein to bind the kni enhancer. To test this possibility, we conducted gel shift assays with both full-length Bcd and Bcd(92-489) on hb and kni enhancers (Fig. 2B for hb and C for kni). Our results show that wild-type Bcd can bind to both enhancer elements cooperatively (Fig. 2B and C, lanes 1 to 4). In particular, as the protein concentration increased, the monomeric complex containing a single Bcd molecule (marked with an arrowhead) was depleted to form oligomeric complexes. Interestingly, contrary to the simple affinity threshold model (see Discussion for further details), the kni enhancer did not exhibit a higher affinity for Bcd than the hb enhancer did (Fig. 2B and C, lanes 1 to 4).
|
We also analyzed systematic Bcd deletion derivatives in yeast cells (Fig. 3). In this experiment, Bcd derivatives fused to the strong activation domain from VP16 (27) were compared for their abilities to activate integrated hb-lacZ and kni-lacZ reporter genes in yeast cells. Both full-length Bcd-VP16 and Bcd(1-246)-VP16 fusion proteins activated transcription from the hb-lacZ and kni-lacZ reporter genes effectively (lines 5 and 12). Therefore, all deletions were derived from the truncated protein Bcd(1-246)-VP16, with its activities on each reporter set to 100. Figure 3 shows that a derivative lacking the amino-terminal domain, Bcd(92-246)-VP16, failed to activate transcription from the kni-lacZ reporter while remaining active on hb-lacZ (line 4). In contrast, Bcd(1-151)-VP16, which lacks the Bcd sequence on the C-terminal side of its homeodomain, had a more significant reduction in activity on hb-lacZ than on kni-lacZ (line 9). The data shown in Fig. 3 further suggest that protein sequences immediately flanking the homeodomain play positive roles on kni-lacZ (compare lines 4 and 8) and hb-lacZ (compare lines 9 and 10; also compare lines 2 through 4). Interestingly, Bcd(92-151)-VP16, which contains only the homeodomain of Bcd fused to VP16, failed to efficiently activate from both reporters (line 1). Since the self-inhibitory function conferred by the amino-terminal domain of Bcd is not operative in yeast cells (38 and unpublished data) and, moreover, since all the derivatives contain the strong activation domain VP16, these results are consistent with the suggestion that different Bcd sequences outside its homeodomain play important but differential roles in cooperative binding to hb and kni enhancers.
|
To further determine whether the self-inhibitory function provided by the amino-terminal domain of Bcd is differentially implemented on different enhancers, we took advantage of these point mutation derivatives. We performed transient transfection assays by using high protein concentrations to specifically compare the effectiveness of the self-inhibitory function on the hb-CAT and kni-CAT reporters. Figure 4 shows that, while Bcd(A52-56) was 25 times more active than wild-type Bcd on hb-CAT, the difference was only 6.25-fold on kni-CAT (lane 3). In addition, while Bcd(A57-61) was virtually inactive on hb-CAT, its activity on kni-CAT was reduced only by 40% (lane 4). As shown in Fig. 1C and D, at high protein concentrations the activity difference between full-length Bcd and the N-terminally truncated derivative Bcd(92-489) was also much greater on hb-CAT (>40-fold) than on kni-CAT (<4-fold) (Fig. 4, lane 2). Together these results demonstrate that the self-inhibitory function is more robust on hb-CAT than on kni-CAT. Interestingly, Bcd(92-489) was more active on hb-CAT than Bcd(A52-56), but the opposite was true on kni-CAT; this difference is consistent with the finding that the amino-terminal domain is preferentially required for cooperative binding to the kni enhancer (Fig. 2).
|
|
On enhancers that contain multiple Bcd binding sites, the Bcd-DNA complexes were poorly resolved in gel shift assays (Fig. 2 and 5). In contrast, the monomeric and dimeric Bcd complexes could be readily resolved on DNA probes containing only two Bcd binding sites (Fig. 6). This permitted a more accurate and quantitative analysis of the role of the amino-terminal domain in cooperative DNA binding. Specifically, we determined Hill coefficient values for full-length Bcd and Bcd(92-489) on single pairs of closely spaced HH (Fig. 6A and C) and TT (Fig. 6B and D) binding sites. Wild-type Bcd bound to both HH and TT pairs cooperatively, with an estimated Hill coefficient of 1.88 ± 0.28 and 1.48 ± 0.25, respectively. In contrast, Bcd(92-489) had a significantly reduced but measurable cooperativity on the HH pair, with an estimated Hill coefficient of 1.40 ± 0.30. This same Bcd derivative became completely noncooperative on the TT pair, with an estimated Hill coefficient of 0.66 ± 0.10, which is similar to the value (0.8) on a single Bcd binding site. These results further demonstrate that the amino-terminal domain plays an important role on closely spaced symmetric pairs of Bcd binding sites, particularly the TT pair.
|
The results shown in Fig. 7 demonstrate that both 3HH-CAT and 3TT-CAT reporters became more sensitive to the self-inhibitory function of Bcd. In particular, Bcd(A57-61) had much lower activities on both 3TT-CAT and 3HH-CAT than on kni-CAT (compare lane 4 of Fig. 7 with Fig. 4). In addition, both Bcd(A52-56) and Bcd(92-489) exhibited higher activities on 3TT-CAT and 3HH-CAT than on the kni-CAT reporter (Fig. 7, lanes 2 and 3). Our results also reveal a difference between 3HH-CAT and 3TT-CAT in their responses to the A52-56 mutation and removal of the entire amino-terminal domain (Fig. 7, lanes 2 and 3). This difference is consistent with the finding that cooperative binding to the TT pair is more dependent on the amino-terminal domain than to the HH pair (Fig. 6). Finally, although the responses to the mutations affecting the self-inhibitory function on the kni(3R)-CAT reporter were not dramatically altered compared to those of kni-CAT, all activities were significantly reduced on knic(3R)-CAT (27%) (see also reference 6). These results are consistent with the finding that Bcd fails to bind cooperatively to kni(3R) (Fig. 5), demonstrating that reversion of the orientations of Bcd binding sites in the kni enhancer alters its characteristics in responding to Bcd.
|
| DISCUSSION |
|---|
|
|
|---|
Both cooperative DNA binding and self inhibition involve protein-protein interactions (19, 35, 38), suggesting that the amino-terminal domain of Bcd is engaged in distinct protein-protein interaction events on different enhancers. Our results suggest that the engagement of the amino-terminal domain in cooperative DNA binding on the kni enhancer likely hinders its availability for self inhibition. Tight spacing between the Bcd binding sites may also reduce the effectiveness of the self-inhibitory function. We presently do not know precisely how the amino-terminal domain of Bcd facilitates cooperative DNA binding. While this domain is strictly required for a single pair of closely spaced TT sites, its role on a single pair of closely spaced HH sites, though important, is somewhat less critical (Fig. 6). It is possible that the amino-terminal domains of two Bcd molecules are engaged in a direct homodimeric interaction. It is also possible that this domain of one Bcd molecule interacts with a different domain of another Bcd molecule on DNA. Our previous protein-protein interaction experiments suggest that the amino-terminal domain alone does not self associate efficiently without the homeodomain (35). In this context it is relevant that cooperative DNA binding can be achieved on the hb enhancer even without the amino-terminal domain. Bcd contains multiple self-association domains (35) which play differential roles on distinct enhancers (Fig. 3). Regardless, our findings suggest that Bcd can accommodate and coordinate multiple interaction events that involve different parts of the protein, particularly its amino-terminal domain. Interestingly, the amino-terminal domain of Bcd also contains a motif that can interact with the translation initiation factor eIF4E and participates in translation repression of caudal (cad) mRNA (23), further underscoring the importance of this evolutionarily conserved domain.
One of the fundamental biological questions concerning the action of Bcd is how different target genes respond to distinct thresholds of the Bcd gradient in embryos. It was previously proposed that the affinity of Bcd binding sites within an enhancer dictates the concentration of Bcd required for activation in embryos (12, 14). However, this simple affinity threshold model has been challenged recently (36). In particular, the assumption of this modelthe consensus A-type sites have higher Bcd affinity than the nonconsensus X-type siteshas been proven incorrect. Biochemical studies have shown that Bcd binds to a consensus A1 site and a nonconsensus X1 site with comparable affinities (7, 19). In addition, the kni enhancer, despite its ability to respond to lower Bcd concentrations, does not exhibit any higher Bcd affinity than the hb enhancer (Fig. 2). A more recent study also suggests that the concentration of Bcd in embryos, though important, is unlikely to be the only determinant for defining the target gene responses to the Bcd gradient (17). Our characterizations of the self-inhibitory function in embryos suggest that the activation potential of Bcd may play an important role in controlling target gene expression. In particular, Bcd(A52-56), which has a stronger activation potential, can exert a dominant effect in wild-type embryos, causing a posterior shift of the target genes (38). These results suggest that a stronger Bcd can activate transcription at a lower concentration (or with fewer DNA-bound molecules) in embryos, a concept that has been proposed previously for various activators (2, 5, 18, 20, 29). We propose that the muted self-inhibitory function of Bcd on the kni enhancer contributes to the protein's ability to activate transcription at low concentrations.
Many homeodomain proteins have been shown to cooperate specifically with their DNA binding partners (16, 21, 22, 31-33). Although no such partner proteins have been identified to specifically help Bcd to select its target genes, this possibility cannot be formally ruled out at this time. If such proteins do exist, they could potentially play important roles in facilitating Bcd to act distinctively on different enhancers. In this context, it is interesting that the Bcd-responsive kni enhancer element used in this study contains binding sites for the Caudal (Cad) protein, although shorter kni enhancer fragments without these Cad binding sites can similarly respond to low concentrations of Bcd in embryos (26). The fact that the experiments described in this report use only Bcd-responsive enhancer elements in Drosophila cells indicates that Bcd may have an intrinsic ability to activate different target genes distinctively. These studies demonstrate that even relatively simple enhancers can have profound effects in influencing activator behaviors. Understanding precisely how complex enhancers orchestrate in vivo the individual actions of different transcription factors represents a rewarding challenge for molecular biologists.
| ACKNOWLEDGMENTS |
|---|
This work was supported in part by NIH and AHA grants (to J.M.).
| FOOTNOTES |
|---|
Present address: Institute of Cellular and Molecular Biology, University of Texas, Austin, TX 78712. ![]()
| REFERENCES |
|---|
|
|
|---|
2. Arnosti, D. N., S. Barolo, M. Levine, and S. Small. 1996. The eve stripe 2 enhancer employs multiple modes of transcriptional synergy. Development 122:205-214.[Abstract]
3. Ausubel, F., R. Brent, R. Kingston, D. Moore, J. Seidman, J. Smith, and K. Struhl. 1994. Current protocols in molecular biology. John Wiley & Sons, Inc., New York, N.Y.
4. Berleth, T., M. Burri, G. Thoma, D. Bopp, S. Richstein, G. Frigerio, M. Noll, and C. Nüsslein-Volhard. 1988. The role of localization of bicoid RNA in organizing the anterior pattern of the Drosophila embryo. EMBO J. 7:1749-1756.[Medline]
5. Bunker, C. A., and R. E. Kingston. 1996. Activation domain-mediated enhancement of activator binding to chromatin in mammalian cells. Proc. Natl. Acad. Sci. USA 93:10820-10825.
6. Burz, D. S., R. Pivera-Pomar, H. Jackle, and S. D. Hanes. 1998. Cooperative DNA-binding by Bicoid provides a mechanism for threshold-dependent gene activation in the Drosophila embryo. EMBO J. 17:5998-6009.[CrossRef][Medline]
7. Dave, V., C. Zhao, F. Yang, C. S. Tung, and J. Ma. 2000. Reprogrammable recognition codes in Bicoid homeodomain-DNA interaction. Mol. Cell. Biol. 20:7673-7684.
8. Driever, W. 1992. The Bicoid morphogen: concentration dependent transcriptional activation of zygotic target genes during early Drosophila development, p. 1221-1250. In S. L. McKnight and K. Yamamoto (ed.), Transcriptional regulation. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, N.Y.
9. Driever, W., J. Ma, C. Nusslein-Volhard, and M. Ptashne. 1989. Rescue of bicoid mutant Drosophila embryos by Bicoid fusion proteins containing heterologous activating sequences. Nature 342:149-154.[CrossRef][Medline]
10. Driever, W., and C. Nüsslein-Volhard. 1989. Bicoid protein is a positive regulator of hunchback transcription in the early Drosophila embryo. Nature 337:138-143.[CrossRef][Medline]
11. Driever, W., and C. Nüsslein-Volhard. 1988. A gradient of bicoid protein in Drosophila embryos. Cell 54:83-93.[CrossRef][Medline]
12. Driever, W., G. Thoma, and C. Nüsslein-Volhard. 1989. Determination of spatial domains of zygotic gene expression in the Drosophila embryo by the affinity of binding site for the bicoid morphogen. Nature 340:363-367.[CrossRef][Medline]
13. Frohnhöfer, H. G., and C. Nüsslein-Volhard. 1986. Organization of anterior pattern in the Drosophila embryo by the maternal gene bicoid. Nature 324:120-125.[CrossRef]
14. Gao, Q., and R. Finkelstein. 1998. Targeting gene expression to the head: the Drosophila orthodenticle gene is a direct target of the Bicoid morphogen. Development 125:4185-4193.[Abstract]
15. Gao, Q., Y. Wang, and R. Finkelstein. 1996. Orthodenticle regulation during embryonic head development in Drosophila. Mech. Dev. 56:3-15.[CrossRef][Medline]
16. Guichet, A., J. W. Copeland, M. Erdelyi, D. Hlousek, P. Zavorszky, J. Ho, S. Brown, A. Percival-Smith, H. M. Krause, and A. Ephrussi. 1997. The nuclear receptor homologue Ftz-F1 and the homeodomain protein Ftz are mutually dependent cofactors. Nature 385:548-552.[CrossRef][Medline]
17. Houchmandzadeh, B., E. Wieschaus, and S. Leibler. 2002. Establishment of developmental precision and proportions in the early Drosophila embryo. Nature 415:798-802.[Medline]
18. Lehman, A. M., K. B. Ellwood, B. E. Middleton, and M. Carey. 1998. Compensatory energetic relationships between upstream activators and the RNA polymerase II general transcription machinery. J. Biol. Chem. 273:932-939.
19. Ma, X., D. Yuan, K. Diepold, T. Scarborough, and J. Ma. 1996. The Drosophila morphogenetic protein Bicoid binds DNA cooperatively. Development 122:1195-1206.[Abstract]
20. Ma, X., D. Yuan, T. Scarborough, and J. Ma. 1999. Contributions to gene activation by multiple functions of Bicoid. Biochem. J. 338:447-455.
21. Mak, A., and A. D. Johnson. 1993. The carboxy-terminal tail of the homeodomain protein
2 is required for function with a second homeodomain protein. Genes Dev. 7:1862-1870.
22. Mann, R. S., and S. K. Chan. 1996. Extra specificity from extradenticle: the partnership between HOX and PBX/EXD homeodomain proteins. Trends Genet. 12:258-262.[CrossRef][Medline]
23. Niessing, D., S. Blanke, and H. Jackle. 2002. Bicoid associates with the 5'-cap-bound complex of caudal mRNA and represses translation. Genes Dev. 16:2576-2582.
24. Nusslein-Volhard, C. 1991. Determination of the embryonic axis of Drosophila. Development 1(Suppl.):1-10.
25. Rivera-Pomar, R., and H. Jackle. 1996. From gradients to stripes in Drosophila embryogenesis: filling in the gaps. Trends Genet. 12:478-483.[CrossRef][Medline]
26. Rivera-Pomar, R., X. Lu, H. Taubert, N. Perrimon, and H. Jackle. 1995. Activation of posterior gap gene expression in the Drosophila blastoderm. Nature 376:253-256.[CrossRef][Medline]
27. Sadowski, I., J. Ma, S. Triezenberg, and M. Ptashne. 1988. GAL4-VP16 is an unusually potent transcriptional activator. Nature 335:563-564.[CrossRef][Medline]
28. Struhl, G., K. Struhl, and P. Macdonald. 1989. The gradient morphogen bicoid is a concentration-dependent transcriptional activator. Cell 57:1259-1273.[CrossRef][Medline]
29. Tanaka, M. 1996. Modulation of promoter occupancy by cooperative DNA binding and activation-domain function is a major determinant of transcriptional regulation by activators in vivo. Proc. Natl. Acad. Sci. USA 93:4311-4315.
30. Tautz, D., R. Lehmann, H. Schnürch, R. Schuh, E. Seifert, A. Kienlin, K. Jones, and H. Jäckle. 1987. Finger protein of novel structure encoded by hunchback, a second member of the gap class of Drosophila segmentation genes. Nature 327:383-389.[CrossRef]
31. Vershon, A. K., and A. D. Johnson. 1993. A short, disordered protein region mediates interactions between the homeodomain of the yeast
2 protein and the MCM1 protein. Cell 72:105-112.[CrossRef][Medline]
32. Xue, D., Y. Tu, and M. Chalfie. 1993. Cooperative interactions between the Caenorhabditis elegans homeoproteins UNC-86 and MEC-3. Science 261:1324-1328.
33. Yu, Y., W. Li, K. Su, M. Yussa, W. Han, N. Perrimon, and L. Pick. 1997. The nuclear hormone receptor Ftz-F1 is a cofactor for the Drosophila homeodomain protein Ftz. Nature 385:552-555.[CrossRef][Medline]
34. Yuan, D., X. Ma, and J. Ma. 1999. Recognition of multiple patterns of DNA sites by Drosophila homeodomain protein Bicoid. J. Biochem. 125:809-817.
35. Yuan, D., X. Ma, and J. Ma. 1996. Sequences outside the homeodomain of Bicoid are required for protein-protein interaction. J. Biol. Chem. 271:21660-21665.
36. Zhao, C., V. Dave, D. Fu, A. York, and J. Ma. Insights into the molecular functions of the Drosophila morphogenetic protein Bicoid. In A. Gayathri (ed.), Recent developments in molecular and cellular biology, vol. 3, in press. Research Signpost, Kerala, India.
37. Zhao, C., V. Dave, F. Yang, T. Scarborough, and J. Ma. 2000. Target selectivity of Bicoid is dependent on non-consensus site recognition and protein-protein interaction. Mol. Cell. Biol. 20:8112-8123.
38. Zhao, C., A. York, F. Yang, D. J. Forsthoefel, V. Dave, D. Fu, D. Zhang, M. S. Corado, S. Small, M. A. Seeger, and J. Ma. 2002. The activity of the Drosophila morphogenetic protein Bicoid is inhibited by a domain located outside its homeodomain. Development 129:1669-1680.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| J. Bacteriol. | J. Virol. | Eukaryot. Cell |
|---|
| Microbiol. Mol. Biol. Rev. | Clin. Vaccine Immunol. | All ASM Journals |
|---|