Comparative genomics of the core and accessory genomes of 48 Sinorhizobium strains comprising five genospecies

Horizontal gene transfer of the Mer operon is associated with large effects on the transcriptome and increased tolerance to mercury in nitrogen-fixing bacteria

BACKGROUND

Mercury (Hg) is highly toxic and has the potential to cause severe health problems for humans and foraging animals when transported into edible plant parts. Soil rhizobia that form symbiosis with legumes may possess mechanisms to prevent heavy metal translocation from roots to shoots in plants by exporting metals from nodules or compartmentalizing metal ions inside nodules. Horizontal gene transfer has potential to confer immediate de novo adaptations to stress. We used comparative genomics of high quality de novo assemblies to identify structural differences in the genomes of nitrogen-fixing rhizobia that were isolated from a mercury (Hg) mine site that show high variation in their tolerance to Hg.

RESULTS

Our analyses identified multiple structurally conserved merA homologs in the genomes of Sinorhizobium medicae and Rhizobium leguminosarum but only the strains that possessed a Mer operon exhibited 10-fold increased tolerance to Hg. RNAseq analysis revealed nearly all genes in the Mer operon were significantly up-regulated in response to Hg stress in free-living conditions and in nodules. In both free-living and nodule environments, we found the Hg-tolerant strains with a Mer operon exhibited the fewest number of differentially expressed genes (DEGs) in the genome, indicating a rapid and efficient detoxification of Hg from the cells that reduced general stress responses to the Hg-treatment. Expression changes in S. medicae while in bacteroids showed that both rhizobia strain and host-plant tolerance affected the number of DEGs. Aside from Mer operon genes, nif genes which are involved in nitrogenase activity in S. medicae showed significant up-regulation in the most Hg-tolerant strain while inside the most Hg-accumulating host-plant. Transfer of a plasmid containing the Mer operon from the most tolerant strain to low-tolerant strains resulted in an immediate increase in Hg tolerance, indicating that the Mer operon is able to confer hyper tolerance to Hg.

CONCLUSIONS

Mer operons have not been previously reported in nitrogen-fixing rhizobia. This study demonstrates a pivotal role of the Mer operon in effective mercury detoxification and hypertolerance in nitrogen-fixing rhizobia. This finding has major implications not only for soil bioremediation, but also host plants growing in mercury contaminated soils.

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.

The impact of Elaeagnus angustifolia root exudates on Parafrankia soli NRRL B-16219 exoproteome

Root exudates from host plant species are known to play a critical role in the establishment and maintenance of symbiotic relationships with soil bacteria. In this study, we investigated the impact of root exudates from compatible host plant species; Elaeagnus angustifolia on the exoproteome of Parafrankia soli strain NRRL B-16219. A total of 565 proteins were evidenced as differentially abundant, with 32 upregulated and 533 downregulated in presence of the plant exudates. Analysis of the function of these proteins suggests that the bacterial strain is undergoing a complex metabolic reprogramming towards a new developmental phase elicited in presence of host plant root exudates. The upregulation of Type II/IV secretion system proteins among the differentially expressed proteins indicates their possible role in infecting the host plant, as shown for some rhizobia. Additionally, EF-Tu, proteins upregulated in this study, may function as an effector for the T4SSs and trigger plant defense responses. These findings suggest that Parafrankia soli may use EF-Tu to infect the actinorhizal host plant and pave the way for further investigations of the molecular mechanisms underlying the establishment of symbiotic relationships.

Omics technology draws a comprehensive heavy metal resistance strategy in bacteria

The rapid industrial revolution significantly increased heavy metal pollution, becoming a major global environmental concern. This pollution is considered as one of the most harmful and toxic threats to all environmental components (air, soil, water, animals, and plants until reaching to human). Therefore, scientists try to find a promising and eco-friendly technique to solve this problem i.e., bacterial bioremediation. Various heavy metal resistance mechanisms were reported. Omics technologies can significantly improve our understanding of heavy metal resistant bacteria and their communities. They are a potent tool for investigating the adaptation processes of microbes in severe conditions. These omics methods provide unique benefits for investigating metabolic alterations, microbial diversity, and mechanisms of resistance of individual strains or communities to harsh conditions. Starting with genome sequencing which provides us with complete and comprehensive insight into the resistance mechanism of heavy metal resistant bacteria. Moreover, genome sequencing facilitates the opportunities to identify specific metal resistance genes, operons, and regulatory elements in the genomes of individual bacteria, understand the genetic mechanisms and variations responsible for heavy metal resistance within and between bacterial species in addition to the transcriptome, proteome that obtain the real expressed genes. Moreover, at the community level, metagenome, meta transcriptome and meta proteome participate in understanding the microbial interactive network potentially novel metabolic pathways, enzymes and gene species can all be found using these methods. This review presents the state of the art and anticipated developments in the use of omics technologies in the investigation of microbes used for heavy metal bioremediation.

Guidelines for the description of rhizobial symbiovars

Rhizobia are bacteria that form nitrogen-fixing nodules in legume plants. The sets of genes responsible for both nodulation and nitrogen fixation are carried in plasmids or genomic islands that are often mobile. Different strains within a species sometimes have different host specificities, while very similar symbiosis genes may be found in strains of different species. These specificity variants are known as symbiovars, and many of them have been given names, but there are no established guidelines for defining or naming them. Here, we discuss the requirements for guidelines to describe symbiovars, propose a set of guidelines, provide a list of all symbiovars for which descriptions have been published so far, and offer a mechanism to maintain a list in the future.

The Type IV Secretion System (T4SS) Mediates Symbiosis between Bradyrhizobium sp. SUTN9-2 and Legumes

There has been little study of the type IV secretion system (T4SS) of bradyrhizobia and its role in legume symbiosis. Therefore, broad host range Bradyrhizobium sp. SUTN9-2 was selected for study. The chromosome of Bradyrhizobium sp. SUTN9-2 contains two copies of the T4SS gene, homologous with the tra/trb operons. A phylogenetic tree of the T4SS gene traG was constructed, which exemplified its horizontal transfer among Bradyrhizobium and Mesorhizobium genera. They also showed similar gene arrangements for the tra/trb operons. However, the virD2 gene was not observed in Mesorhizobium, except M. oppotunistum WSM2075. Interestingly, the orientation of copG, traG, and virD2 cluster was unique to the Bradyrhizobium genus. The phylogenetic tree of copG, traG, and virD2 demonstrated that copies 1 and 2 of these genes were grouped in different clades. In addition, the derived mutant and complementation strains of T4SS were investigated in representative legumes Genistoids, Dalbergioids, and Millettiods. When T4SS copy 1 (T4SS1) was deleted, the nodule number and nitrogenase activity decreased. This supports a positive effect of T4SS1 on symbiosis. In addition, delayed nodulation was observed 7 dpi, which was restored by the complementation of T4SS1. Therefore, T4SS plays an important role in the symbiotic interaction between Bradyrhizobium sp. SUTN9-2 and its leguminous hosts. IMPORTANCE SUTN9-2 is a broad host range strain capable of symbiosis with several legumes. Two copies of T4SS clusters belonging to the tra/trb operon are observed on chromosomes with different gene arrangements. We use phylogenetic tree and gene annotation analysis to predict the evolution of the tra/trb operon of rhizobia. Our finding suggests that the gene encoding the T4SS gene among Bradyrhizobium and Mesorhizobium may have coevolution. In addition, Bradyrhizobium has a uniquely arranged copG, traG, and virD2 gene cluster. The results of T4SS1 gene deletion and complementation revealed its positive effect on nodulation. Therefore, T4SS seems to be another determinant for symbiosis. This is the first report on the role of T4SS in Bradyrhizobium symbiosis.

Scrutiny of NolA and NodD1 Regulatory Roles in Symbiotic Compatibility Unveils New Insights into Bradyrhizobium guangxiense CCBAU53363 Interacting with Peanut (Arachis hypogaea) and Mung Bean (Vigna radiata)

