- Open Access
Phenotypes and gene expression profiles of Saccharopolyspora erythraea rifampicin-resistant (rif) mutants affected in erythromycin production
- Elisabetta Carata†1,
- Clelia Peano†2,
- Salvatore M Tredici1,
- Francesco Ferrari3,
- Adelfia Talà1,
- Giorgio Corti2,
- Silvio Bicciato4,
- Gianluca De Bellis2 and
- Pietro Alifano1Email author
© Carata et al; licensee BioMed Central Ltd. 2009
Received: 31 January 2009
Accepted: 30 March 2009
Published: 30 March 2009
There is evidence from previous works that bacterial secondary metabolism may be stimulated by genetic manipulation of RNA polymerase (RNAP). In this study we have used rifampicin selection as a strategy to genetically improve the erythromycin producer Saccharopolyspora erythraea.
Spontaneous rifampicin-resistant (rif) mutants were isolated from the parental strain NRRL2338 and two rif mutations mapping within rpoB, S444F and Q426R, were characterized. With respect to the parental strain, S444F mutants exhibited higher respiratory performance and up to four-fold higher final erythromycin yields; in contrast, Q426R mutants were slow-growing, developmental-defective and severely impaired in erythromycin production. DNA microarray analysis demonstrated that these rif mutations deeply changed the transcriptional profile of S. erythraea. The expression of genes coding for key enzymes of carbon (and energy) and nitrogen central metabolism was dramatically altered in turn affecting the flux of metabolites through erythromycin feeder pathways. In particular, the valine catabolic pathway that supplies propionyl-CoA for biosynthesis of the erythromycin precursor 6-deoxyerythronolide B was strongly up-regulated in the S444F mutants, while the expression of the biosynthetic gene cluster of erythromycin (ery) was not significantly affected. In contrast, the ery cluster was down-regulated (<2-fold) in the Q426R mutants. These strains also exhibited an impressive stimulation of the nitrogen regulon, which may contribute to lower erythromycin yields as erythromycin production was strongly inhibited by ammonium.
Rifampicin selection is a simple and reliable tool to investigate novel links between primary and secondary metabolism and morphological differentiation in S. erythraea and to improve erythromycin production. At the same time genome-wide analysis of expression profiles using DNA microarrays allowed information to be gained about the mechanisms underlying the stimulatory/inhibitory effects of the rif mutations on erythromycin production.
Actinomycetes are used for fermentative production of a wide range of bioactive molecules including antibiotics, anticancer agents and immune-suppressants. A crucial point of this process is that these microorganisms must often be genetically improved for higher production before they can be used in the industry. Historically, strain improvement has been carried out by multiple rounds of random mutagenesis and selection in close association with process improvements to optimize large-scale industrial fermentations . Since the late 1970s, the availability of molecular genetics tools and useful information about the biosynthetic pathways and genetic control for most of secondary metabolites of commercial interest has opened the way for improving strains by rational engineering [2, 3]. More recently, these rational strain improvement strategies benefit of the support of genomic, transcriptomic, proteomic, and metabolomic technologies [4–11].
The erythromycin fermentation is a classic antibiotic fermentation that has been improved by the traditional mutate-and-screen method over the past 50 years. Erythromycin biosynthesis in the mycelial actinomycete, Saccharopolyspora erythraea, has been widely studied as a model system for antibiotic production [12–15], and erythromycin and its semi-synthetic derivatives are widely used in the clinic; therefore, improved producers are still highly sought after.
Erythromycin A is made by a three-stage pathway : assembly of the 14-membered macrolactone 6-deoxyerythronolide B (6DEB) from one propionyl-CoA and six (2S)-methylmalonyl-CoA units followed by its hydroxylation to erythronolide B (EB), formation of the deoxysugars mycarose and desosamine from glucose and their addition to EB to make erythromycin D, and then C-12 hydroxylation and C-3" O-methylation of the latter compound to produce erythromycin A. The synthesis of 6DEB is catalyzed by multifunctional modular polyketide synthase with at least 28 distinct active sites, in a process similar to that of fatty-acid biosynthesis [17, 18].
Extensive genetic studies have provided some insight into the genes involved in erythromycin biosynthesis [19, 20]. The erythromycin gene cluster contains 20 genes arranged in four major polycistronic units . Evidence for regulatory genes has been missing for a long time hampering efforts to enhance erythromycin production other than by medium manipulation, random mutagenesis and selection. In recent times, the availability of the entire genome sequence of S. erythraea has opened the possibility of defining the mechanisms by which erythromycin is controlled by using global approaches [22, 23]. Very recently, these approaches have led to the discovery that BldD, a key developmental regulator in actinomycetes [24, 25], regulates the synthesis of erythromycin . Meanwhile, there is evidence that increasing the flux through feeder metabolic pathways strongly influences the erythromycin yields. This has been recently obtained by engineering the methylmalonyl-CoA metabolite node in S. erythraea and in Aeromicrobium erythreum, a non-filamentous erythromycin A producer [27–29].
The focus of this study was to explore the possibility to increase the erythromycin production by genetic manipulation of the RNA polymerase (RNAP) of S. erythraea. This working hypothesis relies on well-documented evidence that: i) ppGpp, the effector of the stringent response , triggers antibiotic production in streptomycetes [31–35]; ii) several rifampicin-resistance mutations (rif) in the RNAP beta chain confer ppGpp-independent biosynthesis of the pigmented actinorhodin (Act) and undecylprodigiosin (Red), methylenomycin, and calcium-dependent antibiotic (CDA) in the model actinomycete Streptomyces coelicolor A3(2) ; rif mutations also activate cryptic antibiotic biosynthesis in Streptomyces lividans, a fast-growing close relative of S. coelicolor A3(2), which produces less or no Act, Red and CDA despite the existence of all the required biosynthetic genes [37–41]. Mechanistically, it has been proposed that several rif mutations mimic binding of ppGpp to RNAP [36, 40].
These premises prompted us to investigate: i) the role of the stringent response in the activation of erythromycin biosynthesis; ii) the utility of the rif screening to search for erythromycin over-producing mutants; iii) how transcriptional changes, analyzed at a whole genome level, could shed light on the molecular mechanisms underlying, on one hand, erythromycin over-production by rif1 mutants, and on the other the almost null antibiotic production by rif6 mutants.
Results and Discussion
Isolation and phenotypes of S. erythraea rif mutants
Rif6 and rif12 mutants were conditionally defective in antibiotic production. Indeed, their defect was less apparent when the strains were grown on complex YS solid medium (Figure 1B and data not shown). It should be noted, however, that the erythromycin production by the wild type strain and by rif1 and rif9 mutants was significantly lower on YS than on R3/1. Phenotypical analysis demonstrated that rif6 and rif12 were severely defective in aerial mycelium development and spore formation on R3/1 (Figure 1C and data not shown). In contrast, on the same medium rif1 and rif9 exhibited luxurious growth and more abundant spores with respect to the wild type.
Genotypic characterization of the S. erythraea rif mutants
In prokaryotes more than 90% of rifampicin-resistant isolates have missense mutations, deletions or insertions in the 81-bp rifampicin resistance-determining region (RRDR) of the rpoB gene [42–44]. To determine the location and nature of the S. erythraea rif mutations, this region was amplified by PCR from the wild type and mutant strains and subjected to nucleotide sequencing. The sequence analysis demonstrated: i) a C to T transition at position 1333 of the nucleotide sequence of the coding region of S. erythraea rpoB in the rif1 and rif9 mutants; ii) a C to T transition at position 1326 in the rif2, rif3, rif4 and rif5 mutants; iii) an A to G transition at position 1279 in the rif6 and rif12 mutant strains.
Stimulation of erythromycin production in the rif1 and rif9 mutants (harboring the S444F missense) was thus consistent with the results obtained in S. lividans. Indeed, phenylalanine and tyrosine (replacing, respectively, a serine residue in S. erythraea and an asparagine residue in S. lividans) are structurally similar amino acids. In contrast, the R442W missense was not functionally equivalent to R440H or R440C in S. lividans and thus failed to stimulate erythromycin production. Interestingly, the missense Q426R in the slow-growing and developmental-defective rif6 and rif12 mutants affected a glutamine residue that, when mutated to proline, leucine or arginine (Figure 3) was responsible for transcription termination defects and slow growth in E. coli [42–46].
Microarray analysis of the transcriptome of the S. erythraea rif mutants
To gain information about the mechanisms underlying the stimulatory/inhibitory effects of the rif mutations on erythromycin production, genome-wide analysis of expression profiles using DNA microarrays was performed. To this purpose, the wild type strain and the rif1 (S444F) and rif6 (Q426R) mutants were grown in shake-flasks containing R3/1 medium to either phase a or phase b of the growth curve (i.e. 24 h for wt and rif1 and 48 h for rif6; Figure 2). RNA samples were extracted from two independent cultures, processed and hybridized to custom made GeneChips containing DNA oligonucleotide probes corresponding to all the predicted S. erythaea ORFs.
Expression data of the wild type, rif1 and rif6 mutants during growth phases a and b were compared using Significance Analysis of Microarray (SAM) multiclass analysis. Setting the q-value threshold at 1% allowed identifying 198 and 270 differentially expressed genes (DEG) among wild type, rif1 and rif6 strains in phases a and b, respectively (see Additional file 1, Table S1 and Additional file 2, Figure S1). Among the 198 DEG characterizing the phase a , the five most represented gene-functional classes showing a significant enrichment were: II.6-Posttranslational modification, protein turnover, chaperone; II.8-RNA processing and modification; II.12-Translation, ribosomal structure and biogenesis; III.5-Energy production and conversion; III.8-Nucleotide transport and metabolism. The gene-functional classes showing a significant depletion were: I.3-Cell wall/membrane/envelope biogenesis, III.2-Carbohydrate transport and metabolism.
Among the 270 DEG characterizing the phase b the only functional category showing a significant enrichment was II.12-Translation, ribosomal structure and biogenesis, while a significant depletion was evidenced only for the functional category III.10-Secondary metabolites biosynthesis, transport catabolism.
This cluster includes also many genes involved in nitrogen metabolism: glnB (SACE 6061), encoding the nitrogen regulatory protein PII, amt (SACE 6062), coding for an ammonium transporter, glnA-1 (SACE 1623), coding for the glutamine synthetase, gudB (SACE 4093), coding for the NAD-specific glutamate dehydrogenase, ureA (SACE 0634) and ureC (SACE 0636), coding for alpha and gamma subunits of the urease respectively, nirB (SACE 3801, SACE 3802) and nirD (SACE 3803), encoding the assimilatory nitrite reductase large and small subunits respectively, and narK coding for a nitrite extrusion protein.
Up-regulation of the nitrogen regulon, which in actinomycetes includes genes involved in ammonium assimilation and supply [47, 48], is noteworthy because demonstrates that the erythromycin production is strongly inhibited by ammonium in the wild type strain (see Additional file 2, Figure S2) and is supported by the evidence that, in Amycolatopsis mediterranei, GlnR links rifamycin biosynthesis to nitrogen metabolism .
This cluster comprises 36 genes that were down-regulated in the rif6 mutant and not affected in the rif1 compared to the wild type (Figure 5, right up-panel). The most important members of this cluster are genes encoding NADH dehydrogenase I (nuoA [SACE 6902] and nioL [SACE 6891]), and two genes involved in nucleotide metabolism (purE [SACE 6485] and pyrH [SACE 6036]). More interestingly, all genes of the ery cluster except eryBI (SACE 0732) showed a similar gene expression profile, being down-regulated in the rif6 mutant and not significantly affected in the rif1 compared to the wild type (Figure 4A, left panel). In contrast, eryBI, which is transcribed monocistronically with the transcript start site facing toward that of eryE , exhibited an opposite behavior, e.g. being up-regulated in the rif6 mutant and not significantly affected in the rif1 compared to the wild type. It is worthwhile noting that, although none of the genes of the ery cluster had a statistically significant differential expression at the stringent q-value threshold of 1%, half of them (i.e., eryCVI, eryCI, eryBI, eryCV, eryBV, eryCII, eryBII, eryG, eryBVI, eryBVII) were statistically significant at the still stringent q-value threshold of 5%. Down-regulation of eryCII in the rif6 mutant was supported by reverse transcriptase real-time PCR analysis (Figure 4B).
Unexpectedly, bldD, coding for a key developmental regulator that seems to regulate the ery cluster positively, exhibited an opposite gene expression pattern with respect to most of ery genes. In particular, similarly to the genes of cluster 1, bldD was up-regulated in the rif6 and not significantly affected in the rif1 compared to the wild type (Figure 4A, left panel and Figure 4B). It should be noted, however, that bldD mRNA levels may not parallel BldD protein levels because the regulation of bldD expression in actinomycetes is very complex and involves different mechanisms including transcription repression by BldD [26, 50] and developmental stage-dependent proteolysis . Thus, it is possible that the increased bldD mRNA levels may reflect an enhanced BldD proteolysis in the rif6 mutant. This hypothesis is supported by up-regulation, in this strain, of many genes which are known to be negatively regulated by BldD in streptomycetes  including bldN (SACE 6952), coding for a developmental regulator, cya (SACE 1989), coding for an adenylate cyclase, a two-component system response regulator (SACE 6712) and a HAD-superfamily subfamily IB, PSPase-like protein (SACE 0365) (Figure 4A, left panel). In addition to a possible effect on the expression of the ery cluster, the altered BldD-dependent regulation may also account for the "bald" phenotype of this mutant (Figure 1C). Future work will be aimed at verifying this hypothesis.
This cluster comprises 27 genes that were moderately down-regulated in the rif1 and strongly down-regulated in the rif6 compared to the expression levels into the wild type (Figure 5, right central-panel). This cluster includes genes coding for DNA gyrase subunits A and B (gyrA [SACE 0009] and gyrB [SACE 0008], respectively), ATP-dependent RNA helicases (SACE 5942 and deaD [SACE 6058]) and a large number of ribosomal proteins (rpsB, rpsD, rplD, rplO, rplW, rpmC and rpmD corresponding to SACE 6038, SACE 6804, SACE 6835, SACE 6817, SACE 6834, SACE 6828 and SACE 6818, respectively).
Apparently, down-regulation of genes encoding ribosomal proteins and up-regulation of genes involved in amino acid biosynthesis (cluster 1, see above) was suggestive for an ability of the rif6 mutations to induce a stringent phenotype as well as other well-characterized rif mutations (see Background). However, it should be pointed out that in Corynebacterium glutamicum up-regulation of ribosomal protein operons was also observed in ppGpp-defective strains following treatment with serine hydroxamate, a serine analogue that competitively binds to the seryl-tRNA synthetase and prevents the seryl-tRNA from being charged thus inducing the stringent response .
As the ery cluster was found to be only slightly down-regulated in the rif6 and not significantly affected in the rif1, it is likely that the different yields of erythromycin in the rif1 and the rif6 mutants may be due, at least in part, to a differential activity of this major feeder pathway (Figure 6B). To validate the microarray analysis and obtain more quantitative data, the amounts of mmsA2- and mmsB-specific transcripts were measured by RT real-time PCR. The results confirmed the up-regulation of the cluster SACE 1456-SACE 1459 in the rif1 mutant (Figure 6C).
Moreover, the data analysis identified additional clusters of contiguous genes sharing similar transcriptional modulation. These groups of physically proximal genes contain, in most cases, members of known or putative operons showing similar expression profile (Additional file 2, Figure S3) including: i.) the str locus containing the S10-spc-alpha ribosomal protein operons ; ii.) the rpsB-tsf-pyrH operon coding for ribosomal protein S2, elongation factor Ts and UMP-kinase [54, 55]; iii.) the nuo operon encoding the energy-generating NADH dehydrogenase complex I ; iv.) the ctaE-qcrCA gene locus coding for cytochrome c oxidase subunit III and for subunits of the ubiquinol c reductase including cytochrome cc and a Rieske Fe-S protein, respectively ; v.) the F0 F1 ATP synthase operon containing the eight genes atpBEFHAGDC ; vi.) the ure operon (ureABCFGD) coding for all subunits of the urease and its accessory proteins ; vii.) a putative nar/nir gene cluster coding for a nitrite extrusion protein (NarK), the catalytic subunit of the assimilatory nitrate reductase (NasA/NarB), the large (NirB) and small (NirD) subunits of NAD(P)H-nitrite reductase and two enzymes involved in the biosynthesis of siroheme (HemD and NirE/SirB), the prostetic group of nitrite reductase ; viii.) the amt-glnB operon coding for the ammonium transporter and the nitrogen regulon regulatory protein PII. Up-regulation of the amt-glnB operon in the rif6 mutant may be due to increased glnR (SACE 7101) mRNA levels. Indeed, there is evidence for a central role of GlnR in regulation of nitrogen metabolism in actinomycetes . The expression profiles of several genes located in these operons were confirmed by RT real-time PCR (Additional file 2, Figure S4).
This study demonstrates the usefulness of the rif screening as a tool to search for higher-producer strains and provides new information about the molecular mechanisms underlying the stimulatory effect of several rif mutations on bacterial secondary metabolism. The rif mutations deeply changed the transcriptional profile of S. erythraea suggesting that their effects on erythromycin biosynthesis go beyond the stringent response, as previously reported (see Introduction). In addition to specific effects on clusters of genes coding for secondary metabolites, the expression of genes coding for key enzymes of the carbon and nitrogen central metabolism was dramatically altered affecting, in turn, energy supply, growth rates and fluxes of metabolites through the erythromycin feeder pathways.
The ery cluster was found to be slightly down-regulated in the hypo-producing rif6 mutant possibly as a consequence of perturbed BldD-dependent regulation. In contrast, the expression of the ery cluster was slightly affected in the rif1 mutant. In this strain, the enhanced activity of an erythromycin feeder pathway may account for the hyper-producing phenotype.
The intimate connection between the erythromycin biosynthesis and the central metabolism is consistent with both genomic and expression data. In general, secondary metabolism is believed dispensable for survival, and most of gene clusters coding for secondary metabolites occupy non-core genomic regions and are maximally expressed during late growth phases. In contrast, the ery cluster maps in the core region of the S. erythraea chromosome, and is transcribed during the middle pseudo-exponential growth phase when the activities of the carbon and nitrogen central metabolic pathways are maximal [22, 23, 60]. As shown in figure 6B, these pathways are strictly connected to erythromycin biosynthesis. Indeed, assembly of the 14-membered macrolactone 6DEB requires one propionyl-CoA and six (2S)-methylmalonyl-CoA units.
The analysis of the hyper-producing mutant rif1 suggests that the valine catabolic pathway, which is strongly up-regulated in this strain, may be a major feeder pathway supplying propionyl-CoA. Propionyl-CoA is then transformed into (2S)-methylmalonyl-CoA biotin-dependent carboxylation that in S. erythraea may be accomplished by at least five isoenzymes [22, 23]. On the other hand, there is experimental evidence that the erythromycin biosynthetic pathway is connected to the Krebs cycle via the methylmalonyl-CoA mutase, an adenosylcobalamin-dependent enzyme that catalyzes the reversible isomerization of (2R)-methylmalonyl-CoA and succinyl-CoA [27–29]. This connection may underlie the negative regulation of erythromycin production by ammonium in starch-based fermentations accounting, at least in part, for the phenotype of the hypo-producing mutant rif6 in which the nitrogen regulon, including genes for ammonium uptake and assimilation, was impressively up-regulated.
Bacterial strains and media
S. erythraea wild type strain NRRL2338 was a gift of S. Donadio (KtedoGen, Milano). This strain has been deposited at the American Type Culture Collection. The strain was stored in 1-ml cryotubes at -80°C as frozen mycelium in YS medium containing 15% glycerol at a biomass concentration of approximately 0.25 g dry cell weight (DCW) ml-1, or at -20°C as spores in 20% glycerol (in distilled water) at a title of approximately 5 × 108 ml-1.
Composition of the media used in this study
Composition (per liter)
5 g peptone, 3 g yeast extract, 1 g MgSO4·7 H2O
Seed medium (SM)
4 g peptone, 4 g yeast extract, 2 g KH2PO4, 4 g K2HPO4, 0.5 g MgSO47H20, 10 g glucose
5 g yeast extract, 0.1 g casamino acids, 3 g L-proline, 10 g MgCl2·6H2O, 4 g CaCl2·2H2O, 0.2 g K2SO4, 0.05 g KH2PO4, 5.6 g TES, 10 g glucose
Medium 266 (YS)
2 g yeast extract, 10 g soluble starch
40 g oatmeal, 1 g yeast extract
3 g beef extract, 5 g tryptone, 15 g NaCl
7 g NH4Cl, 3 g KH2PO4, 7 g K2HPO4, 0.25 g MgSO4·7 H2O, 0.0138 g CaCl2·2 H2O, 10 g glucose, 2 ml trace solution element (TSE)a
1 g NH4Cl, 1 g CaCO3, 0.5 g NaCl, 0.4 g MgSO4·7 H2O, 0.15 g KH2PO4, 0.35 g K2HPO4, 10 g glucose, 2 ml trace solution element (TSE)a
Escherichia coli strain DH5α was used in cloning procedures. This strain was grown in Luria Bertani (LB) medium. To allow plasmid selection, LB medium was supplemented with ampicillin (50 μg ml-1). The composition (per liter) of the nutrient broth agar in the microbiological assays with Micrococcus luteus tester strain was: 3 g beef extract, 5 g tryptone, 15 g NaCl, 15 g agar.
Preparation of spores
Concentrated spore suspensions (5 × 108 ml-1) are crucial for purposes like starting reproducible cultures for physiological or fermentation studies. To prepare spores adapted to the conditions of liquid medium, spores were spread on the same medium with agar. Mycelium with spores was strongly attached to the surface agar, which made impossible to collect spores without agar traces. Therefore, the method of growing strains on cellophane disc was used . The cellophane discs were sterilized in distilled water and then placed on agar, and the inoculum was spread on cellophane by glass stick. After two weeks, spores (control in microscope) were easily scraped from cellophane, and stored in 20% glycerol at -20°C.
For shake-flask experiments, spores in frozen aliquots were collected by centrifugation, re-suspended in medium 707 (for rehydration), and readily separated by vortexing. Individual aliquots (about 5 × 108 spores) were used to inoculate each 500 ml buffled Erlenmeyer flask containing 50 ml of the liquid media described above. Cultures were incubated at 30°C with shaking at 250 rpm. Bioreactor cultures were carried out by using Minifors mini-fermenters (Infors AG, Bottmingen, CH) that operated with a working volume of 1.5 l. Stirring was provided by Rushton-type impellors rotating at 250 rpm. Sterile air was supplied through a sparger. The bioreactors were equipped with pH electrode, pO2 electrode (polarographic), antifoam probe and Pt-100. Glucose concentration was monitored during fermentation by coupled glucose oxidase-peroxidase reaction using a commercial kit (distributed by Laboser srl).
Erythromycin production in solid media was assayed by bioassay. To this purpose, S. erythraea strains were grown in solid media (30 ml) in Petri dishes (8.5 cm). After desired time of cultivation, 1.6 cm (diameter) agar discs (with mycelium on the surface) were removed and placed into empty Petri dishes (diameter 8.5 cm) that were filled with soft nutrient agar that was seeded with Micrococcus luteus. Diameters of the zone of inhibition were measured after 2 days incubation at 37°C. Agar discs containing defined amounts of >95% pure erythromycin A (Sigma) were used as a reference. In liquid media, erythromycin was extracted and assayed by both microbiological assay and thin layer chromatography (TLC). Extraction of erythromycin from fermentation broths was performed as described . TLC identification was carried out on silica-gel GF254 plates as previously described .
High molecular weight genomic DNA was extracted from S. erythraea strains grown in 50 ml of SM medium with shaking at 28°C for 5 days (120 h). After centrifugation, the mycelium was re-suspended in 10 ml SET buffer (75 mM NaCl, 25 mM EDTA, 20 mM Tris-Cl pH 7.5) and incubated in the presence of 5 mg ml-1 lysozyme for 30' at 37°C. Samples were sonicated (Sonifer sonicator Model 250/240, Brain Ultrasonic Corporation) 3 times 30 sec, and incubated in the presence of 20 mg ml-1 Proteinase K and 1.2% sodium dodecyl sulfate (SDS) for 2 h at 55°C. Nucleic acids were extracted by fenol-chloroform:isoamylic alcohol (24:1) extraction according to standard procedure , and 15 μg ml-1 ribonuclease A was used to remove RNA. After fenol-chloroform:isoamylic alcohol (24:1) extraction and ethanol-precipitation, high molecular-weight DNA was collected by spooling using Shepherd's crooks .
Oligonucleotides used in RT real-time PCR experiments
DNA sequencing was performed as a service by MWG Biotech. DNA similarity searches were carried out using BLAST at NCBI http://www.ncbi.nlm.nih.gov/. Sequence alignments were performed with Clustal W at EBI http://www.ebi.ac.uk/.
RNA extraction, microarray and reverse transcriptase real-time PCR experiments
For each strain and time point, total RNA was extracted from mycelium pellets deriving from 1-ml culture samples using the GeneElute™ total RNA Purification Kit (SIGMA), recovering it in 50 μl of Elution Solution. After extraction RNAs were quantified with a NanoDrop spectrophotometer (NanoDrop Technologies) and analyzed by capillary electrophoresis on an Agilent Bioanalyzer (Agilent).
The RNA samples showing an RNA Integrity Number (i.e., the quality parameter calculated by the Bioanalyzer software) higher than 7 were processed for microarray hybridization, following the instructions of "Target Labeling for Prokaryotic GeneChip® Antisense Arrays" (Affymetrix Prokaryotic gene Expression Manual). The protocol consists in cDNA synthesis by reverse transcription (starting with 10 μg RNA), followed by cDNA fragmentation with DNase I and labeling with Terminal Deoxynucleotidyl Transferase. The labeled cDNAs were then hybridized for 16 h at 50°C on individual S. erythraea GeneChips. After hybridization, GeneChips were washed and stained with streptavidin-conjugated phycoerythrin by using the Fluidic Station FS450 (Affymetrix) according to the ProkGE-WS2v3_450 Protocol. Fluorescent images of the microarrays were acquired using a GeneChip Scanner 3000 (Affymetrix). All Chip images and files have been deposited in the GEO (Gene Expression Omnibus) http://www.ncbi.nlm.nih.gov/geo/ repository (accession number: GSE12017).
The S. erythraea GeneChip was previously described in Peano et al. . The microarray, targeting the whole set of S. erythraea genes, is composed of 25-mer oligonucleotide probes organized in 7060 probe-sets. The vast majority of probe-sets (i.e. 6494) are composed of at least 11 probes pairs (Perfect Match/MisMatch); among these, 29 probe-sets, targeting genes longer than 5000 bp, include from 13 to 86 probe pairs; the remaining 566 probe-sets are instead composed of less than 11 probe pairs.
Quantitative analysis of eryCII, bldD, SACE 0636, 1456–1459, 1684, 3802, 6038, 6062, 6280, 6818, 6899, 6902 and 7101 specific transcripts, normalized to 16S rRNA, was performed by reverse transcriptase (RT) real-time PCR. Total RNAs (1 μg) from S. erythraea NRRL2338 and several rif derivatives grown in R3/1 medium up to phase a (Figure 2) were reverse-transcribed by using random hexamer (2.5 μM) with Superscript RT (Invitrogen). About 0.1–1% of each RT reaction was used to run real-time PCR on a SmartCycler System (Cepheid) with SYBR® Green JumpStart Taq ReadyMix (Sigma-Aldrich) and the primer pairs indicated in Table 2. Real-time PCR samples were run in triplicate. The real-time PCR conditions were: 10 min initial denaturation at 94°C, 30 sec at 94°C, 30 sec at 55°C, 30 sec at 72°C for 35 cycles; detection of PCR products was performed at 83°C.
The quality of the raw data obtained from microarray hybridization was assessed by MAS5.0 (Microarray Suite/Software, Affymetrix) control parameters after a global scaling at a target intensity of 100. Control parameters, as well as box plot of raw intensities, indicated the overall high quality of the data set and the absence of any outlying sample. Probe level data was converted to expression values using the Robust Multi-array Average (RMA) procedure . Briefly, PM values (Perfect Match) were background-adjusted, normalized using quantile normalization, and log transformed.
Significance Analysis of Microarray (SAM) has been applied to detect differentially expressed genes.
SAM was introduced by Tusher  as a statistical technique for finding significant genes in microarrays while controlling the False Discovery Rate (FDR). SAM computes a statistic for each gene, measuring the strength of the relationship between gene expression and the response variable (e.g. the strain types, the growth phases, etc.). It uses repeated permutations of the data to determine if the expression level of any genes is significantly related to the response. The cutoff for significance is determined controlling the q-value, i.e. the lowest False Discovery Rate at which a gene is called significant . Similarly to the p-value, the q-value measures how significant a gene is differentially expressed in the context of a large number of genes.
Hierarchical clustering and Eisen's maps were used to group modulated genes and samples in the software package dChip . Before clustering, the expression values for a gene across all samples were standardized and these standardized values were used to calculate correlations between genes and samples and served as the basis for merging nodes. In the hierarchical agglomerative clustering Pearson correlation coefficient and centroid were used as distance metric and linkage method, respectively.
This work was partially supported by grant from Progetto MIUR Cofin 2006 to P.A.
- Adrio JL, Demain AL: Genetic improvement of processes yielding microbial products. FEMS Microbiol Rev. 2006, 30 (2): 187-214. 10.1111/j.1574-6976.2005.00009.x.View ArticleGoogle Scholar
- Baltz RH: Genetic methods and strategies for secondary metabolite yield improvement in actinomycetes. Antonie Van Leeuwenhoek. 2001, 79 (3–4): 251-259. 10.1023/A:1012020918624.View ArticleGoogle Scholar
- Vinci VA, Byng G: Manual of Industrial Microbiology and Biotechnology, second ed. Edited by: Demain AL, Davies JE, Atlas RM, Cohen G, Hershberger CL, Hu W-S, Sherman DH, Willson RC, Wu JHD. 1999, 103-113. Washington, DC: American Society for MicrobiologyGoogle Scholar
- Blanchard AP, Hood L: Sequence to array: probing the genome's secrets. Nat Biotechnol. 1996, 14 (13): 1649- 10.1038/nbt1296-1649.View ArticleGoogle Scholar
- Bro C, Nielsen J: Impact of 'ome' analyses on inverse metabolic engineering. Metab Eng. 2004, 6 (3): 204-211. 10.1016/j.ymben.2003.11.005.View ArticleGoogle Scholar
- Donadio S, Monciardini P, Sosio M: Polyketide synthases and nonribosomal peptide synthetases: the emerging view from bacterial genomics. Nat Prod Rep. 2007, 24 (5): 1073-1109. 10.1039/b514050c.View ArticleGoogle Scholar
- Lander ES: The new genomics: global views of biology. Science. 1996, 274 (5287): 536-539. 10.1126/science.274.5287.536.View ArticleGoogle Scholar
- Lander ES: Array of hope. Nat Genet. 1999, 21 (1 Suppl): 3-4. 10.1038/4427.View ArticleGoogle Scholar
- Lee SY, Lee DY, Kim TY: Systems biotechnology for strain improvement. Trends Biotechnol. 2005, 23 (7): 349-358. 10.1016/j.tibtech.2005.05.003.View ArticleGoogle Scholar
- Stafford DE, Stephanopoulos G: Metabolic engineering as an integrating platform for strain development. Curr Opin Microbiol. 2001, 4 (3): 336-340. 10.1016/S1369-5274(00)00214-9.View ArticleGoogle Scholar
- Werf van der MJ, Jellema RH, Hankemeier T: Microbial metabolomics: replacing trial-and-error by the unbiased selection and ranking of targets. J Ind Microbiol Biotechnol. 2005, 32 (6): 234-252. 10.1007/s10295-005-0231-4.View ArticleGoogle Scholar
- Katz L, Donadio S: Macrolides. Biotechnology. 1995, 28: 385-420.Google Scholar
- Leadlay PF: Combinatorial approaches to polyketide biosynthesis. Curr Opin Chem Biol. 1997, 1 (2): 162-168. 10.1016/S1367-5931(97)80005-1.View ArticleGoogle Scholar
- Hutchinson CR, Decker H, Madduri K, Otten SL, Tang L: Genetic control of polyketide biosynthesis in the genus Streptomyces. Antonie Van Leeuwenhoek. 1993, 64 (2): 165-176. 10.1007/BF00873025.View ArticleGoogle Scholar
- McDaniel R, Licari P, Khosla C: Process development and metabolic engineering for the overproduction of natural and unnatural polyketides. Adv Biochem Eng Biotechnol. 2001, 73: 31-52.Google Scholar
- Corcoran JW: Biochemical mechanisms in the biosynthesis of the erythromycins. Antibiotics IV Biosynthesis. Edited by: Corcoran JW. 1981, 132-174. Springer-Verlag edn. New YorkView ArticleGoogle Scholar
- Cortes J, Haydock SF, Roberts GA, Bevitt DJ, Leadlay PF: An unusually large multifunctional polypeptide in the erythromycin-producing polyketide synthase of Saccharopolyspora erythraea. Nature. 1990, 348 (6297): 176-178. 10.1038/348176a0.View ArticleGoogle Scholar
- Donadio S, Staver MJ, McAlpine JB, Swanson SJ, Katz L: Modular organization of genes required for complex polyketide biosynthesis. Science. 1991, 252 (5006): 675-679. 10.1126/science.2024119.View ArticleGoogle Scholar
- Donadio S, McAlpine JB, Sheldon PJ, Jackson M, Katz L: An erythromycin analog produced by reprogramming of polyketide synthesis. Proc Natl Acad Sci USA. 1993, 90 (15): 7119-7123. 10.1073/pnas.90.15.7119.View ArticleGoogle Scholar
- Weber JM, Leung JO, Maine GT, Potenz RH, Paulus TJ, DeWitt JP: Organization of a cluster of erythromycin genes in Saccharopolyspora erythraea. J Bacteriol. 1990, 172 (5): 2372-2383.Google Scholar
- Reeves AR, English RS, Lampel JS, Post DA, Boom Vanden TJ: Transcriptional organization of the erythromycin biosynthetic gene cluster of Saccharopolyspora erythraea. J Bacteriol. 1999, 181 (22): 7098-7106.Google Scholar
- Oliynyk M, Samborskyy M, Lester JB, Mironenko T, Scott N, Dickens S, Haydock SF, Leadlay PF: Complete genome sequence of the erythromycin-producing bacterium Saccharopolyspora erythraea NRRL23338. Nat Biotechnol. 2007, 25 (4): 447-453. 10.1038/nbt1297.View ArticleGoogle Scholar
- Peano C, Bicciato S, Corti G, Ferrari F, Rizzi E, Bonnal RJ, Bordoni R, Albertini A, Bernardi LR, Donadio S, et al: Complete gene expression profiling of Saccharopolyspora erythraea using GeneChip DNA microarrays. Microb Cell Fact. 2007, 6: 37- 10.1186/1475-2859-6-37.View ArticleGoogle Scholar
- Elliot M, Damji F, Passantino R, Chater K, Leskiw B: The bldD gene of Streptomyces coelicolor A3(2). J Bacteriol. 1998, 180 (6): 1549-1555.Google Scholar
- Elliot MA, Bibb MJ, Buttner MJ, Leskiw BK: BldD is a direct regulator of key developmental genes in Streptomyces coelicolor A3(2). Mol Microbiol. 2001, 40 (1): 257-269. 10.1046/j.1365-2958.2001.02387.x.View ArticleGoogle Scholar
- Chng C, Lum AM, Vroom JA, Kao CM: A key developmental regulator controls the synthesis of the antibiotic erythromycin in Saccharopolyspora erythraea. Proc Natl Acad Sci USA. 2008, 105 (32): 11346-11351. 10.1073/pnas.0803622105.View ArticleGoogle Scholar
- Reeves AR, Cernota WH, Brikun IA, Wesley RK, Weber JM: Engineering precursor flow for increased erythromycin production in Aeromicrobium erythreum. Metab Eng. 2004, 6 (4): 300-312. 10.1016/j.ymben.2004.03.003.View ArticleGoogle Scholar
- Reeves AR, Brikun IA, Cernota WH, Leach BI, Gonzalez MC, Weber JM: Effects of methylmalonyl-CoA mutase gene knockouts on erythromycin production in carbohydrate-based and oil-based fermentations of Saccharopolyspora erythraea. J Ind Microbiol Biotechnol. 2006, 33 (7): 600-609. 10.1007/s10295-006-0094-3.View ArticleGoogle Scholar
- Reeves A, Haldenwang WG: Isolation and characterization of dominant mutations in the Bacillus subtilis stressosome components RsbR and RsbS. J Bacteriol. 2007, 189 (5): 1531-1541. 10.1128/JB.01649-06.View ArticleGoogle Scholar
- Cashel M, Gentry DR, Hernandez VJ, Vinella D: The stringent response. Escherichia coli and Salmonella: Cellular and Molecular Biology. Edited by: Neidhart FC, Curtiss III R, Ingraham JL, Lin ECC, Low KB, Magasanik B, Reznikoff WS, Riley M, Schaechter M, Umbarger HE. 1996, 1: 1458-1496. Washington, D.C.: American society for Microbiology, 2Google Scholar
- Artsimovitch I, Patlan V, Sekine S, Vassylyeva MN, Hosaka T, Ochi K, Yokoyama S, Vassylyev DG: Structural basis for transcription regulation by alarmone ppGpp. Cell. 2004, 117 (3): 299-310. 10.1016/S0092-8674(04)00401-5.View ArticleGoogle Scholar
- Chakraburtty R, Bibb M: The ppGpp synthetase gene (relA) of Streptomyces coelicolor A3(2) plays a conditional role in antibiotic production and morphological differentiation. J Bacteriol. 1997, 179 (18): 5854-5861.Google Scholar
- Martinez-Costa OH, Arias P, Romero NM, Parro V, Mellado RP, Malpartida F: A relA/spoT homologous gene from Streptomyces coelicolor A3(2) controls antibiotic biosynthetic genes. J Biol Chem. 1996, 271 (18): 10627-10634. 10.1074/jbc.271.18.10627.View ArticleGoogle Scholar
- Ochi K: Metabolic initiation of differentiation and secondary metabolism by Streptomyces griseus: significance of the stringent response (ppGpp) and GTP content in relation to A factor. J Bacteriol. 1987, 169 (8): 3608-3616.Google Scholar
- Ochi K, Zhang D, Kawamoto S, Hesketh A: Molecular and functional analysis of the ribosomal L11 and S12 protein genes (rplK and rpsL) of Streptomyces coelicolor A3(2). Mol Gen Genet. 1997, 256 (5): 488-498.Google Scholar
- Xu J, Tozawa Y, Lai C, Hayashi H, Ochi K: A rifampicin resistance mutation in the rpoB gene confers ppGpp-independent antibiotic production in Streptomyces coelicolor A3(2). Mol Genet Genomics. 2002, 268 (2): 179-189. 10.1007/s00438-002-0730-1.View ArticleGoogle Scholar
- Hu H, Ochi K: Novel approach for improving the productivity of antibiotic-producing strains by inducing combined resistant mutations. Appl Environ Microbiol. 2001, 67 (4): 1885-1892. 10.1128/AEM.67.4.1885-1892.2001.View ArticleGoogle Scholar
- Hu H, Zhang Q, Ochi K: Activation of antibiotic biosynthesis by specified mutations in the rpoB gene (encoding the RNA polymerase beta subunit) of Streptomyces lividans. J Bacteriol. 2002, 184 (14): 3984-3991. 10.1128/JB.184.14.3984-3991.2002.View ArticleGoogle Scholar
- Lai C, Xu J, Tozawa Y, Okamoto-Hosoya Y, Yao X, Ochi K: Genetic and physiological characterization of rpoB mutations that activate antibiotic production in Streptomyces lividans. Microbiology. 2002, 148 (Pt 11): 3365-3373.View ArticleGoogle Scholar
- Tamehiro N, Hosaka T, Xu J, Hu H, Otake N, Ochi K: Innovative approach for improvement of an antibiotic-overproducing industrial strain of Streptomyces albus. Appl Environ Microbiol. 2003, 69 (11): 6412-6417. 10.1128/AEM.69.11.6412-6417.2003.View ArticleGoogle Scholar
- Talà A, Wang G, Zemanova M, Okamoto S, Ochi K, Alifano P: Activation of dormant bacterial genes by Nonomuraea sp. Strain ATCC 39727 mutant-type RNA polymerase. J Bacteriol. 2009, 191 (3): 805-814. 10.1128/JB.01311-08.View ArticleGoogle Scholar
- Jin DJ, Gross CA: Mapping and sequencing of mutations in the Escherichia coli rpoB gene that lead to rifampicin resistance. J Mol Biol. 1988, 202 (1): 45-58. 10.1016/0022-2836(88)90517-7.View ArticleGoogle Scholar
- Ovchinnikov Yu A, Monastyrskaya GS, Gubanov VV, Lipkin VM, Sverdlov ED, Kiver IF, Bass IA, Mindlin SZ, Danilevskaya ON, Khesin RB: Primary structure of Escherichia coli RNA polymerase nucleotide substitution in the beta subunit gene of the rifampicin resistant rpoB255 mutant. Mol Gen Genet. 1981, 184 (3): 536-538. 10.1007/BF00352535.View ArticleGoogle Scholar
- Ovchinnikov YA, Monastyrskaya GS, Guriev SO, Kalinina NF, Sverdlov ED, Gragerov AI, Bass IA, Kiver IF, Moiseyeva EP, Igumnov VN, et al: RNA polymerase rifampicin resistance mutations in Escherichia coli: sequence changes and dominance. Mol Gen Genet. 1983, 190 (2): 344-348. 10.1007/BF00330662.View ArticleGoogle Scholar
- Landick R, Stewart J, Lee DN: Amino acid changes in conserved regions of the beta-subunit of Escherichia coli RNA polymerase alter transcription pausing and termination. Genes Dev. 1990, 4 (9): 1623-1636. 10.1101/gad.4.9.1623.View ArticleGoogle Scholar
- Severinov K, Soushko M, Goldfarb A, Nikiforov V: Rifampicin region revisited. New rifampicin-resistant and streptolydigin-resistant mutants in the beta subunit of Escherichia coli RNA polymerase. J Biol Chem. 1993, 268 (20): 14820-14825.Google Scholar
- Fink D, Weissschuh N, Reuther J, Wohlleben W, Engels A: Two transcriptional regulators GlnR and GlnRII are involved in regulation of nitrogen metabolism in Streptomyces coelicolor A3(2). Mol Microbiol. 2002, 46 (2): 331-347. 10.1046/j.1365-2958.2002.03150.x.View ArticleGoogle Scholar
- Tiffert Y, Supra P, Wurm R, Wohlleben W, Wagner R, Reuther J: The Streptomyces coelicolor GlnR regulon: identification of new GlnR targets and evidence for a central role of GlnR in nitrogen metabolism in actinomycetes. Mol Microbiol. 2008, 67 (4): 861-880.View ArticleGoogle Scholar
- Yu H, Yao Y, Liu Y, Jiao R, Jiang W, Zhao GP: A complex role of Amycolatopsis mediterranei GlnR in nitrogen metabolism and related antibiotics production. Arch Microbiol. 2007, 188 (1): 89-96. 10.1007/s00203-007-0228-7.View ArticleGoogle Scholar
- Elliot MA, Leskiw BK: The BldD protein from Streptomyces coelicolor is a DNA-binding protein. J Bacteriol. 1999, 181 (21): 6832-6835.Google Scholar
- Lee CJ, Won HS, Kim JM, Lee BJ, Kang SO: Molecular domain organization of BldD, an essential transcriptional regulator for developmental process of Streptomyces coelicolor A3(2). Proteins. 2007, 68 (1): 344-352. 10.1002/prot.21338.View ArticleGoogle Scholar
- Brockmann-Gretza O, Kalinowski J: Global gene expression during stringent response in Corynebacterium glutamicum in presence and absence of the rel gene encoding (p)ppGpp synthase. BMC Genomics. 2006, 7: 230- 10.1186/1471-2164-7-230.View ArticleGoogle Scholar
- Coenye T, Vandamme P: Organisation of the S10, spc and alpha ribosomal protein gene clusters in prokaryotic genomes. FEMS Microbiol Lett. 2005, 242 (1): 117-126. 10.1016/j.femsle.2004.10.050.View ArticleGoogle Scholar
- Weidner U, Geier S, Ptock A, Friedrich T, Leif H, Weiss H: The gene locus of the proton-translocating NADH: ubiquinone oxidoreductase in Escherichia coli. Organization of the 14 genes and relationship between the derived proteins and subunits of mitochondrial complex I. J Mol Biol. 1993, 233 (1): 109-122. 10.1006/jmbi.1993.1488.View ArticleGoogle Scholar
- Serina L, Blondin C, Krin E, Sismeiro O, Danchin A, Sakamoto H, Gilles AM, Barzu O: Escherichia coli UMP-kinase, a member of the aspartokinase family, is a hexamer regulated by guanine nucleotides and UTP. Biochemistry. 1995, 34 (15): 5066-5074. 10.1021/bi00015a018.View ArticleGoogle Scholar
- Sone N, Nagata K, Kojima H, Tajima J, Kodera Y, Kanamaru T, Noguchi S, Sakamoto J: A novel hydrophobic diheme c-type cytochrome. Purification from Corynebacterium glutamicum and analysis of the QcrCBA operon encoding three subunit proteins of a putative cytochrome reductase complex. Biochim Biophys Acta. 2001, 1503 (3): 279-290. 10.1016/S0005-2728(00)00205-X.View ArticleGoogle Scholar
- Gaballo A, Abbrescia A, Palese LL, Micelli L, di Summa R, Alifano P, Papa S: Structure and expression of the atp operon coding for F1F0-ATP synthase from the antibiotic-producing actinomycete Nonomuraea sp. ATCC 39727. Res Microbiol. 2006, 157 (7): 675-683. 10.1016/j.resmic.2006.02.005.View ArticleGoogle Scholar
- Mobley HL, Hausinger RP: Microbial ureases: significance, regulation, and molecular characterization. Microbiol Rev. 1989, 53 (1): 85-108.Google Scholar
- Xiong J, Bauer CE, Pancholy A: Insight into the haem d1 biosynthesis pathway in heliobacteria through bioinformatics analysis. Microbiology. 2007, 153 (Pt 10): 3548-3562. 10.1099/mic.0.2007/007930-0.View ArticleGoogle Scholar
- Katz L, Khosla C: Antibiotic production from the ground up. Nat Biotechnol. 2007, 25 (4): 428-429. 10.1038/nbt0407-428.View ArticleGoogle Scholar
- Kieser T, Bibb MJ, Buttner MJ, Chater KF, Hopwood DA: Practical Streptomyces genetics. 2000, Norwich, England: The John Innes Foundation. CrowesGoogle Scholar
- Le QH SL, Shi YH: Extraction of erythromycin from fermentation broth using salt-induced phase separation processes. Sep Purif Technol. 2001, 24: 85-91. 10.1016/S1383-5866(00)00217-3.View ArticleGoogle Scholar
- Hu CQ, Zou WB, Hu WS, Ma XK, Yang MZ, Zhou SL, Sheng JF, Li Y, Cheng SH, Xue J: Establishment of a Fast Chemical Identification System for screening of counterfeit drugs of macrolide antibiotics. J Pharm Biomed Anal. 2006, 40 (1): 68-74. 10.1016/j.jpba.2005.06.026.View ArticleGoogle Scholar
- Sambrook J, Russell DW: Molecular Cloning. A Laboratory Manual. 2001, Cold Spring Harbor, New York: Cold Spring Harbor Laboratory PressGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.View ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680. 10.1093/nar/22.22.4673.View ArticleGoogle Scholar
- Barrett T, Troup DB, Wilhite SE, Ledoux P, Rudnev D, Evangelista C, Kim IF, Soboleva A, Tomashevsky M, Edgar R: NCBI GEO: mining tens of millions of expression profiles – database and tools update. Nucleic Acids Res. 2007, D760-765. 35 DatabaseView ArticleGoogle Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-264. 10.1093/biostatistics/4.2.249.View ArticleGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.View ArticleGoogle Scholar
- Storey JD: A direct approach to false discovery rates. J Roy Stat Soc. 2002, 64: 479-498. 10.1111/1467-9868.00346. 10.1111/1467-9868.00346.View ArticleGoogle Scholar
- Li C, Wong WH: Model-based analysis of oligonucleotide arrays: expression index computation and outlier detection. Proc Natl Acad Sci USA. 2001, 98 (1): 31-36. 10.1073/pnas.011404098.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.