Previous Article | Next Article ![]()
Molecular and Cellular Biology, March 2007, p. 1859-1867, Vol. 27, No. 5
0270-7306/07/$08.00+0 doi:10.1128/MCB.01395-06
Copyright © 2007, American Society for Microbiology. All Rights Reserved.
,
Molecular Oncology Research Institute, Tufts-New England Medical Center, Boston, Massachusetts 02111,1 Department of Experimental and Diagnostic Medicine and Interdepartmental Center for Cancer Research, University of Ferrara, Ferrara 44100, Italy,2 Comprehensive Cancer Center, Ohio State University, Columbus, Ohio 432103
Received 28 July 2006/ Returned for modification 28 August 2006/ Accepted 10 December 2006
|
|
|---|
|
|
|---|
Moreover, global microRNA expression changes have been described to occur in human cancers and in some cases shown to correlate with the clinico-pathological features of the tumor (2, 13, 29). However, no mechanism has been proposed to date for these profile alterations.
Despite this wealth of data, relatively little is known about microRNA regulation and response to microenvironmental factors. One mechanism involves the activation of specific signal transduction pathways that in turn promote the transcription of certain microRNAs. For example, it was reported that the miR-1 genes are targets of serum response factor, a converging downstream effector for a variety of oncoproteins and growth factors (30). Another transcription factor, the c-myc oncogene product, was also found to activate the expression of a microRNA cluster (25).
Hypoxia is an essential feature of the neoplastic microenvironment. Tumors with widespread low oxygenation tend to exhibit increased invasion and resistance to conventional therapy (9). The molecular mechanisms responsible for the hypoxic survival of neoplastic cells are not fully characterized, and a better understanding of this process may lead to novel strategies for pharmacological intervention.
Our data indicate that hypoxia leaves a specific mark on microRNA profiles in a variety of cell types, with a critical contribution of the hypoxia-inducible factor (HIF). Moreover, at least a subgroup of these hypoxia-regulated microRNAs (HRMs) seem to play a role in cell survival in a low-oxygen environment.
Finally, by comparing hypoxia-associated microRNA spectra with published data from a large number of tumors (28), we propose that cancer-associated microRNA profiles exhibit a hypoxic signature.
|
|
|---|
Plasmids.
The mutant versions of HIF-1
and HIF-2
with double proline-to-alanine substitutions have been described previously (14, 15). The three-HRE-thymidine kinase (tk)-luciferase reporter is a HIF-responsive construct containing a tandem of hypoxia-responsive elements from the erythropoietin promoter. Both plasmids were provided by William Kaelin (Dana-Farber Cancer Institute, Boston, MA).
MicroRNA microarray analysis. RNA was extracted by using TRIzol (Invitrogen, Carlsbad, CA), according to the manufacturer's protocol. RNA was labeled and hybridized on microRNA microarray chips as previously described (20). Briefly, 5 µg of RNA from each sample was biotin labeled during reverse transcription using random hexamers. Hybridization was carried out on the second version of our microRNA chip (ArrayExpress accession number A-MEXP-258), which contains 381 probes for mature and precursor human microRNAs and 393 probes for mouse microRNAs, together with controls. There are four spot replicates for each probe on the chip. Hybridization signals were detected by biotin binding of a Streptavidin-Alexa 647 conjugate using a GenePix 4000B scanner (Axon Instruments). Images were quantified using the GenePix Pro 6.0 apparatus (Axon Instruments). We used the microRNA nomenclature according to the microRNA Registry (miRBase http://microrna.sanger.ac.uk/sequences/) at the Sanger Institute and the Genome Browser (http://genome.ucsc.edu).
Analysis of microarray data. Raw data were normalized and analyzed by GeneSpring GX software version 7.3 (Agilent Technologies). The GeneSpring software generated a unique value for each microRNA, determining the average of the results from the four spots. Samples were normalized using the on-chip median normalization feature of GeneSpring software; the result for each cell line time course experiment was normalized to that at time zero. MicroRNAs showing an altered expression across the time course were identified using the filter of the n-fold-change tool. In detail, a GeneSpring GX filter with the n-fold-change tool was used to point out the microRNAs that, for at least two cell lines, were 1.5-fold upregulated at 24 h compared to expression levels at time zero. The 1.5-fold-change threshold was chosen on the basis of its use in previously published articles employing these particular types of microarrays. Analysis of variance (ANOVA) with the Benjamini and Hochberg correction for false-positive reduction was used to compare the average values obtained with microRNA probes in all cell lines at 24 and 48 h with values at time zero. Hierarchical cluster analysis was performed using average linkage and Pearson correlation as a measure of similarity.
Transfections. Plasmid transfections were performed using Lipofectamine 2000 (Invitrogen) as per the manufacturer's protocol. Fresh medium was added after 5 h of transfection, and RNA and protein were harvested 24/48 h posttransfection. Transfection of mature microRNAs or inhibitors (Ambion, Inc.) was performed using siPORT NeoFX transfection agent (Ambion, Inc.) according to the manufacturer's protocol. Transfection complexes were added to cells at a final oligonucleotide concentration of 20 nM, and medium was replaced 24 h later.
Luciferase assays using reporters based on HRM promoter fragments. We started by amplifying HRM promoter fragments predicted to encompass HIF sites, as shown in Table S1b in the supplemental material. Thus, the genomic region 4.9 kb immediately upstream of miR-24-1 was amplified using primers 5'ATACTCGAGCTGCTAGGCCATGCGTGTCC3' (forward) and 5'ATTAAGCTTCAAGAGAGAGTTCACCGCGC3' (reverse) (underlining indicates restriction sites inserted). For miR-181c, a region of approximately 2.0 kb immediately upstream was PCR amplified using 5'ATAGGTACCCACTCCACAGCCTGAATG3' (forward) and 5'TATAAGCTTGGTGGGGTAGGTGGCAGGGAAC3' (reverse). For miR-26b, a region situated approximately 3 kb upstream of the 5' end (and encompassing a cluster of four predicted HIF sites) was PCR amplified using the following primers: 5'ATAGCTAGCGAGACAGATGTCCCGCTCCCAG3' (forward) and 5'ATCGCTAGCACGCTCTTGAATGGGACGG3' (reverse). Due to the size and CG contents, PCR amplification was done using the Herculase II fusion DNA polymerase (Stratagene). The resulting fragments were cloned in the pGL3 basic vector (Promega) (miR-24-1 and -181c) or pGL3-tk-luciferase vector (from David Fisher, DFCI) (miR-26b). MCF7 or HT29 cells were cotransfected using Lipofectamine 2000 (Invitrogen) with the reporter plasmids and either of the HIF mutants or the empty vector. Luciferase assays were performed 24 h posttransfection by following the manufacturer's protocol (luciferase assay system; Promega). Luciferase activity (measured as relative light units) was scored using a Femtomaster FB12 luminometer (Zylux Corporation) and normalized to that of an equal protein concentration in all samples.
Apoptosis assays. Cells were plated in triplicate in a 96-well plate, transfected with microRNA precursors or inhibitors (Ambion, Inc.) as described above, and scored for caspase activation 72 h posttransfection using an Apo-ONE homogeneous caspase-3/7 assay kit (Promega) according to manufacturer's instructions. The protein concentration was determined using Bradford's method, and caspase activities for all samples were normalized to that of an equal protein amount. The data are presented as means plus standard deviations from three independent experiments performed in triplicate.
Northern hybridization.
Total RNA isolation was performed using the TRIzol reagent (Invitrogen) according to the manufacturer's instructions. RNA samples (20 µg each) were run on 15% acrylamide-denaturing (urea) Criterion precast gels (Bio-Rad Laboratories, Hercules, CA) and then transferred onto a Hybond-N+ membrane (Amersham Pharmacia Biotech). The hybridization with [
-32P]ATP was performed at 42°C in 7% sodium dodecyl sulfate (SDS)-0.2 M Na2PO4 (pH 7.0) overnight. Membranes were washed at 42°C, twice in 2x SSPE (1x SSPE is 0.18 M NaCl, 10 mM NaH2PO4, and 1 mM EDTA [pH 7.7])-0.1% SDS and twice with 0.5x SSPE-0.1% SDS. Blots were stripped by boiling them in 0.1% aqueous SDS-0.1x SSC for 10 min and were reprobed several times. As a loading control, we used 5S rRNA stained with ethidium bromide or probing for the U6 snRNA.
TaqMan RT-PCR. The method was optimized for microRNA, and reagents, primers, and probes were obtained from Applied Biosystems. Human 18S rRNA was used to normalize all RNA samples. Reverse transcriptase (RT) reactions and real-time PCR (PCR) were performed according to manufacturer protocols. RNA concentrations were determined with a NanoDrop apparatus (NanoDrop Technologies, Inc.), and 1 nanogram per sample was used for the assays. All RT reaction mixtures, including no-template controls and RT-minus controls, were run in duplicate in an Applied Biosystems 9700 thermocycler. Gene expression levels were quantified using the ABI Prism 7900HT sequence detection system (Applied Biosystems). Analysis was performed in duplicate, including with no-template controls. Relative expression was calculated using the comparative cycle threshold method.
ChIP.
HT29 cells (60% confluent) were exposed to hypoxia (0.2% O2) or normoxia (21% O2) for 24 h. Cross-linking was performed using 1% (vol/vol) formaldehyde for 10 min, and the reaction was stopped with 125 mM glycine for 5 min. Cells were washed twice with ice-cold 1x phosphate-buffered saline and collected by scraping them in 1 ml of 1x phosphate-buffered saline supplemented with 1 mM phenylmethylsulfonyl fluoride (PMSF), followed by centrifugation and lysis in 400 µl of buffer (1% SDS, 10 mM EDTA, 50 mM Tris-Cl, pH 8.0, protease inhibitor cocktail, and 1 mM PMSF). The resulting lysate was incubated on ice for 10 min, followed by sonication for DNA shearing to fragments between 200 bp and 1 kb. The supernatant was recovered and diluted in chromatin immunoprecipitation (ChIP) dilution buffer (1% Triton X-100, 2 mM EDTA, 150 mM NaCl, 20 mM Tris-Cl, pH 8.0, protease inhibitor cocktail, 1 mM PMSF), followed by immunoprecipitation overnight at 4°C with a polyclonal anti-HIF-1
antiserum (ab2185; Abcam) or whole rabbit antiserum (immunoglobulin G [IgG] control). Immunocomplexes were recovered by the addition of a 50% slurry of salmon sperm DNA-protein A-agarose (Upstate) to the samples and sequentially washed for 4 min each in buffer I (20 mM Tris, pH 8.0, 200 mM NaCl, 0.5% Triton X-100, 0.05% deoxycholate, 0.5% NP-40, 1 mM PMSF), buffer II (20 mM Tris, pH 8.0, 500 mM NaCl, 0.5% Triton X-100, 0.05% deoxycholate, 0.5% NP-40, 1 mM PMSF), buffer III (10 mM Tris, pH 8.0, 250 mM LiCl, 1% Triton X-100, 0.1% deoxycholate, 0.5% NP-40, 0.5 mM EDTA, 1 mM PMSF), and buffer IV (10 mM Tris, pH 8.0, 5 mM EDTA). The immunoprecipitated DNA was retrieved from the beads by incubation in elution buffer (10 mM Tris, pH 8.0, 5 mM EDTA, 1% SDS) at 65°C for 1 h. Cross-linking was reversed using 200 mM NaCl at 65°C overnight followed by proteinase K digestion at 47°C for 2 h. DNA was then purified using a PCR purification kit (QIAGEN), and PCR was performed using primers spaced approximately 15 to 250 bp apart and encompassing predicted HIF binding sites in miR-210 and -26a-2 promoters (see below).
5'GGAGCCTTGACGGTTTGACC 3' (forward) and 5'CGAGGACCAGGGTGACAGTG3' (reverse) were used to PCR amplify the miR-210 promoter fragment (210-A) containing predicted HIF binding sites located at positions 1720 and 1822. For the 210-B fragment containing a predicted HIF consensus at position 1166, the following pair was used: 5'GGTCGTGAAGGAGCCTCTAAG3' (forward) and 5'GGACCATCCCACTCACTACAG3' (reverse). Finally, the miR-26a-2 promoter fragment (26-A), encompassing a predicted HIF binding site at position 2860, was amplified using 5'CCAAGGACTATGCACATCACC3' (forward) and 5'GGAAAGGCAGTGAGATGGCAC3' (reverse). The following primers were used to amplify a region in the promoter of miR-130b (negative control, hypoxia nonresponsive): 5'GCGAAACCCCAGCTCTACTA3' (forward) and 5'ACACTCTCACTCTGTCGCCC3' (reverse). For HIF site localization and sequences, see Table S1b in the supplemental material. The PCR products were resolved on a 1.5% agarose gel and visualized by ethidium bromide staining.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Summary of hypoxia-regulated microRNAsa
|
![]() View larger version (65K): [in a new window] |
FIG. 1. Coordinated hypoxic changes of microRNA expression in colon and breast cancer cell lines. Cluster analysis of four cell lines according to the expression of microRNAs upregulated by hypoxia in at least two cell lines. Expression data were normalized to expression at time zero.
|
![]() View larger version (25K): [in a new window] |
FIG. 2. Confirmation of HRM induction by hypoxia by Northern analysis or quantitative RT-PCR. (a) Northern blotting confirmation of miR-210 induction under hypoxia (the mature form is indicated). Lanes NOR, 6HY, and 30HY show results under normoxia and at 6 h and 30 h under hypoxia, respectively. An ethidium bromide-stained gel picture is shown as a loading control. (b) Quantitative RT-PCR confirmation of HRM induction by 24-hour hypoxia (H) compared with HRM induction for normoxic controls (N). Bars indicate means from two independent experiments. I bars indicate standard deviations.
|
In order to investigate whether HRMs exhibit a significant enrichment in predicted HIF binding sites, we performed an in silico search upstream of the genomic sequences that encode all the known microRNAs (1,039 sequences, experimentally demonstrated and predicted; http://microrna.sanger.ac.uk/sequences/ftp/genomes/hsa.gff).
Since only a few microRNA promoters have been identified experimentally, a 6-kb region (5 kb upstream and 1 kb downstream region of the 5' end of the annotated microRNA) was designated as a putative promoter sequence. Indeed, most of the experimentally confirmed transcription factor binding sites are located within the kb 5 region of the transcription start site. Predicted HIF binding sites were analyzed using the MATCH program and the V$HIF1_Q3 (GNNKACGTGCGGNN; boldfacing indicates the core HIF consensus), V$HIF1_Q5 (NGTACGTGCNGB), and V$HIF1_Q6 (NRCGTGNGN) position weight matrices from the TRANSFAC database (version 9.1) (23). The position weight matrices describe the position preferences of different nucleotides in the HIF binding site. We scanned regions around the transcription start sites of all the microRNAs from kilobase positions 5 to +1 using the "minFP_good91.prf" profile (the profile of cutoff values with a minimum number of false-positive predictions) of MATCH, similarly to searches described previously (16, 26, 27).
We then tested whether these consensus sequences are significantly more abundant in the promoters of the 23 HRMs (target set) than in the rest of microRNA-ome. Thus, we generated 50,000 groups, each consisting of 23 promoters randomly selected from the 1,039 microRNAs and calculated the number of promoters that contain at least one predicted HIF binding site in both target and random sets of promoters. The search was performed separately for the HIF_Q3 and HIF_Q5 consensus sequences. The selection was performed using the function "sample," which is part of the random module in the Python 2.4 programming language (http://www.python.org/), followed by data analysis using the R programming language (http://www.r-project.org/). The results of the analysis are summarized in Fig. 3 and Table S1a in the supplemental material, indicating a highly significant enrichment of the HIF binding consensuses in the HRM group (P = 0.00294 for HIF_Q3; P = 0.011 for the HIF_Q5 consensus). The search based on the HIF_Q6 matrix did not yield a significant P value (not shown), which is not surprising, given the high probability of such a short and degenerate sequence arising by chance very often in the genome.
![]() View larger version (17K): [in a new window] |
FIG. 3. Distribution of random 23-microRNA groups (samples) based on HIF1_Q3 (a) or HIF1_Q5 (b) binding sites. The arrow indicates where the experimental data (HRMs) fall within the random sample population, with the corresponding P value.
|
subunits versus a vector-only control (pcDNA3.1) in HT29 and MCF7 cell lines under normoxia. HIF stabilization was achieved by substituting the two prolines (at positions 564 and 402 in the case of HIF-1) in the alpha subunits that are subject to oxygen-dependent hydroxylation and proteasomal degradation via VHL-dependent ubiquitylation (15, 16, 22). The activity of exogenous HIFs was confirmed by cotransfection with an HRE-tk-luciferase reporter (containing three hypoxia response elements) followed by standard luciferase assay (Fig. S2a in the supplemental material). Expression of miR-103, -210, and -213 was measured following HIF transfection using a modified real-time RT-PCR protocol which revealed a robust and reproducible increase in the expression of all three transcripts (Fig. 4).
![]() View larger version (50K): [in a new window] |
FIG. 4. Effect of HIF on specific microRNA expression. The impact of exogenous constitutively active HIF-1 (HIF1P/A) and HIF-2 (HIF2P/A) subunits on miR-103, -210, and -213 expression was determined by quantitative RT-PCR in HT29 (a) and MCF7 (b) cells. The control was the pcDNA3.1 empty vector (PC). Bars indicate means from two independent experiments. I bars indicate standard deviations.
|
![]() View larger version (24K): [in a new window] |
FIG. 5. Direct effect of HIF in the up-regulation of select HRMs. Relative luciferase activities of HRM promoter reporter constructs in HT29 cells. (a) miR-24-1 and -181c promoters in a pGL3 context; (b) miR-26b fragment in a pGL3-tk context. The constructs were cotransfected with constitutively active HIF-1 (HIF1P/A), HIF-2 (HIF2P/A), or the empty vector pcDNA3.1 (PC) and incubated under normoxia (NOR). The effect of hypoxia (HYP) on the reporter is also shown. Bars indicate means from three independent experiments. I bars indicate standard deviations.
|
antibody (but not the control IgG antibody) immunoprecipitated the miR-210 and miR-26a-2 promoter fragments in hypoxic HT29 cells but very little in the normoxic controls. A similar assay performed for a region upstream of miR-130b (which is not an HRM), did not reveal measurable recruitment of HIF upon hypoxia exposure.
![]() View larger version (40K): [in a new window] |
FIG. 6. Direct recruitment of HIF on the miR-210 and miR-26 promoters under hypoxia. Chromatin was immunoprecipitated from HT29 cells using a HIF-1 antibody or an IgG control, and the enriched genomic fragment was amplified using primers spanning the candidate HREs located at positions 1720 and 1822 (210-A); 1166 (210-B) of the miR-210 promoter; or 2860 for the miR-26a-2 promoter (26-A). A fragment of the miR-130b promoter was used as a negative control.
|
![]() View larger version (41K): [in a new window] |
FIG. 7. Antiapoptotic effect of select microRNAs under hypoxia. Blockade of miR-26, -107, and -210 with antisense inhibitors leads to an increased apoptotic response in three independent experiments (each performed in triplicate). In contrast, an excess of sense microRNAs decreases the apoptotic response. The dotted line represents the apoptotic caspase-3/7 baseline activity in response to negative-control microRNA under hypoxia (P, precursor [sense]; A, antisense). (b) Northern blot confirmation of efficient transduction or blockade of miR-210 in MCF7 cells. MCF7 cells were transfected with the precursor or antisense miR-210 or the scramble control (SCR). U6 snRNA is shown as the loading control.
|
Correlations between cancer and hypoxia-specific microRNA profiles. Recent investigations have dissected a large number of human neoplasms for microRNA expression (13, 21, 28) and identified specific alterations from normal cells; however, the mechanism and biological impact of these changes remain elusive. In order to address a potential correlation between the pattern of microRNAs altered in solid cancers and under hypoxia, we took advantage of the largest genome-wide microarray profiling study published to date, including 540 tumor samples from six types of solid cancer (breast, lung, colon, stomach, and pancreatic endocrine tumors and prostate carcinomas) and corresponding normal samples (28). From the 228 microRNAs, 137 exhibiting expression values above threshold in at least 90% of samples were retained, which generated a "common signature" of abnormally expressed microRNAs (presented in Fig. 2B and supplementary Tables 10 and 11 in reference 28). Of note, the microarray search for HRMs was performed using the same profiling technology.
The vast majority of HRMs identified by our study (Table 2 ) are also overexpressed in at least some types of tumors, suggesting that hypoxia may represent a key contributing "trigger" for microRNA alterations in cancer.
|
View this table: [in a new window] |
TABLE 2. Solid correlations of microRNAs with cancersa
|
|
|
|---|
We are fully aware that the microarray-based strategy leaves open the possibility that other microRNAs may respond to hypoxia and were simply not detected by the screen. Indeed, with consideration of the well-recognized technical limitations of microarrays, the microRNAs' world is still expanding and it is conceivable that real HRMs were not present on the chip at that stage.
On the other hand, it is possible that some of the microRNAs identified as upregulated by microarray analysis, but not validated independently by quantitative PCR, are false positives. However, given the extremely high confirmation rate of the candidates tested, it is conceivable that high proportions of the remaining candidates indeed respond to oxygen deprivation and represent bona fide HRMs.
The direct contribution of HIFs to the upregulation of HRMs was dissected for select microRNAs using a combination of luciferase-based reporters (containing fragments of microRNA promoters) and chromatin immunoprecipitation. A large variety of direct HIF target genes have been reported, and an indirect, microRNA-mediated component could further add to the complexity of the molecular response orchestrated by these transcription factors. Several in silico methods for target gene prediction have been developed and are publicly available, including PicTar (pictar.bio.nyu.edu), TargetScan 3.0. (http://www.targetscan.org/), and miRBase (http://microrna.sanger.ac.uk/sequences/). For a given microRNA, these utilize different algorithms and ranking criteria (8, 18, 19) and are known to produce a partially overlapping set of candidates, making the search for targets a complex endeavor. Using these three programs, we performed in silico searches for HRM targets, which revealed a highly complex spectrum, including genes involved in apoptosis and proliferation. This could indeed be highly significant for the response to low oxygen, since hypoxia is known to have an impact on both processes.
We show experimental evidence that caspase activation is inhibited by several HRMs during hypoxia. Interestingly, and in keeping with this, our searches with the above-mentioned programs predict that several HRMs target core components of apoptosis: BAK1 (miR-26), BIM (miR-24 and -181), BID (miR-23), caspase-7 (miR-23), CASP3 (miR-30), APAF1 (miR-27), and NIX/BNIP3L (miR-23). In the case of miR-26, we demonstrated a direct antiapoptotic effect with oxygen depletion, increasing the possibility that BAK1 (a proapoptotic protein) is a relevant target.
With regard to miR-210, another HRM with antiapoptotic effect, in silico searches did not reveal candidate targets that are part of the apoptotic machinery. However, a large variety or other genes can influence apoptosis in response to specific stresses. For example, one putative target revealed by PicTar is neuronal pentraxin 1 (11), which has been shown to mediate apoptosis in ischemic neurons, and miR-210 could help neutralize such an effect. Whether such a gene could play a role in apoptosis in nonneuronal cells in low oxygen is not known.
Another process known to be affected by hypoxia is proliferation, with many cell types undergoing cell cycle arrest during oxygen deprivation. HRMs could contribute to this process via predicted targets, such as cdc25A (miR-21, miR-103, and miR-107), cyclin D2 (miR-26, miR-103), cyclin E1 (miR-26), cyclin H (miR23), and cdk6 (miR-26, miR-103). Interestingly, and similarly to the case of apoptosis, several cell cycle genes are predicted to be targeted by multiple HRMs, thereby increasing the chance of efficient downregulation.
Since the hypoxic response described in this paper involves a multitude of microRNAs, it is conceivable that manipulation of any individual HRM could fail to fully capture the phenotypic impact of this mechanism in low oxygen. The concerted induction of these HRMs could therefore have a much more robust impact on apoptotic/proliferative behavior when oxygen is limiting for extended periods, such as in cancer. One could speculate that differences between HRM induction in various cell types could contribute to a variability in the response to hypoxia, with important consequences for cancer progression and response to therapy.
Our analysis shows that a surprisingly high proportion of HRMs are overexpressed in human tumors. The alterations of microRNAs in various cancer types is conceivably the sum of a variety of factors (including oncogene signaling, paracrine factors, and pH alterations), but hypoxia could have a significant impact, setting in motion microtranscripts with biological impact on survival and/or proliferation.
We declare that we have no competing financial interests.
Supplemental material for this article may be found at http://mcb.asm.org/. ![]()
Published ahead of print on 28 December 2006. ![]()
|
|
|---|
) and HIF-2
in stem cells. Mol. Cell. Biol. 26:3514-3526.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»