Bradyrhizobium guangxiense CCBAU53363 efficiently nodulates peanut but exhibits incompatible interaction with mung bean. By comparing the common nod region with those of other peanut bradyrhizobia efficiently nodulating these two hosts, distinctive characteristics with a single nodD isoform (nodD1) and a truncated nolA were identified. However, the regulatory roles of NodD1 and NolA and their coordination in legume-bradyrhizobial interactions remain largely unknown in terms of explaining the contrasting symbiotic compatibility. Here, we report that nolA was important for CCBAU53363 symbiosis with peanut but restricted nodulation on mung bean, while nodD1 was dispensable for CCBAU53363 symbiosis with peanut but essential for nodulation on mung bean. Moreover, nolA exerted a cumulative contribution with nodD1 to efficient symbiosis with peanut. Additionally, mutants lacking nolA delayed nodulation on peanut, and both nolA and nodD1 were required for competitive nodule colonization. It is noteworth that most of the nodulation genes and type III secretion system (T3SS)-related genes were significantly downregulated in a strain 53ΔnodD1nolA mutant compared to wild-type strain CCBAU53363, and the downregulated nodulation genes also had a greater impact than T3SS-related genes on the symbiotic defect of 53ΔnodD1nolA on peanut, which was supported by a more severe symbiotic defect induced by 53ΔnodC than that with the 53ΔnodD1nopP, 53ΔnodD1rhcJ, and 53ΔnodD1ttsI mutants. NolA did not regulate nod gene expression but did regulate the T3SS effector gene nopP in an indirect way. Meanwhile, nolA, nodW, and some T3SS-related genes besides nopP were also demonstrated as new "repressors" that seriously impaired CCBAU53363 symbiosis with mung bean. Taken together, the roles and essentiality of nolA and nodD1 in modulating symbiotic compatibility are sophisticated and host dependent. IMPORTANCE The main findings of this study were that we clarified that the roles and essentiality of nodD1 and nolA are host dependent. Importantly, for the first time, NolA was found to positively regulate T3SS effector gene nopP to mediate incompatibility on mung bean. Additionally, NolA does not regulate nod genes, which are activated by NodD1. nolA exerts a cumulative effect with nodD1 on CCBAU53363 symbiosis with peanut. These findings shed new light on our understanding of coordinated regulation of NodD1 and NolA in peanut bradyrhizobia with different hosts.

Diversification of Sinorhizobium populations associated with Medicago polymorpha and Medicago lupulina in purple soil of China

The double selection of environment adaptation and host specificity forced the diversification of rhizobia in nature. In the tropical region of China, Medicago polymorpha and Medicago lupulina are widely distributed, particularly in purple soil. However, the local distribution and diversity of rhizobia associated with these legumes has not been systematically investigated. To this end, root nodules of M. polymorpha and M. lupulina grown in purple soil at seven locations in Yunnan Province of China were collected for rhizobial isolation. The obtained rhizobia were characterized by RFLP of 16S-23S rRNA intergenic spacer, BOXAIR fingerprinting, and phylogeny of housekeeping and symbiosis genes. As result, a total of 91 rhizobial strains were classified into species Sinorhizobium medicae and S. meliloti, while three nodC gene types were identified among them. S. medicae containing nodC of type I was dominant in farmlands associated with M. polymorpha; while S. meliloti harboring nodC of type III was dominant in wild land nodulated by M. lupulina. For both rhizobial species, greater genetic diversity was detected in the populations isolated from their preferred host plant. A high level of genetic differentiation was observed between the two Sinorhizobium species, and gene flow was evident within the populations of the same species derived from different soil types, indicating that rhizobial evolution is likely associated with the soil features. To examine the effects of environmental features on rhizobial distribution, soil physicochemical traits and rhizobial genotypes were applied for constrained analysis of principle coordinates, which demonstrated that soil features like pH, nitrogen and sodium were the principle factors governing the rhizobial geographical distribution. Altogether, both S. medicae and S. meliloti strains could naturally nodulate with M. polymorpha and M. lupulina, but the rhizobium-legume symbiosis compatibility determined by both the host species and soil factors was also highlighted.

Genome-Wide Association Studies across Environmental and Genetic Contexts Reveal Complex Genetic Architecture of Symbiotic Extended Phenotypes

A goal of modern biology is to develop the genotype-phenotype (G→P) map, a predictive understanding of how genomic information generates trait variation that forms the basis of both natural and managed communities. As microbiome research advances, however, it has become clear that many of these traits are symbiotic extended phenotypes, being governed by genetic variation encoded not only by the host's own genome, but also by the genomes of myriad cryptic symbionts. Building a reliable G→P map therefore requires accounting for the multitude of interacting genes and even genomes involved in symbiosis. Here, we use naturally occurring genetic variation in 191 strains of the model microbial symbiont Sinorhizobium meliloti paired with two genotypes of the host Medicago truncatula in four genome-wide association studies (GWAS) to determine the genomic architecture of a key symbiotic extended phenotype-partner quality, or the fitness benefit conferred to a host by a particular symbiont genotype, within and across environmental contexts and host genotypes. We define three novel categories of loci in rhizobium genomes that must be accounted for if we want to build a reliable G→P map of partner quality; namely, (i) loci whose identities depend on the environment, (ii) those that depend on the host genotype with which rhizobia interact, and (iii) universal loci that are likely important in all or most environments. IMPORTANCE Given the rapid rise of research on how microbiomes can be harnessed to improve host health, understanding the contribution of microbial genetic variation to host phenotypic variation is pressing, and will better enable us to predict the evolution of (and select more precisely for) symbiotic extended phenotypes that impact host health. We uncover extensive context-dependency in both the identity and functions of symbiont loci that control host growth, which makes predicting the genes and pathways important for determining symbiotic outcomes under different conditions more challenging. Despite this context-dependency, we also resolve a core set of universal loci that are likely important in all or most environments, and thus, serve as excellent targets both for genetic engineering and future coevolutionary studies of symbiosis.

The Bradyrhizobium Sp. LmicA16 Type VI Secretion System Is Required for Efficient Nodulation of Lupinus Spp

Many bacteria of the genus Bradyrhizobium are capable of inducing nodules in legumes. In this work, the importance of a type VI secretion system (T6SS) in a symbiotic strain of the genus Bradyrhizobium is described. T6SS of Bradyrhizobium sp. LmicA16 (A16) is necessary for efficient nodulation with Lupinus micranthus and Lupinus angustifolius. A mutant in the gene vgrG, coding for a component of the T6SS nanostructure, induced less nodules and smaller plants than the wild-type (wt) strain and was less competitive when co-inoculated with the wt strain. A16 T6SS genes are organized in a 26-kb DNA region in two divergent gene clusters of nine genes each. One of these genes codes for a protein (Tsb1) of unknown function but containing a methyltransferase domain. A tsb1 mutant showed an intermediate symbiotic phenotype regarding vgrG mutant and higher mucoidity than the wt strain in free-living conditions. T6SS promoter fusions to the lacZ reporter indicate expression in nodules but not in free-living cells grown in different media and conditions. The analysis of nodule structure revealed that the level of nodule colonization was significantly reduced in the mutants with respect to the wt strain.

The Rhizobial Type 3 Secretion System: The Dr. Jekyll and Mr. Hyde in the Rhizobium–Legume Symbiosis

Rhizobia are soil bacteria that can establish a symbiotic association with legumes. As a result, plant nodules are formed on the roots of the host plants where rhizobia differentiate to bacteroids capable of fixing atmospheric nitrogen into ammonia. This ammonia is transferred to the plant in exchange of a carbon source and an appropriate environment for bacterial survival. This process is subjected to a tight regulation with several checkpoints to allow the progression of the infection or its restriction. The type 3 secretion system (T3SS) is a secretory system that injects proteins, called effectors (T3E), directly into the cytoplasm of the host cell, altering host pathways or suppressing host defense responses. This secretion system is not present in all rhizobia but its role in symbiosis is crucial for some symbiotic associations, showing two possible faces as Dr. Jekyll and Mr. Hyde: it can be completely necessary for the formation of nodules, or it can block nodulation in different legume species/cultivars. In this review, we compile all the information currently available about the effects of different rhizobial effectors on plant symbiotic phenotypes. These phenotypes are diverse and highlight the importance of the T3SS in certain rhizobium–legume symbioses.

Root exudate chemistry affects soil carbon mobilization via microbial community reassembly

