Análisis genómico y funcional de los sistemas de Quorum Sensing en Rhizobium leguminosarum

Rhizobium leguminosarum (Rl) es una alfa-proteobacteria capaz de establecer una simbiosis diazotrófica con distintas leguminosas. A pesar de la importancia de esta simbiosis en el balance global del ciclo del nitrógeno, muy pocos genomas de rhizobios han sido secuenciados, que aporten nuevos conocimientos relacionados con las características genéticas que contribuyen a importantes procesos simbióticos. Únicamente tres secuencias completas de Rl han sido publicadas: Rl bv. viciae 3841 y dos genomas de Rl bv. trifolii (WSM1325 y WSM2304), ambos simbiontes de trébol. La secuencia genómica de Rlv UPM791 se ha determinado por medio de secuenciación 454. Este genoma tiene un tamaño aproximado de 7.8 Mb, organizado en un cromosoma y 5 replicones extracromosómicos, que incluyen un plásmido simbiótico de 405 kb. Este nuevo genoma se ha analizado en relación a las funciones simbióticas y adaptativas en comparación con los genomas completos de Rlv 3841 y Rl bv. trifolii WSM1325 y WSM2304. Mientras que los plásmidos pUPM791a y b se encuentran conservados, el plásmido simbiótico pUPM791c exhibe un grado de conservación muy bajo comparado con aquellos descritos en las otras cepas de Rl. Uno de los factores implicados en el establecimiento de la simbiosis es el sistema de comunicación intercelular conocido como Quorum Sensing (QS). El análisis del genoma de Rlv UPM791 ha permitido la identificación de dos sistemas tipo LuxRI mediados por señales de tipo N-acyl-homoserina lactonas (AHLs). El análisis mediante HPLC-MS ha permitido asociar las señales C6-HSL, C7-HSL y C8-HSL al sistema rhiRI, codificado en el plásmido simbiótico; mientras que el sistema cinRI, localizado en el cromosoma, produce 3OH-C14:1-HSL. Se ha identificado una tercera sintasa (TraI) codificada en el plásmido simbiótico, pero su regulador correspondiente se encuentra truncado debido a un salto de fase. Adicionalmente, se han encontrado tres reguladores de tipo LuxR-orphan que no presentan una sintasa LuxI asociada. El efecto potencial de las señales tipo AHL se ha estudiado mediante una estrategia de quorum quenching, la cual interfiere con los sistemas de QS de la bacteria. Esta estrategia está basada en la introducción del gen aiiA de Bacillus subtilis, que expresa constitutivamente una enzima lactonasa degradadora de AHLs. Para llevar a cabo el análisis en condiciones simbióticas, se ha desarrollado un sistema de doble marcaje que permite la identificación basado en los marcadores gusA y celB, que codifican para una enzima β–glucuronidasa y una β–galactosidasa termoestable, respectivamente. Los resultados obtenidos indican que Rlv UPM791 predomina sobre la cepa Rlv 3841 para la formación de nódulos en plantas de guisante. La baja estabilidad del plásmido que codifica para aiiA, no ha permitido obtener una conclusión definitiva sobre el efecto de la lactonasa AiiA en competitividad. Con el fin de analizar el significado y la regulación de la producción de moléculas señal tipo AHL, se han generado mutantes defectivos en cada uno de los dos sistemas de QS. Se ha llevado a cabo un análisis detallado sobre la producción de AHLs, formación de biofilm y simbiosis con plantas de guisante, veza y lenteja. El efecto de las deleciones de los genes rhiI y rhiR en Rlv UPM791 es más drástico en ausencia del plásmido pUPM791d. Mutaciones en cinI o cinRIS muestran tanto ausencia de señales, como producción exclusivamente de las de bajo peso molecular, respectivamente, producidas por el sistema rhiRI. Estas mutaciones mostraron un efecto importante en simbiosis. El sistema rhiRI se necesita para un comportamiento simbiótico normal. Además, mutantes cinRIS generaron nódulos blancos e ineficientes, mientras que el mutante cinI fue incapaz de producir nódulos en ninguna de las leguminosas utilizadas. Dicha mutación resultó en la inestabilización del plasmido simbiótico por un mecanismo dependiente de cinI que no ha sido aclarado. En general, los resultados obtenidos indican la existencia de un modelo de regulación dependiente de QS significativamente distinto a los que se han descrito previamente en otras cepas de R. leguminosarum, en las cuales no se había observado ningún fenotipo relevante en simbiosis. La regulación de la producción de AHLs Rlv UPM791 es un proceso complejo que implica genes situados en los plásmidos UPM791c y UPM791d, además de la señal 3-OH-C14:1-HSL. Finalmente, se ha identificado un transportador de tipo RND, homologo a mexAB-oprM de P. aeruginosa e implicado en la extrusión de AHLs de cadena larga. La mutación he dicho transportador no tuvo efectos apreciables sobre la simbiosis. ABSTRACT Rhizobium leguminosarum (Rl) is a soil alpha-proteobacterium that establishes a diazotrophic symbiosis with different legumes. Despite the importance of this symbiosis to the global nitrogen cycling balance, very few rhizobial genomes have been sequenced so far which provide new insights into the genetic features contributing to symbiotically relevant processes. Only three complete sequences of Rl strains have been published: Rl bv. viciae 3841, harboring six plasmids (7.75 Mb) and two Rl bv. trifolii (WSM1325 and WSM2304), both clover symbionts, harboring 5 and 4 plasmids, respectively (7.41 and 6.87 Mb). The genomic sequence of Rlv UPM791 was undertaken by means of 454 sequencing. Illumina and Sanger reads were used to improve the assembly, leading to 17 final contigs. This genome has an estimated size of 7.8 Mb organized in one chromosome and five extrachromosomal replicons, including a 405 kb symbiotic plasmid. Four of these plasmids are already closed, whereas there are still gaps in the smallest one (pUPM791d) due to the presence of insertion elements and repeated sequences, which difficult the assembly. The annotation has been carried out thanks to the Manatee pipeline. This new genome has been analyzed as regarding symbiotic and adaptive functions in comparison to the Rlv 3841 complete genome, and to those from Rl bv. trifolii strains WSM1325 and WSM2304. While plasmids pUPM791a and b are conserved, the symbiotic plasmid pUPM791c exhibited the lowest degree of conservation as compared to those from the other Rl strains. One of the factors involved in the symbiotic process is the intercellular communication system known as Quorum Sensing (QS). This mechanism allows bacteria to carry out diverse biological processes in a coordinate way through the production and detection of extracellular signals that regulate the transcription of different target genes. Analysis of the Rlv UPM791 genome allowed the identification of two LuxRI-like systems mediated by N-acyl-homoserine lactones (AHLs). HPLC-MS analysis allowed the adscription of C6-HSL, C7-HSL and C8-HSL signals to the rhiRI system, encoded in the symbiotic plasmid, whereas the cinRI system, located in the chromosome, produces 3OH-C14:1-HSL, previously described as “bacteriocin small”. A third synthase (TraI) is encoded also in the symbiotic plasmid, but its cognate regulator TraR is not functional due to a fameshift mutation. Three additional LuxR orphans were also found which no associated LuxI-type synthase. The potential effect of AHLs has been studied by means of a quorum quenching approach to interfere with the QS systems of the bacteria. This approach is based upon the introduction into the strains Rl UPM791 and Rl 3841 of the Bacillus subtilis gene aiiA expressing constitutively an AHL-degrading lactonase enzyme which led to virtual absence of AHL even when AiiA-expressing cells were a fraction of the total population. No significant effect of AiiA-mediated AHL removal on competitiveness for growth in solid surface was observed. For analysis under symbiotic conditions we have set up a two-label system to identify nodules produced by two different strains in pea roots, based on the markers gusA and celB, encoding a β–glucuronidase and a thermostable β–galactosidase enzymes, respectively. The results obtained show that Rlv UPM791 outcompetes Rlv 3841 for nodule formation in pea plants, and that the presence of the AiiA plasmid does not significantly affect the relative competitiveness of the two Rlv strains. However, the low stability of the pME6863 plasmid, encoding aiiA, did not lead to a clear conclusion about the AiiA lactonase effect on competitiveness. In order to further analyze the significance and regulation of the production of AHL signal molecules, mutants deficient in each of the two QS systems were constructed. A detailed analysis of the effect of these mutations on AHL production, biofilm formation and symbiosis with pea, vetch and lentil plants has been carried out. The effect of deletions on Rlv UPM791 rhiI and rhiR genes is more pronounced in the absence of plasmid pUPM791d, as no signal is detected in UPM791.1, lacking this plasmid. Mutations in cinI or cinRIS show either no signals, or only the small ones produced by the rhiRI system, suggesting that cinR might be regulating the rhiRI system. These mutations had a strong effect on symbiosis. Analysis of rhi mutants revealed that rhiRI system is required for normal symbiotic performance, as a drastic reduction of symbiotic fitness is observed when rhiI is deleted, and rhiR is essential for nitrogen fixation in the absence of plasmid pUPM791d. Furthermore, cinRIS mutants resulted in white and inefficient nodules, whereas cinI mutant was unable to form nodules on any legume tested. The latter mutation is associated to the instabilization of the symbiotic plasmid through a mechanism still uncovered. Overall, the results obtained indicate the existence of a model of QS-dependent regulation significantly different to that previously described in other R. leguminosarum strains, where no relevant symbiotic phenotype had been observed. The regulation of AHL production in Rlv UPM791 is a complex process involving the symbiotic plasmid (pUPM791c) and the smallest plasmid (pUPM791d), with a key role for the 3-OH-C14:1-HSL signal. Finally, we made a search for potential AHL transporters in Rlv UPM791 genome. These signals diffuse freely across membranes, but in the case of the long-chain AHLs an active efflux system might be required, as it has been described for C12-HSL in the case of Pseudomonas aeruginosa. We have identified a putative AHL transporter of the RND family homologous to P. aeruginosa mexAB-oprM. A mutant strain deficient in this transporter has been generated, and TLC analysis shows absence of 3OH-C14:1-HSL in its supernatant. This deficiency was complemented by the reintroduction of an intact copy of the genes via plasmid transfer. The mutation in mexAB genes had no significant effects on the symbiotic performance of R. leguminosarum bv. viciae.

