Coordinated regulation of core and accessory genes in the multipartite genome of Sinorhizobium fredii

Cofitness network connectivity determines a fuzzy essential zone in open bacterial pangenome

Most in silico evolutionary studies commonly assumed that core genes are essential for cellular function, while accessory genes are dispensable, particularly in nutrient-rich environments. However, this assumption is seldom tested genetically within the pangenome context. In this study, we conducted a robust pangenomic Tn-seq analysis of fitness genes in a nutrient-rich medium for Sinorhizobium strains with a canonical open pangenome. To evaluate the robustness of fitness category assignment, Tn-seq data for three independent mutant libraries per strain were analyzed by three methods, which indicates that the Hidden Markov Model (HMM)-based method is most robust to variations between mutant libraries and not sensitive to data size, outperforming the Bayesian and Monte Carlo simulation-based methods. Consequently, the HMM method was used to classify the fitness category. Fitness genes, categorized as essential (ES), advantage (GA), and disadvantage (GD) genes for growth, are enriched in core genes, while nonessential genes (NE) are over-represented in accessory genes. Accessory ES/GA genes showed a lower fitness effect than core ES/GA genes. Connectivity degrees in the cofitness network decrease in the order of ES, GD, and GA/NE. In addition to accessory genes, 1599 out of 3284 core genes display differential essentiality across test strains. Within the pangenome core, both shared quasi-essential (ES and GA) and strain-dependent fitness genes are enriched in similar functional categories. Our analysis demonstrates a considerable fuzzy essential zone determined by cofitness connectivity degrees in Sinorhizobium pangenome and highlights the power of the cofitness network in understanding the genetic basis of ever-increasing prokaryotic pangenome data.

Fulvic Acid Promotes Legume-Rhizobium Symbiosis by Stimulating Endogenous Flavonoids Synthesis and Secretion

Fulvic acid (FA) promotes symbiosis between legumes and rhizobia. To elucidate from the aspect of symbiosis, the effects of root irrigation of water-soluble humic materials (WSHM) or foliar spraying of its highly active component, FA, on soybean root exudates and on rhizosphere microorganisms were investigated. As a result, WSHM/FA treatments significantly altered root exudate metabolite composition, and isoflavonoids were identified as key contributors in both treatments compared to the control. Increased expression of genes related to the isoflavonoid biosynthesis were validated by RT-qPCR in both treatments, which notably elevated the synthesis of symbiotic signals genistein, daidzin, coumestrol, and biochanin A. Moreover, the WSHM/FA treatments induced a change in rhizosphere microbial community, coupled with an increase in the relative abundance of rhizobia. Our findings showed that WSHM/FA promotes symbiosis by stimulating the endogenous flavonoid synthesis and leads to rhizobia accumulation in the rhizosphere. This study provides new insights into mechanisms underlying the FA-mediated promotion of symbiosis.

Selective regulation of endophytic bacteria and gene expression in soybean by water-soluble humic materials

BACKGROUND

As part of the plant microbiome, endophytic bacteria play an essential role in plant growth and resistance to stress. Water-soluble humic materials (WSHM) is widely used in sustainable agriculture as a natural and non-polluting plant growth regulator to promote the growth of plants and beneficial bacteria. However, the mechanisms of WSHM to promote plant growth and the evidence for commensal endophytic bacteria interaction with their host remain largely unknown. Here, 16S rRNA gene sequencing, transcriptomic analysis, and culture-based methods were used to reveal the underlying mechanisms.

RESULTS

WSHM reduced the alpha diversity of soybean endophytic bacteria, but increased the bacterial interactions and further selectively enriched the potentially beneficial bacteria. Meanwhile, WSHM regulated the expression of various genes related to the MAPK signaling pathway, plant-pathogen interaction, hormone signal transduction, and synthetic pathways in soybean root. Omics integration analysis showed that Sphingobium was the genus closest to the significantly changed genes in WSHM treatment. The inoculation of endophytic Sphingobium sp. TBBS4 isolated from soybean significantly improved soybean nodulation and growth by increasing della gene expression and reducing ethylene release.

CONCLUSION

All the results revealed that WSHM promotes soybean nodulation and growth by selectively regulating soybean gene expression and regulating the endophytic bacterial community, Sphingobium was the key bacterium involved in plant-microbe interaction. These findings refined our understanding of the mechanism of WSHM promoting soybean nodulation and growth and provided novel evidence for plant-endophyte interaction.

Using Machine Learning to Predict Genes Underlying Differentiation of Multipartite and Unipartite Traits in Bacteria

Since the discovery of the second chromosome in the Rhodobacter sphaeroides 2.4.1 by Suwanto and Kaplan in 1989 and the revelation of gene sequences, multipartite genomes have been reported in over three hundred bacterial species under nine different phyla. This phenomenon shattered the dogma of a unipartite genome (a single circular chromosome) in bacteria. Recently, Artificial Intelligence (AI), machine learning (ML), and Deep Learning (DL) have emerged as powerful tools in the investigation of big data in a plethora of disciplines to decipher complex patterns in these data, including the large-scale analysis and interpretation of genomic data. An important inquiry in bacteriology pertains to the genetic factors that underlie the structural evolution of multipartite and unipartite bacterial species. Towards this goal, here we have attempted to leverage machine learning as a means to identify the genetic factors that underlie the differentiation of, in general, bacteria with multipartite genomes and bacteria with unipartite genomes. In this study, deploying ML algorithms yielded two gene lists of interest: one that contains 46 discriminatory genes obtained following an assessment on all gene sets, and another that contains 35 discriminatory genes obtained based on an investigation of genes that are differentially present (or absent) in the genomes of the multipartite bacteria and their respective close relatives. Our study revealed a small pool of genes that discriminate bacteria with multipartite genomes and their close relatives with single-chromosome genomes. Machine learning thus aided in uncovering the genetic factors that underlie the differentiation of bacterial multipartite and unipartite traits.

MicroRNA 4407 modulates nodulation in soybean by repressing a root-specific ISOPENTENYLTRANSFERASE (GmIPT3)

