Saccharopolyspora erythraea’sgenome is organised in high-order transcriptional regions mediated by targeted degradation at the metabolic switch
- Esteban Marcellin†1,
- Tim R Mercer†1,
- Cuauhtemoc Licona-Cassani1,
- Robin W Palfreyman1,
- Marcel E Dinger2,
- Jennifer A Steen1,
- John S Mattick2 and
- Lars K Nielsen1Email author
© Marcellin et al.; licensee BioMed Central Ltd. 2013
Received: 4 October 2012
Accepted: 20 December 2012
Published: 16 January 2013
Actinobacteria form a major bacterial phylum that includes numerous human pathogens. Actinobacteria are primary contributors to carbon cycling and also represent a primary source of industrial high value products such as antibiotics and biopesticides. Consistent with other members of the actinobacterial phylum, Saccharopolyspora erythraea undergo a transitional switch. This switch is characterized by numerous metabolic and morphological changes.
We performed RNA sequencing to analyze the transcriptional changes that occur during growth of Saccharopolyspora erythraea in batch culture. By sequencing RNA across the fermentation time course, at a mean coverage of 4000X, we found the vast majority of genes to be prominently expressed, showing that we attained close to saturating sequencing coverage of the transcriptome. During the metabolic switch, global changes in gene expression influence the metabolic machinery of Saccharopolyspora erythraea, resetting an entirely novel gene expression program. After the switch, global changes include the broad repression of half the genes regulated by complex transcriptional mechanisms. Paralogous transposon clusters, delineate these transcriptional programs. The new transcriptional program is orchestrated by a bottleneck event during which mRNA levels are severely restricted by targeted mRNA degradation.
Our results, which attained close to saturating sequencing coverage of the transcriptome, revealed unanticipated transcriptional complexity with almost one third of transcriptional content originating from un-annotated sequences. We showed that the metabolic switch is a sophisticated mechanism of transcriptional regulation capable of resetting and re-synchronizing gene expression programs at extraordinary speed and scale.
Actinobacteria represent one of the most dominant and distinct bacterial phylum. The evolutionary divergence of actinobacteria is an ancient event, reflected in their unique biology. The majority of actinobacteria are found in soil, where they play a vital role in carbon cycling. In human medicine, they comprise a large number of pathogens and provide a major source of antibiotics . They have a large, high GC content genome that equips them for growth in hostile, nutrient poor environments. They are also able to synthesise secondary metabolites that have important commercial and medical significance .
Saccharopolyspora erythraea (S.erythraea) is a model gram-positive filamentous actinomycete that has served for the study of antibiotic production for several decades. The bacteria produces erythromycin A, the first macrolide antibiotic to be discovered, and the backbone for numerous modern antibiotics . S.erythraea exhibits a distinct and complex lifecycle, in which an initial growth phase is followed by a transition period, known as the metabolic switch, followed by a secondary growth phase. Like in Streptomyces coelicolor (S.coelicolor), the switch in S.erythraea is followed by morphological changes that coincide with potential cell death and the transcription of secondary metabolite gene clusters . A detailed study of transcriptional activity during this transitional period is of major importance to understand the complex, yet little understood, life cycle of this sophisticated bacterial phylum. Studies have profiled actinomycete transcriptomes using microarrays [5, 6]. In fact all annotated coding sequences have been profiled using Affymetrix arrays in S.coelicolor. However, using long probes which are less likely to produce false negatives has proved challenging due to the high GC content of actinobacterial genomes. The recent advent of RNA sequencing provides a new opportunity to profile the global transcriptome of high GC content microorganisms.
Here, we present a detailed transcriptional description of S.erythraea life cycle from in-depth RNA sequencing. Our study shows that the major morphological changes observed at the metabolic switch are accompanied by an extreme transcriptional phenomenon, with the repression of large regions of the genome. In addition, we characterise a restrictive bottleneck event mediated by targeted RNase activity that results in a wholesale and rapid global transcriptional rearrangement event.
RNA sequencing shows pervasive transcription of the genome
Like most actinobacteria, S.erythraea has a 71.1% GC content genome which prevented us designing and performing microarray experiments that contained all ORFs (Additional file 1: Figure S1). After performing microarray experiments in a control sample, we employed RNA sequencing to profile the whole transcriptome. We first assessed the ability of two alternative RNA sequencing approaches, duplex-specific thermostable normalisation (DSN)  and MicrobExpress depletion , to deplete rRNA in GC-rich transcriptomes. Enriched mRNA samples were compared to untreated matched reference controls (Additional file 1: Figure S1B,G,H). We did not observe any systematic divergence in nucleotide composition of sequenced reads between treated and untreated libraries (r 2 = 0.01, Additional file 1: Figure S1B-H). We observed a high correlation between gene abundance estimates by independent validation by qRT-PCR and the expression of low GC content genes represented in the microarray (Additional file 1: Figure S1I,J), thereby confirming the fidelity of both strategies to profile GC rich transcriptomes.
Our comprehensive sequencing showed the S.erythraea genome to be prevalently transcribed, with at least 85.8% of genome sequence from either strand being represented (98.7% of the genome without considering strand). This is similar to recent studies in yeast and other eukaryotes that show pervasive transcription of complex genomes that express intergenic and noncoding RNA [11, 12]. Furthermore, given less than 0.3% of expressed nucleotides were represented singly suggests that close to saturating sequencing coverage of the transcriptome was attained, thereby permitting us to detect rare and transient transcriptional events.
The transcription of genomic macro-regions alternates at the switch
We next applied RNA sequencing to resolve large-scale topological features of the S.erythraea transcriptome. Our sequencing clearly illustrates the strand-specific enrichment for transcription from the leading strand, consistent with co-directional transcription of essential and highly expressed genes with replication  (Figure 1B). Like other actinomycetes, the genome is divided into two large regions, with a core region thought to contain the majority of genes essential for survival, and a non-core region enriched for genes for conditionally adaptive functions [2, 14]. Alignment of sequenced reads across these regions allowed the resolution of genome-scale, contiguous transcriptional regions associated with these core/non-core regions (Figure 1C). These macro-regions undergo alternative enrichment across the switch, with marked repression of the non-core transcriptional loci following the switch (Figure 1C).
The metabolic switch defines a distinct gene expression program
The high temporal resolution and sequencing depth with which we profiled the growth cycle and metabolic switch using RNA-seq (with >4000 fold genome coverage), revealed considerable complexity in gene expression profiles. Although we detected transcription from all genes, for clarity we restricted our analysis to prominently expressed genes (top 50% of the total most prominently expressed genes). Shifts to gene expression invoked at the switch were so broad as to affect 60.5% of the genes that were profiled (>2-fold change). To delineate cogent trends within this complexity, we identified major gene groups that exhibited a normalised enrichment at exponential, transitional or stationary phases across the S.erythraea life cycle (Additional file 1: Figure S3A-C, S4A-C). A comparative analysis delineated two distinct gene expression programs with diverging functional trends before and following the switch (Additional file 2: Table S2). For example, the greater number of genes expressed prior to the switch were enriched for programs associated with central carbon metabolism, such as those involved in carbohydrate synthesis (p = 1.75x10-16), cell wall biogenesis (p = 0.043), energy production (p = 6.80x10-4) and lipid metabolism (p = 0.011). In contrast, genes associated with secretion and transport (p = 0.037), posttranslational modification (p = 8.98x10-6) and ribosomal structure and biogenesis (p = 1.12x10-16) were enriched after the switch. Similarly, the erythromycin gene cluster was enriched after the metabolic switch (Additional file 1: Figure S3). During the metabolic switch the greatest enrichment was for genes associated with signal transduction (p = 0.0011) and transcription regulation (p = 1.94x10-9). These alternate programs were reciprocally enriched for core/non-core and leading/lagging strand organisation (Figure 1D), reflecting the overarching context these macro-regions provide for gene expression. Interestingly we found similar patterns of up or down regulation to previously reported homologous genes in closely related microorganisms [7, 15], suggesting that the transcriptional phenomenon described here comprise a characteristic event that punctuates and defines gene expression in related microorganisms (Additional file 1: Figure S8).
RNA sequencing facilitates operon annotation
By the contiguous transcription across shared and co-expressed genes, we looked at the predicted organization of the 1,134 previously predicted operons that encompass 4,919 genes . Our de novo assembly using Oases, from the strand specific RNA sequencing after trimming all reads to 70 bases and assembling them using a 45 bases k-mer, resulted in 9,170 contigs. We found several of these contigs aligned with previously predicted operons (Additional file 2: Table S3). Moreover, using the strand specific sequencing, we found 2.9% of total transcription to be antisense to these annotations, with an enrichment of antisense transcription at the 3’ end of the genes. Lastly, we employed a range of metrics, including presence, size and structure to annotate novel independent ncRNAs within intergenic regions. By Combining small RNA sequencing and long RNA-seq data, we identified 190 putative ncRNA from transcribed sections in intergenic genomic regions (Additional file 2: Table S7 and ncRNA track in the genome browser). The novel ncRNAs display distinct CPC scores , dynamic transcriptional pattern and 14 of them displayed a distinct ncRNA secondary structure (Additional file 1: Figure S7; Additional file 2: Table S4, S7) .
RNA degradation mediates a transcriptional bottleneck event
Degradation targets specific mRNA classes
Our analysis of RNA integrity across the S.erythraea life cycle revealed extensive RNA degradation coincident with the metabolic switch (Figure 3C and Additional file 1: Figures S4H). Given that samples were treated with RNA-later and that we obtained high RIN numbers (except for the sample corresponding to the metabolic switch), it was unlikely that degradation could occurred post-extraction. Although the Qiagen columns have less binding affinity for RNAs smaller than 200 bp, it does not completely remove sRNAs . Thus, to provide an insight into the degradation pattern, we isolated and sequenced these degraded RNA fractions. We firstly isolated fragments corresponding to 20-50nt in size, thereby depleting the preparation of small regulatory RNAs that were greater than 50nt in length . We also reasoned that, unlike the 5’ monophosphate termini of a fragmented transcript, primary RNA transcripts with a 5’ triphosphate would be refractory to library preparation and sequencing (Additional file 1: Figure S5D) . For comparison we also sequenced matching fractions from before (40 hours) and after (93 hours) the switch, returning approximately 16 million sequenced small RNAs (Additional file 1: Figure S5C). We next omitted any reads that did not align in sense to transcribed regions from matched RNA sequencing libraries, reasoning they may reflect bona-fide small RNAs. Alignment of the remaining small RNA fragments exhibit no apparent size or nucleotide enrichment to annotated coding genes previously resolved.
Our analysis of the S.erythraea transcriptome describes a major transcriptional phenomenon occurring during fermentation in bioreactors. A genome-wide transcriptional non-core region, delineated by genes of markedly lower GC content, sparse gene density and inverted repeats, is broadly repressed following the switch. This region organization is present in related microorganisms, including some with linear genomes, such as the terminal regions of the S.coelicolor linear chromosome where analogous insertion elements mark the inner boundary of terminal inverted repeats . Similar low GC regions have been reported in Streptomyces ambofaciens, where these regions have been associated to genome plasticity and may be used to elucidate mechanisms of evolution . The presence of analogous transposable elements delineating the core/non-core regions is apparent in other related microorganisms anticipating a role in the maintenance and possible regulation of these macro-regions. The alternating expression of these regions, previously observed in S.erythraea using microarrays , is able to mediate a broad coordinated suite of genetic changes, representing a broad architectural mechanism by which a cell can rapidly modulate large changes in the transcriptome in response to environmental cues .
This transcriptional reorganization is complemented by a targeted bottleneck event where mRNAs are degraded. The bottleneck event establishes that the degradation may be targeted to specific mRNA classes, thereby contributing to the stage-specific enrichments following the S.erythraea growth cycle. The importance of these degradation processes in global transcriptomic regulation has been recently dissected in yeast [32, 33] and individual cases described in S.coelicolor, whereby endoribonuclease activity of RNase III has been shown to bind to specific transcript and to regulate antibiotic production [21, 34]. Studies have also described the role of (p)ppGpp in regulating mRNA half-life in S.coelicolor by inhibiting PNPase activity which preserves mRNA for crucial cellular processes [35–37]. Such widespread degradation may also indicate a role for programmed cell death of a distinct cellular subpopulation, in mediating the changes at the metabolic switch. The new transcriptome observed after the switch may be reflective of the transcriptome of the surviving subpopulation. Collectively, these observations of cellular degradation and targeted mRNA degradation supports the emergence of key, yet little studied, mechanism by which gene expression are regulated.
Our results, which attained close to saturating sequencing coverage of the transcriptome, reveal unanticipated additional transcriptional complexity with almost one third of transcriptional content originating from unannotated sequences. Combining RNA sequencing with sequencing of small RNAs showed the importance of non-coding RNAs in bacterial transcription.
Like other related microorganisms, S.erythraea undergoes a distinct transitional switch that is characterized by growth arrest, morphological changes and the production of secondary metabolites. We show this to be accompanied by a massive transcriptional modulation whose speed and scale, revealed here for the first time by RNA sequencing, represents a significant and novel transcriptional phenomenon. Given similarities shared between the life cycle and genome organization with other microorganisms, we suggest that this transcriptional phenomenon may be a widespread feature across the large critically important bacterial group of industrial microorganisms.
Bacterial strain, growth and fermentation conditions
S.erythraea (NRRL2338) was used in this study. All fermentations were conducted in 2L bioreactors (Applikon) in mineral salt medium MM-101 without casaminoacids, as described elsewhere [38, 39]. Medium ISP 2 (yeast extract, 4 g/L; malt extract, 10 g/L; Dextrose, 4 g/L; Agar, 20 g/L) was used as solid media for spore germination and seed cultures. Approximately 0.5 mL of glycerol stock was used to inoculate starting cultures in a 500 mL baffled flask with 100 mL of ISP 2 media incubated at 30°C in a rotary shaker (INFORS HT, Bottmingen, Switzerland) at 220 rpm for 30 h. When the seed culture reached an OD450 of 2.5 (early stationary phase), a second seed culture (1 L baffled flasks with 150 mL of ISP 2) was inoculated to an initial OD450 of 0.3 and incubated under the same culture conditions for 72 h. Cells were then centrifuged at 10,000 rpm at room temperature (Allegra X-15R, Beckman Coulter, USA), washed and resuspended in MM-101 prior to inoculation. Temperature and pH remained constant at 30°C and pH = 7.0 respectively by addition of 20% NaOH or 10.9% HCl. Dissolved oxygen was maintained between 45 and 60% saturation. Oxygen uptake rate (OUR) and CO2 production were measured using a DASGIP Off-Gas Analyzer GA4 (DASGIP AG, Jülich, Germany) or a mass spectrometer (Hiden HPR20-QIC B, Warrington, England) attached to the bioreactor condensers. Cells were harvested from two biological replicates at seven time points for RNA extraction (see below). Erythromycin A was analysed by LCMS as described in .
Total RNA isolation and mRNA enrichment
To protect RNA from endogenous RNA degradation, RNA isolation was performed after 8 h incubation at 4°C in RNA-later (Ambion). Total RNA was extracted using two cycles of cellular lysis in RNAse-free zirconia beads. RNA was purified using RNeasy midi kit (Qiagen) with on column DNase treatment and a second off column DNase treatment. Although the Qiagen columns have less binding affinity for RNAs smaller than 200 bp, it does not completely remove sRNAs  and are therefore suitable for sRNA analysis. RNA quality was evaluated using BioAnalyzer (Agilent) and Nanodrop prior analysis. Ribosomal RNA was removed with MicrobExpress Bacterial mRNA Enrichment kit (Ambion) or duplex-specific thermostable nuclease enzyme from Kamchatka crab (DSN) as described in the Illumina’s protocol. RNA samples quality and integrity was analyzed using 2100 BioAnalyzer (Agilent) according to manufacturer’s instructions.
Microarrays and real-time quantitative PCR
Probes for microarrays, were designed using eArray and OligoRankPick  and evaluated using Agilent’s Base Composition (BC) algorithm. RNA was amplified using MessageAmp II-bacterial RNA Amplification kit (Life Technologies), reversed transcribed to cDNA, labeled using Agilent QuickAmp labeling kit (Life Technologies), hybridized, washed and scanned according to Agilent’s protocols. RT-PCR was performed from cDNA amplified from 600 ng of total RNA using Superscript III (Invitrogen) as per the manufacturer’s instructions for high GC% transcripts (50°C, 1 h). Diluted cDNA (1 in 50) was used for qRT PCRs performed in a Corbett Rotorgene 3000 using SYBRgreen UDG master-mix (Invitrogen). Ct values were determined using the default values on the Rotor Gene- 6 software. A correction for DNA contamination was performed using Ct values from RT- controls. The concentration of cDNA (ng/reaction) for each sample was determined using corrected values. These values were then log2 transformed and normalized against eryBV (Additional file 1: Figure 1-i).
Library preparation and sequencing
Deep sequencing was performed from libraries prepared as follow: Two libraries were prepared for total RNA, seven after DSN treatment  and the second round of sequencing (from an independent fermentation), following mRNA enrichment using MicrobExpress (Ambion), was sequenced strand specifically. Additionally, we ran a trial DSN-sequencing sample sequenced after preparation of three independent 90-nt, paired-end libraries. To achieve optimized clustering density, the trial run was sequenced with two 400mer libraries sequenced at two different concentrations following DSN normalization [41, 42]. The third library was size selected from a gel after the final amplification to remove adapter multimer product and sequenced at 300mer. This optimized protocol for GC-rich microorganisms was used for subsequent DSN-sequencing. Libraries prepared after MicrobExpress mRNA enrichment were directionally sequenced at 90bp as described in the direction mRNA-seq prep pre-realised protocol from Illumina . RNAs corresponding from 15-50nt fractions relative to ladder were excised for from the PAGE gel and purified for sequencing using the Illumina small RNA sequencing protocol with minor modifications as previously described . All sequencing was performed at Geneworks (Adelaide, Australia) on the Illumina GAII. All Data is available at GEO GSE39722.
Sequence read alignment and analysis
Reads were independently aligned to the S.erythraea genome using multiple software to assess for possible alignment bias caused by the high GC content of the genome. Reads were aligned using Burrows-Wheeler Aligner BWA  Bowtie  and Bowtie2  requiring no more than 2 mismatches. Only reads that aligned to one genomic location were retained unless specified for the analysis. We employed CuffDiff  to determine gene expression using SAM files resulting from alignment, annotated gene file upper quantile normalization and masking of rRNA sequences. For expression of non-standard features, such as transcriptional macro-regions, read alignments were summed over genomic annotations. Hierarchal clustering of gene expression was performed by Pearson’s pair-wise complete linkage using Cluster3 . Gene expression was normalized within libraries and between libraries as indicated in . To determine significance of gene ontology enrichments, we performed Fisher’s exact test with multiple hypothesis correction. Additional statistical tests and figure generation were conducted using the Prism 5 (http://www.graphpad.com/prism/). Additional in-house perl scripts provided at http://matticklab.com/index.php?title=Marcel_Dinger was also employed. We also employed for figure generation Circos to generate circular genomic figures . Operon prediction was performed from de novo assembly of all reads resulting from the strand specific biological replicate using Oases . Reads were trimmed to 70 bases and assembled using a 45 bases kmer, from which 23% aligned with previously annotated operons (Additional file 2: Table S3, operon track in Gbrowser and de novo track in Gbrowser).
The authors would like to thank the following funding sources: The University of Queensland (UQ ECR Grant 2010002194); Australian Research Council/University of Queensland co-sponsored Federation Fellowship (FF0561986; to J.S.M. and T.R.M.); Australian National Health and Medical Research Council Australia Fellowship (631668; to J.S.M) and Career Development Award (CDA631542; to M.E.D.); CONACYT scholarship (to C.L.C)
- Bibb MJ: Regulation of secondary metabolism in streptomycetes. Curr Opin Microbiol. 2005, 8 (2): 208-215.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle Scholar
- Wu J, Zhang Q, Deng W, Qian J, Zhang S, Liu W: Toward improvement of erythromycin A production in an industrial Saccharopolyspora erythraea strain via facilitation of genetic manipulation with an artificial attB site for specific recombination. Appl Environ Microbiol. 2011, 77 (21): 7508-7516.PubMed CentralView ArticlePubMedGoogle Scholar
- Manteca A, Sanchez J, Jung HR, Schwammle V, Jensen ON: Quantitative proteomics analysis of Streptomyces coelicolor development demonstrates that onset of secondary metabolism coincides with hypha differentiation. Mol Cell Proteomics. 2010, 9 (7): 1423-1436.PubMed CentralView ArticlePubMedGoogle Scholar
- Nieselt K, Battke F, Herbig A, Bruheim P, Wentzel A, Jakobsen OM, Sletta H, Alam MT, Merlo ME, Moore J, et al: The dynamic architecture of the metabolic switch in Streptomyces coelicolor. BMC Genomics. 2010, 11: 10-PubMed CentralView ArticlePubMedGoogle Scholar
- Peano C, Bicciato S, Corti G, Ferrari F, Rizzi E, Bonnal R, Bordoni R, Albertini A, Bernardi L, Donadio S, et al: Complete gene expression profiling of Saccharopolyspora erythraea using GeneChip DNA microarrays. Microb Cell Fact. 2007, 6 (1): 37-PubMed CentralView ArticlePubMedGoogle Scholar
- Battke F, Herbig A, Wentzel A, Jakobsen OM, Bonin M, Hodgson DA, Wohlleben W, Ellingsen TE, Nieselt K: A technical platform for generating reproducible expression data from Streptomyces coelicolor batch cultivations. Adv Exp Med Biol. 2011, 696: 3-15.View ArticlePubMedGoogle Scholar
- Zhulidov PA, Bogdanova EA, Shcheglov AS, Vagner LL, Khaspekov GL, Kozhemyako VB, Matz MV, Meleshkevitch E, Moroz LL, Lukyanov SA, et al: Simple cDNA normalization using kamchatka crab duplex-specific nuclease. Nucleic Acids Res. 2004, 32 (3): e37-PubMed CentralView ArticlePubMedGoogle Scholar
- He S, Wurtzel O, Singh K, Froula JL, Yilmaz S, Tringe SG, Wang Z, Chen F, Lindquist EA, Sorek R, et al: Validation of two ribosomal RNA removal methods for microbial metatranscriptomics. Nat Meth. 2010, 7 (10): 807-812.View ArticleGoogle Scholar
- Shagina I, Bogdanova E, Mamedov IZ, Lebedev Y, Lukyanov S, Shagin D: Normalization of genomic DNA using duplex-specific nuclease. Biotechniques. 2010, 48 (6): 455-459.View ArticlePubMedGoogle Scholar
- Kapranov P, Cheng J, Dike S, Nix DA, Duttagupta R, Willingham AT, Stadler PF, Hertel J, Hackermuller J, Hofacker IL, et al: RNA Maps Reveal New RNA Classes and a Possible Function for Pervasive Transcription. Science. 2007, 316 (5830): 1484-1488.View ArticlePubMedGoogle Scholar
- Nagalakshmi U, Wang Z, Waern K, Shou C, Raha D, Gerstein M, Snyder M: The Transcriptional Landscape of the Yeast Genome Defined by RNA Sequencing. Science. 2008, 320 (5881): 1344-1349.PubMed CentralView ArticlePubMedGoogle Scholar
- Brewer BJ: When polymerases collide: Replication and the transcriptional organization of the E. coli chromosome. Cell. 1988, 53 (5): 679-686.View ArticlePubMedGoogle Scholar
- Hopwood DA: Soil to genomics: the Streptomyces chromosome. Annu Rev Genet. 2006, 40: 1-23.View ArticlePubMedGoogle Scholar
- Martin JF, Santos-Beneit F, Rodriguez-Garcia A, Sola-Landa A, Smith MC, Ellingsen TE, Nieselt K, Burroughs NJ, Wellington EM: Transcriptomic studies of phosphate control of primary and secondary metabolism in Streptomyces coelicolor. Appl Microbiol Biotechnol. 2012, 95 (1): 61-75.View ArticlePubMedGoogle Scholar
- Mao F, Dam P, Chou J, Olman V, Xu Y: DOOR: a database for prokaryotic operons. Nucleic Acids Res. 2009, 37 (suppl 1): D459-D463.PubMed CentralView ArticlePubMedGoogle Scholar
- Schulz MH, Zerbino DR, Vingron M, Birney E: Oases: robust de novo RNA-seq assembly across the dynamic range of expression levels. Bioinformatics. 2012, 28 (8): 1086-1092.PubMed CentralView ArticlePubMedGoogle Scholar
- Kong L, Zhang Y, Ye Z-Q, Liu X-Q, Zhao S-Q, Wei L, Gao G: CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007, 35 (suppl 2): W345-W349.PubMed CentralView ArticlePubMedGoogle Scholar
- Passalacqua KD, Varadarajan A, Ondov BD, Okou DT, Zwick ME, Bergman NH: The Structure and Complexity of a Bacterial Transcriptome. J Bacteriol. 2009, 191 (10): 3203-3211.PubMed CentralView ArticlePubMedGoogle Scholar
- Manteca A, Fernandez M, Sanchez J: Cytological and biochemical evidence for an early cell dismantling event in surface cultures of Streptomyces antibioticus. Res Microbiol. 2006, 157 (2): 143-152.View ArticlePubMedGoogle Scholar
- Gatewood ML, Bralley P, Weil MR, Jones GH: The transcriptome of Streptomyces coelicolor: RNA-seq and RNA immunoprecipitation identify substrates for RNase III. J Bacteriol. 2012, 10.1128/JB.06541-11. Published ahead of printGoogle Scholar
- Jones GH: RNA degradation and the regulation of antibiotic synthesis in Streptomyces. Future Microbiol. 2010, 5 (3): 419-429.View ArticlePubMedGoogle Scholar
- Qiu Y, Cho BK, Park YS, Lovley D, Palsson BO, Zengler K: Structural and operational complexity of the Geobacter sulfurreducens genome. Genome Res. 2010, 20 (9): 1304-1311.PubMed CentralView ArticlePubMedGoogle Scholar
- Hatoum-Aslan A, Maniv I, Marraffini LA: Mature clustered, regularly interspaced, short palindromic repeats RNA (crRNA) length is measured by a ruler mechanism anchored at the precursor processing site. Proc Natl Acad Sci U S A. 2011, 108 (52): 21218-21222.PubMed CentralView ArticlePubMedGoogle Scholar
- Sharma CM, Hoffmann S, Darfeuille F, Reignier J, Findeisz S, Sittka A, Chabas S, Reiche K, Hackermuller J, Reinhardt R, et al: The primary transcriptome of the major human pathogen Helicobacter pylori. Nature. 2010, 464 (7286): 250-255.View ArticlePubMedGoogle Scholar
- McCormick JR, Flärdh K: Signals and regulators that govern Streptomyces development. FEMS Microbiol Rev. 2012, 36: 1-View ArticleGoogle Scholar
- Hindra , Pak P, Elliot MA: Regulation of a novel gene cluster involved in secondary metabolite production in Streptomyces coelicolor. J Bacteriol. 2010, 192 (19): 4973-4982.PubMed CentralView ArticlePubMedGoogle Scholar
- Xu W, Huang J, Lin R, Shi J, Cohen SN: Regulation of morphological differentiation in S. coelicolor by RNase III (AbsB) cleavage of mRNA encoding the AdpA transcription factor. Mol Microbiol. 2010, 75 (3): 781-791.PubMed CentralView ArticlePubMedGoogle Scholar
- Bentley SD, Chater KF, Cerdeno-Tarraga AM, Challis GL, Thomson NR, James KD, Harris DE, Quail MA, Kieser H, Harper D, et al: Complete genome sequence of the model actinomycete Streptomyces coelicolor A3(2). Nature. 2002, 417 (6885): 141-147.View ArticlePubMedGoogle Scholar
- Choulet F, Gallois A, Aigle B, Mangenot S, Gerbaud C, Truong C, Francou F-X, Borges F, Fourrier C, Guerineau M, et al: Intraspecific Variability of the Terminal Inverted Repeats of the Linear Chromosome of Streptomyces ambofaciens. J Bacteriol. 2006, 188 (18): 6599-6610.PubMed CentralView ArticlePubMedGoogle Scholar
- Umbarger MA, Toro E, Wright MA, Porreca GJ, Bau D, Hong SH, Fero MJ, Zhu LJ, Marti-Renom MA, McAdams HH, et al: The three-dimensional architecture of a bacterial genome and its alteration by genetic perturbation. Mol Cell. 2011, 44 (2): 252-264.View ArticlePubMedGoogle Scholar
- Harigaya Y, Parker R: Global analysis of mRNA decay intermediates in Saccharomyces cerevisiae. Proc Natl Acad Sci U S A. 2012, 109 (29): 11764-11769.PubMed CentralView ArticlePubMedGoogle Scholar
- Miller C, Schwalb B, Maier K, Schulz D, Dumcke S, Zacher B, Mayer A, Sydow J, Marcinowski L, Dolken L, et al: Dynamic transcriptome analysis measures rates of mRNA synthesis and decay in yeast. Mol Syst Biol. 2011, 7: 458-PubMed CentralView ArticlePubMedGoogle Scholar
- Gravenbeek ML, Jones GH: The endonuclease activity of RNase III is required for the regulation of antibiotic production by Streptomyces coelicolor. Microbiology. 2008, 154 (Pt 11): 3547-3555.View ArticlePubMedGoogle Scholar
- Dalebroux ZD, Swanson MS: ppGpp magic beyond RNA polymerase. Nat Rev Microbiol. 2012, 10 (3): 203-212.View ArticlePubMedGoogle Scholar
- Gatewood ML, Jones GH: (p)ppGpp inhibits polynucleotide phosphorylase from streptomyces but not from Escherichia coli and increases the stability of bulk mRNA in Streptomyces coelicolor. J Bacteriol. 2010, 192 (17): 4275-4280.PubMed CentralView ArticlePubMedGoogle Scholar
- Siculella L, Damiano F, di Summa R, Tredici SM, Alduina R, Gnoni GV, Alifano P: Guanosine 5'-diphosphate 3'-diphosphate (ppGpp) as a negative modulator of polynucleotide phosphorylase activity in a 'rare' actinomycete. Mol Microbiol. 2010, 77 (3): 716-729.View ArticlePubMedGoogle Scholar
- Carata E, Peano C, Tredici SM, Ferrari F, Tala A, Corti G, Bicciato S, De Bellis G, Alifano P: Phenotypes and gene expression profiles of Saccharopolyspora erythraea rifampicin-resistant (rif) mutants affected in erythromycin production. Microb Cell Fact. 2009, 8: 18-PubMed CentralView ArticlePubMedGoogle Scholar
- Licona-Cassani C, Marcellin E, Quek LE, Jacob S, Nielsen LK: Reconstruction of the Saccharopolyspora erythraea genome-scale model and its use for enhancing erythromycin production. Antonie Van Leeuwenhoek. 2012, 102 (3): 493-502.View ArticlePubMedGoogle Scholar
- Hu G, Llinas M, Li J, Preiser PR, Bozdech Z: Selection of long oligonucleotides for gene expression microarrays using weighted rank-sum strategy. BMC Bioinformatics. 2007, 8: 350-PubMed CentralView ArticlePubMedGoogle Scholar
- DSN Normalization Sample Prep Guide: Pilot Release.Google Scholar
- Storz G, Vogel J, Wassarman Karen M: Regulation by Small RNAs in Bacteria: Expanding Frontiers. Mol Cell. 2011, 43 (6): 880-891.PubMed CentralView ArticlePubMedGoogle Scholar
- Directional mRNA-Seq Library Prep: Pre-Release Protocol.Google Scholar
- Taft RJ, Simons C, Nahkuri S, Oey H, Korbie DJ, Mercer TR, Holst J, Ritchie W, Wong JJ, Rasko JE, et al: Nuclear-localized tiny RNAs are associated with transcription initiation and splice sites in metazoans. Nat Struct Mol Biol. 2010, 17 (8): 1030-1034.View ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760.PubMed CentralView ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg S: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-PubMed CentralView ArticlePubMedGoogle Scholar
- Langmead B, Salzberg SL: Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012, 9 (4): 357-359.PubMed CentralView ArticlePubMedGoogle Scholar
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotech. 2010, 28 (5): 511-515.View ArticleGoogle Scholar
- de Hoon MJL, Imoto S, Nolan J, Miyano S: Open source clustering software. Bioinformatics. 2004, 20 (9): 1453-1454.View ArticlePubMedGoogle Scholar
- Bullard JH, Purdom E, Hansen KD, Dudoit S: Evaluation of statistical methods for normalization and differential expression in mRNA-Seq experiments. BMC Bioinformatics. 2010, 11: 94-PubMed CentralView ArticlePubMedGoogle Scholar
- Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, Jones SJ, Marra MA: Circos: An information aesthetic for comparative genomics. Genome Res. 2009, 19 (9): 1639-1645.PubMed CentralView ArticlePubMedGoogle 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.