The genome structure of Pseudomonas putida: high-resolution mapping and microarray analysis

Differential transcriptional response to antibiotics by Pseudomonas putida DOT-T1E

Multi-drug resistant bacteria are a major threat to humanity, especially because the current battery of known antibiotics is not sufficient to combat infections produced by these microbes. Therefore, the study of how current antibiotics act and how bacteria defend themselves against antibiotics is of critical importance. Pseudomonas putida DOT-T1E exhibits an impressive array of RND efflux pumps, which confer this microorganism high resistance to organic solvents and antibiotics that would kill most other microorganisms. We have chosen DOT-T1E as a model microbe to study the microbial responses to a wide battery of antibiotics (chloramphenicol, rifampicin, tetracycline, ciprofloxacin, ampicillin, kanamycin, spectinomycin and gentamicin). Ribonucleic acid sequencing (RNA)-seq analyses revealed that each antibiotic provokes a unique transcriptional response profile in DOT-T1E. While many of the genes identified were related to known antibiotic targets, others were unrelated or encoded hypothetical proteins. These results indicate that our knowledge of antibiotic resistance mechanisms is still partial. We also identified 138 new small RNAs (sRNAs) in DOT-T1E, dramatically adding to the 16 that have been previously described. Importantly, our results reveal that a correlation exists between the expression of messenger RNA and sRNA, indicating that some of these sRNAs are likely involved in fine tuning the expression of antibiotic resistance genes. Taken together, these findings open new frontiers in the fight against multi-drug resistant bacteria and point to the potential use of sRNAs as novel antimicrobial targets.

Regulation of glucose metabolism in Pseudomonas: the phosphorylative branch and entner-doudoroff enzymes are regulated by a repressor containing a sugar isomerase domain

In Pseudomonas putida, genes for the glucose phosphorylative pathway and the Entner-Doudoroff pathway are organized in two operons; one made up of the zwf, pgl, and eda genes and another consisting of the edd, glk, gltR2, and gltS genes. Divergently with respect to the edd gene is the gap-1 gene. Expression from P(zwf), P(edd), and P(gap) is modulated by HexR in response to the availability of glucose in the medium. To study the regulatory process in greater detail we purified HexR and showed that it is a monomer in solution. Electrophoretic mobility shift assays and isothermal titration calorimetry assays were done showing that HexR recognizes the P(edd), P(zwf), and P(gap-1) promoters with affinity in the nanomolar range. DNA footprinting assays identified the binding site between +30 and +1 at P(zwf), between +16 and +41 at P(edd), and between -6 and +18 at P(gap-1). Based on DNA sequence alignment of the target sites and isothermal titration calorimetry data, two monomers of HexR bind to a pseudopalindrome with a consensus sequence of 5'-TTGTN(7-8)ACAA-3'. Binding of the Entner-Doudoroff pathway intermediate 2-keto-3-deoxy-6-phosphogluconate to HexR released the repressor from its target operators, whereas other chemicals such as glucose, glucose 6-phosphate, and 6-phosphogluconate did not induce complex dissociation. The phosphorylated effector is likely to be recognized by a sugar isomerase domain located at the C-terminal end of HexR, whereas the helix-turn-helix DNA binding domain of HexR exhibits high similarity to proteins of the RpiR family of regulators.

Use of complex DNA and antibody microarrays as tools in functional analyses

While the deciphering of basic sequence information on a genomic scale is yielding complete genomic sequences in ever-shorter intervals, experimental procedures for elucidating the cellular effects and consequences of the DNA-encoded information become critical for further analyses. In recent years, DNA microarray technology has emerged as a prime candidate for the performance of many such functional assays. Technically, array technology has come a long way since its conception some 15 years ago, initially designed as a means for large-scale mapping and sequencing.

The basic arrangement, however, could be adapted readily to serve eventually as an analytical tool in a large variety of applications. On their own or in combination with other methods, microarrays open up many new avenues of functional analysis.

Growth phase-dependent expression of the Pseudomonas putida KT2440 transcriptional machinery analysed with a genome-wide DNA microarray

Bacterial transcriptional networks are built on a hierarchy of regulators, on top of which lie the components of the RNA polymerase (in particular the sigma factors) and the global control elements, which play a pivotal role. We have designed a genome-wide oligonucleotide-based DNA microarray for Pseudomonas putida KT2440. In combination with real-time reverse transcription polymerase chain reaction (RT-PCR), we have used it to analyse the expression pattern of the genes encoding the RNA polymerase subunits (the core enzyme and the 24 sigma factors), and various proteins involved in global regulation (Crc, Lrp, Fur, Anr, Fis, CsrA, IHF, HupA, HupB, HupN, BipA and several MvaT-like proteins), during the shift from exponential growth in rich medium into starvation and stress brought about by the entry into stationary phase. Expression of the genes encoding the RNA polymerase core and the vegetative sigma factor decreased in stationary phase, while that of sigma(S) increased. Data obtained for sigma(N), sigma(H), FliA and for the 19 extracytoplasmic function (ECF)-like sigma factors suggested that their mRNA levels change little upon entry into stationary phase. Expression of Crc, BipA, Fis, HupB, HupN and the MvaT-like protein PP3693 decreased in stationary phase, while that of HupA and the MvaT-like protein PP3765 increased significantly. Expression of IHF was indicative of post-transcriptional control. These results provide the first global study of the expression of the transcriptional machinery through the exponential stationary-phase shift in P. putida.

Functional genomics of stress response in Pseudomonas putida KT2440

