Functional and evolutionary significance of unknown genes from uncultivated taxa.

Many of the Earth's microbes remain uncultured and understudied, limiting our understanding of the functional and evolutionary aspects of their genetic material, which remain largely overlooked in most metagenomic studies1. Here we analysed 149,842 environmental genomes from multiple habitats2-6 and compiled a curated catalogue of 404,085 functionally and evolutionarily significant novel (FESNov) gene families exclusive to uncultivated prokaryotic taxa. All FESNov families span multiple species, exhibit strong signals of purifying selection and qualify as new orthologous groups, thus nearly tripling the number of bacterial and archaeal gene families described to date. The FESNov catalogue is enriched in clade-specific traits, including 1,034 novel families that can distinguish entire uncultivated phyla, classes and orders, probably representing synapomorphies that facilitated their evolutionary divergence. Using genomic context analysis and structural alignments we predicted functional associations for 32.4% of FESNov families, including 4,349 high-confidence associations with important biological processes. These predictions provide a valuable hypothesis-driven framework that we used for experimental validatation of a new gene family involved in cell motility and a novel set of antimicrobial peptides. We also demonstrate that the relative abundance profiles of novel families can discriminate between environments and clinical conditions, leading to the discovery of potentially new biomarkers associated with colorectal cancer. We expect this work to enhance future metagenomics studies and expand our knowledge of the genetic repertory of uncultivated organisms.

Prevalence and Specificity of Chemoreceptor Profiles in Plant-Associated Bacteria.

Chemosensory pathways are among the most abundant prokaryotic signal transduction systems, allowing bacteria to sense and respond to environmental stimuli. Signaling is typically initiated by the binding of specific molecules to the ligand binding domain (LBD) of chemoreceptor proteins (CRs). Although CRs play a central role in plant-microbiome interactions such as colonization and infection, little is known about their phylogenetic and ecological specificity. Here, we analyzed 82,277 CR sequences from 11,806 representative microbial species covering the whole prokaryotic phylogeny, and we classified them according to their LBD type using a de novo homology clustering method. Through phylogenomic analysis, we identified hundreds of LBDs that are found predominantly in plant-associated bacteria, including several LBDs specific to phytopathogens and plant symbionts. Functional annotation of our catalogue showed that many of the LBD clusters identified might constitute unknown types of LBDs. Moreover, we found that the taxonomic distribution of most LBD types that are specific to plant-associated bacteria is only partially explained by phylogeny, suggesting that lifestyle and niche adaptation are important factors in their selection. Finally, our results show that the profile of LBD types in a given genome is related to the lifestyle specialization, with plant symbionts and phytopathogens showing the highest number of niche-specific LBDs. The LBD catalogue and information on how to profile novel genomes are available at https://github.com/compgenomicslab/CRs. IMPORTANCE Considering the enormous variety of LBDs at sensor proteins, an important question resides in establishing the forces that have driven their evolution and selection. We present here the first clear demonstration that environmental factors play an important role in the selection and evolution of LBDs. We were able to demonstrate the existence of LBD families that are highly enriched in plant-associated bacteria but show a wide phylogenetic spread. These findings offer a number of research opportunities in the field of single transduction, such as the exploration of similar relationships in chemoreceptors of bacteria with a different lifestyle, like those inhabiting or infecting the human intestine. Similarly, our results raise the question whether similar LBD types might be shared by members of different sensor protein families. Lastly, we provide a comprehensive catalogue of CRs classified by their LBD region that includes a large number of putative new LBD types.

Pseudomonas syringae pv. tomato exploits light signals to optimize virulence and colonization of leaves.

Light is pervasive in the leaf environment, creating opportunities for both plants and pathogens to cue into light as a signal to regulate plant-microbe interactions. Light enhances plant defences and regulates opening of stomata, an entry point for foliar bacterial pathogens such as Pseudomonas syringae pv. tomato DC3000 (PsPto). The effect of light perception on gene expression and virulence was investigated in PsPto. Light induced genetic reprogramming in PsPto that entailed significant changes in stress tolerance and virulence. Blue light-mediated up-regulation of type three secretion system genes and red light-mediated down-regulation of coronatine biosynthesis genes. Cells exposed to white light, blue light or darkness before inoculation were more virulent when inoculated at dawn than dusk probably due to an enhanced entry through open stomata. Exposure to red light repressed coronatine biosynthesis genes which could lead to a reduced stomatal re-opening and PsPto entry. Photoreceptor were required for the greater virulence of light-treated and dark-treated PsPto inoculated at dawn as compared to dusk, indicating that these proteins sense the absence of light and contribute to priming of virulence in the dark. These results support a model in which PsPto exploits light changes to maximize survival, entry and virulence on plants.