Plant roots are one of the major mediators that allocate carbon captured from the atmosphere to soils as rhizodeposits, including root exudates. Although rhizodeposition regulates both microbial activity and the biogeochemical cycling of nutrients, the effects of particular exudate species on soil carbon fluxes and key rhizosphere microorganisms remain unclear. By combining high-throughput sequencing, q-PCR, and NanoSIMS analyses, we characterized the bacterial community structure, quantified total bacteria depending on root exudate chemistry, and analyzed the consequences on the mobility of mineral-protected carbon. Using well-controlled incubation experiments, we showed that the three most abundant groups of root exudates (amino acids, carboxylic acids, and sugars) have contrasting effects on the release of dissolved organic carbon (DOC) and bioavailable Fe in an Ultisol through the disruption of organo-mineral associations and the alteration of bacterial communities, thus priming organic matter decomposition in the rhizosphere. High resolution (down to 50 nm) NanoSIMS images of mineral particles indicated that iron and silicon co-localized significantly more organic carbon following amino acid inputs than treatments without exudates or with carboxylic acids. The application of sugar strongly reduced microbial diversity without impacting soil carbon mobilization. Carboxylic acids increased the prevalence of Actinobacteria and facilitated carbon mobilization, whereas amino acid addition increased the abundances of Proteobacteria that prevented DOC release. In summary, root exudate functions are defined by their chemical composition that regulates bacterial community composition and, consequently, the biogeochemical cycling of carbon in the rhizosphere.

A newly evolved chimeric lysin motif receptor-like kinase in Medicago truncatula spp. tricycla R108 extends its Rhizobia symbiotic partnership

Rhizobial lipochitooligosaccharidic Nod factors (NFs), specified by nod genes, are the primary determinants of host specificity in the legume-Rhizobia symbiosis. We examined the nodulation ability of Medicago truncatula cv Jemalong A17 and M. truncatula ssp. tricycla R108 with the Sinorhizobium meliloti nodF/nodL mutant, which produces modified NFs. We then applied genetic and functional approaches to study the genetic basis and mechanism of nodulation of R108 by this mutant. We show that the nodF/nodL mutant can nodulate R108 but not A17. Using genomics and reverse genetics, we identified a newly evolved, chimeric LysM receptor-like kinase gene in R108, LYK2bis, which is responsible for the phenotype and can allow A17 to gain nodulation with the nodF/nodL mutant. We found that LYK2bis is involved in nodulation by mutants producing nonO-acetylated NFs and interacts with the key receptor protein NFP. Many, but not all, natural S. meliloti and S. medicae strains tested require LYK2bis for efficient nodulation of R108. Our findings reveal that a newly evolved gene in R108, LYK2bis, extends nodulation specificity to mutants producing nonO-acetylated NFs and is important for nodulation by many natural Sinorhizobia. Evolution of this gene may present an adaptive advantage to allow nodulation by a greater variety of strains.

The changing paradigm of rhizobial taxonomy and its systematic growth upto postgenomic technologies

Rhizobia are a diazotrophic group of bacteria that are usually isolated form the nodules in roots, stem of leguminous plants and are able to form nodules in the host plant owing to the presence of symbiotic genes. The rhizobial community is highly diverse, and therefore, the taxonomy and genera-wise classification of rhizobia has been constantly changing since the last three decades. This is mainly due to technical advancements, and shifts in definitions, resulting in a changing paradigm of rhizobia taxonomy. Initially, the taxonomic definitions at the species and sub species level were based on phylogenetic analysis of 16S rRNA sequence, followed by polyphasic approach to have phenotypic, biochemical, and genetic analysis including multilocus sequence analysis. Rhizobia mainly belonging to α- and β-proteobacteria, and recently new additions from γ-proteobacteria had been classified. Nowadays rhizobial taxonomy has been replaced by genome-based taxonomy that allows gaining more insights of genomic characteristics. These omics-technologies provide genome specific information that considers nodulation and symbiotic genes, along with molecular markers as taxonomic traits. Taxonomy based on complete genome sequence (genotaxonomy), average nucleotide identity, is now being considered as primary approach, resulting in an ongoing paradigm shift in rhizobial taxonomy. Also, pairwise whole-genome comparisons, phylogenomic analyses offer correlations between DNA and DNA re-association values that have delineated biologically important species. This review elaborates the present classification and taxonomy of rhizobia, vis-a-vis development of technical advancements, parameters and controversies associated with it, and describe the updated information on evolutionary lineages of rhizobia.

Characterization of symbiotic and nitrogen fixing bacteria

Symbiotic nitrogen fixing bacteria comprise of diverse species associated with the root nodules of leguminous plants. Using an appropriate taxonomic method to confirm the identity of superior and elite strains to fix nitrogen in legume crops can improve sustainable global food and nutrition security. The current review describes taxonomic methods preferred and commonly used to characterize symbiotic bacteria in the rhizosphere. Peer reviewed, published and unpublished articles on techniques used for detection, classification and identification of symbiotic bacteria were evaluated by exploring their advantages and limitations. The findings showed that phenotypic and cultural techniques are still affordable and remain the primary basis of species classification despite their challenges. Development of new, robust and informative taxonomic techniques has really improved characterization and identification of symbiotic bacteria and discovery of novel and new species that are effective in biological nitrogen fixation (BNF) in diverse conditions and environments.

Combining GWAS and population genomic analyses to characterize coevolution in a legume-rhizobia symbiosis

The mutualism between legumes and rhizobia is clearly the product of past coevolution. However, the nature of ongoing evolution between these partners is less clear. To characterize the nature of recent coevolution between legumes and rhizobia, we used population genomic analysis to characterize selection on functionally annotated symbiosis genes as well as on symbiosis gene candidates identified through a two-species association analysis. For the association analysis, we inoculated each of 202 accessions of the legume host Medicago truncatula with a community of 88 Sinorhizobia (Ensifer) meliloti strains. Multistrain inoculation, which better reflects the ecological reality of rhizobial selection in nature than single-strain inoculation, allows strains to compete for nodulation opportunities and host resources and for hosts to preferentially form nodules and provide resources to some strains. We found extensive host by symbiont, that is, genotype-by-genotype, effects on rhizobial fitness and some annotated rhizobial genes bear signatures of recent positive selection. However, neither genes responsible for this variation nor annotated host symbiosis genes are enriched for signatures of either positive or balancing selection. This result suggests that stabilizing selection dominates selection acting on symbiotic traits and that variation in these traits is under mutation-selection balance. Consistent with the lack of positive selection acting on host genes, we found that among-host variation in growth was similar whether plants were grown with rhizobia or N-fertilizer, suggesting that the symbiosis may not be a major driver of variation in plant growth in multistrain contexts.

Varietas Delectat: Exploring Natural Variations in Nitrogen-Fixing Symbiosis Research

The nitrogen-fixing symbiosis between leguminous plants and soil bacteria collectively called rhizobia plays an important role in the global nitrogen cycle and is an essential component of sustainable agriculture. Genetic determinants directing the development and functioning of the interaction have been identified with the help of a very limited number of model plants and bacterial strains. Most of the information obtained from the study of model systems could be validated on crop plants and their partners. The investigation of soybean cultivars and different rhizobia, however, has revealed the existence of ineffective interactions between otherwise effective partners that resemble gene-for-gene interactions described for pathogenic systems. Since then, incompatible interactions between natural isolates of model plants, called ecotypes, and different bacterial partner strains have been reported. Moreover, diverse phenotypes of both bacterial mutants on different host plants and plant mutants with different bacterial strains have been described. Identification of the genetic factors behind the phenotypic differences did already and will reveal novel functions of known genes/proteins, the role of certain proteins in some interactions, and the fine regulation of the steps during nodule development.

The Versatile Roles of Type III Secretion Systems in Rhizobia-Legume Symbioses

To suppress plant immunity and promote the intracellular infection required for fixing nitrogen for the benefit of their legume hosts, many rhizobia use type III secretion systems (T3SSs) that deliver effector proteins (T3Es) inside host cells. As reported for interactions between pathogens and host plants, the immune system of legume hosts and the cocktail of T3Es secreted by rhizobia determine the symbiotic outcome. If they remain undetected, T3Es may reduce plant immunity and thus promote infection of legumes by rhizobia. If one or more of the secreted T3Es are recognized by the cognate plant receptors, defense responses are triggered and rhizobial infection may abort. However, some rhizobial T3Es can also circumvent the need for nodulation (Nod) factors to trigger nodule formation. Here we review the multifaceted roles played by rhizobial T3Es during symbiotic interactions with legumes. Expected final online publication date for the Annual Review of Microbiology, Volume 76 is September 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

An Overview of the Mechanisms Against "Candidatus Liberibacter asiaticus": Virulence Targets, Citrus Defenses, and Microbiome