The metabolically versatile soil bacterium Pseudomonas putida has to cope with numerous abiotic stresses in its habitats. The stress responses of P. putida KT2440 to 4 degrees C, pH 4.5, 0.8 M urea, and 45 mM sodium benzoate were analyzed by determining the global mRNA expression profiles and screening for stress-intolerant nonauxotrophic Tn5 transposon mutants. In 392 regulated genes or operons, 36 gene regions were differentially expressed by more than 2.5-fold, and 32 genes in 23 operons were found to be indispensable for growth during exposure to one of the abiotic stresses. The transcriptomes of the responses to urea, benzoate, and 4 degrees C correlated positively with each other but negatively with the transcriptome of the mineral acid response. The CbrAB sensor kinase, the cysteine synthase CysM, PcnB and VacB, which control mRNA stability, and BipA, which exerts transcript-specific translational control, were essential to cope with cold stress. The cyo operon was required to cope with acid stress. A functional PhoP, PtsP, RelA/SpoT modulon, and adhesion protein LapA were necessary for growth in the presence of urea, and the outer membrane proteins OmlA and FepA and the phosphate transporter PstBACS were indispensable for growth in the presence of benzoate. A lipid A acyltransferase (PP0063) was a mandatory component of the stress responses to cold, mineral acid, and benzoate. Adaptation of the membrane barrier, uptake of phosphate, maintenance of the intracellular pH and redox status, and translational control of metabolism are key mechanisms of the response of P. putida to abiotic stresses.

Regulation of superoxide stress in Pseudomonas putida KT2440 is different from the SoxR paradigm in Escherichia coli

In Escherichia coli, the SoxR regulon orchestrates genes for defense against certain types of oxidative stress through the SoxR-regulated synthesis of the SoxS transcription activator. The Pseudomonas putida genome did not reveal a clear soxS homolog. The P. putida SoxR protein appears to be functional: its expression in an E. coli DeltasoxR strain restored the paraquat inducibility of soxS. Of nine candidate P. putida oxidative stress genes, which are known to be SoxR regulon in E. coli, tested for response to superoxide or nitric oxide, fumC-1, sodA, zwf-1, and particularly fpr, encoding ferredoxin:NADP(+) reductase, were induced, all independent of P. putida soxR. Disruption of the fpr and finR, a regulatory protein that is required for paraquat-dependent expression of the fpr, resulted in more oxidative stress sensitivity. However, a P. putida soxR-deletion strain had normal resistance to the superoxide-generating agent paraquat. The data presented here show that the genetic responses to superoxide stress in P. putida differ markedly from those seen in E. coli and Salmonella, and the role of P. putida soxR remains to be established.

Projector 2: contig mapping for efficient gap-closure of prokaryotic genome sequence assemblies

With genome sequencing efforts increasing exponentially, valuable information accumulates on genomic content of the various organisms sequenced. Projector 2 uses (un)finished genomic sequences of an organism as a template to infer linkage information for a genome sequence assembly of a related organism being sequenced. The remaining gaps between contigs for which no linkage information is present can subsequently be closed with direct PCR strategies. Compared with other implementations, Projector 2 has several distinctive features: a user-friendly web interface, automatic removal of repetitive elements (repeat-masking) and automated primer design for gap-closure purposes. Moreover, when using multiple fragments of a template genome, primers for multiplex PCR strategies can also be designed. Primer design takes into account that, in many cases, contig ends contain unreliable DNA sequences and repetitive sequences. Closing the remaining gaps in prokaryotic genome sequence assemblies is thereby made very efficient and virtually effortless. We demonstrate that the use of single or multiple fragments of a template genome (i.e. unfinished genome sequences) in combination with repeat-masking results in mapping success rates close to 100%. The web interface is freely accessible at http://molgen.biol.rug.nl/websoftware/projector2.

Creation of a minimal tiling path of genomic clones for Drosophila: provision of a common resource

On the basis of shotgun subclone libraries used in the sequencing of the Drosophila melanogaster genome, a minimal tiling path of subclones across much of the genome was determined. About 320,000 shotgun clones for chromosomes X(12–20), 2R, 2L, 3R, and 4 were available from the Berkeley Drosophila Genome Project. The clone inserts have an average length of 3.4 kb and are amenable to standard PCR amplification. The resulting tiling path covers 86.2% of chromosome X(12–20), 86.2% of chromosomal arm 2R, 79.0% of 2L, 89.6% of 3R, and 80.5% of chromosome 4. In total, the 25,135 clones represent 76.7 Mb—equivalent to about 67% of the genome —and would be suitable for producing a microarray on a single slide.

Projector: automatic contig mapping for gap closure purposes

Projector was designed for automatic positioning of contigs from an unfinished prokaryotic genome onto a template genome of a closely related strain or species. Projector mapped 84 contigs of Lactococcus lactis MG1363 (corresponding to 81% of the assembly nucleotides) against the genome of L.lactis IL1403. Ninety three percent of subsequent gap closure PCRs were successful. Moreover, a significant improvement in the N50 and N80 values (describing the assembly quality) was observed after the use of Projector. Because increasing numbers of bacterial genomes are being sequenced, Projector provides an efficient method to close a significant number of remaining gaps in the late stages of a genome sequencing project.

Analysis of nutrient-dependent transcript variations in Neurospora crassa

Nutrient-dependent variations in transcript levels of the filamentous fungus Neurospora crassa were studied on a microarray containing some 4700 cDNAs. Cells were grown in minimal and acetate medium. The isolated RNA was analyzed in comparison to the results obtained upon the hybridization of samples prepared from the RNA of cells grown in full medium. Altogether, 160 cDNA clones exhibited significant variations, falling into five distinct subgroups of very similar transcription profiles. This is indicative of the occurrence of a high degree of co-regulation of genes in N. crassa. Especially the regulation of the expression of proteins involved in metabolic pathways was found to be strongly regulated at the RNA level.