Next-generation annotation of prokaryotic genomes with EuGene-P: application to Sinorhizobium meliloti 2011

Comparative Proteomic Analysis Revealing ActJ-Regulated Proteins in Sinorhizobium meliloti

To adapt to different environmental conditions, Sinorhizobium meliloti relies on finely tuned regulatory networks, most of which are unexplored to date. We recently demonstrated that deletion of the two-component system ActJK renders an acid-vulnerable phenotype in S. meliloti and negatively impacts bacteroid development and nodule occupancy as well. To fully understand the role of ActJ in acid tolerance, S. meliloti wild-type and S. meliloti ΔactJ proteomes were compared in the presence or absence of acid stress by nanoflow ultrahigh-performance liquid chromatography coupled to mass spectrometry. The analysis demonstrated that proteins involved in the synthesis of exopolysaccharides (EPSs) were notably enriched in ΔactJ cells in acid pH. Total EPS quantification further revealed that although EPS production was augmented at pH 5.6 in both the ΔactJ and the parental strain, the lack of ActJ significantly enhanced this difference. Moreover, several efflux pumps were found to be downregulated in the ΔactJ strain. Promoter fusion assays suggested that ActJ positively modulated its own expression in an acid medium but not at under neutral conditions. The results presented here identify several ActJ-regulated genes in S. meliloti, highlighting key components associated with ActJK regulation that will contribute to a better understanding of rhizobia adaptation to acid stress.

Rhizobial migration toward roots mediated by FadL-ExoFQP modulation of extracellular long-chain AHLs

Migration from rhizosphere to rhizoplane is a key selecting process in root microbiome assembly, but not fully understood. Rhizobiales members are overrepresented in the core root microbiome of terrestrial plants, and here we report a genome-wide transposon-sequencing of rhizoplane fitness genes of beneficial Sinorhizobium fredii on wild soybean, cultivated soybean, rice, and maize. There were few genes involved in broad-host-range rhizoplane colonization. The fadL mutant lacking a fatty acid transporter exhibited high colonization rates, while mutations in exoFQP (encoding membrane proteins directing exopolysaccharide polymerization and secretion), but not those in exo genes essential for exopolysaccharide biosynthesis, led to severely impaired colonization rates. This variation was not explainable by their rhizosphere and rhizoplane survivability, and associated biofilm and exopolysaccharide production, but consistent with their migration ability toward rhizoplane, and associated surface motility and the mixture of quorum-sensing AHLs (N-acylated-L-homoserine lactones). Genetics and physiology evidences suggested that FadL mediated long-chain AHL uptake while ExoF mediated the secretion of short-chain AHLs which negatively affected long-chain AHL biosynthesis. The fadL and exoF mutants had elevated and depleted extracellular long-chain AHLs, respectively. A synthetic mixture of long-chain AHLs mimicking that of the fadL mutant can improve rhizobial surface motility. When this AHL mixture was spotted into rhizosphere, the migration toward roots and rhizoplane colonization of S. fredii were enhanced in a diffusible way. This work adds novel parts managing extracellular AHLs, which modulate bacterial migration toward rhizoplane. The FadL-ExoFQP system is conserved in Alphaproteobacteria and may shape the "home life" of diverse keystone rhizobacteria.

Symbiotic Nodule Development and Efficiency in the Medicago truncatula Mtefd-1 Mutant Is Highly Dependent on Sinorhizobium Strains

Symbiotic nitrogen fixation (SNF) can play a key role in agroecosystems to reduce the negative impact of nitrogen fertilizers. Its efficiency is strongly affected by the combination of bacterial and plant genotypes, but the mechanisms responsible for the differences in the efficiency of rhizobium strains are not well documented. In Medicago truncatula, SNF has been mostly studied using model systems, such as M. truncatula A17 in interaction with Sinorhizobium meliloti Sm2011. Here we analyzed both the wild-type (wt) A17 and the Mtefd-1 mutant in interaction with five S. meliloti and two Sinorhizobium medicae strains. ETHYLENE RESPONSE FACTOR REQUIRED FOR NODULE DIFFERENTIATION (MtEFD) encodes a transcription factor, which contributes to the control of nodule number and differentiation in M. truncatula. We found that, in contrast to Sm2011, four strains induce functional (Fix+) nodules in Mtefd-1, although less efficient for SNF than in wt A17. In contrast, the Mtefd-1 hypernodulation phenotype is not strain-dependent. We compared the plant nodule transcriptomes in response to SmBL225C, a highly efficient strain with A17, versus Sm2011, in wt and Mtefd-1 backgrounds. This revealed faster nodule development with SmBL225C and early nodule senescence with Sm2011. These RNA sequencing analyses allowed us to identify candidate plant factors that could drive the differential nodule phenotype. In conclusion, this work shows the value of having a set of rhizobium strains to fully evaluate the biological importance of a plant symbiotic gene.

A dual legume-rhizobium transcriptome of symbiotic nodule senescence reveals coordinated plant and bacterial responses

Senescence determines plant organ lifespan depending on aging and environmental cues. During the endosymbiotic interaction with rhizobia, legume plants develop a specific organ, the root nodule, which houses nitrogen (N)-fixing bacteria. Unlike earlier processes of the legume-rhizobium interaction (nodule formation, N fixation), mechanisms controlling nodule senescence remain poorly understood. To identify nodule senescence-associated genes, we performed a dual plant-bacteria RNA sequencing approach on Medicago truncatula-Sinorhizobium meliloti nodules having initiated senescence either naturally (aging) or following an environmental trigger (nitrate treatment or salt stress). The resulting data allowed the identification of hundreds of plant and bacterial genes differentially regulated during nodule senescence, thus providing an unprecedented comprehensive resource of new candidate genes associated with this process. Remarkably, several plant and bacterial genes related to the cell cycle and stress responses were regulated in senescent nodules, including the rhizobial RpoE2-dependent general stress response. Analysis of selected core nodule senescence plant genes allowed showing that MtNAC969 and MtS40, both homologous to leaf senescence-associated genes, negatively regulate the transition between N fixation and senescence. In contrast, overexpression of a gene involved in the biosynthesis of cytokinins, well-known negative regulators of leaf senescence, may promote the transition from N fixation to senescence in nodules.

Genome Sequences of 17 Strains from Eight Races of Xanthomonas campestris pv. campestris

Xanthomonas campestris pv. campestris is a group of phytopathogenic bacteria causing black rot disease on Brassicaceae crops. Here, we report on draft genome sequences of 17 strains representing eight of nine known races of this pathogen, including the pathotype strain CFBP 6865.

Reference nodule transcriptomes for Melilotus officinalis and Medicago sativa cv. Algonquin

Host/symbiont compatibility is a hallmark of the symbiotic nitrogen-fixing interaction between rhizobia and legumes, mediated in part by plant-produced nodule-specific cysteine-rich (NCR) peptides and the bacterial BacA membrane protein that can act as a NCR peptide transporter. In addition, the genetic and metabolic properties supporting symbiotic nitrogen fixation often differ between compatible partners, including those sharing a common partner, highlighting the need for multiple study systems. Here, we report high-quality nodule transcriptome assemblies for Medicago sativa cv. Algonquin and Melilotus officinalis, two legumes able to form compatible symbioses with Sinorhizobium meliloti. The compressed M. sativa and M. officinalis assemblies consisted of 79,978 and 64,593 contigs, respectively, of which 33,341 and 28,278 were assigned putative annotations, respectively. As expected, the two transcriptomes showed broad similarity at a global level. We were particularly interested in the NCR peptide profiles of these plants, as these peptides drive bacterial differentiation during the symbiosis. A total of 412 and 308 NCR peptides were predicted from the M. sativa and M. officinalis transcriptomes, respectively, with approximately 9% of the transcriptome of both species consisting of NCR transcripts. Notably, transcripts encoding highly cationic NCR peptides (isoelectric point > 9.5), which are known to have antimicrobial properties, were ∼2-fold more abundant in M. sativa than in M. officinalis, and ∼27-fold more abundant when considering only NCR peptides in the six-cysteine class. We hypothesize that the difference in abundance of highly cationic NCR peptides explains our previous observation that some rhizobial bacA alleles which can support symbiosis with M. officinalis are unable to support symbiosis with M. sativa.