Citrus Huanglongbing (HLB) or citrus greening, is the most destructive disease for citrus worldwide. It is caused by the psyllid-transmitted, phloem-limited bacteria "Candidatus Liberibacter asiaticus" (CLas). To date, there are still no effective practical strategies for curing citrus HLB. Understanding the mechanisms against CLas can contribute to the development of effective approaches for combatting HLB. However, the unculturable nature of CLas has hindered elucidating mechanisms against CLas. In this review, we summarize the main aspects that contribute to the understanding about the mechanisms against CLas, including (1) CLas virulence targets, focusing on inhibition of virulence genes; (2) activation of citrus host defense genes and metabolites of HLB-tolerant citrus triggered by CLas, and by agents; and (3) we also review the role of citrus microbiome in combatting CLas. Finally, we discuss novel strategies to continue studying mechanisms against CLas and the relationship of above aspects.

Comparative Genomic Analyses of the Genus Nesterenkonia Unravels the Genomic Adaptation to Polar Extreme Environments

The members of the Nesterenkonia genus have been isolated from various habitats, like saline soil, salt lake, sponge-associated and the human gut, some of which are even located in polar areas. To identify their stress resistance mechanisms and draw a genomic profile across this genus, we isolated four Nesterenkonia strains from the lakes in the Tibetan Plateau, referred to as the third pole, and compared them with all other 30 high-quality Nesterenkonia genomes that are deposited in NCBI. The Heaps’ law model estimated that the pan-genome of this genus is open and the number of core, shell, cloud, and singleton genes were 993 (6.61%), 2782 (18.52%), 4117 (27.40%), and 7132 (47.47%), respectively. Phylogenomic and ANI/AAI analysis indicated that all genomes can be divided into three main clades, named NES-1, NES-2, and NES-3. The strains isolated from lakes in the Tibetan Plateau were clustered with four strains from different sources in the Antarctic and formed a subclade within NES-2, described as NES-AT. Genome features of this subclade, including GC (guanine + cytosine) content, tRNA number, carbon/nitrogen atoms per residue side chain (C/N-ARSC), and amino acid composition, in NES-AT individuals were significantly different from other strains, indicating genomic adaptation to cold, nutrient-limited, osmotic, and ultraviolet conditions in polar areas. Functional analysis revealed the enrichment of specific genes involved in bacteriorhodopsin synthesis, biofilm formation, and more diverse nutrient substance metabolism genes in the NES-AT clade, suggesting potential adaptation strategies for energy metabolism in polar environments. This study provides a comprehensive profile of the genomic features of the Nesterenkonia genus and reveals the possible mechanism for the survival of Nesterenkonia isolates in polar areas.

Specialization in a Nitrogen-Fixing Symbiosis: Proteome Differences Between Sinorhizobium medicae Bacteria and Bacteroids

Using tandem mass spectrometry (MS/MS), we analyzed the proteome of Sinorhizobium medicae WSM419 growing as free-living cells and in symbiosis with Medicago truncatula. In all, 3,215 proteins were identified, over half of the open reading frames predicted from the genomic sequence. The abundance of 1,361 proteins displayed strong lifestyle bias. In total, 1,131 proteins had similar levels in bacteroids and free-living cells, and the low levels of 723 proteins prevented statistically significant assignments. Nitrogenase subunits comprised approximately 12% of quantified bacteroid proteins. Other major bacteroid proteins included symbiosis-specific cytochromes and FixABCX, which transfer electrons to nitrogenase. Bacteroids had normal levels of proteins involved in amino acid biosynthesis, glycolysis or gluconeogenesis, and the pentose phosphate pathway; however, several amino acid degradation pathways were repressed. This suggests that bacteroids maintain a relatively independent anabolic metabolism. Tricarboxylic acid cycle proteins were highly expressed in bacteroids and no other catabolic pathway emerged as an obvious candidate to supply energy and reductant to nitrogen fixation. Bacterial stress response proteins were induced in bacteroids. Many WSM419 proteins that are not encoded in S. meliloti Rm1021 were detected, and understanding the functions of these proteins might clarify why S. medicae WSM419 forms a more effective symbiosis with M. truncatula than S. meliloti Rm1021.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Comparative Genomics across Three Ensifer Species Using a New Complete Genome Sequence of the Medicago Symbiont Sinorhizobium (Ensifer) meliloti WSM1022

Here, we report an improved and complete genome sequence of Sinorhizobium (Ensifer) meliloti strain WSM1022, a microsymbiont of Medicago species, revealing its tripartite structure. This improved genome sequence was generated combining Illumina and Oxford nanopore sequencing technologies to better understand the symbiotic properties of the bacterium. The 6.75 Mb WSM1022 genome consists of three scaffolds, corresponding to a chromosome (3.70 Mb) and the pSymA (1.38 Mb) and pSymB (1.66 Mb) megaplasmids. The assembly has an average GC content of 62.2% and a mean coverage of 77X. Genome annotation of WSM1022 predicted 6058 protein coding sequences (CDSs), 202 pseudogenes, 9 rRNAs (3 each of 5S, 16S, and 23S), 55 tRNAs, and 4 ncRNAs. We compared the genome of WSM1022 to two other rhizobial strains, closely related Sinorhizobium (Ensifer) meliloti Sm1021 and Sinorhizobium (Ensifer) medicae WSM419. Both WSM1022 and WSM419 species are high-efficiency rhizobial strains when in symbiosis with Medicago truncatula, whereas Sm1021 is ineffective. Our findings report significant genomic differences across the three strains with some similarities between the meliloti strains and some others between the high efficiency strains WSM1022 and WSM419. The addition of this high-quality rhizobial genome sequence in conjunction with comparative analyses will help to unravel the features that make a rhizobial symbiont highly efficient for nitrogen fixation.

Isolation and Characterization of High-Efficiency Rhizobia From Western Kenya Nodulating With Common Bean

Common bean is one of the primary protein sources in third-world countries. They form nodules with nitrogen-fixing rhizobia, which have to be adapted to the local soils. Commercial rhizobial strains such as Rhizobium tropici CIAT899 are often used in agriculture. However, this strain failed to significantly increase the common bean yield in many places, including Kenya, due to the local soils' low pH. We isolated two indigenous rhizobial strains from the nodules of common bean from two fields in Western Kenya that have never been exposed to commercial inocula. We then determined their ability to fix nitrogen in common beans, solubilize phosphorus, and produce indole acetic acid. In greenhouse experiments, common bean plants inoculated with two isolates, B3 and S2 in sterile vermiculite, performed better than those inoculated with CIAT899 or plants grown with nitrogen fertilizer alone. In contrast to CIAT899, both isolates grew in the media with pH 4.8. Furthermore, isolate B3 had higher phosphate solubilization ability and produced more indole acetic acid than the other two rhizobia. Genome analyses revealed that B3 and S2 are different strains of Rhizobium phaseoli. We recommend fieldwork studies in Kenyan soils to test the efficacy of the two isolates in the natural environment in an effort to produce inoculants specific for these soils.

Genomic Diversity of Pigeon Pea (Cajanus cajan L. Millsp.) Endosymbionts in India and Selection of Potential Strains for Use as Agricultural Inoculants

Pigeon pea (Cajanus cajan L. Millsp. ) is a legume crop resilient to climate change due to its tolerance to drought. It is grown by millions of resource-poor farmers in semiarid and tropical subregions of Asia and Africa and is a major contributor to their nutritional food security. Pigeon pea is the sixth most important legume in the world, with India contributing more than 70% of the total production and harbouring a wide variety of cultivars. Nevertheless, the low yield of pigeon pea grown under dry land conditions and its yield instability need to be improved. This may be done by enhancing crop nodulation and, hence, biological nitrogen fixation (BNF) by supplying effective symbiotic rhizobia through the application of "elite" inoculants. Therefore, the main aim in this study was the isolation and genomic analysis of effective rhizobial strains potentially adapted to drought conditions. Accordingly, pigeon pea endosymbionts were isolated from different soil types in Southern, Central, and Northern India. After functional characterisation of the isolated strains in terms of their ability to nodulate and promote the growth of pigeon pea, 19 were selected for full genome sequencing, along with eight commercial inoculant strains obtained from the ICRISAT culture collection. The phylogenomic analysis [Average nucleotide identity MUMmer (ANIm)] revealed that the pigeon pea endosymbionts were members of the genera Bradyrhizobium and Ensifer. Based on nodC phylogeny and nod cluster synteny, Bradyrhizobium yuanmingense was revealed as the most common endosymbiont, harbouring nod genes similar to those of Bradyrhizobium cajani and Bradyrhizobium zhanjiangense. This symbiont type (e.g., strain BRP05 from Madhya Pradesh) also outperformed all other strains tested on pigeon pea, with the notable exception of an Ensifer alkalisoli strain from North India (NBAIM29). The results provide the basis for the development of pigeon pea inoculants to increase the yield of this legume through the use of effective nitrogen-fixing rhizobia, tailored for the different agroclimatic regions of India.