MicroRNAs (miRNAs) are important regulators of plant biological processes, including soybean nodulation. One miRNA, miR4407, was identified in soybean roots and nodules. However, the function of miR4407 in soybean is still unknown. MiR4407, unique to soybean, positively regulates lateral root emergence and root structures and represses a root-specific ISOPENTENYLTRANSFERASE (GmIPT3). By altering the expression of miR4407 and GmIPT3, we investigated the role of miR4407 in lateral root and nodule development. Both miR4407 and GmIPT3 are expressed in the inner root cortex and nodule primordia. Upon rhizobial inoculation, miR4407 was downregulated while GmIPT3 was upregulated. Overexpressing miR4407 reduced the number of nodules in transgenic soybean hairy roots while overexpressing the wild-type GmIPT3 or a miR4407-resistant GmIPT3 mutant (mGmIPT3) significantly increased the nodule number. The mechanism of miR4407 and GmIPT3 functions was also linked to autoregulation of nodulation (AON), where miR4407 overexpression repressed miR172c and activated its target, GmNNC1, turning on AON. Exogenous CK mimicked the effects of GmIPT3 overexpression on miR172c, supporting the notion that GmIPT3 regulates nodulation by enhancing root-derived CK. Overall, our data revealed a new miRNA-mediated regulatory mechanism of nodulation in soybean. MiR4407 showed a dual role in lateral root and nodule development.

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.

Intracellular common gardens reveal niche differentiation in transposable element community during bacterial adaptive evolution

The distribution and abundance of transposable elements across the tree of life have significantly shaped the evolution of cellular organisms, but the underlying mechanisms shaping these ecological patterns remain elusive. Here we establish a "common garden" approach to study causal ecological interactions between a xenogeneic conditional lethal sacB gene and the community of transposable insertion sequences (ISs) in a multipartite prokaryote genome. Xenogeneic sacB of low, medium, or high GC content was individually inserted into three replicons of a model bacterium Sinorhizobium fredii, and exhibited replicon- and GC-dependent variation in genetic stability. This variation was largely attributable to multidimensional niche differentiation for IS community members. The transposition efficiency of major active ISs depended on the nucleoid-associated xenogeneic silencer MucR. Experimentally eliminating insertion activity of specific ISs by deleting MucR strongly demonstrated a dominant role of niche differentiation among ISs. This intracellular common garden approach in the experimental evolution context allows not only for evaluating genetic stability of natural and synthetic xenogeneic genes of different sequence signatures in host cells but also for tracking and testing causal relationships in unifying ecological principles in genome ecology.

Adaptive Evolution of Rhizobial Symbiosis beyond Horizontal Gene Transfer: From Genome Innovation to Regulation Reconstruction

There are ubiquitous variations in symbiotic performance of different rhizobial strains associated with the same legume host in agricultural practices. This is due to polymorphisms of symbiosis genes and/or largely unexplored variations in integration efficiency of symbiotic function. Here, we reviewed cumulative evidence on integration mechanisms of symbiosis genes. Experimental evolution, in concert with reverse genetic studies based on pangenomics, suggests that gain of the same circuit of key symbiosis genes through horizontal gene transfer is necessary but sometimes insufficient for bacteria to establish an effective symbiosis with legumes. An intact genomic background of the recipient may not support the proper expression or functioning of newly acquired key symbiosis genes. Further adaptive evolution, through genome innovation and reconstruction of regulation networks, may confer the recipient of nascent nodulation and nitrogen fixation ability. Other accessory genes, either co-transferred with key symbiosis genes or stochastically transferred, may provide the recipient with additional adaptability in ever-fluctuating host and soil niches. Successful integrations of these accessory genes with the rewired core network, regarding both symbiotic and edaphic fitness, can optimize symbiotic efficiency in various natural and agricultural ecosystems. This progress also sheds light on the development of elite rhizobial inoculants using synthetic biology procedures.

Analysis of multipartite bacterial genomes using alignment free and alignment-based pipelines

Since the discovery of second chromosome in Rhodobacter sphaeroides 2.4.1 in 1989, multipartite genomes have been reported in over three hundred bacterial species under nine different phyla. This has shattered the unipartite (single chromosome) genome dogma in bacteria. Since then, many questions on various aspects of multipartite genomes in bacteria have been addressed. However, our understanding of how multipartite genomes emerge and evolve is still lacking. Importantly, the knowledge of genetic factors underlying the differences in multipartite and single-chromosome genomes is lacking. In this work, we have performed comparative evolutionary and functional genomics analyses to identify molecular factors that discriminate multipartite from unipartite bacteria, with the goal to decipher taxon-specific factors, and those that are prevalent across the taxa, underlying these traits. We assessed the roles of evolutionary mechanisms, specifically gene gain, in driving the divergence of bacteria with single and multiple chromosomes. In addition, we performed functional genomic analysis to garner support for our findings from comparative evolutionary analysis. We found genes such as those encoding conserved hypothetical proteins in Deinococcus radiodurans R1, and putative phage phi-C31 gp36 major capsid like and hypothetical proteins in Rhodobacter sphaeroides 2.4.1, which are located on accessory chromosomes in these bacteria but were not found in the inferred ancestral sequences, and on the primary chromosomes, as well as were not found in their closest relatives with single chromosome within the same clade. Our study shines a new light on the potential roles of the secondary chromosomes in helping bacteria with multipartite genomes to adapt to specialized environments or growth conditions.

The Ros/MucR Zinc-Finger Protein Family in Bacteria: Structure and Functions