Xanthomonas transcriptome inside cauliflower hydathodes reveals bacterial virulence strategies and physiological adaptations at early infection stages

Xanthomonas campestris pv. campestris (Xcc) is a seed-transmitted vascular pathogen causing black rot disease on cultivated and wild Brassicaceae. Xcc enters the plant tissues preferentially via hydathodes, which are organs localized at leaf margins. To decipher both physiological and virulence strategies deployed by Xcc during early stages of infection, the transcriptomic profile of Xcc was analysed 3 days after entry into cauliflower hydathodes. Despite the absence of visible plant tissue alterations and despite a biotrophic lifestyle, 18% of Xcc genes were differentially expressed, including a striking repression of chemotaxis and motility functions. The Xcc full repertoire of virulence factors had not yet been activated but the expression of the HrpG regulon composed of 95 genes, including genes coding for the type III secretion machinery important for suppression of plant immunity, was induced. The expression of genes involved in metabolic adaptations such as catabolism of plant compounds, transport functions, sulphur and phosphate metabolism was upregulated while limited stress responses were observed 3 days postinfection. We confirmed experimentally that high-affinity phosphate transport is needed for bacterial fitness inside hydathodes. This analysis provides information about the nutritional and stress status of bacteria during the early biotrophic infection stages and helps to decipher the adaptive strategy of Xcc to the hydathode environment.

Prokaryotic Genome Annotation

In the last decade, the high-throughput and relatively low cost of short-read sequencing technologies have revolutionized prokaryotic genomics. This has led to an exponential increase in the number of bacterial and archaeal genome sequences available, as well as corresponding increase of genome assembly and annotation tools developed. Together, these hardware and software technologies have given scientists unprecedented options to study their chosen microbial systems without the need for large teams of bioinformaticists or supercomputing facilities. While these analysis tools largely fall into only a few categories, each may have different requirements, caveats and file formats, and some may be rarely updated or even abandoned. And so, despite the apparent ease in sequencing and analyzing a prokaryotic genome, it is no wonder that the budding genomicist may quickly find oneself overwhelmed. Here, we aim to provide the reader with an overview of genome annotation and its most important considerations, as well as an easy-to-follow protocol to get started with annotating a prokaryotic genome.

Easy identification of insertion sequence mobilization events in related bacterial strains with ISCompare

Bacterial genomes are composed of core and accessory genomes. The first is composed of housekeeping and essential genes, while the second is highly enriched in mobile genetic elements, including transposable elements (TEs). Insertion sequences (ISs), the smallest TEs, have an important role in genome evolution, and contribute to bacterial genome plasticity and adaptability. ISs can spread in a genome, presenting different locations in nearly related strains, and producing phenotypic variations. Few tools are available which can identify differentially located ISs (DLISs) on assembled genomes. Here, we introduce ISCompare, a new program to profile IS mobilization events in related bacterial strains using complete or draft genome assemblies. ISCompare was validated using artificial genomes with simulated random IS insertions and real sequences, achieving the same or better results than other available tools, with the advantage that ISCompare can analyze multiple ISs at the same time and outputs a list of candidate DLISs. ISCompare provides an easy and straightforward approach to look for differentially located ISs on bacterial genomes.

The two-component system ActJK is involved in acid stress tolerance and symbiosis in Sinorhizobium meliloti

The nitrogen-fixing α-proteobacterium Sinorhizobium meliloti genome codifies at least 50 response regulator (RR) proteins mediating different and, in many cases, unknown processes. RR-mutant library screening allowed us to identify genes potentially implicated in survival to acid conditions. actJ mutation resulted in a strain with reduced growth rate under mildly acidic conditions as well as a lower capacity to tolerate a sudden shift to lethal acidic conditions compared with the parental strain. Mutation of the downstream gene actK, which encodes for a histidine kinase, showed a similar phenotype in acidic environments suggesting a functional two-component system. Interestingly, even though nodulation kinetics, quantity, and macroscopic morphology of Medicago sativa nodules were not affected in actJ and actK mutants, ActK was required to express the wild-type nitrogen fixation phenotype and ActJK was necessary for full bacteroid development and nodule occupancy. The actJK regulatory system presented here provides insights into an evolutionary process in rhizobium adaptation to acidic environments and suggests that actJK-controlled functions are crucial for optimal symbiosis development.

Nonadditive Transcriptomic Signatures of Genotype-by-Genotype Interactions during the Initiation of Plant-Rhizobium Symbiosis

Rhizobia are ecologically important, facultative plant-symbiotic microbes. In nature, there is a large variability in the association of rhizobial strains and host plants of the same species. Here, we evaluated whether plant and rhizobial genotypes influence the initial transcriptional response of rhizobium following perception of a host plant. RNA sequencing of the model rhizobium Sinorhizobium meliloti exposed to root exudates or luteolin (an inducer of nod genes, involved in the early steps of symbiotic interaction) was performed on a combination of three S. meliloti strains and three alfalfa varieties as host plants. The response to root exudates involved hundreds of changes in the rhizobium transcriptome. Of the differentially expressed genes, 35% were influenced by the strain genotype, 16% were influenced by the plant genotype, and 29% were influenced by strain-by-host plant genotype interactions. We also examined the response of a hybrid S. meliloti strain in which the symbiotic megaplasmid (∼20% of the genome) was mobilized between two of the above-mentioned strains. Dozens of genes were upregulated in the hybrid strain, indicative of nonadditive variation in the transcriptome. In conclusion, this study demonstrated that transcriptional responses of rhizobia upon perception of legumes are influenced by the genotypes of both symbiotic partners and their interaction, suggesting a wide spectrum of genetic determinants involved in the phenotypic variation of plant-rhizobium symbiosis.IMPORTANCE A sustainable way for meeting the need of an increased global food demand should be based on a holobiont perspective, viewing crop plants as intimately associated with their microbiome, which helps improve plant nutrition, tolerance to pests, and adverse climate conditions. However, the genetic repertoire needed for efficient association with plants by the microbial symbionts is still poorly understood. The rhizobia are an exemplary model of facultative plant symbiotic microbes. Here, we evaluated whether genotype-by-genotype interactions could be identified in the initial transcriptional response of rhizobium perception of a host plant. We performed an RNA sequencing study to analyze the transcriptomes of different rhizobial strains elicited by root exudates of three alfalfa varieties as a proxy of an early step of the symbiotic interaction. The results indicated strain- and plant variety-dependent variability in the observed transcriptional changes, providing fundamentally novel insights into the genetic basis of rhizobium-plant interactions. Our results provide genetic insights and perspective to aid in the exploitation of natural rhizobium variation for improvement of legume growth in agricultural ecosystems.

Nonadditive Transcriptomic Signatures of Genotype-by-Genotype Interactions during the Initiation of Plant-Rhizobium Symbiosis

Rhizobia are ecologically important, facultative plant-symbiotic microbes. In nature, there is a large variability in the association of rhizobial strains and host plants of the same species. Here, we evaluated whether plant and rhizobial genotypes influence the initial transcriptional response of rhizobium following perception of a host plant. RNA sequencing of the model rhizobium Sinorhizobium meliloti exposed to root exudates or luteolin (an inducer of nod genes, involved in the early steps of symbiotic interaction) was performed on a combination of three S. meliloti strains and three alfalfa varieties as host plants. The response to root exudates involved hundreds of changes in the rhizobium transcriptome. Of the differentially expressed genes, 35% were influenced by the strain genotype, 16% were influenced by the plant genotype, and 29% were influenced by strain-by-host plant genotype interactions. We also examined the response of a hybrid S. meliloti strain in which the symbiotic megaplasmid (∼20% of the genome) was mobilized between two of the above-mentioned strains. Dozens of genes were upregulated in the hybrid strain, indicative of nonadditive variation in the transcriptome. In conclusion, this study demonstrated that transcriptional responses of rhizobia upon perception of legumes are influenced by the genotypes of both symbiotic partners and their interaction, suggesting a wide spectrum of genetic determinants involved in the phenotypic variation of plant-rhizobium symbiosis.IMPORTANCE A sustainable way for meeting the need of an increased global food demand should be based on a holobiont perspective, viewing crop plants as intimately associated with their microbiome, which helps improve plant nutrition, tolerance to pests, and adverse climate conditions. However, the genetic repertoire needed for efficient association with plants by the microbial symbionts is still poorly understood. The rhizobia are an exemplary model of facultative plant symbiotic microbes. Here, we evaluated whether genotype-by-genotype interactions could be identified in the initial transcriptional response of rhizobium perception of a host plant. We performed an RNA sequencing study to analyze the transcriptomes of different rhizobial strains elicited by root exudates of three alfalfa varieties as a proxy of an early step of the symbiotic interaction. The results indicated strain- and plant variety-dependent variability in the observed transcriptional changes, providing fundamentally novel insights into the genetic basis of rhizobium-plant interactions. Our results provide genetic insights and perspective to aid in the exploitation of natural rhizobium variation for improvement of legume growth in agricultural ecosystems.