Competitiveness for Nodule Colonization in Sinorhizobium meliloti: Combined In Vitro-Tagged Strain Competition and Genome-Wide Association Analysis

Associations between leguminous plants and symbiotic nitrogen-fixing rhizobia are a classic example of mutualism between a eukaryotic host and a specific group of prokaryotic microbes. Although this symbiosis is in part species specific, different rhizobial strains may colonize the same nodule. Some rhizobial strains are commonly known as better competitors than others, but detailed analyses that aim to predict rhizobial competitive abilities based on genomes are still scarce. Here, we performed a bacterial genome-wide association (GWAS) analysis to define the genomic determinants related to the competitive capabilities in the model rhizobial species Sinorhizobium meliloti. For this, 13 tester strains were green fluorescent protein (GFP) tagged and assayed versus 3 red fluorescent protein (RFP)-tagged reference competitor strains (Rm1021, AK83, and BL225C) in a Medicago sativa nodule occupancy test. Competition data and strain genomic sequences were employed to build a model for GWAS based on k-mers. Among the k-mers with the highest scores, 51 k-mers mapped on the genomes of four strains showing the highest competition phenotypes (>60% single strain nodule occupancy; GR4, KH35c, KH46, and SM11) versus BL225C. These k-mers were mainly located on the symbiosis-related megaplasmid pSymA, specifically on genes coding for transporters, proteins involved in the biosynthesis of cofactors, and proteins related to metabolism (e.g., fatty acids). The same analysis was performed considering the sum of single and mixed nodules obtained in the competition assays versus BL225C, retrieving k-mers mapped on the genes previously found and on vir genes. Therefore, the competition abilities seem to be linked to multiple genetic determinants and comprise several cellular components. IMPORTANCE Decoding the competitive pattern that occurs in the rhizosphere is challenging in the study of bacterial social interaction strategies. To date, the single-gene approach has mainly been used to uncover the bases of nodulation, but there is still a knowledge gap regarding the main features that a priori characterize rhizobial strains able to outcompete indigenous rhizobia. Therefore, tracking down which traits make different rhizobial strains able to win the competition for plant infection over other indigenous rhizobia will improve the strain selection process and, consequently, plant yield in sustainable agricultural production systems. We proved that a k-mer-based GWAS approach can efficiently identify the competition determinants of a panel of strains previously analyzed for their plant tissue occupancy using double fluorescent labeling. The reported strategy will be useful for detailed studies on the genomic aspects of the evolution of bacterial symbiosis and for an extensive evaluation of rhizobial inoculants.

Competitiveness for Nodule Colonization in Sinorhizobium meliloti: Combined In Vitro-Tagged Strain Competition and Genome-Wide Association Analysis

Associations between leguminous plants and symbiotic nitrogen-fixing rhizobia are a classic example of mutualism between a eukaryotic host and a specific group of prokaryotic microbes. Although this symbiosis is in part species specific, different rhizobial strains may colonize the same nodule. Some rhizobial strains are commonly known as better competitors than others, but detailed analyses that aim to predict rhizobial competitive abilities based on genomes are still scarce. Here, we performed a bacterial genome-wide association (GWAS) analysis to define the genomic determinants related to the competitive capabilities in the model rhizobial species Sinorhizobium meliloti. For this, 13 tester strains were green fluorescent protein (GFP) tagged and assayed versus 3 red fluorescent protein (RFP)-tagged reference competitor strains (Rm1021, AK83, and BL225C) in a Medicago sativa nodule occupancy test. Competition data and strain genomic sequences were employed to build a model for GWAS based on k-mers. Among the k-mers with the highest scores, 51 k-mers mapped on the genomes of four strains showing the highest competition phenotypes (>60% single strain nodule occupancy; GR4, KH35c, KH46, and SM11) versus BL225C. These k-mers were mainly located on the symbiosis-related megaplasmid pSymA, specifically on genes coding for transporters, proteins involved in the biosynthesis of cofactors, and proteins related to metabolism (e.g., fatty acids). The same analysis was performed considering the sum of single and mixed nodules obtained in the competition assays versus BL225C, retrieving k-mers mapped on the genes previously found and on vir genes. Therefore, the competition abilities seem to be linked to multiple genetic determinants and comprise several cellular components.

IMPORTANCE Decoding the competitive pattern that occurs in the rhizosphere is challenging in the study of bacterial social interaction strategies. To date, the single-gene approach has mainly been used to uncover the bases of nodulation, but there is still a knowledge gap regarding the main features that a priori characterize rhizobial strains able to outcompete indigenous rhizobia. Therefore, tracking down which traits make different rhizobial strains able to win the competition for plant infection over other indigenous rhizobia will improve the strain selection process and, consequently, plant yield in sustainable agricultural production systems. We proved that a k-mer-based GWAS approach can efficiently identify the competition determinants of a panel of strains previously analyzed for their plant tissue occupancy using double fluorescent labeling. The reported strategy will be useful for detailed studies on the genomic aspects of the evolution of bacterial symbiosis and for an extensive evaluation of rhizobial inoculants.

Spatial Virome Analysis of Zanthoxylum armatum Trees Affected With the Flower Yellowing Disease

Zanthoxylum armatum is an important woody crop with multiple applications in pharmaceutics, cosmetics, and food industries. With continuous increases in the plantation area, integrated pest management is required for scale production when diseases caused by biotic factors such as pests and pathogens have become new problems, one of which is the infectious flower yellowing disease (FYD). Here, isolates of a new illarvirus (3) and a new nepovirus-associated subviral satellite RNA (12) were identified in Z. armatum, in addition to 38 new isolates of four previously reported RNA viruses. Sequence variation can be observed in viral/subviral quasispecies and among predominant isolates from the same or different samples and geographic origins. Intriguingly, RNA sequencing of different diseased trees invariably showed an extraordinary pattern of particularly high reads accumulation of the green Sichuan pepper-nepovirus (GSPNeV) and the satellite RNA in symptomatic tissues. In addition, we also examined small RNAs of the satellite RNA, which show similar patterns to those of coinfecting viruses. This study provides further evidence to support association of the FYD with viral/subviral infections and deepens our understanding of the diversity and molecular characteristics of the viruses and satellite, as well as their interactions with the host.

Minimal gene set from Sinorhizobium (Ensifer) meliloti pSymA required for efficient symbiosis with Medicago

Reduction of N₂ gas to ammonia in legume root nodules is a key component of sustainable agricultural systems. Root nodules are the result of a symbiosis between leguminous plants and bacteria called rhizobia. Both symbiotic partners play active roles in establishing successful symbiosis and nitrogen fixation: while root nodule development is mostly controlled by the plant, the rhizobia induce nodule formation, invade, and perform N₂ fixation once inside the plant cells. Many bacterial genes involved in the rhizobia-legume symbiosis are known, and there is much interest in engineering the symbiosis to include major nonlegume crops such as corn, wheat, and rice. We sought to identify and combine a minimal bacterial gene complement necessary and sufficient for symbiosis. We analyzed a model rhizobium, Sinorhizobium (Ensifer) meliloti, using a background strain in which the 1.35-Mb symbiotic megaplasmid pSymA was removed. Three regions representing 162 kb of pSymA were sufficient to recover a complete N₂-fixing symbiosis with alfalfa, and a targeted assembly of this gene complement achieved high levels of symbiotic N₂ fixation. The resulting gene set contained just 58 of 1,290 pSymA protein-coding genes. To generate a platform for future synthetic manipulation, the minimal symbiotic genes were reorganized into three discrete nod, nif, and fix modules. These constructs will facilitate directed studies toward expanding the symbiosis to other plant partners. They also enable forward-type approaches to identifying genetic components that may not be essential for symbiosis, but which modulate the rhizobium's competitiveness for nodulation and the effectiveness of particular rhizobia-plant symbioses.

Biodiversity-based development and evolution: the emerging research systems in model and non-model organisms

Evolutionary developmental biology, or Evo-Devo for short, has become an established field that, broadly speaking, seeks to understand how changes in development drive major transitions and innovation in organismal evolution. It does so via integrating the principles and methods of many subdisciplines of biology. Although we have gained unprecedented knowledge from the studies on model organisms in the past decades, many fundamental and crucially essential processes remain a mystery. Considering the tremendous biodiversity of our planet, the current model organisms seem insufficient for us to understand the evolutionary and physiological processes of life and its adaptation to exterior environments. The currently increasing genomic data and the recently available gene-editing tools make it possible to extend our studies to non-model organisms. In this review, we review the recent work on the regulatory signaling of developmental and regeneration processes, environmental adaptation, and evolutionary mechanisms using both the existing model animals such as zebrafish and Drosophila, and the emerging nonstandard model organisms including amphioxus, ascidian, ciliates, single-celled phytoplankton, and marine nematode. In addition, the challenging questions and new directions in these systems are outlined as well.

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.

