Protein domains and architectural innovation in plant-associated Proteobacteria

BACKGROUND: Evolution of new complex biological behaviour tends to arise by novel combinations of existing building blocks. The functional and evolutionary building blocks of the proteome are protein domains, the function of a protein being dependent on its constituent domains. We clustered completely-sequenced proteomes of prokaryotes on the basis of their protein domain content, as defined by Pfam (release 16.0). This revealed that, although there was a correlation between phylogeny and domain content, other factors also have an influence. This observation motivated an investigation of the relationship between an organism's lifestyle and the complement of domains and domain architectures found within its proteome. RESULTS: We took a census of all protein domains and domain combinations (architectures) encoded in the completely-sequenced proteobacterial genomes. Nine protein domain families were identified that are found in phylogenetically disparate plant-associated bacteria but are absent from non-plant-associated bacteria. Most of these are known to play a role in the plant-associated lifestyle, but they also included domain of unknown function DUF1427, which is found in plant symbionts and pathogens of the alpha-, beta- and gamma-Proteobacteria, but not known in any other organism. Further, several domains were identified as being restricted to phytobacteria and Eukaryotes. One example is the RolB/RolC glucosidase family, which is found only in Agrobacterium species and in plants. We identified the 0.5% of Pfam protein domain families that were most significantly over-represented in the plant-associated Proteobacteria with respect to the background frequencies in the whole set of available proteobacterial proteomes. These included guanylate cyclase, domains implicated in aromatic catabolism, cellulase and several domains of unknown function. We identified 459 unique domain architectures found in phylogenetically diverse plant pathogens and symbionts that were absent from non-pathogenic and non-symbiotic relatives. The vast majority of these were restricted to a single species or several closely related species and so their distributions could be better explained by phylogeny than by lifestyle. However, several architectures were found in two or more very distantly related phytobacteria but absent from non-plant-associated bacteria. Many of the proteins with these unique architectures are predicted to be secreted. In Pseudomonas syringae pathovar tomato, those genes encoding genes with novel domain architectures tended to have atypical GC contents and were adjacent to insertion sequence elements and phage-like sequences, suggesting acquisition by horizontal transfer. CONCLUSIONS: By identifying domains and architectures unique to plant pathogens and symbionts, we highlighted candidate proteins for involvement in plant-associated bacterial lifestyles. Given that characterisation of novel gene products in vivo and in vitro is time-consuming and expensive, this computational approach may be useful for reducing experimental search space. Furthermore we discuss the biological significance of novel proteins highlighted by this study in the context of plant-associated lifestyles

Adaptation of Rhizobium leguminosarum to pea, alfalfa and sugar beet rhizospheres investigated by comparative transcriptomics

BACKGROUND: The rhizosphere is the microbe-rich zone around plant roots and is a key determinant of the biosphere's productivity. Comparative transcriptomics was used to investigate general and plant-specific adaptations during rhizosphere colonization. Rhizobium leguminosarum biovar viciae was grown in the rhizospheres of pea (its legume nodulation host), alfalfa (a non-host legume) and sugar beet (non-legume). Gene expression data were compared to metabolic and transportome maps to understand adaptation to the rhizosphere. RESULTS: Carbon metabolism was dominated by organic acids, with a strong bias towards aromatic amino acids, C1 and C2 compounds. This was confirmed by induction of the glyoxylate cycle required for C2 metabolism and gluconeogenesis in all rhizospheres. Gluconeogenesis is repressed in R. leguminosarum by sugars, suggesting that although numerous sugar and putative complex carbohydrate transport systems are induced in the rhizosphere, they are less important carbon sources than organic acids. A common core of rhizosphere-induced genes was identified, of which 66% are of unknown function. Many genes were induced in the rhizosphere of the legumes, but not sugar beet, and several were plant specific. The plasmid pRL8 can be considered pea rhizosphere specific, enabling adaptation of R. leguminosarum to its host. Mutation of many of the up-regulated genes reduced competitiveness for pea rhizosphere colonization, while two genes specifically up-regulated in the pea rhizosphere reduced colonization of the pea but not alfalfa rhizosphere. CONCLUSIONS: Comparative transcriptome analysis has enabled differentiation between factors conserved across plants for rhizosphere colonization as well as identification of exquisite specific adaptation to host plants

SCARN a Novel Class of SCAR Protein That Is Required for Root-Hair Infection during Legume Nodulation