Whole-Genome Sequencing of Bradyrhizobium diazoefficiens 113-2 and Comparative Genomic Analysis Provide Molecular Insights Into Species Specificity and Host Specificity

In the present study, we sequenced the complete genome of Bradyrhizobium diazoefficiens 113-2. The genomic characteristics of six selected rhizobial strains (two fast-growing rhizobia, two medium-slow-growing rhizobia and two slow-growing rhizobia) with four different legume hosts were analyzed by comparative genomic analysis. Genomes of B. diazoefficiens 113-2 and B. diazoefficiens USDA110 were found to share a large synteny blocks and a high ANI value, supporting 113-2 as a strain of B. diazoefficiens. 5,455 singletons and 11,656 clusters were identified among the six rhizobia genomes, and most of the pair-wise comparisons clusters were shared by the two genomes of strains in the same genus. Similar genus-specific gene numbers in the assigned COG functional terms were present in the two strains of the same genus, while the numbers were decreased with the increase of growth rate in most of the COG terms. KEGG pathway analysis of B. diazoefficiens 113-2 suggested that the rhizobial genes in ABC transporters and Two-Component system were mainly species-specific. Besides, the candidate genes related to secretion system and surface polysaccharides biosynthesis in the genomes of the six strains were explored and compared. 39 nodulation gene families, 12 nif gene families and 10 fix gene families in the genomes of these six strains were identified, and gene classes in most of gene families and the types and total gene numbers of gene families were substantially different among these six genomes. We also performed synteny analyses for above-mentioned nod, nif, and fix gene groupings, and selected NodW, NolK, NoeJ, NifB, FixK, and FixJ gene families to perform phylogeny analyses. Our results provided valuable molecular insights into species specificity and host specificity. The genetic information responsible for host specificity will play important roles in expanding the host range of rhizobia among legumes, which might provide new clues for the understanding of the genetic determinants of non-legume-rhizobium symbiosis.

Tyrosine Nitration of Flagellins: a Response of Sinorhizobium meliloti to Nitrosative Stress

Rhizobia are bacteria which can either live as free organisms in the soil or interact with plants of the legume family with, as a result, the formation of root organs called nodules in which differentiated endosymbiotic bacteria fix atmospheric nitrogen to the plant's benefit. In both lifestyles, rhizobia are exposed to nitric oxide (NO) which can be perceived as a signaling or toxic molecule. NO can act at the transcriptional level but can also modify proteins by S-nitrosylation of cysteine or nitration of tyrosine residues. However, only a few molecular targets of NO have been described in bacteria and none of them have been characterized in rhizobia. Here, we examined tyrosine nitration of Sinorhizobium meliloti proteins induced by NO. We found three tyrosine-nitrated proteins in S. meliloti grown under free-living conditions, in response to an NO donor. Two nitroproteins were identified by mass spectrometry and correspond to flagellins A and B. We showed that one of the nitratable tyrosines is essential to flagellin function in motility.IMPORTANCE Rhizobia are found as free-living bacteria in the soil or in interaction with plants and are exposed to nitric oxide (NO) in both environments. NO is known to have many effects on animals, plants, and bacteria where only a few molecular targets of NO have been described so far. We identified flagellin A and B by mass spectrometry as tyrosine-nitrated proteins in Sinorhizobium meliloti in vivo We also showed that one of the nitratable tyrosines is essential to flagellin function in motility. The results enhanced our understanding of NO effects on rhizobia. Identification of bacterial flagellin nitration opens a new possible role of NO in plant-microbe interactions.

Glutaredoxins Play an Important Role in the Redox Homeostasis and Symbiotic Capacity of Azorhizobium caulinodans ORS571

Glutaredoxin (GRX) plays an essential role in the control of the cellular redox state and related pathways in many organisms. There is limited information on GRXs from the model nitrogen (N₂)-fixing bacterium Azorhizobium caulinodans. In the present work, we identified and performed functional analyses of monothiol and dithiol GRXs in A. caulinodans in the free-living state and during symbiosis with Sesbania rostrata. Our data show that monothiol GRXs may be very important for bacterial growth under normal conditions and in response to oxidative stress due to imbalance of the redox state in grx mutants of A. caulinodans. Functional redundancies were also observed within monothiol and dithiol GRXs in terms of different physiological functions. The changes in catalase activity and iron content in grx mutants were assumed to favor the maintenance of bacterial resistance against oxidants, nodulation, and N₂ fixation efficiency in this bacterium. Furthermore, the monothiol GRX12 and dithiol GRX34 play a collective role in symbiotic associations between A. caulinodans and Sesbania rostrata. Our study provided systematic evidence that further investigations are required to understand the importance of glutaredoxins in A. caulinodans and other rhizobia.

Responses of mature symbiotic nodules to the whole-plant systemic nitrogen signaling

In symbiotic root nodules of legumes, terminally differentiated rhizobia fix atmospheric N2 producing an NH4+ influx that is assimilated by the plant. The plant, in return, provides photosynthates that fuel the symbiotic nitrogen acquisition. Mechanisms responsible for the adjustment of the symbiotic capacity to the plant N demand remain poorly understood. We have investigated the role of systemic signaling of whole-plant N demand on the mature N2-fixing nodules of the model symbiotic association Medicago truncatula/Sinorhizobium using split-root systems. The whole-plant N-satiety signaling rapidly triggers reductions of both N2 fixation and allocation of sugars to the nodule. These responses are associated with the induction of nodule senescence and the activation of plant defenses against microbes, as well as variations in sugars transport and nodule metabolism. The whole-plant N-deficit responses mirror these changes: a rapid increase of sucrose allocation in response to N-deficit is associated with a stimulation of nodule functioning and development resulting in nodule expansion in the long term. Physiological, transcriptomic, and metabolomic data together provide evidence for strong integration of symbiotic nodules into whole-plant nitrogen demand by systemic signaling and suggest roles for sugar allocation and hormones in the signaling mechanisms.

Riboregulation in Nitrogen-Fixing Endosymbiotic Bacteria

Small non-coding RNAs (sRNAs) are ubiquitous components of bacterial adaptive regulatory networks underlying stress responses and chronic intracellular infection of eukaryotic hosts. Thus, sRNA-mediated regulation of gene expression is expected to play a major role in the establishment of mutualistic root nodule endosymbiosis between nitrogen-fixing rhizobia and legume plants. However, knowledge about this level of genetic regulation in this group of plant-interacting bacteria is still rather scarce. Here, we review insights into the rhizobial non-coding transcriptome and sRNA-mediated post-transcriptional regulation of symbiotic relevant traits such as nutrient uptake, cell cycle, quorum sensing, or nodule development. We provide details about the transcriptional control and protein-assisted activity mechanisms of the functionally characterized sRNAs involved in these processes. Finally, we discuss the forthcoming research on riboregulation in legume symbionts.

A partner-switching system controls activation of mixed-linkage β-glucan synthesis by c-di-GMP in Sinorhizobium meliloti

Sinorhizobium meliloti synthesizes a linear mixed-linkage (1 → 3)(1 → 4)-β-d-glucan (ML β-glucan, MLG) in response to high levels of cyclic diguanylate (c-di-GMP). Two proteins BgsA and BgsB are required for MLG synthesis, BgsA being the glucan synthase which is activated upon c-di-GMP binding to its C-terminal domain. Here we report that the product of bgrR (SMb20447) is a diguanylate cyclase (DGC) that provides c-di-GMP for the synthesis of MLG by BgsA. bgrR is the first gene of a hexacistronic bgrRSTUWV operon, likely encoding a partner-switching regulatory network where BgrR is the final target. Using different approaches, we have determined that the products of genes bgrU (containing a putative PP2C serine phosphatase domain) and bgrW (with predicted kinase effector domain), modulate the phosphorylation status and the activity of the STAS domain protein BgrV. We propose that unphosphorylated BgrV inhibits BgrR DGC activity, perhaps through direct protein-protein interactions as established for other partner switchers. A bgrRSTUWV operon coexists with MLG structural bgsBA genes in many rhizobial genomes but is also present in some MLG non-producers, suggesting a role of this partner-switching system in other processes besides MLG biosynthesis.