Ros/MucR is a widespread family of bacterial zinc-finger-containing proteins that integrate multiple functions, such as symbiosis, virulence, transcription regulation, motility, production of surface components, and various other physiological processes in cells. This regulatory protein family is conserved in bacteria and is characterized by its zinc-finger motif, which has been proposed as the ancestral domain from which the eukaryotic C₂H₂ zinc-finger structure has evolved. The first prokaryotic zinc-finger domain found in the transcription regulator Ros was identified in Agrobacterium tumefaciens. In the past decades, a large body of evidence revealed Ros/MucR as pleiotropic transcriptional regulators that mainly act as repressors through oligomerization and binding to AT-rich target promoters. The N-terminal domain and the zinc-finger-bearing C-terminal region of these regulatory proteins are engaged in oligomerization and DNA binding, respectively. These properties of the Ros/MucR proteins are similar to those of xenogeneic silencers, such as H-NS, MvaT, and Lsr2, which are mainly found in other lineages. In fact, a novel functional model recently proposed for this protein family suggests that they act as H-NS-'like' gene silencers. The prokaryotic zinc-finger domain exhibits interesting structural and functional features that are different from that of its eukaryotic counterpart (a βββα topology), as it folds in a significantly larger zinc-binding globular domain (a βββαα topology). Phylogenetic analysis of Ros/MucR homologs suggests an ancestral origin of this type of protein in α-Proteobacteria. Furthermore, multiple duplications and lateral gene transfer events contributing to the diversity and phyletic distribution of these regulatory proteins were found in bacterial genomes.

Coordinated regulation of symbiotic adaptation by NodD proteins and NolA in the type I peanut bradyrhizobial strain Bradyrhizobium zhanjiangense CCBAU51778

Type I peanut bradyrhizobial strains can establish efficient symbiosis in contrast to symbiotic incompatibility induced by type II strains with mung bean. The notable distinction in the two kinds of key symbiosis-related regulators nolA and nodD close to the nodABCSUIJ operon region between these two types of peanut bradyrhizobia was found. Therefore, we determined whether NolA and NodD proteins regulate the symbiotic adaptations of type I strains to different hosts. We found that NodD1-NolA synergistically regulated the symbiosis between the type I strain Bradyrhizobium zhanjiangense CCBAU51778 and mung bean, and NodD1-NodD2 jointly regulated nodulation ability. In contrast, NodD1-NolA coordinately regulated nodulation ability in the CCBAU51778-peanut symbiosis. Meanwhile, NodD1 and NolA collectively contributes to competitive nodule colonization of CCBAU51778 on both hosts. The Fucosylated Nod factors and intact type 3 secretion system (T3SS), rather than extra nodD2 and full-length nolA, were critical for effective symbiosis with mung bean. Unexpectedly, T3SS-related genes were activated by NodD2 but not NodD1. Compared to NodD1 and NodD2, NolA predominantly inhibits exopolysaccharide production by promoting exoR expression. Importantly, this is the first report that NolA regulates rhizobial T3SS-related genes. The coordinated regulation and integration of different gene networks to fine-tune the expression of symbiosis-related genes and other accessory genes by NodD1-NolA might be required for CCBAU51778 to efficiently nodulate peanut. This study shed new light on our understanding of the regulatory roles of NolA and NodD proteins in symbiotic adaptation, highlighting the sophisticated gene networks dominated by NodD1-NolA.

The convergent xenogeneic silencer MucR predisposes α-proteobacteria to integrate AT-rich symbiosis genes

Bacterial adaptation is largely shaped by horizontal gene transfer, xenogeneic silencing mediated by lineage-specific DNA bridgers (H-NS, Lsr2, MvaT and Rok), and various anti-silencing mechanisms. No xenogeneic silencing DNA bridger is known for α-proteobacteria, from which mitochondria evolved. By investigating α-proteobacterium Sinorhizobium fredii, a facultative legume microsymbiont, here we report the conserved zinc-finger bearing MucR as a novel xenogeneic silencing DNA bridger. Self-association mediated by its N-terminal domain (NTD) is required for DNA–MucR–DNA bridging complex formation, maximizing MucR stability, transcriptional silencing, and efficient symbiosis in legume nodules. Essential roles of NTD, CTD (C-terminal DNA-binding domain), or full-length MucR in symbiosis can be replaced by non-homologous NTD, CTD, or full-length protein of H-NS from γ-proteobacterium Escherichia coli, while NTD rather than CTD of Lsr2 from Gram-positive Mycobacterium tuberculosis can replace the corresponding domain of MucR in symbiosis. Chromatin immunoprecipitation sequencing reveals similar recruitment profiles of H-NS, MucR and various functional chimeric xenogeneic silencers across the multipartite genome of S. fredii, i.e. preferring AT-rich genomic islands and symbiosis plasmid with key symbiosis genes as shared targets. Collectively, the convergently evolved DNA bridger MucR predisposed α-proteobacteria to integrate AT-rich foreign DNA including symbiosis genes, horizontal transfer of which is strongly selected in nature.

The zinc-finger bearing xenogeneic silencer MucR in α-proteobacteria balances adaptation and regulatory integrity

Foreign AT-rich genes drive bacterial adaptation to new niches while challenging the existing regulation network. Here we report that MucR, a conserved regulator in α-proteobacteria, balances adaptation and regulatory integrity in Sinorhizobium fredii, a facultative microsymbiont of legumes. Chromatin immunoprecipitation sequencing coupled with transcriptomic data reveal that average transcription levels of both target and non-target genes, under free-living and symbiotic conditions, increase with their conservation levels. Targets involved in environmental adaptation and symbiosis belong to genus or species core and can be repressed or activated by MucR in a condition-dependent manner, implying regulatory integrations. However, most targets are enriched in strain-specific genes of lower expression levels and higher AT%. Within each conservation levels, targets have higher AT% and average transcription levels than non-target genes and can be further up-regulated in the mucR mutant. This is consistent with higher AT% of spacers between -35 and -10 elements of promoters for target genes, which enhances transcription. The MucR recruitment level linearly increases with AT% and the number of a flexible pattern (with periodic repeats of Ts) of target sequences. Collectively, MucR directly represses AT-rich foreign genes with predisposed high transcription potential while progressive erosions of its target sites facilitate regulatory integrations of foreign genes.

Genome-Scale Metabolic Modelling of Lifestyle Changes in Rhizobium leguminosarum