Exploring new roles for the rpoS gene in the survival and virulence of the fire blight pathogen Erwinia amylovora.

Erwinia amylovora causes fire blight in economically important plants of the family Rosaceae. This bacterial pathogen spends part of its life cycle coping with starvation and other fluctuating environmental conditions. In many Gram-negative bacteria, starvation and other stress responses are regulated by the sigma factor RpoS. We obtained an E. amylovora rpoS mutant to explore the role of this gene in starvation responses and its potential implication in other processes not yet studied in this pathogen. Results showed that E. amylovora needs rpoS to develop normal starvation survival and viable but nonculturable (VBNC) responses. Furthermore, this gene contributed to stationary phase cross-protection against oxidative, osmotic, and acid stresses and was essential for cross-protection against heat shock, but nonessential against acid shock. RpoS also mediated regulation of motility, exopolysaccharide synthesis, and virulence in immature loquats, but not in pear plantlets, and contributed to E. amylovora survival in nonhost tissues during incompatible interactions. Our results reveal some unique roles for the rpoS gene in E. amylovora and provide new knowledge on the regulation of different processes related to its ecology, including survival in different environments and virulence in immature fruits.

The type II secretion system (Xcp) of Pseudomonas putida is active and involved in the secretion of phosphatases.

The genome of the Gram-negative bacterium Pseudomonas putida harbours a complete set of xcp genes for a type II protein secretion system (T2SS). This study shows that expression of these genes is induced under inorganic phosphate (Pi ) limitation and that the system enables the utilization of various organic phosphate sources. A phosphatase of the PhoX family, previously designated UxpB, was identified, which was produced under low Pi conditions and transported across the cell envelope in an Xcp-dependent manner demonstrating that the xcp genes encode an active T2SS. The signal sequence of UxpB contains a twin-arginine translocation (Tat) motif as well as a lipobox, and both processing by leader peptidase II and Tat dependency were experimentally confirmed. Two different tat gene clusters were detected in the P. putida genome, of which one, named tat-1, is located adjacent to the uxpB and xcp genes. Both Tat systems appeared to be capable of transporting the UxpB protein. However, expression of the tat-1 genes was strongly induced by low Pi levels, indicating a function of this system in survival during Pi starvation.

A bacterial cysteine protease effector protein interferes with photosynthesis to suppress plant innate immune responses.

The bacterial pathogen Pseudomonas syringae pv tomato DC3000 suppresses plant innate immunity with effector proteins injected by a type III secretion system (T3SS). The cysteine protease effector HopN1, which reduces the ability of DC3000 to elicit programmed cell death in non-host tobacco, was found to also suppress the production of defence-associated reactive oxygen species (ROS) and callose when delivered by Pseudomonas fluorescens heterologously expressing a P. syringae T3SS. Purified His(6) -tagged HopN1 was used to identify tomato PsbQ, a member of the oxygen evolving complex of photosystem II (PSII), as an interacting protein. HopN1 localized to chloroplasts and both degraded PsbQ and inhibited PSII activity in chloroplast preparations, whereas a HopN1(D299A) non-catalytic mutant lost these abilities. Gene silencing of NtPsbQ in tobacco compromised ROS production and programmed cell death by DC3000. Our data reveal PsbQ as a contributor to plant immunity responses and a target for pathogen suppression.

Genome-wide analysis of the response of Dickeya dadantii 3937 to plant antimicrobial peptides.

Antimicrobial peptides constitute an important factor in the defense of plants against pathogens, and bacterial resistance to these peptides have previously been shown to be an important virulence factor in Dickeya dadantii, the causal agent of soft-rot disease of vegetables. In order to understand the bacterial response to antimicrobial peptides, a transcriptional microarray analysis was performed upon treatment with sub-lethal concentration of thionins, a widespread plant peptide. In all, 36 genes were found to be overexpressed, and were classified according to their deduced function as i) transcriptional regulators, ii) transport, and iii) modification of the bacterial membrane. One gene encoding a uricase was found to be repressed. The majority of these genes are known to be under the control of the PhoP/PhoQ system. Five genes representing the different functions induced were selected for further analysis. The results obtained indicate that the presence of antimicrobial peptides induces a complex response which includes peptide-specific elements and general stress-response elements contributing differentially to the virulence in different hosts.

Identification and characterization of the PhhR regulon in Pseudomonas putida.