Stress-inducible NHEJ in bacteria: function in DNA repair and acquisition of heterologous DNA

DNA double-strand breaks (DSB) in bacteria can be repaired by non-homologous end-joining (NHEJ), a two-component system relying on Ku and LigD. While performing a genetic characterization of NHEJ in Sinorhizobium meliloti, a representative of bacterial species encoding several Ku and LigD orthologues, we found that at least two distinct functional NHEJ repair pathways co-exist: one is dependent on Ku2 and LigD2, while the other depends on Ku3, Ku4 and LigD4. Whereas Ku2 likely acts as canonical bacterial Ku homodimers, genetic evidences suggest that Ku3-Ku4 form eukaryotic-like heterodimers. Strikingly, we found that the efficiency of both NHEJ systems increases under stress conditions, including heat and nutrient starvation. We found that this stimulation results from the transcriptional up-regulation of the ku and/or ligD genes, and that some of these genes are controlled by the general stress response regulator RpoE2. Finally, we provided evidence that NHEJ not only repairs DSBs, but can also capture heterologous DNA fragments into genomic breaks. Our data therefore suggest that NHEJ could participate to horizontal gene transfer from distantly related species, bypassing the need of homology to integrate exogenous DNA. This supports the hypothesis that NHEJ contributes to evolution and adaptation of bacteria under adverse environmental conditions.

Characterization of the Sinorhizobium meliloti HslUV and ClpXP Protease Systems in Free-Living and Symbiotic States

Symbiotic nitrogen fixation (SNF) in the interaction between the soil bacteria Sinorhizobium meliloti and legume plant Medicago sativa is carried out in specialized root organs called nodules. During nodule development, each symbiont must drastically alter their proteins, transcripts, and metabolites in order to support nitrogen fixation. Moreover, bacteria within the nodules are under stress, including challenges by plant antimicrobial peptides, low pH, limited oxygen availability, and strongly reducing conditions, all of which challenge proteome integrity. S. meliloti stress adaptation, proteome remodeling, and quality control are controlled in part by the large oligomeric protease complexes HslUV and ClpXP1. To improve understanding of the roles of S. meliloti HslUV and ClpXP1 under free-living conditions and in symbiosis with M. sativa, we generated ΔhslU, ΔhslV, ΔhslUV, and ΔclpP1 knockout mutants. The shoot dry weight of M. sativa plants inoculated with each deletion mutant was significantly reduced, suggesting a role in symbiosis. Further, slower free-living growth of the ΔhslUV and ΔclpP1 mutants suggests that HslUV and ClpP1 were involved in adapting to heat stress, the while ΔhslU and ΔclpP1 mutants were sensitive to kanamycin. All deletion mutants produced less exopolysaccharide and succinoglycan, as shown by replicate spot plating and calcofluor binding. We also generated endogenous C-terminal enhanced green fluorescent protein (eGFP) fusions to HslU, HslV, ClpX, and ClpP1 in S. meliloti Using anti-eGFP antibodies, native coimmunoprecipitation experiments with proteins from free-living and nodule tissues were performed and analyzed by mass spectrometry. The results suggest that HslUV and ClpXP were closely associated with ribosomal and proteome quality control proteins, and they identified several novel putative protein-protein interactions.IMPORTANCE Symbiotic nitrogen fixation (SNF) is the primary means by which biologically available nitrogen enters the biosphere, and it is therefore a critical component of the global nitrogen cycle and modern agriculture. SNF is the result of highly coordinated interactions between legume plants and soil bacteria collectively referred to as rhizobia, e.g., Medicago sativa and S. meliloti, respectively. Accomplishing SNF requires significant proteome changes in both organisms to create a microaerobic environment suitable for high-level bacterial nitrogenase activity. The bacterial protease systems HslUV and ClpXP are important in proteome quality control, in metabolic remodeling, and in adapting to stress. This work shows that S. meliloti HslUV and ClpXP are involved in SNF, in exopolysaccharide production, and in free-living stress adaptation.

Characterization of l-Carnitine Metabolism in Sinorhizobium meliloti

l-Carnitine is a trimethylammonium compound mostly known for its contribution to fatty acid transport into mitochondria. In bacteria, it is synthesized from γ-butyrobetaine (GBB) and can be used as a carbon source. l-Carnitine can be formed directly by GBB hydroxylation or synthesized via a biosynthetic route analogous to fatty acid degradation. However, this multistep pathway has not been experimentally characterized. In this work, we identified by gene context analysis a cluster of l-carnitine anabolic genes next to those involved in its catabolism and proceeded to the complete in vitro characterization of l-carnitine biosynthesis and degradation in Sinorhizobium meliloti The five enzymes catalyzing the seven steps that convert GBB to glycine betaine are described. Metabolomic analysis confirmed the multistage synthesis of l-carnitine in GBB-grown cells but also revealed that GBB is synthesized by S. meliloti To our knowledge, this is the first report of aerobic GBB synthesis in bacteria. The conservation of l-carnitine metabolism genes in different bacterial taxonomic classes underscores the role of l-carnitine as a ubiquitous nutrient.IMPORTANCE The experimental characterization of novel metabolic pathways is essential for realizing the value of genome sequences and improving our knowledge of the enzymatic capabilities of the bacterial world. However, 30% to 40% of genes of a typical genome remain unannotated or associated with a putative function. We used enzyme kinetics, liquid chromatography-mass spectroscopy (LC-MS)-based metabolomics, and mutant phenotyping for the characterization of the metabolism of l-carnitine in Sinorhizobium meliloti to provide an accurate annotation of the corresponding genes. The occurrence of conserved gene clusters for carnitine metabolism in soil, plant-associated, and marine bacteria underlines the environmental abundance of carnitine and suggests this molecule might make a significant contribution to ecosystem nitrogen and carbon cycling.

Characterization of HicAB toxin-antitoxin module of Sinorhizobium meliloti

Background

Toxin-antitoxin (TA) systems are little genetic units generally composed of two genes encoding antitoxin and toxin. These systems are known to be involved in many functions that can lead to growth arrest and cell death. Among the different types of TA systems, the type II gathers together systems where the antitoxin directly binds and inhibits the toxin. Among these type II TA systems, the HicAB module is widely distributed in free-living Bacteria and Archaea and the toxin HicA functions via RNA binding and cleavage. The genome of the symbiotic Sinorhizobium meliloti encodes numerous TA systems and only a few of them are functional. Among the predicted TA systems, there is one homologous to HicAB modules.

Results

In this study, we characterize the HicAB toxin-antitoxin module of S. meliloti. The production of the HicA of S. meliloti in Escherichia coli cells abolishes growth and decreases cell viability. We show that expression of the HicB of S. meliloti counteracts HicA toxicity. The results of double hybrid assays and co-purification experiments allow demonstrating the interaction of HicB with the toxin HicA. Purified HicA, but not HicAB complex, is able to degrade ribosomal RNA in vitro. The analysis of separated domains of HicB protein permits us to define the antitoxin activity and the operator-binding domain.

Conclusions

This study points out the first characterization of the HicAB system of the symbiotic S. meliloti whereas HicA is a toxin with ribonuclease activity and HicB has two domains: the COOH-terminal one that binds the operator and the NH2-terminal one that inhibits the toxin.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1382-6) contains supplementary material, which is available to authorized users.

EuGene: An Automated Integrative Gene Finder for Eukaryotes and Prokaryotes