Biological nitrogen fixation in rhizobium-legume symbioses is of major importance for sustainable agricultural practices. To establish a mutualistic relationship with their plant host, rhizobia transition from free-living bacteria in soil to growth down infection threads inside plant roots and finally differentiate into nitrogen-fixing bacteroids. We reconstructed a genome-scale metabolic model for Rhizobium leguminosarum and integrated the model with transcriptome, proteome, metabolome, and gene essentiality data to investigate nutrient uptake and metabolic fluxes characteristic of these different lifestyles. Synthesis of leucine, polyphosphate, and AICAR is predicted to be important in the rhizosphere, while myo-inositol catabolism is active in undifferentiated nodule bacteria in agreement with experimental evidence. The model indicates that bacteroids utilize xylose and glycolate in addition to dicarboxylates, which could explain previously described gene expression patterns. Histidine is predicted to be actively synthesized in bacteroids, consistent with transcriptome and proteome data for several rhizobial species. These results provide the basis for targeted experimental investigation of metabolic processes specific to the different stages of the rhizobium-legume symbioses. IMPORTANCE Rhizobia are soil bacteria that induce nodule formation on plant roots and differentiate into nitrogen-fixing bacteroids. A detailed understanding of this complex symbiosis is essential for advancing ongoing efforts to engineer novel symbioses with cereal crops for sustainable agriculture. Here, we reconstruct and validate a genome-scale metabolic model for Rhizobium leguminosarum bv. viciae 3841. By integrating the model with various experimental data sets specific to different stages of symbiosis formation, we elucidate the metabolic characteristics of rhizosphere bacteria, undifferentiated bacteria inside root nodules, and nitrogen-fixing bacteroids. Our model predicts metabolic flux patterns for these three distinct lifestyles, thus providing a framework for the interpretation of genome-scale experimental data sets and identifying targets for future experimental studies.

Global Transcriptional Repression of Diguanylate Cyclases by MucR1 Is Essential for Sinorhizobium-Soybean Symbiosis

The ubiquitous bacterial second messenger c-di-GMP is intensively studied in pathogens but less so in mutualistic bacteria. Here, we report a genome-wide investigation of functional diguanylate cyclases (DGCs) synthesizing c-di-GMP from two molecules of GTP in Sinorhizobium fredii CCBAU45436, a facultative microsymbiont fixing nitrogen in nodules of diverse legumes, including soybean. Among 25 proteins harboring a putative GGDEF domain catalyzing the biosynthesis of c-di-GMP, eight functional DGCs were identified by heterogenous expression in Escherichia coli in a Congo red binding assay. This screening result was further verified by in vitro enzymatic assay with purified full proteins or the GGDEF domains from representative functional and nonfunctional DGCs. In the same in vitro assay, a functional EAL domain catalyzing the degradation of c-di-GMP into pGpG was identified in a protein that has an inactive GGDEF domain but with an active phosphodiesterase (PDE) function. The identified functional DGCs generally exhibited low transcription levels in soybean nodules compared to free-living cultures, as revealed in transcriptomes. An engineered upregulation of a functional DGC in nodules led to a significant increase of c-di-GMP level and symbiotic defects, which were not observed when a functional EAL domain was upregulated at the same level. Further transcriptional analysis and gel shift assay demonstrated that these functional DGCs were all transcriptionally repressed in nodules by a global pleiotropic regulator, MucR1, that is essential in Sinorhizobium-soybean symbiosis. These findings shed novel insights onto the systematic regulation of c-di-GMP biosynthesis in mutualistic symbiosis. IMPORTANCE The ubiquitous second messenger c-di-GMP is well-known for its role in biofilm formation and host adaptation of pathogens, whereas it is less investigated in mutualistic symbioses. Here, we reveal a cocktail of eight functional diguanylate cyclases (DGCs) catalyzing the biosynthesis of c-di-GMP in a broad-host-range Sinorhizobium that can establish nitrogen-fixing nodules on soybean and many other legumes. These functional DGCs are generally transcribed at low levels in soybean nodules compared to free-living conditions. The engineered nodule-specific upregulation of DGC can elevate the c-di-GMP level and cause symbiotic defects, while the upregulation of a phosphodiesterase that quenches c-di-GMP has no detectable symbiotic defects. Moreover, eight functional DGCs located on two different replicons are all directly repressed in nodules by a global silencer, MucR1, that is essential for Sinorhizobium-soybean symbiosis. These findings represent a novel mechanism of a strategic regulation of the c-di-GMP biosynthesis arsenal in prokaryote-eukaryote interactions.

Water-Soluble Humic Materials Modulating Metabolism and Triggering Stress Defense in Sinorhizobium fredii

Bacteria have evolved a series of mechanisms to maintain their survival and reproduction in changeable and stressful environments. In-depth understanding of these mechanisms can allow for better developing and utilizing of bacteria with various biological functions. In this study, we found that water-soluble humic materials (WSHM), a well-known environment-friendly plant growth biostimulant, significantly promoted the free-living growth and survival of Sinorhizobium fredii CCBAU45436 in a bell-shaped, dose-dependent manner, along with more-efficient carbon source consumption and relief of medium acidification. By using RNA-Seq analysis, a total of 1,136 genes significantly up-/downregulated by external addition of WSHM were identified under test conditions. These differentially expressed genes (DEGs) were enriched in functional categories related to carbon/nitrogen metabolism, cellular stress response, and genetic information processing. Further protein-protein interaction (PPI) network analysis and reverse genetic engineering indicated that WSHM might reprogram the transcriptome through inhibiting the expression of key hub gene rsh, which encodes a bifunctional enzyme catalyzing synthesis and hydrolysis of the "magic spot" (p)ppGpp. In addition, the root colonization and viability in soil of S. fredii CCBAU45436 were increased by WSHM. These findings provide us with new insights into how WSHM benefit bacterial adaptations and demonstrate great application value to be a unique inoculant additive. IMPORTANCE Sinorhizobium fredii CCBAU45436 is a highly effective, fast-growing rhizobium that can establish symbiosis with multiple soybean cultivars. However, it is difficult to maintain the high-density effective viable cells in the rhizobial inoculant for the stressful conditions during production, storage, transport, and application. Here, we showed that WSHM greatly increased the viable cells of S. fredii CCBAU45436 in culture, modulating metabolism and triggering stress defense. The root colonization and viability in soil of S. fredii CCBAU45436 were also increased by WSHM. Our results shed new insights into the effects of WSHM on bacteria and the importance of metabolism and stress defense during the bacteria's whole life. In addition, the functional mechanism of WSHM may provide candidate genes for improving environmental adaptability and application potential of bacteria through genetic engineering.