Iron distribution through the developmental stages of Medicago truncatula nodules.

Paramount to symbiotic nitrogen fixation (SNF) is the synthesis of a number of metalloenzymes that use iron as a critical component of their catalytical core. Since this process is carried out by endosymbiotic rhizobia living in legume root nodules, the mechanisms involved in iron delivery to the rhizobia-containing cells are critical for SNF. In order to gain insight into iron transport to the nodule, we have used synchrotron-based X-ray fluorescence to determine the spatio-temporal distribution of this metal in nodules of the legume Medicago truncatula with hitherto unattained sensitivity and resolution. The data support a model in which iron is released from the vasculature into the apoplast of the infection/differentiation zone of the nodule (zone II). The infected cell subsequently takes up this apoplastic iron and delivers it to the symbiosome and the secretory system to synthesize ferroproteins. Upon senescence, iron is relocated to the vasculature to be reused by the shoot. These observations highlight the important role of yet to be discovered metal transporters in iron compartmentalization in the nodule and in the recovery of an essential and scarce nutrient for flowering and seed production.

MtNramp1 mediates iron import in rhizobia-infected Medicago truncatula cells.

Symbiotic nitrogen fixation is a process that requires relatively high quantities of iron provided by the host legume. Using synchrotron-based X-ray fluorescence, we have determined that this iron is released from the vasculature into the apoplast of zone II of M. truncatula nodules. This overlaps with the distribution of MtNramp1, a plasma membrane iron importer. The importance of MtNramp1 in iron transport for nitrogen fixation is indicated by the 60% reduction of nitrogenase activity observed in knock-down lines, most likely due to deficient incorporation of this essential metal cofactor at the necessary levels.