Pseudomonas putida is a soil microorganism that utilizes aromatic amino acids present in root exudates as a nitrogen source. We have previously shown that the PhhR transcriptional regulator induces phhAB genes encoding a phenylalanine hydroxylase. In this study we show, using microarray assays and promoter fusions, that PhhR is a global regulator responsible for the activation of genes essential for phenylalanine degradation, phenylalanine homeostasis and other genes of unknown function. Recently, it has been shown that phenylalanine catabolism occurs through more than one pathway. One of these possible pathways involves the metabolism of phenylalanine via tyrosine, p-hydroxyphenylpyruvate, and homogentisate. We identified two genes within this pathway that encode an acyl-CoA transferase involved in the metabolism of acetoacetate. All genes in this pathway were induced in response to phenylalanine in a PhhR-proficient background. The second potential degradative pathway involves the degradation of phenylalanine to produce phenylpyruvate, which seems to be degraded via phenylacetyl-CoA. A number of mutants in the paa genes encoding phenylacetyl-CoA degradation enzymes fail to grow on phenylpyruvate or phenylacetate, further supporting the existence of this second pathway. We found that the PhhR regulon also includes genes involved in the biosynthesis of aromatic amino acids that are repressed in the presence of phenylalanine, suggesting the possibility of feedback at the transcriptional level. In addition, we found that PhhR modulates the level of expression of the broad-substrate-specificity MexEF/OprN efflux pump. Expression from this pump is under the control of mexT gene product because phenylalanine-dependent transcription from the mexE promoter does not occur in a mexT mutant background. These results place PhhR as an important regulator in the control of bacterial responses to aromatic amino acids.

Redundancy of enzymes for formaldehyde detoxification in Pseudomonas putida.

Pseudomonas putida KT2440 exhibits redundant formaldehyde dehydrogenases and formate dehydrogenases that contribute to the detoxification of formaldehyde, a highly toxic compound. Physical and transcriptional analyses showed that the open reading frame (ORF) PP0328, encoding one of the formaldehyde dehydrogenases, is self-sufficient, whereas the other functional formaldehyde dehydrogenase gene (ORF PP3970) forms an operon with another gene of unknown function. Two formate dehydrogenase gene clusters (PP0489 to PP0492 and PP2183 to PP2186) were identified, and genes in these clusters were found to form operons. All four transcriptional promoters were mapped by primer extension and revealed the presence of noncanonical promoters expressed at basal level in the exponential growth phase and at a higher level in the stationary phase regardless of the presence of extracellular formaldehyde or formate. These promoters were characterized by a 5'-AG-CCA-C/A-CT-3' conserved region between -7 and -16. To determine the contribution of the different gene products to formaldehyde and formate mineralization, mutants with single and double mutations of formaldehyde dehydrogenases were generated, and the effect of the mutations on formaldehyde catabolism was tested by measuring (14)CO(2) evolution from (14)C-labeled formaldehyde. The results showed that both enzymes contributed to formaldehyde catabolism. A double mutant lacking these two enzymes still evolved CO(2) from formaldehyde, suggesting the presence of one or more still-unidentified formaldehyde dehydrogenases. Mutants with single and double mutations in the clusters for formate dehydrogenases were also generated, and all of them were able to metabolize [(14)C]formate to (14)CO(2), suggesting a redundancy of functions that was not limited to only the annotated genes. Single and double mutants deficient in formaldehyde dehydrogenases and formate dehydrogenases exhibited longer lag phases than did the parental strain when confronted with concentrations of formaldehyde close to the MICs. This suggests a role for the detoxification system in tolerance to sublethal concentrations of formaldehyde.

Pseudomonas savastanoi pv. savastanoi contains two iaaL paralogs, one of which exhibits a variable number of a trinucleotide (TAC) tandem repeat.

In this study, Pseudomonas savastanoi pv. savastanoi isolates were demonstrated to contain two iaaL paralogs, which are both chromosomally located in most strains. Comparative analysis of iaaL nucleotide sequences amplified from these two paralogs revealed that one paralog, iaaL(Psn), is 100% identical to iaaL from P. savastanoi pv. nerii, while the other paralog, iaaL(Psv), exhibited 93% identity to iaaL from Pseudomonas syringae pv. tomato (iaaL(Pto)). A 3-nucleotide motif (TAC) comprised of 3 to 15 repeats, which remained stable after propagation of the strains in olive plants, was found in iaaL(Psv). Based on the observed nucleotide sequence variations, a restriction fragment length polymorphism assay was developed that allowed differentiation among iaaL(Psn), iaaL(Psv), and iaaL(Pto)(.) In addition, reverse transcriptase PCR on total RNA from P. savastanoi pv. savastanoi strains demonstrated that both iaaL(Psv) and iaaL(Psn) containing 14 or fewer TAC repeats are transcribed. Capillary electrophoresis analysis of PCR-amplified DNA fragments containing the TAC repeats from iaaL(Psv) allowed the differentiation of P. savastanoi pv. savastanoi isolates.

