Native bradyrhizobial symbionts of Lupinus mariae-josephae, a unique endemism thriving in alkaline soils in Eastern Spain

Lupinus mariae-josephae (Lmj) es una especie de lupino endémica de una pequeña y específica área de Comunidad Valenciana (Este de España), donde prospera en suelos alcalinoscalcáreos, un hábitat singular para los altramuces, que crecen preferentemente en suelos ácidos o neutros. Esto hace de Lmj una especie de lupino única. Cuando se inició este trabajo, la extensión conocida de este endemismo abarcaba unos 700 kilómetros cuadrados, confinados en la provincia de Valencia. En esta área, Lmj prospera en pequeñas poblaciones aisladas que contienen un número reducido de plantas por lo que se la consideró una especie en peligro de extinción. Todos los esfuerzos, utilizando estrategias clásicas dirigidas a ampliar el área de crecimiento de Lmj y garantizar su conservación, han tenido un éxito limitado. El trabajo que se presenta está dirigido a mejorar el conocimiento de la ecología de Lmj, en particular la interacción simbiótica que establece con bacterias del suelo denominadas rizobios y se centra en la caracterización fenotípica, filogenética y genómica de esos rizobios. También se investiga la posible contribución de la simbiosis en mejorar la conservación de Lmj. Para este fin, se han estudiado diferentes aspectos que se describen a continuación. El primero objetivo se centró en aislar y estudiar de la diversidad genética de las bacterias endosimbióticas de Lmj. . Se realizó un análisis filogenético de genes esenciales que mostró que las cepas de Lmj pertenecen al género Bradyrhizobium y que presentan una gran diversidad con características fenotípicas y simbióticas diferentes de cepas de Bradyrhizobium que nodulan otras especies de lupinos nativos de España (cepas ISLU). Las cepas estudiadas se dividieron en dos grupos (Clado I y Clado II). El Clado I, incluye a las cepas Lmj, definiendo un nuevo linaje, filogenéticamente relacionado con otras especies de Bradyrhizobium, como B. jicamae y B. elkanii. El Clado II contiene cepas ISLU relacionadas con cepas de B. canariense y B. japonicum que establecen simbiosis con lupinos de suelos ácidos. Otro análisis filogenético basado en genes simbióticos, distribuyó las cepas de Lmj en sólo dos grupos diferentes. La singularidad y gran diversidad de estas cepas en una pequeña área geográfica, hacen de este, un atractivo sistema para el estudio de la evolución y adaptación de las bacterias simbióticas a su respectiva planta huésped. Adicionalmente, se estudio la presencia de bacterias capaces de nodular Lmj en suelos básicos de Chiapas, México. Sorprendentemente, estos suelos contienen bacterias capaces establecer interacciones simbióticas eficientes con Lmj en ensayos de invernadero. A continuación se investigó la taxonomía de los endosimbiontes de Lmj analizando la secuencia de cuatro genes esenciales (16S rRNA, recA, glnII y atpD) y el promedio de identidad de nucleótidos de genomas completos de algunas cepas representativas de la diversidad (ANIm). Se identificaron nuevas especies de Bradyrhizobium dentro del Clado I y se definió una de ellas: 'Bradyrhizobium valentinum' sp. nov (cepa tipo LmjM3T = CECT 8364T, LMG 2761T). También se abordó cómo conservar Lmj en su hábitat natural mediante inoculación con alguna de las cepas aisladas. Se demostró la ausencia de bacterias capaces de nodular Lmj en suelos rojos alcalinos o ‘‘terra rossa’’ de la Península Ibérica y Baleares. Dos cepas, altamente eficientes en cuanto a la fijación de nitrógeno, LmjC y LmjM3T, fueron seleccionadas para ser empleadas como inoculantes. Dos experimentos de campo llevados a cabo en años consecutivos en áreas con características edafoclimáticas similares a las que presentan las poblaciones de Lmj, lograron la reproducción exitosa de la planta. Se concluyó que un ciclo reproductivo exitoso de Lmj es absolutamente dependiente de la inoculación con sus simbiontes naturales y que la simbiosis debe ser considerada un factor esencial en estrategias de conservación de leguminosas en peligro. La obtención de varias secuencias genómicas de cepas aisladas de Lmj y de otras cepas de Bradyrhizobium reveló una alta similitud entre los genomas de las cepas del Clado I, y permitió la identificación de cinco posibles nuevas especies. Además, se estudiaron tres agrupaciones de genes relacionados con la simbiosis (nod, nif y fix) definiendo un nuevo linaje para las cepas de Lmj, diferente del symbiovar “genistearum” de B. canariense y B. japonicum. La baja diversidad encontrada en el análisis filogenético de los genes simbióticos contrasta con la gran diversidad asociada a genes esenciales. La presencia de plásmidos en cepas del género Bradyrhizobium ha sido descrita en muy pocas ocasiones, sin embargo el análisis de la secuencia genómica de la cepa ISLU101, aislada de Lupinus angustifolius, reveló la presencia de un origen de replicación extracromosómico homólogo al operón repABC, presente en el plásmido de Bradyrhizobium sp BTAi1. Gracias a esta secuencia se identificaron genes homólogos en 19 de 72 cepas ISLU. Filogenéticamente, las secuencias de repABC se agruparon en un grupo monofilético con las de pBTAi1 y separadas de los rizobios de crecimiento rápido. Finalmente, se identificaron sistemas de secreción de proteínas de tipo III (T3SS) en nueve genomas de cepas de Lmj. Los T3SS pueden inyectar proteínas efectoras al interior de células vegetales. Su presencia en rizobios se ha relacionado con la gama de hospedador que pueden nodular y puede tener un efecto beneficioso, neutro o perjudicial en la simbiosis. Los T3SS de las cepas de Lmj codifican para una proteína efectora similar a NopE, un efector dependiente de T3SS descrito en B. diazoefficiens USDA 110T. La proteína NopE de la cepa LmjC se ha caracterizado bioquímicamente. ABSTRACT Lupinus mariae-josephae (Lmj) is a lupine species endemic of a unique small area in Valencia region (Eastern Spain) where the lupine plants thrive in alkaline-limed soils, which preferentially grow in acid or neutral soils. This is the type of soils native lupines of Spain. When this work was initiated, the extension of the endemic area of Lmj was of about 700 squared kilometers confined to the Valencia province. In this area, Lmj thrives in small, isolated patches containing a reduced number of plants, and points to an endemism that can easily became endangered or extinct. Consequently, the Valencia Community authorities gave a ‘‘microreserve” status for conservation of the species. All efforts, using classical strategies directed to extend the area of Lmj growth and ensure its conservation have been so far unsuccessful. The work presented here is directed to improve our knowledge of Lmj ecology and it is centered in the characterization of the rhizobial symbiosis by phenotypic, phylogenetic and genomic analysis as well as in investigate the potential contribution of the symbiosis to improve its conservation. To this end, five different topics have been studied, and results are briefly described here. Extensive details can be followed en the attached, published articles. The first topic deals with the indigenous rhizobial symbionts of the Lmj endemism, and its genetic diversity was investigated. The Lmj root symbionts belong to the Bradyrhizobium genus, and phylogenetic analysis based on core genes identified a large diversity of Bradyrhizobium strains with phenotypic and symbiotic characteristics different from rhizobia nodulating other Lupinus spp. native of Spain. The strains were split in two clades. Clade II contained strains close to classical B. canariense and B. japonicum lineages that establish symbioses with lupines in acid soils of the Mediterranean area. Clade I included Lmj strains that define a new lineage, close to other Bradyrhizobium species as B. jicamae and B. elkanii. The phylogenetic analysis based on symbiotic genes identified only two distinct clusters. The singularity and large diversity of these strains in such a small geographical area makes this an attractive system for studying the evolution and adaptation of the rhizobial symbiont to the plant host. Additionally, the presence of bacteria able to nodulate Lmj in basic soils from Chiapas, Mexico was investigated. Surprisingly, these soils contain bacteria able to effectively nodulate and fix nitrogen with Lmj plants in greenhouse assays. In the second topic, the taxonomic status of the endosymbiotic bacteria of Lmj from Valencia endemism and Chiapas was investigated. Results from phylogenetic analysis of core genes and Average Nucleotide Identity (ANIm) using draft genomic sequences identified new Bradyrhizobium species within strains of Clade I of Lmj endosymbiotic bacteria. Only one of these potentially new species has been defined, meanwhile the others are under process of characterization. The name ‘Bradyrhizobium valentinum’ sp. nov. was proposed for the defined species (type strain LmjM3T= CECT 8364T, LMG 2761T). The third topic was directed to conservation of endangered Lmj in its natural habitat. The relevant conclusion of this experimentation is that the symbiosis should be considered as a relevant factor in the conservation strategies for endangered legumes. First, we showed absence of bacteria able to nodulate Lmj in all the inspected ‘‘terra rossa’’ or alkaline red soils of the Iberian Peninsula and Balearic Islands. Then, two efficient nitrogen fixing strains with Lmj plants, LmjC and LmjM3T, were selected as inoculum for seed coating. Two planting experiments were carried out in consecutive years under natural conditions in areas with edapho-climatic characteristics identical to those sustaining natural Lmj populations, and successful reproduction of the plant was achieved. The relevant conclusion from these assays was that the successful reproductive cycle was absolutely dependent on seedling inoculation with effective bradyrhizobia The forth topic deep into the analysis of the genomic of Lmj representative strains. To this end, draft genomic sequences of selected Lmj strains and type strains of Bradyrhizobium spp. were assembled. The comparison analysis of the draft genomic sequences of Lmj strains and related Bradyrhizobium species grouped in Clade I, revealed a high genomic homology among them, and allowed the definition of five potentially new species of Lmj nodulating bacteria. Also, based on the available draft genomic sequences, only three clusters of nod, fix and nif genes from Lmj strains were identified and showed to define a new symbiotic lineage, distant from that of B. canariense and B. japonicum bv. genistearum. The low diversity exhibited by the phylogenetic analysis of symbiotic genes contrast with the large diversity of strains as regards the housekeeping genes analyzed. Besides, the genomic analysis of a Lupinus angustifolius strain ISLU101, revealed the presence of an extrachromosomal replication origin homologous to repABC cluster from plasmid present in Bradyrhizobium spp BTAi1. This repABC cluster gene sequence allowed the identification of extrachromosomic replication origin in 19 out of 72 Bradyrhizobium strains from Lupinus spp., a highly significant result since the absence of plasmids in the Bradyrhizobium genus was traditionally assumed. The repABC gene sequences of these strains grouped them in a unique monophyletic group, related to B. sp. BTAi1 plasmid, but differentiated from the repABC gene cluster of plasmids in fast growing rhizobium strains. The last topic was focused on characterization of type III secreted effectors present in Lmj endosymbiotic bacteria. Type III secretion systems (T3SS) are specialized protein export machineries which can deliver effector proteins into plant cells. The presence of T3SS in rhizobia has frequently been related to the symbiotic nodulation host-range and may have a beneficial or detrimental effect on the symbiosis with legumes. In this context, the presence of T3SS in genomes of nine Lmj strains was investigated, and it was shown the presence of clusters encoding NopE type III-secreted protein similar to the NopE1 and NopE2 of B. diazoefficiens USDA 110T. The putative NopE protein of LmjC strain is at present being characterized regarding its structure and function.