Functional characterization of Rhizobium leguminosarum bv. viciae DmeRF, a cation diffusion facilitator system involved in nickel and cobalt resistance

In prokaryotes, nickel is an essential element participating in the structure of enzymes involved in multiple cellular processes. Nickel transport is a challenge for microorganisms since, although essential, high levels of this metal inside the cell are toxic. For this reason, bacteria have developed high-affinity nickel transporters as well as nickel-specific detoxification systems. Ultramafic soils, and soils contaminated with heavy metals are excellent sources of nickel resistant bacteria. Molecular analysis of strains isolated in the habitats has revealed novel genetic systems involved in adaptation to such hostile conditions.

Identification and functional characterization of Rhizobium leguminosarum bv. viciae genetic systems involved in nickel homeostasis.

A collection of Rhizobium leguminosarum bv. viciae strains isolated from ultramafic and contaminated soils in Italy and Germany, respectively, was analyzed for resistance to nickel and cobalt ions. These assays led to the identification of strain UPM1137, which is able to grow at high concentrations of nickel and cobalt. In order to identify genetic systems involved in the homeostasis to these metals, a random mutagenesis was carried out in UPM1137 by inserting a Tn5-derivative minitransposon. As a result 4313 transconjugants were obtained, being 39 of them (0.90%) unable to grow at 1.5 mM NiCl2. The identification of the transposon insertion site in these mutants showed that the disrupted genes encode proteins belonging to different functional categories, where the secreted and membrane proteins were the most numerous. The analysis of heavy metal resistance and phenotypes in symbiotic and free –living cells will define the contribution of these genes to metal homeostasis.