Differentially expressed microRNAs that target functional genes in mature soybean nodules

MicroRNAs (miRNAs) are important regulators of biological functions in plants. To find out what roles miRNAs play in regulating symbiotic nitrogen fixation (SNF) in soybean [Glycine max (L.) Merr.], we identified high-confidence differentially expressed (DE) miRNAs from uninoculated roots (UR), rhizobium-inoculated roots (IR), and nodules (NODs) of soybean by robust small RNA sequencing (sRNA-seq). Based on their predicted target messenger RNAs (mRNAs), the expression profiles of some of these DE miRNAs could be linked to nodule functions. In particular, several miRNAs associated with nutrient transportation genes were differentially expressed in IRs and mature NODs. MiR399b, specifically, was highly induced in IRs and NODs, as well as by inorganic phosphate (Pi) starvation. In composite soybean plants overexpressing miR399b, PHOSPHATE2 (PHO2), a known target of miR399b that inhibits the activities of high-affinity Pi transporters, was strongly repressed. In addition, the overexpression of miR399b in the roots of transgenic composite plants significantly improved whole-plant Pi and ureide concentrations and the overall growth in terms of leaf node numbers and whole-plant dry weight. Our findings suggest that the induction of miR399b in NODs could enhance nitrogen fixation and soybean growth, possibly via improving Pi uptake to achieve a better Pi-nitrogen balance to promote SNF in nodules.

Lineage-Specific Rewiring of Core Pathways Predating Innovation of Legume Nodules Shapes Symbiotic Efficiency

The interkingdom coevolution innovated the rhizobium-legume symbiosis. The application of this nitrogen-fixing system in sustainable agriculture is usually impeded by incompatible interactions between partners. However, the progressive evolution of rhizobium-legume compatibility remains elusive. In this work, deletions of rhcV encoding a structural component of the type three secretion system allow related Sinorhizobium strains to nodulate a previously incompatible soybean cultivar (Glycine max). These rhcV mutants show low to medium to high symbiotic efficiency on the same cultivated soybean while being indistinguishable on wild soybean plants (Glycine soja). The dual pantranscriptomics reveals nodule-specific activation of core symbiosis genes of Sinorhizobium and Glycine genes associated with genome duplication events along the chronogram. Unexpectedly, symbiotic efficiency is in line with lineage-dependent transcriptional profiles of core pathways which predate the diversification of Fabaceae and Sinorhizobium. This is supported by further physiological and biochemical experiments. Particularly, low-efficiency nodules show disordered antioxidant activity and low-energy status, which restrict nitrogen fixation activity. Collectively, the ancient core pathways play a crucial role in optimizing the function of later-evolved mutualistic arsenals in the rhizobium-legume coevolution.

IMPORTANCE Significant roles of complex extracellular microbiota in environmental adaptation of eukaryotes in ever-changing circumstances have been revealed. Given the intracellular infection ability, facultative endosymbionts can be considered pioneers within complex extracellular microbiota and are ideal organisms for understanding the early stage of interkingdom adaptation. This work reveals that the later innovation of key symbiotic arsenals and the lineage-specific network rewiring in ancient core pathways, predating the divergence of legumes and rhizobia, underline the progressive evolution of rhizobium-legume compatibility. This insight not only is significant for improving the application benefits of rhizobial inoculants in sustainable agriculture but also advances our general understanding of the interkingdom coevolution which is theoretically explored by all host-microbiota interactions.

Rhizobiales-Specific RirA Represses a Naturally "Synthetic" Foreign Siderophore Gene Cluster To Maintain Sinorhizobium-Legume Mutualism

Iron homeostasis is strictly regulated in cellular organisms. The Rhizobiales order enriched with symbiotic and pathogenic bacteria has evolved a lineage-specific regulator, RirA, responding to iron fluctuations. However, the regulatory role of RirA in bacterium-host interactions remains largely unknown. Here, we report that RirA is essential for mutualistic interactions of Sinorhizobium fredii with its legume hosts by repressing a gene cluster directing biosynthesis and transport of petrobactin siderophore. Genes encoding an inner membrane ABC transporter (fat) and the biosynthetic machinery (asb) of petrobactin siderophore are sporadically distributed in Gram-positive and Gram-negative bacteria. An outer membrane siderophore receptor gene (fprA) was naturally assembled with asb and fat, forming a long polycistron in S. fredii. An indigenous regulation cascade harboring an inner membrane protease (RseP), a sigma factor (FecI), and its anti-sigma protein (FecR) were involved in direct activation of the fprA-asb-fat polycistron. Operons harboring fecI and fprA-asb-fat, and those encoding the indigenous TonB-ExbB-ExbD complex delivering energy to the outer membrane transport activity, were directly repressed by RirA under iron-replete conditions. The rirA deletion led to upregulation of these operons and iron overload in nodules, impaired intracellular persistence, and symbiotic nitrogen fixation of rhizobia. Mutualistic defects of the rirA mutant can be rescued by blocking activities of this naturally "synthetic" circuit for siderophore biosynthesis and transport. These findings not only are significant for understanding iron homeostasis of mutualistic interactions but also provide insights into assembly and integration of foreign machineries for biosynthesis and transport of siderophores, horizontal transfer of which is selected in microbiota. IMPORTANCE Iron is a public good explored by both eukaryotes and prokaryotes. The abundant ferric form is insoluble under neutral and basic pH conditions, and many bacteria secrete siderophores forming soluble ferric siderophore complexes, which can be then taken up by specific receptors and transporters. Siderophore biosynthesis and uptake machineries can be horizontally transferred among bacteria in nature. Despite increasing attention on the importance of siderophores in host-microbiota interactions, the regulatory integration process of transferred siderophore biosynthesis and transport genes is poorly understood in an evolutionary context. By focusing on the mutualistic rhizobium-legume symbiosis, here, we report how a naturally synthetic foreign siderophore gene cluster was integrated with the rhizobial indigenous regulation cascade, which is essential for maintaining mutualistic interactions.