EuGene is an integrative gene finder applicable to both prokaryotic and eukaryotic genomes. EuGene annotated its first genome in 1999. Starting from genomic DNA sequences representing a complete genome, EuGene is able to predict the major transcript units in the genome from a variety of sources of information: statistical information, similarities with known transcripts and proteins, but also any GFF3 structured information supporting the presence or absence of specific types of elements. EuGene has been used to find genes in the plants Arabidopsis thaliana, Medicago truncatula, and Theobroma cacao; tomato, sunflower, and Rosa genomes; and in the nematode Meloidogyne incognita genome, among many others. The large fraction of plant in this list probably influenced EuGene development, especially in its capacities to withstand a genome with a large number of repeated regions and transposable elements.Depending on the sources of information used for prediction, EuGene can be considered as purely ab initio, purely similarity based, or hybrid. With the general availability of NGS-transcribed sequence data in genome projects, EuGene adopts a default hybrid behavior that strongly relies on similarity information. Initially targeted at eukaryotic genomes, EuGene has also been extended to offer integrative gene prediction for bacteria, allowing for richer and robust predictions than either purely statistical or homology-based prokaryotic gene finders.This text has been written as a practical guide that will give you the capacity to train and execute EuGene on your favorite eukaryotic genome. As the prokaryotic case is simpler and has already been described, only the main differences with the eukaryotic version were reported.

Creation and Characterization of a Genomically Hybrid Strain in the Nitrogen-Fixing Symbiotic Bacterium Sinorhizobium meliloti

Many bacteria, often associated with eukaryotic hosts and of relevance for biotechnological applications, harbor a multipartite genome composed of more than one replicon. Biotechnologically relevant phenotypes are often encoded by genes residing on the secondary replicons. A synthetic biology approach to developing enhanced strains for biotechnological purposes could therefore involve merging pieces or entire replicons from multiple strains into a single genome. Here we report the creation of a genomic hybrid strain in a model multipartite genome species, the plant-symbiotic bacterium Sinorhizobium meliloti. We term this strain as cis-hybrid, since it is produced by genomic material coming from the same species' pangenome. In particular, we moved the secondary replicon pSymA (accounting for nearly 20% of total genome content) from a donor S. meliloti strain to an acceptor strain. The cis-hybrid strain was screened for a panel of complex phenotypes (carbon/nitrogen utilization phenotypes, intra- and extracellular metabolomes, symbiosis, and various microbiological tests). Additionally, metabolic network reconstruction and constraint-based modeling were employed for in silico prediction of metabolic flux reorganization. Phenotypes of the cis-hybrid strain were in good agreement with those of both parental strains. Interestingly, the symbiotic phenotype showed a marked cultivar-specific improvement with the cis-hybrid strains compared to both parental strains. These results provide a proof-of-principle for the feasibility of genome-wide replicon-based remodelling of bacterial strains for improved biotechnological applications in precision agriculture.

Most Sinorhizobium meliloti Extracytoplasmic Function Sigma Factors Control Accessory Functions

Fixed (reduced) soil nitrogen plays a critical role in soil fertility and successful food growth. Much soil fertility relies on symbiotic nitrogen fixation: the bacterial partner infects the host plant roots and reduces atmospheric dinitrogen in exchange for host metabolic fuel, a process that involves complex interactions between the partners mediated by changes in gene expression in each partner. Here we test the roles of a family of 11 extracytoplasmic function (ECF) gene regulatory proteins (sigma factors [σs]) that interact with RNA polymerase to determine if they play a significant role in establishing a nitrogen-fixing symbiosis or in responding to various stresses, including cell envelope stress. We discovered that symbiotic nitrogen fixation occurs even when all 11 of these regulatory genes are deleted, that most ECF sigma factors control accessory functions, and that none of the ECF sigma factors are required to survive envelope stress.

Bacteria must sense alterations in their environment and respond with changes in function and/or structure in order to cope. Extracytoplasmic function sigma factors (ECF σs) modulate transcription in response to cellular and environmental signals. The symbiotic nitrogen-fixing alphaproteobacterium Sinorhizobium meliloti carries genes for 11 ECF-like σs (RpoE1 to -E10 and FecI). We hypothesized that some of these play a role in mediating the interaction between the bacterium and its plant symbiotic partner. The bacterium senses changes in its immediate environment as it establishes contact with the plant root, initiates invasion of the plant as the root nodule is formed, traverses several root cell layers, and enters plant cortical cells via endocytosis. We used genetics, transcriptomics, and functionality to characterize the entire S. meliloti cohort of ECF σs. We discovered new targets for individual σs, confirmed others by overexpressing individual ECF σs, and identified or confirmed putative promoter motifs for nine of them. We constructed precise deletions of each ECF σ gene and its demonstrated or putative anti-σ gene and also a strain in which all 11 ECF σ and anti-σ genes were deleted. This all-ECF σ deletion strain showed no major defects in free-living growth, in Biolog Phenotype MicroArray assays, or in response to multiple stresses. None of the ECF σs were required for symbiosis on the host plants Medicago sativa and Medicago truncatula: the strain deleted for all ECF σ and anti-σ genes was symbiotically normal.

IMPORTANCE Fixed (reduced) soil nitrogen plays a critical role in soil fertility and successful food growth. Much soil fertility relies on symbiotic nitrogen fixation: the bacterial partner infects the host plant roots and reduces atmospheric dinitrogen in exchange for host metabolic fuel, a process that involves complex interactions between the partners mediated by changes in gene expression in each partner. Here we test the roles of a family of 11 extracytoplasmic function (ECF) gene regulatory proteins (sigma factors [σs]) that interact with RNA polymerase to determine if they play a significant role in establishing a nitrogen-fixing symbiosis or in responding to various stresses, including cell envelope stress. We discovered that symbiotic nitrogen fixation occurs even when all 11 of these regulatory genes are deleted, that most ECF sigma factors control accessory functions, and that none of the ECF sigma factors are required to survive envelope stress.

Functional Genomics Approaches to Studying Symbioses between Legumes and Nitrogen-Fixing Rhizobia

Biological nitrogen fixation gives legumes a pronounced growth advantage in nitrogen-deprived soils and is of considerable ecological and economic interest. In exchange for reduced atmospheric nitrogen, typically given to the plant in the form of amides or ureides, the legume provides nitrogen-fixing rhizobia with nutrients and highly specialised root structures called nodules. To elucidate the molecular basis underlying physiological adaptations on a genome-wide scale, functional genomics approaches, such as transcriptomics, proteomics, and metabolomics, have been used. This review presents an overview of the different functional genomics approaches that have been performed on rhizobial symbiosis, with a focus on studies investigating the molecular mechanisms used by the bacterial partner to interact with the legume. While rhizobia belonging to the alpha-proteobacterial group (alpha-rhizobia) have been well studied, few studies to date have investigated this process in beta-proteobacteria (beta-rhizobia).

Inter-replicon Gene Flow Contributes to Transcriptional Integration in the Sinorhizobium meliloti Multipartite Genome

Integration of newly acquired genes into existing regulatory networks is necessary for successful horizontal gene transfer (HGT). Ten percent of bacterial species contain at least two DNA replicons over 300 kilobases in size, with the secondary replicons derived predominately through HGT. The Sinorhizobium meliloti genome is split between a 3.7 Mb chromosome, a 1.7 Mb chromid consisting largely of genes acquired through ancient HGT, and a 1.4 Mb megaplasmid consisting primarily of recently acquired genes. Here, RNA-sequencing is used to examine the transcriptional consequences of massive, synthetic genome reduction produced through the removal of the megaplasmid and/or the chromid. Removal of the pSymA megaplasmid influenced the transcription of only six genes. In contrast, removal of the chromid influenced expression of ∼8% of chromosomal genes and ∼4% of megaplasmid genes. This was mediated in part by the loss of the ETR DNA region whose presence on pSymB is due to a translocation from the chromosome. No obvious functional bias among the up-regulated genes was detected, although genes with putative homologs on the chromid were enriched. Down-regulated genes were enriched in motility and sensory transduction pathways. Four transcripts were examined further, and in each case the transcriptional change could be traced to loss of specific pSymB regions. In particularly, a chromosomal transporter was induced due to deletion of bdhA likely mediated through 3-hydroxybutyrate accumulation. These data provide new insights into the evolution of the multipartite bacterial genome, and more generally into the integration of horizontally acquired genes into the transcriptome.

The complete replicons of 16 Ensifer meliloti strains offer insights into intra- and inter-replicon gene transfer, transposon-associated loci, and repeat elements