Identification and analysis of a nickel-inducible cation diffusion facilitator-NiCoT combined system in Rhizobium leguminosarum bv. viciae

Nickel, like other transition metals, can be toxic to cells even at moderate concentration (low microM range) by displacing essential metals from their native binding sites or by generating reactive oxygen species that cause oxidative DNA damage. For this reason, cells have evolved mechanisms to deal with excess nickel. Efflux systems include members of the Resistance-Nodulation-cell Division (RND) protein family, P-type ATPases, cation diffusion facilitators (CDF) and other resistance factors. Nickel-specific exporters have been characterized in Cupravidus metallidurans, Helicobacter pylori, Achromobacter xylosoxidans, Serratia marcenses and Escherichia coli.

Azotobacter vinelandii nitrogenase: “Kinetics of nif gene expression and insights into the roles of FdxN and NifQ in FeMo-co biosynthesis”

The Molybdenum-nitrogenase is responsible for most biological nitrogen fixation activity (BNF) in the biosphere. Due to its great agronomical importance, it has been the subject of profound genetic and biochemical studies. The Mo nitrogenase carries at its active site a unique iron-molybdenum cofactor (FeMoco) that consists of an inorganic 7 Fe, 1 Mo, 1 C, 9 S core coordinated to the organic acid homocitrate. Biosynthesis of FeMo-co occurs outside nitrogenase through a complex and highly regulated pathway involving proteins acting as molecular scaffolds, metallocluster carriers or enzymes that provide substrates in appropriate chemical forms. Specific expression regulatory factors tightly control the accumulation levels of all these other components. Insertion of FeMo-co into a P-cluster containing apo-NifDK polypeptide results in nitrogenase reconstitution. Investigation of FeMo-co biosynthesis has uncovered new radical chemistry reactions and new roles for Fe-S clusters in biology.

Metagenomic Anlaysis of microsymbiont selection by the legume plant host

Rhizobium leguminosarum bv.viciae is able to establish nitrogen-fixing symbioses with legumes of the genera Pisum, Lens, Lathyrus and Vicia. Classic studies using trap plants (Laguerre et al., Young et al.) provided evidence that different plant hosts are able to select different rhizobial genotypes among those available in a given soil. However, these studies were necessarily limited by the paucity of relevant biodiversity markers. We have now reappraised this problem with the help of genomic tools. A well-characterized agricultural soil (INRA Bretennieres) was used as source of rhizobia. Plants of Pisum sativum, Lens culinaris, Vicia sativa and V. faba were used as traps. Isolates from 100 nodules were pooled, and DNA from each pool was sequenced (BGI-Hong Kong; Illumina Hiseq 2000, 500 bp PE libraries, 100 bp reads, 12 Mreads). Reads were quality filtered (FastQC, Trimmomatic), mapped against reference R. leguminosarum genomes (Bowtie2, Samtools), and visualized (IGV). An important fraction of the filtered reads were not recruited by reference genomes, suggesting that plant isolates contain genes that are not present in the reference genomes. For this study, we focused on three conserved genomic regions: 16S-23S rDNA, atpD and nodDABC, and a Single Nucleotide Polymorphism (SNP) analysis was carried out with meta / multigenomes from each plant. Although the level of polymorphism varied (lowest in the rRNA region), polymorphic sites could be identified that define the specific soil population vs. reference genomes. More importantly, a plant-specific SNP distribution was observed. This could be confirmed with many other regions extracted from the reference genomes (data not shown). Our results confirm at the genomic level previous observations regarding plant selection of specific genotypes. We expect that further, ongoing comparative studies on differential meta / multigenomic sequences will identify specific gene components of the plant-selected genotypes