Comparative Analysis of Core and Accessory Genes in Coexpression Network

Prokaryotes harbor a various proportion of accessory genes in their genomes. The integration of accessory functions with the core regulation network is critical for environmental adaptation, particularly considering a theoretically unlimited number of niches on the earth for microorganisms. Comparative genomics can reveal a co-occurrence pattern between a subset of accessory genes (or variations in core genes) and an adaptation trait, while comparative transcriptomics can further uncover whether a coordinated regulation of gene expression is involved. In this chapter, we introduce a protocol for weighted gene coexpression network construction by using well-developed open source tools, and a further application of such a network in comparative analysis of bacterial core and accessory genes.

The Location of Substitutions and Bacterial Genome Arrangements

Increasing evidence supports the notion that different regions of a genome have unique rates of molecular change. This variation is particularly evident in bacterial genomes where previous studies have reported gene expression and essentiality tend to decrease, whereas substitution rates usually increase with increasing distance from the origin of replication. Genomic reorganization such as rearrangements occur frequently in bacteria and allow for the introduction and restructuring of genetic content, creating gradients of molecular traits along genomes. Here, we explore the interplay of these phenomena by mapping substitutions to the genomes of Escherichia coli, Bacillus subtilis, Streptomyces, and Sinorhizobium meliloti, quantifying how many substitutions have occurred at each position in the genome. Preceding work indicates that substitution rate significantly increases with distance from the origin. Using a larger sample size and accounting for genome rearrangements through ancestral reconstruction, our analysis demonstrates that the correlation between the number of substitutions and the distance from the origin of replication is significant but small and inconsistent in direction. Some replicons had a significantly decreasing trend (E. coli and the chromosome of S. meliloti), whereas others showed the opposite significant trend (B. subtilis, Streptomyces, pSymA and pSymB in S. meliloti). dN, dS, and ω were examined across all genes and there was no significant correlation between those values and distance from the origin. This study highlights the impact that genomic rearrangements and location have on molecular trends in some bacteria, illustrating the importance of considering spatial trends in molecular evolutionary analysis. Assuming that molecular trends are exclusively in one direction can be problematic.

Ancestral zinc-finger bearing protein MucR in alpha-proteobacteria: A novel xenogeneic silencer?

The MucR/Ros family protein is conserved in alpha-proteobacteria and characterized by its zinc-finger motif that has been proposed as the ancestral domain from which the eukaryotic C2H2 zinc-finger structure evolved. In the past decades, accumulated evidences have revealed MucR as a pleiotropic transcriptional regulator that integrating multiple functions such as virulence, symbiosis, cell cycle and various physiological processes. Scattered reports indicate that MucR mainly acts as a repressor, through oligomerization and binding to multiple sites of AT-rich target promoters. The N-terminal region and zinc-finger bearing C-terminal region of MucR mediate oligomerization and DNA-binding, respectively. These features are convergent to those of xenogeneic silencers such as H-NS, MvaT, Lsr2 and Rok, which are mainly found in other lineages. Phylogenetic analysis of MucR homologs suggests an ancestral origin of MucR in alpha- and delta-proteobacteria. Multiple independent duplication and lateral gene transfer events contribute to the diversity and phyletic distribution of MucR. Finally, we posed questions which remain unexplored regarding the putative roles of MucR as a xenogeneic silencer and a general manager in balancing adaptation and regulatory integration in the pangenome context.

The histone modification H3K4me3 marks functional genes in soybean nodules

Nitrogen fixation in legumes requires the development of specialized organs called root nodules. Here we characterized the high-confidence transcriptome and genome-wide patterns of H3K4me3 marks in soybean roots and mature nodules symbiotic with Sinorhizobium fredii. Changes in H3K4me3 levels were positively associated with the transcription levels of functional genes in the nodules. The up-regulation of H3K4me3 levels was not only present in leghaemoglobin and nodulin-related genes, but also in genes involved in nitrogen and carbon metabolic pathways. In addition, genes regulating the transmembrane transport of metal ions, phosphates, sulphates, peptides, and sugars were differentially modified. On the contrary, a loss of H3K4me3 marks was found in several key transcription factor genes and was correlated with the down-regulation of the defense-related network in nodules, which could contribute to nodule maintenance. All these findings demonstrate massive reprogramming of gene expressions via alterations in H3K4me3 levels in the genes in mature soybean nodules.

Experimental Evolution of Legume Symbionts: What Have We Learnt?

Rhizobia, the nitrogen-fixing symbionts of legumes, are polyphyletic bacteria distributed in many alpha- and beta-proteobacterial genera. They likely emerged and diversified through independent horizontal transfers of key symbiotic genes. To replay the evolution of a new rhizobium genus under laboratory conditions, the symbiotic plasmid of Cupriavidus taiwanensis was introduced in the plant pathogen Ralstonia solanacearum, and the generated proto-rhizobium was submitted to repeated inoculations to the C. taiwanensis host, Mimosa pudica L.. This experiment validated a two-step evolutionary scenario of key symbiotic gene acquisition followed by genome remodeling under plant selection. Nodulation and nodule cell infection were obtained and optimized mainly via the rewiring of regulatory circuits of the recipient bacterium. Symbiotic adaptation was shown to be accelerated by the activity of a mutagenesis cassette conserved in most rhizobia. Investigating mutated genes led us to identify new components of R. solanacearum virulence and C. taiwanensis symbiosis. Nitrogen fixation was not acquired in our short experiment. However, we showed that post-infection sanctions allowed the increase in frequency of nitrogen-fixing variants among a non-fixing population in the M. pudica-C. taiwanensis system and likely allowed the spread of this trait in natura. Experimental evolution thus provided new insights into rhizobium biology and evolution.

Modulation of Symbiotic Compatibility by Rhizobial Zinc Starvation Machinery