Rhizobial infection of legume root hairs requires a rearrangement of the actin cytoskeleton to enable the establishment of plant-made infection structures called infection threads. In the SCAR/WAVE (Suppressor of cAMP receptor defect/WASP family verpolin homologous protein) actin regulatory complex, the conserved N-terminal domains of SCAR proteins interact with other components of the SCAR/WAVE complex. The conserved C-terminal domains of SCAR proteins bind to and activate the actin-related protein 2/3 (ARP2/3) complex, which can bind to actin filaments catalyzing new actin filament formation by nucleating actin branching. We have identified, SCARN (SCAR-Nodulation),a gene required for root hair infection of Lotus japonicus by Mesorhizobium loti. Although the SCARN protein is related to Arabidopsis thaliana SCAR2 and SCAR4, it belongs to a distinct legume-sub clade. We identified other SCARN-like proteins in legumes and phylogeny analyses suggested that SCARN may have arisen from a gene duplication and acquired specialized functions in root nodule symbiosis. Mutation of SCARN reduced formation of infection-threads and their extension into the root cortex and slightly reduced root-hair length. Surprisingly two of the scarn mutants showed constitutive branching of root hairs in uninoculated plants. However we observed no effect of scarn mutations on trichome development or on the early actin cytoskeletal accumulation that is normally seen in root hair tips shortly after M. loti inoculation, distinguishing them from other symbiosis mutations affecting actin nucleation. The C-terminal domain of SCARN binds to ARPC3 and ectopic expression of the N-terminal SCAR-homology domain (but not the full length protein) inhibited nodulation. In addition, we found that SCARN expression is enhanced by M. loti in epidermal cells and that this is directly regulated by the NODULE INCEPTION (NIN) transcription factor

MtLAX2, a functional homologue of the Arabidopsis auxin influx transporter AUX1, is required for nodule organogenesis

Most legume plants can form nodules, specialized lateral organs that form on roots, and house nitrogen-fixing bacteria collectively called rhizobia. The uptake of the phytohormone auxin into cells is known to be crucial for development of lateral roots. To test the role of auxin influx in nodulation we used the auxin influx inhibitors 1-naphthoxyacetic acid (1-NOA) and 2-NOA, which we found reduced nodulation of Medicago truncatula. This suggested the possible involvement of the AUX/LAX family of auxin influx transporters in nodulation. Gene expression studies identified MtLAX2, a paralogue of Arabidopsis (Arabidopsis thaliana) AUX1, as being induced at early stages of nodule development. MtLAX2 is expressed in nodule primordia, the vasculature of developing nodules, and at the apex of mature nodules. The MtLAX2 promoter contains several auxin response elements, and treatment with indole-acetic acid strongly induces MtLAX2 expression in roots. mtlax2 mutants displayed root phenotypes similar to Arabidopsis aux1 mutants, including altered root gravitropism, fewer lateral roots, shorter root hairs, and auxin resistance. In addition, the activity of the synthetic DR5-GUS auxin reporter was strongly reduced in mtlax2 roots. Following inoculation with rhizobia, mtlax2 roots developed fewer nodules, had decreased DR5-GUS activity associated with infection sites, and had decreased expression of the early auxin responsive gene ARF16a. Our data indicate that MtLAX2 is a functional analog of Arabidopsis AUX1 and is required for the accumulation of auxin during nodule formation in tissues underlying sites of rhizobial infection

Widely conserved AHL transcription factors are essential for NCR gene expression and nodule development in Medicago

Symbiotic nitrogen fixation by Rhizobium bacteria in the cells of legume root nodules alleviates the need for nitrogen fertilizers. Nitrogen fixation requires the endosymbionts to differentiate into bacteroids which can be reversible or terminal. The latter is controlled by the plant, it is more beneficial and has evolved in multiple clades of the Leguminosae family. The plant effectors of terminal differentiation in inverted repeat-lacking clade legumes (IRLC) are nodule-specific cysteine-rich (NCR) peptides, which are absent in legumes such as soybean where there is no terminal differentiation of rhizobia. It was assumed that NCR s co-evolved with specific transcription factors, but our work demonstrates that expression of NCR genes does not require NCR -specific transcription factors. Introduction of the Medicago truncatula NCR169 gene under its own promoter into soybean roots resulted in its nodule-specific expression, leading to bacteroid changes associated with terminal differentiation. We identified two AT-Hook Motif Nuclear Localized (AHL) transcription factors from both M. truncatula and soybean nodules that bound to AT-rich sequences in the NCR169 promoter inducing its expression. Whereas mutation of NCR169 arrested bacteroid development at a late stage, the absence of MtAHL1 or MtAHL2 completely blocked bacteroid differentiation indicating that they also regulate other NCR genes required for the development of nitrogen-fixing nodules. Regulation of NCR s by orthologous transcription factors in non-IRLC legumes opens up the possibility of increasing the efficiency of nitrogen fixation in legumes lacking NCR s

Morphotype of bacteroids in different legumes correlates with the number and type of symbiotic NCR peptides