Population genomics of host specificity in Rhizobium leguminosarum bv. viciae

Legumes establish a root-nodule symbiosis with soil bacteria collectively known as rhizobia. This symbiosis allows legumes to benefit from the nitrogen fixation capabilities of rhizobia and thus to grow in the absence of any fixed nitrogen source. This is especially relevant for Agriculture, where intensive plant growth depletes soils of useable, fixed nitrogen sources. One of the main features of the root nodule symbiosis is its specificity. Different rhizobia are able to nodulate different legumes. Rhizobium leguminosarum bv. viciae is able to establish an effective symbiosis with four different plant genera (Pisum, Lens, Vicia, Lathyrus), and any given isolate will nodulate any of the four plant genera. A population genomics study with rhizobia isolated from P. sativum, L. culinaris, V. sativa or V. faba, all originating in the same soil, showed that plants select specific genotypes from those available in that soil. This was demonstrated at the genome-wide level, but also for specific genes. Accelerated mesocosm studies with successive plant cultures provided additional evidence on this plant selection and on the nature of the genotypes selected. Finally, representatives from the major rhizobial genotypes isolated from these plants allowed characterization of the size and nature of the respective pangenome and specific genome compartments. These were compared to the different genotypes ?symbiotic and non-symbiotic?present in rhizobial populations isolated directly from the soil without plant intervention.

Medicago truncatula Natural Resistance Associated Macrophage Protein1 is required for iron uptake by rhizobia infected nodule cells

Iron is critical for symbiotic nitrogen fixation (SNF) as a key component ofmultiple ferroproteins involved in this biological process. In the model legume Medicago truncatula, iron is delivered by the vasculature to the infection/maturation zone (zone II) of the nodule, where it is released to the apoplast. From there, plasma membrane iron transporters move it into rhizobia-containing cells, where iron is used as the cofactor of multiple plant and rhizobial proteins (e.g. plant leghemoglobin and bacterial nitrogenase). MtNramp1 (Medtr3g088460) is the M. truncatula Natural Resistance-Associated Macrophage Protein family member, with the highest expression levels in roots and nodules. Immunolocalization studies indicate that MtNramp1 is mainly targeted to the plasma membrane. A loss-of-function nramp1 mutant exhibited reduced growth compared with the wild type under symbiotic conditions, but not when fertilized with mineral nitrogen. Nitrogenase activity was low in the mutant, whereas exogenous iron and expression of wild-type MtNramp1 in mutant nodules increased nitrogen fixation to normal levels. These data are consistent with a model in which MtNramp1 is the main transporter responsible for apoplastic iron uptake by rhizobia-infected cells in zone II.

Symbiosis-specific expression of Rhizobium etli casA encoding a secreted calmodulin-related protein

Symbiosis between Rhizobium and its leguminous host requires elaborate communication between the partners throughout the interaction process. A calmodulin-like protein, termed calsymin, was identified in Rhizobium etli; a calmodulin-related protein in a Gram-negative bacterium had not been described previously. Calsymin possesses three repeated homologous domains. Each domain contains two predicted EF-hand Ca2+-binding motifs. Ca2+-binding activity of calsymin was demonstrated on purified protein. R. etli efficiently secretes calsymin without N-terminal cleavage of the protein. The gene encoding calsymin, casA, is exclusively expressed during colonization and infection of R. etli with the host. Expression of casA is controlled by a repressor protein, termed CasR, belonging to the TetR family of regulatory proteins. Mutation of the casA gene affects the development of bacteroids during symbiosis and symbiotic nitrogen fixation.