Genomic analysis reveals the major driving forces of bacterial life in the rhizosphere.

BACKGROUND: Mutualistic interactions less well known than those between rhizobia and legumes are commonly found between plants and bacteria, frequently pseudomonads, which colonize roots and adjacent soil areas (the rhizosphere). RESULTS: A global analysis of Pseudomonas putida genes expressed during their interaction with maize roots revealed how a bacterial population adjusts its genetic program to this lifestyle. Differentially expressed genes were identified by comparing rhizosphere-colonizing populations with three distinct controls covering a variety of nutrients, growth phases and life styles (planktonic and sessile). Ninety rhizosphere up-regulated (rup) genes, which were induced relative to all three controls, were identified, whereas there was no repressed gene in common between the experiments. Genes involved in amino acid uptake and metabolism of aromatic compounds were preferentially expressed in the rhizosphere, which reflects the availability of particular nutrients in root exudates. The induction of efflux pumps and enzymes for glutathione metabolism indicates that adaptation to adverse conditions and stress (oxidative) response are crucial for bacterial life in this environment. The finding of a GGDEF/EAL domain response regulator among the induced genes suggests a role for the turnover of the secondary messenger c-diGMP in root colonization. Several mutants in rup genes showed reduced fitness in competitive root colonization. CONCLUSION: Our results show the importance of two selective forces of different nature to colonize the rhizosphere: stress adaptation and availability of particular nutrients. We also identify new traits conferring bacterial survival in this niche and open a way to the characterization of specific signalling and regulatory processes governing the plant-Pseudomonas association.

The ttgGHI solvent efflux pump operon of Pseudomonas putida DOT-T1E is located on a large self-transmissible plasmid.

Pseudomonas putida DOT-T1E is a solvent-tolerant strain able to grow in the presence of > 1% (v/v) toluene in the culture medium. A set of multidrug efflux pumps have been found to play a major role in the tolerance of this bacterium to organic solvents (Rojas et al., J Bacteriol 183: 3967-3973). In the course of studies of the mechanisms underlying solvent tolerance in DOT-T1E, we isolated a spontaneous solvent-sensitive mutant derivative which had lost the genes encoding the TtgGHI efflux pump, the most important extrusion element in quantitative terms. Genomic comparisons between the mutant and its parental strain by microarray analysis revealed that in addition to the ttgVW-ttgGHI gene cluster, another group of genes, highly similar to those found in the Tn4653A and ISPpu12 transposable elements of the TOL plasmid pWW0 from P. putida mt-2, were also absent from this strain. Further analysis demonstrated that strain DOT-T1E harboured a large plasmid (named pGRT1) that was lost from the solvent-sensitive mutant. Mapping analysis revealed that the ttgVW-ttgGHI genes and the Tn4653A-like transposon are borne by the pGRT1 plasmid. Plasmid pGRT1 is highly stable and its frequency of loss is below 10(-8) per cell per generation under a variety of growth conditions, including nutritional and physical stresses. The pGRT1 plasmid is self-transmissible, and its acquisition by the toluene-sensitive P. putida KT2440 and Pseudomonas aeruginosa PAO1 increased the recipient's tolerance to toluene up to levels similar to those exhibited by P. putida DOT-T1E. We discuss the importance and potential benefits of this plasmid for the development of bacteria with enhanced solvent tolerance, and its potential impact for bioremediation and whole-cell biotransformations.

Convergent peripheral pathways catalyze initial glucose catabolism in Pseudomonas putida: genomic and flux analysis.

