The Complete Genome of Rhizobium favelukesii LPU83T: Insights into Plastic pSym and Its Symbiotic Incompatibility with a Broad Range of Legume Hosts

Achieving completeness of multipartite bacterial genomes has been a difficult task, especially in rhizobia. In this study, we performed a deep bioinformatic analysis of the newly re-sequenced genome of Rhizobium favelukesii LPU83T. This strain was isolated from acid soils in Argentina and is capable of nodulating several leguminous plants, although it is unable to fix nitrogen efficiently in any of them. Oxford Nanopore sequencing allowed us to completely assemble the symbiotic plasmid of the strain, pRfaLPU83b, and we discovered that it harbors three intact prophages and a high density of insertion sequences (ISs). These characteristics show why it is often so difficult to complete the symbiotic plasmids of rhizobial strains and the importance of having long-read sequencing methods. Upon detailed analysis of this replicon, we identified a complete conjugation system with gene structure consistent with quorum sensing-associated systems that may have contributed to the genetic mosaic structure of the strain. Furthermore, we identified in the symbiotic plasmid of R. favelukesii LPU83T a large proportion of the symbiotic genes previously identified as essential for Biological Nitrogen Fixation (BNF) in symbiosis with alfalfa, with a high percentage of identity with respect to those of Sinorhizobium meliloti 2011. Among the determinants related to BNF, we found genes encoding the HrrP and SapA peptidases in the LPU83 genome, previously described and related to the degradation of nodule-specific cysteine-rich peptides. These peptides are essential for bacteroid differentiation and, therefore, efficient BNF. Our results show that despite having these genes, they are not directly responsible for the inefficient BNF phenotype of LPU83.