- Research article
- Open Access
Lateral gene transfer of an ABC transporter complex between major constituents of the human gut microbiome
BMC Microbiologyvolume 12, Article number: 248 (2012)
Several links have been established between the human gut microbiome and conditions such as obesity and inflammatory bowel syndrome. This highlights the importance of understanding what properties of the gut microbiome can affect the health of the human host. Studies have been undertaken to determine the species composition of this microbiome and infer functional profiles associated with such host properties. However, lateral gene transfer (LGT) between community members may result in misleading taxonomic attributions for the recipient organisms, thus making species-function links difficult to establish.
We identified a peptides/nickel transport complex whose components differed in abundance based upon levels of host obesity, and assigned the encoded proteins to members of the microbial community. Each protein was assigned to several distinct taxonomic groups, with moderate levels of agreement observed among different proteins in the complex. Phylogenetic trees of these proteins produced clusters that differed greatly from taxonomic attributions and indicated that habitat-directed LGT of this complex is likely to have occurred, though not always between the same partners.
These findings demonstrate that certain membrane transport systems may be an important factor within an obese-associated gut microbiome and that such complexes may be acquired several times by different strains of the same species. Additionally, an example of individual proteins from different organisms being transferred into one operon was observed, potentially demonstrating a functional complex despite the donors of the subunits being taxonomically disparate. Our results also highlight the potential impact of habitat-directed LGT on the resident microbiota.
A vast array of bacteria, archaea, viruses and eukaryotes inhabit the tract of the human gut and form its microbiome [1, 2]. Investigation into the composition of this densely packed community and its effect on the host have revealed several benefits derived from the microorganisms such as plant polysaccharide processing and amino acid synthesis [1, 3]. The species structure of the community has also been linked to several health problems such as inflammatory bowel disease  and obesity [5–7].
Initial studies of the human gut microbiome involved sequencing of the 16S ribosomal RNA gene to determine the main constituents of the community. Although many organisms observed in these studies were previously uncharacterised , members of the phyla Firmicutes and Bacteroidetes comprised over 90% of the population of known bacterial species within the gut . The Human Microbiome Project (HMP) utilised both a 16S-based approach and a large-scale study of obese and lean twin pairs, and found that the species composition of the gut microbiome was more similar in related individuals than unrelated individuals . However no core species group was observed in all studied individuals. A preliminary investigation of full genome sequences was also performed on a subset of samples in this study, revealing that similar taxonomic profiles were linked to similar metabolic profiles between individuals . Each of the two main phyla (Firmicutes and Bacteroidetes) was associated with enrichment of different metabolic pathways (transporters and carbohydrate metabolism respectively) and although the species composition differed between individuals, there was a relatively high level of functional conservation in the majority of gut microbiomes studied.
Associative studies between the human gut microbiome and host factors such as inflammatory bowel disease (IBD) and weight have revealed close ties between the composition of the microorganism community and human health [4, 6, 9, 10]. Metagenomic sequencing of faecal samples from 124 European individuals was performed in order to study multiple portions of the community gene pool and link variation in community to IBD . A core gut microbiome gene pool was reported along with a proposed list of possible core species. These species were primarily from the two main phyla identified previously, and taxonomic rank abundances were used to distinguish between IBD and non-IBD individuals. Taxonomic differences have also been linked to obesity, especially based upon relative abundances of different phyla. Turnbaugh et al. found that obese twins had a lower proportion of Bacteroidetes than lean twins . This relationship between weight and the reduction of Bacteroidetes species has also been supported by other studies [5, 10]. However, additional studies have found either no significant change in the Firmicutes: Bacteroidetes ratio [6, 11] or even an increase in Bacteroidetes in obese individuals .
The aim of our study was to investigate whether links could be made between an individual’s body mass index (BMI) and metabolic functions within the microbiome. Findings indicate that multiple components of the peptides/nickel transport system show consistent differences in abundance based upon levels of obesity within the sampled individuals. This transporter is comprised of five proteins and is likely used to transport nickel into cells and regulate its intracellular concentration , or potentially regulate the expression of cell surface molecules through selective uptake of short peptides . Despite significant differences in the abundance of complex members, the taxonomic distribution of these proteins did not differ between obese and lean individuals. However, phylogenetic analysis of abundant species, regardless of BMI, revealed that these proteins were likely laterally acquired from other gut-associated microbes, indicating that habitat-directed LGT can influence microbial metabolic systems that are linked to human health.
Results and discussion
Prediction of open reading frames (ORFs) from the dataset of 124 patients presented in  revealed an average of 203,300 potential ORFs per sample. Use of BLAST sequence matching resulted in predicted protein functions for, on average, 46% of the ORFs per sample. Subsequent characterisation of these putative protein sequence fragments using the KEGG database allowed for metabolic classification of 39% of the ORFs with BLAST hits (18% of the original predicted ORF set). Each microbiome sample had an average of 2,400 KO groupings containing at least one sequence fragment with a total of 4,849 KOs being present in at least one sample in the dataset.
Distributions of predicted metabolic functions between low and high-BMI groups
Sequence counts for all 4,849 KOs were compared across patients in order to identify metabolic functions that differ in abundance between low BMI (18 to 22) and high BMI (30+) associated samples. Present KEGG Orthology groups ranged in relative abundance from 4 × 10-5 (i.e. one copy of the protein in the largest sample) to 0.8% of the total assigned proteins, with K06147 (bacterial ATP-binding cassette, subfamily B) as the most abundant KO across all patients, regardless of BMI. Fifty-two KOs were found to differ significantly (Bonferroni-corrected p value <0.01) in abundance levels between lean- and obese-related samples. The majority of these KOs were low in frequency in both BMI categories; apart from the ABC transporter mentioned above, only five of the 52 KOs had a mean proportion in both BMI sets of 0.2% or higher (Figure 1). K06147, in addition to being the most abundant protein in all patients, was 46% more abundant in low-BMI samples. The other four KOs that were found to have significant differences in abundances all belong to the peptides/nickel transport system module (KEGG module M00239). This module contains five ABC transporter proteins (K02031-K02035), four of which were found to be significantly more abundant in low-BMI patients (K02031-K02034; ratios ranging between 42 and 44%; corrected p-values < 0.01) (Figure 1). This transport system contains two ATP-binding proteins (K02031 and K02032), two permeases (K02033 and K02034) and one substrate-binding protein (K02035). Variation in abundances of each KO between patients in the same BMI group (lean or obese) was found to be low, with mean proportions at most 0.2%. Although differences in abundance of K02035 were not found to be as statistically supported as the other subunits (p-value 0.021) it was found at similar levels of abundance between patients as the other four members of the transport system. Thus K02035 was included alongside the other subunits in the module in order to identify if specific species are associated with the complex as a whole.
Taxonomic assignment of metagenomic fragments associated with nickel transporters
Reference phylogenetic trees were constructed for each of the five KOs within the peptides/nickel transport complex using proteins from 3,181 sequenced genomes retrieved from IMG  (Additional file 1: Figure S1). Habitat metadata from the IMG database  was used to assign species to the human gastrointestinal tract resulting in 472 gut-associated species. It was found that these species were spread throughout the trees and did not appear to cluster based upon habitat (Additional file 1: Figure S1). We constructed subtrees containing only gut-associated species and assessed the cohesion of taxonomic groups using the consistency index (CI): CIs close to 1.0 indicate perfect clustering of all taxonomic groups at a particular rank, while low CIs indicate intermingling of organisms from different groups and are suggestive of LGT, especially if organisms in the same cluster are from very disparate groups. The CIs of all trees were less than 0.5 when evaluated at the ranks of family, class, order and phylum (Additional file 2: Table S1), suggesting a lack of cohesion of major lineages. CIs at the genus (0.60 to 0.64) and species (0.93 to 0.96) levels were higher, indicating less disruption of these groups. Examples of disrupted species include Faecalibacterium prausnitzii and Clostridium difficile in the tree of K02031 sequences from gut-associated species (Additional file 3: Figure S2); in this case, large evolutionary distances separated sequences associated with strains of the same species. However as such disparities were also observed within the trees containing all species, not just gut-associated strains, further analysis was required to discover whether LGT events were directed by environment.
Pplacer  was used to place metagenomic fragments onto expanded reference trees for each of the KOs of interest. Not all fragments were mapped down to species level and thus a proportion was assigned only to a rank of genus or higher. The quantity of reads that were unclassified at different levels due either to lack of placement confidence of the read below a certain taxonomic level or lack of NCBI taxonomy information varied between KOs (Table 1). Taxonomic assignment was above 75% at all levels of classification with an average of 93% per rank. Fragments that were not mapped below a certain level were labelled as ‘unclassified’ and disregarded in further abundance analysis at that level. In general, Firmicutes were the dominant phylum associated with each KO, as is to be expected by their abundance within the gut , with the class Clostridia and order Clostridiales making up the largest proportion of classified reads in each sample. Several Firmicute genera, including Clostridium, Blautia, Ruminococcus and Faecalibacterium, were found to be in relatively high abundance in almost every protein set (up to 15%). Members of other phyla such as Proteobacteria and Actinobacteria also contributed to the species composition of proteins within this complex though these signals were less abundant and consistent than the Firmicute members. Thus, although correlation of assignments at higher taxonomic ranks was found between KOs, this did not extend to the genus level. This could be due to incorrect taxonomic assignments as a result of a deficiency in relevant reference genomes or lack of predictive power from the metagenomic ORFs. Inconsistencies could also be due to recent LGT events between members of different genera, which would result in discordant taxonomic assignments associated with the recipient species. Thus it is possible that this protein complex is present in a smaller, more consistent, set of genera with the human gut microbiome than is observed here.
Mapping of species classifications revealed further disparate signals between the KOs. Within each of the proteins K02031-K02035, no single species was represented in more than 9% of taxonomic attributions (Table 2). Collectively, the top four contributing species did not comprise more than 25% of the taxonomic groups associated with any of these KOs. As many of the fragments were not classified to the species level (average of 17.12%), it is difficult to determine exactly what species are most commonly associated with each protein. Analysis of the peptides/nickel transport system revealed very little overlap in species composition between the individual proteins of the complex. Only Faecalibacterium prausnitzii was found in relatively high abundance in all five KO phylogenies, with most other highly abundant species only being highly associated with at most three components. However, all of the most abundantly associated species are resident within either the gut or the oral cavity of the human microbiome. Thus, despite low overlap of species composition, fragments were found to be derived from microbes associated with the human alimentary canal as is to be expected.
Analysis of Faecalibacterium prausnitziistrains within reference protein phylogenetic trees
The probable origin of each subunit of the peptides/nickel transport system within F. prausnitzii was examined using full-length protein trees derived from 3,181 sequenced species. It was found that the five sequenced strains of this species (M21/2, A2-165, KLE1255, SL3/3 and L2-6) contained up to 6 copies of each gene, which were spread across up to six operons with an average of 2.8 per strain (Figure 2). Operons were classified based upon whether the strains formed a closely related group within the full protein tree of the constituent KOs. Up to six such groups were found within each protein tree for K02031-K02035, resulting in the postulation of six operon types, each with a potential separate origin. Each operon type appeared to be derived from an LGT event from strains of various taxonomically spread species (Additional file 4: Figure S3). However, most of these species are associated with the human gut microbiome, suggesting that habitat-direct LGT occurred. Operon 3, which is complete only in strain A2-165, appears to have been potentially acquired from multiple bacterial species with the ATP-binding proteins (K02031 and K02032) separately acquired from the remaining proteins (Additional file 4: Figure S3). Gene neighbourhood analysis revealed preservation of operon organisation between F. prausnitzii strains and potential donors of operons, though not similarity in flanking regions, adding credence to the possibility of LGT resulting in acquisition of this function. Although multiple strains of F. prausnitzii contain each type of operon, suggesting acquisition prior to strain separation, rearrangement of the gene constituents appears to be frequent with inversions observed in operon types 2 and 5 and potential loss of components in operons 3, 4, 5 and 6 (although sequence similarity between missing sections of operon 5 in strains A2-165 and L2-6 and K02035 indicate this gene is present, though not annotated correctly).
Although high abundance of F. prausnitzii was found in association with the peptides/nickel transport complex, regardless of BMI, analysis of the species abundance associated with changes in BMI revealed no noticeable difference between low and high BMI patients. This could be due to the high numbers of unclassified reads, several cases of LGT confusing the species abundance signals or the difference in gene copy numbers between strains of F. prausnitzii.
The investigation into function-species relationships undertaken here highlights some important aspects of microbiome studies and the possible inferences that can be made from such information. Although there are potential pitfalls with analysis of abundance of functions within a microbiome as has been done here such as insufficient sampling depth or incomplete sequencing of all DNA fragments, such approaches have revealed marked differences previously [5, 17]. It was found that the abundance of components of the peptides/nickel transport system differed between low and high BMI related samples, likely indicating a link between this system and obesity although such a correlation would require validation on other datasets. Taxonomic assignment of KO-associated reads showed that within the peptides/nickel transport system, there are multiple species associated with each KO, with dominance by one species being rare (Table 2). There are numerous possible reasons for this inconsistency of dominant species between KOs. As it is highly implausible that each protein is being created by different species and somehow incorporated separately into the transport systems, it is more likely LGT has resulted in operon or single gene transfers between organisms. This would result in conflicting phylogenetic relationships as observed here and makes determination of the true species involved in pathways difficult. This situation is likely due to the high degree of LGT known to occur in the human gut [18–20]. Although the presence of F. prausnitzii in all five KO sets may indicate that this species is one of the dominant organisms involved in this pathway, such extrapolation cannot be confirmed without knowing the exact history of LGT events within the microbiome, or with much deeper sequencing that allows for assembly of large genomic fragments as was performed in . Therefore further insight into detecting lateral gene transfer within the microbiome and determining the true species involved in each pathway is required before accurate relationships between species, functions and host properties such as disease be made with confidence.
Analysis of the peptides/nickel transport complex with F. prausnitzii revealed multiple operons associated with this function, each of which appeared to have been acquired through lateral gene transfer. Previous work on Fusobacterium nucleatum found an iron transport complex within the genome that resulted both from LGT of an entire operon and separate LGT events of single genes from multiple strains of other species resulting in two other operons of heterogeneous origins . Within F. prausnitzii it appears that a similar scenario has occurred within the peptides/nickel transporter with six operons types discovered. It was determined that each operon arose from separate LGT events through analysis of congruent gene trees within the operon (Additional file 4: Figure S3), which is a strong indicator of LGT [22, 23]. Five of the six operon types appear to be derived from the transfer of the whole operon into strains of F. prausnitzii, though the presence of the same operon type in some but not all strains suggests such transfers occurred prior to the divergence of certain strains. The remaining operon which was only found in a complete form within strain A2-165 appears to have been acquired from multiple sources, with the majority of the genes derived from Lachnospiraceae bacterium 3_1_57FAA_CT1 with the two ATP-binding related genes derived from other sources (Additional file 4: Figure S3). This may be due to a whole operon transfer followed by subsequent orthologous replacement and demonstrates that although the complexity hypothesis suggests such interactions between a new protein and the pre-existing complex would fail , heterogeneous integration can occur and may result in loss of fitness [25, 26], if this operon is active. Thus if multiple acquisitions did take place, this could point to a system of gradual gain of novel functions from multiple sources. However, functional assays (such as those performed in ) would be required to determine if this operon is transcribed and translated into a complex within this strain.
It may be that all five strains of F. prausnitzii acquired this transport system from independent sources within their environment (or across habitats from strains of closely related species) via gain-of-function LGT or already possessed the operon which was subsequently overwritten by multiple orthologous replacements, making the history of the lateral gene transfers difficult to trace. The relevance of nickel or short peptide transport within this species is difficult to interpret. Several enzymes such as ureases, hydrogenases, methane reductases and carbon monoxide dehydrogenases use nickel as a cofactor  though F. prausnitzii is not known to have urease activity or many hydrolases . However, a relationship between nickel concentration and butyrate production, a product of F. prausnitzii, has been postulated, and demonstrated in cattle . This could indicate that these strains are influencing the levels of butyrate within the surrounding environment. Concentrations of butyrate and butyrate-producing bacteria have been associated with lower carbohydrate intake  and also reduced obesity in mice . This suggests that a subset of the enzymatic functions associated with nickel , specifically links to butyrate production and may be connected to levels of obesity with the host, possibly through influence of butyrate production. Additionally, as this transport system can also be involved in more general transport of peptide from two to five amino acid residues in length it could be another unknown function being utilised by this species within the human digestive tract habitat. This module was characterised based upon the Opp complex in Salmonella typhimurium, which has been shown to be involved in modulating expression of surface-exposed proteins . These proteins may be involved in functions such as sporulation and virulence, both of which have been shown to be important in the human gut microbiome [19, 33]. Thus it is possible that this transporter is not involved in nickel regulation but actually modulating the cell surface responses to the digestive tract environment. As it has been shown that low levels of F. prausnitzii are associated with Crohn’s disease  and we have shown here that F. prausnitzii may also be associated with obesity, it is likely that LGT of systems such as peptides/nickel transport may contribute to host adaptation of this species, as has been observed with LGT in other species [35, 36], or play a role in determining the importance of the species within the microbiome. However, further experimental analysis would be required to confirm the link between this membrane transport system and host obesity and also determine is precise function.
Understanding the effect of habitat-directed LGT is a difficult problem. Microbiome data can be utilised to address this as has been shown here. We have found that although an overall signal for clustering of gut-associated organisms was not observed, this is not indicative of a lack of LGT. Each protein tree did not correlate exactly with a species tree as would be usually derived from single-gene studies based on 16S or other marker genes. Subsequent analysis revealed that some species that were clustered together in the protein trees were from taxonomically distant groups (Additional file 4: Figure S3). These species were usually found to be occupying similar environmental niches and were possibly associated with influencing the habitat, in this case the BMI of the host. Thus these findings signify that subsets of species may share genetic information within the environment and such LGT may impact how the habitat as a whole is shaped.
The dataset of  derived from 124 European individuals using Illumina sequencing was used for this analysis. Deep sequencing of samples from these individuals resulted in an average of 4.5 Gb of data per patient, which was further assembled into contigs as described in reference . Associated with these sequences is a range of metadata including BMI, an indicator of the level of obesity of the patient. Low BMI (18 to 22) indicates underweight/healthy patients and a BMI of 30 and above indicates an obese individual. Only lean (low BMI; 34 samples) and obese (high BMI; 33 samples) patients were selected for further analysis to maximise any differences in the microbiome that may be associated with weight.
Functional assignment of proteins and estimation of abundances within the microbiome metabolic profile
Assembled contigs from each patient were used as input into Orphelia  for prediction of open reading frames (ORFs). Any predicted ORFs of length < 150 nucleotides were removed to ensure greater coverage for prediction of function. Prediction of protein function for each ORF was undertaken using UBLAST as implemented in USEARCH version 4.0.38  against a protein dataset derived from 3,181 completed and draft reference genomes obtained from IMG on 4th September 2012. An expectation value cut-off of 10-30 was utilised to ensure a high confidence level for the assigned functions. Metabolic functions were linked to a sample’s protein sequence fragments using the KEGG database (v58)  with annotations as listed in the IMG database for each genome . If the top hit for an ORF within the reference genome dataset had an associated KEGG Orthologous (KO) group that KO was assigned to the ORF.
A count of each KO within each of the 67 samples was compiled and input to STAMP version 2  in order to detect significant differences in abundances between lean and obese patients, including those that are absent in one but present in the other. Each sample was compared between these two groups using the Welch two-sided t-test with Bonferroni multiple test correction. A cut-off p-value of 0.01 was used to identify KOs whose mean abundance differed significantly between low and high BMI samples.
Phylogenetic reconstruction and taxonomic assignment
Sequences assigned to the same KO set were aligned using ClustalOmega  and then trimmed using BMGE  with an entropy score of 0.7 and a BLOSUM30 matrix. A hidden Markov model was built from this alignment and all metagenome ORF sequences that were assigned a particular KO were aligned to the reference alignment for that KO using hmmalign. Phylogenetic trees were built for each reference KO alignment using FastTree 2.1 with the JTT substitution model and a gamma distribution . In order to calculate bootstrap support, 100 resampled alignments were built per KO using SEQBOOT of the phylip package . FastTree was then used to create a tree per resampled alignment and the original tree was subsequently compared to these 100 resampled trees to infer bootstrap support per node. Subtrees containing only gut-associated species (as listed in the IMG database ) were created and tested for consistency with taxonomy using Chameleon, a visualisation and analysis environment for phylogenetic diversity currently in development.
Classification of metagenomic fragments was undertaken using the Pplacer package v1.1 alpha11 . The taxonomic assignment of each reference sequence was retrieved from the NCBI taxonomy database using Taxtastic (http://fhcrc.github.com/taxtastic) and a Pplacer reference package was created for each KO of interest. Metagenomic sequence fragments were then placed on the tree using Pplacer. This allowed for assignment of each ORF to a taxonomic attribution with a high level of confidence. These classifications were then retrieved using the guppy classification method of Pplacer, which reports the closest taxonomic attribution for each phylogenetically placed read. Differences in abundances of species between lean and obese patients were examined using STAMP version 2 employing the Welch two-sided t-test with Bonferroni multiple test correction and a 0.05 p-value cut-off.
Bäckhed F, Ley RE, Sonnenburg JL, Peterson DA, Gordon JI: Host-bacterial mutualism in the human intestine. Science. 2005, 307: 1915-1920. 10.1126/science.1104816.
Gill SR, Pop M, Deboy RT, Eckburg PB, Turnbaugh PJ, Samuel BS, Gordon JI, Relman DA, Fraser-Liggett CM, Nelson KE: Metagenomic analysis of the human distal gut microbiome. Science. 2006, 312: 1355-1359. 10.1126/science.1124234.
Metges CC: Contribution of Microbial Amino Acids to Amino Acid Homeostasis of the Host. J Nutr. 2000, 130: 1857-1864.
Qin J, Li R, Raes J, Arumugam M, Burgdorf KS, Manichanh C, Nielsen T, Pons N, Levenez F, Yamada T, Mende DR, Li J, Xu J, Li S, Li D, Cao J, Wang B, Liang H, Zheng H, Xie Y, Tap J, Lepage P, Bertalan M, Batto J-M, Hansen T, Le Paslier D, Linneberg A, Nielsen HB, Pelletier E, Renault P, Sicheritz-Ponten T, Turner K, Zhu H, Yu C, Li S, Jian M, Zhou Y, Li Y, Zhang X, Li S, Qin N, Yang H, Wang J, Brunak S, Doré J, Guarner F, Kristiansen K, Pedersen O, Parkhill J, Weissenbach J, Bork P, Ehrlich SD, Wang J: A human gut microbial gene catalogue established by metagenomic sequencing. Nature. 2010, 464: 59-65. 10.1038/nature08821.
Ley RE, Turnbaugh PJ, Klein S, Gordon JI: Human gut microbes associated with obesity. Nature. 2006, 444: 1022-1023. 10.1038/4441022a.
Duncan SH, Lobley GE, Holtrop G, Ince J, Johnstone AM, Louis P, Flint HJ: Human colonic microbiota associated with diet, obesity and weight loss. Int J Obes. 2008, 32: 1720-1724. 10.1038/ijo.2008.155.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, Egholm M, Henrissat B, Heath AC, Knight R, Gordon JI: A core gut microbiome in obese and lean twins. Nature. 2008, 457: 480-484.
Eckburg PB, Bik EM, Bernstein CN, Purdom E, Dethlefsen L, Sargent M, Gill SR, Nelson KE, Relman DA: Diversity of the human intestinal microbial flora. Science. 2005, 308: 1635-1638. 10.1126/science.1110591.
Million M, Maraninchi M, Henry M, Armougom F, Richet H, Carrieri P, Valero R, Raccah D, Vialettes B, Raoult D: Obesity-associated gut microbiota is enriched in Lactobacillus reuteri and depleted in Bifidobacterium animalis and Methanobrevibacter smithii. Int J Obes. 2005, 2011: 1-9.
Armougom F, Henry M, Vialettes B, Raccah D, Raoult D: Monitoring bacterial community of human gut microbiota reveals an increase in Lactobacillus in obese patients and Methanogens in anorexic patients. PLoS One. 2009, 4: e7125-10.1371/journal.pone.0007125.
Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto J-M, Bertalan M, Borruel N, Casellas F, Fernandez L, Gautier L, Hansen T, Hattori M, Hayashi T, Kleerebezem M, Kurokawa K, Leclerc M, Levenez F, Manichanh C, Nielsen HB, Nielsen T, Pons N, Poulain J, Qin J, Sicheritz-Ponten T, Tims S, Torrents D, Ugarte E, Zoetendal EG, Wang J, Guarner F, Pedersen O, de Vos WM, Brunak S, Doré J, Consortium M, Weissenbach J, Ehrlich SD, Bork P, Antolín M, Artiguenave F, Blottiere HM, Almeida M, Brechot C, Cara C, Chervaux C, Cultrone A, Delorme C, Denariaz G, Dervyn R, Foerstner KU, Friss C, van de Guchte M, Guedon E, Haimet F, Huber W, van Hylckama-Vlieg J, Jamet A, Juste C, Kaci G, Knol J, Lakhdari O, Layec S, Le Roux K, Maguin E, Mérieux A, Melo Minardi R, M’rini C, Muller J, Oozeer R, Parkhill J, Renault P, Rescigno M, Sanchez N, Sunagawa S, Torrejon A, Turner K, Vandemeulebrouck G, Varela E, Winogradsky Y, Zeller G: Enterotypes of the Human Gut Microbiome. Nature. 2011, 473: 174-180. 10.1038/nature09944.
Schwiertz A, Taras D, Schäfer K, Beijer S, Bos NA, Donus C, Hardt PD: Microbiota and SCFA in lean and overweight healthy subjects. Obesity. 2010, 18: 190-195. 10.1038/oby.2009.167.
Navarro C, Wu LF, Mandrand-Berthelot MA: The nik operon of Escherichia coli encodes a periplasmic binding-protein-dependent transport system for nickel. Mol Microbiol. 1993, 9: 1181-1191. 10.1111/j.1365-2958.1993.tb01247.x.
Flores-Valdez MA, Morris RP, Laval F, Daffé M, Schoolnik GK: Mycobacterium tuberculosis modulates its cell surface via an oligopeptide permease (Opp) transport system. FASEB J. 2009, 23: 4091-4104. 10.1096/fj.09-132407.
Markowitz VM, Chen I-M A, Palaniappan K, Chu K, Szeto E, Grechkin Y, Ratner A, Jacob B, Huang J, Williams P, Huntemann M, Anderson I, Mavromatis K, Ivanova NN, Kyrpides NC: IMG: the integrated microbial genomes database and comparative analysis system. Nucleic Acids Res. 2012, 40: D115-D122. 10.1093/nar/gkr1044.
Matsen FA, Kodner RB, Armbrust EV: pplacer: linear time maximum-likelihood and Bayesian phylogenetic placement of sequences onto a fixed reference tree. BMC Bioinforma. 2010, 11: 538-10.1186/1471-2105-11-538.
Dinsdale EA, Edwards RA, Hall D, Angly F, Breitbart M, Brulc JM, Furlan M, Desnues C, Haynes M, Li L, McDaniel L, Moran MA, Nelson KE, Nilsson C, Olson R, Paul J, Brito BR, Ruan Y, Swan BK, Stevens R, Valentine DL, Thurber RV, Wegley L, White BA, Rohwer F: Functional metagenomic profiling of nine biomes. Nature. 2008, 452: 629-632. 10.1038/nature06810.
Langille MGI, Meehan CJ, Beiko RG: Human Microbiome: A Genetic Bazaar for Microbes?. Curr Biol. 2012, 22: R20-R22. 10.1016/j.cub.2011.11.023.
Smillie CS, Smith MB, Friedman J, Cordero OX, David LA, Alm EJ: Ecology drives a global network of gene exchange connecting the human microbiome. Nature. 2011, 480: 241-244. 10.1038/nature10571.
Kurokawa K, Itoh T, Kuwahara T, Oshima K, Toh H, Toyoda A, Takami H, Morita H, Sharma VK, Srivastava TP, Taylor TD, Noguchi H, Mori H, Ogura Y, Ehrlich DS, Itoh K, Takagi T, Sakaki Y, Hayashi T, Hattori M: Comparative metagenomics revealed commonly enriched gene sets in human gut microbiomes. DNA Res. 2007, 14: 169-181. 10.1093/dnares/dsm018.
Hess M, Sczyrba A, Egan R, Kim T-W, Chokhawala H, Schroth G, Luo S, Clark DS, Chen F, Zhang T, Mackie RI, Pennacchio LA, Tringe SG, Visel A, Woyke T, Wang Z, Rubin EM: Metagenomic discovery of biomass-degrading genes and genomes from cow rumen. Science. 2011, 331: 463-467. 10.1126/science.1200387.
Mira A, Pushker R, Legault BA, Moreira D, Rodríguez-Valera F: Evolutionary relationships of Fusobacterium nucleatum based on phylogenetic analysis and comparative genomics. BMC Evol Biol. 2004, 4: 50-10.1186/1471-2148-4-50.
Yap WH, Zhang Z, Wang Y: Distinct types of rRNA operons exist in the genome of the actinomycete Thermomonospora chromogena and evidence for horizontal transfer of an entire rRNA operon. J Bacteriol. 1999, 181: 5201-5209.
Jain R, Rivera MC, Lake JA: Horizontal gene transfer among genomes: the complexity hypothesis. Proc Natl Acad Sci U S A. 1999, 96: 3801-3806. 10.1073/pnas.96.7.3801.
Wellner A, Gophna U: Neutrality of foreign complex subunits in an experimental model of lateral gene transfer. Mol Biol Evol. 2008, 25: 1835-1840. 10.1093/molbev/msn131.
Omer S, Kovacs A, Mazor Y, Gophna U: Integration of a foreign gene into a native complex does not impair fitness in an experimental model of lateral gene transfer. Mol Biol Evol. 2010, 27: 2441-2445. 10.1093/molbev/msq145.
Hausinger RP: Nickel utilization by microorganisms. Microbiol Rev. 1987, 51: 22-42.
Duncan SH, Hold GL, Harmsen HJM, Stewart CS, Flint HJ: Growth requirements and fermentation products of Fusobacterium prausnitzii, and a proposal to reclassify it as Faecalibacterium prausnitzii gen. nov., comb. nov. Int J Syst Evol Microbiol. 2002, 52: 2141-2146. 10.1099/ijs.0.02241-0.
O’Dell GD, Miller WJ, King WA, Moore SL, Blackmon DM: Nickel toxicity in the young bovine. J Nutr. 1970, 100: 1447-1453.
Duncan SH, Belenguer A, Holtrop G, Johnstone AM, Flint HJ, Lobley GE: Reduced dietary intake of carbohydrates by obese subjects results in decreased concentrations of butyrate and butyrate-producing bacteria in feces. Appl Environ Microbiol. 2007, 73: 1073-1078. 10.1128/AEM.02340-06.
Gao Z, Yin J, Zhang J, Ward RE, Martin RJ, Lefevre M, Cefalu WT, Ye J: Butyrate Improves Insulin Sensitivity and Increases Energy Expenditure in Mice. Diabetes. 2009, 58: 1509-1517. 10.2337/db08-1637.
Hiles ID, Gallagher MP, Jamieson DJ, Higgins CF: Molecular characterization of the oligopeptide permease of Salmonella typhimurium. J Mol Biol. 1987, 195: 125-142. 10.1016/0022-2836(87)90332-9.
Turnbaugh PJ, Ley RE, Hamady M, Fraser-Liggett CM, Knight R, Gordon JI: The human microbiome project. Nature. 2007, 449: 804-810. 10.1038/nature06244.
Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermúdez-Humarán LG, Gratadoux J-J, Blugeon S, Bridonneau C, Furet J-P, Corthier G, Grangette C, Vasquez N, Pochart P, Trugnan G, Thomas G, Blottière HM, Doré J, Marteau P, Seksik P, Langella P: Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc Natl Acad Sci U S A. 2008, 105: 16731-16736. 10.1073/pnas.0804812105.
Richards VP, Lang P, Pavinski Bitar PD, Lefébure T, Schukken YH, Zadoks RN, Stanhope MJ: Comparative genomics and the role of lateral gene transfer in the evolution of bovine adapted Streptococcus agalactiae. Infect Genet Evol: J Mol Epidemiol Evol Genet Infect Dis. 2011, 11: 1263-1275.
Lurie-Weinberger MN, Peeri M, Gophna U: Contribution of lateral gene transfer to the gene repertoire of a gut-adapted methanogen. Genomics. 2011, 99: 52-58.
Hoff KJ, Lingner T, Meinicke P, Tech M: Orphelia: predicting genes in metagenomic sequencing reads. Nucleic Acids Res. 2009, 37: W101-W105. 10.1093/nar/gkp327.
Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26: 2460-2461. 10.1093/bioinformatics/btq461.
Kanehisa M, Goto S, Kawashima S, Okuno Y, Hattori M: The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004, 32: D277-D280. 10.1093/nar/gkh063.
Parks DH, Beiko RG: Identifying biologically relevant differences between metagenomic communities. Bioinformatics. 2010, 26: 715-721. 10.1093/bioinformatics/btq041.
Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Söding J, Thompson JD, Higgins DG: Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011, 7: 539-
Criscuolo A, Gribaldo S: BMGE (Block Mapping and Gathering with Entropy): a new software for selection of phylogenetic informative regions from multiple sequence alignments. BMC Evol Biol. 2010, 10: 210-10.1186/1471-2148-10-210.
Price MN, Dehal PS, Arkin AP: FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS One. 2010, 5: e9490-10.1371/journal.pone.0009490.
Felsenstein J: PHYLIP - Phylogeny Inference Package (Version 3.2). Cladistics. 1989, 5: 164-166.
We would like to thank Donovan Parks, Robert Eveleigh, Morgan Langille and Erick Matsen for assistance with statistical analysis, alignment processing, phylogenetic clustering and taxonomic assignments.
This work is supported by CIHR grant number CMF-108026. RGB acknowledges the support of Genome Atlantic and the Canada Research Chairs program.
The authors declare that they have no competing interests.
CJM carried out the study design, analysis, and manuscript preparation and editing. RGB contributed to study design, and manuscript preparation and editing. Both authors read and approved the final manuscript.