Integrating pan-genome with metagenome for microbial community profiling

Advances in sequencing technology have led to the increased availability of genomes and metagenomes, which has greatly facilitated microbial pan-genome and metagenome analysis in the community. In line with this trend, studies on microbial genomes and phenotypes have gradually shifted from individuals to environmental communities.

Pan-genomics and metagenomics are powerful strategies for in-depth profiling study of microbial communities. Pan-genomics focuses on genetic diversity, dynamics, and phylogeny at the multi-genome level, while metagenomics profiles the distribution and function of culture-free microbial communities in special environments. Combining pan-genome and metagenome analysis can reveal the microbial complicated connections from an individual complete genome to a mixture of genomes, thereby extending the catalog of traditional individual genomic profile to community microbial profile. Therefore, the combination of pan-genome and metagenome approaches has become a promising method to track the sources of various microbes and decipher the population-level evolution and ecosystem functions.

This review summarized the pan-genome and metagenome approaches, the combined strategies of pan-genome and metagenome, and applications of these combined strategies in studies of microbial dynamics, evolution, and function in communities. We discussed emerging strategies for the study of microbial communities that integrate information in both pan-genome and metagenome. We emphasized studies in which the integrating pan-genome with metagenome approach improved the understanding of models of microbial community profiles, both structural and functional. Finally, we illustrated future perspectives of microbial community profile: more advanced analytical techniques, including big-data based artificial intelligence, will lead to an even better understanding of the patterns of microbial communities.

Comparative Genomics, from the Annotated Genome to Valuable Biological Information: A Case Study

High availability of fast, cheap, and high-throughput next generation sequencing techniques resulted in acquisition of numerous de novo sequenced and assembled bacterial genomes. It rapidly became clear that digging out useful biological information from such a huge amount of data presents a considerable challenge. In this chapter we share our experience with utilization of several handy open source comparative genomic tools. All of them were applied in the studies focused on revealing inter- and intraspecies variation in pectinolytic plant pathogenic bacteria classified to Dickeya solani and Pectobacterium parmentieri. As the described software performed well on the species within the Pectobacteriaceae family, it presumably may be readily utilized on some closely related taxa from the Enterobacteriaceae family. First of all, implementation of various annotation software is discussed and compared. Then, tools computing whole genome comparisons including generation of circular juxtapositions of multiple sequences, revealing the order of synteny blocks or calculation of ANI or Tetra values are presented. Besides, web servers intended either for functional annotation of the genes of interest or for detection of genomic islands, plasmids, prophages, CRISPR/Cas are described. Last but not least, utilization of the software designed for pangenome studies and the further downstream analyses is explained. The presented work not only summarizes broad possibilities assured by the comparative genomic approach but also provides a user-friendly guide that might be easily followed by nonbioinformaticians interested in undertaking similar studies.

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.

Massive rhizobial genomic variation associated with partner quality in Lotus-Mesorhizobium symbiosis

Variation in partner quality is commonly observed in diverse cooperative relationships, despite the theoretical prediction that selection favoring high-quality partners should eliminate such variation. Here, we investigated how genetic variation in partner quality could be maintained in the nitrogen-fixing mutualism between Lotus japonicus and Mesorhizobium bacteria. We reconstructed de novo assembled full-genome sequences from nine rhizobial symbionts, finding massive variation in the core genome and the similar symbiotic islands, indicating recent horizontal gene transfer (HGT) of the symbiosis islands into diverse Mesorhizobium lineages. A cross-inoculation experiment using 9 sequenced rhizobial symbionts and 15 L. japonicus accessions revealed extensive quality variation represented by plant growth phenotypes, including genotype-by-genotype interactions. Variation in quality was not associated with the presence/absence variation in known symbiosis-related genes in the symbiosis island; rather, it showed significant correlation with the core genome variation. Given the recurrent HGT of the symbiosis islands into diverse Mesorhizobium strains, local Mesorhizobium communities could serve as a major source of variation for core genomes, which might prevent variation in partner quality from fixing, even in the presence of selection favoring high-quality partners. These findings highlight the novel role of HGT of symbiosis islands in maintaining partner quality variation in the legume-rhizobia symbiosis.

Symbiotic and Nonsymbiotic Members of the Genus Ensifer (syn. Sinorhizobium) Are Separated into Two Clades Based on Comparative Genomics and High-Throughput Phenotyping

Rhizobium–legume symbioses serve as paradigmatic examples for the study of mutualism evolution. The genus Ensifer (syn. Sinorhizobium) contains diverse plant-associated bacteria, a subset of which can fix nitrogen in symbiosis with legumes. To gain insights into the evolution of symbiotic nitrogen fixation (SNF), and interkingdom mutualisms more generally, we performed extensive phenotypic, genomic, and phylogenetic analyses of the genus Ensifer. The data suggest that SNF likely emerged several times within the genus Ensifer through independent horizontal gene transfer events. Yet, the majority (105 of 106) of the Ensifer strains with the nodABC and nifHDK nodulation and nitrogen fixation genes were found within a single, monophyletic clade. Comparative genomics highlighted several differences between the “symbiotic” and “nonsymbiotic” clades, including divergences in their pangenome content. Additionally, strains of the symbiotic clade carried 325 fewer genes, on average, and appeared to have fewer rRNA operons than strains of the nonsymbiotic clade. Initial characterization of a subset of ten Ensifer strains identified several putative phenotypic differences between the clades. Tested strains of the nonsymbiotic clade could catabolize 25% more carbon sources, on average, than strains of the symbiotic clade, and they were better able to grow in LB medium and tolerate alkaline conditions. On the other hand, the tested strains of the symbiotic clade were better able to tolerate heat stress and acidic conditions. We suggest that these data support the division of the genus Ensifer into two main subgroups, as well as the hypothesis that pre-existing genetic features are required to facilitate the evolution of SNF in bacteria.

Rhizobial Microsymbionts of Kamchatka Oxytropis Species Possess Genes of the Type III and VI Secretion Systems, Which Can Affect the Development of Symbiosis

A collection of rhizobial strains isolated from root nodules of the narrowly endemic legume species Oxytropis erecta, O. anadyrensis, O. kamtschatica, and O. pumilio originating from the Kamchatka Peninsula (Russian Federation) was obtained. Analysis of the 16S ribosomal RNA gene sequence showed a significant diversity of isolates belonging to families Rhizobiaceae (genus Rhizobium), Phyllobacteriaceae (genera Mesorhizobium, Phyllobacterium), and Bradyrhizobiaceae (genera Bosea, Tardiphaga). A plant nodulation assay showed that only strains belonging to genus Mesorhizobium could form nitrogen-fixing nodules on Oxytropis plants. The strains M. loti 582 and M. huakuii 583, in addition to symbiotic clusters, possessed genes of the type III and type VI secretion systems (T3SS and T6SS, respectively), which can influence the host specificity of strains. These strains formed nodules of two types (elongated and rounded) on O. kamtschatica roots. We suggest this phenomenon may result from Nod factor-dependent and -independent nodulation strategies. The obtained strains are of interest for further study of the T3SS and T6SS gene function and their role in the development of rhizobium-legume symbiosis. The prospects of using rhizobia having both gene systems related to symbiotic and nonsymbiotic nodulation strategies to enhance the efficiency of plant-microbe interactions by expanding the host specificity and increasing nodulation efficiency are discussed.

Polyploid plants obtain greater fitness benefits from a nutrient acquisition mutualism

Polyploidy is a key driver of ecological and evolutionary processes in plants, yet little is known about its effects on biotic interactions. This gap in knowledge is especially profound for nutrient acquisition mutualisms, despite the fact that they regulate global nutrient cycles and structure ecosystems. Generalism in mutualistic interactions depends on the range of potential partners (niche breadth), the benefits obtained and ability to maintain benefits across a variety of partners (fitness plasticity). Here, we determine how each of these is influenced by polyploidy in the legume-rhizobium mutualism. We inoculated a broad geographic sample of natural diploid and autotetraploid alfalfa (Medicago sativa) lineages with a diverse panel of Sinorhizobium bacterial symbionts. To analyze the extent and mechanism of generalism, we measured host growth benefits and functional traits. Autotetraploid plants obtained greater fitness enhancement from mutualistic interactions and were better able to maintain this across diverse rhizobial partners (i.e. low plasticity in fitness) relative to diploids. These benefits were not attributed to increases in niche breadth, but instead reflect increased rewards from investment in the mutualism. Polyploid plants displayed greater generalization in bacterial mutualisms relative to diploids, illustrating another axis of advantage for polyploids over diploids.