Ensifer meliloti (formerly Rhizobium meliloti and Sinorhizobium meliloti) is a model bacterium for understanding legume–rhizobial symbioses. The tripartite genome of E. meliloti consists of a chromosome, pSymA and pSymB, and in some instances strain-specific accessory plasmids. The majority of previous sequencing studies have relied on the use of assemblies generated from short read sequencing, which leads to gaps and assembly errors. Here we used PacBio-based, long-read assemblies and were able to assemble, de novo, complete circular replicons. In this study, we sequenced, de novo-assembled and analysed 10 E. meliloti strains. Sequence comparisons were also done with data from six previously published genomes. We identified genome differences between the replicons, including mol% G+C and gene content, nucleotide repeats, and transposon-associated loci. Additionally, genomic rearrangements both within and between replicons were identified, providing insight into evolutionary processes at the structural level. There were few cases of inter-replicon gene transfer of core genes between the main replicons. Accessory plasmids were more similar to pSymA than to either pSymB or the chromosome, with respect to gene content, transposon content and G+C content. In our population, the accessory plasmids appeared to share an open genome with pSymA, which contains many nodulation- and nitrogen fixation-related genes. This may explain previous observations that horizontal gene transfer has a greater effect on the content of pSymA than pSymB, or the chromosome, and why some rhizobia show unstable nodulation phenotypes on legume hosts.

Comparative analysis of rhizobial chromosomes and plasmids to estimate their evolutionary relationships

In the present study, complete genomic sequences retrieved from 57 rhizobial strains that covered four genera including 11 species were analyzed comprehensively. The four types of replicons: chromosomes, chromids, nonsymbiotic plasmids, and symbiotic plasmids were investigated and compared among these strains. Results showed that co-evolution occurred among these four replicons based on the similarities in average nucleotide identity. High correlation coefficient r values were observed between chromosomes and chromids, as well as between chromosomes and nonsymbiotic plasmids. Chromosomes and symbiotic plasmids showed different phylogenetic topology based on their core genes. Population structure analyses were performed to extrapolate the evolutionary histories of the test strains based on their chromosomal and symbiotic plasmid background. This resulted in seven ancestral types for chromosomal genes and three ancestral types for symbiotic plasmid genes. Rhizobial strains containing chromosome genes with ancestral type E tend to contain symbiotic plasmid genes with ancestral type II, while rhizobial strains containing chromosome genes with ancestral type G tend to contain symbiotic plasmid genes with ancestral type III. Seventeen strains associated with different host plant species which harbored the symbiotic genes with ancestral type I, exhibited high genetic diversity. In addition, Fu's test of the symbiotic plasmid genes with ancestral type III had undergone an expansion event, implying the influence of negative selection on these symbiotic plasmid genes.

Transcriptomic Studies of the Effect of nod Gene-Inducing Molecules in Rhizobia: Different Weapons, One Purpose

Simultaneous quantification of transcripts of the whole bacterial genome allows the analysis of the global transcriptional response under changing conditions. RNA-seq and microarrays are the most used techniques to measure these transcriptomic changes, and both complement each other in transcriptome profiling. In this review, we exhaustively compiled the symbiosis-related transcriptomic reports (microarrays and RNA sequencing) carried out hitherto in rhizobia. This review is specially focused on transcriptomic changes that takes place when five rhizobial species, Bradyrhizobium japonicum (=diazoefficiens) USDA 110, Rhizobium leguminosarum biovar viciae 3841, Rhizobium tropici CIAT 899, Sinorhizobium (=Ensifer) meliloti 1021 and S. fredii HH103, recognize inducing flavonoids, plant-exuded phenolic compounds that activate the biosynthesis and export of Nod factors (NF) in all analysed rhizobia. Interestingly, our global transcriptomic comparison also indicates that each rhizobial species possesses its own arsenal of molecular weapons accompanying the set of NF in order to establish a successful interaction with host legumes.

Is there sufficient Ensifer and Rhizobium species diversity in UK farmland soils to support red clover (Trifolium pratense), white clover (T. repens), lucerne (Medicago sativa) and black medic (M. lupulina)?

Rhizobia play important roles in agriculture owing to their ability to fix nitrogen through a symbiosis with legumes. The specificity of rhizobia-legume associations means that underused legume species may depend on seed inoculation with their rhizobial partners. For black medic (Medicago lupulina) and lucerne (Medicago sativa) little is known about the natural prevalence of their rhizobial partner Ensifer meliloti in UK soils, so that the need for inoculating them is unclear. We analysed the site-dependence of rhizobial seed inoculation effects on the subsequent ability of rhizobial communities to form symbioses with four legume species (Medicago lupulina, M. sativa, Trifolium repens and T. pratense). At ten organic farms across the UK, a species-diverse legume based mixture (LBM) which included these four species was grown. The LBM seed was inoculated with a mix of commercial inocula specific for clover and lucerne. At each site, soil from the LBM treatment was compared to the soil sampled prior to the sowing of the LBM (the control). From each site and each of the two treatments, a suspension of soils was applied to seedlings of the four legume species and grown in axenic conditions for six weeks. Root nodules were counted and their rhizobia isolated. PCR and sequencing of a fragment of the gyrB gene from rhizobial isolates allowed identification of strains. The number of nodules on each of the four legume species was significantly increased when inoculated with soil from the LBM treatment compared to the control. Both the proportion of plants forming nodules and the number of nodules formed varied significantly by site, with sites significantly affecting the Medicago species but not the Trifolium species. These differences in nodulation were broadly reflected in plant biomass where site and treatment interacted; at some sites there was a significant advantage from inoculation with the commercial inoculum but not at others. In particular, this study has demonstrated the commercial merit of inoculation of lucerne with compatible rhizobia.

RNA silencing in plant symbiotic bacteria: Insights from a protein-centric view

Extensive work in model enterobacteria has evidenced that the RNA chaperone Hfq and several endoribonucleases, such as RNase E or RNase III, serve pivotal roles in small RNA-mediated post-transcriptional silencing of gene expression. Characterization of these protein hubs commonly provide global functional and mechanistic insights into complex sRNA regulatory networks. The legume endosymbiont Sinorhizobium meliloti is a non-classical model bacterium with a very complex lifestyle in which riboregulation is expected to play important adaptive functions. Here, we discuss current knowledge about RNA silencing in S. meliloti from the perspective of the activity of Hfq and a recently discovered endoribonuclease (YbeY) exhibiting unprecedented catalytic versatility for the cleavage of single- and double-stranded RNA molecules.

A Key Regulator of the Glycolytic and Gluconeogenic Central Metabolic Pathways in Sinorhizobium meliloti

The order Rhizobiales contains numerous agriculturally, biotechnologically, and medically important bacteria, including the rhizobia, and the genera Agrobacterium, Brucella, and Methylobacterium, among others. These organisms tend to be metabolically versatile, but there has been relatively little investigation into the regulation of their central carbon metabolic pathways. Here, RNA-sequencing and promoter fusion data are presented to show that the PckR protein is a key regulator of central carbon metabolism in Sinorhizobium meliloti; during growth with gluconeogenic substrates, PckR represses expression of the complete Entner-Doudoroff glycolytic pathway and induces expression of the pckA and fbaB gluconeogenic genes. Electrophoretic mobility shift assays indicate that PckR binds an imperfect palindromic sequence that overlaps the promoter or transcriptional start site in the negatively regulated promoters, or is present in tandem upstream the promoter motifs in the positively regulated promoters. Genetic and in vitro electrophoretic mobility shift assay experiments suggest that elevated concentrations of a PckR effector ligand results in the dissociation of PckR from its target binding site, and evidence is presented that suggests phosphoenolpyruvate may function as the effector. Characterization of missense pckR alleles identified three conserved residues important for increasing the affinity of PckR for its cognate effector molecule. Bioinformatics analyses illustrates that PckR is limited to a narrow phylogenetic range consisting of the Rhizobiaceae, Phyllobacteriaceae, Brucellaceae, and Bartonellaceae families. These data provide novel insights into the regulation of the core carbon metabolic pathways of this pertinent group of α-proteobacteria.

AraC-like transcriptional activator CuxR binds c-di-GMP by a PilZ-like mechanism to regulate extracellular polysaccharide production

