Liangzhi Li1,2
Yongjun Liu3*
Qinzhi Xiao4Zhipeng Xiao5
Delong Meng1,2Zhaoyue Yang1,2Wenqiao Deng6*
Huaqun Yin1,2
Zhenghua Liu1,2*- 1School of Minerals Processing and Bioengineering, Central South University, Changsha, China
- 2Key Laboratory of Biometallurgy of Ministry of Education, Central South University, Changsha, China
- 3Hunan Tobacco Science Institute, Changsha, China
- 4Yongzhou Tobacco Company of Hunan Province, Yongzhou, China
- 5Hengyang Tobacco Company of Hunan Province, Hengyang, China
- 6Changsha Institute of Agricultural Science, Changsha, China
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.
1. Introduction
Microorganisms are the foundation of the Earth’s biosphere and play an integral and unique role in various ecological processes and functions, where they interact to form complex functional networks (Zhou et al., 2010). Soil is an important component of the global carbon cycle and is critical to climate change mitigation. Almost all life on Earth cannot leave the living soils, which emphasizes the significance of soil-borne microbiota for essential life-supporting processes (e.g., carbon and nitrogen cycling). The soil organic carbon (SOC) reservoir (about 1,500 Gt) is supposed to be greater than the sum of the carbon stocks in the air and global flora (about 560 Gt) (van Elsas, 2019; Jassey et al., 2022). Microorganisms have much higher growth rates and carbon turnover rates than plants. The size of the soil organic carbon pool depends to a large extent on microorganisms, as their growth and activity balance the accumulation of organic carbon and its release through the decomposition of plant die-offs. Soil-borne microbes control the kinetics of soil carbon transformation by converting carbon from plants, incorporating carbon resources to increase biomass, and breaking down terrestrial organic compounds. For instance, strong biological methane-oxidizing activities in agricultural soils can lead to the emissions of biogenic CO2 linked to CH4 oxidation by a large biodiversity of methanotrophs (Cappelletti et al., 2016). Globally, soil algae absorb about 3.6 Pg of carbon per year, which is 6.4% of annual terrestrial net primary productivity (NPP) and equivalent to 31% of global anthropogenic carbon emissions (Jassey et al., 2022). Soil-borne microflora is not only dynamic temporally but also varied geographically. As a result, fluctuations in the activities and quantity of the communities that constitute such soil-borne microflora are frequently seen (Fierer, 2017). The heterogeneity in soil-borne microbiota composition is primarily caused by the spatial heterogeneity of amounts and concentrations of, for example, nutrition, mineral resources, pH, and moisture in the soil mass (Fierer, 2017; van Elsas, 2019), as well as shaped by genetic recombination and gene-specific selection processes (Crits-Christoph et al., 2020). Besides, the carbon resource abundance and diversity in soil have been proven to correlate the ecological certainty during the bio-control of microbe-induced plant disease (Zhou et al., 2023).
Horizontal gene transfer (HGT), or the interchange of genetic material across phylogenetic clades, is thought to be an efficient strategy for dispersing reproductive fitness for prokaryotic and eukaryotic cells (Huang, 2013; Daubin and Szöllősi, 2016; Li et al., 2022). Employing transmitting mobile genetic elements (MGEs) like plasmids, viruses, transposons, and gene transfer agents (GTA, tailed phage-like entities capable of packing and transferring random pieces of the host genome) as well as by direct absorption and assimilation of naked DNA by homologous or unauthorized recombination, new genes are transferred through HGT (Thomas and Nielsen, 2005; Lang and Beatty, 2007). As a result, the origins of specific genes within a particular species can vary, and patterns of gene exchange between near and remotely affiliated lineages can be seen on various genomes (Gogarten et al., 2002; Koonin, 2005; Kunin et al., 2005). HGT is commonly described based on contradicting phylogenetic trees upon genomic comparison (Delsuc et al., 2005; Li et al., 2019). In recent decades, the evolution and HGT processes of genetic traits have been investigated in a variety of ecosystems, where multiple factors (e.g., environmental conditions, ancestral genome sizes) might have influenced the frequency of HGT during evolution (Huang, 2013; Daubin and Szöllősi, 2016; Li et al., 2022). Stable HGT flux induced under selection was also suggested to enhance microbial interplay’s structural stability and thereby maintain microbial communities’ equilibrium (Fan et al., 2018). For instance, it was reported that 9.6% of the genes within a prokaryotic genome were recently acquired on average (Kloesges et al., 2011), while in E. coli 18% of genes were recently acquired via HGT (Lawrence and Ochman, 2002), and in Rickettsiales 25% of core genes were recently transferred (Hernández-López et al., 2013).
Even though the terrestrial area is a continually fluctuating and difficult place (Fierer, 2017), factors like limited nutritional and ion intensities that endorse competency, clay deposits that sustain the perseverance of bacteriophages and naked DNA, and the capacity of subsoil microbes to aggregate for genetic exchange suggest that HGT processes could be performing in this globalized context. Consistently, the transmission of antibiotic-resistance genes in the soil via MGEs (e.g., bacteriophage) has been highlighted as a public concern (Gallo et al., 2019; Zhang et al., 2022). Besides, soil is stratified geographically and varied physiologically and biologically on various dimensions, which offers a unique habitat for microbiota of diverse background. Many agronomic factors, such as pH, saltiness, warmth, and humidity, influence the organization of the microbiota inside the complicated soil-borne biomass (Frindte et al., 2019). The microscopic society’s shape and function in soils could also be influenced by temporal changes such as meteorological conditions, rhizosphere exudation, as well as other periodical supplies of plant organic matter (Zheng et al., 2019; Chernov and Zhelezova, 2020). In natural soil settings, the first investigative research on HGT among microorganisms was conducted in the 1970s (Weinberg and Stotzky, 1972; Graham and Istock, 1978). Since then, further investigations have employed fieldwork and microcosm experiments to evaluate ecological HGT (van Elsas et al., 2006). Plant-associated soils and biofilm communities are well-known hotspots for HGT due to the great genetic variety on such a restricted geographical level (Fan et al., 2018). It was found that the genetic factors in HGT are more abundant in the rhizosphere than in bulk soil, therefore HGT may aid in the development of rhizosphere competence (Vieira et al., 2020). Previous bipartite network analyses have provided evidence of genetic exchange, plasmid fusion and fission, exogenetic plasmid transfer, and environmental adaptation of the soil bacteria such as Rhizobium (Corel et al., 2018; Li L. et al., 2020; Li X. et al., 2020). In addition, plasmid transmission from transplanted P. fluorescens to native gram-negative rhizobacteria in soil has been proven to happen in natural settings (Balthazar et al., 2021).
Many terrestrial subterranean microbes have also been found to contain MGEs, and some of them have already shown conjugative activity in a lab environment or even in bulk soil. Solitary (Feng et al., 2007) or multiple (Brockman et al., 1989) plasmids could be present in subterranean isolates, with big plasmids that are more apt to own recombinant functionality seeming to predominate (Fredrickson et al., 1988). Microorganisms buried deeper underground have been shown to harbor extremely large (>150 kb) conjugative plasmids at a higher rate than microorganisms from shallower subsurface soils (Ogunseitan et al., 1987; Fredrickson et al., 1988). Degradative genetic determinants are commonly found in subsurface large plasmids, which significantly increase the host’s metabolism adaptability (Romine et al., 1999; Basta et al., 2004). Elements similar to recombinases are also present in several accessible plasmids, implying the possibility of incorporation and excision from the soil microbial genome (Romine et al., 1999). Among the sequenced soil subsoil genomes, G. thermodenitrificans was discovered to match 11 and 3% of its genetic makeup with Bacillus sp. and other Firmicutes, respectively. Around 2.7% of the genes in G. thermodenitrificans NG80-2 exclusively occur in distant relatives and may have undergone HGT. This includes two groups of proteins associated with nitrogen use that seem to have originated separately. The addition of such genes and the catabolic plasmid pLW1071 has significantly improved the metabolic adaptability in nutrient-poor conditions (Feng et al., 2007). Furthermore, the introduction of MGEs into synthesized microbial communities (SynComs) was considered an effective way for manipulating the microbial community for applicable purposes (e.g., carbon storage) (Song et al., 2014; Tsoi et al., 2019; de Lorenzo, 2022).
Yet, unlike simple bacterial mixes in a lab, the complexities of HGT mechanisms in wild areas, including the heavily heterogeneous populations, still need to be fully deciphered. In situ critical cellular events, additional chemical triggers for genetic exchange or absorption, and natural factors like cell-MGE interactions are examples of the complexities. In addition, there needs to be more focus on the environmental implications, such as how these HGT activities affect microbial communities’ capacity to adapt to wild areas, and how these HGT mechanisms reflect the features of the particular ecosystem. Though there is a growing interest in elucidating the underlying microbial mechanisms driving soil carbon transformation, stabilization, and release processes, many unknowns remain. Nonetheless, many earlier studies mainly concentrated on HGT episodes of particular lineages and on significant issues such as adaption-associated processes (Li et al., 2019; Li L. et al., 2020; Li X. et al., 2020; Li et al., 2021). Regardless of these efforts, there is so much to be learned in a holistic mode about the quantity and effects of HGT-driven genetic makeup in soil-borne microbiota that can confer a variety of carbon assimilation and dissimilation functions (e.g., central carbon cycles, fatty acid biosynthesis). This poses significant challenges to understanding gene flow in terrestrial microbial communities. As a result, it remains to be seen whether there are inter-microbe HGT(s) that have promoted the development of soil-borne microbiota in a variety of conditions. If these HGT markers exist, how do they work to counteract environmental stressors, and how widely disseminated are these markers throughout the genetic material of soil-borne microbiota?
A variety of natural systems, including protein expression (Alon, 2007), biochemical processes (Pál et al., 2005), biomolecule interplay (Jeong et al., 2001), contradicting evolutionary indications (Huson and Bryant, 2006), and ecosystem dynamics (Rezende et al., 2007), are being modeled using network infrastructure (Proulx et al., 2005). In theory, the network is capable of better displaying the patterns of microbial genomic evolution (Doolittle and Bapteste, 2007). HGT networks are a distinctive form of sharing gene networks. They are intended to investigate trends in genetic dispersion brought on by HGT throughout ecological evolution.
Under this context, we performed BLASTP-driven screenings and HGT network constructions to discover and investigate potential inter-microbe HGT genes linked to carbon metabolism throughout genomic sequences of all accessible strains isolated from soils, followed by network constructions. Our results suggested that the HGT episodes may significantly contribute to genetic imports that contribute to the soil-borne microbiota’s increased carbon metabolic versatility and adaptivity. Our knowledge of microbial interplay and the adaptable development of microbes to deal with various situations have therefore been enhanced as a result of this work. However, further experimental work is still necessary to evaluate the occurrence of HGT in the vast and unexplored environment of the terrestrial subsurface.
2. Results
2.1. Overall statues
We first queried and retrieved all the items in the public database (Genbank/IMG) with the keyword “soil” within the “biosample” or “habitat” regions. We then chose and downloaded from the public database (Genbank/IMG) the out-coming high-quality genomes of soil microbial isolates (n = 764, see Supplementary Table S1 for details at https://doi.org/10.6084/m9.figshare.22154828.v1) for downstream analyses. Putative microbial horizontally transferred genes (HTGs) in these genomes were identified employing the BLASTP-based IMGAP pipeline under default mode (Markowitz et al., 2010), which resulted in a total of 37,481 HTG entries, followed by manual pickup of HTGs that are related to carbon metabolism. These processes eventually produce a total of 4,554 HTG entries related to carbon metabolism, as listed in Supplementary Table S2 at https://doi.org/10.6084/m9.figshare.22154828.v1.
The significance of HGT to biological diversification could be better understood through a network study of gene patterns amongst genomes. A directed phylogenomic network can be built using HGT inference techniques that consider the provider and receptor of the lateral gene delivery event (Popa et al., 2011). Similar to the methods used in previous work (Beiko et al., 2005), we established the phylogenomic HGT networks of carbon metabolic genes in soil-borne microbiota using HGT occurrences predicted to examine the HGT component within microbiological diversification (Figure 1). A baseline species (tested genomes) evolutionary tree’s exterior and internal nodes make up the HGT network’s vertex. Between the nodes, edges in the network refer to alleged gene transfer occurrences. If a gene relative is transferred across the nodes through a suspected HGT episode, these nodes are interconnected in the HGT networks. The HGT network’s edge weight is correlated with the amount of horizontally shared genes. Clusters are formed when specific subgroups of organisms have stronger connections among themselves than with other factions in the networks (Guimerà and Nunes Amaral, 2005; Gallos et al., 2007).
Figure 1. An overall directed horizontal gene transfer (HGT) network of carbon metabolic genes harbored by the soil-borne microbiota. Nodes represent genomes of species. The colors are marked according to the node’s taxonomy (left, phylum; right, class). The edges represent HGT events and are directed from the donor to the recipient.
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 (Figure 1). Among the nodes, phyla Proteobacteria (34.65%), Actinobacteriota (26.68%), Bacteroidota (9.28%), and Firmicutes (7.89%) are predominant, followed by Cyanobacteria (1.54%) and Ascomycota (1.30%). Regarding specific classes, Actinobacteria (26.26%), Gammaproteobacteria (16.1%), Alphaproteobacteria (14.15%), Bacteroidia (7.42%), Actinomycetia (6%), Betaproteobacteria (4.08%), Clostridia (3.31%), and Bacili (3.25%) are dominant in the overall HGT network of carbon metabolic genes. The average path length (the average of the shortest route between any two nodes in the network) of the total HGT network is 1.55, indicating that the related carbon metabolic genes could be easily and rapidly transported through tested soil-borne microbiota. The HGT network’s thickest (weight > = 5) edges are mainly formed by cross-domain HGT events from archaeal phylum Haloarchaea to Bacteria (26/36 entries, 74.3%), with HGT from Methanosarcina acetivorans C2A (Haloarchaea) to Streptomyces yangpuensis fd2-tb (Actinobacteriota) having the highest HGT frequency (n = 14) (see Supplementary Table S3 at https://doi.org/10.6084/m9.figshare.22154828.v1). The chimeric nature and high externalization rate (especially on carbohydrate biosynthesis genes) of Haloarchaea have been highlighted previously (Galagan et al., 2002; Nelson-Sathi et al., 2012; Garushyants et al., 2015).
In Kyoto Encyclopedia of Genes and Genomes (KEGG) functional assignments, glycosyltransferases (18.5%), glycolysis/gluconeogenesis (8.8%), carbohydrate-related transporters (7.9%), fatty acid biosynthesis (6.5%), benzoate degradation (3.1%), and butanoate metabolism (3.0%) are most frequently identified (Figure 2). These results suggest that soil-borne microbiota is involved in a range of carbon metabolic processes, particularly those related to the breakdown and transportation of carbohydrates and fatty acids. This may have important implications for soil health and ecosystem function.
Figure 2. Sankey diagram showing the distribution of donor taxa, recipient taxa, and function categories of carbon metabolic genes within the soil-borne microbiota HGT network, with edge thickness and proportion indicating the relative frequency of transfer events.
We then used two classical network centralities, degree and betweenness, to examine what taxonomic groups are holding and shaping the HGT network (see Supplementary Tables S4, S5 at https://doi.org/10.6084/m9.figshare.22154828.v1). The most commonly used centrality measurement, degree, counts the edges that link a specific node. In a directed graph, the arrows are directional, pointing from one node (HGT donor) to another (HGT acceptor), so the number of arrows pointing to each node is its entry degree. In contrast, the number of arrows pointing away from this node is its outgoing degree. Of the top 30 nodes having the highest degree values (>25), class Actinomycetia occupies the most (43.2%) of these nodes (see Supplementary Table S4 at https://doi.org/10.6084/m9.figshare.22154828.v1), followed by classes Deltaproteobacteria (13.3%) and Betaproteobacteria (13.3%). In comparison, most other nodes (92.5%) in the overall HGT network of carbon metabolic genes in soil-borne microbiota display a degree value <5. In particular, the highest degree values in correspondent classes were seen in the genomes of Kribbella soli VKM Ac-2,540 (98, Actinomycetia), Anaeromyxobacter sp. Red232 (55, Deltaproteobacteria), Paenibacillus piri MS74 (41, Bacilli), Pseudoduganella violaceinigra DSM 15887(36, Betaproteobacteria), and Pyrinomonas methylaliphatogenes K22 (35, Blastocatellia).
Another centrality measurement, betweenness, measures how frequently a specific node is located on the shortest path possible between any pair of nodes within the network. Nodes with large betweenness centralities act as essential reservoirs or communicators for the common-good genetic codes. Of the top 30 nodes having the highest betweenness centrality values (>20) in the overall HGT network of carbon metabolic genes in soil-borne microbiota (see Supplementary Table S5 at https://doi.org/10.6084/m9.figshare.22154828.v1), class Actinomycetia also occupies most (60.0%) of these nodes, followed by class Deltaproteobacteria (10.0%). Primarily exhibiting the highest betweenness centrality values in correspondent classes are genomes of Jiangella anatolica GTF31 (558, Actinomycetia), Pyrinomonas methylaliphatogenes K22 (195.5, Blastocatellia), Saccharicrinis fermentans DSM 9555 (168, Bacteroidia), Methylovulum miyakonense HT12 (89, Gammaproteobacteria), and Nitrosospira multiformis ATCC 25196 (60.5, Betaproteobacteria). In comparison, most other nodes (98.2%) in the overall HGT network of carbon metabolic genes in soil-borne microbiota display a betweenness centrality value <5. On the other hand, nodes that display lower betweenness centralities either occupy peripheral positions in the network or form their linked constituents.
We further investigated the network’s community composition by segmenting it into sub-network according to the major phyla present in the overall HGT network (Figure 3). The sub-networks were created by extracting corresponding taxonomic (e.g., Proteobacteria) nodes acting as HGT donors and acceptors, respectively, as well as extracting their direct edges (Figure 3A). The donor and acceptor sub-networks of the same phyla generally exhibited similar network properties in comparison (Table 1). The same possesses for sub-network construction were applied to nodes of phyla Actinobacteria and Bacteroidota (Figures 3B,C).
Figure 3. The sub-networks created by extracting corresponding taxonomic (e.g., Proteobacteria) nodes acting as HGT donors and acceptors, respectively, as well as extracting their direct edges according to the major phyla present in the overall HGT network: (A) Proteobacteria; (B) Actinobacteriota; (C) Bacteroidota. Nodes represent genomes of species. The edges represent HGT events and are directed from the donor to the recipient. The colors are marked according to the node’s taxonomy.
Table 1. Network properties of the overall and sub-networks of carbon metabolic HGT detected in soil-borne microbiota.
In the sub-network that Proteobacteria members act as HGT donors (Figure 3A, left), Actinomycetia (4.51%) and Bacili (1.13%) are the most frequently observed classes outside the phyla Proteobacteria, indicating that classes Actinomycetia and Bacili are the major cross-phylum recipients of the genetic goods carried by Proteobacteria; whereas in the sub-network where Proteobacteria acting as HGT acceptor (Figure 3A, right), Actinobacteria (7.39%) and Bacteroidia (3.5%) are the most frequently observed classes outside the phyla Proteobacteria, suggesting that classes Actinobacteria and Bacteroidia are major cross-phylum donors of carbon metabolic genes to Proteobacteria.
Secondly, in the sub-network that Actinobacteriota members act as HGT donors (Figure 3B, left), Betaproteobacteria (1.61%), Deltaproteobacteria (1.44%), and Alphaproteobacteria (1.02%) the most frequently observed classes outside the phyla Actinobacteriota, indicating that these classes are the primary cross-phylum recipients of the genetic goods carried by Actinobacteriota; whereas in the sub-network where Actinobacteriota acting as HGT acceptor (Figure 3B, right), Alphaproteobacteria (7.16%), superior to Gammaproteobacteria (6.18%), and Betaproteobacteria (1.21%), stands for the most frequently observed classes outside the phyla Actinobacteriota, which indicates that classes from Proteobacteria like Alphaproteobacteria are major cross-phylum donors of carbon metabolic genes to Actinobacteriota.
Lastly, in the sub-network where Bacteroidota members act as HGT donors (Figure 3C, left), Deltaproteobacteria (2.84%), Betaproteobacteria (1.99%), and Alphaproteobacteria (1.99%) are the most frequently observed classes outside the phyla Bacteroidota, indicating that these classes are the major cross-phylum recipients of the genetic goods carried by Bacteroidota; whereas in the sub-network where Bacteroidota acting as HGT acceptor (Figure 3C, right), Gammaproteobacteria (6.73%), Alphaproteobacteria (3.51%), and Bacilli (2.5%) are the most frequently observed classes outside the phyla Bacteroidota, suggesting that these classes are major cross-phylum donors of carbon metabolic genes to Bacteroidota.
Network properties such as the node number, edge number, average (weighted) degree, as well as average path length were higher in the “acceptor” network than in the “donor” one in the corresponding networks of Proteobacteria and Actinobacteria, respectively. While the connector component (number of subgraphs in which each pair of nodes is connected) shows the opposite trend (higher in the “acceptor”). Average eigenvector centralities (that measure the transitive influence of nodes) decline from the Proteobacteria network to that of Actinobacteria, and, finally, Bacteroidota. Other network properties like modularity (measuring the structural strength and density of a network community) show only slight variations among groups (Table 1).
2.2. Carbon metabolic functions
We further created sub-networks that extracted HGT events (edges) corresponding to specific carbon metabolic functions (e.g., core carbon metabolic pathway, fatty acid biosynthesis) and nodes acting as HGT donors and acceptors, respectively, as well as extracting their direct edges (Figure 4).
Figure 4. The sub-networks created by extracting HGT events (edges) corresponding to specific carbon metabolic functions: (A) core carbon metabolism; (B) butanoate metabolism; (C) lipid metabolism; (D) xenobiotics biodegradation and metabolism; (E) glycosyltransferase; (F) transport protein family. Nodes represent genomes of species. The edges represent HGT events and are directed from the donor to the recipient.
In the sub-network reflecting HGT episodes of core carbon metabolism (700 edges/genes, 882 nodes/genomes), genes attributed to pathways of glycolysis/gluconeogenesis (10.14%), butanoate metabolism (7.29%), pyruvate metabolism (7.29%), and galactose metabolism (7.00%) are predominant. Microbial classes Actinobacteria (27.66%), Gammaproteobacteria (13.83%), and Alphaproteobacteria (12.59%) are mostly seen (Figure 4A).
Alcohol dehydrogenase that metabolizes alcoholic substances (with 230 entries) is dominant in the current sub-network of core carbon metabolism (Figure 4A, marked with red rectangles) and ranks third in the list of most frequently identified carbon-metabolism-related horizontal transferred gene (HTG) from the soil-borne microbiota (see Supplementary Table S2 at https://doi.org/10.6084/m9.figshare.22154828.v1). Alcohol dehydrogenase exhibits a wide range of substrate specificity additionally, oxidizing mainly primary and secondary aliphatic alcohols when utilizingNAD+ as a co-substrate. It is also able to reduce aldehydes and ketones (Tsigos et al., 1998). Accordingly, HGT of alcohol dehydrogenase genes from prokaryotes and unicellular eukaryotes was also reported previously (Espinosa and Paz-Y-Miño-C, 2012; McCarthy and Fitzpatrick, 2016).
We further extracted the sub-networks of butanoate metabolism (53 edges, 88 nodes). In the sub-network, HGT genes encoding acetolactate synthase (28.3%), succinate-semialdehyde dehydrogenase (22.64%), and 3-hydroxybutyryl-CoA dehydrogenase (9.43%) are predominant. Microbial classes Actinobacteria (26.14%), Actinomycetia (21.59%), and Alphaproteobacteria (13.64%) are mostly seen (Figure 4B).
Likewise, in the sub-network reflecting HGT episodes of lipid metabolism (374 edges, 543 nodes), genes attributed to fatty acid biosynthesis (28.88%), primary bile acid biosynthesis (2.94%), and fatty acid degradation (2.67%) are predominant. Microbial classes Actinobacteria (29.65%), Alphaproteobacteria (15.47%), Actinomycetia (14.0%), Gammaproteobacteria (9.02%), and Betaproteobacteria (7.73%) are mostly seen at the same time (Figure 4C).
Unexpectedly, the pimeloyl-ACP methyl ester carboxylesterase (with 231 entries) is dominant in the current sub-network of lipid metabolism (Figure 4C, marked with red rectangles) and ranks second in the list of most frequently identified carbon-metabolism-related HTGs from the soil-borne microbiota (see Supplementary Table S2 at https://doi.org/10.6084/m9.figshare.22154828.v1). Pimeloyl-ACP methyl ester carboxylesterase catalyzes the hydrolysis of the ester bonds of pimeloyl-ACP esters that allow the synthesis of pimeloyl-ACP via the fatty acid synthetic pathway (Lin et al., 2010). Previous reports have also found pimeloyl-ACP methyl ester carboxylesterase located in MGEs, such as the plasmids of the psychrotolerant Polaromonas spp. isolated from Arctic and Antarctic glaciers (Ciok et al., 2018). It follows that 3-oxoacyl-ACP reductase (FabG) (with 74 entries) ranks third in the list of most frequently identified carbon-metabolism-related HTG from the soil-borne microbiota (see Supplementary Table S2 at https://doi.org/10.6084/m9.figshare.22154828.v1). FabG is the key enzyme directly responsible for the synthesis of 3-hydroxyacyl-ACPs in the fatty acid synthesis elongation cycle (Yu et al., 2019). As shown in previous studies, FabG has also been frequently observed in mobile genetic elements, like the plasmids of plant growth-promoting bacteria Azospirillum brasilense (Yu et al., 2019), and the plastid genomes of Eustigmatophyte Algae (Ševcíková et al., 2019).
Similarly, in the sub-network reflecting HGT episodes of xenobiotics biodegradation and metabolism (374 edges, 543 nodes), genes attributed to benzoate degradation (47.71%), chloroalkane/chloroalkene degradation (11.93%), steroid degradation (9.17%), fluorobenzoate degradation (7.34%), dioxin degradation (6.42%), and toluene degradation (6.42%) are predominant. At the same time, microbial classes Actinomycetia (19.07%), Actinobacteria (18.56%), Gammaproteobacteria (15.98%), Alphaproteobacteria (14.43%), and Betaproteobacteria (10.31%) are mostly observed (Figure 4D).
Regarding specific protein families, glycosyltransferase family (with 363 HGT entries) ranks first overwhelmingly in the list of most frequently identified carbon-metabolism-related HGT enzymes from the soil-borne microbiota (see Supplementary Table S2 at https://doi.org/10.6084/m9.figshare.22154828.v1). While microbial classes Gammaproteobacteria (14.37%), Actinobacteria (9.58%), Alphaproteobacteria (9.58%), Actinomycetia (8.54%), Bacteroidia, (8.12%), Betaproteobacteria (6.04%), and Bacilli (5.42%) are mostly observed to encode glycosyltransferase (Figure 4E). Glycosyltransferase is capable of catalyzing the sequential transfer of glycosyl moieties to the undecaprenyl phosphate carrier lipid during the early steps of polysaccharide synthesis, as well as various other acceptors like proteins, nucleic acid, and secondary metabolites (Garrison-Schilling et al., 2014; Zhang et al., 2020). Glycosyltransferase complexes are known to be involved in cell wall biosynthesis, protein glycosylation, and primary and secondary metabolism in both microbes and plants (Chang et al., 2003; Harholt et al., 2012), whose vital roles in enhancing metabolic flexibility and adaptive advantages might have endorsed the frequent HGT-driven expansions of glycosyltransferase across microbial species. Bacterial and viral glycosyltransferases determine their surface chemistry, influencing the interplay with their hosts (Brew et al., 2010). Our finding is also in line with previous investigations, which have found glycosyltransferase in Rickettsial plasmids (El et al., 2016). Also, the vascular wilt fungus Verticillium acquired glycosyltransferase involved in the synthesis of extracellular glucans by HGT from Proteobacteria (Lacroix and Citovsky, 2016). Besides, glycosyltransferase operons acquired via HGT were found to modify the O-antigen in Salmonella enterica serovar Typhi (Kintz et al., 2017).
Transport protein families (124 edges, 191 nodes) of organic carbon compounds were also visualized by sub-network creation (Figure 4F). In the sub-network, genes attributed to multiple sugar transport system (38.71%), raffinose/stachyose/melibiose transport system (16.13%), and lipopolysaccharide transport system (12.91%) are predominant. At the same time, microbial classes Actinobacteria (24.08%), Actinomycetia (13.61%), Gammaproteobacteria (9.42%), Alphaproteobacteria (8.90%), Betaproteobacteria (7.33%), and Clostridia (6.81%) are mostly observed to encode the transporters.
Also, there are other specific carbon metabolic enzymes enriched in the HGT network (see Supplementary Table S2 at https://doi.org/10.6084/m9.figshare.22154828.v1), such as beta-lactamase family enzymes (159 entries) that degrade multiple metabolites (e.g., antibiotics) for the use of the substrates as nutrients or carbon sources (Graham et al., 2016; Crofts et al., 2018; Maatouk et al., 2022), catechol 2,3-dioxygenase-like lactoylglutathione lyase (54 entries,COG0346) that could enable the degradation of phenol or lignin (Hupert-Kocurek et al., 2012), 2-polyprenyl-6-methoxyphenol hydroxylase (53 entries,COG0654) involved in the aerobic degradation of aromatic compounds (Trias et al., 2017) as well as quinol monooxygenase YgiN (34 entries) related to aerobic redox and physiology. Mapping these HTGs of carbon metabolism onto the correspondent KEGG pathways showed that most HTGs locate on the fringe of the metabolic cycles, and the core carbon metabolic sections rarely transfer (see Supplementary Figures S1–S5 at https://doi.org/10.6084/m9.figshare.22154828.v1).
3. Discussion
An editorial in the journal Science has advocated that more thorough investigations be conducted into the appearance and mechanics of HGT episodes (“Why does lateral transfer occur in so many species and how?”) (American Association for the Advancement of Science, 2005). We thoughtfully evaluated that the current study could have provided several suggestions for this subject.
A variety of soil-borne carbon-cycle HGT events of cross-class and even cross-phylum levels were discovered with considerably high transfer frequency (Figure 1). The constructed HGT networks demonstrate that genetic exchange across microbial genera is a significant contributor to microorganisms’ biodiversity and adaptivity. The characteristic of a “small-world system” found in previous HGT networks was also observed in our study, which refers to a network that has a small diameter in terms of the number of nodes and a handful of strongly linked nodes that allow for flux to flow across the system (Proulx et al., 2005). A small-world architecture in the HGT network indicates that significantly advantageous genes that arise in any microbe could transcend taxonomic boundaries and stretch another microbe through a limited amount of HGT episodes, where genomes with high betweenness can act as a link between otherwise unconnected parts of the network and transfer genes to numerous other genomes in the ecosystem with a small number of genetic transactions. Similarly, a recent study also found that the HGT rate was increased in organisms with similar ecological distributions (Zhou et al., 2021). Besides, the observed variations in network properties (Table 1) may have multiple implications for HGT events, highlighting the importance of understanding network dynamics in studying HGT events (Kunin et al., 2005; Kloesges et al., 2011; Popa et al., 2011; Li L. et al., 2020; Li X. et al., 2020). For example, a higher number of nodes and edges in the “acceptor” network may indicate a greater potential for gene transfer due to more opportunities for contact between bacteria. Similarly, higher degrees or connector components may create more interconnectedness amongst bacterial populations and facilitate the spread of genetic material. On the other hand, a longer average path length between nodes in a network may reduce the likelihood of gene transfer, as the distance between two bacteria would be higher and the chance of interaction would be lower. Moreover, differences in eigenvector centralities among network groups may illustrate the presence of influential nodes--nodes with high centrality values that may play significant roles in mediating gene transfer. Modularity can also impact the strength and density of community structures within a network, potentially creating barriers or pathways for gene transfer.
The phylum Proteobacteria was the most abundant taxon in soil samples, accounting for an average of 30% of metagenomic sequences (Guo et al., 2018). In keeping with this, our study found that nodes of classes, including Actinobacteria as well as Proteobacteria classes, such as Alphaproteobacteria and Gammaproteobacteria are extensively present in our overall HGT network of carbon metabolism (Figure 1), implying that HGT is widespread during the evolution of species in these classes. However, the sampling density of sequenced microbial genomes might be biased toward Proteobacteria, since the predominance of Proteobacteria in the genome dataset might be responsible for the enormous rate of HGTs within this group. Previous studies have confirmed that Proteobacterial MGEs constitute the major connected component in the virulence network, and extensive gene sharing exists among Actinobacteria and Gammaproteobacteria (Tamminen et al., 2012; Jiang et al., 2017). Previous studies have also validated experimentally the possibility of MGE-mediated HGT in soil samples, in which plasmids from the donor strains Psudomonas putida KT2440, Escherichia coli MG1655, and Kluyvera sp. can be transferred to a wide range of bacterial phyla from agricultural soils (Klümper et al., 2015), including (α–ε) Proteobacteria, Acidobacteria. Actinobacteria, Bacteroidetes, Firmicutes, Fusobacter, Gemmatimonadetes, Planctomyces, spirochaetes, Candidate division TM7, Verrucomicrobia as well as Eukaryote (Zhang et al., 2015).
Externalized genes carried by MGEs could act as containers for the shared genetic pool (Norman et al., 2009). In our study, uneven functional distributions of the externalized gene are also discovered, in which HTGs are mostly incorporated into the peripheral functions of the carbon metabolic pathway (e.g., nutrient transport and dispensable reaction). In contrast, the core metabolic sections (e.g., intermediate reactions and biomass production) of putative competent significance are mostly evolutionarily unvaried (see Supplementary Figures S1–S5 at https://doi.org/10.6084/m9.figshare.22154828.v1). Similar findings were also reported previously. For example, according to research on the horizontally gained genes within the E. coli metabolic pathways (Pál et al., 2005), HGT is more common among proteins engaged in the absorption and consumption of resources than those responsible for the generation of biomass, which suggested that their function influences the HGT probability of metabolic genes in the internal metabolic pathway. This uneven functional distribution could also be explained by the “complexity hypothesis” (Jain et al., 1999; Muller et al., 2018), which proposed that proteins in a complicated system, like ribosomal systems or core metabolic cycles, are specialized to work. Decreased adaptation of the microbial recipients will arise from an HGT incidence that leads to substituting such a gene with a less-suited counterpart. The proportional influence of functionality type and the number of interactive participants on HGT incidence was also investigated, demonstrating that the “complexity theory” still holds up in the genome-scale analyses (Cohen et al., 2011).
On the other hand, soil multi-functionality is affected by the environment and by microbial community composition and diversity (Zheng et al., 2019). Environmental factors in the soil also impact the HGT processes and lead to uneven distribution. The microbial HGT rate in soil relies on environmental stress variables like surface temperatures, pH (Rochelle et al., 1989), soil composition (Richaume et al., 1992), and wetness (Richaume et al., 1989). The gene-sharing network revealed strong correlations between gene connectivity and the trailed soil variables (Zhou et al., 2010). Besides, the persistence and transportation of genetically modified bacteria were impacted by fluidity in subsoil (Trevors et al., 1990). In terms of biological variables, the existence of fungi (Sengeløv et al., 2000), protozoa (Johannes Sørensen et al., 1999), and roundworms (Daane et al., 1996, 1997) could also impact sexual plasmid transmission in soil. Moreover, the generation of leachate and root elongation looked to be the main contributors to the incidence of HGT inside the rhizosphere (Mølbak et al., 2007). Transmission rates were roughly 10 times lower in the bean and cereal rhizospheres than in the control soils (Musovic et al., 2006). Lastly, the other key factor contributing to the rise in local MGE quantities and HGT incidence in this environment is the implementation of manure, pesticides, and antimicrobials to the land (Gotz and Smalla, 1997). Correspondently, various organic degradative genes were found mobilized among soil-borne microbiota via HGT in our study (Figure 4D), whose catabolic activities might be further applied for bio-remediation of polluted environments (Top et al., 2002; Nojiri et al., 2004). It is recently confirmed that genes of microalgal origin have conferred Caenorhabditis elegans the ability to degrade cyanogenic toxins (Wang et al., 2022). Previous studies also found that functional categories “biosynthesis and degradation of surface polysaccharides and lipopolysaccharides” and “DNA regulation and modification” tend to be enriched in the HGT entries. In consistence, we found glycosyltransferase in our study as the most abundant protein encoded by HTG related to carbon metabolism. The glycosyltransferase family was reported to be involved in the glycosylation and modifications of biomolecules, including the bases in DNA, which might alter the host’s gene expression pattern (Iyer et al., 2013). Also, it was reported that glycosyltransferases were significantly enriched in horizontally transferred genes in the human gut, while soil microbiota has similar expansions of the glycosyltransferase repertoires as the gut (Lozupone et al., 2008). Another generally enriched functional category in the HGT profile is metabolite transporter (Cordero and Hogeweg, 2009; Paquola et al., 2018), which is also reflected by our results (Figure 4F).
In conclusion, the inter-microbe HGT genetic traits in soil-borne microbiota genetic sequences that we recognized through our assessments, as well as their involvement in carbon metabolism and resilience to various environmental stressors typically found in territory ecosystems, suggest a pervasive and substantial effect of HGT on the evolvement of microbes. Nevertheless, the information we have provided here is not thorough. The examples of inter-microbe HGT documented so far represent just the tip of a giant biological iceberg. Upcoming studies are expected to offer a more intriguing view of both the scope and the biological importance of HGT in microbes as increasing numbers of genomic data of higher quality are becoming accessible for yet more branches in the tree of life for microbes. Further work will be necessary to evaluate the occurrence of HGT in the vast, heterogeneous, and isolated environment of the terrestrial subsurface and to assess the full impact of gene transfer on terrestrial subsurface microbial evolution.
4. Materials and methods
We first queried and retrieved all the items in the public database (Genbank/IMG) with the keyword “soil” within the “biosample” or “habitat” regions. We have then chosen and downloaded from the public database (Genbank/IMG) the resultant high-quality 764 genomes of soil microbial isolates as listed in Supplementary Table S1 at https://doi.org/10.6084/m9.figshare.22154828.v1 for downstream analyses. The Integrated Microbial Genomes Annotation Pipeline (IMGAP) v.5.0 (Markowitz et al., 2010) under default mode was used to identify horizontally transferred genes (HTGs) in the genetic sequences of individual collected soil-borne microbiota, as conducted in previous studies (Li L. et al., 2020; Li X. et al., 2020; Li et al., 2021, 2022). It used the following criteria to determine which genes inside the trialed genomes were horizontally transferred from remote descendants: genes with the finest BLASTP matches (most significant bit scores) or over 90% of the best hits discovered beyond the phylogenetic clade of the trialed genome (i.e., from remote phylum, class, etc.) and with lower-scoring matches or no hits within the original phylogenetic clade of the trialed genome. The HTG sequences were annotated with the eggnog mapper v.2.0 (Huerta-Cepas et al., 2019). This resulted in a total of 37,481 HTG entries, followed by the pickup of HTGs that are related to carbon metabolism according to the Kyoto Encyclopedia of Genes and Genomes (KEGG) classification. These processes eventually produce a total of 4,554 HTG entries related to carbon metabolism. We then used the Gephi network visualization package1 for network visualization and exploration of carbon metabolism HTGs. Fruchterman Reingold and Openord layout approaches were used for network visualization in the Gephi. The graphical user interface of the KEGG pathway map color tool was applied for coloring map objects, including split-coloring and gradation2.
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.
Author contributions
LL and ZL conceived and designed the research. YL, QX, and ZX provided the computing sources. LL, DM, ZY, and HY analyzed the data. LL wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by the key project of Science and Technology of Hunan Branch of China National Tobacco Corporation (CD2022KJ01, 2022431021240240, HN2021KJ05), Natural Science Foundation of Changsha (No. kq2202089), Fundamental Research Funds for the Central Universities of Central South University (No. 2022ZZTS0420), Hunan International Scientific and Technological Cooperation Base of Environmental Microbiome and Application (No. 2018WK4019), and key project of Science and Technology of China National Tobacco Corporation (No. 110202201004) (JY-04).
Acknowledgments
We are grateful for resources from the High-Performance Computing Center of Central South University.
Conflict of interest
YL, QX, and ZX were employed by Zhangjiajie Tobacco Company of Hunan Province.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
Alon, U. (2007). Network motifs: theory and experimental approaches. Nat. Rev. Genet. 8, 450–461. doi: 10.1038/nrg2102
American Association for the Advancement of Science (2005). So much more to know. Science 309, 78–102. doi: 10.1126/science.309.5731.78b
Balthazar, C., Joly, D. L., and Filion, M. (2021). Exploiting beneficial Pseudomonas spp. for Cannabis production. Front. Virol. 12:833172. doi: 10.3389/fmicb.2021.833172
Basta, T., Keck, A., Klein, J., and Stolz, A. (2004). Detection and characterization of conjugative degradative plasmids in xenobiotic-degrading sphingomonas strains. J. Bacteriol. 186, 3862–3872. doi: 10.1128/JB.186.12.3862-3872.2004
Beiko, R. G., Harlow, T. J., and Ragan, M. A. (2005). Highways of gene sharing in prokaryotes. Proc. Natl. Acad. Sci. 102, 14332–14337. doi: 10.1073/pnas.0504068102
Brew, K., Tumbale, P., and Acharya, K. R. (2010). Family 6 glycosyltransferases in vertebrates and bacteria: inactivation and horizontal gene transfer may enhance mutualism between vertebrates and bacteria. J. Biol. Chem. 285, 37121–37127. doi: 10.1074/jbc.R110.176248
Brockman, F. J., Denovan, B. A., Hicks, R. J., and Fredrickson, J. K. (1989). Isolation and characterization of quinoline-degrading bacteria from subsurface sediments. Appl. Environ. Microbiol. 55, 1029–1032. doi: 10.1128/aem.55.4.1029-1032.1989
Cappelletti, M., Ghezzi, D., Zannoni, D., Capaccioni, B., and Fedi, S. (2016). Diversity of methane-oxidizing Bacteria in soils from “hot lands of Medolla” (Italy) featured by anomalous high-temperatures and biogenic CO2 emission. Microbes Environ. 31, 369–377. doi: 10.1264/jsme2.ME16087
Chang, R., Yeager, A. R., and Finney, N. S. (2003). Probing the mechanism of a fungal glycosyltransferase essential for cell wall biosynthesis. UDP-chitobiose is not a substrate for chitin synthase. Org. Biomol. Chem. 1, 39–41. doi: 10.1039/b208953j
Chernov, T. I., and Zhelezova, A. D. (2020). The dynamics of soil microbial communities on different timescales: a review. Eurasian Soil Sci. 53, 643–652. doi: 10.1134/S106422932005004X
Ciok, A., Budzik, K., Zdanowski, M. K., Gawor, J., Grzesiak, J., Decewicz, P., et al. (2018). Plasmids of psychrotolerant Polaromonas spp. isolated from arctic and Antarctic glaciers – diversity and role in adaptation to polar environments. Front. Microbiol. 9:1285. doi: 10.3389/fmicb.2018.01285
Cohen, O., Gophna, U., and Pupko, T. (2011). The complexity hypothesis revisited: connectivity rather than function constitutes a barrier to horizontal gene transfer. Mol. Biol. Evol. 28, 1481–1489. doi: 10.1093/molbev/msq333
Cordero, O. X., and Hogeweg, P. (2009). The impact of long-distance horizontal gene transfer on prokaryotic genome size. Proc. Natl. Acad. Sci. 106, 21748–21753. doi: 10.1073/pnas.0907584106
Corel, E., Méheust, R., Watson, A. K., McInerney, J. O., Lopez, P., and Bapteste, E. (2018). Bipartite network analysis of gene sharings in the microbial world. Mol. Biol. Evol. 35, 899–913. doi: 10.1093/molbev/msy001
Crits-Christoph, A., Olm, M. R., Diamond, S., Bouma-Gregson, K., and Banfield, J. F. (2020). Soil bacterial populations are shaped by recombination and gene-specific selection across a grassland meadow. ISME J. 14, 1834–1846. doi: 10.1038/s41396-020-0655-x
Crofts, T. S., Wang, B., Spivak, A., Gianoulis, T. A., Forsberg, K. J., Gibson, M. K., et al. (2018). Shared strategies for β-lactam catabolism in the soil microbiome. Nat. Chem. Biol. 14, 556–564. doi: 10.1038/s41589-018-0052-1
Daane, L. L., Molina, J. A., Berry, E. C., and Sadowsky, M. J. (1996). Influence of earthworm activity on gene transfer from Pseudomonas fluorescens to indigenous soil bacteria. Appl. Environ. Microbiol. 62, 515–521. doi: 10.1128/aem.62.2.515-521.1996
Daane, L. L., Molina, J., and Sadowsky, M. J. (1997). Plasmid transfer between spatially separated donor and recipient bacteria in earthworm-containing soil microcosms. Appl. Environ. Microbiol. 63, 679–686. doi: 10.1128/aem.63.2.679-686.1997
Daubin, V., and Szöllősi, G. J. (2016). Horizontal gene transfer and the history of life. Cold Spring Harb. Perspect. Biol. 8:a18036. doi: 10.1101/cshperspect.a018036
de Lorenzo, V. (2022). Environmental Galenics: large-scale fortification of extant microbiomes with engineered bioremediation agents. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 377:20210395. doi: 10.1098/rstb.2021.0395
Delsuc, F., Brinkmann, H., and Philippe, H. (2005). Phylogenomics and the reconstruction of the tree of life. Nat. Rev. Genet. 6, 361–375. doi: 10.1038/nrg1603
Doolittle, W. F., and Bapteste, E. (2007). Pattern pluralism and the tree of life hypothesis. Proc. Natl. Acad. Sci. U. S. A. 104, 2043–2049. doi: 10.1073/pnas.0610699104
El, K. K., Pontarotti, P., Raoult, D., and Fournier, P. E. (2016). Origin and evolution of rickettsial plasmids. PLoS One 11:e147492. doi: 10.1371/journal.pone.0147492
Espinosa, A., and Paz-Y-Miño-C, G. (2012). Discrimination, crypticity, and incipient taxa in entamoeba. J. Eukaryot. Microbiol. 59, 105–110. doi: 10.1111/j.1550-7408.2011.00606.x
Fan, Y., Xiao, Y., Momeni, B., and Liu, Y. (2018). Horizontal gene transfer can help maintain the equilibrium of microbial communities. J. Theor. Biol. 454, 53–59. doi: 10.1016/j.jtbi.2018.05.036
Feng, L., Wang, W., Cheng, J., Ren, Y., Zhao, G., Gao, C., et al. (2007). Genome and proteome of long-chain alkane degrading Geobacillus thermodenitrificans NG80-2 isolated from a deep-subsurface oil reservoir. Proc. Natl. Acad. Sci. U. S. A. 104, 5602–5607. doi: 10.1073/pnas.0609650104
Fierer, N. (2017). Embracing the unknown: disentangling the complexities of the soil microbiome. Nat. Rev. Microbiol. 15, 579–590. doi: 10.1038/nrmicro.2017.87
Fredrickson, J. K., Hicks, R. J., Li, S. W., and Brockman, F. J. (1988). Plasmid incidence in bacteria from deep subsurface sediments. Appl. Environ. Microbiol. 54, 2916–2923. doi: 10.1128/aem.54.12.2916-2923.1988
Frindte, K., Pape, R., Werner, K., Löffler, J., and Knief, C. (2019). Temperature and soil moisture control microbial community composition in an arctic–alpine ecosystem along elevational and micro-topographic gradients. ISME J. 13, 2031–2043. doi: 10.1038/s41396-019-0409-9
Galagan, J. E., Nusbaum, C., Roy, A., Endrizzi, M. G., Macdonald, P., FitzHugh, W., et al. (2002). The genome of M. acetivorans reveals extensive metabolic and physiological diversity. Genome Res. 12, 532–542. doi: 10.1101/gr.223902
Gallo, I., Furlan, J., Sanchez, D. G., and Stehling, E. G. (2019). Heavy metal resistance genes and plasmid-mediated quinolone resistance genes in Arthrobacter sp. isolated from Brazilian soils. Antonie Van Leeuwenhoek 112, 1553–1558. doi: 10.1007/s10482-019-01281-9
Gallos, L. K., Song, C., Havlin, S., and Makse, H. A. (2007). Scaling theory of transport in complex biological networks. Proc. Natl. Acad. Sci. 104, 7746–7751. doi: 10.1073/pnas.0700250104
Garrison-Schilling, K. L., Kaluskar, Z. M., Lambert, B., and Pettis, G. S. (2014). Genetic analysis and prevalence studies of the brp exopolysaccharide locus of Vibrio vulnificus. PLoS One 9:e100890. doi: 10.1371/journal.pone.0100890
Garushyants, S. K., Kazanov, M. D., and Gelfand, M. S. (2015). Horizontal gene transfer and genome evolution in Methanosarcina. BMC Evol. Biol. 15:102. doi: 10.1186/s12862-015-0393-2
Gogarten, J. P., Doolittle, W. F., and Lawrence, J. G. (2002). Prokaryotic evolution in light of gene transfer. Mol. Biol. Evol. 19, 2226–2238. doi: 10.1093/oxfordjournals.molbev.a004046
Gotz, A., and Smalla, K. (1997). Manure enhances plasmid mobilization and survival of Pseudomonas putida introduced into field soil. Appl. Environ. Microbiol. 63, 1980–1986. doi: 10.1128/aem.63.5.1980-1986.1997
Graham, J. B., and Istock, C. A. (1978). Genetic exchange in Bacillus subtilis in soil. Mol. Gen. Genet. MGG 166, 287–290. doi: 10.1007/BF00267620
Graham, D. W., Knapp, C. W., Christensen, B. T., McCluskey, S., and Dolfing, J. (2016). Appearance of β-lactam resistance genes in agricultural soils and clinical isolates over the 20th century. Sci. Rep. 6:21550. doi: 10.1038/srep21550
Guimerà, R., and Nunes Amaral, L. A. (2005). Functional cartography of complex metabolic networks. Nature 433, 895–900. doi: 10.1038/nature03288
Guo, X., Feng, J., Shi, Z., Zhou, X., Yuan, M., Tao, X., et al. (2018). Climate warming leads to divergent succession of grassland microbial communities. Nat. Clim. Chang. 8, 813–818. doi: 10.1038/s41558-018-0254-2
Harholt, J., Sørensen, I., Fangel, J., Roberts, A., Willats, W. G., Scheller, H. V., et al. (2012). The glycosyltransferase repertoire of the spikemoss Selaginella moellendorffii and a comparative study of its cell wall. PLoS One 7:e35846. doi: 10.1371/journal.pone.0035846
Hernández-López, A., Chabrol, O., Royer-Carenzi, M., Merhej, V., Pontarotti, P., and Raoult, D. (2013). To tree or not to tree? Genome-wide quantification of recombination and reticulate evolution during the diversification of strict intracellular bacteria. Genome Biol. Evol. 5, 2305–2317. doi: 10.1093/gbe/evt178
Huang, J. (2013). Horizontal gene transfer in eukaryotes: the weak-link model. BioEssays 35, 868–875. doi: 10.1002/bies.201300007
Huerta-Cepas, J., Szklarczyk, D., Heller, D., Hernández-Plaza, A., Forslund, S. K., Cook, H., et al. (2019). EggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses. Nucleic Acids Res. 47, D309–D314. doi: 10.1093/nar/gky1085
Hupert-Kocurek, K., Guzik, U., and Wojcieszyńska, D. (2012). Characterization of catechol 2,3-dioxygenase from Planococcus sp. strain S5 induced by high phenol concentration. Acta Biochim. Pol. 59, 345–351. doi: 10.18388/abp.2012_2119
Huson, D. H., and Bryant, D. (2006). Application of phylogenetic networks in evolutionary studies. Mol. Biol. Evol. 23, 254–267. doi: 10.1093/molbev/msj030
Iyer, L. M., Zhang, D., Burroughs, A. M., and Aravind, L. (2013). Computational identification of novel biochemical systems involved in oxidation, glycosylation and other complex modifications of bases in DNA. Nucleic Acids Res. 41, 7635–7655. doi: 10.1093/nar/gkt573
Jain, R., Rivera, M. C., and Lake, J. A. (1999). Horizontal gene transfer among genomes: the complexity hypothesis. Proc. Natl. Acad. Sci. U. S. A. 96, 3801–3806. doi: 10.1073/pnas.96.7.3801
Jassey, V. E. J., Walcker, R., Kardol, P., Geisen, S., Heger, T., Lamentowicz, M., et al. (2022). Contribution of soil algae to the global carbon cycle. New Phytol. 234, 64–76. doi: 10.1111/nph.17950
Jeong, H., Mason, S. P., Barabási, A. L., and Oltvai, Z. N. (2001). Lethality and centrality in protein networks. Nature 411, 41–42. doi: 10.1038/35075138
Jiang, X., Ellabaan, M., Charusanti, P., Munck, C., Blin, K., Tong, Y., et al. (2017). Dissemination of antibiotic resistance genes from antibiotic producers to pathogens. Nat. Commun. 8:15784. doi: 10.1038/ncomms15784
Johannes Sørensen, S., Schyberg, T., and Rønn, R. (1999). Predation by protozoa on Escherichia coli K12 in soil and transfer of resistance plasmid RP4 to indigenous bacteria in soil. Appl. Soil Ecol. 11, 79–90. doi: 10.1016/S0929-1393(98)00117-6
Kintz, E., Heiss, C., Black, I., Donohue, N., Brown, N., Davies, M. R., et al. (2017). Salmonella enterica Serovar Typhi lipopolysaccharide O-antigen modification impact on serum resistance and antibody recognition. Infect. Immun. 85. doi: 10.1128/IAI.01021-16
Kloesges, T., Popa, O., Martin, W., and Dagan, T. (2011). Networks of gene sharing among 329 proteobacterial genomes reveal differences in lateral gene transfer frequency at different phylogenetic depths. Mol. Biol. Evol. 28, 1057–1074. doi: 10.1093/molbev/msq297
Klümper, U., Riber, L., Dechesne, A., Sannazzarro, A., Hansen, L. H., Sørensen, S. J., et al. (2015). Broad host range plasmids can invade an unexpectedly diverse fraction of a soil bacterial community. ISME J. 9, 934–945. doi: 10.1038/ismej.2014.191
Koonin, E. V. (2005). Orthologs, paralogs, and evolutionary genomics. Annu. Rev. Genet. 39, 309–338. doi: 10.1146/annurev.genet.39.073003.114725
Kunin, V., Goldovsky, L., Darzentas, N., and Ouzounis, C. A. (2005). The net of life: reconstructing the microbial phylogenetic network. Genome Res. 15, 954–959. doi: 10.1101/gr.3666505
Lacroix, B., and Citovsky, V. (2016). Transfer of DNA from Bacteria to eukaryotes. MBio 7. doi: 10.1128/mBio.00863-16
Lang, A. S., and Beatty, J. T. (2007). Importance of widespread gene transfer agent genes in α-proteobacteria. Trends Microbiol. 15, 54–62. doi: 10.1016/j.tim.2006.12.001
Lawrence, J. G., and Ochman, H. (2002). Reconciling the many faces of lateral gene transfer. Trends Microbiol. 10, 1–4. doi: 10.1016/s0966-842x(01)02282-x
Li, L., Liu, Z., Meng, D., Liu, X., Li, X., Zhang, M., et al. (2019). Comparative genomic analysis reveals the distribution, organization, and evolution of metal resistance genes in the genus acidithiobacillus. Appl. Environ. Microbiol. 85, e2118–e2153. doi: 10.1128/AEM.02153-18
Li, L., Liu, Z., Zhang, M., Meng, D., Liu, X., Wang, P., et al. (2020). Insights into the metabolism and evolution of the genus Acidiphilium, a typical acidophile in acid mine drainage. MSystems 5, e820–e867. doi: 10.1128/mSystems.00867-20
Li, L., Liu, Z., Zhou, Z., Zhang, M., Meng, D., Liu, X., et al. (2021). Comparative genomics provides insights into the genetic diversity and evolution of the DPANN superphylum. MSystems 6, e0060221–e0060621. doi: 10.1128/mSystems.00602-21
Li, L., Peng, S., Wang, Z., Zhang, T., Li, H., Xiao, Y., et al. (2022). Genome mining reveals abiotic stress resistance genes in plant genomes acquired from microbes via HGT. Front. Plant Sci. 13:1025122. doi: 10.3389/fpls.2022.1025122
Li, X., Wang, H., Tong, W., Feng, L., Wang, L., Rahman, S. U., et al. (2020). Exploring the evolutionary dynamics of Rhizobium plasmids through bipartite network analysis. Environ. Microbiol. 22, 934–951. doi: 10.1111/1462-2920.14762
Lin, S., Hanson, R. E., and Cronan, J. E. (2010). Biotin synthesis begins by hijacking the fatty acid synthetic pathway. Nat. Chem. Biol. 6, 682–688. doi: 10.1038/nchembio.420
Lozupone, C. A., Hamady, M., Cantarel, B. L., Coutinho, P. M., Henrissat, B., Gordon, J. I., et al. (2008). The convergence of carbohydrate active gene repertoires in human gut microbes. Proc. Natl. Acad. Sci. U. S. A. 105, 15076–15081. doi: 10.1073/pnas.0807339105
Maatouk, M., Ibrahim, A., Pinault, L., Armstrong, N., Azza, S., Rolain, J. M., et al. (2022). New beta-lactamases in candidate phyla radiation: owning pleiotropic enzymes is a smart paradigm for microorganisms with a reduced genome. Int. J. Mol. Sci. 23. doi: 10.3390/ijms23105446
Markowitz, V. M., Chen, I. A., Palaniappan, K., Chu, K., Szeto, E., Grechkin, Y., et al. (2010). The integrated microbial genomes system: an expanding comparative analysis resource. Nucleic Acids Res. 38, D382–D390. doi: 10.1093/nar/gkp887
McCarthy, C. G., and Fitzpatrick, D. A. (2016). Systematic search for evidence of interdomain horizontal gene transfer from prokaryotes to oomycete lineages. MSphere 1. doi: 10.1128/mSphere.00195-16
Mølbak, L., Molin, S., and Kroer, N. (2007). Root growth and exudate production define the frequency of horizontal plasmid transfer in the Rhizosphere. FEMS Microbiol. Ecol. 59, 167–176. doi: 10.1111/j.1574-6941.2006.00229.x
Muller, E. E. L., Faust, K., Widder, S., Herold, M., Martínez Arbas, S., and Wilmes, P. (2018). Using metabolic networks to resolve ecological properties of microbiomes. Curr. Opin. Syst. Biol. 8, 73–80. doi: 10.1016/j.coisb.2017.12.004
Musovic, S., Oregaard, G., Kroer, N., and Sørensen, S. J. (2006). Cultivation-independent examination of horizontal transfer and host range of an IncP-1 plasmid among gram-positive and gram-negative bacteria indigenous to the barley rhizosphere. Appl. Environ. Microbiol. 72, 6687–6692. doi: 10.1128/AEM.00013-06
Nelson-Sathi, S., Dagan, T., Landan, G., Janssen, A., Steel, M., McInerney, J. O., et al. (2012). Acquisition of 1,000 eubacterial genes physiologically transformed a methanogen at the origin of Haloarchaea. Proc. Natl. Acad. Sci. U. S. A. 109, 20537–20542. doi: 10.1073/pnas.1209119109
Nojiri, H., Shintani, M., and Omori, T. (2004). Divergence of mobile genetic elements involved in the distribution of xenobiotic-catabolic capacity. Appl. Microbiol. Biotechnol. 64, 154–174. doi: 10.1007/s00253-003-1509-y
Norman, A., Hansen, L. H., and Sørensen, S. J. (2009). Conjugative plasmids: vessels of the communal gene pool. Philos. Trans. R. Soc. B Biol. Sci. 364, 2275–2289. doi: 10.1098/rstb.2009.0037
Ogunseitan, O. A., Tedford, E. T., Pacia, D., Sirotkin, K. M., and Sayler, G. S. (1987). Distribution of plasmids in groundwater bacteria. J. Ind. Microbiol. 1, 311–317. doi: 10.1007/BF01569309
Pál, C., Papp, B., and Lercher, M. J. (2005). Adaptive evolution of bacterial metabolic networks by horizontal gene transfer. Nat. Genet. 37, 1372–1375. doi: 10.1038/ng1686
Paquola, A. C. M., Asif, H., Pereira, C. A. D. B., Feltes, B. C., Bonatto, D., Lima, W. C., et al. (2018). Horizontal gene transfer building prokaryote genomes: genes related to exchange between cell and environment are frequently transferred. J. Mol. Evol. 86, 190–203. doi: 10.1007/s00239-018-9836-x
Popa, O., Hazkani-Covo, E., Landan, G., Martin, W., and Dagan, T. (2011). Directed networks reveal genomic barriers and DNA repair bypasses to lateral gene transfer among prokaryotes. Genome Res. 21, 599–609. doi: 10.1101/gr.115592.110
Proulx, S. R., Promislow, D. E. L., and Phillips, P. C. (2005). Network thinking in ecology and evolution. Trends Ecol. Evol. 20, 345–353. doi: 10.1016/j.tree.2005.04.004
Rezende, E. L., Lavabre, J. E., Guimarães, P. R., Jordano, P., and Bascompte, J. (2007). Non-random coextinctions in phylogenetically structured mutualistic networks. Nature 448, 925–928. doi: 10.1038/nature05956
Richaume, A., Angle, J. S., and Sadowsky, M. J. (1989). Influence of soil variables on in situ plasmid transfer from Escherichia coli to Rhizobium fredii. Appl. Environ. Microbiol. 55, 1730–1734. doi: 10.1128/aem.55.7.1730-1734.1989
Richaume, A., Smit, E., Faurie, G., and van Elsas, J. D. (1992). Influence of soil type on the transfer of plasmid RP4p from Pseudomonas fluorescens to introduced recipient and to indigenous bacteria. FEMS Microbiol. Lett. 101, 281–291. doi: 10.1016/0378-1097(92)90825-9
Rochelle, P. A., Fry, J. C., and Day, M. J. (1989). Factors affecting conjugal transfer of plasmids encoding mercury resistance from pure cultures and mixed natural suspensions of epilithic bacteria. J. Gen. Microbiol. 135, 409–424. doi: 10.1099/00221287-135-2-409
Romine, M. F., Stillwell, L. C., Wong, K. K., Thurston, S. J., Sisk, E. C., Sensen, C., et al. (1999). Complete sequence of a 184-kilobase catabolic plasmid from Sphingomonas aromaticivorans F199. J. Bacteriol. 181, 1585–1602. doi: 10.1128/JB.181.5.1585-1602.1999
Sengeløv, G., Kowalchuk, G. A., and Sørensen, S. J. (2000). Influence of fungal-bacterial interactions on bacterial conjugation in the residuesphere. FEMS Microbiol. Ecol. 31, 39–45. doi: 10.1016/s0168-6496(99)00079-3
Ševcíková, T., Yurchenko, T., Fawley, K. P., Amaral, R., Strnad, H., Santos, L., et al. (2019). Plastid genomes and proteins illuminate the evolution of eustigmatophyte algae and their bacterial endosymbionts. Genome Biol. Evol. 11, 362–379. doi: 10.1093/gbe/evz004
Song, H., Ding, M., Jia, X., Ma, Q., and Yuan, Y. (2014). Synthetic microbial consortia: from systematic analysis to construction and applications. Chem. Soc. Rev. 43, 6954–6981. doi: 10.1039/C4CS00114A
Tamminen, M., Virta, M., Fani, R., and Fondi, M. (2012). Large-scale analysis of plasmid relationships through gene-sharing networks. Mol. Biol. Evol. 29, 1225–1240. doi: 10.1093/molbev/msr292
Thomas, C. M., and Nielsen, K. M. (2005). Mechanisms of, and barriers to, horizontal gene transfer between bacteria. Nat. Rev. Microbiol. 3, 711–721. doi: 10.1038/nrmicro1234
Top, E. M., Springael, D., and Boon, N. (2002). Catabolic mobile genetic elements and their potential use in bioaugmentation of polluted soils and waters. FEMS Microbiol. Ecol. 42, 199–208. doi: 10.1111/j.1574-6941.2002.tb01009.x
Trevors, J. T., van Elsas, J. D., van Overbeek, L. S., and Starodub, M. E. (1990). Transport of a genetically engineered Pseudomonas fluorescens strain through a soil microcosm. Appl. Environ. Microbiol. 56, 401–408. doi: 10.1128/aem.56.2.401-408.1990
Trias, R., Ménez, B., le Campion, P., Zivanovic, Y., Lecourt, L., Lecoeuvre, A., et al. (2017). High reactivity of deep biota under anthropogenic CO(2) injection into basalt. Nat. Commun. 8:1063. doi: 10.1038/s41467-017-01288-8
Tsigos, I., Velonia, K., Smonou, I., and Bouriotis, V. (1998). Purification and characterization of an alcohol dehydrogenase from the Antarctic psychrophile Moraxella sp. TAE123. Eur. J. Biochem. 254, 356–362. doi: 10.1046/j.1432-1327.1998.2540356.x
Tsoi, R., Dai, Z., and You, L. (2019). Emerging strategies for engineering microbial communities. Biotechnol. Adv. 37:107372. doi: 10.1016/j.biotechadv.2019.03.011
van Elsas, J. D., Turner, S., and Trevors, J. T. (2006). “Bacterial conjugation in soil” in eds. P. Nannipieri and K. Smalla (Berlin, Heidelberg: Springer Berlin Heidelberg), 331–353.
Vieira, S., Sikorski, J., Dietz, S., Herz, K., Schrumpf, M., Bruelheide, H., et al. (2020). Drivers of the composition of active rhizosphere bacterial communities in temperate grasslands. ISME J. 14, 463–475. doi: 10.1038/s41396-019-0543-4
Wang, B., Pandey, T., Long, Y., Delgado-Rodriguez, S. E., Daugherty, M. D., and Ma, D. K. (2022). Co-opted genes of algal origin protect C. elegans against cyanogenic toxins. Curr. Biol. 32, 4941–4948.e3. doi: 10.1016/j.cub.2022.09.041
Weinberg, S. R., and Stotzky, G. (1972). Conjugation and genetic recombination of Escherichia coli in soil. Soil Biol. Biochem. 4, 171–180. doi: 10.1016/0038-0717(72)90008-9
Yu, Y., Ma, J., Guo, Q., Ma, J., and Wang, H. (2019). A novel 3-oxoacyl-ACP reductase (FabG3) is involved in the xanthomonadin biosynthesis of Xanthomonas campestris pv. Campestris. Mol Plant Pathol 20, 1696–1709. doi: 10.1111/mpp.12871
Zhang, Y., Guo, Y., Qiu, T., Gao, M., and Wang, X. (2022). Bacteriophages: underestimated vehicles of antibiotic resistance genes in the soil. Front. Microbiol. 13:936267. doi: 10.3389/fmicb.2022.936267
Zhang, M., Warmink, J., Pereira, E. S. M., Brons, J., Smalla, K., and van Elsas, J. D. (2015). IncP-1β plasmids are important carriers of fitness traits for Variovorax species in the Mycosphere—two novel plasmids, pHB44 and pBS64, with differential effects unveiled. Microb. Ecol. 70, 141–153. doi: 10.1007/s00248-014-0550-y
Zhang, P., Zhang, Z., Zhang, L., Wang, J., and Wu, C. (2020). Glycosyltransferase GT1 family: phylogenetic distribution, substrates coverage, and representative structural features. Comput. Struct. Biotechnol. J. 18, 1383–1390. doi: 10.1016/j.csbj.2020.06.003
Zheng, Q., Hu, Y., Zhang, S., Noll, L., Böckle, T., Dietrich, M., et al. (2019). Soil multifunctionality is affected by the soil environment and by microbial community composition and diversity. Soil Biol. Biochem. 136:107521. doi: 10.1016/j.soilbio.2019.107521
Zhou, H., Beltrán, J. F., and Brito, I. L. (2021). Functions predict horizontal gene transfer and the emergence of antibiotic resistance. Sci. Adv. 7:j5056:eabj5056. doi: 10.1126/sciadv.abj5056
Zhou, J., Deng, Y., Luo, F., He, Z., Tu, Q., and Zhi, X. (2010). Functional molecular ecological networks. MBio 1, e110–e169. doi: 10.1128/mBio.00169-10
Keywords: soil, microbiota, horizontal gene transfer, network, carbon metabolic gene
Citation: Li L, Liu Y, Xiao Q, Xiao Z, Meng D, Yang Z, Deng W, Yin H and Liu Z (2023) Dissecting the HGT network of carbon metabolic genes in soil-borne microbiota. Front. Microbiol. 14:1173748. doi: 10.3389/fmicb.2023.1173748
Edited by:
Ye Deng, Chinese Academy of Sciences (CAS), ChinaReviewed by:
Zheng Gao, Shandong Agricultural University, ChinaZiyan Wei, Chinese Academy of Sciences (CAS), China
Copyright © 2023 Li, Liu, Xiao, Xiao, Meng, Yang, Deng, Yin and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yongjun Liu, VmluY2VudGxpdTIwMjBAMTYzLmNvbQ==; Zhenghua Liu, bGl1emhlbmdodWEyMDE3Y3N1QDE2My5jb20=; Wenqiao Deng, MjM3MzgzMjBAcXEuY29t