Complete genome sequencing and comparative genomic analyses of Bacillus sp. S3, a novel hyper Sb(III)-oxidizing bacterium

BACKGROUND

Antimonite [Sb(III)]-oxidizing bacterium has great potential in the environmental bioremediation of Sb-polluted sites. Bacillus sp. S3 that was previously isolated from antimony-contaminated soil displayed high Sb(III) resistance and Sb(III) oxidation efficiency. However, the genomic information and evolutionary feature of Bacillus sp. S3 are very scarce.

RESULTS

Here, we identified a 5,436,472 bp chromosome with 40.30% GC content and a 241,339 bp plasmid with 36.74% GC content in the complete genome of Bacillus sp. S3. Genomic annotation showed that Bacillus sp. S3 contained a key aioB gene potentially encoding As (III)/Sb(III) oxidase, which was not shared with other Bacillus strains. Furthermore, a wide variety of genes associated with Sb(III) and other heavy metal (loid) s were also ascertained in Bacillus sp. S3, reflecting its adaptive advantage for growth in the harsh eco-environment. Based on the analysis of phylogenetic relationship and the average nucleotide identities (ANI), Bacillus sp. S3 was proved to a novel species within the Bacillus genus. The majority of mobile genetic elements (MGEs) mainly distributed on chromosomes within the Bacillus genus. Pan-genome analysis showed that the 45 genomes contained 554 core genes and many unique genes were dissected in analyzed genomes. Whole genomic alignment showed that Bacillus genus underwent frequently large-scale evolutionary events. In addition, the origin and evolution analysis of Sb(III)-resistance genes revealed the evolutionary relationships and horizontal gene transfer (HGT) events among the Bacillus genus. The assessment of functionality of heavy metal (loid) s resistance genes emphasized its indispensable role in the harsh eco-environment of Bacillus genus. Real-time quantitative PCR (RT-qPCR) analysis indicated that Sb(III)-related genes were all induced under the Sb(III) stress, while arsC gene was down-regulated.

CONCLUSIONS

The results in this study shed light on the molecular mechanisms of Bacillus sp. S3 coping with Sb(III), extended our understanding on the evolutionary relationships between Bacillus sp. S3 and other closely related species, and further enriched the Sb(III) resistance genetic data sources.

Diversity of Bradyrhizobium in Non-Leguminous Sorghum Plants: B. ottawaense Isolates Unique in Genes for N2O Reductase and Lack of the Type VI Secretion System

Diverse members of Bradyrhizobium diazoefficiens, B. japonicum, and B. ottawaense were isolated from the roots of field-grown sorghum plants in Fukushima, and classified into “Rhizobia” with nodulated soybeans, “Free-living diazotrophs”, and “Non-diazotrophs” by nitrogen fixation and nodulation assays. Genome analyses revealed that B. ottawaense members possessed genes for N₂O reduction, but lacked those for the Type VI secretion system (T6SS). T6SS is a new bacterial weapon against microbial competitors. Since T6SS-possessing B. diazoefficiens and B. japonicum have mainly been isolated from soybean nodules in Japan, T6SS-lacking B. ottawaense members may be a cryptic lineage of soybean bradyrhizobia in Japan.

Gene Expression in Nitrogen-Fixing Symbiotic Nodule Cells in Medicago truncatula and Other Nodulating Plants

Root nodules formed by plants of the nitrogen-fixing clade (NFC) are symbiotic organs that function in the maintenance and metabolic integration of large populations of nitrogen-fixing bacteria. These organs feature unique characteristics and processes, including their tissue organization, the presence of specific infection structures called infection threads, endocytotic uptake of bacteria, symbiotic cells carrying thousands of intracellular bacteria without signs of immune responses, and the integration of symbiont and host metabolism. The early stages of nodulation are governed by a few well-defined functions, which together constitute the common symbiosis-signaling pathway (CSSP). The CSSP activates a set of transcription factors (TFs) that orchestrate nodule organogenesis and infection. The later stages of nodule development require the activation of hundreds to thousands of genes, mostly expressed in symbiotic cells. Many of these genes are only active in symbiotic cells, reflecting the unique nature of nodules as plant structures. Although how the nodule-specific transcriptome is activated and connected to early CSSP-signaling is poorly understood, candidate TFs have been identified using transcriptomic approaches, and the importance of epigenetic and chromatin-based regulation has been demonstrated. We discuss how gene regulation analyses have advanced our understanding of nodule organogenesis, the functioning of symbiotic cells, and the evolution of symbiosis in the NFC.

An Alkane Sulfonate Monooxygenase Is Required for Symbiotic Nitrogen Fixation by Bradyrhizobium diazoefficiens (syn. Bradyrhizobium japonicum) USDA110T

Sulfur (S)-containing molecules play an important role in symbiotic nitrogen fixation and are critical components of nitrogenase and other iron-S proteins. S deficiency inhibits symbiotic nitrogen fixation by rhizobia. However, despite its importance, little is known about the sources of S that rhizobia utilize during symbiosis. We previously showed that Bradyrhizobium diazoefficiens USDA110T can assimilate both inorganic and organic S and that genes involved in organic S utilization are expressed during symbiosis. Here, we show that a B. diazoefficiens USDA110T mutant with a sulfonate monooxygenase (ssuD) insertion is defective in nitrogen fixation. Microscopy analyses revealed that the ΔssuD mutant was defective in root hair infection and that ΔssuD mutant bacteroids showed degradation compared to the wild-type strain. Moreover, the ΔssuD mutant was significantly more sensitive to hydrogen peroxide-mediated oxidative stress than the wild-type strain. Taken together, these results show that the ability of rhizobia to utilize organic S plays an important role in symbiotic nitrogen fixation. Since nodules have been reported to be an important source of reduced S used during symbiosis and nitrogen fixation, further research will be needed to determine the mechanisms involved in the regulation of S assimilation by rhizobia.IMPORTANCE Rhizobia form symbiotic associations with legumes that lead to the formation of nitrogen-fixing nodules. Sulfur-containing molecules play a crucial role in nitrogen fixation; thus, the rhizobia inside nodules require large amounts of sulfur. Rhizobia can assimilate both inorganic (sulfate) and organic (sulfonates) sources of sulfur. However, very little is known about rhizobial sulfur metabolism during symbiosis. In this report, we show that sulfonate utilization by Bradyrhizobium diazoefficiens is important for symbiotic nitrogen fixation in both soybean and cowpea. The symbiotic defect is probably due to increased sensitivity to oxidative stress from sulfur deficiency in the mutant strain defective for sulfonate utilization. The results of this study can be extended to other rhizobium-legume symbioses, as sulfonate utilization genes are widespread in these bacteria.

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.

The genome of Ensifer alkalisoli YIC4027 provides insights for host specificity and environmental adaptations

Background

Ensifer alkalisoli YIC4027, a recently characterized nitrogen-fixing bacterium of the genus Ensifer, has been isolated from root nodules of the host plant Sesbania cannabina. This plant is widely used as green manure and for soil remediation. E. alkalisoli YIC4027 can grow in saline-alkaline soils and is a narrow-host-range strain that establishes a symbiotic relationship with S. cannabina. The complete genome of this strain was sequenced to better understand the genetic basis of host specificity and adaptation to saline-alkaline soils.

Results

E. alkalisoli YIC4027 was found to possess a 6.1-Mb genome consisting of three circular replicons: one chromosome (3.7 Mb), a chromid (1.9 Mb) and a plasmid (0.46 Mb). Genome comparisons showed that strain YIC4027 is phylogenetically related to broad-host-range Ensifer fredii strains. Synteny analysis revealed a strong collinearity between chromosomes of E. alkalisoli YIC4027 and those of the E. fredii NGR234 (3.9 Mb), HH103 (4.3 Mb) and USDA257 (6.48 Mb) strains. Notable differences were found for genes required for biosynthesis of nodulation factors and protein secretion systems, suggesting a role of these genes in host-specific nodulation. In addition, the genome analysis led to the identification of YIC4027 genes that are presumably related to adaptation to saline-alkaline soils, rhizosphere colonization and nodulation competitiveness. Analysis of chemotaxis cluster genes and nodulation tests with constructed che gene mutants indicated a role of chemotaxis and flagella-mediated motility in the symbiotic association between YIC4027 and S. cannabina.