The rhizobium-legume symbiosis contributes around 65% of biological nitrogen fixation in agriculture systems and is critical for sustainable agriculture by reducing the amount of chemical nitrogen fertilizer being used. Rhizobial inocula have been commercialized for more than 100 years, but the efficiency of inoculation can vary among legume cultivars, field sites, and years. These long-lasting challenging problems impede the establishment of a sustainable agriculture, particularly in developing countries. Here, we report that rhizobial zinc starvation machinery containing a conserved high-affinity zinc transporter and accessory components makes cumulative contributions to modulating rhizobial symbiotic compatibility. This work highlights a critical role of largely unexplored nutritional immunity in the rhizobium-legume symbiosis, which makes zinc starvation machinery an attractive target for improving rhizobial symbiotic compatibility.

Pathogenic bacteria need high-affinity zinc uptake systems to counteract the nutritional immunity exerted by infected hosts. However, our understanding of zinc homeostasis in mutualistic systems such as the rhizobium-legume symbiosis is limited. Here, we show that the conserved high-affinity zinc transporter ZnuABC and accessory transporter proteins (Zip1, Zip2, and c06450) made cumulative contributions to nodulation of the broad-host-range strain Sinorhizobium fredii CCBAU45436. Zur acted as a zinc-dependent repressor for the znuC-znuB-zur operon, znuA, and c06450 by binding to the associated Zur box, but did not regulate zip1 and zip2. ZnuABC was the major zinc transporter. Combined mutants lacking znuA and one of the three accessory genes had more severe defects in nodulation and growth under zinc starvation conditions than the znuA mutant, though rhizoplane colonization by these mutants was not impaired. In contrast to the elite strain CCBAU45436, more drastic symbiotic defects were observed for the znuA mutants of other Sinorhizobium strains, which lack at least one of the three accessory genes in their genomes and are characterized by their limited host range and geographical distribution. The znu-derived mutants showed a higher expression level of nod genes involved in Nod factor biosynthesis and a reduced expression of genes encoding a type three secretion system and its effector NopP, which can interfere with the host immune system. Application of exogenous zinc restored the nodulation ability of these znu-derived mutants. Therefore, the conserved ZnuABC and accessory components in the zinc starvation machinery play an important role in modulating symbiotic compatibility.

Metabolic Analyses of Nitrogen Fixation in the Soybean Microsymbiont Sinorhizobium fredii Using Constraint-Based Modeling

Nitrogen is the most limiting macronutrient for plant growth, and rhizobia are important bacteria for agriculture because they can fix atmospheric nitrogen and make it available to legumes through the establishment of a symbiotic relationship with their host plants. In this work, we studied the nitrogen fixation process in the microsymbiont Sinorhizobium fredii at the genome level. A metabolic model was built using genome annotation and literature to reconstruct the symbiotic form of S. fredii. Genes controlling the nitrogen fixation process were identified by simulating gene knockouts. Additionally, the nitrogen-fixing capacities of S. fredii CCBAU45436 in symbiosis with cultivated and wild soybeans were evaluated. The predictions suggested an outperformance of S. fredii with cultivated soybean, consistent with published experimental evidence. The reconstruction presented here will help to understand and improve nitrogen fixation capabilities of S. fredii and will be beneficial for agriculture by reducing the reliance on fertilizer applications.

Rhizobia are soil bacteria able to establish symbiosis with diverse host plants. Specifically, Sinorhizobium fredii is a soil bacterium that forms nitrogen-fixing root nodules in diverse legumes, including soybean. The strain S. fredii CCBAU45436 is a dominant sublineage of S. fredii that nodulates soybeans in alkaline-saline soils in the Huang-Huai-Hai Plain region of China. Here, we present a manually curated metabolic model of the symbiotic form of Sinorhizobium fredii CCBAU45436. A symbiosis reaction was defined to describe the specific soybean-microsymbiont association. The performance and quality of the reconstruction had a 70% score when assessed using a standardized genome-scale metabolic model test suite. The model was used to evaluate in silico single-gene knockouts to determine the genes controlling the nitrogen fixation process. One hundred forty-one of 541 genes (26%) were found to influence the symbiotic process, wherein 121 genes were predicted as essential and 20 others as having a partial effect. Transcriptomic profiles of CCBAU45436 were used to evaluate the nitrogen fixation capacity in cultivated versus in wild soybean inoculated with the microsymbiont. The model quantified the nitrogen fixation activities of the strain in these two hosts and predicted a higher nitrogen fixation capacity in cultivated soybean. Our results are consistent with published data demonstrating larger amounts of ureides and total nitrogen in cultivated soybean than in wild soybean. This work presents the first metabolic network reconstruction of S. fredii as an example of a useful tool for exploring the potential benefits of microsymbionts to sustainable agriculture and the ecosystem.

IMPORTANCE Nitrogen is the most limiting macronutrient for plant growth, and rhizobia are important bacteria for agriculture because they can fix atmospheric nitrogen and make it available to legumes through the establishment of a symbiotic relationship with their host plants. In this work, we studied the nitrogen fixation process in the microsymbiont Sinorhizobium fredii at the genome level. A metabolic model was built using genome annotation and literature to reconstruct the symbiotic form of S. fredii. Genes controlling the nitrogen fixation process were identified by simulating gene knockouts. Additionally, the nitrogen-fixing capacities of S. fredii CCBAU45436 in symbiosis with cultivated and wild soybeans were evaluated. The predictions suggested an outperformance of S. fredii with cultivated soybean, consistent with published experimental evidence. The reconstruction presented here will help to understand and improve nitrogen fixation capabilities of S. fredii and will be beneficial for agriculture by reducing the reliance on fertilizer applications.

Pulse Crop Genetics for a Sustainable Future: Where We Are Now and Where We Should Be Heading

