BACKGROUND
n. gonorrhoeae was long considered to be an obligate aerobe until it was discovered that anaerobic growth was possible when nitrite or nitric oxide  was used as a terminal electron acceptor  <cit> .
anaerobic growth is accomplished through utilization of a truncated denitrification pathway, which is encoded within the gonococcal genome as a pair of divergently transcribed genes, ania, encoding a nitrite reductase, and norb, encoding a nitric oxide reductase  <cit> .
anaerobiosis is presumed to be a physiologically significant state during infection, as the gonococcus is often recovered from infected individuals in co-culture with obligate anaerobes such as peptococcus and bacteroides spp.
<cit> .
furthermore, gonococci have been shown to induce and repress the expression of several genes in response to anaerobiosis, and antibody to ania, the major anaerobically induced outer membrane protein, can be found in sera from infected women, demonstrating that this protein is expressed in vivo  <cit> .
the ability of gonococci to utilize this denitrification pathway to reduce no may have immunomodulary effects during the course of infection.
some evidence suggests that the gonococcal reduction of host-produced no may be responsible, at least in part, for the high incidence of asymptomatic disease  <cit> .
in vitro, n. gonorrhoeae was shown to be capable of setting a no steady-state level in the anti-inflammatory range  <cit> .
recently, in an in vitro cervical cell model of infection, it was shown that gonococcal activation of inos promoted bacterial survival.
in this study, it was suggested that host-derived nitric oxide is not protective against gonococci, rather, nitric oxide may actually be required to sustain cervical bacterial disease  <cit> .
previous research has shown that the genes involved in denitrification and/or adaptation to anaerobic growth in neisseria spp.
are subject to transcriptional regulation by the oxygen-sensitive regulator, fnr, the no-sensitive repressor, nsrr, and the nitrite-insensitive two component system, narqp  <cit> .
earlier data from a microarray-based approach suggested that the gonococcal fnr-regulon was composed of fourteen activated and six repressed transcripts, making the gonococcal fnr regulon much smaller than that of e. coli, where fnr was shown to regulate more than  <dig> operons  <cit> .
discovery of nsrr and narqp-regulated genes in neisseria spp.
has relied mainly on in silico analysis of intergenic chromosomal regions to identify similarity to previously defined regulator binding sites.
as is the case with fnr, the neisserial nsrr and narqp regulons, at least the currently defined members, are much smaller than their e. coli counterparts  <cit> .
in this study we employ a powerful whole-genome approach, rna-seq, to quantitatively sequence the complete gonococcal transcriptome.
using this method we were able to define global changes in gene expression in response to anaerobiosis.
we show that the gonococcal anaerobic stimulon is not small, and that  <dig> chromosomal open reading frames  are differentially expressed.
we present several novel findings that, taken together, support the view that anaerobic growth is an important facet of life for this organism, and should be considered when studying the host/pathogen interaction.
RESULTS
sequencing the gonococcal transcriptome
whole genome mrna sequencing is an attractive method of monitoring global changes in gene expression while overcoming many of the pitfalls of traditional dna microarrays  <cit> .
for the purposes of this study, anaerobic conditions were defined as anoxia, the presence of nitrite  and the concomitant presence of nitric oxide .
in order to define the gonococcal anaerobic stimulon, rna-seq was performed on two biological replicates of aerobically or anaerobically plate-grown gonococci using the abi solid™ system .
unique sequence reads from the rna-seq data were mapped to the annotated fa <dig> genome, and gene expression was quantified as reads per kilobase of coding sequence per million reads  .
when expression data for each replicate were plotted against each other, rpkm values were observed to be adequately reproducible , though expectedly, the extent of reproducibility was slightly less in genes with low expression  <cit> .
for this reason, we required mapped genes to have an rpkm ≥  <dig> in order to be considered as a candidate for differential expression.
however, genes that had an rpmk <  <dig> under one growth condition, but were highly induced in the other condition were accepted as candidates.
genes were considered to be differentially expressed if there was a three-fold or greater difference in rpkm between the two growth conditions.
analysis of rna-seq data revealed that  <dig> chromosomal genes met the criteria for being differentially expressed, with  <dig> genes being anaerobically induced and  <dig> genes being anaerobically repressed .
figure 1b and 1c illustrate a reasonable measure of reproducibility in expression level between replicates in the  <dig> differentially expressed genes.
a gonococci grown in an anaerobe chamber on plates containing  <dig> mm nitrite.
b the fold change for each replicate was calculated by comparing rpkm values of anaerobically grown gonococci to rpkm values of aerobically grown gonococci.
a positive value represents an induction of gene expression anaerobically, and a negative value represents a repression of gene expression anaerobically
c protein definitions are derived from the annotated fa <dig> genome , and/or the ncbi conserved domain database  <cit> .
d  genes that have been shown to be directly regulated by fur in sebastian et al.,   <cit>  and/or jackson et al.,   <cit> .
genes reported to be upregulated during rpoh overexpression in gunesekere et al.,   <cit> .
genes reported to be in the lexa regulon by schook et al.,   <cit> .
genes found to be differentially expressed in a n. meningitidis Δlrp mutant  by microarray analysis in ren et al.,   <cit> .
genes found to be differentially expressed in a gonococcal Δmtrr mutant by microarray analysis in folsterΔ et al.,   <cit> .
genes reported to be in the farr regulon by friedrich et al.,   <cit> .
genes reported to be in the nmlr regulon by kidd et al.,   <cit> .
genes found to be differentially expressed microaerobically in a Δfnr mutant with and without nitrite by microarray analysis in whitehead et al.,   <cit> .
genes reported to be in the nsrr regulon by isabella et al.,   <cit>  and isabella,   <cit> .
e fur-repressed small regulatory rna described by ducey et al.,   <cit> .
as expected, anaerobic expression levels cluster at higher rpkm values in the set of anaerobically induced genes , while the opposite is true in anaerobically repressed genes .
secondary verification of differentially expressed genes
many of the genes found to have been differentially expressed by rna-seq have previously been determined to be members of the neisserial fnr or nsrr regulons, validating the use of this technique to monitor transcriptome changes   <cit> .
to further corroborate the rna-seq data, several genes shown to be differentially expressed were selected for secondary confirmation.
analysis of translational promoter-lacz fusions or rt-pcr was used to accomplish this task .
analysis of translational promoter-lacz fusions demonstrated that the genes encoding heat shock protein, clpb , dna repair enzyme, recn , nitropropane dioxygenase, npd , filamentous hemagglutinin, heca , glutaredoxin, grx <dig> , bacterioferritin, bfra  and repressor protein lexa  were significantly upregulated anaerobically .
rt-pcr confirmed that 2-isopropylmalate synthase, leu <dig> , and adhesin, leca , were also increased in expression anaerobically .
translational promoter-lacz fusions to the genes encoding the iron-binding protein, fbpa , and the solute binding protein, oxia  confirmed that these genes are anaerobically repressed .
rt-pcr confirmed that ngo <dig>  l-lactate dehydrogenase, lldd , thioredoxin, trx <dig> , heme oxygenase, hemo , and nadh dehydrogenase component, nqra  were also anaerobically repressed .
analysis of raw rna-seq data allows for the prediction of transcription start sites, - <dig>  and - <dig> elements, as has been described in rna-seq projects in both helicobacter pylori and listeria monocytogenes  <cit> .
compared to e. coli, the use of alternative sigma factors in the gonococcus is very limited.
gonococci have no sigma- <dig> homolog, and only utilize sigma factors in the sigma- <dig> family.
this class of sigma factors recognize - <dig> and - <dig> elements within bacterial promoters  <cit> .
the predictions of transcriptional elements for the individual genes selected for secondary analysis are shown in figure  <dig>  these predictions were generally a high match to the e. coli consensus for transcriptional elements  or the α-proteobacteria consensus for rpoh-dependent transcriptional elements   <cit> .
the ability to predict the location of these transcriptional elements in a large scale manner will prove very useful in the search for potential regulatory sites, and will aid in future work to define gonococcal transcriptional networks.
differentially expressed genes involved in macromolecular biosynthesis
the genes found to be differentially expressed in response to anaerobic growth were broadly characterized according to their putative function .
several genes encoding proteins involved in macromolecular synthesis displayed decreased expression under anaerobic conditions.
the synthesis of ribosomal proteins is strongly related to growth rate, and ribosomal protein synthesis has been documented to decrease during times of energy deficiency  <cit> .
for facultative anaerobes, growth in an oxygen-limited environment is more energetically-deprived than in aerobic conditions.
the reduced synthesis of ribosomal proteins is indicative of adaptation to a slower growth rate following recovery from a nutritional shift-down.
the reduced expression of rna polymerase β an β' subunits is also indicative of this , and has been observed in e. coli as well  <cit> .
genes involved in transport and binding
many genes encoding proteins involved in transport and binding were found to be anaerobically repressed .
most of these proteins are involved in the transport of amino acids, sugars, or other growth metabolites into the cell.
this result comes as no surprise, as a slower anaerobic growth rate would result in decreased demand and slower utilization of such substrates.
proteins of this type were also shown to be downregulated anaerobically in e. coli, p. aeruginosa, and s. enterica  <cit> .
three repressed genes, fbpa, hemo, and ompu, are involved in the acquisition and/or transport of iron into the cell.
in s. enterica, expression of the gene encoding iron acquisition protein, sita, was also found to be repressed anaerobically  <cit> .
while genes involved in iron acquisition and transport are repressed anaerobically, the genes encoding the iron-storage bacterioferritin proteins, bfra and bfrb are induced .
three other proteins implicated in binding that were shown to be anaerobically induced, leca , heca , and fhbp, may be involved directly in the host-pathogen interaction.
the role of leca in gonococci is unclear, but a leca homolog in another member of the neisseriaceae family, eikenella corrodens, has been shown to be a part of an adhesin complex important in adherence to a variety of host tissues  <cit> .
the gonococcal heca protein contains a hemagglutination domain that, in other bacteria, has been shown to promote bacterial aggregation or host attachment  <cit> .
the role of factor h binding protein  as it relates to gonococcal pathogenesis has not been examined, however, in n. meningitidis, fhbp expression on the cell surface was shown to be a means of immune evasion that rendered meningococci resistant to complement-mediated killing  <cit> .
genes involved in energy production/oxidation-reduction
unsurprisingly, numerous genes encoding proteins involved in energy metabolism were found to be anaerobically repressed .
gonococci use a combination of the entner-douderoff and pentose phosphate pathways in the utilization of glucose  <cit> .
key enzymes involved in these pathways were found to be anaerobically repressed, as well as several enzymatic components of the krebs cycle.
anaerobic repression of genes involved in central intermediary metabolism has also been observed in other bacteria, including e. coli, p. aeruginosa, and b. subtilus  <cit> .
interestingly, gonococci contain two major nadh dehydrogenases capable of coupling the oxidation of nadh to the reduction of ubiquinone.
one of these dehydrogenases is the prototypical h+-translocating 14-subunit complex encoded by the nuo operon .
the second is a seven-subunit na+-translocating complex encoded by the nqr operon .
though genes in the nuo operon did not meet the criteria used in this study to define differential expression, there was an approximate two-fold level of anaerobic repression across its entire length .
the nqr operon, however, demonstrated a significant decrease in anaerobic expression .
na+-translocating nadh dehydrogenases were originally discovered and characterized in marine bacteria, though it has now been revealed that several pathogenic species encode such complexes  <cit> .
transport of sodium ions into the periplasmic space through the nqr complex generates an electrochemical membrane potential that plays a role in solute import, atp synthesis, and, for organisms that are motile, flagellar rotation  <cit> .
the flow of electrons through the active portions of the electron transport chain would be expected to decrease anaerobically when:  <cit>  oxygen is unavailable, rendering the terminal cytochrome oxidase  inactive,  <cit>  enzymatic components of the krebs cycle are repressed, resulting in decreased production of reducing equivalents, and  <cit>  growth is dependent on high potential alternative electron acceptors such as no2- and no.
furthermore, neither the nitrite or nitric oxide reductases are thought to be capable of pumping protons across the membrane  <cit> .
downregulation of nqr would decrease the electrochemical potential across the membrane , thus decreasing the proton motive force  <cit> .
because expression of the nuo operon is not affected to the same degree anaerobically, this would now result in a greater proportion of nadh oxidation being coupled to atp production through h+ translocation.
also, many proteins involved in the transport of metabolites into the cell are repressed anaerobically.
downregulation of the nqr complex may add another level of regulation to cellular transport.
it stands to reason that the activity of solute transporters that rely on the sodium motive force would decrease when fewer na+ ions are pumped across the membrane  <cit> .
the function and regulation of the most highly induced gene involved in anaerobic growth, the nitrite reductase ania, has been thoroughly characterized, and has been shown to be required for anaerobic growth on nitrite  <cit> .
interestingly, gonococci contain a gene encoding a 2-nitropropane dioxygenase-like protein that was found to be induced anaerobically .
this would seem unusual, as gonococci would not be expected to encounter the toxic nitroalkanes utilized by npd proteins.
furthermore, molecular oxygen is consumed in the npd reaction mechanism  <cit> .
further work should determine if this protein plays a novel role in gonococci, possibly microaerobically when oxygen would be present.
gonococci also showed increased anaerobic expression of a gene encoding electron transfer flavoprotein-ubiquinone oxidoreductase, eftd , a protein containing 57% identity  to human eftd.
the presence of eftd in bacteria is rare, and its role has only been studied in plants and mammals.
in arabidopsis thaliana, eftd was shown to be upregulated in times of carbon starvation, when genes involved in glycolysis and the krebs cycle were reduced in expression, as is observed when gonococci are grown anaerobically  <cit> .
in a. thaliana, it was reported that eftd can act as an electron acceptor from various dehydrogenases involved in amino acid catabolism and fatty acid degradation.
the reduction of ubiquinone by etfd was proposed to be a mechanism of utilizing alternative respiratory substrates to feed into the electron transport chain  <cit> .
eftd upregulation in n. gonorrhoeae could indicate a shift to the use of alternative growth substrates under anaerobic conditions.
genes involved in small molecule biosynthesis
anaerobically, several genes involved in the synthesis of amino acids were found to be repressed , which has also been observed in other bacteria  <cit> .
genes involved in the synthesis of iron-sulfur clusters were also repressed.
this result is not entirely unexpected, as many iron-sulfur cluster containing proteins involved in respiration were shown to be downregulated anaerobically.
the demand for iron-sulfur clusters may therefore be reduced under anaerobic conditions.
also, because it is common for the fe-s clusters in fe-s cluster-containing proteins to be oxygen-sensitive, a lower rate of oxygen-induced fe-s cluster turnover would be expected under anaerobic conditions  <cit> .
unlike other amino acid biosynthetic genes, the gene encoding leu <dig>  the first committed step in leucine biosynthesis, was found to be induced anaerobically .
this could be indicative of a decrease in intracellular leucine concentration, as leu <dig>  is repressed by leucine in e. coli and other organisms  <cit> .
in e. coli, the anaerobic stimulon is intimately associated with global regulation through the leucine responsive regulatory protein, lrp.
<cit> .
in n. meningitidis, the reported lrp regulon demonstrated striking overlap to many of the genes in the gonococcal anaerobic stimulon , particularly in regards to proteins involved in glucose utilization, the krebs cycle, and ribosomal assembly  <cit> .
future work will examine whether gonococcal lrp plays a role in control of the anaerobic stimulon.
genes involved in transcriptional regulation
the expression of several regulatory proteins was found to be induced under anaerobic conditions, though no regulatory proteins appeared to be repressed .
the gene encoding arsr , a regulator involved in norb repression in response to iron, was shown to be moderately upregulated  <cit> .
it is unclear if arsr regulates any other target genes.
the fnr-regulated nmlr protein has been shown to act as both an activator and repressor of gene expression, and genes in the nmlr regulon are purported to be involved in resistance to nitrosative stress  <cit> .
though there was an anaerobic increase in nmlr expression, none of the genes in the nmlr regulon, except for nmlr itself, were found to be upregulated.
one member of the regulon, the thioredoxin reductase trxb , was shown to be significantly downregulated .
surprisingly, in a recent study, gonococcal trxb mutants were found to be susceptible to killing by no, though this property was correlated with low expression of ania and norb  <cit> .
it stands to reason, therefore, that genes in the nmlr regulon may be important in adaptation to anaerobic conditions or in certain environments where reactive oxygen or nitrogen species are present, but they are possibly not as important once the anaerobic respiratory chain becomes highly induced.
in this circumstance, norb is likely able to reduce no levels below a threshold needed to induce the activation function of nmlr   <cit> .
under anaerobic conditions, and in the absence of oxygen or reactive oxygen species , no may not need to be detoxified, as toxic no-derived reactive nitrogen compounds will not be produced.
in this instance no would preferentially be used as an energy source.
the gene encoding lexa was found to be upregulated anaerobically .
this transcriptional regulator was proposed to be important in regulating genes involved in defense against ros  <cit> .
lexa controls a small regulon in gonococci that includes itself, dna repair protein, recn, and a gene of unknown function, ngo <dig>  lexa has been demonstrated to sense hydrogen peroxide through thiol modification of a single cysteine residue  <cit> , though the anaerobic induction of the lexa regulon observed here suggests that lexa may also sense no.
interestingly, analysis of the raw rna-seq data and rt-pcr showed that under conditions where the lexa regulon was induced, significant transcriptional read-through occurred beyond lexa into the adjacent erpa gene , as well as past the divergently transcribed ngo <dig> gene into the adjacent dnak gene .
both erpa and dnak are members of the rpoh regulon  <cit> , suggesting that there may be crosstalk between the rpoh and lexa regulons.
though transcriptional read-through may account for part of the anaerobic induction of erpa and dnak, our data suggest that rpoh activity itself was likely modulated anaerobically.
although there was no significant increase in the quantity of rpoh transcripts at the  <dig> hr time point at which the cells were harvested, control of rpoh at the transcriptional level plays a relatively minor role in maintaining rpoh levels in the cell  <cit> .
a majority of rpoh regulation is posttranscriptional.
translation of rpoh is inhibited by secondary structure in the rpoh mrna, and activity of rpoh can be further modulated by the chaperone dnak  <cit> .
including erpa and dnak, eight out of the twelve genes previously demonstrated to encompass the gonococcal rpoh regulon were found to be significantly upregulated anaerobically, including genes encoding a transcriptional regulator of unknown function, marr , dihydrofolate reductase, fola , putative dna-binding protein, crea , and chaperone proteins clpb , grpe , and dnaj   <cit> .
the anaerobic induction of genes in the rpoh regulon may be of significance in vivo, as these genes have been shown to contribute to epithelial cell invasion  <cit> .
the gene encoding farr  was found to be anaerobically induced .
farr was originally described for its role in repression of the farab operon, which is involved in fatty acid resistance  <cit> .
analysis of the rna-seq data demonstrated that farab expression was not altered anaerobically, suggesting that this operon may already be fully repressed under standard laboratory conditions.
a recent study established that farr, possibly indirectly, was also responsible for the repression of two abc transporter cassettes  as well as the gene encoding the multiple transferable resistance protein, mtrf   <cit> .
all of these genes were found to be anaerobically repressed by rna-seq analysis , suggesting that farr is yet another player in the anaerobic lifestyle of n. gonorrhoeae.
interestingly, expression of farr itself was determined to be controlled by repression through mtrr, a regulator of the mtrcde efflux pump involved in resistance to antimicrobial agents  <cit> .
mtrr was recently shown to act in a global fashion, regulating many genes, including rpoh  <cit> .
many of the genes in the reported mtrr regulon also appear to be part of the anaerobic stimulon in gonococci, suggesting a potential role for this regulator in anaerobic growth as well  <cit>  .
genes involved in adaptation and stress response
in addition to the rpoh regulon, many genes involved in response and adaptation to stress were found to be anaerobically induced.
several genes encoding proteins involved in dna repair were found to be upregulated anaerobically .
recn is involved in recombinational repair of dna damage, and has been shown to protect gonococci from both oxidative and nonoxidative mechanisms of killing  <cit> .
recx  has been demonstrated to be involved in enhancement of reca-mediated processes, including repair of dna damage  <cit> .
a gene of unknown function, ngo <dig>  was also upregulated anaerobically, and encodes a deddh-family protein with a high level of homology to dnaq , the epsilon subunit of dna polymerase iii  <cit> .
as the epsilon subunit of dna polymerase iii contains 3' → 5' exonuclease activity and is involved in dna mismatch repair, this dnaq-like protein may also function in the repair of dna damage  <cit> .
several genes typically involved in resistance to oxidative stress were found to be induced anaerobically, including those encoding glutaredoxin, grx <dig>  cytochrome c peroxidase, ccp , tellurite resistance protein, tehb , and nitroreductase-like protein, ntra .
ccp was shown to protect gonococci from hydrogen peroxide mediated damage, but not to be involved in protection from reactive nitrogen species  <cit> .
activation of ccp by fnr suggests that gonococci may typically be exposed to ros under anaerobic or microaerobic conditions in vivo.
proteins similar to tehb were originally described for their role in resistance to tellurite, however, gonococci do not inhabit an environment where tellurite would be encountered.
in haemophilus influenzae, it was recently shown that a deletion of tehb resulted in increased sensitivity to oxidizing agents including hydrogen peroxide, suggesting a dual role for this protein  <cit> .
the protein encoded by ngo <dig> was shown to contain a nitroreductase-like domain  <cit> , and to contain 54% identity  to a novel nitroreductase in staphylococcus aureus, ntra.
in s. aureus, ntra was shown to contribute to nitrosative stress resistance through its s-nitrosoglutathione  reductase activity  <cit> .
two other proteins, dnrn and erpa, were also demonstrated to be upregulated anaerobically.
dnrn has been implicated in the repair of iron sulfur clusters damaged by oxidative or nitrosative stress  <cit> , and erpa may be important in maintaining the iron-sulfur cluster status of specific proteins in cells growing under stressful conditions  <cit> .
interestingly, genes encoding two nudix-family hydrolases, nudh  and ngo <dig>  as well as a histidine triad family protein, hit , were all found to be upregulated anaerobically .
analysis of the conserved domains of these open reading frames reveals the greatest level of homology to proteins involved in metabolism/hydrolysis of diadenosine tetraphosphate   <cit> .
ap4a is an important and ubiquitous signaling molecule in nature, and has been implicated in the maintenance and regulation of vital cellular functions  <cit> .
as a bacterial second messenger, ap4a has been shown to be involved in modulating chaperone and heat shock protein activity, coupling dna replication to cell division, and altering the ability of cells to use alternate carbohydrate sources, though these may be just a few of its functions  <cit> .
if allowed to accumulate unchecked, high ap4a concentration can interfere with a number of atp-dependent reactions  <cit> .
in e. coli, salmonella typhimurium, and bartonella bacilliformis, deletion of the nudh gene resulted in gross deficiencies in the ability of these organisms to invade mammalian cells  <cit> .
the role of ap4a as a second messenger in gonococci has yet to be evaluated thus far, though it appears that the metabolism of this molecule may be modulated anaerobically, making it an interesting target for future study.
hypothetical proteins regulated anaerobically
many genes encoding hypothetical proteins were found to be induced anaerobically .
one gene, listed as duf <dig> in table  <dig>  was not found in the fa <dig> genome annotation, rather it was discovered in proximity to the iron sulfur cluster assembly operon  using raw rna-seq data.
using this data, it was discovered that the isc operon spanned the iscr gene  through hesb , making it longer than the current annotation would suggest.
ngo <dig> probably does not encode a protein, and no reads were mapped on the positive strand at this location .
a high level of transcription was found at an unannotated putative open reading frame with an atg start codon on the positive strand at chromosomal coordinate  <dig>  .
there was a high match to consensus - <dig>  - <dig>  and rbs sequences upstream of this open reading frame, and transcription of this gene partially overlapped the 3' end of the isc operon.
this anaerobically induced gene encodes a protein with a conserved domain of unknown function .
this is just an example of how rna-seq data can be used to locate and correct annotation errors in the gonococcal genome.
the large quantity of hypothetical proteins that are differentially expressed in response to anaerobiosis may make appealing targets for future study.
many of these annotated proteins are small and neisseria-specific.
future work should determine if these are akin to small proteins in e. coli that have been shown to accumulate in response to stress  <cit> .
alternatively, these small open reading frames could be mis-annotations, and may be involved in transcription of regulatory rna.
current work involves integrating raw rna-seq data in a genome-wide search for regulatory rna.
n. gonorrhoeae encodes an anaerobically induced srna
in e. coli and several other enterobacterial species, recent studies have described the anaerobic induction and function of a small regulatory rna termed fnrs  <cit> .
in e. coli, the regulation of this srna was determined to be relatively complex.
maximal fnrs expression occurred anaerobically, however, the available carbon source and terminal electron acceptor present also impacted levels to a lesser extent.
anaerobic induction of fnrs was shown to be fnr-dependent and mediated by fnr binding to a class-ii activation site centered at - <dig>  with respect to the transcription start site  <cit> .
no fnrs homolog has been observed in bacteria outside of the enterobacteriaceae family, however, previous microarray analysis of the gonococcal fnr regulon identified a small fnr-activated transcript of unknown function.
coincidently, this transcript was also determined to be the most highly induced by fnr  <cit> .
this small fnr-induced transcript is located within the coding region and on the opposite strand of ngo <dig>  flanked on the 5' end by the bacterioferritin genes, brfa and bfrb, and on the 3' end by an xre-family repressor  end by an .
though not included the fa <dig> genome annotation, analysis of the raw rna-seq data confirms that this transcript is highly induced anaerobically .
the region upstream of this transcript also contains a perfect match to the e. coli fnr consensus sequence  located at - <dig>  with respect to the transcription start site.
secondary structure prediction of this  <dig> nt transcript was performed by the mfold program http://mfold.rna.albany.edu/, and the most stable predicted structure is shown in figure 5b.
this predicted structure is typical of a small regulatory rna, with a 5' stem loop followed by an unstructured region and a 3' rho-independent terminator .
to confirm that this small rna is expressed and induced anaerobically, a transcriptional fusion to lacz was constructed.
in this construct, the promoter region from the induced transcript  was cloned in front of the leader region of lacz, which provided spacing and an rbs for efficient translation.
aerobic β-galactosidase activity from this transcriptional fusion was low, but was induced 260-fold anaerobically , validating the rna-seq data and the work of whitehead, et al.
<cit> .
to show that transcription was halted at the predicted rho-independent terminator, a transcriptional fusion was generated in which the entire small rna, including the promoter region , was fused to the leader region of lacz.
no β-galactosidase activity was observed from this fusion gene aerobically or anaerobically, confirming the terminator prediction and the raw rna-seq data .
the first three nucleotides of this small transcript are aug. to ensure that this transcript is indeed a small rna and not a leaderless protein-coding message, a translational fusion was constructed.
a leaderless lacz gene was fused in-frame past nucleotide + <dig> in the transcript.
no β-galactosidase activity was observed in this translational fusion under aerobic or anaerobic conditions .
for these reasons, we feel that this small anaerobically induced transcript should be termed fnrs.
much work needs to be completed to define the role, if any, of this gonococcal fnrs transcript in gene regulation, however, homologs of genes in the e. coli fnrs regulon that were found to be anaerobically repressed in gonococci, including gpma  and mqo , will make interesting candidates to study for potential interactions  <cit> .
future work will employ a similar rna-seq approach to define the prospective fnrs regulon in gonococci.
overlap of the anaerobic, iron, and hydrogen peroxide stimulons
in vivo, n. gonorrhoeae will be faced with a multitude of environmental stresses, including oxygen and iron limitation, as well as exposure to reactive oxygen species.
previous studies have employed genome-wide approaches to define the gonococcal iron and hydrogen peroxide responsive stimulons  <cit> .
genes differentially expressed in response to these stimuli were compared to the anaerobic stimulon to search for similarities in gene expression between these environmental signals .
a average of rpkm ratios from replicates in table  <dig>
b data from jackson et al.,   <cit> , taken from time point where iron depletion had the biggest effect on gene expression.
a positive value represents an induction of gene expression under iron depleted conditions, and a negative value represents a repression of gene expression under iron depleted conditions.
c hydrogen peroxide ; data from stohl et al.,   <cit> .
a positive value represents induction of gene expression in the presence of hydrogen peroxide, and a negative value represents a repression in gene expression in the presence of hydrogen peroxide.
† data from ducey, et al.,   <cit> .
it should be noted that there are caveats in directly comparing these data sets due to the differences in experimental design that exist between each study.
the use of microarrays in the iron and hydrogen peroxide studies may be less sensitive than the use of transcriptome sequencing, and therefore it is possible that some additional iron- or hydrogen peroxide-regulated genes may have been missed  <cit> .
differences in the cutoff level used to define differential expression, the growth medium utilized, and the growth phase of the cells used for rna isolation between these studies may also impact the results.
though the iron and hydrogen peroxide studies utilized rna from broth-grown cells in mid-log phase, the growth phase of anaerobically grown gonococcal cells is difficult to determine due to the requirement for growth on solid media, and consequently, this may also affect the comparison of these results with other studies  <cit> .
with these factors in mind, a significant overlap does exist between the members of these stimulons, and it should be noted that of the common members, the response is not always predictable.
that is, genes regulated in response to these stimuli do not follow a common expression pattern between gene sets, underscoring the complexity of gonococcal gene regulation .
for example, the nqr operon and bfra are both members of the iron stimulon and are both repressed under iron-depleted conditions, however, nqr and bfra are repressed and induced, respectively, during anaerobic growth.
from this comparison we hope to discover neisserial genes that may play an important role in survival under a broad range of environmental stress.
analysis of this gene set may also help future work to define transcriptional networks at the regulator level.
gonococcal prophage genes are highly induced anaerobically
somewhat unexpectedly, the rna-seq data showed that a large subset of the anaerobically induced genes were bacteriophage in origin .
the dsdna filamentous prophage, ngoΦ <dig>  is believed to contain all of the coding regions necessary to produce functionally active phage particles.
in a previous study, phage dna from ngoΦ <dig>  as well as phage particles that were likely ngoΦ <dig> -derived, were observed in gonococci, though plaques were not observed  <cit> .
in this study, rna-seq data showed dramatic anaerobic induction of the ngoΦ <dig> coding region.
though no plaques have been observed in anaerobically grown gonococci, this is not the first report of anaerobic prophage induction.
in anaerobically grown pseudomonas aeruginosa, prophage genes were shown to be among the most highly induced, and it was suggested that, more specifically, nitric oxide was likely involved in this apparent induction  <cit> .
evidence has also been provided to show that these anaerobically induced phage are involved in biofilm development and differentiation in p. aeruginosa  <cit> .
current work proposes that gonococcal biofilms are important in pathogenesis.
it was recently shown that anaerobic metabolism occurs in the substratum of gonococcal biofilms, and that no plays an important role in biofilm maintenance  <cit> .
it is interesting to speculate that in the gonococcus, ngoΦ <dig> may play a role similar to that of anaerobically induced phage in p. aeruginosa.
ngoΦ <dig>  ngoΦ <dig>  and ngoΦ <dig> encode incomplete phage genomes and are probably derived from ngoΦ <dig> or ngoΦ <dig>  <cit> .
ngoΦ <dig>  ngoΦ <dig>  and ngoΦ <dig> were all shown to contain coding regions that were upregulated anaerobically .
the functions of these genes are unknown, but the presence of prophage sequences in bacterial genomes could have a profound effect on host cell fitness or pathogenicity.
a diverse group of bacterial virulence factors, including toxins, have been found to be encoded by genes of phage origin.
some prophage can also contain regulatory proteins that affect expression of genes not encoded by phage .
it is certainly possible that none of the anaerobically upregulated gonococcal prophage genes are involved in making active phage particles; alternatively, these genes may have acquired a novel function.
regardless, further research will be required to determine the role of these genes as they relate to anaerobiosis.
expression from the cryptic plasmid is induced anaerobically
the gonococcal cryptic plasmid, pjd <dig>  is  <dig>  bp long and is found in approximately 96% of gonococcal strains.
despite its predominance, the function and replication mechanism of this plasmid remains largely unknown  <cit> .
this plasmid consists of two divergent transcripts, each consisting of five open reading frames .
surprisingly, pjd <dig> transcripts were greatly increased in prevalence anaerobically .
it appears as though the increase in transcript quantity was not due to an increase in anaerobic copy number, as there was no significant difference in the quantity of plasmid recovered from equal numbers of aerobically or anaerobically grown cells .
this result implies that the two divergent cryptic plasmid promoters are subject to transcriptional regulation.
a protein definitions are derived from korch et al.,   <cit>  or the ncbi conserved domain database  <cit> .
b the fold change for each replicate was calculated by comparing total number of reads for each gene in anaerobically grown gonococci to number of reads for each gene in aerobically grown gonococci.
values were normalized to total number of chromosomal reads.
a positive value represents an induction of gene expression anaerobically, and a negative value represents a repression of gene expression anaerobically.
in other plasmids, the essential repa protein has been shown to be involved in regulating gene expression as well as the initiation of plasmid replication   <cit> .
in several organisms, dnak has been shown to facilitate the dissociation of repa dimers into monomers, resulting in derepression  <cit> .
therefore it is possible that the anaerobic increase in dnak concentration could result in a decrease in repa-mediated autorepression .
unfortunately, little is known about the replication cycle of pjd <dig> or the mechanism of repa binding in gonococci, so this potential regulatory mechanism remains only speculative.
an alternate possibility is that an as-of-yet unidentified chromosomally expressed transcription factor or factors acts at the divergent plasmid promoters to regulate their expression.
the effect of increased cryptic plasmid expression on gonococcal physiology is unknown.
previous studies involving pjd <dig> have provided evidence that orf <dig>  orf <dig>  orf <dig>  orf <dig>  cppc, cppb, and cppa are capable of forming protein products, and several of these protein products have no homology to any other proteins in the ncbi conserved domain database   <cit> .
it is uncertain if orf <dig>  orf <dig>  or orf <dig> are translated.
the high frequency and sequence conservation of this plasmid among gonococcal isolates suggests that the plasmid may be important in pathogenesis.
the fact that transcription from this plasmid is induced anaerobically is perhaps suggestive that the gene products of this plasmid are important in an anaerobic niche.
in transcript- <dig>  orf <dig> and orf <dig> encode a putative toxin/antitoxin system analogous to the vapdx systems of haemophilus influenzae and several species of actinobacillus  <cit> .
it is assumed that the vapdx system is involved in stable maintenance of pjd <dig>  <cit> .
interestingly, in h. influenzae, the vapdx locus is located on the chromosome, and vapd has been characterized as a virulence factor.
it has been proposed that the ability of vapd to cause an arrest of bacterial translation was responsible for an observed enhancement of bacterial survival within human cells compared to survival in a Δvapd mutant  <cit> .
vapd was postulated to facilitate persistent or chronic mucosal infections under stressful conditions.
the role of vapd as it pertains to gonococcal virulence, if any, is unknown, though it may make an attractive target for future study.
CONCLUSIONS
the goal of this study was to define the gonococcal transcriptome in response to anaerobiosis.
unlike previous studies that have looked only at the contributions of individual transcriptional regulators, we used a powerful and sensitive methodology, rna-seq  <cit> , to look at the environmental condition of anaerobiosis as a whole.
by using this method, we were able to show that the anaerobic stimulon in n. gonorrhoeae was larger than the fnr and nsrr regulons indentified in previous studies  <cit> .
the induction of prophage and cryptic plasmid genes was previously unknown to occur during anaerobiosis.
many new targets for future study have been uncovered, and we hope that other investigators will take an interest in working on these potential projects.
because anaerobic growth is assumed to be a physiologically significant state in vivo  <cit> , the results of this study may help to elucidate factors or mechanisms of virulence that have previously been overlooked.

 n. gonorrhoeae was long considered to be an obligate aerobe until it was discovered that anaerobic growth was possible when nitrite or nitric oxide was used as a terminal electron acceptor. anaerobic growth is accomplished through utilization of a truncated denitrification pathway, which is encoded within the gonococcal genome as a pair of divergently transcribed genes, ania, encoding a nitrite reductase, and norb, encoding a nitric oxide reductase. anaerobiosis is presumed to be a physiologically significant state during infection, as the gonococcus is often recovered from infected individuals in co-culture with obligate anaerobes such as peptococcus and bacteroides spp.. furthermore, gonococci have been shown to induce and repress the expression of several genes in response to anaerobiosis, and antibody to ania, the major anaerobically induced outer membrane protein, can be found in sera from infected women, demonstrating that this protein is expressed in vivo. the ability of gonococci to utilize this denitrification pathway to reduce no may have immunomodulary effects during the course of infection. some evidence suggests that the gonococcal reduction of host-produced no may be responsible, at least in part, for the high incidence of asymptomatic disease. in vitro, n. gonorrhoeae was shown to be capable of setting a no steady-state level in the anti-inflammatory range. recently, in an in vitro cervical cell model of infection, it was shown that gonococcal activation of inos promoted bacterial survival. in this study, it was suggested that host-derived nitric oxide is not protective against gonococci, rather, nitric oxide may actually be required to sustain cervical bacterial disease. previous research has shown that the genes involved in denitrification and/or adaptation to anaerobic growth in neisseria spp. are subject to transcriptional regulation by the oxygen-sensitive regulator, fnr, the no-sensitive repressor, nsrr, and the nitrite-insensitive two component system, narqp. earlier data from a microarray-based approach suggested that the gonococcal fnr-regulon was composed of fourteen activated and six repressed transcripts, making the gonococcal fnr regulon much smaller than that of e. coli, where fnr was shown to regulate more than operons. discovery of nsrr and narqp-regulated genes in neisseria spp. has relied mainly on in silico analysis of intergenic chromosomal regions to identify similarity to previously defined regulator binding sites. as is the case with fnr, the neisserial nsrr and narqp regulons, at least the currently defined members, are much smaller than their e. coli counterparts. in this study we employ a powerful whole-genome approach, rna-seq, to quantitatively sequence the complete gonococcal transcriptome. using this method we were able to define global changes in gene expression in response to anaerobiosis. we show that the gonococcal anaerobic stimulon is not small, and that chromosomal open reading frames are differentially expressed. we present several novel findings that, taken together, support the view that anaerobic growth is an important facet of life for this organism, and should be considered when studying the host/pathogen interaction.  sequencing the gonococcal transcriptome whole genome mrna sequencing is an attractive method of monitoring global changes in gene expression while overcoming many of the pitfalls of traditional dna microarrays. for the purposes of this study, anaerobic conditions were defined as anoxia, the presence of nitrite and the concomitant presence of nitric oxide. in order to define the gonococcal anaerobic stimulon, rna-seq was performed on two biological replicates of aerobically or anaerobically plate-grown gonococci using the abi solid™ system. unique sequence reads from the rna-seq data were mapped to the annotated fa genome, and gene expression was quantified as reads per kilobase of coding sequence per million reads. when expression data for each replicate were plotted against each other, rpkm values were observed to be adequately reproducible, though expectedly, the extent of reproducibility was slightly less in genes with low expression. for this reason, we required mapped genes to have an rpkm ≥ in order to be considered as a candidate for differential expression. however, genes that had an rpmk < under one growth condition, but were highly induced in the other condition were accepted as candidates. genes were considered to be differentially expressed if there was a three-fold or greater difference in rpkm between the two growth conditions. analysis of rna-seq data revealed that chromosomal genes met the criteria for being differentially expressed, with genes being anaerobically induced and genes being anaerobically repressed. figure 1b and 1c illustrate a reasonable measure of reproducibility in expression level between replicates in the differentially expressed genes. a gonococci grown in an anaerobe chamber on plates containing mm nitrite. b the fold change for each replicate was calculated by comparing rpkm values of anaerobically grown gonococci to rpkm values of aerobically grown gonococci. a positive value represents an induction of gene expression anaerobically, and a negative value represents a repression of gene expression anaerobically c protein definitions are derived from the annotated fa genome, and/or the ncbi conserved domain database. d genes that have been shown to be directly regulated by fur in sebastian et al., and/or jackson et al.,. genes reported to be upregulated during rpoh overexpression in gunesekere et al.,. genes reported to be in the lexa regulon by schook et al.,. genes found to be differentially expressed in a n. meningitidis Δlrp mutant by microarray analysis in ren et al.,. genes found to be differentially expressed in a gonococcal Δmtrr mutant by microarray analysis in folsterΔ et al.,. genes reported to be in the farr regulon by friedrich et al.,. genes reported to be in the nmlr regulon by kidd et al.,. genes found to be differentially expressed microaerobically in a Δfnr mutant with and without nitrite by microarray analysis in whitehead et al.,. genes reported to be in the nsrr regulon by isabella et al., and isabella,. e fur-repressed small regulatory rna described by ducey et al.,. as expected, anaerobic expression levels cluster at higher rpkm values in the set of anaerobically induced genes, while the opposite is true in anaerobically repressed genes. secondary verification of differentially expressed genes many of the genes found to have been differentially expressed by rna-seq have previously been determined to be members of the neisserial fnr or nsrr regulons, validating the use of this technique to monitor transcriptome changes. to further corroborate the rna-seq data, several genes shown to be differentially expressed were selected for secondary confirmation. analysis of translational promoter-lacz fusions or rt-pcr was used to accomplish this task. analysis of translational promoter-lacz fusions demonstrated that the genes encoding heat shock protein, clpb, dna repair enzyme, recn, nitropropane dioxygenase, npd, filamentous hemagglutinin, heca, glutaredoxin, grx, bacterioferritin, bfra and repressor protein lexa were significantly upregulated anaerobically. rt-pcr confirmed that 2-isopropylmalate synthase, leu, and adhesin, leca, were also increased in expression anaerobically. translational promoter-lacz fusions to the genes encoding the iron-binding protein, fbpa, and the solute binding protein, oxia confirmed that these genes are anaerobically repressed. rt-pcr confirmed that ngo l-lactate dehydrogenase, lldd, thioredoxin, trx, heme oxygenase, hemo, and nadh dehydrogenase component, nqra were also anaerobically repressed. analysis of raw rna-seq data allows for the prediction of transcription start sites, - and - elements, as has been described in rna-seq projects in both helicobacter pylori and listeria monocytogenes. compared to e. coli, the use of alternative sigma factors in the gonococcus is very limited. gonococci have no sigma- homolog, and only utilize sigma factors in the sigma- family. this class of sigma factors recognize - and - elements within bacterial promoters. the predictions of transcriptional elements for the individual genes selected for secondary analysis are shown in figure these predictions were generally a high match to the e. coli consensus for transcriptional elements or the α-proteobacteria consensus for rpoh-dependent transcriptional elements. the ability to predict the location of these transcriptional elements in a large scale manner will prove very useful in the search for potential regulatory sites, and will aid in future work to define gonococcal transcriptional networks. differentially expressed genes involved in macromolecular biosynthesis the genes found to be differentially expressed in response to anaerobic growth were broadly characterized according to their putative function. several genes encoding proteins involved in macromolecular synthesis displayed decreased expression under anaerobic conditions. the synthesis of ribosomal proteins is strongly related to growth rate, and ribosomal protein synthesis has been documented to decrease during times of energy deficiency. for facultative anaerobes, growth in an oxygen-limited environment is more energetically-deprived than in aerobic conditions. the reduced synthesis of ribosomal proteins is indicative of adaptation to a slower growth rate following recovery from a nutritional shift-down. the reduced expression of rna polymerase β an β' subunits is also indicative of this, and has been observed in e. coli as well. genes involved in transport and binding many genes encoding proteins involved in transport and binding were found to be anaerobically repressed. most of these proteins are involved in the transport of amino acids, sugars, or other growth metabolites into the cell. this result comes as no surprise, as a slower anaerobic growth rate would result in decreased demand and slower utilization of such substrates. proteins of this type were also shown to be downregulated anaerobically in e. coli, p. aeruginosa, and s. enterica. three repressed genes, fbpa, hemo, and ompu, are involved in the acquisition and/or transport of iron into the cell. in s. enterica, expression of the gene encoding iron acquisition protein, sita, was also found to be repressed anaerobically. while genes involved in iron acquisition and transport are repressed anaerobically, the genes encoding the iron-storage bacterioferritin proteins, bfra and bfrb are induced. three other proteins implicated in binding that were shown to be anaerobically induced, leca, heca, and fhbp, may be involved directly in the host-pathogen interaction. the role of leca in gonococci is unclear, but a leca homolog in another member of the neisseriaceae family, eikenella corrodens, has been shown to be a part of an adhesin complex important in adherence to a variety of host tissues. the gonococcal heca protein contains a hemagglutination domain that, in other bacteria, has been shown to promote bacterial aggregation or host attachment. the role of factor h binding protein as it relates to gonococcal pathogenesis has not been examined, however, in n. meningitidis, fhbp expression on the cell surface was shown to be a means of immune evasion that rendered meningococci resistant to complement-mediated killing. genes involved in energy production/oxidation-reduction unsurprisingly, numerous genes encoding proteins involved in energy metabolism were found to be anaerobically repressed. gonococci use a combination of the entner-douderoff and pentose phosphate pathways in the utilization of glucose. key enzymes involved in these pathways were found to be anaerobically repressed, as well as several enzymatic components of the krebs cycle. anaerobic repression of genes involved in central intermediary metabolism has also been observed in other bacteria, including e. coli, p. aeruginosa, and b. subtilus. interestingly, gonococci contain two major nadh dehydrogenases capable of coupling the oxidation of nadh to the reduction of ubiquinone. one of these dehydrogenases is the prototypical h+-translocating 14-subunit complex encoded by the nuo operon. the second is a seven-subunit na+-translocating complex encoded by the nqr operon. though genes in the nuo operon did not meet the criteria used in this study to define differential expression, there was an approximate two-fold level of anaerobic repression across its entire length. the nqr operon, however, demonstrated a significant decrease in anaerobic expression. na+-translocating nadh dehydrogenases were originally discovered and characterized in marine bacteria, though it has now been revealed that several pathogenic species encode such complexes. transport of sodium ions into the periplasmic space through the nqr complex generates an electrochemical membrane potential that plays a role in solute import, atp synthesis, and, for organisms that are motile, flagellar rotation. the flow of electrons through the active portions of the electron transport chain would be expected to decrease anaerobically when: oxygen is unavailable, rendering the terminal cytochrome oxidase inactive, enzymatic components of the krebs cycle are repressed, resulting in decreased production of reducing equivalents, and growth is dependent on high potential alternative electron acceptors such as no2- and no. furthermore, neither the nitrite or nitric oxide reductases are thought to be capable of pumping protons across the membrane. downregulation of nqr would decrease the electrochemical potential across the membrane, thus decreasing the proton motive force. because expression of the nuo operon is not affected to the same degree anaerobically, this would now result in a greater proportion of nadh oxidation being coupled to atp production through h+ translocation. also, many proteins involved in the transport of metabolites into the cell are repressed anaerobically. downregulation of the nqr complex may add another level of regulation to cellular transport. it stands to reason that the activity of solute transporters that rely on the sodium motive force would decrease when fewer na+ ions are pumped across the membrane. the function and regulation of the most highly induced gene involved in anaerobic growth, the nitrite reductase ania, has been thoroughly characterized, and has been shown to be required for anaerobic growth on nitrite. interestingly, gonococci contain a gene encoding a 2-nitropropane dioxygenase-like protein that was found to be induced anaerobically. this would seem unusual, as gonococci would not be expected to encounter the toxic nitroalkanes utilized by npd proteins. furthermore, molecular oxygen is consumed in the npd reaction mechanism. further work should determine if this protein plays a novel role in gonococci, possibly microaerobically when oxygen would be present. gonococci also showed increased anaerobic expression of a gene encoding electron transfer flavoprotein-ubiquinone oxidoreductase, eftd, a protein containing 57% identity to human eftd. the presence of eftd in bacteria is rare, and its role has only been studied in plants and mammals. in arabidopsis thaliana, eftd was shown to be upregulated in times of carbon starvation, when genes involved in glycolysis and the krebs cycle were reduced in expression, as is observed when gonococci are grown anaerobically. in a. thaliana, it was reported that eftd can act as an electron acceptor from various dehydrogenases involved in amino acid catabolism and fatty acid degradation. the reduction of ubiquinone by etfd was proposed to be a mechanism of utilizing alternative respiratory substrates to feed into the electron transport chain. eftd upregulation in n. gonorrhoeae could indicate a shift to the use of alternative growth substrates under anaerobic conditions. genes involved in small molecule biosynthesis anaerobically, several genes involved in the synthesis of amino acids were found to be repressed, which has also been observed in other bacteria. genes involved in the synthesis of iron-sulfur clusters were also repressed. this result is not entirely unexpected, as many iron-sulfur cluster containing proteins involved in respiration were shown to be downregulated anaerobically. the demand for iron-sulfur clusters may therefore be reduced under anaerobic conditions. also, because it is common for the fe-s clusters in fe-s cluster-containing proteins to be oxygen-sensitive, a lower rate of oxygen-induced fe-s cluster turnover would be expected under anaerobic conditions. unlike other amino acid biosynthetic genes, the gene encoding leu the first committed step in leucine biosynthesis, was found to be induced anaerobically. this could be indicative of a decrease in intracellular leucine concentration, as leu is repressed by leucine in e. coli and other organisms. in e. coli, the anaerobic stimulon is intimately associated with global regulation through the leucine responsive regulatory protein, lrp.. in n. meningitidis, the reported lrp regulon demonstrated striking overlap to many of the genes in the gonococcal anaerobic stimulon, particularly in regards to proteins involved in glucose utilization, the krebs cycle, and ribosomal assembly. future work will examine whether gonococcal lrp plays a role in control of the anaerobic stimulon. genes involved in transcriptional regulation the expression of several regulatory proteins was found to be induced under anaerobic conditions, though no regulatory proteins appeared to be repressed. the gene encoding arsr, a regulator involved in norb repression in response to iron, was shown to be moderately upregulated. it is unclear if arsr regulates any other target genes. the fnr-regulated nmlr protein has been shown to act as both an activator and repressor of gene expression, and genes in the nmlr regulon are purported to be involved in resistance to nitrosative stress. though there was an anaerobic increase in nmlr expression, none of the genes in the nmlr regulon, except for nmlr itself, were found to be upregulated. one member of the regulon, the thioredoxin reductase trxb, was shown to be significantly downregulated. surprisingly, in a recent study, gonococcal trxb mutants were found to be susceptible to killing by no, though this property was correlated with low expression of ania and norb. it stands to reason, therefore, that genes in the nmlr regulon may be important in adaptation to anaerobic conditions or in certain environments where reactive oxygen or nitrogen species are present, but they are possibly not as important once the anaerobic respiratory chain becomes highly induced. in this circumstance, norb is likely able to reduce no levels below a threshold needed to induce the activation function of nmlr. under anaerobic conditions, and in the absence of oxygen or reactive oxygen species, no may not need to be detoxified, as toxic no-derived reactive nitrogen compounds will not be produced. in this instance no would preferentially be used as an energy source. the gene encoding lexa was found to be upregulated anaerobically. this transcriptional regulator was proposed to be important in regulating genes involved in defense against ros. lexa controls a small regulon in gonococci that includes itself, dna repair protein, recn, and a gene of unknown function, ngo lexa has been demonstrated to sense hydrogen peroxide through thiol modification of a single cysteine residue, though the anaerobic induction of the lexa regulon observed here suggests that lexa may also sense no. interestingly, analysis of the raw rna-seq data and rt-pcr showed that under conditions where the lexa regulon was induced, significant transcriptional read-through occurred beyond lexa into the adjacent erpa gene, as well as past the divergently transcribed ngo gene into the adjacent dnak gene. both erpa and dnak are members of the rpoh regulon, suggesting that there may be crosstalk between the rpoh and lexa regulons. though transcriptional read-through may account for part of the anaerobic induction of erpa and dnak, our data suggest that rpoh activity itself was likely modulated anaerobically. although there was no significant increase in the quantity of rpoh transcripts at the hr time point at which the cells were harvested, control of rpoh at the transcriptional level plays a relatively minor role in maintaining rpoh levels in the cell. a majority of rpoh regulation is posttranscriptional. translation of rpoh is inhibited by secondary structure in the rpoh mrna, and activity of rpoh can be further modulated by the chaperone dnak. including erpa and dnak, eight out of the twelve genes previously demonstrated to encompass the gonococcal rpoh regulon were found to be significantly upregulated anaerobically, including genes encoding a transcriptional regulator of unknown function, marr, dihydrofolate reductase, fola, putative dna-binding protein, crea, and chaperone proteins clpb, grpe, and dnaj. the anaerobic induction of genes in the rpoh regulon may be of significance in vivo, as these genes have been shown to contribute to epithelial cell invasion. the gene encoding farr was found to be anaerobically induced. farr was originally described for its role in repression of the farab operon, which is involved in fatty acid resistance. analysis of the rna-seq data demonstrated that farab expression was not altered anaerobically, suggesting that this operon may already be fully repressed under standard laboratory conditions. a recent study established that farr, possibly indirectly, was also responsible for the repression of two abc transporter cassettes as well as the gene encoding the multiple transferable resistance protein, mtrf. all of these genes were found to be anaerobically repressed by rna-seq analysis, suggesting that farr is yet another player in the anaerobic lifestyle of n. gonorrhoeae. interestingly, expression of farr itself was determined to be controlled by repression through mtrr, a regulator of the mtrcde efflux pump involved in resistance to antimicrobial agents. mtrr was recently shown to act in a global fashion, regulating many genes, including rpoh. many of the genes in the reported mtrr regulon also appear to be part of the anaerobic stimulon in gonococci, suggesting a potential role for this regulator in anaerobic growth as well. genes involved in adaptation and stress response in addition to the rpoh regulon, many genes involved in response and adaptation to stress were found to be anaerobically induced. several genes encoding proteins involved in dna repair were found to be upregulated anaerobically. recn is involved in recombinational repair of dna damage, and has been shown to protect gonococci from both oxidative and nonoxidative mechanisms of killing. recx has been demonstrated to be involved in enhancement of reca-mediated processes, including repair of dna damage. a gene of unknown function, ngo was also upregulated anaerobically, and encodes a deddh-family protein with a high level of homology to dnaq, the epsilon subunit of dna polymerase iii. as the epsilon subunit of dna polymerase iii contains 3' → 5' exonuclease activity and is involved in dna mismatch repair, this dnaq-like protein may also function in the repair of dna damage. several genes typically involved in resistance to oxidative stress were found to be induced anaerobically, including those encoding glutaredoxin, grx cytochrome c peroxidase, ccp, tellurite resistance protein, tehb, and nitroreductase-like protein, ntra. ccp was shown to protect gonococci from hydrogen peroxide mediated damage, but not to be involved in protection from reactive nitrogen species. activation of ccp by fnr suggests that gonococci may typically be exposed to ros under anaerobic or microaerobic conditions in vivo. proteins similar to tehb were originally described for their role in resistance to tellurite, however, gonococci do not inhabit an environment where tellurite would be encountered. in haemophilus influenzae, it was recently shown that a deletion of tehb resulted in increased sensitivity to oxidizing agents including hydrogen peroxide, suggesting a dual role for this protein. the protein encoded by ngo was shown to contain a nitroreductase-like domain, and to contain 54% identity to a novel nitroreductase in staphylococcus aureus, ntra. in s. aureus, ntra was shown to contribute to nitrosative stress resistance through its s-nitrosoglutathione reductase activity. two other proteins, dnrn and erpa, were also demonstrated to be upregulated anaerobically. dnrn has been implicated in the repair of iron sulfur clusters damaged by oxidative or nitrosative stress, and erpa may be important in maintaining the iron-sulfur cluster status of specific proteins in cells growing under stressful conditions. interestingly, genes encoding two nudix-family hydrolases, nudh and ngo as well as a histidine triad family protein, hit, were all found to be upregulated anaerobically. analysis of the conserved domains of these open reading frames reveals the greatest level of homology to proteins involved in metabolism/hydrolysis of diadenosine tetraphosphate. ap4a is an important and ubiquitous signaling molecule in nature, and has been implicated in the maintenance and regulation of vital cellular functions. as a bacterial second messenger, ap4a has been shown to be involved in modulating chaperone and heat shock protein activity, coupling dna replication to cell division, and altering the ability of cells to use alternate carbohydrate sources, though these may be just a few of its functions. if allowed to accumulate unchecked, high ap4a concentration can interfere with a number of atp-dependent reactions. in e. coli, salmonella typhimurium, and bartonella bacilliformis, deletion of the nudh gene resulted in gross deficiencies in the ability of these organisms to invade mammalian cells. the role of ap4a as a second messenger in gonococci has yet to be evaluated thus far, though it appears that the metabolism of this molecule may be modulated anaerobically, making it an interesting target for future study. hypothetical proteins regulated anaerobically many genes encoding hypothetical proteins were found to be induced anaerobically. one gene, listed as duf in table was not found in the fa genome annotation, rather it was discovered in proximity to the iron sulfur cluster assembly operon using raw rna-seq data. using this data, it was discovered that the isc operon spanned the iscr gene through hesb, making it longer than the current annotation would suggest. ngo probably does not encode a protein, and no reads were mapped on the positive strand at this location. a high level of transcription was found at an unannotated putative open reading frame with an atg start codon on the positive strand at chromosomal coordinate. there was a high match to consensus - - and rbs sequences upstream of this open reading frame, and transcription of this gene partially overlapped the 3' end of the isc operon. this anaerobically induced gene encodes a protein with a conserved domain of unknown function. this is just an example of how rna-seq data can be used to locate and correct annotation errors in the gonococcal genome. the large quantity of hypothetical proteins that are differentially expressed in response to anaerobiosis may make appealing targets for future study. many of these annotated proteins are small and neisseria-specific. future work should determine if these are akin to small proteins in e. coli that have been shown to accumulate in response to stress. alternatively, these small open reading frames could be mis-annotations, and may be involved in transcription of regulatory rna. current work involves integrating raw rna-seq data in a genome-wide search for regulatory rna. n. gonorrhoeae encodes an anaerobically induced srna in e. coli and several other enterobacterial species, recent studies have described the anaerobic induction and function of a small regulatory rna termed fnrs. in e. coli, the regulation of this srna was determined to be relatively complex. maximal fnrs expression occurred anaerobically, however, the available carbon source and terminal electron acceptor present also impacted levels to a lesser extent. anaerobic induction of fnrs was shown to be fnr-dependent and mediated by fnr binding to a class-ii activation site centered at - with respect to the transcription start site. no fnrs homolog has been observed in bacteria outside of the enterobacteriaceae family, however, previous microarray analysis of the gonococcal fnr regulon identified a small fnr-activated transcript of unknown function. coincidently, this transcript was also determined to be the most highly induced by fnr. this small fnr-induced transcript is located within the coding region and on the opposite strand of ngo flanked on the 5' end by the bacterioferritin genes, brfa and bfrb, and on the 3' end by an xre-family repressor end by an. though not included the fa genome annotation, analysis of the raw rna-seq data confirms that this transcript is highly induced anaerobically. the region upstream of this transcript also contains a perfect match to the e. coli fnr consensus sequence located at - with respect to the transcription start site. secondary structure prediction of this nt transcript was performed by the mfold program http://mfold.rna.albany.edu/, and the most stable predicted structure is shown in figure 5b. this predicted structure is typical of a small regulatory rna, with a 5' stem loop followed by an unstructured region and a 3' rho-independent terminator. to confirm that this small rna is expressed and induced anaerobically, a transcriptional fusion to lacz was constructed. in this construct, the promoter region from the induced transcript was cloned in front of the leader region of lacz, which provided spacing and an rbs for efficient translation. aerobic β-galactosidase activity from this transcriptional fusion was low, but was induced 260-fold anaerobically, validating the rna-seq data and the work of whitehead, et al.. to show that transcription was halted at the predicted rho-independent terminator, a transcriptional fusion was generated in which the entire small rna, including the promoter region, was fused to the leader region of lacz. no β-galactosidase activity was observed from this fusion gene aerobically or anaerobically, confirming the terminator prediction and the raw rna-seq data. the first three nucleotides of this small transcript are aug. to ensure that this transcript is indeed a small rna and not a leaderless protein-coding message, a translational fusion was constructed. a leaderless lacz gene was fused in-frame past nucleotide + in the transcript. no β-galactosidase activity was observed in this translational fusion under aerobic or anaerobic conditions. for these reasons, we feel that this small anaerobically induced transcript should be termed fnrs. much work needs to be completed to define the role, if any, of this gonococcal fnrs transcript in gene regulation, however, homologs of genes in the e. coli fnrs regulon that were found to be anaerobically repressed in gonococci, including gpma and mqo, will make interesting candidates to study for potential interactions. future work will employ a similar rna-seq approach to define the prospective fnrs regulon in gonococci. overlap of the anaerobic, iron, and hydrogen peroxide stimulons in vivo, n. gonorrhoeae will be faced with a multitude of environmental stresses, including oxygen and iron limitation, as well as exposure to reactive oxygen species. previous studies have employed genome-wide approaches to define the gonococcal iron and hydrogen peroxide responsive stimulons. genes differentially expressed in response to these stimuli were compared to the anaerobic stimulon to search for similarities in gene expression between these environmental signals. a average of rpkm ratios from replicates in table b data from jackson et al.,, taken from time point where iron depletion had the biggest effect on gene expression. a positive value represents an induction of gene expression under iron depleted conditions, and a negative value represents a repression of gene expression under iron depleted conditions. c hydrogen peroxide; data from stohl et al.,. a positive value represents induction of gene expression in the presence of hydrogen peroxide, and a negative value represents a repression in gene expression in the presence of hydrogen peroxide. † data from ducey, et al.,. it should be noted that there are caveats in directly comparing these data sets due to the differences in experimental design that exist between each study. the use of microarrays in the iron and hydrogen peroxide studies may be less sensitive than the use of transcriptome sequencing, and therefore it is possible that some additional iron- or hydrogen peroxide-regulated genes may have been missed. differences in the cutoff level used to define differential expression, the growth medium utilized, and the growth phase of the cells used for rna isolation between these studies may also impact the results. though the iron and hydrogen peroxide studies utilized rna from broth-grown cells in mid-log phase, the growth phase of anaerobically grown gonococcal cells is difficult to determine due to the requirement for growth on solid media, and consequently, this may also affect the comparison of these results with other studies. with these factors in mind, a significant overlap does exist between the members of these stimulons, and it should be noted that of the common members, the response is not always predictable. that is, genes regulated in response to these stimuli do not follow a common expression pattern between gene sets, underscoring the complexity of gonococcal gene regulation. for example, the nqr operon and bfra are both members of the iron stimulon and are both repressed under iron-depleted conditions, however, nqr and bfra are repressed and induced, respectively, during anaerobic growth. from this comparison we hope to discover neisserial genes that may play an important role in survival under a broad range of environmental stress. analysis of this gene set may also help future work to define transcriptional networks at the regulator level. gonococcal prophage genes are highly induced anaerobically somewhat unexpectedly, the rna-seq data showed that a large subset of the anaerobically induced genes were bacteriophage in origin. the dsdna filamentous prophage, ngoΦ is believed to contain all of the coding regions necessary to produce functionally active phage particles. in a previous study, phage dna from ngoΦ as well as phage particles that were likely ngoΦ -derived, were observed in gonococci, though plaques were not observed. in this study, rna-seq data showed dramatic anaerobic induction of the ngoΦ coding region. though no plaques have been observed in anaerobically grown gonococci, this is not the first report of anaerobic prophage induction. in anaerobically grown pseudomonas aeruginosa, prophage genes were shown to be among the most highly induced, and it was suggested that, more specifically, nitric oxide was likely involved in this apparent induction. evidence has also been provided to show that these anaerobically induced phage are involved in biofilm development and differentiation in p. aeruginosa. current work proposes that gonococcal biofilms are important in pathogenesis. it was recently shown that anaerobic metabolism occurs in the substratum of gonococcal biofilms, and that no plays an important role in biofilm maintenance. it is interesting to speculate that in the gonococcus, ngoΦ may play a role similar to that of anaerobically induced phage in p. aeruginosa. ngoΦ ngoΦ and ngoΦ encode incomplete phage genomes and are probably derived from ngoΦ or ngoΦ. ngoΦ ngoΦ and ngoΦ were all shown to contain coding regions that were upregulated anaerobically. the functions of these genes are unknown, but the presence of prophage sequences in bacterial genomes could have a profound effect on host cell fitness or pathogenicity. a diverse group of bacterial virulence factors, including toxins, have been found to be encoded by genes of phage origin. some prophage can also contain regulatory proteins that affect expression of genes not encoded by phage. it is certainly possible that none of the anaerobically upregulated gonococcal prophage genes are involved in making active phage particles; alternatively, these genes may have acquired a novel function. regardless, further research will be required to determine the role of these genes as they relate to anaerobiosis. expression from the cryptic plasmid is induced anaerobically the gonococcal cryptic plasmid, pjd is bp long and is found in approximately 96% of gonococcal strains. despite its predominance, the function and replication mechanism of this plasmid remains largely unknown. this plasmid consists of two divergent transcripts, each consisting of five open reading frames. surprisingly, pjd transcripts were greatly increased in prevalence anaerobically. it appears as though the increase in transcript quantity was not due to an increase in anaerobic copy number, as there was no significant difference in the quantity of plasmid recovered from equal numbers of aerobically or anaerobically grown cells. this result implies that the two divergent cryptic plasmid promoters are subject to transcriptional regulation. a protein definitions are derived from korch et al., or the ncbi conserved domain database. b the fold change for each replicate was calculated by comparing total number of reads for each gene in anaerobically grown gonococci to number of reads for each gene in aerobically grown gonococci. values were normalized to total number of chromosomal reads. a positive value represents an induction of gene expression anaerobically, and a negative value represents a repression of gene expression anaerobically. in other plasmids, the essential repa protein has been shown to be involved in regulating gene expression as well as the initiation of plasmid replication. in several organisms, dnak has been shown to facilitate the dissociation of repa dimers into monomers, resulting in derepression. therefore it is possible that the anaerobic increase in dnak concentration could result in a decrease in repa-mediated autorepression. unfortunately, little is known about the replication cycle of pjd or the mechanism of repa binding in gonococci, so this potential regulatory mechanism remains only speculative. an alternate possibility is that an as-of-yet unidentified chromosomally expressed transcription factor or factors acts at the divergent plasmid promoters to regulate their expression. the effect of increased cryptic plasmid expression on gonococcal physiology is unknown. previous studies involving pjd have provided evidence that orf orf orf orf cppc, cppb, and cppa are capable of forming protein products, and several of these protein products have no homology to any other proteins in the ncbi conserved domain database. it is uncertain if orf orf or orf are translated. the high frequency and sequence conservation of this plasmid among gonococcal isolates suggests that the plasmid may be important in pathogenesis. the fact that transcription from this plasmid is induced anaerobically is perhaps suggestive that the gene products of this plasmid are important in an anaerobic niche. in transcript- orf and orf encode a putative toxin/antitoxin system analogous to the vapdx systems of haemophilus influenzae and several species of actinobacillus. it is assumed that the vapdx system is involved in stable maintenance of pjd. interestingly, in h. influenzae, the vapdx locus is located on the chromosome, and vapd has been characterized as a virulence factor. it has been proposed that the ability of vapd to cause an arrest of bacterial translation was responsible for an observed enhancement of bacterial survival within human cells compared to survival in a Δvapd mutant. vapd was postulated to facilitate persistent or chronic mucosal infections under stressful conditions. the role of vapd as it pertains to gonococcal virulence, if any, is unknown, though it may make an attractive target for future study.  the goal of this study was to define the gonococcal transcriptome in response to anaerobiosis. unlike previous studies that have looked only at the contributions of individual transcriptional regulators, we used a powerful and sensitive methodology, rna-seq, to look at the environmental condition of anaerobiosis as a whole. by using this method, we were able to show that the anaerobic stimulon in n. gonorrhoeae was larger than the fnr and nsrr regulons indentified in previous studies. the induction of prophage and cryptic plasmid genes was previously unknown to occur during anaerobiosis. many new targets for future study have been uncovered, and we hope that other investigators will take an interest in working on these potential projects. because anaerobic growth is assumed to be a physiologically significant state in vivo, the results of this study may help to elucidate factors or mechanisms of virulence that have previously been overlooked. 