Cyclic dimeric GMP (c-di-GMP) has emerged as a key regulatory player in the transition between planktonic and sedentary biofilm-associated bacterial lifestyles. It controls a multitude of processes including production of extracellular polysaccharides (EPSs). The PilZ domain, consisting of an N-terminal "RxxxR" motif and a β-barrel domain, represents a prototype c-di-GMP receptor. We identified a class of c-di-GMP-responsive proteins, represented by the AraC-like transcription factor CuxR in plant symbiotic α-proteobacteria. In Sinorhizobium meliloti, CuxR stimulates transcription of an EPS biosynthesis gene cluster at elevated c-di-GMP levels. CuxR consists of a Cupin domain, a helical hairpin, and bipartite helix-turn-helix motif. Although unrelated in sequence, the mode of c-di-GMP binding to CuxR is highly reminiscent to that of PilZ domains. c-di-GMP interacts with a conserved N-terminal RxxxR motif and the Cupin domain, thereby promoting CuxR dimerization and DNA binding. We unravel structure and mechanism of a previously unrecognized c-di-GMP-responsive transcription factor and provide insights into the molecular evolution of c-di-GMP binding to proteins.

A conserved α-proteobacterial small RNA contributes to osmoadaptation and symbiotic efficiency of rhizobia on legume roots

Small non-coding RNAs (sRNAs) are expected to have pivotal roles in the adaptive responses underlying symbiosis of nitrogen-fixing rhizobia with legumes. Here, we provide primary insights into the function and activity mechanism of the Sinorhizobium meliloti trans-sRNA NfeR1 (Nodule Formation Efficiency RNA). Northern blot probing and transcription tracking with fluorescent promoter-reporter fusions unveiled high nfeR1 expression in response to salt stress and throughout the symbiotic interaction. The strength and differential regulation of nfeR1 transcription are conferred by a motif, which is conserved in nfeR1 promoter regions in α-proteobacteria. NfeR1 loss-of-function compromised osmoadaptation of free-living bacteria, whilst causing misregulation of salt-responsive genes related to stress adaptation, osmolytes catabolism and membrane trafficking. Nodulation tests revealed that lack of NfeR1 affected competitiveness, infectivity, nodule development and symbiotic efficiency of S. meliloti on alfalfa roots. Comparative computer predictions and a genetic reporter assay evidenced a redundant role of three identical NfeR1 unpaired anti Shine-Dalgarno motifs for targeting and downregulation of translation of multiple mRNAs from transporter genes. Our data provide genetic evidence of the hyperosmotic conditions of the endosymbiotic compartments. NfeR1-mediated gene regulation in response to this cue could contribute to coordinate nutrient uptake with the metabolic reprogramming concomitant to symbiotic transitions.

RNase E and RNase J are needed for S-adenosylmethionine homeostasis in Sinorhizobium meliloti

The ribonucleases (RNases) E and J play major roles in E. coli and Bacillus subtilis, respectively, and co-exist in Sinorhizobium meliloti. We analysed S. meliloti 2011 mutants with mini-Tn5 insertions in the corresponding genes rne and rnj and found many overlapping effects. We observed similar changes in mRNA levels, including lower mRNA levels of the motility and chemotaxis related genes flaA, flgB and cheR and higher levels of ndvA (important for glucan export). The acyl-homoserine lactone (AHL) levels were also higher during exponential growth in both RNase mutants, despite no increase in the expression of the sinI AHL synthase gene. Furthermore, several RNAs from both mutants migrated aberrantly in denaturing gels at 300 V but not under stronger denaturing conditions at 1300 V. The similarities between the two mutants could be explained by increased levels of the key methyl donor S-adenosylmethionine (SAM), since this may result in faster AHL synthesis leading to higher AHL accumulation as well as in uncontrolled methylation of macromolecules including RNA, which may strengthen RNA secondary structures. Indeed, we found that in both mutants the N6-methyladenosine content was increased almost threefold and the SAM level was increased at least sevenfold. Complementation by induced ectopic expression of the respective RNase restored the AHL and SAM levels in each of the mutants. In summary, our data show that both RNase E and RNase J are needed for SAM homeostasis in S. meliloti.

DNA double-strand break repair is involved in desiccation resistance of Sinorhizobium meliloti, but is not essential for its symbiotic interaction with Medicago truncatula

The soil bacterium Sinorhizobium meliloti, a nitrogen-fixing symbiont of legume plants, is exposed to numerous stress conditions in nature, some of which cause the formation of harmful DNA double-strand breaks (DSBs). In particular, the reactive oxygen species (ROS) and the reactive nitrogen species (RNS) produced during symbiosis, and the desiccation occurring in dry soils, are conditions which induce DSBs. Two major systems of DSB repair are known in S. meliloti: homologous recombination (HR) and non-homologous end-joining (NHEJ). However, their role in the resistance to ROS, RNS and desiccation has never been examined in this bacterial species, and the importance of DSB repair in the symbiotic interaction has not been properly evaluated. Here, we constructed S. meliloti strains deficient in HR (by deleting the recA gene) or in NHEJ (by deleting the four ku genes) or both. Interestingly, we observed that ku and/or recA genes are involved in S. meliloti resistance to ROS and RNS. Nevertheless, an S. meliloti strain deficient in both HR and NHEJ was not altered in its ability to establish and maintain an efficient nitrogen-fixing symbiosis with Medicago truncatula, showing that rhizobial DSB repair is not essential for this process. This result suggests either that DSB formation in S. meliloti is efficiently prevented during symbiosis or that DSBs are not detrimental for symbiosis efficiency. In contrast, we found for the first time that both recA and ku genes are involved in S. meliloti resistance to desiccation, suggesting that DSB repair could be important for rhizobium persistence in the soil.

High-Quality Draft Genome Sequences of Two Xanthomonas Pathotype Strains Infecting Aroid Plants

We present here the draft genome sequences of bacterial pathogens of the Araceae family, Xanthomonas axonopodis pv. dieffenbachiae LMG 695 and Xanthomonas campestris pv. syngonii LMG 9055, differing in host range. A comparison between genome sequences will help understand the mechanisms involved in tissue specificity and adaptation to host plants.

Organic acid mediated repression of sugar utilization in rhizobia

Rhizobia are a class of symbiotic diazotrophic bacteria which utilize C4 acids in preference to sugars and the sugar utilization is repressed as long as C4 acids are present. This can be manifested as a diauxie when rhizobia are grown in the presence of a sugar and a C4 acid together. Succinate, a C4 acid is known to repress utilization of sugars, sugar alcohols, hydrocarbons, etc by a mechanism termed as Succinate Mediated Catabolite Repression (SMCR). Mechanism of catabolite repression determines the hierarchy of carbon source utilization in bacteria. Though the mechanism of catabolite repression has been well studied in model organisms like E. coli, B. subtilis and Pseudomonas sp., mechanism of SMCR in rhizobia has not been well elucidated. C4 acid uptake is important for effective symbioses while mutation in the sugar transport and utilization genes does not affect symbioses. Deletion of hpr and sma0113 resulted in the partial relief of SMCR of utilization of galactosides like lactose, raffinose and maltose in the presence of succinate. However, no such regulators governing SMCR of glucoside utilization have been identified till date. Though rhizobia can utilize multitude of sugars, high affinity transporters for many sugars are yet to be identified. Identifying high affinity sugar transporters and studying the mechanism of catabolite repression in rhizobia is important to understand the level of regulation of SMCR and the key regulators involved in SMCR.

Draft genome sequences of five Pseudomonas syringae pv. actinidifoliorum strains isolated in France

Pseudomonas syringae pv. actinidifoliorum causes necrotic spots on the leaves of Actinidia deliciosa and Actinidia chinensis. P. syringae pv. actinidifoliorum has been detected in New Zealand, Australia, France and Spain. Four lineages were previously identified within the P. syringae pv. actinidifoliorum species group. Here, we report the draft genome sequences of five strains of P. syringae pv. actinidifoliorum representative of lineages 1, 2 and 4, isolated in France. The whole genomes of strains isolated in New Zealand, representative of P. syringae pv. actinidifoliorum lineages 1 and 3, were previously sequenced. The availability of supplementary P. syringae pv. actinidifoliorum genome sequences will be useful for developing molecular tools for pathogen detection and for performing comparative genomic analyses to study the relationship between P. syringae pv. actinidifoliorum and other kiwifruit pathogens, such as P. syringae pv. actinidiae.

The rise of genomics

