Comparative analysis of integrative and conjugative mobile genetic elements in the genus Mesorhizobium

Beav: a bacterial genome and mobile element annotation pipeline

Comprehensive and accurate genome annotation is crucial for inferring the predicted functions of an organism. Numerous tools exist to annotate genes, gene clusters, mobile genetic elements, and other diverse features. However, these tools and pipelines can be difficult to install and run, be specialized for a particular element or feature, or lack annotations for larger elements that provide important genomic context. Integrating results across analyses is also important for understanding gene function. To address these challenges, we present the Beav annotation pipeline. Beav is a command-line tool that automates the annotation of bacterial genome sequences, mobile genetic elements, molecular systems and gene clusters, key regulatory features, and other elements. Beav uses existing tools in addition to custom models, scripts, and databases to annotate diverse elements, systems, and sequence features. Custom databases for plant-associated microbes are incorporated to improve annotation of key virulence and symbiosis genes in agriculturally important pathogens and mutualists. Beav includes an optional Agrobacterium-specific pipeline that identifies and classifies oncogenic plasmids and annotates plasmid-specific features. Following the completion of all analyses, annotations are consolidated to produce a single comprehensive output. Finally, Beav generates publication-quality genome and plasmid maps. Beav is on Bioconda and is available for download at https://github.com/weisberglab/beav.

IMPORTANCE

Annotation of genome features, such as the presence of genes and their predicted function, or larger loci encoding secretion systems or biosynthetic gene clusters, is necessary for understanding the functions encoded by an organism. Genomes can also host diverse mobile genetic elements, such as integrative and conjugative elements and/or phages, that are often not annotated by existing pipelines. These elements can horizontally mobilize genes encoding for virulence, antimicrobial resistance, or other adaptive functions and alter the phenotype of an organism. We developed a software pipeline, called Beav, that combines new and existing tools for the comprehensive annotation of these and other major features. Existing pipelines often misannotate loci important for virulence or mutualism in plant-associated bacteria. Beav includes custom databases and optional workflows for the improved annotation of plant-associated bacteria. Beav is designed to be easy to install and run, making comprehensive genome annotation broadly available to the research community.

Identifying functional multi-host shuttle plasmids to advance synthetic biology applications in Mesorhizobium and Bradyrhizobium

Ammonia availability has a crucial role in agriculture as it ensures healthy plant growth and increased crop yields. Since diazotrophs are the only organisms capable of reducing dinitrogen to ammonia, they have great ecological importance and potential to mitigate the environmental and economic costs of synthetic fertilizer use. Rhizobia are especially valuable being that they can engage in nitrogen-fixing symbiotic relationships with legumes, and they demonstrate great diversity and plasticity in genomic and phenotypic traits. However, few rhizobial species have sufficient genetic tractability for synthetic biology applications. This study established a basic genetic toolbox with antibiotic resistance markers, multi-host shuttle plasmids and a streamlined protocol for biparental conjugation with Mesorhizobium and Bradyrhizobium species. We identified two repABC origins of replication from Sinorhizobium meliloti (pSymB) and Rhizobium etli (p42d) that were stable across all three strains of interest. Furthermore, the NZP2235 genome was sequenced and phylogenetic analysis determined its reclassification to Mesorhizobium huakuii. These tools will enable the use of plasmid-based strategies for more advanced genetic engineering projects and ultimately contribute towards the development of more sustainable agriculture practices by means of novel nitrogen-fixing organelles, elite bioinoculants, or symbiotic association with nonlegumes.

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.

Integrative Conjugative Elements (ICEs) of the SXT/R391 family drive adaptation and evolution in γ-Proteobacteria

Integrative Conjugative Elements (ICEs) are mosaics containing functional modules allowing maintenance by site-specific integration and excision into and from the host genome and conjugative transfer to a specific host range. Many ICEs encode a range of adaptive functions that aid bacterial survival and evolution in a range of niches. ICEs from the SXT/R391 family are found in γ-Proteobacteria. Over 100 members have undergone epidemiological and molecular characterization allowing insight into their diversity and function. Comparative analysis of SXT/R391 elements from a wide geographic distribution has revealed conservation of key functions, and the accumulation and evolution of adaptive genes. This evolution is associated with gene acquisition in conserved hotspots and variable regions within the SXT/R391 ICEs catalysed via element-encoded recombinases. The elements can carry IS elements and transposons, and a mutagenic DNA polymerase, PolV, which are associated with their evolution. SXT/R391 ICEs isolated from different niches appear to have retained adaptive functions related to that specific niche; phage resistance determinants in ICEs carried by wastewater bacteria, antibiotic resistance determinants in clinical isolates and metal resistance determinants in bacteria recovered from polluted environments/ocean sediments. Many genes found in the element hotspots are undetermined and have few homologs in the nucleotide databases.

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.

