This Article Abstract Full Text (PDF) Supplemental material Other Versions of this Article: JB.00360-07v1 189/14/5247    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Reprints and Permissions Copyright Information Books from ASM Press MicrobeWorld Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Campbell, E. L. Articles by Meeks, J. C. Search for Related Content PubMed PubMed Citation Articles by Campbell, E. L. Articles by Meeks, J. C.

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 ^,

Section of Microbiology, University of California, Davis, California 95616,¹ Department of Biology, California State University, Northridge, California 91330²

Received 9 March 2007/ Accepted 26 April 2007

	ABSTRACT

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

 
The vegetative cells of the filamentous cyanobacterium Nostoc punctiforme can differentiate into three mutually exclusive cell types: nitrogen-fixing heterocysts, spore-like akinetes, and motile hormogomium filaments. A DNA microarray consisting of 6,893 N. punctiforme genes was used to identify the global transcription patterns at single time points in the three developmental states, compared to those in ammonium-grown time zero cultures. Analysis of ammonium-grown cultures yielded a transcriptome of 2,935 genes, which is nearly twice the size of a soluble proteome. The NH₄⁺-grown transcriptome was enriched in genes encoding core metabolic functions. A steady-state N₂-grown (heterocyst-containing) culture showed differential transcription of 495 genes, 373 of which were up-regulated. The majority of the up-regulated genes were predicted from studies of heterocyst differentiation and N₂ fixation; other genes are candidates for more detailed genetic analysis. Three days into the developmental process, akinetes showed a similar number of differentially expressed genes (497 genes), which were equally up- and down-regulated. The down-regulated genes were enriched in core metabolic functions, consistent with entry into a nongrowth state. There were relatively few adaptive genes up-regulated in 3-day akinetes, and there was little overlap with putative heterocyst developmental genes. There were 1,827 differentially transcribed genes in 24-h hormogonia, which was nearly fivefold greater than the number in akinete-forming or N₂-fixing cultures. The majority of the up-regulated adaptive genes were genes encoding proteins for signal transduction and transcriptional regulation, which is characteristic of a motile filament that is poised to sense and respond to the environment. The greatest fraction of the 883 down-regulated genes was involved in core metabolism, also consistent with entry into a nongrowth state. The differentiation of heterocysts (steady state, N₂ grown), akinetes, and hormogonia appears to involve the up-regulation of genes distinct for each state.

	INTRODUCTION

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

 
The vegetative cells of the filamentous, oxygenic, phototrophic cyanobacterium Nostoc punctiforme have the potential for four developmental alternatives (32). In the presence of replete nutrients, including a source of combined nitrogen such as ammonium, light, and CO₂ or in the dark with glucose or fructose as a source of reduced carbon, a permissive vegetative cell cycle consisting of growth and division is maintained. Nutrient depletion or change elicits three mutually exclusive differentiation responses. The differentiated cells are morphologically and physiologically distinct from vegetative cells. The question to be asked is, what is the extent of differential gene expression in the developmental states of N. punctiforme?

Depletion of combined nitrogen induces the transcription of genes encoding the nitrogenase enzyme complex. Nitrogenase expression is preceded by the differentiation of heterocysts, which are specialized microoxic cells in which oxygen-sensitive nitrogen fixation takes place. The microoxic environment of heterocysts is achieved by inactivation of the photosystem II oxygen-evolving reaction, an increased rate of aerobic respiration, and the synthesis of a bilayer envelope that retards diffusion of gases (44). Heterocysts are terminally differentiated; they constitute 10% or less of the cells in a filament and are singly spaced in a nonrandom pattern (33). They are the source of organic nitrogen for vegetative cells and a sink for photosynthate from vegetative cells. Heterocyst differentiation allows continued vegetative cell growth in a habitat depleted of combined nitrogen. The formation and maintenance of heterocysts are variously estimated to involve the differential expression of 140 (43) to 600 (12) to more than 1,000 (28) genes, and this is the most highly studied cyanobacterial differentiation event (for reviews, see references 1, 17, 22, 33, 43, 44, and 47).

When limited in cellular energy, typically in response to phosphate deprivation, vegetative cells differentiate into spore-like cells termed akinetes. Although akinetes retain some metabolic activity (40), they do not grow. Akinetes are resistant to cold and desiccation, thereby allowing perennation or longer periods of survival. This is a transient differentiation process; when environmental conditions are appropriate for growth, akinetes germinate, and the emerging filaments continue a vegetative cell cycle. Akinete differentiation is confined to the heterocyst-forming cyanobacteria (36). Based on this taxonomic clustering and some biochemical and genetic similarities, there is speculation that akinetes served as the evolutionary progenitors of heterocysts (44). Although a priori there is no obvious reason that heterocyst differentiation and akinete differentiation should be mutually exclusive, continued heterocyst differentiation is generally not observed when environmental conditions favor akinete differentiation. In heterocyst-containing filaments of N. punctiforme, akinetes first appear midway in the interval between two heterocysts, with subsequent differentiation events radiating in both directions away from the first akinete and towards the adjacent heterocysts (32). However, in a glucose-6-phosphate dehydrogenase (zwf) mutant of N. punctiforme that is unable to fix nitrogen or grow in the dark, vegetative cells synchronously differentiate into akinetes when they are placed in the dark (4). Five genes that are differentially transcribed during akinete differentiation have recently been identified (6, 48).

Vegetative filaments of N. punctiforme are nonmotile. Gliding and buoyant mobility can be induced following the transient differentiation of vegetative filaments into small-celled filaments called hormogonia. No single unique environmental factor induces hormogonium differentiation; rather, changes in a multiplicity of factors, including depleted and replete nutrients, light quality, the presence of symbiotic plant partners, or dilution into fresh medium, stimulate hormogonium differentiation to various degrees (33, 39). Once induced, N. punctiforme filaments proceed through a non-N₂-fixing hormogonium differentiation cycle that lasts for 72 to 96 h (9). Upon entry into the cycle, there is a cessation of net increases in biomass, DNA, protein, and pigments (9, 10, 20), while cell division continues, resulting in an increase in the number of cells, all of which are smaller and shaped differently than the parental vegetative cells (9, 10, 32, 36, 39). Since filamentous cyanobacteria contain multiple copies of the chromosome (19), there is a high probability that hormogonium cells contain at least one genome copy. In N. punctiforme, gliding motility is maximal between 18 and 72 h after induction (9). Once the filaments become sessile, biomass increase resumes and the filaments exit the hormogonium cycle to reenter the vegetative cell cycle in the presence or absence of combined nitrogen. There have been few molecular studies of events in the hormogonium cycle. It is known that the transcription of genes encoding gas vesicle proteins, pili, (39), an alternative sigma subunit of RNA polymerase, and a Clp protease (8) is up-regulated during entry into the cycle, while synthesis of phycobiliproteins is down-regulated (39).