we also observed a large induction of genes encoded within the cryptic plasmid, pjd validation of rna-seq data using translational-lacz fusions or rt-pcr demonstrated the rna-seq results to be very reproducible. we also confirmed expression and regulation of a small rna, likely a functional equivalent of fnrs in the enterobacteriaceae family. here we determined that the anaerobic stimulon in gonococci was larger than previous studies would suggest. we determined that chromosomal genes were differentially expressed in response to anaerobic conditions. gonococci will be subject to many forms of environmental stress, including oxygen-limitation, during the course of infection. many new targets for future research have been uncovered, and the results derived from this study may have helped to elucidate factors or mechanisms of virulence that may have otherwise been overlooked. maintenance of an anaerobic denitrification system in the obligate human pathogen, neisseria gonorrhoeae, suggests that an anaerobic lifestyle may be important during the course of infection. using the information derived from this sequencing, we discuss the implications of the robust transcriptional response to anaerobic growth. surprisingly, many genes of prophage origin were induced anaerobically, as well as several transcriptional regulators previously unknown to be involved in anaerobic growth. furthermore, mounting evidence suggests that reduction of host-produced nitric oxide has several immunomodulary effects on the host.
 maintenance of an anaerobic denitrification system in the obligate human pathogen, neisseria gonorrhoeae, suggests that an anaerobic lifestyle may be important during the course of infection. furthermore, mounting evidence suggests that reduction of host-produced nitric oxide has several immunomodulary effects on the host. however, at this point there have been no studies analyzing the complete gonococcal transcriptome response to anaerobiosis. here we performed deep sequencing to compare the gonococcal transcriptomes of aerobically and anaerobically grown cells. using the information derived from this sequencing, we discuss the implications of the robust transcriptional response to anaerobic growth.  we determined that chromosomal genes were differentially expressed in response to anaerobic conditions. we also observed a large induction of genes encoded within the cryptic plasmid, pjd validation of rna-seq data using translational-lacz fusions or rt-pcr demonstrated the rna-seq results to be very reproducible. surprisingly, many genes of prophage origin were induced anaerobically, as well as several transcriptional regulators previously unknown to be involved in anaerobic growth. we also confirmed expression and regulation of a small rna, likely a functional equivalent of fnrs in the enterobacteriaceae family. we also determined that many genes found to be responsive to anaerobiosis have also been shown to be responsive to iron and/or oxidative stress.  gonococci will be subject to many forms of environmental stress, including oxygen-limitation, during the course of infection. here we determined that the anaerobic stimulon in gonococci was larger than previous studies would suggest. many new targets for future research have been uncovered, and the results derived from this study may have helped to elucidate factors or mechanisms of virulence that may have otherwise been overlooked. 