Comparative genome analysis of Sesbania cannabina-nodulating Rhizobium spp. revealing the symbiotic and transferrable characteristics of symbiosis plasmids

Symbiotic nitrogen fixation between legumes and rhizobia makes a great contribution to the terrestrial ecosystem. The successful symbiosis between the partners mainly depends on the nod and nif genes in rhizobia, while the specific symbiosis is mainly determined by the structure of Nod factors and the corresponding secretion systems (type III secretion system; T3SS), etc. These symbiosis genes are usually located on symbiotic plasmids or a chromosomal symbiotic island, both could be transferred interspecies. In our previous studies, Sesbania cannabina-nodulating rhizobia across the world were classified into 16 species of four genera and all the strains, especially those of Rhizobium spp., harboured extraordinarily highly conserved symbiosis genes, suggesting that horizontal transfer of symbiosis genes might have happened among them. In order to learn the genomic basis of diversification of rhizobia under the selection of host specificity, we performed this study to compare the complete genome sequences of four Rhizobium strains associated with S. cannabina, YTUBH007, YTUZZ027, YTUHZ044 and YTUHZ045. Their complete genomes were sequenced and assembled at the replicon level. Each strain represents a different species according to the average nucleotide identity (ANI) values calculated using the whole-genome sequences; furthermore, except for YTUBH007, which was classified as Rhizobium binae, the remaining three strains were identified as new candidate species. A single symbiotic plasmid sized 345-402 kb containing complete nod, nif, fix, T3SS and conjugal transfer genes was detected in each strain. The high ANI and amino acid identity (AAI) values, as well as the close phylogenetic relationships among the entire symbiotic plasmid sequences, indicate that they have the same origin and the entire plasmid has been transferred among different Rhizobium species. These results indicate that S. cannabina stringently selects a certain symbiosis gene background of the rhizobia for nodulation, which might have forced the symbiosis genes to transfer from some introduced rhizobia to the related native or local-condition-adapted bacteria. The existence of almost complete conjugal transfer related elements, but not the gene virD, indicated that the self-transfer of the symbiotic plasmid in these rhizobial strains may be realized via a virD-independent pathway or through another unidentified gene. This study provides insight for the better understanding of high-frequency symbiotic plasmid transfer, host-specific nodulation and the host shift for rhizobia.

Population genomics of Australian indigenous Mesorhizobium reveals diverse nonsymbiotic genospecies capable of nitrogen-fixing symbioses following horizontal gene transfer

Mesorhizobia are soil bacteria that establish nitrogen-fixing symbioses with various legumes. Novel symbiotic mesorhizobia frequently evolve following horizontal transfer of symbiosis-gene-carrying integrative and conjugative elements (ICESyms) to indigenous mesorhizobia in soils. Evolved symbionts exhibit a wide range in symbiotic effectiveness, with some fixing nitrogen poorly or not at all. Little is known about the genetic diversity and symbiotic potential of indigenous soil mesorhizobia prior to ICESym acquisition. Here we sequenced genomes of 144 Mesorhizobium spp. strains cultured directly from cultivated and uncultivated Australian soils. Of these, 126 lacked symbiosis genes. The only isolated symbiotic strains were either exotic strains used previously as legume inoculants, or indigenous mesorhizobia that had acquired exotic ICESyms. No native symbiotic strains were identified. Indigenous nonsymbiotic strains formed 22 genospecies with phylogenomic diversity overlapping the diversity of internationally isolated symbiotic Mesorhizobium spp. The genomes of indigenous mesorhizobia exhibited no evidence of prior involvement in nitrogen-fixing symbiosis, yet their core genomes were similar to symbiotic strains and they generally lacked genes for synthesis of biotin, nicotinate and thiamine. Genomes of nonsymbiotic mesorhizobia harboured similar mobile elements to those of symbiotic mesorhizobia, including ICESym-like elements carrying aforementioned vitamin-synthesis genes but lacking symbiosis genes. Diverse indigenous isolates receiving ICESyms through horizontal gene transfer formed effective symbioses with Lotus and Biserrula legumes, indicating most nonsymbiotic mesorhizobia have an innate capacity for nitrogen-fixing symbiosis following ICESym acquisition. Non-fixing ICESym-harbouring strains were isolated sporadically within species alongside effective symbionts, indicating chromosomal lineage does not predict symbiotic potential. Our observations suggest previously observed genomic diversity amongst symbiotic Mesorhizobium spp. represents a fraction of the extant diversity of nonsymbiotic strains. The overlapping phylogeny of symbiotic and nonsymbiotic clades suggests major clades of Mesorhizobium diverged prior to introduction of symbiosis genes and therefore chromosomal genes involved in symbiosis have evolved largely independent of nitrogen-fixing symbiosis.