The complete 9.059-Mbp genome sequence of N. punctiforme strain ATCC 29133 has been determined. A preliminary report has been published (34), and the sequence can be accessed at the Integrated Microbial Genomes site (http://img.jgi.doe.gov/cgi-bin/pub/main.cgi). The complete genome sequence allows entry into systems-level approaches to identification of gene expression profiles. We report here the construction of a DNA microarray of the N. punctiforme genome. This array was used to identify genes that are transcribed differentially during heterocyst, akinete, and hormogonium differentiation or function by reference to transcription during growth of vegetative cells on NH₄⁺.

	MATERIALS AND METHODS

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

 
Cultures and culture conditions. N. punctiforme strain ATCC 29133 (= PCC 73102) and its derivatives were grown with shaking at 100 to 120 rpm and at 23 to 25°C under ca. 19 to 46 µmol photons/m²/s of cool white fluorescent light in the standard minimal salts medium of Allan and Arnon, diluted fourfold, as described previously (14). Since heterocyst differentiation results in perpetuation of a modified vegetative cell cycle, comparisons were made between steady-state NH₄⁺- and N₂-grown cultures. A spontaneously derived strain of N. punctiforme, UCD 102, which does not form hormogonia, was used for these and subsequent N₂ growth experiments so that heterocyst differentiation was not complicated by formation of hormogonia. Duplicate 50-ml cultures from the same N₂-grown inoculum (200 µl) were incubated for 7 days under the same environmental conditions, with one culture supplemented with 2.5 mM NH₄Cl and 5.0 mM morpholinepropanesulfonic acid (MOPS) and the other culture without these compounds (N₂ steady state). Cultures grown with NH₄⁺ were devoid of heterocysts, as determined by microscopy, and were used as the reference. Under these low-cell-density culture conditions, the heterocyst frequency of the N₂ steady-state culture was, on average, 8.7%. The cultures were harvested at densities of 2 to 3 µg chlorophyll a per ml by centrifugation at 1,000 x g for 5 min and stored at –80°C for subsequent RNA extraction.

The gliding motility of N. punctiforme hormogonia is pronounced 24 h after induction, and this time point was chosen for these experiments, using the symbiotically competent strain UCD 153. Cultures were grown in 500 ml of minimal medium supplemented with 2.5 mM NH₄Cl and 5.0 mM MOPS for 7 days, reaching densities of 2 to 3 µg chlorophyll a/ml; the cultures were then synchronized for differentiation of hormogonia by centrifugation as described above and suspended in fresh medium without supplementation with combined nitrogen. Time zero (T₀) samples were collected at this time. After 24 h under standard incubation conditions, cultures were examined by microscopy for verification of hormogonium formation, and samples showing essentially 100% differentiation were collected by centrifugation and stored as described above.

Akinete differentiation occurs over a period of up to 6 days (4). A gene with enhanced expression in akinetes, avaK, shows strong transcription 3 days following induction (4), as well as does the onset of resistance to desiccation, cold, or lysozyme treatments (5). Therefore, cells in the process of differentiation into akinetes were induced by incubation in the dark for 3 days of a 50-ml culture of the zwf::P_psbA-npt insertion mutant CSUN 466 (a reconstruction of UCD 466) (18) previously grown for 5 days in minimal medium supplemented with 2.5 mM NH₄Cl and 5.0 mM MOPS plus 50 mM fructose. Under these conditions, strain CSUN 466 shows synchronized differentiation into akinetes, and approximately 50% of the cells are resistant to cold and desiccation at day 3 (4, 5). Cells were collected at T₀ and at 3 days by centrifugation and stored as described above.

RNA isolation, cDNA synthesis, and dye labeling. Frozen samples were subjected to breakage using a mini bead beater (Biospec) in the presence of phenol, 0.5-mm zirconium/silicate beads, and Celite as described previously (38). RNA in the aqueous phase was precipitated using 4 M lithium chloride, 20 mM Tris (pH 7.4), 10 mM EDTA (pH 8). The samples were then further purified using QIAGEN RNeasy columns by following the manufacturer's specifications, including the on-column DNase digestion to eliminate DNA contamination. mRNA was enriched from 10 µg of total RNA using a cyanobacterial capture module in a MicrobExpress rRNA subtraction kit (Ambion); the average yield was 1.05 µg of RNA. Fifteen micrograms of total RNA was used in each cDNA synthesis reaction in the presence of 4 µg of a random 18-mer, 10 µg of random hexamers, and a deoxynucleoside triphosphate mixture containing aminoallyl-modified nucleotides or 1 µg of purified mRNA and 10 µg of random hexamers; otherwise, the manufacturer's instructions (Invitrogen) were used. Alexa 647 and Alexa 555 dyes (Invitrogen) were coupled to samples prior to hybridizations. The fluorescently labeled probes were purified using QIAGEN PCR "Qiaquick" columns and were quantified according to the manufacturer's instructions.

Array construction and hybridization. The N. punctiforme array consists of internal gene fragments ranging in size from 53 to 725 bp, with the great majority between 200 and 700 bp. Unique primer pairs for each gene in the array were designed and synthesized under a contract with Operon. PCRs were performed under standard conditions using genomic DNA as the template. The products were then purified using QIAGEN columns in a 96-well format (Qiaquick 96), and the concentrations were adjusted to between 0.1 and 0.3 µM. There are 6,893 individual gene fragments in the array; in addition, five genes are spotted in duplicate as controls (rbcS [NpF4197], rbcL [NpF4195], rpoC [NpF4987], glnA [NpR5387], and hemL [NpR4096]), four concentrations of sheared genomic DNA, and salmon sperm DNA as a negative control. Absent in the array are 351 of the 416 transposase genes with multiple members in various families, three psbA and one psbD family member genes, and 75 hypothetical genes less than 200 bp long. DNA fragments were spotted in duplicate onto aminosilane A+-coated slides (Schott Nexterion) by the ArrayCore Facility at UC Davis (array.ucdavis.edu). All hybridizations were done with an automated Tecan HS4800 hybridization station at 50°C. Slides were scanned using a GenePix 4000A scanner (Axon Instruments).

Data analysis. Raw data in a .gpr file format were generated using the GenePix Pro 3.0 software. Three biological replicates were used in the analysis for each of the differentiation states, and each of the replicates was analyzed with dye swaps. Thus, data from six hybridizations, with T₀ as the reference, for each experiment (three biological replicates plus the reciprocal dye swaps) for a total of 12 technical replicates (since each array was printed in duplicate on the slide) were used in each analysis. To determine which genes showed statistically significant differential expression in the three developmental states relative to T₀ expression, the free source software package LIMMA GUI (Linear model for microarray analysis graphical user interface) (37), written in the R-programming language (www.R-project.org) and available through Bioconductor (www.Bioconductor.org), was employed. The data were first normalized within arrays using Print-Tip Loess and across arrays in the analysis using Aquantile methods. The Linear model fit was computed using the duplicates on the arrays and the least-squares method, and all P values were adjusted using the Benjamini method (reference 37 and references therein) for controlling for the false discovery rate. Statistically significantly differentially expressed genes were defined for N. punctiforme arrays as the genes with B values of >0; this corresponded to adjusted P values of <1.49 x 10^–2 for the N₂ steady state, 4.0 x 10^–3 for hormogonia, and 1.02 x 10^–2 for akinetes. The raw data, as .gpr files, are available at http://microbiology.ucdavis.edu/faculty/jcMeeks/arrays/.

Functional gene assignments. Open reading frames were initially assigned to a functional category based on Clusters of Orthologous Groups (COG) identities from the genome computational annotation. The COG identities were manually verified or corrected by computational sequence comparisons to protein databases; the proteins with Interpro (TIGRfam, Pfam, and SwissPro) identities were accepted as accurate. Gene products were assigned to subcategories based on their physiological roles from metabolic tables and pathways (EcoCyc) or, in some instances, genetic phenotypes.

As a means to organize large data sets in a physiological framework, the NH₄⁺-grown and differentially transcribed genes were grouped into the following five functional categories based on their gene products (2): (i) core metabolism, which includes proteins required for vegetative cell growth and involved in the synthesis of precursors and monomers, ion incorporation, cofactors, polymers, assembly of cell structures, and cell division and genome segregation; (ii) transport metabolism, which involves import and export protein complexes; (iii) adaptive metabolism, including alternative C, N, P, and S macronutrient sources, signal transduction, transcriptional regulation, stress, secondary products, and developmental processes; (iv) selfish metabolism, including phage and transposition functions; and (v) unassigned metabolism, including proteins with known domains or motifs but unknown physiological functions, conserved hypothetical and hypothetical proteins, and pseudogenes (Table 1) (a detailed breakdown is provided in File S1 in the supplemental material). Some category assignments, such as proteases and DNA/RNA metabolic enzymes, are subjective and were defined based on the presence or absence in an NH₄⁺-grown soluble proteome of N. punctiforme (2) and differential expression profiles in this study. This is a dynamic database, in which category assignments can change as experimental identification of the functional role of a gene product occurs.

	RESULTS AND DISCUSSION

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

The results of a single biological and two technical replicates for NH₄⁺-grown cells and each developmental state are shown in the scatter plots in Fig. 1. When the sets of differently labeled cDNAs from NH₄⁺-grown cultures were self-hybridized, few gene spots showed more than a 1.5-fold departure from a linear function (Fig. 1A), and plots of replicate self-hybridizations showed different low-level outliers. These results verified a lack of dye bias in labeling of the cDNAs and established that in replicate experiments departures from linearity are experimentally significant and not technical artifacts.

View larger version (56K): [in this window] [in a new window]   FIG. 1. Scatter plots of a single biological and two technical replicates of two-color hybridization intensities. The y axes represent experimental (exp) intensities, and the x axes represent the reference (ref) (T0 or steady-state NH4+-grown) intensities.

There are no distinct and extended genetic islands of developmentally up- or down-regulated genes in any of the three differentiation states; the genes appear to be randomly distributed in the N. punctiforme chromosome (Fig. 2). The only appearance of island distribution in Fig. 2 is the nondifferentially transcribed blocks of genes at approximately 6:00 o'clock that putatively encode proteins for the synthesis of cyclic peptide secondary metabolites.

View larger version (75K): [in this window] [in a new window]   FIG. 2. Circular map of the chromosome. The outside circle depicts the annotated genome with color-coded COG categories of genes. The next circle shows differentially expressed genes in 24-h hormogonia, with up-regulated genes indicated by red and down-regulated genes indicated by green; this is followed by circles showing differentially expressed genes in steady-state N2-grown cultures and then 3-day akinetes.

In general, the fraction of expressed genes clustered similar to the genome categories (Table 1), although the fraction of genes encoding core metabolic functions was proportionally higher, while the fractions of genes encoding adaptive and selfish metabolism were proportionally lower. This result is reasonable since the array contains few gene fragments encoding putative transposases, and many genes encoding adaptive functions are expected to be more highly expressed under adaptive conditions. The increased proportion of expression of genes involved in core metabolism is similar to the proportion observed in the soluble proteome (2) and would be physiologically expected in a photoautotrophic NH₄⁺-grown transcriptome, as these conditions are the most permissive culture conditions for N. punctiforme. The core metabolic category contains more than 60 gene families with two to eight members, 9 of which are involved in precursor carbon metabolism. One hundred fifteen genes were expressed from the five N. punctiforme plasmids (691 total open reading frames). This set of genes is a reasonable NH₄⁺-grown transcriptome of N. punctiforme, but we do not suggest that it is definitive. It provides an internal control for genes that are differentially transcribed in cellular differentiation.

Functional categories of dinitrogen-grown cultures. The most striking result obtained with N₂-grown cultures is that the majority of the statistically significantly regulated genes were up-regulated (373 genes). Of these genes, 43% were genes involved in unassigned metabolism, 28% were genes involved in adaptive metabolism, and 24% were genes involved in core metabolism (Table 1). The unassigned metabolic genes in all differentiation states are clearly a reservoir of undescribed developmental, stress, and metabolic potential in N. punctiforme. The majority of the N₂-grown up-regulated adaptive and core genes can be predicted from the results of many studies of heterocyst differentiation and function (17, 22, 31, 33, 43, 44, 47). The presence of such known genes (e.g., het, dev, pat, hgl, hep, nif, hup, cyo, zwf, and gnd [see below]) provides an internal control to assess genome coverage.

The adaptive metabolic category includes known nitrogenase (nif), hydrogenase (hup), and nitrate assimilation structural genes and genes encoding enzymes for synthesis of the unique polysaccharide (hep) and glycolipid (hgl) heterocyst envelope. It is worth noting that steady-state expression of nitrate assimilation genes was relatively high in the absence of nitrate or nitrite. In Anabaena sp. strain PCC 7120, these genes are derepressed upon combined-nitrogen starvation but are repressed again as cultures grow with N₂ as the nitrogen source (21). Six genes encoding regulatory proteins involved in heterocyst differentiation and nitrogen assimilation are in the adaptive category; hetR (NpR1722), hetF (NpF3553), nrrA (NpR3907), and patB (NpR0334) are in the development subgroup, while ntcA (NpF5511) and ntcB (NpF1534) plus two genes encoding alternative sigma subunits of RNA polymerase (NpF0307 and NpF4153, neither of which has yet been shown to be essential for differentiation) are in the transcriptional regulation group. Genes encoding other identified heterocyst regulatory proteins, such as devH (NpR6193), hepK (NpR1012), patA (NpF5682), patN (NpF6624), and patU (NpR1831), were not present in the differentially expressed N₂ transcriptome, although devH, patA, patN, and patU are in the NH₄⁺-grown transcriptome and DevH, HepK, and PatN are in the NH₄⁺-grown proteome (2). Expression of five genes encoding serine/threonine protein kinases (PK) (NpF1910, NpF5105, NpF6098, NpF6345, and NpR5250), two PK phosphatases (NpF1519 and NpR4927), three histidine kinases (HK) (NpF0957, NpR1236, and NpR4835), and four response regulators (RR) (NpF1800, NpF5044, NpF5659, and NpR1903) was enhanced in the N₂ transcriptome, implying that these genes have roles in heterocyst differentiation or steady-state N₂ growth. These signal transduction genes are obvious targets for genetic analyses, especially analyses to identify the cognate HKs of NrrA (13) and PatA (27).

The differentially up-regulated genes in the core metabolic category include two putative operons encoding terminal cytochrome c oxidases (NpF0336 [cyoA], NpF0337 [cyoB], and NpF0338 [cyoC]; NpR3535 [cyoC], NpR3536 [cyoB], and NpR3537 [cyoA]), none of which is present in the proteome or expressed in the NH₄⁺ transcriptome, and genes encoding 15 carbon catabolism enzymes that are likely to be involved in reductant generation in heterocysts. Enhanced expression of two cytochrome oxidase complexes (41) and of oxidative pentose phosphate proteins (42) is known to occur in heterocysts. Thirteen genes encoding amino acid metabolic enzymes, including predicted heterocyst enhanced glutamine synthetase (NpR5387) (35), were up-regulated, while the level of the gene encoding vegetative cell-localized glutamate synthase (NpR3877) (29) was not different from the NH₄⁺-grown level. A glutamate dehydrogenase gene (NpF3641) was expressed in the NH₄⁺-grown transcriptome, but it was not differentially transcribed in steady-state N₂-grown cultures. The sequential activities of glutamine synthetase and glutamate synthase constitute the dominant route of NH₄⁺ assimilation in cyanobacteria, and glutamate dehydrogenase does not have a significant role (15). The remaining genes do not appear to encode enzymes specific for any complete amino acid biosynthetic pathways. The expression of five genes for synthesis of proteins for metabolism of reactive oxygen species was enhanced in the N₂-grown differential transcriptome; three of these genes (NpF2660 [ggt], NpF6212 [dps], and NpR0901 [gor], but not NpR5469 [dps] and NpR6212 [katE]), were also present in the NH₄⁺-grown proteome. The NH₄⁺-grown proteome contains another 25 proteins for metabolism of reactive oxygen (2). Whether the up-regulated gene products are localized in the heterocyst is unknown, but metabolism of reactive oxygen is clearly an important process in N. punctiforme. Nineteen genes encoding import and export transport functions are up-regulated. Their substrates are unknown, but several of the gene products are annotated as Fe and Mo ion transporters, consistent with the demands of nitrogenase. Neither of the two genes encoding an NH₄⁺ transport channel (NpF4661 and NpR3288) was significantly up- or down-regulated compared with the NH₄⁺-grown transcriptome level.

The fact that 122 genes were down-regulated in N₂-grown cultures implies that at least some of the gene products may enhance growth in the presence of NH₄⁺. A change from nitrogen-replete (NH₄⁺) to nitrogen-limited (N₂) growth is predicted to involve down-regulation of proteins for low-affinity NH₄⁺ assimilation. This was not the result observed, as exemplified by low-affinity glutamate dehydrogenase (see above). The majority of the down-regulated genes are in the unassigned (51%) and core (24%) categories, and most of the core genes encode enzymes of cell envelope and cofactor metabolism. The functionally identified core genes do not exhibit an obvious physiological pattern so that one could predict the basis of their differential expression. Genes encoding three HKs (NpF3675, NpF3677, and NpR3054), two PKs (NpF01289 and NpF1616), and two RRs (NpF3676 and NpR2902) were significantly down-regulated. The NpF3675, NpF3676, and NpF3677 genes may be transcribed as an operon, and the NpF3677 gene product is a hybrid HK with an N-terminal RR receiver domain. Since the collocated RR (NpF3676) has no DNA binding output domain, there may be additional components of this putative phosphorelay signal transduction pathway, or the regulator may modulate activity of another protein.

A recent analysis (46; see also http://www.people.vcu.edu/elhaij/microarray/) of the transcriptome data generated by Ehira and Ohmori (13) using Anabaena sp. strain PCC 7120 yielded 396 differentially transcribed genes 24 h following deprivation of combined nitrogen; 189 genes were up-regulated, and 207 genes were down-regulated. There are 144 and 177 N. punctiforme orthologs in the up- and down-regulated clusters of Anabaena sp. strain PCC 7120 differentially transcribed genes, respectively. However, when the Anabaena sp. strain PCC 7120 24-h differential expression database was compared to the N. punctiforme differential steady-state N₂-grown organism database, only 42 of the up-regulated and 6 of the down-regulated genes were in both databases (see File S2 in the supplemental material). Most of the common up-regulated genes fell into three categories: 11 genes for heterocyst development (including nrrA but not hetR, hetF, patB, or ntcA); 12 genes for nitrogenase catalytic proteins, assembly, or function; and 12 genes encoding either unassigned or conserved hypothetical proteins. Three of the down-regulated genes also encode conserved hypothetical proteins, plus a cofactor (NpF3770), a cell wall catabolic enzyme (NpF1408), and the ATP binding protein of a transporter (NpR2354). We note that many genes predicted to be transcriptionally enhanced in developing or functional steady-state heterocysts are present in the N. punctiforme steady-state N₂-grown database but are not present in the Anabaena sp. strain PCC 7120 24-h database. These genes include the above-mentioned developmental genes, as well as the metabolic genes glnA, zwf, opcA, pgl, and idh, the hydrogen recycling genes hupS and hupL, and the second cyo operon. Moreover, the Anabaena sp. strain PCC 7120 24-h database contains a marked number of up-regulated genes encoding ribosomal proteins and down-regulated genes encoding phycobiliproteins and photosynthetic reaction center and electron transport proteins that are not seen as differentially transcribed genes in N₂-grown N. punctiforme. A difference in the two databases could result from the fact that the N. punctiforme data are data from a steady-state culture, while the data for Anabaena sp. strain PCC 7120 were obtained 24 h following nitrogen deprivation. This rationale implies that the Anabaena sp. strain PCC 7120 cultures had not reached steady-state N₂ growth conditions by 24 h.

Functional categories of differentiating akinetes. Unlike N₂-grown cultures, 3-day akinete-forming cultures had relatively equal numbers of up- and down-regulated genes (Table 1). The majority of the 255 genes that are up-regulated as akinetes differentiate are in the unassigned category (59%). Akinete differentiation appears to involve relatively few differentially expressed adaptive proteins (11%), and the low number of these proteins may facilitate analysis of their developmental role. Transcription of avaK (NpF5452), encoding an akinete marker protein (6, 48), was enhanced, as was transcription of patA and hetF, whose products are assigned to heterocyst differentiation. The transcription of 12 signal transduction proteins was up-regulated, but only one of these proteins (encoded by NpF6364; a hybrid HK with an RR receiver domain in the N-terminal region) was from a not detectably expressed level in the NH₄⁺-grown state. Four putative transcriptional regulatory proteins were up-regulated: two DNA binding proteins (encoded by NpR3597 and NpR1565, with a Crp family binding domain and a PAS sensor domain) and two sigma subunits (encoded by NpF4153 and NpR4091, the latter of which appears to be specifically expressed in developing akinetes).

The up-regulated core metabolic genes fall into all functional subgroups. The genes involved in cell envelope and cofactor metabolism are the functional counterparts of down-regulated genes in N₂-grown cultures, but they are not homologous. Genes encoding polymer synthetic or modification activity do not appear to be essential for growth and survival and, in most cases, appear to be functionally duplicated by other genes. There was not a dramatic increase in the transcription of genes encoding proteins for carbon, nitrogen, or phosphate storage, although akinetes contain glycogen and cyanophycin granules that are visible by light microscopy (1). Presumably, the constitutive level of expression of these genes is sufficient to ensure accumulation of the polymers. None of the four genes determined by differential display to be up-regulated during akinete development (NpF0062, NpF5999, NpF6000, and NpR4070) (6) was present in the akinete transcriptome or in the NH₄⁺-grown transcriptome or proteome. NpR4070 was up-regulated in steady-state N₂-grown cultures at a significant but low level. The reasons for this discrepancy are not known, but the discrepancy may be a consequence of the particular time point chosen for analysis; the differential display results were confirmed by highly sensitive real-time quantitative PCR and reporter gene localization (6).

Also in contrast to N₂-grown cultures, the majority (64%) of the 242 genes down-regulated in akinete differentiation are involved in core metabolism. These genes include genes encoding proteins essential for energy, precursor, monomer, and polymer metabolism. The down-regulation of such genes is consistent with a cell entering a metabolically quiescent state. Nevertheless, few genes for constitutively transcribed transport proteins were down-regulated during akinete differentiation, implying that there was maintenance of an external assimilatory capacity. Nearly one-half of the down-regulated adaptive genes encode proteins for production of secondary metabolites. The most conspicuously down-regulated adaptive gene is the gene encoding HetR. HetR has been suggested to also regulate akinete differentiation in Nostoc ellipsosporum (26). However, a hetR mutant of N. punctiforme differentiates akinete-like cells that lack granulation (45). It is possible that HetR is involved in very early steps of akinete differentiation and that the gene is down-regulated by the 3-day time point. Time course transcription profiles of akinete differentiation are currently being analyzed and should resolve this possibility for N. punctiforme (M. Summers, E. Campbell, and J. Meeks, unpublished data).

Functional categories of differentiating hormogonia. This is the first systems-level assessment of differential gene expression at any stage of hormogonium development. The most striking feature of the data is the number of differentially transcribed genes in hormogonia 24 h after induction, which is nearly five times higher than the number of differentially transcribed genes in akinete-forming or N₂-fixing cultures (Table 1). This was an unpredicted result, especially the up-regulation of 944 genes in cells that have entered a nongrowth state. The impression that one gets from the hormogonium database is that the hormogonium is a clearly metabolically active, but nongrowing, motile filament that is highly tuned to sense the environment and respond to it.

Whereas most up-regulated genes were again in the unassigned category (52%), the second largest cluster was the adaptive metabolism cluster, which contained 204 genes. Slightly more than 50% of the adaptive genes encode proteins for signal transduction (85 genes) and transcriptional regulation (20 genes). On average, 85% of these genes are up-regulated uniquely in the hormogonium state. The products of the genes encoding transcriptional regulators likely account for part of the extensive differential gene expression in hormogonia at 24 h. The N. punctiforme genome contains genes encoding 156 HK, 103 RR, 33 PK, and 119 transcriptional regulators (31, 34). Some of the sensory proteins are organized with a multiple-domain hybrid kinase and/or regulator architecture. In the cluster of up-regulated hormogonium sensory proteins there are four hybrid HK-PKs and 24 HK-RRs, and in 13 of the HK-RRs the receiver domain is in the N terminus of the protein, where it likely functions as a sensor domain. An alternative sigma subunit (sigH [NpR1771]), which is induced during coculture of N. punctiforme with a symbiotic plant partner and whose mutation results in an altered infection frequency (8), was not up-regulated in the combined-nitrogen starvation-induced 24-h hormogonia. This result implies that hormogonia may be subtly different depending on the environmental signal for their induction.

In addition, 18 genes encoding putative chemotaxis proteins from five separate genomic loci, five gas vesicle proteins, and five pilus formation proteins were up-regulated. Up-regulation of the gas vesicle and pilus genes confirms results for hormogonia of Calothrix sp. strain PCC 7601 (10). Cyanobacterial gliding motility has been reported to be dependent on type IV pilus retraction (7) and/or polysaccharide secretion (23). There are several gene clusters encoding proteins of polysaccharide synthesis assigned to core metabolism (cell envelope synthesis) and development (heterocyst envelope synthesis); eight genes in the latter group were also up-regulated in hormogonia. A cluster of genes that encode polysaccharide biosynthetic and type II secretion proteins was up-regulated uniquely in hormogonia (NpF0066 to NpF0073; NpF0073 [pulG] was also up-regulated in akinetes). If polysaccharide extrusion is a mechanism of gliding motility in hormogonia, there appears to be sufficient biosynthetic capacity to support it, although there is currently no experimental evidence for either the motility mechanisms or specific polysaccharide synthetic genes.

Since hormogonia in this study were induced by transfer from NH₄⁺ into combined-nitrogen-free medium, we anticipated some overlap in the hormogonium differentiation and N₂-grown databases. Seventeen developmental and 18 alternative nitrogen (N) source genes were up-regulated in 24-h hormogonia. The developmental genes include the akinete marker avaK and the heterocyst-specific genes devH, hetF, and patU. Hormogonia are the infective units of Nostoc sp.-plant symbioses (30), and an N. punctiforme hetF mutant infects a symbiotic plant partner (45), which implies that HetF, at least, is not essential for hormogonium differentiation. The alternative N source genes include the genes encoding nitrate (NpR1526) and nitrite (NpF1528) reductases and a nitrate/nitrite permease (NpR1527), as well as 11 nitrogenase-associated proteins, including the catalytic subunit proteins NifH (NpR0415), NifD (NpR0414), and NifK (NpR0390). As in essentially all heterocyst-forming cyanobacteria, the N. punctiforme nifD gene is interrupted by an insertion element (31, 34). Northern hybridization established that the nifD element was excised in 24-h hormogonia and that nifK can also be transcribed in the absence of the excision rearrangement (data not shown). There is evidence, however, that hormogonia do not fix N₂ in normal atmospheric conditions (9), implying that an active protein complex is not assembled. We suggest that these alternative N source genes and some of the developmental genes are expressed in response to the NH₄⁺ starvation method of hormogonium induction. Heterocysts were not present in the 24-h hormogonia, and their differentiation does not occur until 48 to 56 h later. Moreover, HetR, which is essential for heterocyst differentiation, was strongly down-regulated during hormogonium development (see below), indicating that a complete developmental pathway was not operational in 24-h hormogonia.

A substantial number of core metabolic genes (179 genes) were up-regulated in 24-h hormogonia. Most of these genes are involved in cell envelope, polymer, monomer, cofactor, and precursor metabolism. The up-regulated cell envelope proteins include the cell shape-determining proteins MreC (NpR1840) and MreD (NpR1839); MreB (NpR1841) is constitutively expressed in the NH₄⁺-grown transcriptome and proteome. The presence of these proteins may account for the change in shape of hormogonium cells. The polymer proteins are enriched in DNA metabolism and protein processing. The 16 adaptive and core DNA metabolic genes include the genes encoding five helicases, four nucleases, three proteins for modification, two proteins for replication (polA [NpF1050] and priA [NpR0306]), and two proteins for replication repair (sbcC [NpF4484] and ssb [NpF5942]). The presumptive presence of these proteins indicates that there is some DNA metabolism in hormogonia, perhaps chromosome stabilization, although there is no net accumulation of DNA. There are 13 genes encoding protein processing factors and 12 proteases in the adaptive and core categories. The chaperones include two different DnaKs (NpF0122 and NpF5907), two DnaJs (NpF0150 and NpR3872), and two IbpAs (Hsp20) (NpF0198 and NpF0784). A protease whose structural gene was significantly up-regulated is NblA (NpF0291), which degrades phycobilisomes and accounts for the decrease in phycobiliprotein-induced fluorescence in hormogonia. The hormogonium cultures were incubated in the absence of an environmental source of combined nitrogen. Therefore, the products of phycobiliprotein degradation most likely supply amino nitrogen for synthesis of adaptive proteins, aided by the action of the products of the up-regulated amino acid biosynthetic genes. Thus, substantial protein metabolism by reorganization and processing may occur in hormogonia.

Nearly 50% of the 883 genes down-regulated in 24-h hormogonia encode proteins of core metabolism, and most of the remainder cluster in the unassigned metabolism category (38%). Except for the cell division genes, every subgroup of core metabolic genes was substantially down-regulated, and the down-regulation was much more extensive than that of akinetes 3 days into the differentiation process. These results are consistent with a stringent nongrowth state of hormogonia. Notably down-regulated were two genes encoding the chaperone GroEL (NpF1230 and NpR0829) and one gene encoding GroES (NpR0830). Many genes for energy metabolism (e.g., ATP synthase) and photosynthetic and respiratory electron transport were down-regulated, although their transcripts remained detectable. A specific Fe-S protein (putative PsaC; encoded by NpF5213) has been shown to be down-regulated in hormogonia using a reporter gene (5). While nongrowing, hormogonia maintain a measurable rate of photosynthetic CO₂ assimilation (9). The photosynthetic rate may be sufficient to provide photosynthate for polysaccharide synthesis and extrusion or energy for pilus assembly and retraction for hormogonium gliding.

Almost as many genes encoding transport proteins were down-regulated as were up-regulated. This observation may imply that hormogonia acquire different classes of compounds from the environment than NH₄⁺-grown filaments acquire. There is a smaller, but still numerically substantial, fraction of adaptive genes that were down-regulated in hormogonia. Approximately 35% as many genes encoding signal transduction proteins were down-regulated as were up-regulated. This also suggests that there are major specific changes in the signal transduction pathways as hormogonia differentiate from NH₄⁺-grown filaments. The notable down-regulated genes include patA and hetR involved in heterocyst differentiation and glnB (NpF4466) encoding P_II, which has an intracellular nitrogen- and carbon-sensing role in unicellular cyanobacteria (16). The down-regulation of 10 transcriptional regulators, coupled with the up-regulation of 20 similar proteins, indicates that there is a shift in the cellular composition of transcription factors concurrent with a change in physiology. Included in the down-regulated adaptive category are four genes (NpF2182, NpF2183, NpF2186, and NpF2187) located in a cluster of genes encoding proteins for the synthesis of a cyclic peptide (nostopeptide) (nosB, -C, -E, and -F). The differential expression of these genes in hormogonia and other genes in akinetes and N₂-grown cultures may lead to an understanding of the physiological role(s) of such secondary metabolites.

The differential expression in hormogonia of seven genes (up-regulated genes NpF0152 [tonB], NpF2507 [pilT], NpF3829 [atoC], NpF4811 [fliA], NpF5963 [che, encoding methyl-accepting chemotaxis protein], and NpR0414 [nifD] and down-regulated gene NpR0830 [groES]) was verified by Northern hybridization (data not shown).

Overlap of genes differentially expressed in the developmental states. The extent of overlap of up- and down-regulated genes between developmental states is depicted in the Venn diagrams in Fig. 3. This depiction does not include the genes that were up-regulated in one developmental state in any of the two- or three-way comparisons and down-regulated in another developmental state and the converse. A detailed list of the up- and down-regulated genes shared by all three developmental states and by the akinete-forming state and the N₂ steady state is presented in Table 2.

View larger version (20K): [in this window] [in a new window]   FIG. 3. Numerical relationships between up- and down-regulated genes in the N2-grown, hormogonium, and akinete developmental states.

Conclusions. The collective data indicate that the three alternative developmental states of N. punctiforme do have some genes in common, but their differentiation largely involves up-regulation of genes distinct for each state. The most remarkable result is the extent of differential gene expression by 24-h hormogonia. These filaments were once characterized as growth precursor cells (9, 11). Our data support such a conclusion; hormogonia have a distinctly down-regulated capacity for core metabolic functions, similar to the swarmer cells of Caulobacter crescentus (24, 25). Nevertheless, hormogonia are transcriptionally very active, especially in genes encoding adaptive metabolic environmental sense and response functions. Hormogonia appear to be poised to move in specific directions, perhaps depending on the light quality and quantity or the presence of a chemoattractant, to a microhabitat which is somehow perceived as permissive for subsequent growth and reproduction.

Including the hetR results, there is little evidence in the differential expression profiles of N₂-grown cultures and the profiles for 3 days into the differentiation of akinetes to support an evolutionary sequence of akinete to heterocyst in the ancestors of N. punctiforme by, for example, gene duplication. It is possible that such evidence will emerge from comparative time course experiments during differentiation, but it is also possible that too much evolutionary divergence has occurred in developmental pathways of the two cell types.

	ACKNOWLEDGMENTS

 
This work was supported by grant EF0317104 from the NSF to J.C.M. and by grant GM048680 from the NIH to M.L.S.

We thank Ian Campbell and Michael Goodson for assistance in database and statistical analyses, Jeff Elhai for sharing data, and Michael Schultz for help in the development of hybridization procedures. We also thank Dawei Lin and Brad Sickler of the UC Davis Genome Center and Bioinformatics Program for collaboration in database analysis and visualization tools.

	FOOTNOTES

 
* Corresponding author. Mailing address: Section of Microbiology, University of California, Davis, CA 95616. Phone: (530) 752-3346. Fax: (530) 752-9014. E-mail: jcmeeks{at}ucdavis.edu

Published ahead of print on 4 May 2007.

Supplemental material for this article may be found at http://jb.asm.org/.

	REFERENCES

Top ABSTRACT INTRODUCTION MATERIALS AND METHODS RESULTS AND DISCUSSION REFERENCES

 

1. Adams, D. G., and P. S. Duggan. 1999. Heterocyst and akinete differentiation in cyanobacteria. New Phytol. 144:1-33.[CrossRef]
2. Anderson, D. C., E. L. Campbell, and J. C. Meeks. 2006. 3D LC/MS/MS analysis of the soluble proteome of the cyanobacterium Nostoc punctiforme ATCC 29133. J. Proteome Res. 5:3096-3104.[CrossRef][Medline]
3. Aravind, L., and E. V. Koonin. 2002. Classification of the caspase-hemoglobinase fold: detection of new families and implications for the origin of eukaryotic separins. Proteins Struct. Funct. Genet. 46:355-367.[CrossRef][Medline]
4. Argueta, C., and M. L. Summers. 2005. Characterization of a model system for the study of Nostoc punctiforme akinetes. Arch. Microbiol. 183:338-346.[CrossRef][Medline]
5. Argueta, C., K. Yuksek, and M. L. Summers. 2004. Construction and use of GFP reporter vectors for analysis of cell-type-specific gene expression in Nostoc punctiforme. J. Microbiol. Methods 59:181-188.[CrossRef][Medline]
6. Argueta, C., K. Yuksek, R. Patel, and M. L. Summers. 2006. Identification of Nostoc punctiforme akinete-expressed genes using differential display. Mol. Microbiol. 61:748-757.[CrossRef][Medline]
7. Bhaya, D., N. R. Bianco, D. Bryant, and A. Grossman. 2000. Type IV pilus biogenesis and motility in the cyanobacterium Synechocystis sp. PCC6803. Mol. Microbiol. 37:941-951.[CrossRef][Medline]
8. Campbell, E. L., B. Brahamsha, and J. C. Meeks. 1998. Mutation of an alternative sigma factor in the cyanobacterium Nostoc punctiforme results in increased infection of its symbiotic plant partner, Anthoceros punctatus. J. Bacteriol. 180:4938-4941.[Abstract/Free Full Text]
9. Campbell, E. L., and J. C. Meeks. 1989. Characteristics of hormogonia formation by symbiotic Nostoc spp. in response to the presence of Anthoceros punctatus or its extracellular products. Appl. Environ. Microbiol. 55:125-131.[Abstract/Free Full Text]
10. Damerval, T., G. Guglielmi, J. Houmard, and N. Tandeau de Marsac. 1991. Hormogonium differentiation in the cyanobacterium Calothrix: a photoregulated developmental process. Plant Cell 3:191-201.[Abstract/Free Full Text]
11. Dow, C. S., R. Whittenbury, and N. G. Carr. 1983. The "shut down" or "growth precursor" cell—an adaptation for survival in a potentially hostile environment, p. 187-247. In J. H. Slater, R. Whittenbury, and J. W. T. Wimpenny (ed.), Microbes in their natural environments. Cambridge University Press, Cambridge, United Kingdom.
12. Ehira, S., M. Ohmori, and N. Sato. 2003. Genome wide expression of the responses to nitrogen deprivation in the heterocyst-forming cyanobacterium Anabaena sp. strain PCC 7120. DNA Res. 10:97-113.[Abstract]
13. Ehira, S., and M. Ohmori. 2006. NrrA, a nitrogen-responsive response regulator facilitates heterocyst development in the cyanobacterium Anabaena sp. strain PCC 7120. Mol. Microbiol. 59:1692-1703.[CrossRef][Medline]
14. Enderlin, C. S., and J. C. Meeks. 1983. Pure culture and reconstitution of the Anthoceros-Nostoc symbiotic association. Planta 158:157-165.[CrossRef]
15. Flores, E., and A. Herrero. 1994. Assimilatory nitrogen metabolism, p. 487-517. In D. A. Bryant (ed.), The molecular biology of cyanobacteria. Kluwer Academic Publishers, Boston, MA.
16. Forchhammer, K. 2004. Global carbon/nitrogen control by P_II signal transduction in cyanobacteria: from signals to targets. FEMS Microbiol. Rev. 28:319-333,[CrossRef][Medline]
17. Golden, J. W., and H. S. Yoon. 2003. Heterocyst development in Anabaena. Curr. Opin. Microbiol. 6:557-563.[CrossRef][Medline]
18. Hagen, K. D., and J. C. Meeks. 2001. The unique cyanobacterial protein OpcA is an allosteric effector of glucose-6-phosphate dehydrogenase in Nostoc punctiforme ATCC 29133. J. Biol. Chem. 276:11477-11486.[Abstract/Free Full Text]
19. Herdmann, M., M. Janvier, R. Rippka, and R. Y. Stanier. 1979. Genome size in cyanobacteria. J. Gen. Microbiol. 111:73-85.
20. Herdmann, M., and R. Rippka. 1988. Cellular differentiation: hormogonia and baeocytes. Methods Enzymol. 167:232-242.
21. Herrero, A., A. M. Muro-Pastor, and E. Flores. 2001. Nitrogen control in cyanobacteria. J. Bacteriol. 183:411-425.[Free Full Text]
22. Herrero, A., A. M. Muro-Pastor, A. Valladares, and E. Flores. 2004. Cellular differentiation and the NtcA transcription factor in filamentous cyanobacteria. FEMS Microbiol. Rev. 28:469-487.[Medline]
23. Hoiczyk, E., and W. Baumeister. 1998. The junctional pore complex, a prokaryotic secretion organelle, is the molecular motor underlying gliding motility in cyanobacteria. Curr. Biol. 8:1161-1168.[CrossRef][Medline]
24. Jacobs-Wagner, C. 2004. Regulatory proteins with a sense of direction: cell cycle signaling network in Caulobacter. Mol. Microbiol. 51:7-13.[CrossRef][Medline]
25. Laub, M., H. H. McAdams, T. Feldblyum, C. M. Fraser, and L. Shapiro. 2000. Global analysis of the genetic network controlling a bacterial cell cycle. Science 290:2144-2148.[Abstract/Free Full Text]
26. Leganés, F., F. Fernandéz-Piñas, and C. P. Wolk. 1994. Two mutations that block heterocyst differentiation have different effects on akinete differentiation in Nostoc ellipsosporum. Mol. Microbiol. 12:679-684.[CrossRef][Medline]
27. Liang, J., L. Scappino, and R. Haselkorn. 1992. The patA gene product, which contains a region similar to CheY of Escherichia coli, controls heterocyst pattern formation in the cyanobacterium Anabaena 7120. Proc. Natl. Acad. Sci. USA 89:5655-5659.[Abstract/Free Full Text]
28. Lynn, M. E., J. A. Bantle, and J. D. Ownby. 1986. Estimation of gene expression in heterocysts of Anabaena variabilis by using DNA-RNA hybridization. J. Bacteriol. 167:940-946.[Abstract/Free Full Text]
29. Martin-Figueroa, E., F. Navaro, and F. J. Florencio. 2000. The GS-GOGAT pathway is not operative in the heterocysts. Cloning and expression of the glsF gene from the cyanobacterium Anabaena sp. PCC 7120. FEBS Lett. 476:282-286.[CrossRef][Medline]
30. Meeks, J. C. 1998. Symbiosis between nitrogen-fixing cyanobacteria and plants. BioScience 48:266-276.[CrossRef]
31. Meeks, J. C. 2005. The genome of the filamentous cyanobacterium Nostoc punctiforme. What can we learn from it about free-living and symbiotic nitrogen fixation?, p. 27-70. In R. Palacios and W. E. Newton (ed.), Nitrogen fixation: 1888-2001, vol. VI. Genomes and genomics of nitrogen-fixing organisms. Springer, Dordrecht, The Netherlands.
32. Meeks, J. C., E. L. Campbell, M. L. Summers, and F. C. Wong. 2002. Cellular differentiation in the cyanobacterium Nostoc punctiforme. Arch. Microbiol. 178:395-403.[CrossRef][Medline]
33. Meeks, J. C., and J. Elhai. 2002. Regulation of cellular differentiation in filamentous cyanobacteria in free-living and plant-associated symbiotic growth states. Microbiol. Mol. Biol. Rev. 65:94-121.
34. Meeks, J. C., J. Elhai, T. Thiel, M. Potts, F. Larimer, J. Lamerdin, P. Predki, and R. Atlas. 2001. An overview of the genome of Nostoc punctiforme, a multicellular, symbiotic cyanobacterium. Photosynth. Res. 70:85-106.[CrossRef][Medline]
35. Renström-Kellner, E., A. N. Rai, and B. Bergmann. 1990. Correlation between nitrogenase and glutamine synthetase expression in the cyanobacterium Anabaena cylindrica. Physiol. Plant. 80:12-19.[CrossRef]
36. Rippka, R., J. Deruelles, J. B. Waterbury, M. Herdmann, and R. Y. Stanier. 1979. Generic assignments, strain histories and properties of pure cultures of cyanobacteria. J. Gen. Microbial. 111:1-61.
37. Smyth, G. K. 2005. Limma: linear models for microarray data, p. 397-420. In R. Gentleman, V. Carry, S. Dudoit, R. Irizarry, and W. Huber (ed.), Bioinformatics and computational biology solutions using R and Bioconductor. Springer, New York, NY.
38. Summers, M. L., and J. C. Meeks. 1996. Transcriptional regulation of zwf, encoding glucose-6-phosphate dehydrogenase, from the cyanobacterium Nostoc punctiforme ATCC 29133. Mol. Microbiol. 22:473-480.[CrossRef][Medline]
39. Tandeau de Marsac, N. 1994. Differentiation of hormogonia and relationships with other biological processes, p. 825-842. In D. A. Bryant (ed.), The molecular biology of cyanobacteria. Kluwer Academic Publishers, Boston, MA.
40. Thiel, T., and C. P. Wolk. 1983. Metabolic activities of spores of Nostoc spongiaeforme. J. Bacteriol. 156:369-374.[Abstract/Free Full Text]
41. Valladares, A., A. Herrero, D. Pils, G. Schmetterer, and E. Flores. 2003. Cytochrome c oxidase genes required for nitrogenase activity and diazotrophic growth in Anabaena sp. PCC 7120. Mol. Microbiol. 47:1239-1249.[CrossRef][Medline]
42. Winkenbach, F., and C. P. Wolk. 1973. Activities of enzymes of the oxidative and the reductive pentose phosphate pathways in heterocysts of a blue-green alga. Plant Physiol. 52:480-483.[Abstract/Free Full Text]
43. Wolk, C. P. 2000. Heterocyst formation in Anabaena, p. 83-104. In Y. V. Brun and L. J. Shimkets (ed.), Prokaryotic development. American Society for Microbiology, Washington, DC.
44. Wolk, C. P., A. Ernst, and J. Elhai. 1994. Heterocyst metabolism and development, p. 769-823. In D. A. Bryant (ed.), The molecular biology of cyanobacteria. Kluwer Academic Publishers, Boston, MA.
45. Wong, F. C., and J. C. Meeks. 2002. Establishment of a functional symbiosis between the cyanobacterium Nostoc punctiforme and the bryophyte hornwort Anthoceros punctatus requires genes involved in nitrogen control and initiation of heterocyst differentiation. Microbiology 148:315-323.[Abstract/Free Full Text]
46. Xu, X., J. Elhai, and C. P. Wolk. 2007. Transcriptional and developmental responses by Anabaena to deprivation of fixed nitrogen, p. 383-422. In A. Herrero and E. Flores (ed.), The cyanobacteria: molecular biology, genomics and evolution. Horizon Scientific Press, Norwich, United Kingdom.
47. Zhang, C.-C., S. Laurent, S. Sakr, L. Peng, and S. Bédu. 2006. Heterocyst differentiation and pattern formation in cyanobacteria: a chorus of signals. Mol. Microbiol. 59:367-375.[CrossRef][Medline]
48. Zhou, R., and C. P. Wolk. 2002. Identification of an akinete marker gene in Anabaena variabilis. J. Bacteriol. 184:2529-2532.[Abstract/Free Full Text]

This article has been cited by other articles:

• Stockel, J., Welsh, E. A., Liberton, M., Kunnvakkam, R., Aurora, R., Pakrasi, H. B. (2008). Global transcriptomic analysis of Cyanothece 51142 reveals robust diurnal oscillation of central metabolic processes. Proc. Natl. Acad. Sci. USA 105: 6156-6161 [Abstract] [Full Text]  
• Chapman, K. E., Duggan, P. S., Billington, N. A. née, Adams, D. G. (2008). Mutation at Different Sites in the Nostoc punctiforme cyaC Gene, Encoding the Multiple-Domain Enzyme Adenylate Cyclase, Results in Different Levels of Infection of the Host Plant Blasia pusilla. J. Bacteriol. 190: 1843-1847 [Abstract] [Full Text]  
• Adams, D. G., Duggan, P. S. (2008). Cyanobacteria-bryophyte symbioses. J Exp Bot 59: 1047-1058 [Abstract] [Full Text]  

This Article Abstract Full Text (PDF) Supplemental material Other Versions of this Article: JB.00360-07v1 189/14/5247    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Reprints and Permissions Copyright Information Books from ASM Press MicrobeWorld Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Campbell, E. L. Articles by Meeks, J. C.