Exploring the evolutionary dynamics of Rhizobium plasmids through bipartite network analysis

Snowflake: visualizing microbiome abundance tables as multivariate bipartite graphs

Current visualizations in microbiome research rely on aggregations in taxonomic classifications or do not show less abundant taxa. We introduce Snowflake: a new visualization method that creates a clear overview of the microbiome composition in collected samples without losing any information due to classification or neglecting less abundant reads. Snowflake displays every observed OTU/ASV in the microbiome abundance table and provides a solution to include the data's hierarchical structure and additional information obtained from downstream analysis (e.g., alpha- and beta-diversity) and metadata. Based on the value-driven ICE-T evaluation methodology, Snowflake was positively received. Experts in microbiome research found the visualizations to be user-friendly and detailed and liked the possibility of including and relating additional information to the microbiome's composition. Exploring the topological structure of the microbiome abundance table allows them to quickly identify which taxa are unique to specific samples and which are shared among multiple samples (i.e., separating sample-specific taxa from the core microbiome), and see the compositional differences between samples. An R package for constructing and visualizing Snowflake microbiome composition graphs is available at https://gitlab.com/vda-lab/snowflake.

Dissecting the HGT network of carbon metabolic genes in soil-borne microbiota

The microbiota inhabiting soil plays a significant role in essential life-supporting element cycles. Here, we investigated the occurrence of horizontal gene transfer (HGT) and established the HGT network of carbon metabolic genes in 764 soil-borne microbiota genomes. Our study sheds light on the crucial role of HGT components in microbiological diversification that could have far-reaching implications in understanding how these microbial communities adapt to changing environments, ultimately impacting agricultural practices. In the overall HGT network of carbon metabolic genes in soil-borne microbiota, a total of 6,770 nodes and 3,812 edges are present. Among these nodes, phyla Proteobacteria, Actinobacteriota, Bacteroidota, and Firmicutes are predominant. Regarding specific classes, Actinobacteria, Gammaproteobacteria, Alphaproteobacteria, Bacteroidia, Actinomycetia, Betaproteobacteria, and Clostridia are dominant. The Kyoto Encyclopedia of Genes and Genomes (KEGG) functional assignments of glycosyltransferase (18.5%), glycolysis/gluconeogenesis (8.8%), carbohydrate-related transporter (7.9%), fatty acid biosynthesis (6.5%), benzoate degradation (3.1%) and butanoate metabolism (3.0%) are primarily identified. Glycosyltransferase involved in cell wall biosynthesis, glycosylation, and primary/secondary metabolism (with 363 HGT entries), ranks first overwhelmingly in the list of most frequently identified carbon metabolic HGT enzymes, followed by pimeloyl-ACP methyl ester carboxylesterase, alcohol dehydrogenase, and 3-oxoacyl-ACP reductase. Such HGT events mainly occur in the peripheral functions of the carbon metabolic pathway instead of the core section. The inter-microbe HGT genetic traits in soil-borne microbiota genetic sequences that we recognized, as well as their involvement in the metabolism and regulation processes of carbon organic, suggest a pervasive and substantial effect of HGT on the evolution of microbes.

Nodules of wild legumes as unique natural hotspots of antibiotic resistance genes

Root nodules (RN) of legumes have distinct microenvironment from their symbiotic roots and surrounding soils. The rhizobia can withstand the host-produced phytoalexins and antimicrobial compounds. We thus hypothesize that the wild legume RN may develop unique natural resistome and be antibiotic resistance gene (ARG) hotspots. In this study, in comparison with rhizosphere soil (RS) and bulk soil (BS), we characterized the feature of antibiotic resistance in the RN of two wild legumes, Medicago polymorpha and Astragalus sinicus, by metagenomics. It was shown that the total relative abundance of ARGs followed the order of RN > RS > BS for both legumes. ARGs encoding antibiotic efflux pump predominated in all samples with increased proportion from BS to RN samples for both legumes. Totally 275 ARG subtypes were detected, and diversity of ARGs in RN was significantly lower than in BS samples for both legumes. 32 and 25 unique ARGs subtypes were detected in RN of both legumes. Bacterial community played a key role in shaping nodule-associated resistome because both ARG profiles and bacterial community differed greatly among BS, RS and RN. Rhizobia potentially hosted 10 and 15 ARGs subtypes for both legumes. The number and proportion of plasmid- and ARG-carrying contigs (ACCs) were higher in RN than in BS. Host tracking analysis of plasmid-ACCs suggests that proportion of rhizobial bacteria identified as their hosts decreased from BS to RN samples. No plasmid-ACCs with multiple ARGs were observed in BS samples, whereas they were detected in RN samples of both legumes. Our study showed that even wild legume nodules are unique natural ARG hotspots and enough attention should be paid to the dissemination risk of ARGs posed by globally produced legume crops.