Crystallographic and X-ray scattering study of RdfS, a recombination directionality factor from an integrative and conjugative element

The recombination directionality factors from Mesorhizobium spp. (RdfS) are involved in regulating the excision and transfer of integrative and conjugative elements. Here, solution small-angle X-ray scattering, and crystallization and preliminary structure solution of RdfS from Mesorhizobium japonicum R7A are presented. RdfS crystallizes in space group P212121, with evidence of eightfold rotational crystallographic/noncrystallographic symmetry. Initial structure determination by molecular replacement using ab initio models yielded a partial model (three molecules), which was completed after manual inspection revealed unmodelled electron density. The finalized crystal structure of RdfS reveals a head-to-tail polymer forming left-handed superhelices with large solvent channels. Additionally, RdfS has significant disorder in the C-terminal region of the protein, which is supported by the solution scattering data and the crystal structure. The steps taken to finalize structure determination, as well as the scattering and crystallographic characteristics of RdfS, are discussed.

Bradyrhizobium cenepequi sp. nov., Bradyrhizobium semiaridum sp. nov., Bradyrhizobium hereditatis sp. nov. and Bradyrhizobium australafricanum sp. nov., symbionts of different leguminous plants of Western Australia and South Africa and definition of three novel symbiovars

Bradyrhizobium is a heterogeneous bacterial genus capable of establishing symbiotic associations with a broad range of legume hosts, including species of economic and environmental importance. This study was focused on the taxonomic and symbiovar definition of four strains – CNPSo 4026T, WSM 1704T, WSM 1738T and WSM 4400T – previously isolated from nodules of legumes in Western Australia and South Africa. The 16S rRNA gene phylogenetic tree allocated the strains to the Bradyrhizobium elkanii supergroup. The multilocus sequence analysis (MLSA) with partial sequences of six housekeeping genes – atpD, dnaK, glnII, gyrB, recA and rpoB – did not cluster the strains under study as conspecific to any described Bradyrhizobium species. Average nucleotide identity and digital DNA–DNA hybridization values were calculated for the four strains of this study and the closest species according to the MLSA phylogeny with the highest values being 95.46 and 62.20 %, respectively; therefore, both being lower than the species delineation cut-off values. The nodC and nifH phylogenies included strains WSM 1738T and WSM 4400T in the symbiovars retamae and vignae respectively, and also allowed the definition of three new symbiovars, sv. cenepequi, sv. glycinis, and sv. cajani. Analysis of morphophysiological characterization reinforced the identification of four novel proposed Bradyrhizobium species that are accordingly named as follows: Bradyrhizobium cenepequi sp. nov. (CNPSo 4026T=WSM 4798T=LMG 31653T), isolated from Vigna unguiculata; Bradyrhizobium semiaridum sp. nov. (WSM 1704T=CNPSo 4028T=LMG 31654T), isolated from Tephrosia gardneri; Bradyrhizobium hereditatis sp. nov. (WSM 1738T=CNPSo 4025T=LMG 31652T), isolated from Indigofera sp.; and Bradyrhizobium australafricanum sp. nov. (WSM 4400T=CNPSo 4015T=LMG 31648T) isolated from Glycine sp.