A brief history of the development of genomics is provided. Complete sequencing of genomes of uni- and multicellular organisms is based on important progress in sequencing and bioinformatics. Evolution of these methods is ongoing and has triggered an explosion in data production and analysis. Initial analyses focused on the inventory of genes encoding proteins. Completeness and quality of gene prediction remains crucial. Genome analyses profoundly modified our views on evolution, biodiversity and contributed to the detection of new functions, yet to be fully elucidated, such as those fulfilled by non-coding RNAs. Genomics has become the basis for the study of biology and provides the molecular support for a bunch of large-scale studies, the omics.

Genomic resources for identification of the minimal N2 -fixing symbiotic genome

The lack of an appropriate genomic platform has precluded the use of gain-of-function approaches to study the rhizobium-legume symbiosis, preventing the establishment of the genes necessary and sufficient for symbiotic nitrogen fixation (SNF) and potentially hindering synthetic biology approaches aimed at engineering this process. Here, we describe the development of an appropriate system by reverse engineering Sinorhizobium meliloti. Using a novel in vivo cloning procedure, the engA-tRNA-rmlC (ETR) region, essential for cell viability and symbiosis, was transferred from Sinorhizobium fredii to the ancestral location on the S. meliloti chromosome, rendering the ETR region on pSymB redundant. A derivative of this strain lacking both the large symbiotic replicons (pSymA and pSymB) was constructed. Transfer of pSymA and pSymB back into this strain restored symbiotic capabilities with alfalfa. To delineate the location of the single-copy genes essential for SNF on these replicons, we screened a S. meliloti deletion library, representing > 95% of the 2900 genes of the symbiotic replicons, for their phenotypes with alfalfa. Only four loci, accounting for < 12% of pSymA and pSymB, were essential for SNF. These regions will serve as our preliminary target of the minimal set of horizontally acquired genes necessary and sufficient for SNF.

Genomics and transcriptomics of Xanthomonas campestris species challenge the concept of core type III effectome

BACKGROUND

The bacterial species Xanthomonas campestris infects a wide range of Brassicaceae. Specific pathovars of this species cause black rot (pv. campestris), bacterial blight of stock (pv. incanae) or bacterial leaf spot (pv. raphani).

RESULTS

In this study, we extended the genomic coverage of the species by sequencing and annotating the genomes of strains from pathovar incanae (CFBP 1606R and CFBP 2527R), pathovar raphani (CFBP 5828R) and a pathovar formerly named barbareae (CFBP 5825R). While comparative analyses identified a large core ORFeome at the species level, the core type III effectome was limited to only three putative type III effectors (XopP, XopF1 and XopAL1). In Xanthomonas, these effector proteins are injected inside the plant cells by the type III secretion system and contribute collectively to virulence. A deep and strand-specific RNA sequencing strategy was adopted in order to experimentally refine genome annotation for strain CFBP 5828R. This approach also allowed the experimental definition of novel ORFs and non-coding RNA transcripts. Using a constitutively active allele of hrpG, a master regulator of the type III secretion system, a HrpG-dependent regulon of 141 genes co-regulated with the type III secretion system was identified. Importantly, all these genes but seven are positively regulated by HrpG and 56 of those encode components of the Hrp type III secretion system and putative effector proteins.

CONCLUSIONS

This dataset is an important resource to mine for novel type III effector proteins as well as for bacterial genes which could contribute to pathogenicity of X. campestris.

OperomeDB: A Database of Condition-Specific Transcription Units in Prokaryotic Genomes

Background. In prokaryotic organisms, a substantial fraction of adjacent genes are organized into operons—codirectionally organized genes in prokaryotic genomes with the presence of a common promoter and terminator. Although several available operon databases provide information with varying levels of reliability, very few resources provide experimentally supported results. Therefore, we believe that the biological community could benefit from having a new operon prediction database with operons predicted using next-generation RNA-seq datasets. Description. We present operomeDB, a database which provides an ensemble of all the predicted operons for bacterial genomes using available RNA-sequencing datasets across a wide range of experimental conditions. Although several studies have recently confirmed that prokaryotic operon structure is dynamic with significant alterations across environmental and experimental conditions, there are no comprehensive databases for studying such variations across prokaryotic transcriptomes. Currently our database contains nine bacterial organisms and 168 transcriptomes for which we predicted operons. User interface is simple and easy to use, in terms of visualization, downloading, and querying of data. In addition, because of its ability to load custom datasets, users can also compare their datasets with publicly available transcriptomic data of an organism. Conclusion. OperomeDB as a database should not only aid experimental groups working on transcriptome analysis of specific organisms but also enable studies related to computational and comparative operomics.

Genome Sequences of the Race 1 and Race 4 Xanthomonas campestris pv. campestris Strains CFBP 1869 and CFBP 5817

Xanthomonas campestris pv. campestris is the causal agent of black rot on Brassicaceae. The draft genome sequences of strains CFBP 1869 and CFBP 5817 have been determined and are the first ones corresponding to race 1 and race 4 strains, which have a predominant agronomic and economic impact on cabbage cultures worldwide.

Evolution of Intra-specific Regulatory Networks in a Multipartite Bacterial Genome

Reconstruction of the regulatory network is an important step in understanding how organisms control the expression of gene products and therefore phenotypes. Recent studies have pointed out the importance of regulatory network plasticity in bacterial adaptation and evolution. The evolution of such networks within and outside the species boundary is however still obscure. Sinorhizobium meliloti is an ideal species for such study, having three large replicons, many genomes available and a significant knowledge of its transcription factors (TF). Each replicon has a specific functional and evolutionary mark; which might also emerge from the analysis of their regulatory signatures. Here we have studied the plasticity of the regulatory network within and outside the S. meliloti species, looking for the presence of 41 TFs binding motifs in 51 strains and 5 related rhizobial species. We have detected a preference of several TFs for one of the three replicons, and the function of regulated genes was found to be in accordance with the overall replicon functional signature: house-keeping functions for the chromosome, metabolism for the chromid, symbiosis for the megaplasmid. This therefore suggests a replicon-specific wiring of the regulatory network in the S. meliloti species. At the same time a significant part of the predicted regulatory network is shared between the chromosome and the chromid, thus adding an additional layer by which the chromid integrates itself in the core genome. Furthermore, the regulatory network distance was found to be correlated with both promoter regions and accessory genome evolution inside the species, indicating that both pangenome compartments are involved in the regulatory network evolution. We also observed that genes which are not included in the species regulatory network are more likely to belong to the accessory genome, indicating that regulatory interactions should also be considered to predict gene conservation in bacterial pangenomes.

The influence of transcriptional regulatory networks on the evolution of bacterial pangenomes has not yet been elucidated, even though the role of transcriptional regulation is widely recognized. Using the model symbiont Sinorhizobium meliloti we have predicted the regulatory targets of 41 transcription factors in 51 strains and 5 other rhizobial species, showing a correlation between regulon diversity and pangenome evolution, through upstream sequence diversity and accessory genome composition. We have also shown that genes not wired to the regulatory network are more likely to belong to the accessory genome, thus suggesting that inclusion in the regulatory circuits may be an indicator of gene conservation. We have also highlighted a series of transcription factors that preferentially regulate genes belonging to one of the three replicons of this species, indicating the presence of replicon-specific regulatory modules, with peculiar functional signatures. At the same time the chromid shares a significant part of the regulatory network with the chromosome, indicating an additional way by which this replicon integrates itself in the pangenome.

Inactivation of group II intron RmInt1 in the Sinorhizobium meliloti genome

Group II introns are self-splicing catalytic RNAs that probably originated in bacteria and act as mobile retroelements. The dispersal and dynamics of group II intron spread within a bacterial genome are thought to follow a selection-driven extinction model. Likewise, various studies on the evolution of group II introns have suggested that they are evolving toward an inactive form by fragmentation, with the loss of the intron 3′-terminus, but with some intron fragments remaining and continuing to evolve in the genome. RmInt1 is a mobile group II intron that is widespread in natural populations of Sinorhizobium meliloti, but some strains of this species have no RmInt1 introns. We studied the splicing ability and mobility of the three full-length RmInt1 copies harbored by S. meliloti 1021, and obtained evidence suggesting that specific mutations may lead to the impairment of intron splicing and retrohoming. Our data suggest that the RmInt1 copies in this strain are undergoing a process of inactivation.