Previous Article | Next Article ![]()
Journal of Virology, October 2008, p. 9689-9699, Vol. 82, No. 19
0022-538X/08/$08.00+0 doi:10.1128/JVI.00995-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Institute of Molecular Biology, Friedrich-Loeffler-Institut, Südufer 10, 17493 Greifswald-Insel Riems, Germany
Received 13 May 2008/ Accepted 16 July 2008
|
|
|---|
|
|
|---|
|
|
|---|
Stable isotope labeling. The original stable isotope labeling procedure (38) was adapted as follows. Dulbecco's modified Eagle (DME)/F12 medium (D-9785; Sigma-Aldrich, Taufkirchen, Germany) was supplemented with 5% dialyzed fetal calf serum and all missing amino acids (Sigma-Aldrich) except L-leucine. Medium was then divided and supplemented with conventional or deuterated L-leucine (L-leucine-5,5,5-D3 [99 atom% D]) (catalog number 486825; Sigma-Aldrich) to produce PROLeu-DME/F12 or DEULeu-DME/F12 medium, respectively. MDBK cells were passaged in parallel in both media at a 1:10 ratio every 3 days. After four passages, aliquots of the cell cultures were lysed, and proteins were separated by gel electrophoresis. The efficiency of the exchange of normal by deuterated leucine was controlled by mass spectrometry. The incorporation of L-leucine-5,5,5-D3 increases the mass of leucine-containing peptides by 3 mass units per leucine residue.
Infection. Deuterium-labeled cells were used as mock-infected controls, and cells grown on the conventional amino acid source were infected. Cell batches passaged in the two media were seeded in 75-cm2 or 150-cm2 cell culture flasks and inoculated with virus stock corresponding to a multiplicity of infection of 10 (cells grown in PROLeu-DME/F12 medium) or mock inoculated (cells grown in DEULeu-DME/F12 medium) on ice for 1 h and then incubated at 37°C for 4 h.
Extraction of phosphoproteins. Phosphoproteins of infected and mock-infected cells were purified with the PhosphoProtein purification kit (catalog number 37101; Qiagen, Hilden, Germany). In short, 1 x 107 cells were extracted with the buffer supplied with the kit, containing 0.24% 3-[(3-cholamidopropyl)-dimethylammonio]-1-propansulfonate (CHAPS), and unsolubilized material was removed by centrifugation. The protein concentration in the clarified extract was determined with the BCA protein assay (Pierce, Rockford, IL), and equal amounts of proteins were mixed and applied onto the column supplied with the kit. The eluate containing mainly phosphoproteins and the flowthrough containing mainly the nonphosphorylated proteins were concentrated by ultrafiltration in Nanosep tubes equipped with a 10-kDa-cutoff membrane. Protein concentrations were assayed with the BCA protein assay (Pierce). Total eluates, usually containing approximately 200 µg of protein and aliquots of 1 mg of the concentrated flowthrough, were precipitated (2D clean-up kit, catalog number 80-6484-51; GE Healthcare, Braunschweig, Germany), and the pellets were dissolved in 200 µl of rehydration buffer (RHB) (7) with mild sonication and either used immediately or stored at –20°C. The resulting fractions are referred to as the phosphoprotein fraction and nonphosphoprotein fraction (flowthrough). Both fractions were tested for the phosphorylation status of the included proteins by Western blotting with antibodies directed against phosphoserine, phosphothreonine, or phosphotyrosine, and the efficiency of the separation for HeLa cells as given by the manufacturer (PhosphoProtein purification kit user manual) was confirmed for the MDBK cells used.
ASPE with cibacron blue F3G-A-Sepharose and heparin-Sepharose. Two 150-cm2 flasks each of labeled and unlabeled cells corresponding to approximately 8 x 107 cells in total were harvested in 10 ml extraction buffer (10 mM Na2HPO4-KH2PO4, 150 mM NaCl, 1% CHAPS [pH 7.0]) supplemented with protease inhibitors (Complete Mini, 1 tablet/10 ml; Roche), extracted for 60 min on ice with occasional shaking, and centrifuged (15 min at 4°C at 4,000 x g) to remove unsoluble material. Extracts from unlabeled and isotope-labeled cells were mixed at a 1:1 protein ratio. HiTrap columns (1-ml column volume) (catalog numbers 17-0412-01 and 17-0406-01; GE Healthcare) were washed with 10 column volumes each of water, elution buffer (10 mM Na2HPO4-KH2PO4, 2 M NaCl, 1% CHAPS [pH 7.0]), and extraction buffer before the extract was applied with a peristaltic pump and recirculated for 30 min. The columns were washed with 10 column volumes of extraction buffer, and bound material was eluted with 5 ml elution buffer. Samples containing 1 mg of the eluates or the flowthrough were precipitated with trichloroacetic acid, and precipitates were resuspended in RHB and stored at –20°C. The resulting fractions are referred to as the cibacron fraction, the heparin fraction, and the affinity solid-phase extraction (ASPE) flowthrough in the text. Protein yields were approximately 11% in the cibacron fraction and 17% in the heparin fraction, and the rest was found in the flowthrough. Protein recoveries approximated 100%.
Two-dimensional (2D) gel electrophoresis. Phosphoproteins were analyzed on 11-cm ReadyStrips (Bio-Rad, Munich, Germany) with a nonlinear pH range of 3 to 10, nonphosphoproteins were analyzed on 24-cm strips with a nonlinear pH range of 3 to 10, and all other fractions were analyzed on 24-cm strips with linear pH ranges of 3 to 6, 4 to 7, and 6 to 9. Precipitated samples were resuspended in RHB, briefly sonicated on ice, and extracted for 2 h at 20°C with intensive shaking. Undissolved material was removed by centrifugation (20°C at 10 min at 14,000 x g), and sample proteins were allowed to diffuse into ReadyStrips and focused in an IEF cell (Bio-Rad) according to the guidelines provided by the manufacturer. Focused strips were frozen at –80°C, thawed, and sequentially equilibrated in buffers containing dithiothreitol and iodoacetamide as recommended by the manufacturer. The second dimension was run on hand-cast full-size 12% acrylamide gels in a Dodeca cell electrophoresis chamber (Bio-Rad) with two 11-cm strips or one 24-cm strip loaded per gel. After the electrophoretic run, gels were fixed, stained overnight with colloidal Coomassie brilliant blue (36), and scanned.
Peptide mass fingerprint (PMF) analysis. After evaluation of the gel scans with Delta2D software (version 3.4; Decodon, Greifswald, Germany), lists of protein spots to be picked were used to operate a Proteineer SPII Spotpicker (Bruker Daltonics, Bremen, Germany). Tryptic digestion (46) was carried out in 96-well V-bottom polypropylene microtiter plates for 3 h with 30 ng trypsin (catalog number V5111; Promega, Mannheim, Germany) per sample at 37°C. Fingerprint and up to four tandem mass (MS/MS) spectra per sample were registered on a Bruker Ultraflex I tandem time-of-flight instrument (Bruker) and processed by flexAnalysis 2.0 software (Bruker). For the quantitation of mass-tagged peptide peak pairs, it was crucial to choose the "SNAP" option as the peak detection algorithm in the flexAnalysis software, which is robust with respect to overlapping isotope patterns. A batch database search (MASCOT Server 2.0.0 software; Matrix Science Ltd., London, United Kingdom) (42) was launched by Biotools 2.2 (Bruker) using the bovine International Protein Index (IPI) database (www.ebi.ac.uk) (24) or an in-house database covering the PrV proteome as compiled from the Swiss-Prot database (www.expasy.org) (13). Carbamidomethylation was set as a fixed modification for cysteine residues, the significance level was set to 95%, and mass tolerance was set to 50 ppm for the fingerprint spectra. Proteins were considered as identified and selected for quantitative evaluation if significant molecular weight search (MOWSE) scores (39) were obtained with the fingerprint spectrum or in a combined search of the fingerprint and MS/MS spectra.
Quantitation and data processing. Quantitation was carried out by in-house software (AMaDEuS) based on Visual Basic for applications (Microsoft, Redmond, WA) macros. For identified samples, peaks representing the protein with the highest MOWSE score were selected, and the masses of the expected isotope-labeled peaks within error margins of 15 ppm were calculated on the basis of the number of leucine residues present in the peptide and the mass shift of +3 Da per leucine residue. If the spectrum contained the mass of the expected deuterated peak, the peak pair was selected for quantitative evaluation, and the intensity ratio of the peaks was calculated. If more than four peak pairs were obtained from one spectrum, outliers were eliminated by the symmetrical truncation of intensity ratios. Truncation was limited to a maximum of three rounds or the number of rounds that would leave a minimum of three values. Truncated mean values and standard deviations of the intensity ratios of qualified peak pairs representing the identified protein were calculated. Relative standard deviations rarely exceeded 15%. The database query was then repeated with a fixed modification of +3 Da for leucine residues to ensure the identification of proteins that were strongly downregulated after infection. Quantitative evaluation was again carried out as described above, with the only difference being that a –3-Da shift was used to calculate the expected mass of the unlabeled peptide peak in the AMaDEuS software. Results of both queries and calculations were combined and compared. Data sets resulting from the same sample (protein spot) yielding higher MOWSE scores were preferred to those with lower scores. The final output for every sample was the truncated mean of the ratios (unlabeled over labeled, that is, infected over mock [IOM]) of qualified peak pairs, which will be abbreviated as IOM ratios. IOM ratios greater than 1.0 indicate higher relative levels, and IOM ratios below 1.0 indicate lower relative levels of the respective proteins after infection with PrV. The average number of qualified peptides in identified proteins was 3.3, and the maximum number was 18.
|
|
|---|
![]() View larger version (18K): [in a new window] |
FIG. 1. Flowchart of sample preparation and analysis. Note that all separation and analytical procedures (affinity extractions, 2D gel electrophoresis, tryptic digestion, and mass spectrometric analysis) were carried out with a 1:1 protein mixture of extracts of conventional and heavy-isotope-labeled cells. After extraction with three different affinity matrices and separation in different pH ranges, a total of 11 2D electrophoretic gels were analyzed for each experiment. FT, flowthrough.
|
In the course of this study, over 4,000 proteins were identified, and 2,374 proteins were quantified from a minimum number of three qualified peptide pairs. Of the 1,490 spots quantified in infection experiments, IOM ratios of 109 samples (7.3%) representing 55 genes (Table 1) were beyond the empirical 1% and 99% quantiles of 0.63 and 1.63. Thus, approximately 30 of the 109 samples were expected to result from the statistical variance of the experiment itself, and approximately 79 samples were expected to result from virus-induced changes.
|
View this table: [in a new window] |
TABLE 1. Proteins up- or downregulated after infection with PrV-Kad
|
![]() View larger version (58K): [in a new window] |
FIG. 2. Scans of 2D electrophoretic gels representing the two affinity-purified fractions (BLUE, cibacron fraction; HEP, heparin fraction) and the flowthrough (FT) of the ASPE separation procedure. Differing protein patterns reflect an efficient separation into three well-defined protein fractions. Boxes indicate gel regions that were analyzed in more detail (Fig. 6 and 7) for modifications of eIF-4B and SNX-9 (box A), lamin B2 (box B), 60S acidic ribosomal protein P0 (box C), Hsp27 (box D), and hnRNP K (box in panel FT).
|
![]() View larger version (10K): [in a new window] |
FIG. 3. Graphic representation of separation results. Overlaps were calculated on the basis of the number of entries of identified proteins found in the IPI database. Shown are overlaps between the different fractions of the ASPE (BLUE, cibacron fraction; HEP, heparin fraction; FT, flowthrough) (A) and overlaps between the cibacron, heparin, and phosphoprotein (PP) fractions (B) demonstrating a high selectivity of the ASPE procedure (A) and a good complementarity of the ASPE with phosphoprotein extraction (B).
|
![]() View larger version (24K): [in a new window] |
FIG. 4. Distribution of isotope ratios of proteins identified in a control experiment with labeled and unlabeled mock-infected cells. Isotope ratios were calculated as described in Materials and Methods. The frequency of isotope ratios was registered in steps of 0.01. As expected, the distribution centers around 1.0. Quantiles are given in percentages, and the corresponding isotope ratios are shown as plain numbers.
|
![]() View larger version (20K): [in a new window] |
FIG. 5. Distribution of relative cellular protein levels 4 h after infection with PrV-Ka. Error bars indicate the standard deviations of the isotope ratios of the peptides used for the calculation of the IOM ratio. PP, phosphoprotein fraction; PP-FT, flowthrough of the PhosphoProtein purification kit (nonphosphoproteins); BLUE, cibacron fraction; HEP, heparin fraction; FT, ASPE flowthrough. Cutoff values of 0.63 and 1.63 are given as horizontal lines.
|
![]() View larger version (36K): [in a new window] |
FIG. 6. (A) Detailed view from a 2D electrophoretic gel of a heparin fraction analyzed in the pH range 4 to 7 (box D in Fig. 2). Five protein spots were identified as being Hsp27, and the relative protein levels, given in numbers, were calculated as described in the text. Details of mass spectra originating from spot A (most acidic isoform) and spot E (least acidic isoform) show two peptide peaks containing one (1,413 Da) or two (1,163 Da) leucine residues. IOM ratios were calculated from intensity ratios of peak pairs with a mass distance corresponding to the number of leucine residues present, in this case 3 Da (1,413 Da) and 6 Da (1,163 Da). The respective deuterium-labeled peaks at 1,416 Da and 1,169 Da representing peptides from mock-infected cells are marked with asterisks. Both labeled and unlabeled peptides occur in their natural isotope patterns, giving rise to multiple peaks with descending intensities in 1-Da distances. Note that reliable relative quantitation was restricted to ratios between 1:10 and 10:1 so that IOM ratios of the most acidic spots of Hsp27 and lamin C (Fig. 7A) could be assessed only as "greater than 10".
|
![]() View larger version (100K): [in a new window] |
FIG. 7. Proteins that were differentially regulated in different isoforms. Relative abundances of the respective isoforms after infection are given in plain numbers.
|
Lamins. Lamin A/C was found mainly in the phosphoprotein fraction but was also present in the nonphosphoprotein and heparin fractions of the ASPE in two strings of poorly resolved protein spots representing the unprocessed lamin A and the approximately 20-kDa-smaller proteolytic cleavage product, lamin C (Fig. 7A). For both proteins, a shift to more acidic isoforms was observed after infection. As found for the most acidic modification of Hsp27, the protein spot representing the most acidic variant of lamin C contained hardly any material from mock-infected cells. Lamin B2 was found in one major and one minor more acidic spot (Fig. 7B). Considering the relative protein levels and the distribution of the total mass between the two spots, gains and losses are not nearly balanced so that under the condition that no other lamin B2 isoforms have been missed, a total loss of protein occurred after infection. Lamin B1 was also identified in three protein spots with quantitative variations within the isoforms being less pronounced than those for lamin B2, although for lamin B1 increasing abundances were also found with increasing acidity (IOM ratios from the most acidic spot to the least acidic spot were 1.77, 1.08, and 0.80).
Proteins involved in translation. Although numerous ribosomal proteins and translation factors were identified and quantified, only one ribosomal protein (60S acidic ribosomal protein P0) (Fig. 7C) and two translation initiation factors showed a significant (eukaryotic initiation factor 4, subunit B [eIF-4B]) (Fig. 7D) or a moderate response (eIF-3F) (Table 1) to PrV infection. eIF-4B was present in seven charge variants most likely representing different phosphorylation levels.
hnRNPs. hnRNPs were highly enriched in the cibacron fraction with the exception of hnRNP K, which reliably separated into the flowthrough of the ASPE (Fig. 2). A number of hnRNP representatives like hnRNP A3 (IOM of 2.80 ± 0.23), hnRNP D (IOM of 0.42 ± 0.09), and hnRNP A2/B1 (IOM of 2.69 ± 1.50) showed significant variation after infection, while others (hnRNP A1 and hnRNP F) remained constant or were only slightly changed (IOM of 0.62 ± 0.1 for hnRNP H and IOM of 0.60 ± 0.16 for hnRNP A/B). Of the six spots that were identified as being hnRNP K (Fig. 7E), two migrated with a slightly smaller apparent molecular weight and decreased in abundance after infection, whereas relative levels of the slightly larger proteins increased, with gains being more pronounced for the more acidic isoforms.
Proteins related to intracellular transport and the cytoskeleton. Sorting nexin 9 (SNX-9) (Fig. 7D), enriched in the heparin fraction, appeared in four distinct spots, which followed the observed shift to more acidic isoforms. However, in this case, the mass gain was not highest for the most acidic form. A functionally related protein, ezrin-radixin-moesin-binding phosphoprotein 50 (EBP50) was found with decreased relative levels (IOM of 0.48 ± 0.06) in the heparin fraction. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which, apart from its function in the glycolytic pathway, is involved in intracellular membrane trafficking (55, 57), was present in the cibacron fraction and found to be markedly upregulated (IOM of 2.32 ± 0.30). Relative levels of the other glycolytic enzymes that were identified and quantified (fructose-bisphosphate aldolase, triosephosphate isomerase, phosphoglycerate kinase, enolase, and pyruvate kinase) did not differ significantly between infected and noninfected cells.
|
|
|---|
Overall, the cellular proteome appeared to be very stable. The relative levels of the vast majority of proteins were unaffected by infection with PrV, and the calculated IOM ratios were close to 1.0 despite the known detrimental effects of PrV infection on the stability of mRNA and cell protein synthesis due to viral host cell shutoff functions (29, 49). Thus, the delayed character of PrV-induced host cell shutoff indicated by microarray studies in porcine (12), rat (44), and human (6) cells was confirmed on the protein level for the bovine cells used here. From time course studies with unfractioned material and analysis with 2D electrophoresis in the pH range of 3 to 10 (data not shown), there was no indication for a global decrease of protein levels up to 8 h after infection with PrV, indicating that the degradation of cellular mRNA by viral pUL41 and other mechanisms of host cell shutoff seem to exclusively target the RNA metabolism but do not interfere with the physiological steady-state levels of most proteins, at least not at early times after infection.
Nevertheless, the IOM ratios of a considerable number of protein spots did change significantly in either direction. Moreover, in several cases, e.g., Hsp27, lamin A/C, lamin B2, hnRNP K, SNX9, eIF-4B, or 60S acidic ribosomal protein P0, inverse modulations of different charge variants of the same protein were observed.
Stress response proteins. The most acidic of the five isoforms of Hsp27 was detected almost exclusively in infected cells, indicating that PrV infection caused a sharp rise in these isoforms. Most probably, the charge variants are the result of differential phosphorylation, which was previously described for human (37, 48) and bovine (28) Hsp27. A similar shift to higher phosphorylated forms of Hsp27 has been observed after treatment of bovine cells with cadmium (28), indicating that this reaction is not specific for infections with herpesviruses but rather reflects phosphorylation of Hsp27 during cellular stress (27). However, other stress-related proteins that were identified and quantified showed no significant changes in relative protein levels (data not shown). Hsp27 is a multifunctional protein (1) with strong antiapoptotic properties. It is unclear so far if the changes in levels of Hsp27 are related to the pUS3-mediated suppression of apoptosis (14), which is observed after infection with PrV.
Lamins. The observed shift to more acidic isoforms of lamin A/C parallels the reported changes of 2D gel electrophoretic patterns of lamin A/C caused by the viral kinase pUS3 after infection with HSV-1 (35). In an ongoing study with a US3-deleted PrV mutant, the impact of pUS3(PrV) on the phosphorylation pattern of lamin A/C will be addressed. Infection with HSV-1 is also accompanied by an intracellular rearrangement (5) and loss (50) of B-type lamins, which were identified and quantified from the heparin fraction. A shift to more acidic isoforms was observed for both, which corresponds to the previously reported phosphorylation of B-type lamins following HSV-1 infection (40). The hyperphosphorylation of the lamins precedes the disintegration of the nuclear lamina in many physiological processes, which is presumably a prerequisite for the transfer of PrV capsids from the nucleus into the cytoplasm for secondary envelopment (34).
Proteins related to translation. After infection with HSV-1, the synthesis of several ribosomal proteins and their assembly into ribosomes continue in spite of a general inhibition of cellular protein synthesis and a concomitant loss of the mRNA encoding the respective ribosomal proteins (15, 52). With the exception of the redistribution seen within the isoforms of the 60S acidic ribosomal protein P0, we also found stable levels of ribosomal proteins after PrV infection. Likewise, significant changes in relative levels of translation factors were not observed, with the exception of eIF-4B and the F subunit of eIF-3.
Proteins related to intracellular transport and the cytoskeleton. So far, none of the three proteins related to intracellular transport and the cytoskeleton (SNX9, EBP50, and GAPDH), which we found to be significantly altered in abundance, have been correlated with any step during herpesvirus replication. However, they exhibit interesting features, which may have functions related to herpesvirus infections. Members of the SNX family of proteins contain a Phox domain that mediates binding to membrane-anchored phospholipids. They are involved in intracellular membrane trafficking (51), which is important for the replication of enveloped viruses. SNX9 shares a relative promiscuity with other SNXs with respect to the different phosphatidylinositols that are bound by the Phox domain, the presence of an SH3 domain, and the presence of a BAR dimerization domain. SNX9 is activated by phosphorylation (30), which correlates with the observed shift to more negatively charged variants, which we observed after infection (Fig. 7D). After phosphorylation, SNX9 is translocated to the plasma membrane, where it can interact with transmembrane proteins, but also recruits dynamin to the plasma membrane (30). It is also required for effective clathrin-mediated endocytosis (53). By virtue of its BAR domain, it may be involved in the stabilization but also in the sensoring of strongly curved membranes (30), which occur during budding.
EBP50 is a peripheral membrane protein that resides at the plasma membrane and supports functions of the ezrin-radixin-moesin proteins in connecting the actin cytoskeleton to the plasma membrane. Cytomorphological alterations of PrV infection have been attributed to actin stress fiber breakdown mediated by the pUS3 kinase (58), which may occur via the phosphorylation of EBP50.
Apart from its function in glycolysis, GAPDH plays a major role in the early steps of the secretory pathway during the transport of cargo from the endoplasmic reticulum to the Golgi apparatus (55). As was exemplarily demonstrated with the G glycoprotein of vesicular stomatitis virus, the transport function of GAPDH was independent of its enzymatic activity but depended on the phosphorylation of GAPDH mediated by the cellular Src protein (56, 57). In the context of PrV infection, GAPDH may play a role in the transport of viral glycoproteins to late secretory compartments like the trans-Golgi network, where the PrV particles acquire their final envelope in a budding process (33).
hnRNPs. Four types of hnRNP were significantly (hnRNPs A3, A2/B1, D, and K) and two were moderately (hnRNPs H and A/B) modulated after infection with PrV. Since infection with alphaherpesviruses is accompanied by drastic changes in RNA metabolism, which is, at least in part, caused by herpesviral shutoff mechanisms, the modulation of hnRNPs, which are involved in RNA processing and turnover, was not unexpected. The host shutoff protein pUL41 of HSV-1 acts as an endoribonuclease, preferentially degrading mRNA containing AU-rich elements (11). In this context, the loss of hnRNP D, which is also known as AU-rich element binding protein 1 (AUF-1), might be important for the function of pUL41.
Another viral protein that is involved in HSV-1-mediated host cell shutoff, ICP27, interacts with hnRNPs K and CK2 (59), which phosphorylates hnRNP K in an ICP27-dependent manner (25). The impact of this phosphorylation, which is probably reflected by the observed shift to more acidic variants of hnRNP K, in the progress of the infection is still unclear but might be of importance across the herpesvirus family, since a similar interaction between hnRNP K and CK2 and the ORF57 protein has been found in Kaposi's sarcoma-associated herpesvirus-infected cells (32).
The proteomic screen described here resulted in the identification of a number of cellular proteins that are potentially of significance for herpesvirus infection. The fact that several of them have previously been linked to the replication cycle of alphaherpesviruses (lamin A/C, the B-type lamins, and hnRNP K) serves to validate our experimental system. The newly established protocol can be easily adapted to study virus-host cell interactions of any virus that can be propagated in cultured cells of a species with a sequenced and annotated genome. It can be further refined by the application of additional affinity matrices or the use of narrow-range isoelectric focusing strips to improve yields of less abundant proteins. The demonstration that most of the observed changes affected the distribution of protein in different posttranslationally modified isoforms rather than absolute expression levels underlines the importance of analysis at the protein level for a comprehensive understanding of the molecular effects of virus infections on the cellular metabolism.
We thank Barbara Bettin for expert technical assistance.
Published ahead of print on 23 July 2008. ![]()
|
|
|---|
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»