Directional RNA deep sequencing sheds new light on the transcriptional response of Anabaena sp. strain PCC 7120 to combined-nitrogen deprivation
© Flaherty et al; licensee BioMed Central Ltd. 2011
Received: 18 February 2011
Accepted: 28 June 2011
Published: 28 June 2011
Cyanobacteria are potential sources of renewable chemicals and biofuels and serve as model organisms for bacterial photosynthesis, nitrogen fixation, and responses to environmental changes. Anabaena (Nostoc) sp. strain PCC 7120 (hereafter Anabaena) is a multicellular filamentous cyanobacterium that can "fix" atmospheric nitrogen into ammonia when grown in the absence of a source of combined nitrogen. Because the nitrogenase enzyme is oxygen sensitive, Anabaena forms specialized cells called heterocysts that create a microoxic environment for nitrogen fixation. We have employed directional RNA-seq to map the Anabaena transcriptome during vegetative cell growth and in response to combined-nitrogen deprivation, which induces filaments to undergo heterocyst development. Our data provide an unprecedented view of transcriptional changes in Anabaena filaments during the induction of heterocyst development and transition to diazotrophic growth.
Using the Illumina short read platform and a directional RNA-seq protocol, we obtained deep sequencing data for RNA extracted from filaments at 0, 6, 12, and 21 hours after the removal of combined nitrogen. The RNA-seq data provided information on transcript abundance and boundaries for the entire transcriptome. From these data, we detected novel antisense transcripts within the UTRs (untranslated regions) and coding regions of key genes involved in heterocyst development, suggesting that antisense RNAs may be important regulators of the nitrogen response. In addition, many 5' UTRs were longer than anticipated, sometimes extending into upstream open reading frames (ORFs), and operons often showed complex structure and regulation. Finally, many genes that had not been previously identified as being involved in heterocyst development showed regulation, providing new candidates for future studies in this model organism.
Directional RNA-seq data were obtained that provide comprehensive mapping of transcript boundaries and abundance for all transcribed RNAs in Anabaena filaments during the response to nitrogen deprivation. We have identified genes and noncoding RNAs that are transcriptionally regulated during heterocyst development. These data provide detailed information on the Anabaena transcriptome as filaments undergo heterocyst development and begin nitrogen fixation.
Cyanobacteria are photosynthetic prokaryotes that have evolved a wide array of metabolic capabilities . Because of their high photosynthetic efficiency, variety of metabolic pathways, and genetic manipulability, they are a potential source of "green" chemicals and fuels [2, 3]. Some cyanobacteria reduce atmospheric nitrogen to ammonia to support growth in nitrogen-deficient conditions . Because nitrogen is often a limiting resource for growth, this gives nitrogen-fixing strains a competitive edge in some environments. Understanding the response to nitrogen deprivation, nitrogen fixation, and diazotrophic growth in cyanobacteria will shed light on basic mechanisms of bacterial genetic regulation and physiology. In addition, it may help to develop better strains of cyanobacteria for the production of renewable chemicals and biofuels.
The cyanobacterium Anabaena (Nostoc) sp. strain PCC 7120 grows as long filaments of photosynthetic vegetative cells in the presence of combined nitrogen. In an environment lacking combined nitrogen, about 7 to 10% of the cells terminally differentiate into nitrogen-fixing heterocysts. Heterocysts provide a microoxic environment for the expression of the oxygen-sensitive nitrogenase enzyme [5, 6]. Single heterocysts are spaced about every 10-15 cells along filaments and they supply fixed nitrogen, probably in the form of amino acids, to neighboring vegetative cells . Vegetative cells provide heterocysts with products of carbon fixation, probably as sucrose [7, 8], thus creating a multicellular organism with two mutually dependent cell types. Heterocyst development involves the response of vegetative cells to nitrogen deprivation, the formation and maintenance of the pattern of the two cell types, differentiation of heterocysts from vegetative cells, and the adaptations made by vegetative cells to adjust to diazotrophic growth.
The differentiation of a vegetative cell into a heterocyst involves substantial changes in cell morphology and physiology [5, 6]. Heterocysts deposit glycolipid and polysaccharide layers outside of their cell wall to limit the entry of atmospheric oxygen [9–11]. They lack photosystem II activity, which normally produces O2, and increase respiration to consume O2 that enters the cell. Heterocyst differentiation requires dramatic changes in transcription and some of the key components of this regulation are known. Nitrogen limitation is sensed by accumulation of 2-oxoglutarate (2-OG), the backbone for nitrogen assimilation. 2-OG enhances the DNA-binding activity of the transcription factor ntcA , which regulates expression of the response regulator nrrA, which is partially responsible for upregulation of hetR [13, 14]. HetR, deemed the master regulator of heterocyst development, regulates the expression of many genes, including the glycolipid genes (hgl), exopolysaccharide genes (hep), and the patS gene, which encodes a peptide involved in heterocyst pattern formation .
Factors other than those described above are known to be involved in heterocyst development and have been identified through microarrays and genetic screens [16–22]. While these methods are powerful, microarrays and screens often overlook unannotated regions of the genome and antisense or noncoding transcripts. In addition, they lack sensitivity and do not provide information on UTR length or operon structure. Therefore, we have employed directional RNA-seq to analyze the transcriptome of Anabaena filaments during nitrogen step-down to identify and map all transcripts during heterocyst development [23–27].
Our RNA-seq data provide information on the UTR lengths of each mRNA transcript, on the transcription of sense and antisense or noncoding RNAs, and on the changes in expression of all transcripts whether or not they carry an annotation. The data show long 5' UTRs for many genes, likely with multiple transcriptional start or processing sites. In addition, our study identifies antisense transcription in the coding region or UTR of many genes known to be involved in heterocyst development. Finally, we detected new genes that are significantly upregulated in response to nitrogen deprivation.
Results and Discussion
Analysis of the transcriptional response to nitrogen step-down
We obtained RNA-seq data for total RNA isolated from Anabaena filaments grown with ammonium as a nitrogen source and at three times after nitrogen step-down from ammonium to dinitrogen in air. The nitrogen step-down produces relatively synchronous induction of heterocyst development. RNA-seq data were acquired from filaments at 0, 6, 12, and 21 hours after nitrogen step-down, which provides detailed transcriptome data at important stages of heterocyst development (GEO accession #GSE26633). At 0 hours, all cells in the filaments are actively growing photosynthetic vegetative cells. At 6 hours, the cells have responded to the nitrogen step-down and are expressing early heterocyst differentiation genes such as hetR, the master regulator of heterocyst development, and patS, which is involved in pattern formation. By 12 hours, proheterocysts are committed to complete differentiation and are expressing genes required for altering their morphology and physiology to become microoxic. By 21 hours, nearly all heterocysts appear fully formed, contain polar cyanophycin granules, and are actively fixing nitrogen.
RNA-seq expression data are presented as RPKM, or r eads p er k ilobase (kb) of CDS (coding sequence) model per m illion mapped reads in the sample , with the CDS model defined as the CDS plus 100 bp of 5' UTR. RPKM values and changes in RPKM value for the chromosome and six plasmids are presented as additional files 1, 2, 3, 4, 5, 6, and 7. These data can be examined and filtered in many ways. For example, for genes on the chromosome (additional file 1: Chromosome.xlsx) with a RPKM value of at least 2 (which includes only those genes with good read coverage) and a fold change of at least 5, there are 22 genes with increased expression by 6 hours after nitrogen step-down, 434 genes upregulated by 12 hours (including many known heterocyst morphogenesis genes), and 396 genes upregulated by 21 hours (including the nitrogen fixation genes). For genes on the chromosome with decreased expression, there are 6 genes downregulated at 6 hours after nitrogen step-down, 32 genes downregulated at 12 hours, and 35 genes downregulated at 21 hours.
Temporal response of the nifHDK operon and selected heterocyst glycolipid (hgl) and polysaccharide (hep) genes to nitrogen step-down
Formation of the heterocyst envelope involves deposition of a polysaccharide outer layer followed by deposition of an inner glycolipid layer [9, 11]. The RNA-seq data show that the genes responsible for heterocyst exopolysaccharide synthesis (hep genes) were upregulated by 12 hours after nitrogen step-down (Table 1). However, strong upregulation of the genes required for heterocyst glycolipid synthesis (hgl genes) did not occur until the 21-hour sample. These data show that during heterocyst morphogenesis, the polysaccharide genes are expressed first, likely depositing the stabilizing exopolysaccharide later; subsequently, hgl genes are expressed to produce the underlying glycolipid envelope layer, which together are required to help create a microoxic environment within the heterocyst.
Identification of unstudied genes regulated in response to nitrogen deprivation
Unstudied regulatory genes with 5-fold or greater increase in expression by 6, 12, or 21 hours after nitrogen step down1
RPKM fold change
GTP-binding protein TypA/BipA
probable GTP-binding protein
two-component sensor histidine kinase
two-component sensor histidine kinase
two-component hybrid sensor and regulator
two-component response regulator
Transposase gene families with a 2-fold or greater increase in expression by 6 hours after nitrogen step-down1
Average Fold Change in RPKM
0 to 6 h
0 to 12 h
0 to 21 h
However, for other genes such as atp1 (ATP synthase) and rbcL (carbon fixation), the pattern of reads showed evidence for transcripts extending upstream from published 5' ends (Figure 2). For both atp1 and rbcL, there is a considerable drop in read coverage at the previously identified 5' end, but there is still significant coverage upstream of these sites. This suggests that the previously identified 5' end may be a processing site or a site of RNA secondary structure that affects primer extension results on full-length transcripts, and that some of these transcripts originate at upstream start sites. Our analysis of transcript 5' ends suggests that transcription initiation and transcript processing in Anabaena is often complex, and that many transcripts have long 5' UTRs with unclear transcriptional start sites or processing sites. These "trailing" 5' ends make identification of transcriptional start sites via primer extension or RACE difficult for some genes, and RNA-seq will be the method of choice for mapping transcripts.
The RNA-seq data can be used to map transcripts and examine the coordinated expression of genes that form operons. For example, the genes in the nifB-fdxN-nifS-nifU operon showed very low numbers of reads at 0, 6, and 12 hours during nitrogen step-down, when proheterocysts have yet to create a microoxic environment for nitrogen fixation (Figure 3). By 21 hours, there was a large increase in the number of reads for all four genes, indicating that a single transcript for all four genes originates from a promoter upstream of nifB. The RNA-seq data showed a clear 5' end at position -282 upstream of the nifB start codon, which validates previously published results . Further research will be required to fully describe the promoters and operon structure for the entire set of nitrogen fixation genes, which may be unexpectedly complex . Additional RNA-seq data or other types of approaches will be required to clarify transcription in certain regions because of fluctuations in read coverage present in our data set. These fluctuations in reads across ORFs and a general increase in coverage around the 5' ends of ORFs are consistent with other RNA-seq datasets [28, 37–39]. It is likely that RNA stability and secondary structure contribute to these fluctuations in coverage.
Identification of antisense RNAs
The small RNA prep protocol maintains information on the direction of transcription by adding different adaptors to the 3' and 5' ends of each RNA molecule in the sample prior to cDNA synthesis. Therefore, we were able to identify antisense RNAs in ORFs or 5' UTRs of annotated genes and also identify transcripts from unannotated regions of the genome (data available at GEO accession #GSE26633). The antisense transcripts would not have been identified with standard microarray or RNA deep sequencing methods because these methods do not normally distinguish between sense and antisense transcripts.
Antisense RNAs transcribed within the ORF or 5' UTR of genes involved in heterocyst differentiation
Gene Expression After Nitrogen Step-down
Antisense RNA Expression After Nitrogen Step-down
increased by 12 h
narB ORF and 5' UTR
increased by 12 h
slight decrease by 21 h
heterocyst specific ABC transporter
increased by 12 h
5' end of alr3649 ORF into upstream gene
increased by 12 h
similar to nitrogen regulation protein NtrR
decreased slightly at 6 h only
increased by 21 h
glnB 5' UTR
heterocyst differentiation regulator
increased by 6 h
hetR ORF and 5' UTR
decreased by 21 h
increased by 12 h
two RNAs, one in hetC ORF and one in hetC-hetP intergenic region
increased by 12 h
hetF ORF and 5' UTR; NsiR1*
increased by 6 h
nitrogen assimilation regulation
increased by 21 h
increased by 12 h
hepK ORF; contains repeat sequence at 5' end
increased by 12 h; largest increase at 21 h
nblA ORF and 5' and 3' UTRs
antisense transcription until 12 h only
increased by 12 h; largest increase at 21 h
hglE 3' end of ORF
increased by 12 h
increased by 21 h
hglD ORF and 5' UTR
heterocyst specific glycolipid gene
increased by 6 h; largest increase at 21 h
increased at 6 h only
Some genes showed striking changes in antisense RNAs. For example, nblA, which is involved in degradation of phycobilisome proteins in response to nitrogen deprivation , showed extensive antisense reads at 0 hours (Figure 4). After nitrogen step-down, the antisense reads decreased as nblA sense reads increased (Figure 4). We hypothesize that it is critical for vegetative cells to avoid expression of even small amounts of NblA protein and that the antisense RNA ensures no expression from the nblA gene. The ratio of antisense to sense nblA RNA decreases after nitrogen step-down, likely because heterocysts begin to express the coding transcript. Future studies will be required to determine if sense and antisense RNAs for nblA are differentially expressed in heterocysts and vegetative cells.
Antisense transcription within genes involved in heterocyst development (Table 4) suggests that antisense RNAs may be an important mechanism of regulation during the response to combined nitrogen deprivation. Future analysis of these noncoding RNAs will shed light on the regulation of genes required for heterocyst development and diazotrophic growth.
Overall, our data confirm directional RNA deep sequencing as a more thorough method for analyzing transcriptional regulation in cyanobacteria and indicate that further studies using different environmental conditions and mutant strains will yield novel information about Anabaena gene regulation. Our RNA-seq data can be used to improve gene annotation and map RNA ends and operon structure. In addition, directional RNA-seq data provide superior information compared to 5' RACE and primer extension experiments for mapping RNA transcripts and identifying potential promoter regions. Because we used a strand-specific sequencing protocol, the dataset can be used to identify antisense and other noncoding RNAs potentially involved in gene regulation. Finally, our work has provided a systems-level view of the Anabaena transcriptome during the response to nitrogen step-down. Together, these features of the directional RNA-seq data can be used to define future directions for studying heterocyst development.
Preparation of RNA samples
For deep sequencing, total RNA was prepared from Anabaena (Nostoc) sp. strain PCC 7120 cultures grown in 100 ml of liquid medium in 250-ml flasks with cotton plugs as previously described with slight modifications . Briefly, 100-ml liquid cultures were grown to an OD750 of 0.5 in BG-11(NH4) medium, which lacked sodium nitrate and contained 2.5 mM ammonium chloride and 5 mM MOPS (pH 8.0). For the 0 hour sample, cells were collected before deprivation of combined nitrogen. For nitrogen-deprived samples, cells were collected by centrifugation and washed 3 times in BG-110, which lacks sodium nitrate, resuspended in BG-110 to an OD750 of approximately 0.05, and incubated for 6, 12, or 21 hours. These times were chosen because at 6 hours cells are beginning to establish a pattern of differentiating cells, at 12 hours proheterocysts are committed to becoming heterocysts , and at 21 hours nitrogen-fixing heterocysts are fully differentiated. For each sample, cells were rapidly cooled and collected by pouring 50 ml of culture over 100 g crushed ice and centrifugation at 4,000 × g for 10 minutes at 4°C. Cell pellets were transferred to a 2-ml tube and collected by centrifugation at 11,500 × g for 2 minutes at 4°C. Supernatant was removed and cell pellets were flash frozen in liquid nitrogen for storage at -80°C. RNA was isolated from filaments with the Ambion RiboPure RNA isolation kit according to the manufacturer's protocol and total RNA was used for deep sequencing.
Preparation of a cDNA library
Ribosomal RNA was not removed from total RNA to avoid any depletion of coding transcripts. The presence of ribosomal and transfer RNAs in our sequencing sample did result in a decreased yield of transcript information from coding RNAs (only ~10% of the sample at each time point was non-rRNA and non-tRNA transcripts); however, this protocol avoided the known biases introduced by current rRNA depletion methods . Furthermore, samples were not multiplexed, which allows the sequencing of multiple samples in the same lane, to avoid biases that result from using different adapters for different samples. A directional RNA sequencing library was prepared from 1 μg of total RNA from each sample (0, 6, 12, and 21 h). The RNA samples were purified using a Qiagen RNeasy MinElute Cleanup Kit, fragmented for 180 seconds using a Covaris S2 sonicator set at 10% duty cycle, 5 Intensity, 200 cycles/burst in 120 μl of 1 mM Tris-EDTA pH 8.0 and purified again with the cleanup kit. The fragmented RNA (100 ng) was dephosphorylated using Antarctic phosphatase (2 units, 37°C, 30 min), 5' phosphorylated using T4 polynucleotide kinase (20 units, 37°C, 60 min), and purified with the cleanup kit. cDNA library preparation steps were performed as described in the Illumina Small RNA v1.5 Sample Preparation Guide except that during size selection on 4% agarose gels fragments of approximately 250 bp were obtained, which are suitable for up to 100 base sequencing reads. We saw no evidence of genomic DNA contamination of the RNA samples prior to cDNA synthesis, as there were many regions with no read coverage across the genome.
Illumina sequencing and analysis
For sequencing, the library was denatured and diluted following standard Illumina-recommended protocols to a final concentration of 9 pM before being loaded onto an Illumina single read flow-cell for massively parallel sequencing on an Illumina GAIIx. Raw sequences were obtained from GA Pipeline software using CASAVA v1.7. The reads were further processed to remove any adapter sequence using the Illumina Flicker add-on. Flicker v2.7 trims the adaptor sequence from each read and does iterative alignment to the reference genome using ELAND.
One sample was sequenced per lane, yielding an average of 17 million high quality 36-bp reads per sample for the 0-, 12-, and 21-hour time points, 1.7 million of which were from non-rRNA and non-tRNA transcripts. The 6-hour sample was sequenced with longer 100-bp reads; however, these reads were trimmed to 40 bp for our analysis. RNA-seq data were aligned and analyzed with CLC Genomics Workbench 4 to create SAM files that were further analyzed with the Cufflinks software suite . CLC genomics workbench 4 was used to generate expression profiles and clustering. Additional files 1, 2, 3, 4, 5, 6, and 7 contain annotated read data, RPKM values, and fold-change calculations for the time-course expression data from the chromosome and six plasmids.
RNA-seq data are available through the NCBI Gene Expression Omnibus (GEO) database, accession number #GSE26633. The raw sequence reads as well BAM files of reads at each time point aligned to NCBI's current build of the Anabaena sp. strain PCC 7120 genome are included in the accession. Raw reads are .txt files and can be opened with a FASTA viewer. Aligned reads are in .BAM format and can be analyzed with free or commercial software suites.
Acknowledgements and Funding
We thank Rodrigo Mella-Herrera for help with RNA purification, Arnaud Taton for bioinformatics help, and Eric Allen for help with experimental design, software, and use of his facilities. This work was supported by National Science Foundation Grant 0925126 (JWG). BLF was supported by National Institutes of Health Training Grant T32GM007240. FVN was supported by the Fund for Scientific Research-Flanders, Belgium (F.W.O.-Vlaanderen).
- Herrero A, Flores E, (Eds): Cyanobacteria: Molecular Biology, Genomics and Evolution. 2008, Caister Academic PressGoogle Scholar
- Dismukes GC, Carrieri D, Bennette N, Ananyev GM, Posewitz MC: Aquatic phototrophs: efficient alternatives to land-based crops for biofuels. Curr Opin Biotechnol. 2008, 19: 235-240. 10.1016/j.copbio.2008.05.007.View ArticlePubMedGoogle Scholar
- Ducat DC, Way JC, Silver PA: Engineering cyanobacteria to generate high-value products. Trends Biotechnol. 2011, 29: 95-103. 10.1016/j.tibtech.2010.12.003.View ArticlePubMedGoogle Scholar
- Bothe H, Schmitz O, Yates MG, Newton WE: Nitrogen fixation and hydrogen metabolism in cyanobacteria. Microbiol Mol Biol Rev. 2010, 74: 529-551. 10.1128/MMBR.00033-10.View ArticlePubMedPubMed CentralGoogle Scholar
- Flores E, Herrero A: Compartmentalized function through cell differentiation in filamentous cyanobacteria. Nat Rev Micro. 2010, 8: 39-50. 10.1038/nrmicro2242.View ArticleGoogle Scholar
- Kumar K, Mella-Herrera RA, Golden JW: Cyanobacterial heterocysts. Cold Spring Harb Perspect Biol 2009. 2010, 2: a000315-10.1101/cshperspect.a000315.Google Scholar
- Vargas WA, Nishi CN, Giarrocco LE, Salerno GL: Differential roles of alkaline/neutral invertases in Nostoc sp. PCC 7120: Inv-B isoform is essential for diazotrophic growth. Planta. 2010, 233: 153-162.View ArticlePubMedGoogle Scholar
- Lopez-Igual R, Flores E, Herrero A: Inactivation of a heterocyst-specific invertase indicates a principal role of sucrose catabolism in heterocysts of Anabaena sp. J Bacteriol. 2010, 192: 5526-5533. 10.1128/JB.00776-10.View ArticlePubMedPubMed CentralGoogle Scholar
- Awai K, Lechno-Yossef S, Wolk CP: Heterocyst envelope glycolipids. Lipids in Photosynthesis: Essential and Regulatory Functions. Edited by: Wada H, Murata N. 2009, New York: Springer, 179-202.View ArticleGoogle Scholar
- Nicolaisen K, Hahn A, Schleiff E: The cell wall in heterocyst formation by Anabaena sp. PCC 7120. J Basic Microbiol. 2009, 49: 5-24. 10.1002/jobm.200800300.View ArticlePubMedGoogle Scholar
- Huang G, Fan Q, Lechno-Yossef S, Wojciuch E, Wolk CP, Kaneko T, Tabata S: Clustered genes required for the synthesis of heterocyst envelope polysaccharide in Anabaena sp. strain PCC 7120. J Bacteriol. 2005, 187: 1114-1123. 10.1128/JB.187.3.1114-1123.2005.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao MX, Jiang YL, He YX, Chen YF, Teng YB, Chen Y, Zhang CC, Zhou CZ: Structural basis for the allosteric control of the global transcription factor NtcA by the nitrogen starvation signal 2-oxoglutarate. Proc Natl Acad Sci USA. 2010, 107: 12487-12492. 10.1073/pnas.1001556107.View ArticlePubMedPubMed CentralGoogle Scholar
- Ehira S, Ohmori M: NrrA, a nitrogen responsive response regulator facilitates heterocyst development in the cyanobacterium Anabaena sp. strain PCC 7120. Mol Microbiol. 2006, 59: 1692-1703. 10.1111/j.1365-2958.2006.05049.x.View ArticlePubMedGoogle Scholar
- Ehira S, Ohmori M: NrrA directly regulates expression of hetR during heterocyst differentiation in the cyanobacterium Anabaena sp. strain PCC 7120. J Bacteriol. 2006, 188: 8520-10.1128/JB.01314-06.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang CC, Laurent S, Sakr S, Peng L, Bedu S: Heterocyst differentiation and pattern formation in cyanobacteria: a chorus of signals. Mol Microbiol. 2006, 59: 367-375. 10.1111/j.1365-2958.2005.04979.x.View ArticlePubMedGoogle Scholar
- Ehira S, Ohmori M, Sato N: Genome-wide expression analysis of the responses to nitrogen deprivation in the heterocyst-forming cyanobacterium Anabaena sp. strain PCC 7120. DNA Res. 2003, 10: 97-113. 10.1093/dnares/10.3.97.View ArticlePubMedGoogle Scholar
- Ernst A, Black TA, Cai Y, Panoff JM, Tiwari DN, Wolk CP: Synthesis of nitrogenase in mutants of the cyanobacterium Anabaena sp. strain PCC 7120 affected in heterocyst development or metabolism. J Bacteriol. 1992, 174: 6025-PubMedPubMed CentralGoogle Scholar
- Nayar AS, Yamaura H, Rajagopalan R, Risser DD, Callahan SM: FraG is necessary for filament integrity and heterocyst maturation in the cyanobacterium Anabaena sp. strain PCC 7120. Microbiology. 2007, 153: 601-10.1099/mic.0.2006/002535-0.View ArticlePubMedGoogle Scholar
- Liu D, Golden JW: hetL overexpression stimulates heterocyst formation in Anabaena sp. strain PCC 7120. J Bacteriol. 2002, 184: 6873-10.1128/JB.184.24.6873-6881.2002.View ArticlePubMedPubMed CentralGoogle Scholar
- Higa KC, Callahan SM: Ectopic expression of hetP can partially bypass the need for hetR in heterocyst differentiation by Anabaena sp. strain PCC 7120. Mol Microbiol. 2010, 77: 562-574. 10.1111/j.1365-2958.2010.07257.x.View ArticlePubMedGoogle Scholar
- Mella-Herrera RA, Neunuebel MR, Golden JW: Anabaena sp. strain PCC 7120 conR contains a LytR-CpsA-Psr domain, is developmentally regulated, and is essential for diazotrophic growth and heterocyst morphogenesis. Microbiology. 2010, 157: 617-626.View ArticlePubMedGoogle Scholar
- Campbell EL, Summers ML, Christman H, Martin ME, Meeks JC: Global gene expression patterns of Nostoc punctiforme in steady state dinitrogen-grown heterocyst-containing cultures, and at single time points during the differentiation of akinetes and hormogonia. J Bacteriol. 2007, 189: 5247-5256. 10.1128/JB.00360-07.View ArticlePubMedPubMed CentralGoogle Scholar
- Filiatrault MJ, Stodghill PV, Bronstein PA, Moll S, Lindeberg M, Grills G, Schweitzer P, Wang W, Schroth GP, Luo S: Transcriptome analysis of Pseudomonas syringae identifies new genes, noncoding RNAs, and antisense activity. J Bacteriol. 2010, 192: 2359-2372. 10.1128/JB.01445-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Otto TD, Wilinski D, Assefa S, Keane TM, Sarry LR, Böhme U, Lemieux J, Barrell B, Pain A, Berriman M: New insights into the blood-stage transcriptome of Plasmodium falciparum using RNA-Seq. Mol Microbiol. 2010, 76: 12-24. 10.1111/j.1365-2958.2009.07026.x.View ArticlePubMedPubMed CentralGoogle Scholar
- Passalacqua KD, Varadarajan A, Ondov BD, Okou DT, Zwick ME, Bergman NH: Structure and complexity of a bacterial transcriptome. J Bacteriol. 2009, 191: 3203-3211. 10.1128/JB.00122-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Sorek R, Cossart P: Prokaryotic transcriptomics: a new view on regulation, physiology and pathogenicity. Nat Rev Genet. 2009, 11: 9-16.View 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: 511-515. 10.1038/nbt.1621.View ArticleGoogle Scholar
- Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Meth. 2008, 5: 621-628. 10.1038/nmeth.1226.View ArticleGoogle Scholar
- Wolk CP, Lechno-Yossef S, Jager KM: The insertion sequences of Anabaena sp. strain PCC 7120 and their effects on its open reading frames. J Bacteriol. 2010, 192: 5289-5303. 10.1128/JB.00460-10.View ArticlePubMedPubMed CentralGoogle Scholar
- Mulligan ME, Haselkorn R: Nitrogen fixation (nif) genes of the cyanobacterium Anabaena species strain PCC 7120. The nifB-fdxN-nifS-nifU operon. J Biol Chem. 1989, 264: 19200-19207.PubMedGoogle Scholar
- Valladares A, Muro-Pastor AM, Herrero A, Flores E: The NtcA-dependent P1 promoter is utilized for glnA expression in N2-fixing heterocysts of Anabaena sp. strain PCC 7120. J Bacteriol. 2004, 186: 7337-7343. 10.1128/JB.186.21.7337-7343.2004.View ArticlePubMedPubMed CentralGoogle Scholar
- Vrba JM, Curtis SE: Characterization of a four-member psbA gene family from the cyanobacterium Anabaena PCC 7120. Plant Mol Biol. 1990, 14: 81-92. 10.1007/BF00015657.View ArticlePubMedGoogle Scholar
- Alam J, Whitaker R, Krogmann D, Curtis S: Isolation and sequence of the gene for ferredoxin I from the cyanobacterium Anabaena sp. strain PCC 7120. J Bacteriol. 1986, 168: 1265-PubMedPubMed CentralGoogle Scholar
- Lang JD, Haselkorn R: Isolation, sequence and transcription of the gene encoding the photosystem II chlorophyll-binding protein, CP-47, in the cyanobacterium Anabaena 7120. Plant Mol Bio. 1989, 13: 441-456. 10.1007/BF00015556.View ArticleGoogle Scholar
- Muro-Pastor AM, Olmedo-Verd E, Flores E: All4312, an NtcA-regulated two-component response regulator in Anabaena sp. strain PCC 7120. FEMS Microbiology Letters. 2006, 256: 171-177. 10.1111/j.1574-6968.2006.00136.x.View ArticlePubMedGoogle Scholar
- Ungerer JL, Pratte BS, Thiel T: RNA processing of nitrogenase transcripts in the cyanobacterium Anabaena variabilis. J Bacteriol. 2010, 192: 3311-3320. 10.1128/JB.00278-10.View ArticlePubMedPubMed CentralGoogle Scholar
- Jacob NM, Kantardjieff A, Yusufi FNK, Retzel EF, Mulukutla BC, Chuah SH, Yap M, Hu WS: Reaching the depth of the Chinese hamster ovary cell transcriptome. Biotechnol Bioeng. 2009, 105: 1002-1009.Google Scholar
- Mizrachi E, Hefer C, Ranik M, Joubert F, Myburg A: De novo assembled expressed gene catalog of a fast-growing Eucalyptus tree produced by Illumina mRNA-Seq. BMC Genomics. 2010, 11: 681-693. 10.1186/1471-2164-11-681.View ArticlePubMedPubMed CentralGoogle Scholar
- Ozsolak F, Milos PM: RNA sequencing: advances, challenges and opportunities. Nat Rev Genet. 2011, 12: 87-98. 10.1038/nrg2934.View ArticlePubMedGoogle Scholar
- Ionescu D, Voss B, Oren A, Hess WR, Muro-Pastor AM: Heterocyst-specific transcription of NsiR1, a non-coding RNA encoded in a tandem array of direct repeats in cyanobacteria. J Mol Biol. 2010, 398: 177-188. 10.1016/j.jmb.2010.03.010.View ArticlePubMedGoogle Scholar
- Buikema WJ, Haselkorn R: Characterization of a gene controlling heterocyst differentiation in the cyanobacterium Anabaena 7120. Genes Dev. 1991, 5: 321-330. 10.1101/gad.5.2.321.View ArticlePubMedGoogle Scholar
- Khudyakov I, Wolk CP: hetC, a gene coding for a protein similar to bacterial ABC protein exporters, is involved in early regulation of heterocyst differentiation in Anabaena sp. strain PCC 7120. J Bacteriol. 1997, 179: 6971-6978.PubMedPubMed CentralGoogle Scholar
- Muro-Pastor AM, Valladares A, Flores E, Herrero A: The hetC gene is a direct target of the NtcA transcriptional regulator in cyanobacterial heterocyst development. J Bacteriol. 1999, 181: 6664-6669.PubMedPubMed CentralGoogle Scholar
- Baier K, Lehmann H, Stephan DP, Lockau W: NblA is essential for phycobilisome degradation in Anabaena sp. strain PCC 7120 but not for development of functional heterocysts. Microbiology. 2004, 150: 2739-2749. 10.1099/mic.0.27153-0.View ArticlePubMedGoogle Scholar
- Golden JW, Whorff LL, Wiest DR: Independent regulation of nifHDK operon transcription and DNA rearrangement during heterocyst differentiation in the cyanobacterium Anabaena sp. strain PCC 7120. J Bacteriol. 1991, 173: 7098-7105.PubMedPubMed CentralGoogle Scholar
- Yoon HS, Golden JW: PatS and products of nitrogen fixation control heterocyst pattern. J Bacteriol. 2001, 183: 2605-2613. 10.1128/JB.183.8.2605-2613.2001.View ArticlePubMedPubMed CentralGoogle Scholar
- He S, Wurtzel O, Singh K, Froula JL, Yilmaz S, Tringe SG, Wang Z, Chen F, Lindquist EA, Sorek R, Hugenholtz P: Validation of two ribosomal RNA removal methods for microbial metatranscriptomics. Nat Meth. 2010, 7: 807-812. 10.1038/nmeth.1507.View ArticleGoogle Scholar
- Reich M, Liefeld T, Gould J, Lerner J, Tamayo P, Mesirov JP: GenePattern 2.0. Nat Genetics. 2006, 38: 500-501. 10.1038/ng0506-500.View ArticlePubMedGoogle Scholar
- Xu X, Elhai J, Wolk CP: Transcriptional and developmental responses by Anabaena to deprivation of fixed nitrogen. The cyanobacteria: molecular biology, genomics and evolution. Edited by: Herrero A, Flores E. Norwich. 2008, United Kingdom: Horizon Scientific Press, 383-422.Google 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.