In this study, we show that glucose catabolism in Pseudomonas putida occurs through the simultaneous operation of three pathways that converge at the level of 6-phosphogluconate, which is metabolized by the Edd and Eda Entner/Doudoroff enzymes to central metabolites. When glucose enters the periplasmic space through specific OprB porins, it can either be internalized into the cytoplasm or be oxidized to gluconate. Glucose is transported to the cytoplasm in a process mediated by an ABC uptake system encoded by open reading frames PP1015 to PP1018 and is then phosphorylated by glucokinase (encoded by the glk gene) and converted by glucose-6-phosphate dehydrogenase (encoded by the zwf genes) to 6-phosphogluconate. Gluconate in the periplasm can be transported into the cytoplasm and subsequently phosphorylated by gluconokinase to 6-phosphogluconate or oxidized to 2-ketogluconate, which is transported to the cytoplasm, and subsequently phosphorylated and reduced to 6-phosphogluconate. In the wild-type strain, glucose was consumed at a rate of around 6 mmol g(-1) h(-1), which allowed a growth rate of 0.58 h(-1) and a biomass yield of 0.44 g/g carbon used. Flux analysis of (13)C-labeled glucose revealed that, in the Krebs cycle, most of the oxalacetate fraction was produced by the pyruvate shunt rather than by the direct oxidation of malate by malate dehydrogenase. Enzymatic and microarray assays revealed that the enzymes, regulators, and transport systems of the three peripheral glucose pathways were induced in response to glucose in the outer medium. We generated a series of isogenic mutants in one or more of the steps of all three pathways and found that, although all three functioned simultaneously, the glucokinase pathway and the 2-ketogluconate loop were quantitatively more important than the direct phosphorylation of gluconate. In physical terms, glucose catabolism genes were organized in a series of clusters scattered along the chromosome. Within each of the clusters, genes encoding porins, transporters, enzymes, and regulators formed operons, suggesting that genes in each cluster coevolved. The glk gene encoding glucokinase was located in an operon with the edd gene, whereas the zwf-1 gene, encoding glucose-6-phosphate dehydrogenase, formed an operon with the eda gene. Therefore, the enzymes of the glucokinase pathway and those of the Entner-Doudoroff pathway are physically linked and induced simultaneously. It can therefore be concluded that the glucokinase pathway is a sine qua non condition for P. putida to grow with glucose.

Role of Pseudomonas putida tol-oprL gene products in uptake of solutes through the cytoplasmic membrane.

Proteins of the Tol-Pal (Tol-OprL) system play a key role in the maintenance of outer membrane integrity and cell morphology in gram-negative bacteria. Here we describe an additional role for this system in the transport of various carbon sources across the cytoplasmic membrane. Growth of Pseudomonas putida tol-oprL mutant strains in minimal medium with glycerol, fructose, or arginine was impaired, and the growth rate with succinate, proline, or sucrose as the carbon source was lower than the growth rate of the parental strain. Assays with radiolabeled substrates revealed that the rates of uptake of these compounds by mutant cells were lower than the rates of uptake by the wild-type strain. The pattern and amount of outer membrane protein in the P. putida tol-oprL mutants were not changed, suggesting that the transport defect was not in the outer membrane. Consistently, the uptake of radiolabeled glucose and glycerol in spheroplasts was defective in the P. putida tol-oprL mutant strains, suggesting that there was a defect at the cytoplasmic membrane level. Generation of a proton motive force appeared to be unaffected in these mutants. To rule out the possibility that the uptake defect was due to a lack of specific transporter proteins, the PutP symporter was overproduced, but this overproduction did not enhance proline uptake in the tol-oprL mutants. These results suggest that the Tol-OprL system is necessary for appropriate functioning of certain uptake systems at the level of the cytoplasmic membrane.

Transcriptional organization of the Pseudomonas putida tol-oprL genes.

Proteins of the Tol system play a key role in the maintenance of outer membrane integrity and cell morphology in gram-negative bacteria. In Pseudomonas putida, the seven genes, orf1, tolQ, tolR, tolA, tolB, oprL, and orf2, which encode the proteins of this complex, are clustered in a 5.8-kb region of chromosomal DNA. Analysis of polar mutations, reverse transcriptase PCR assays, and transcriptional fusion constructs with a promoterless lacZ gene revealed that the genes are arranged in two operons: orf1 tolQ tolR tolA tolB and oprL orf2. We were also able to find a transcript that was initiated at the orf1 promoter and covered the two operons in a single mRNA. On the basis of the OprL protein level, we surmised that this transcript contributed only about 10 to 15% of the total OprL protein. Primer extension analysis identified the oprL orf2 operon promoter within the tolB gene, and the -10 and -35 regions exhibited some similarity to those of sigma(70)-recognized promoters. The transcription start point of orf1 was located 91 bp upstream of the orf1 start codon, and the -10/-35 region also exhibited sigma(70) -10/-35 recognition sequences. The expression from both promoters in rich and minimal media was constitutive and was very little influenced by the growth phase or iron-deficient conditions. In addition, analyses of the beta-galactosidase activities of different translational fusion constructs revealed that translation of tolA and orf2 genes was dependent on the translation of their corresponding upstream genes (tolR and oprL, respectively).