Previous Article | Next Article ![]()
Journal of Bacteriology, May 2006, p. 3463-3469, Vol. 188, No. 10
0021-9193/06/$08.00+0 doi:10.1128/JB.188.10.3463-3469.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Catherine S. Ryan,1
Lori A. S. Snyder,2
Nigel J. Saunders,2
Julian I. Rood,1 and
John K. Davies1*
Australian Bacterial Pathogenesis Program, Department of Microbiology, Monash University, Clayton, VIC 3800, Australia,1 Bacterial Pathogenesis and Functional Genomics Group, Sir William Dunn School of Pathology, University of Oxford, South Parks Road, Oxford OX1 3RE, United Kingdom2
Received 17 October 2005/ Accepted 29 December 2005
|
|
|---|
|
|
|---|
We have searched the genome sequence of Neisseria gonorrhoeae strain FA1090 (GenBank accession number AE004969) for the presence of genes encoding sigma factors. As expected, there is a gene, rpoD, encoding sigma-70 (NGO0999). Only two intact genes encoding alternative sigma factors, rpoH (NGO0288) and ecf (NGO1944), were found. Laskos et al. (11) have shown that there is also an inactive RpoN-like sigma factor, RLS, although this sequence feature (NGO1766) has not been annotated as such in the publicly available annotations of the N. gonorrhoeae strain FA1090 genome sequence (GenBank accession number AE004969; annotations at www.stdgen.lanl.gov,cmr.tigr.org/tigr-scripts/CMR/CmrHomePage.cgi, and www.ncbi.nlm.nih.gov/genomes/lproks.cgi). The ecf gene encodes a member of the extracytoplasmic function (ECF) family of sigma factors. As the name suggests, ECF sigma factors from different bacterial species appear to respond specifically to a variety of extracytoplasmic stimuli. Characteristically, their activity is controlled by anti-sigma factors, and they control relatively small regulons (1). In most examples of this system investigated to date, the ECF sigma factors regulate not only their own expression but also that of the genes encoding the cognate anti-sigma factor, which are located in the same operon.
We have used DNA microarrays to assess the role of the sigma factor Ecf in gonococci and show that it controls a small regulon which contains the msrAB gene, which encodes an unusual methionine sulfoxide reductase.
|
|
|---|
[F endA1 thi-1 hsdR17 supE44 relA1
lacU169 (
80
lacZM15)] was used for genetic manipulation of constructs. E. coli strains were grown in Luria-Bertani (LB) broth (Difco) at 37°C with agitation or on LB agar plates supplemented with 1.5% agar, and where appropriate, with the following antibiotics: 100 µg/ml ampicillin, 25 µg/ml kanamycin, 150 µg/ml erythromycin, or 10 µg/ml tetracycline. Gonococcal strains were grown on GC agar base (Oxoid) or GC broth as described previously (7). Where appropriate, media were supplemented with 7 µg/ml erythromycin, 12.5 µg/ml tetracycline, or 40 µg/ml kanamycin. Transformation and conjugation experiments with N. gonorrhoeae were performed as described previously (12) with the exception that conjugation experiments involved mixing 5 x 108 donor cells and 1.5 x 109 recipient cells. |
View this table: [in a new window] |
TABLE 1. Bacterial strains and plasmids
|
. Ampicillin- and kanamycin-resistant colonies were screened by PCR for the presence of the insert. The resultant plasmid, pJKD2603 (Table 1), was linearized with SphI, and the mutated gene was recombined into the genome of N. gonorrhoeae strain FA1090 by homologous recombination. Transformants were isolated by selection for kanamycin resistance, and the presence of the aphA1 cassette in ecf was confirmed by PCR. The ecf mutant was designated JKD5069. |
View this table: [in a new window] |
TABLE 2. Oligonucleotide primers
|
For investigation of the Ecf regulon, gonococcal strain JKD5062 was grown to mid-exponential phase (optical density at 560 nm = 0.6), and the culture was then split between two prewarmed flasks. One of the cultures was induced with 2 mM IPTG for 30 min. Samples from cultures were harvested and stored in RNAlater RNA stabilization solution (Ambion). The RNA from these samples was extracted using the RNeasy midi kit (QIAGEN) with an optional on-column DNase I digestion step according to the manufacturer's instructions. The quality and quantity of RNA were determined by gel electrophoresis and spectrophotometry.
cDNA synthesis and fluorescent labeling. Total RNA (30 µg) was mixed with 30 µg of random hexamers, heated to 70°C for 10 min, and then rapidly chilled on ice. To this mixture 0.5 µl (20 units) of RNasin (Promega), 6 µl of Superscript II buffer (Life Technologies, Inc.), 3 µl dithiothreitol (DTT), 1 mM dATP, 1 mM dCTP, 1 mM dGTP, 0.4 mM dTTP (Promega), 0.6 mM amino allyl-dUTP (Sigma), and 2 µl (400 units) Superscript II reverse transcriptase (Life Technologies, Inc) were added, and the mixture was incubated for 2.5 h at 42°C. Reactions were terminated by the addition of 10 µl of 1 M NaOH and 10 µl of 0.5 M EDTA, and the mixtures were heated to 65°C for 15 min and neutralized by the addition of 25 µl of 1 M Tris, pH 7.4. Unincorporated amino allyl-dUTP was removed with Microcon 30 columns (Millipore) according to the manufacturer's instructions. The purified cDNA was concentrated to 12 µl, and the quality of labeling and quantity of cDNA were determined using a NanoDrop ND-100 spectrophotometer (NanoDrop Technologies). Cy3 and Cy5 dyes were dissolved in dimethyl sulfoxide (Sigma) to a final concentration of 14 µg/µl, and a 4-µl aliquot was added to each cDNA sample in 0.1 M sodium bicarbonate buffer (pH 9). Reaction mixtures were incubated for 1 h at room temperature in the dark, labeled cDNA was purified using Microcon columns, and the eluted samples were concentrated to 10 µl using a SpeedVac SVC100 (Savant). Cy3 and Cy5 dye incorporation was estimated by spectrophotometry.
Microarray hybridization and data analysis. The labeled cDNA was hybridized to a pan-Neisseria microarray which contains 2,704 PCR products, spotted in triplicate, corresponding to potential coding sequences from N. gonorrhoeae strain FA1090, N. meningitidis strains MC58 and Z2491, the N. gonorrhoeae strain MS11 gonococcal genetic island (8), and various controls. Details of the construction of this array will be published elsewhere, but details and arrays can be obtained from the corresponding author. Prior to prehybridization, the slide was plunged into a 95°C water bath for 2 min, centrifuged for 5 min at 2,000 x g, and used immediately. Prehybridization was carried out in a volume of 30 µl containing 25% formamide, 5x SSC (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate), 0.1% sodium dodecyl sulfate (SDS), 10 mg/ml bovine serum albumin (BSA; fraction V), and 30 µg of herring sperm DNA (Promega) under a coverslip in a humidified Corning CMT hybridization chamber for 45 min at 42°C. The slide was rinsed in water, dried by centrifugation, and used immediately for hybridization. The labeled cDNA was added to the hybridization fluid (prehybridization solution without BSA) in a total volume of 30 µl and denatured at 95°C for 5 min. The hybridization fluid was placed on the surface of the microarray under a coverslip. The slide was enclosed in the hybridization chamber and submerged in a 42°C water bath overnight. After hybridization, slides were washed once in 2x SSC-0.1% SDS for 5 min at 42°C, once in 0.1x SSC-0.1% SDS for 10 min at room temperature, and four times in 0.1x SSC for 1 min at room temperature. The slides were rinsed with water for 10 seconds, and excess fluid on the surface of the slide was removed by centrifugation before scanning. A GMS 418 array scanner (Genetic Microsystem) was used to scan the images. The Cy3 and Cy5 images were combined, and fluorescent and background intensity for each spot was determined using ImaGene version 5 software (BioDiscovery) as outlined previously (2). Data from poor spots that were manually or automatically flagged in ImaGene were removed from further analysis.
The individual ImaGene data files were uploaded to a web site created with BASE (19) and converted to a common BASE format using a series of custom-made applications. The spot intensities were found to be most reliable when no background correction was performed. Analysis was done using Bioconductor (5) and Limma (23). The normalization of the data to remove various biases involved two steps. Firstly, each array was normalized independently using print-tip loess normalization (Y. H. Yang, S. Dudoit, P. Luu, and T. P. Speed, presented at SPIE BiOS 2001, San Jose, Calif.). Secondly, diagnostic plots suggested a variation in scale between arrays, so the log ratios were scaled in such a way that each array had the same median-absolute-deviation. The normalized data were then used to fit a linear model (23) for each gene using a generalized least-squares method which took into account the correlation between replicate spots (24). The coefficient of the fitted model for each gene describes the inferred difference in RNA expression between the two strains. Empirical Bayes method was then used to calculate the moderated t statistics and associated P values. The P values were adjusted for multiple testing using the false-discovery rate. Genes with an absolute ratio greater than 1.5-fold and significant at the 0.001 level were selected as differentially expressed.
Quantitative real-time reverse transcriptase PCR (RT-PCR). cDNA was generated from 5 µg of the total RNA samples used for the microarray experiments. The RNA preparation was subjected to a second DNase I treatment before cDNA synthesis. The cDNA synthesis was performed as described above with the exception that a 1 mM concentration of each dNTP and 7.8 µg of random hexamers were used. Oligonucleotide primer pairs specific for each gene of interest were designed using ABI PRISM Primer Express software (Applied Biosystems) and are listed in Table 2. To quantitate the cDNA, a gene-specific standard curve method was employed using serial dilutions of strain FA1090 genomic DNA as the template. All assays included 12.5 µl of SYBR Green PCR master mix (Applied Biosystems), 2 µl of each primer (0.5 nM final), 2 µl of template, and diethyl pyrocarbonate-treated water made up to a final volume of 25 µl. Controls lacked reverse transcriptase or template. Reactions were run on an ABI 7700 sequence detection system (Applied Biosystems), and recA was used as the normalizer for all reactions. All reactions amplified a single product, as determined by melting curve analysis (Applied Biosystems).
Northern blots. Northern blotting was performed on 10 µg of the total RNA preparations used in the microarray experiments. RNA was electrophoresed on formaldehyde-MOPS (morpholinepropanesulfonic acid) 1% agarose gels and transferred to Hybond N nylon membranes (Amersham) as previously described (20). The transferred RNA was fixed to the nylon membrane with the auto-cross-link setting on a UV Stratalinker 1800 (Stratagene). Probe labeling and detection were performed using a digoxigenin nonradioactive DNA labeling and detection kit (Roche Diagnostics) according to the manufacturer's instructions. Probes were amplified by PCR from each gene using oligonucleotides listed in Table 2.
Identification of transcription start points.
Oligonucleotide primer 23762 (Table 2), which is complementary to the msrAB sequence, was 5'-end labeled with 30 µCi of [
-32P]dATP (Amersham) and T4 polynucleotide kinase. Primer extension reactions were performed as described in the instructions for the use of the primer extension system (Promega) using avian myeloblastosis virus reverse transcriptase. The precipitated primer extension products were subjected to electrophoresis in an 8% polyacrylamide gel containing 8 M urea, next to a sequencing ladder generated from plasmid pJKD2623 (Table 1) with oligonucleotide primer 23762 (Table 2). Plasmid pJKD2623 was constructed by amplifying a PCR product encompassing the promoter region of msrAB with primers 23764 and 23763 (Table 2), T4 filling, and ligation to HincII-cut pUC18 vector. Primer extension gels were dried onto chromatography paper and subjected to autoradiography.
To confirm the transcription start points, 5' RACE (rapid amplification of cDNA ends) experiments were performed (27) with the same RNA preparations used in the microarray experiments. cDNA synthesis was performed as described previously (27) using a specific oligonucleotide primer, 23111 (Table 2), located 150 to 200 bp downstream from the initiation codon of NGO1948. Briefly, 5'-phosphorylated, 3'-blocked primer DT88 (27) (Geneworks) (Table 2) was anchored to the 5' end of the cDNA, and a first-round PCR product was amplified using primer DT89 (Table 2), which was complementary to DT88 (27), and oligonucleotide 23111 (Table 2). A second round of PCR used oligonucleotides DT89 and a nested primer, 23489 (Table 2), that was internal to 23111. Controls with and without reverse transcriptase and T4 RNA ligase were included. Purified PCR product was sequenced to determine the 5' end of the transcript, which represents the transcription start point.
Western blots. Western blotting was performed as described previously (11). The samples were prepared from whole-cell extracts of induced and noninduced gonococcal strains JKD5062 and JKD5064 (Table 1) grown in liquid medium. Prior to electrophoresis, 1 x 109 cells were harvested by centrifugation at 5,000 x g at 4°C for 5 min, resuspended in 500 µl of phosphate-buffered saline and SDS sample buffer (0.2 M Tris, 20% glycerol, 25% SDS, 10% ß-mercaptoethanol, 0.1% bromophenol blue; pH 6.8), and boiled for 5 min. For equivalent loading, total protein concentration was determined by the microplate procedure using the bicinchoninic acid protein assay kit (Pierce). Primary polyclonal anti-MsrAB rabbit antiserum was used at a 1:20,000 dilution (22). Peroxidase-conjugated anti-rabbit immunoglobulin was used as the secondary antibody (Chemicon) at a 1:3,000 dilution. The Western blot was completed using the protocol provided with the ECL detection reagent kit (Amersham).
|
|
|---|
In the absence of any knowledge of the physiological signals to which Ecf responds, we reasoned that activation of the regulon could be obtained by overexpression of the sigma factor. To this end we used an IPTG-inducible Ptrc promoter in the E. coli/N. gonorrhoeae Hermes shuttle vector system (10), in the background of the ecf mutant JKD5069. The strain constructed, JKD5062, allows exclusive and inducible expression of a plasmid-borne copy of the ecf gene (see Materials and Methods). We again used microarrays to measure changes in gene expression when the ecf gene in strain JKD5062 was induced for 30 min compared to an uninduced control.
Four genes, NGO1946, NGO1947, NGO1948, and NGO2059, were up-regulated along with ecf (NGO1944) itself when ecf was overexpressed (Table 3). As might be expected with the overexpression of an alternative sigma factor, no down-regulated genes were detected. Quantitative real-time RT-PCR was used to validate the microarray data for three of the up-regulated genes (Table 3). NGO1946, NGO1947, and NGO1948 are annotated as encoding proteins with no known function. The presence of possible transmembrane domains, as predicted by PSORT (15), suggests that NGO1948 might encode an integral membrane protein. A conserved-domain search (13) with the derived amino acid sequence from NGO1948 also revealed the presence of a motif that is found in the DoxD-like family of proteins. DoxD is a subunit of quinone oxidoreductase, which is involved in the oxidation of sulfur (14, 17). NGO1947 is predicted by PSORT (15) to encode a periplasmic protein. Genes encoding proteins closely related to Ecf and the genes NGO1946, NGO1947, and NGO1948 can be found clustered in other bacterial genomes, where the gene order also seems to have been at least partially conserved (data not shown). The expression ratio for NGO1945, the gene located between ecf and NGO1946, was just below the 1.5-fold cutoff used in these experiments.
|
View this table: [in a new window] |
TABLE 3. Genes up-regulated in strain JKD5062 as a result of overexpression of Ecf
|
One component of the anti-sigma factors controlling ECF activity is an inner membrane protein with at least one transmembrane domain (18). NGO1948 encodes a protein with a predicted transmembrane domain. However, as is the case with the sigma E proteins from E. coli and Streptomyces coelicolor, multiple genes can have a role in the anti-sigma function (3, 29). Taking into account the observation that Ecf appears to be relatively inactive in strain FA1090 during the exponential phase of growth, the possibility exists that NGO1946 to NGO1948 encode proteins that function to inhibit Ecf activity. In addition, the genomic arrangement of NGO1944 (ecf), NGO1945, NGO1946, NGO1947, and NGO1948 suggests the possibility of cotranscription of at least some of these genes. Genes encoding several ECF sigma factors have been shown to be cotranscribed with a gene(s) encoding a specific anti-sigma factor(s) (3, 9, 21, 25). Northern hybridization was performed with total RNA isolated from strain JKD5062, with or without induction of ecf. A strong signal was observed from RNA prepared after induction of ecf with probes derived from both msrAB and NGO1948, suggesting the presence of Ecf-controlled promoters upstream of these genes (Fig. 1). The probe specific for NGO1948 hybridized to a major 0.9-kb band and a minor 1.9-kb mRNA transcript. The presence of the latter suggests that at least NGO1946, NGO1947, and NGO1948 can be cotranscribed. Larger transcripts are unlikely to have been detected in this experiment. The major 0.9-kb transcript corresponds in nucleotide length to NGO1948 and NGO1947, which suggests the presence of a transcriptional attenuator at the end of NGO1947. An inverted repeat sequence can be found in this intergenic region (data not shown). The minimal spacing between the remaining genes in the cluster (Fig. 1) suggests that they may be cotranscribed. An inverted repeat sequence that may act as a terminator is found downstream of NGO1943. With this strain, we were unable to look directly for cotranscription of the other genes with ecf, as this gene is transcribed from the multicopy vector in strain JKD5062. As expected, given the orientation of the genes neighboring msrAB (Fig. 1), this probe hybridized to a band of the correct size for a single gene transcript.
![]() View larger version (19K): [in a new window] |
FIG. 1. Detection of msrAB and NGO1948 transcripts by Northern hybridization. Total RNA was isolated from N. gonorrhoeae strain JKD5062, both with (In) and without (Non) overexpression of ecf, and separated on a 1.2% denaturing agarose gel. The separated RNA species were probed with an msrAB-specific (A) or an NGO1948-specific (B) probe. The msrAB-specific probe detected a 1.5-kb transcript (A, thin arrow). The NGO1948-specific probe detected two transcripts (B, thin arrows) of 0.9 kb and 1.9 kb. Nonspecific bands (bold arrows) were artifacts comigrating with rRNA. Below the gels are diagrammatic representations of the putative genes in the region, with approximate sizes of transcripts (B), in nucleotides.
|
![]() View larger version (16K): [in a new window] |
FIG. 2. Identification of transcription start points of Ecf-regulated genes. (A) Primer extension mapping of the msrAB TSP. The 5' terminus of the msrAB transcript was determined using RNA prepared from cultures induced for ecf overexpression (In) and noninduced cultures (Non). The lanes marked T, G, C, and A are the sequencing ladder that was generated with the primer used for primer extension and pJKD2623 (Table 1), which contains the putative promoter. In the sequence on the left, the arrow indicates the transcription start point. (B) 5' RACE mapping of the NGO1948 TSP. The sequence upstream of NGO1948 is shown; the start codon for the putative gene is in bold. (C) Alignment of the promoter regions of NGO1948 and msrAB with the putative ecf consensus sequence shown below. The 10 and 35 regions are in bold.
|
![]() View larger version (16K): [in a new window] |
FIG. 3. Western blot of cell extracts from N. gonorrhoeae using a polyclonal antiserum to detect MsrAB. Marker protein sizes are in kDa. Lanes 1, 3, 5, and 7 contain extracts from ecf overexpression strain JKD5062, while lanes 2, 4, 6, and 8 contain extracts from control strain JKD5064, carrying an empty vector. Lanes 1, 2, 5, and 6 contain extracts from cultures that were induced with IPTG, while lanes 3, 4, 7, and 8 contain extracts from noninduced cultures. Lanes 1 to 4 contain extracts that were sampled after 30 min, while lanes 5 to 8 contain extracts that were sampled after 60 min.
|
This work was supported by an NHMRC Program Grant to the Australian Bacterial Pathogenesis Program. L.A.S.S. was supported by a Wellcome Trust project grant awarded to N.J.S. I.C.G. was the recipient of a Monash University Postgraduate Scholarship. We have made use of data from the Gonococcal Genome Sequencing Project, which was supported by USPHS/NIH grant AI-38399.
Present address: Department of Microbiology and Immunology, School of Biomedical, Biomolecular and Chemical Sciences, The University of Western Australia, QE11 Medical Centre, Crawley, WA 6009, Australia. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»