Why are rhizobial symbiosis genes mobile?

Rhizobia are one of the most important and best studied groups of bacterial symbionts. They are defined by their ability to establish nitrogen-fixing intracellular infections within plant hosts. One surprising feature of this symbiosis is that the bacterial genes required for this complex trait are not fixed within the chromosome, but are encoded on mobile genetic elements (MGEs), namely plasmids or integrative and conjugative elements. Evidence suggests that many of these elements are actively mobilizing within rhizobial populations, suggesting that regular symbiosis gene transfer is part of the ecology of rhizobial symbionts. At first glance, this is counterintuitive. The symbiosis trait is highly complex, multipartite and tightly coevolved with the legume hosts, while transfer of genes can be costly and disrupt coadaptation between the chromosome and the symbiosis genes. However, horizontal gene transfer is a process driven not only by the interests of the host bacterium, but also, and perhaps predominantly, by the interests of the MGEs that facilitate it. Thus understanding the role of horizontal gene transfer in the rhizobium-legume symbiosis requires a 'mobile genetic element's-eye view' on the ecology and evolution of this important symbiosis. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

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

Members of the Mesorhizobium genus are soil bacteria that often form nitrogen-fixing symbioses with legumes. Most characterised Mesorhizobium spp. genomes are ~8 Mb in size and harbour extensive pangenomes including large integrative and conjugative elements (ICEs) carrying genes required for symbiosis (ICESyms). Here, we document and compare the conjugative mobilome of 41 complete Mesorhizobium genomes. We delineated 56 ICEs and 24 integrative and mobilizable elements (IMEs) collectively occupying 16 distinct integration sites, along with 24 plasmids. We also demonstrated horizontal transfer of the largest (853,775 bp) documented ICE, the tripartite ICEMspSym^AA22. The conjugation systems of all identified ICEs and several plasmids were related to those of the paradigm ICESym ICEMlSym^R7A, with each carrying conserved genes for conjugative pilus formation (trb), excision (rdfS), DNA transfer (rlxS) and regulation (fseA). ICESyms have likely evolved from a common ancestor, despite occupying a variety of distinct integration sites and specifying symbiosis with diverse legumes. We found extensive evidence for recombination between ICEs and particularly ICESyms, which all uniquely lack the conjugation entry-exclusion factor gene trbK. Frequent duplication, replacement and pseudogenization of genes for quorum-sensing-mediated activation and antiactivation of ICE transfer suggests ICE transfer regulation is constantly evolving. Pangenome-wide association analysis of the ICE identified genes potentially involved in symbiosis, rhizosphere colonisation and/or adaptation to distinct legume hosts. In summary, the Mesorhizobium genus has accumulated a large and dynamic pangenome that evolves through ongoing horizontal gene transfer of large conjugative elements related to ICEMlSym^R7A.

Genetic Variation in Host-Specific Competitiveness of the Symbiont Rhizobium leguminosarum Symbiovar viciae

Legumes of the Fabeae tribe form nitrogen-fixing root nodules resulting from symbiotic interaction with the soil bacteria Rhizobium leguminosarum symbiovar viciae (Rlv). These bacteria are all potential symbionts of the Fabeae hosts but display variable partner choice when co-inoculated in mixture. Because partner choice and symbiotic nitrogen fixation mostly behave as genetically independent traits, the efficiency of symbiosis is often suboptimal when Fabeae legumes are exposed to natural Rlv populations present in soil. A core collection of 32 Rlv bacteria was constituted based on the genomic comparison of a collection of 121 genome sequences, representative of known worldwide diversity of Rlv. A variable part of the nodD gene sequence was used as a DNA barcode to discriminate and quantify each of the 32 bacteria in mixture. This core collection was co-inoculated on a panel of nine genetically diverse Pisum sativum, Vicia faba, and Lens culinaris genotypes. We estimated the relative Early Partner Choice (EPC) of the bacteria with the Fabeae hosts by DNA metabarcoding on the nodulated root systems. Comparative genomic analyses within the bacterial core collection identified molecular markers associated with host-dependent symbiotic partner choice. The results revealed emergent properties of rhizobial populations. They pave the way to identify genes related to important symbiotic traits operating at this level.

Population Genomics of Bacterial Plant Pathogens

Population genomics is transforming our understanding of pathogen biology and evolution, and contributing to the prevention and management of disease in diverse crops. We provide an overview of key methods in bacterial population genomics and describe recent work focusing on three topics of critical importance to plant pathology: (i) resolving pathogen origins and transmission pathways during outbreak events, (ii) identifying the genetic basis of host specificity and virulence, and (iii) understanding how pathogens evolve in response to changing agricultural practices.[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.