Conclusions

This study provides a basis for a better understanding of host specific nodulation and of adaptation to a saline-alkaline rhizosphere. This information offers the perspective to prepare optimal E. alkalisoli inocula for agriculture use and soil remediation.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6004-7) contains supplementary material, which is available to authorized users.

Role of antimicrobial peptides in controlling symbiotic bacterial populations

Covering: up to 2018 Antimicrobial peptides (AMPs) have been known for well over three decades as crucial mediators of the innate immune response in animals and plants, where they are involved in the killing of infecting microbes. However, AMPs have now also been found to be produced by eukaryotic hosts during symbiotic interactions with bacteria. These symbiotic AMPs target the symbionts and therefore have a more subtle biological role: not eliminating the microbial symbiont population but rather keeping it in check. The arsenal of AMPs and the symbionts' adaptations to resist them are in a careful balance, which contributes to the establishment of the host-microbe homeostasis. Although in many cases the biological roles of symbiotic AMPs remain elusive, for a number of symbiotic interactions, precise functions have been assigned or proposed to the AMPs, which are discussed here. The microbiota living on epithelia in animals, from the most primitive ones to the mammals, are challenged by a cocktail of AMPs that determine the specific composition of the bacterial community as well as its spatial organization. In the symbiosis of legume plants with nitrogen-fixing rhizobium bacteria, the host deploys an extremely large panel of AMPs - called nodule-specific cysteine-rich (NCR) peptides - that drive the bacteria into a terminally differentiated state and manipulate the symbiont physiology to maximize the benefit for the host. The NCR peptides are used as tools to enslave the bacterial symbionts, limiting their reproduction but keeping them metabolically active for nitrogen fixation. In the nutritional symbiotic interactions of insects and protists that have vertically transmitted bacterial symbionts with reduced genomes, symbiotic AMPs could facilitate the integration of the endosymbiont and host metabolism by favouring the flow of metabolites across the symbiont membrane through membrane permeabilization.

Expression of the Arabidopsis thaliana immune receptor EFR in Medicago truncatula reduces infection by a root pathogenic bacterium, but not nitrogen-fixing rhizobial symbiosis

Interfamily transfer of plant pattern recognition receptors (PRRs) represents a promising biotechnological approach to engineer broad-spectrum, and potentially durable, disease resistance in crops. It is however unclear whether new recognition specificities to given pathogen-associated molecular patterns (PAMPs) affect the interaction of the recipient plant with beneficial microbes. To test this in a direct reductionist approach, we transferred the Brassicaceae-specific PRR ELONGATION FACTOR-THERMO UNSTABLE RECEPTOR (EFR), conferring recognition of the bacterial EF-Tu protein, from Arabidopsis thaliana to the legume Medicago truncatula. Constitutive EFR expression led to EFR accumulation and activation of immune responses upon treatment with the EF-Tu-derived elf18 peptide in leaves and roots. The interaction of M. truncatula with the bacterial symbiont Sinorhizobium meliloti is characterized by the formation of root nodules that fix atmospheric nitrogen. Although nodule numbers were slightly reduced at an early stage of the infection in EFR-Medicago when compared to control lines, nodulation was similar in all lines at later stages. Furthermore, nodule colonization by rhizobia, and nitrogen fixation were not compromised by EFR expression. Importantly, the M. truncatula lines expressing EFR were substantially more resistant to the root bacterial pathogen Ralstonia solanacearum. Our data suggest that the transfer of EFR to M. truncatula does not impede root nodule symbiosis, but has a positive impact on disease resistance against a bacterial pathogen. In addition, our results indicate that Rhizobium can either avoid PAMP recognition during the infection process, or is able to actively suppress immune signaling.

Flowering plant immune repertoires expand under mycorrhizal symbiosis

Immune perception in flowering plants is mediated by a repertoire of cytoplasmic and cell-surface receptors that detect invading microbes and their effects on cells. Here, we show that several large families of immune receptors exhibit size variations related to a plant's competence to host symbiotic root fungi (mycorrhiza). Plants that do not participate in mycorrhizal associations have significantly smaller immune repertoires, while the most promiscuous symbiotic hosts (ectomycorrhizal plant species) have significantly larger immune repertoires. By contrast, we find no significant increase in immune repertoire size among legumes competent to form a symbiosis with nitrogen-fixing bacteria (rhizobia). To explain these observations, we hypothesize that plant immune repertoire size expands with symbiote species diversity.

traG Gene Is Conserved across Mesorhizobium spp. Able to Nodulate the Same Host Plant and Expressed in Response to Root Exudates

Evidences for an involvement of the bacterial type IV secretion system (T4SS) in the symbiotic relationship between rhizobia and legumes have been pointed out by several recent studies. However, information regarding this secretion system in Mesorhizobium is still very scarce. The aim of the present study was to investigate the phylogeny and expression of the traG gene, which encodes a substrate receptor of the T4SS. In addition, the occurrence and genomic context of this and other T4SS genes, namely, genes from tra/trb and virB/virD4 complexes, were also analyzed in order to unveil the structural and functional organization of T4SS in mesorhizobia. The location of the T4SS genes in the symbiotic region of the analyzed rhizobial genomes, along with the traG phylogeny, suggests that T4SS genes could be horizontally transferred together with the symbiosis genes. Regarding the T4SS structural organization in Mesorhizobium, the virB/virD4 genes were absent in all chickpea (Cicer arietinum L.) microsymbionts and in the Lotus symbiont Mesorhizobium japonicum MAFF303099^T. Interestingly, the presence of genes belonging to another secretion system (T3SS) was restricted to these strains lacking the virB/virD4 genes. The traG gene expression was detected in M. mediterraneum Ca36^T and M. ciceri LMS-1 strains when exposed to chickpea root exudates and also in the early nodules formed by M. mediterraneum Ca36^T, but not in older nodules. This study contributes to a better understanding of the importance of T4SS in mutualistic symbiotic bacteria.

Comparative Genomic Analysis Reveals the Distribution, Organization, and Evolution of Metal Resistance Genes in the Genus Acidithiobacillus

Members of the genus Acidithiobacillus, which can adapt to extremely high concentrations of heavy metals, are universally found at acid mine drainage (AMD) sites. Here, we performed a comparative genomic analysis of 37 strains within the genus Acidithiobacillus to answer the untouched questions as to the mechanisms and the evolutionary history of metal resistance genes in Acidithiobacillus spp. The results showed that the evolutionary history of metal resistance genes in Acidithiobacillus spp. involved a combination of gene gains and losses, horizontal gene transfer (HGT), and gene duplication. Phylogenetic analyses revealed that metal resistance genes in Acidithiobacillus spp. were acquired by early HGT events from species that shared habitats with Acidithiobacillus spp., such as Acidihalobacter, Thiobacillus, Acidiferrobacter, and Thiomonas species. Multicopper oxidase genes involved in copper detoxification were lost in iron-oxidizing Acidithiobacillus ferridurans, Acidithiobacillus ferrivorans, and Acidithiobacillus ferrooxidans and were replaced by rusticyanin genes during evolution. In addition, widespread purifying selection and the predicted high expression levels emphasized the indispensable roles of metal resistance genes in the ability of Acidithiobacillus spp. to adapt to harsh environments. Altogether, the results suggested that Acidithiobacillus spp. recruited and consolidated additional novel functionalities during the adaption to challenging environments via HGT, gene duplication, and purifying selection. This study sheds light on the distribution, organization, functionality, and complex evolutionary history of metal resistance genes in Acidithiobacillus spp.IMPORTANCE Horizontal gene transfer (HGT), natural selection, and gene duplication are three main engines that drive the adaptive evolution of microbial genomes. Previous studies indicated that HGT was a main adaptive mechanism in acidophiles to cope with heavy-metal-rich environments. However, evidences of HGT in Acidithiobacillus species in response to challenging metal-rich environments and the mechanisms addressing how metal resistance genes originated and evolved in Acidithiobacillus are still lacking. The findings of this study revealed a fascinating phenomenon of putative cross-phylum HGT, suggesting that Acidithiobacillus spp. recruited and consolidated additional novel functionalities during the adaption to challenging environments via HGT, gene duplication, and purifying selection. Altogether, the insights gained in this study have improved our understanding of the metal resistance strategies of Acidithiobacillus spp.