Endosymbiotic bacteria nodulating a new endemic lupine Lupinus mariae-josephi from alkaline soils in Eastern Spain represent a new lineage within the Bradyrhizobium genus

Lupinus mariae-josephi is a recently described endemic Lupinus species from a small area in Eastern Spain where it thrives in soils with active lime and high pH. The L. mariae-josephi root symbionts were shown to be very slow-growing bacteria with different phenotypic and symbiotic characteristics from those of Bradyrhizobium strains nodulating other Lupinus. Their phylogenetic status was examined by multilocus sequence analyses of four housekeeping genes (16S rRNA, glnII, recA, and atpD) and showed the existence of a distinct evolutionary lineage for L. mariae-josephi that also included Bradyrhizobium jicamae. Within this lineage, the tested isolates clustered in three different sub-groups that might correspond to novel sister Bradyrhizobium species. These core gene analyses consistently showed that all the endosymbiotic bacteria isolated from other Lupinus species of the Iberian Peninsula were related to strains of the B. canariense or B. japonicum lineages and were separate from the L. mariae-josephi isolates. Phylogenetic analysis based on nodC symbiotic gene sequences showed that L. mariae-josephi bacteria also constituted a new symbiotic lineage distant from those previously defined in the genus Bradyrhizobium. In contrast, the nodC genes of isolates from other Lupinus spp. from the Iberian Peninsula were again clearly related to the B. canariense and B. japonicum bv. genistearum lineages. Speciation of L. mariae-josephi bradyrhizobia may result from the colonization of a singular habitat by their unique legume host.

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.