In legume nodules, rhizobia differentiate into nitrogen-fixing forms called bacteroids, which are enclosed by a plant membrane in an organelle-like structure called the symbiosome. In the Inverted Repeat-Lacking Clade (IRLC) of legumes, this differentiation is terminal due to irreversible loss of cell division ability and is associated with genome amplification and different morphologies of the bacteroids that can be swollen, elongated, spherical, and elongatedbranched, depending on the host plant. In Medicago truncatula, this process is orchestrated by nodule-specific cysteine-rich peptides (NCRs) delivered into developing bacteroids. Here, we identified the predicted NCR proteins in 10 legumes representing different subclades of the IRLC with distinct bacteroid morphotypes. Analysis of their expression and predicted sequences establishes correlations between the composition of the NCR family and the morphotypes of bacteroids. Although NCRs have a single origin, their evolution has followed different routes in individual lineages, and enrichment and diversification of cationic peptides has resulted in the ability to impose major morphological changes on the endosymbionts. The wide range of effects provoked by NCRs such as cell enlargement, membrane alterations and permeabilization, and biofilm and vesicle formation is dependent on the amino acid composition and charge of the peptides. These effects are strongly influenced by the rhizobial surface polysaccharides that affect NCR-induced differentiation and survival of rhizobia in nodule cells

Mutation of praR in Rhizobium leguminosarum enhances root biofilms, improving nodulation competitiveness by increased expression of attachment proteins

Summary In Rhizobium leguminosarum bv. viciae, quorumsensing is regulated by CinR, which induces the cinIS operon. CinI synthesizes an AHL, whereas CinS inactivates PraR, a repressor. Mutation of praR enhanced biofilms in vitro. We developed a light (lux)-dependent assay of rhizobial attachment to roots and demonstrated that mutation of praR increased biofilms on pea roots. The praR mutant out-competed wild-type for infection of pea nodules in mixed inoculations. Analysis of gene expression by microarrays and promoter fusions revealed that PraR represses its own transcription and mutation of praR increased expression of several genes including those encoding secreted proteins (the adhesins RapA2, RapB and RapC, two cadherins and the glycanase PlyB), the polysaccharide regulator RosR, and another protein similar to PraR. PraR bound to the promoters of several of these genes indicating direct repression. Mutations in rapA2, rapB, rapC, plyB, the cadherins or rosR did not affect the enhanced root attachment or nodule competitiveness of the praR mutant. However combinations of mutations in rapA, rapB and rapC abolished the enhanced attachment and nodule competitiveness. We conclude that relief of PraR-mediated repression determines a lifestyle switch allowing the expression of genes that are important for biofilm formation on roots and the subsequent initiation of infection of legume roots

Arabinose and protocatechuate catabolism genes are important for growth of Rhizobium leguminosarum biovar viciae in the pea rhizosphere

Background and aims: To form nitrogen-fixing nodules on pea roots, Rhizobium leguminosarum biovar viciae must be competitive in the rhizosphere. Our aim was to identify genes important for rhizosphere fitness. Methods: Signature-tagged mutants were screened using microarrays to identify mutants reduced for growth in pea rhizospheres. Candidate mutants were assessed relative to controls for growth in minimal medium, growth in pea rhizospheres and for infection of peas in mixed inoculants. Mutated genes were identified by DNA sequencing and confirmed by transduction. Results: Of 5508 signature-tagged mutants, microarrays implicated 50 as having decreased rhizosphere fitness. Growth tests identified six mutants with rhizosphere-specific phenotypes. The mutation in one of the genes (araE) was in an arabinose catabolism operon and blocked growth on arabinose. The mutation in another gene (pcaM), encoding a predicted solute binding protein for protocatechuate and hydroxybenzoate uptake, decreased growth on protocatechuate. Both mutants were decreased for nodule infection competitiveness with mixed inoculants, but nodulated peas normally when inoculated alone. Other mutants with similar phenotypes had mutations predicted to affect secondary metabolism. Conclusions: Catabolism of arabinose and protocatechuate in the pea rhizosphere is important for competitiveness of R.l. viciae. Other genes predicted to be involved in secondary metabolism are also important

Nonlinear Time Series Analysis of Nodulation Factor Induced Calcium Oscillations: Evidence for Deterministic Chaos?

Legume plants form beneficial symbiotic interactions with nitrogen fixing bacteria (called rhizobia), with the rhizobia being accommodated in unique structures on the roots of the host plant. The legume/rhizobial symbiosis is responsible for a significant proportion of the global biologically available nitrogen. The initiation of this symbiosis is governed by a characteristic calcium oscillation within the plant root hair cells and this signal is activated by the rhizobia. Recent analyses on calcium time series data have suggested that stochastic effects have a large role to play in defining the nature of the oscillations. The use of multiple nonlinear time series techniques, however, suggests an alternative interpretation, namely deterministic chaos. We provide an extensive, nonlinear time series analysis on the nature of this calcium oscillation response. We build up evidence through a series of techniques that test for determinism, quantify linear and nonlinear components, and measure the local divergence of the system. Chaos is common in nature and it seems plausible that properties of chaotic dynamics might be exploited by biological systems to control processes within the cell. Systems possessing chaotic control mechanisms are more robust in the sense that the enhanced flexibility allows more rapid response to environmental changes with less energetic costs. The desired behaviour could be most efficiently targeted in this manner, supporting some intriguing speculations about nonlinear mechanisms in biological signaling