The last decade has witnessed dramatic changes in global food consumption patterns mainly because of population growth and economic development. Food substitutions for healthier eating, such as swapping regular servings of meat for protein-rich crops, is an emerging diet trend that may shape the future of food systems and the environment worldwide. To meet the erratic consumer demand in a rapidly changing world where resources become increasingly scarce due largely to anthropogenic activity, the need to develop crops that benefit both human health and the environment has become urgent. Legumes are often considered to be affordable plant-based sources of dietary proteins. Growing legumes provides significant benefits to cropping systems and the environment because of their natural ability to perform symbiotic nitrogen fixation, which enhances both soil fertility and water-use efficiency. In recent years, the focus in legume research has seen a transition from merely improving economically important species such as soybeans to increasingly turning attention to some promising underutilized species whose genetic resources hold the potential to address global challenges such as food security and climate change. Pulse crops have gained in popularity as an affordable source of food or feed; in fact, the United Nations designated 2016 as the International Year of Pulses, proclaiming their critical role in enhancing global food security. Given that many studies have been conducted on numerous underutilized pulse crops across the world, we provide a systematic review of the related literature to identify gaps and opportunities in pulse crop genetics research. We then discuss plausible strategies for developing and using pulse crops to strengthen food and nutrition security in the face of climate and anthropogenic changes.

Pangenomics Comes of Age: From Bacteria to Plant and Animal Applications

The pangenome refers to a collection of genomic sequence found in the entire species or population rather than in a single individual; the sequence can be core, present in all individuals, or accessory (variable or dispensable), found in a subset of individuals only. While pangenomic studies were first undertaken in bacterial species, developments in genome sequencing and assembly approaches have allowed construction of pangenomes for eukaryotic organisms, fungi, plants, and animals, including two large-scale human pangenome projects. Analysis of the these pangenomes revealed key differences, most likely stemming from divergent evolutionary histories, but also surprising similarities.

High-Throughput Mass Spectrometric Analysis of the Whole Proteome and Secretome From Sinorhizobium fredii Strains CCBAU25509 and CCBAU45436

Sinorhizobium fredii is a dominant rhizobium on alkaline-saline land that can induce nitrogen-fixing symbiotic root nodules in soybean. Two S. fredii strains, CCBAU25509 and CCBAU45436, were used in this study to facilitate in-depth analyses of this species and its interactions with soybean. We have previously completed the full assembly of the genomes and detailed transcriptomic analyses for these two S. fredii strains, CCBAU25509 and CCBAU45436, that exhibit differential compatibility toward some soybean hosts. In this work, we performed high-throughput Orbitrap analyses of the whole proteomes and secretomes of CCBAU25509 and CCBAU45436 at different growth stages. Our proteomic data cover coding sequences in the chromosome, chromid, symbiotic plasmid, and other accessory plasmids. In general, we found higher levels of protein expression by genes in the chromosomal genome, whereas proteins encoded by the symbiotic plasmid were differentially accumulated in bacteroids. We identified secreted proteins from the extracellular medium, including seven and eight Nodulation Outer Proteins (Nops) encoded by the symbiotic plasmid of CCBAU25509 and CCBAU45436, respectively. Differential host restriction of CCBAU25509 and CCBAU45436 is regulated by the allelic type of the soybean Rj2(Rfg1) protein. Using sequencing data from this work and available in public databases, our analysis confirmed that the soybean Rj2(Rfg1) protein has three major allelic types (Rj2/rfg1, rj2/Rfg1, rj2/rfg1) that determine the host restriction of some Bradyrhizobium diazoefficiens and S. fredii strains. A mutant defective in the type 3 protein secretion system (T3SS) in CCBAU25509 allowed this strain to nodulate otherwise-incompatible soybeans carrying the rj2/Rfg1 allelic type, probably by disrupting Nops secretion. The allelic forms of NopP and NopI in S. fredii might be associated with the restriction imposed by Rfg1. By swapping the NopP between CCBAU25509 and CCBAU45436, we found that only the strains carrying NopP from CCBAU45436 could nodulate soybeans carrying the rj2/Rfg1 allelic type. However, no direct interaction between either forms of NopP and Rfg1 could be observed.

Coordinated Regulation of the Size and Number of Polyhydroxybutyrate Granules by Core and Accessory Phasins in the Facultative Microsymbiont Sinorhizobium fredii NGR234

The exact roles of various granule-associated proteins (GAPs) of polyhydroxybutyrate (PHB) are poorly investigated, particularly for bacteria associated with plants. In this study, four structural GAPs, named phasins PhaP1 to PhaP4, were identified and demonstrated as true phasins colocalized with PHB granules in Sinorhizobium fredii NGR234, a facultative microsymbiont of Vigna unguiculata and many other legumes. The conserved PhaP2 dominated in regulation of granule size under both free-living and symbiotic conditions. PhaP1, another conserved phasin, made a higher contribution than accessory phasins PhaP4 and PhaP3 to PHB biosynthesis at stationary phase. PhaP3, with limited phyletic distribution on the symbiosis plasmid of Sinorhizobium, was more important than PhaP1 in regulating PHB biosynthesis in V. unguiculata nodules. Under the test conditions, no significant symbiotic defects were observed for mutants lacking individual or multiple phaP genes. The mutant lacking two PHB synthases showed impaired symbiotic performance, while mutations in individual PHB synthases or a PHB depolymerase yielded no symbiotic defects. This phenomenon is not related to either the number or size of PHB granules in test mutants within nodules. Distinct metabolic profiles and cocktail pools of GAPs of different phaP mutants imply that core and accessory phasins can be differentially involved in regulating other cellular processes in the facultative microsymbiont S. fredii NGR234.IMPORTANCE Polyhydroxybutyrate (PHB) granules are a store of carbon and energy in bacteria and archaea and play an important role in stress adaptation. Recent studies have highlighted distinct roles of several granule-associated proteins (GAPs) in regulating the size, number, and localization of PHB granules in free-living bacteria, though our knowledge of the role of GAPs in bacteria associated with plants is still limited. Here we report distinct roles of core and accessory phasins associated with PHB granules of Sinorhizobium fredii NGR234, a broad-host-range microsymbiont of diverse legumes. Core phasins PhaP2 and PhaP1 are conserved major phasins in free-living cells. PhaP2 and accessory phasin PhaP3, encoded by an auxiliary gene on the symbiosis plasmid, are major phasins in nitrogen-fixing bacteroids in cowpea nodules. GAPs and metabolic profiles can vary in different phaP mutants. Contrasting symbiotic performances between mutants lacking PHB synthases, depolymerase, or phasins were revealed.