<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2016.01058</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Variations in the Post-weaning Human Gut Metagenome Profile As Result of <italic>Bifidobacterium</italic> Acquisition in the Western Microbiome</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Soverini</surname> <given-names>Matteo</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Rampelli</surname> <given-names>Simone</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/191006/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Turroni</surname> <given-names>Silvia</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/158027/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Schnorr</surname> <given-names>Stephanie L.</given-names></name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/355582/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Quercia</surname> <given-names>Sara</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/155197/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Castagnetti</surname> <given-names>Andrea</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Biagi</surname> <given-names>Elena</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/359632/overview"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Brigidi</surname> <given-names>Patrizia</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/190984/overview"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>Candela</surname> <given-names>Marco</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>&#x002A;</sup></xref>
<uri xlink:href="http://loop.frontiersin.org/people/157974/overview"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>Department of Pharmacy and Biotechnology, University of Bologna</institution> <country>Bologna, Italy</country></aff>
<aff id="aff2"><sup>2</sup><institution>Department of Anthropology, University of Oklahoma</institution> <country>Norman, OK, USA</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: <italic>Francesca Turroni, University College Cork, Ireland</italic></p></fn>
<fn fn-type="edited-by"><p>Reviewed by: <italic>Abelardo Margolles, Consejo Superior de Investigaciones Cient&#x00ED;ficas, Spain; Jeffrey David Galley, Baylor College of Medicine, USA</italic></p></fn>
<fn fn-type="corresp" id="fn001"><p>&#x002A;Correspondence: <italic>Simone Rampelli, <email>simone.rampelli@unibo.it</email> Marco Candela, <email>marco.candela@unibo.it</email></italic></p></fn>
<fn fn-type="other" id="fn002"><p>This article was submitted to Microbial Symbioses, a section of the journal Frontiers in Microbiology</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>12</day>
<month>07</month>
<year>2016</year>
</pub-date>
<pub-date pub-type="collection">
<year>2016</year>
</pub-date>
<volume>07</volume>
<elocation-id>1058</elocation-id>
<history>
<date date-type="received">
<day>29</day>
<month>04</month>
<year>2016</year>
</date>
<date date-type="accepted">
<day>23</day>
<month>06</month>
<year>2016</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2016 Soverini, Rampelli, Turroni, Schnorr, Quercia, Castagnetti, Biagi, Brigidi and Candela.</copyright-statement>
<copyright-year>2016</copyright-year>
<copyright-holder>Soverini, Rampelli, Turroni, Schnorr, Quercia, Castagnetti, Biagi, Brigidi and Candela</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>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) or licensor 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.</p></license>
</permissions>
<abstract>
<p>Studies of the gut microbiome variation among human populations revealed the existence of robust compositional and functional layouts matching the three subsistence strategies that describe a trajectory of changes across our recent evolutionary history: hunting and gathering, rural agriculture, and urban post-industrialized agriculture. In particular, beside the overall reduction of ecosystem diversity, the gut microbiome of Western industrial populations is typically characterized by the loss of <italic>Treponema</italic> and the acquisition of <italic>Bifidobacterium</italic> as an abundant inhabitant of the post-weaning gut microbial ecosystem. In order to advance the hypothesis about the possible adaptive nature of this exchange, here we explore specific functional attributes that correspond to the mutually exclusive presence of <italic>Treponema</italic> and <italic>Bifidobacterium</italic> using publically available gut metagenomic data from Hadza hunter-gatherers and urban industrial Italians. According to our findings, <italic>Bifidobacterium</italic> provides the enteric ecosystem with a diverse panel of saccharolytic functions, well suited to the array of gluco- and galacto-based saccharides that abound in the Western diet. On the other hand, the metagenomic functions assigned to <italic>Treponema</italic> are more predictive of a capacity to incorporate complex polysaccharides, such as those found in unrefined plant foods, which are consistently incorporated in the Hadza diet. Finally, unlike <italic>Treponema</italic>, the <italic>Bifidobacterium</italic> metagenome functions include genes that permit the establishment of microbe&#x2013;host immunological cross-talk, suggesting recent co-evolutionary events between the human immune system and <italic>Bifidobacterium</italic> that are adaptive in the context of agricultural subsistence and sedentary societies.</p>
</abstract>
<kwd-group>
<kwd><italic>Bifidobacterium</italic></kwd>
<kwd><italic>Treponema</italic></kwd>
<kwd>gut microbiota</kwd>
<kwd>microbiome</kwd>
<kwd>co-evolution</kwd>
<kwd>Westernization</kwd>
</kwd-group>
<counts>
<fig-count count="3"/>
<table-count count="0"/>
<equation-count count="0"/>
<ref-count count="34"/>
<page-count count="9"/>
<word-count count="0"/>
</counts>
</article-meta>
</front>
<body>
<sec><title>Introduction</title>
<p>The gut microbiota (GM) exerts a vital role in human physiology, being strategic for human nutrition, immune protection, and for the preservation of essential metabolic functions (<xref ref-type="bibr" rid="B13">Kau et al., 2011</xref>; <xref ref-type="bibr" rid="B20">Nicholson et al., 2012</xref>). This raises questions of whether, and to what degree, the human GM can provide the host with extra physiological and metabolic flexibility for the adaptation to different lifestyles and environments. Since 2010, several studies have been conducted with the specific aim to explore GM variation across human populations with different subsistence practices, lifestyles and geographical origin. These human GM surveys revealed the existence of robust GM compositional and functional subgroups, so far generally reflective of the variations in subsistence strategy: hunter-gatherer, rural agricultural, and urban industrial Western lifestyle (<xref ref-type="bibr" rid="B33">Yatsunenko et al., 2012</xref>; <xref ref-type="bibr" rid="B25">Schnorr et al., 2014</xref>; <xref ref-type="bibr" rid="B16">Mart&#x00ED;nez et al., 2015</xref>; <xref ref-type="bibr" rid="B21">Obregon-Tito et al., 2015</xref>; <xref ref-type="bibr" rid="B23">Rampelli et al., 2015</xref>).</p>
<p>The findings from this new and emerging field of research, which combines human microbiology ecology and anthropology, resulted in two main conclusions with important implications for both human evolutionary history and human health: first, humans co-evolved with symbiont microbial ecosystems, which have co-adapted along the trajectory of subsistence change across human evolutionary history, from hunter-gatherers to rural agricultural to the most recent development of completely industrialized societies (<xref ref-type="bibr" rid="B22">Quercia et al., 2014</xref>; <xref ref-type="bibr" rid="B21">Obregon-Tito et al., 2015</xref>); second, despite the considerable variation in rural and traditional life-ways, urban industrial populations stand apart as having a distinctly altered GM profile. Indeed, the GM of urban industrial populations seems to universally share certain compositional qualities, such as: (i) an overall compression of microbial diversity as measured by phylogeny and the number of unique taxa (<xref ref-type="bibr" rid="B27">Segata, 2015</xref>), (ii) the loss of the so-called microorganisms &#x201C;old friends&#x201D;, <italic>Treponema</italic> and <italic>Succinivibrio</italic> (<xref ref-type="bibr" rid="B3">Blaser and Falkow, 2009</xref>; <xref ref-type="bibr" rid="B32">Warinner et al., 2015</xref>), and (iii) the acquisition of <italic>Bifidobacterium</italic> as typical inhabitant of the adult gut (<xref ref-type="bibr" rid="B25">Schnorr et al., 2014</xref>; <xref ref-type="bibr" rid="B16">Mart&#x00ED;nez et al., 2015</xref>; <xref ref-type="bibr" rid="B21">Obregon-Tito et al., 2015</xref>).</p>
<p>In agreement with the multiple hit hypothesis (<xref ref-type="bibr" rid="B29">Sonnenburg and Sonnenburg, 2014</xref>), dietary changes, sanitization and antibiotic usage are all potential triggers that would explain the reduction of ecosystem diversity and the loss of co-adapted microbial communities commonly observed in the GM profile of urban industrial Western populations. Conversely, the factors that favored the correspondent acquisition of <italic>Bifidobacterium</italic> in the Western microbiome still remain to be defined. Showing a relative abundance that ranges from 3 to 10% of the total ecosystem, bifidobacteria are an abundant bacterial component in the GM of urban industrial, namely Westernized, adults, and also dominates the GM ecosystem of breast-fed infants, where this bacterial family accounts on average for 80% of the total community (<xref ref-type="bibr" rid="B30">Turroni et al., 2012</xref>; <xref ref-type="bibr" rid="B4">Bottacini et al., 2014</xref>; <xref ref-type="bibr" rid="B18">Milani et al., 2015b</xref>). The recent characterization of the bifidobacterial pangenome &#x2013; 18,181 <italic>Bifidobacterium</italic> specific Cluster of Orthologous Genes (BifCOGs) from 47 sequenced type strains &#x2013; revealed the saccharolytic functions of this microorganism, and indicated a strong adaptation to the human gut environment (<xref ref-type="bibr" rid="B17">Milani et al., 2015a</xref>). Furthermore, recent work has helped advance a possible hypothesis to explain the adoption of <italic>Bifidobacterium</italic> into the Western adult GM. Through comparisons of the GM of Hadza hunter-gatherers and urban industrial Italians, <xref ref-type="bibr" rid="B25">Schnorr et al. (2014)</xref> highlighted &#x2013; for the first time to our knowledge &#x2013; the substantial lack of <italic>Bifidobacterium</italic> from the GM of some traditional populations. The authors proposed that the lack of bifidobacteria in adult Hadza hunter-gatherers may be &#x201C;a consequence of the post-weaning GM composition in the absence of agro-pastoral-derived foods&#x201D;, while the continued consumption of dairy into adulthood is one of the possible vectors by which many Westernized populations maintain a relatively large bifidobacterial presence. <xref ref-type="bibr" rid="B21">Obregon-Tito et al. (2015)</xref> and <xref ref-type="bibr" rid="B16">Mart&#x00ED;nez et al. (2015)</xref> reached similar conclusions through comparative research on the GM of hunter-gatherers from Peru and rural highlanders from Papua New Guinea, respectively, and also found a low proportion of bifidobacteria relative to individuals from the USA. However, these comparative gut metagenome surveys do not specifically explore the impact of <italic>Bifidobacterium</italic> acquisition on the functional configuration of the GM of Western adults.</p>
<p>Taken together, these findings indicate a scenario of exchange, whereby through Westernization, human populations have lost long-standing commensal microorganisms, in particular <italic>Treponema</italic> and <italic>Succinivibrio</italic>, but have potentially compensated through the adult acquisition of bifidobacteria. In order to examine changes to the GM as a result of these community shifts, we investigate how the loss of <italic>Treponema</italic> and the acquisition of <italic>Bifidobacterium</italic> influenced the human gut metagenome profile. To this aim, we compare gut metagenome functions assigned to <italic>Treponema</italic> and <italic>Bifidobacterium</italic> retrieved from downloadable GM metagenomic data for both Hadza hunter-gatherers and urban Italians. Our findings reveal interesting functional gains in the Western microbiome corresponding to the post-weaning retention of <italic>Bifidobacterium</italic> as a symbiont microorganism, suggesting an opportunistic yet important role of this taxon in our recent history.</p>
</sec>
<sec id="s1" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec><title>Sample Collection and Shotgun Sequencing</title>
<p>The Illumina shotgun sequences used in this study were previously generated (<xref ref-type="bibr" rid="B23">Rampelli et al., 2015</xref>), and are publically available at the National Center for Biotechnology Information &#x2013; Sequence Read Archive (NCBI SRA; SRP056480, Bioproject ID PRJNA278393. Leipzig Ethik-kommission review board, reference number 164-12-21052012).</p>
</sec>
<sec><title><italic>Bifidobacterium</italic> and <italic>Treponema</italic> Species Identification within Italian and Hadza Metagenomes</title>
<p>In order to identify the <italic>Bifidobacterium</italic> and <italic>Treponema</italic> species in Italian and Hadza populations, respectively, the 16S rDNA sequences within the assembled metagenomes were taxonomically selected using the assign_taxonomy.py script of the Qiime pipeline (<xref ref-type="bibr" rid="B5">Caporaso et al., 2010</xref>), against the Greengenes database<sup><xref ref-type="fn" rid="fn01">1</xref></sup>. The assignment at species level was performed by blastn of the <italic>Treponema</italic> and <italic>Bifidobacterium</italic> 16S rDNA sequences against the entire NCBI nucleotide database, in particular the top hit results for each sequence were considered.</p>
</sec>
<sec><title>Characterization of the CAZyme Repertoire Assigned to <italic>Bifidobacterium</italic> and <italic>Treponema</italic> in the Gut Metagenome</title>
<p>Reads from a total of 38 individual GM metagenomes, 27 Hadza and 11 Italians from <xref ref-type="bibr" rid="B23">Rampelli et al. (2015)</xref>, were downloaded and included in this study. Reads were assembled into contigs using MetaVelvet (<xref ref-type="bibr" rid="B19">Namiki et al., 2012</xref>) with 350 bp as insert length. Predicted open reading frames (ORFs) were determined by FragGeneScan (<xref ref-type="bibr" rid="B24">Rho et al., 2010</xref>) on assembled contigs, using the &#x201C;&#x2013;w 0&#x201D; option for the fragmented genomic sequences and the parameter &#x201C;&#x2013;t complete&#x201D;. From the translated ORFs we detected the CAZymes-coding sequences using hmmscan tool from the HMMER software package (<xref ref-type="bibr" rid="B9">Eddy, 2011</xref>) and the dbCAN database (<xref ref-type="bibr" rid="B34">Yin et al., 2012</xref>). The outputs were processed by the script hmmscan-parser.sh<sup><xref ref-type="fn" rid="fn02">2</xref></sup>, selecting only the sequences that showed a minimum identity of 30% to the query sequences and an alignment length of at least 100 residues. In order to identify CAZymes derived from <italic>Bifidobacterium</italic> and <italic>Treponema</italic>, we retrieved the nucleotide sequences of the CAZymes detected with hmmscan from the FragGeneScan output, and we blasted them against the NCBI nucleotide database. Only the sequences that showed as best hit an assignment to <italic>Bifidobacterium</italic> for the Italian samples or <italic>Treponema</italic> for the Hadza samples were retained for further analysis. On the basis of the coverage of the contigs, we obtained information concerning the abundance of CAZymes. To compare the data among samples, we obtained a normalized CAZyme abundance by dividing the CAZyme coverages of every correspondent contig for the giga-bases of every correspondent sample. Heatmaps and graphs were generated in R using the packages made4 (<xref ref-type="bibr" rid="B7">Culhane et al., 2005</xref>) and stats<sup><xref ref-type="fn" rid="fn03">3</xref></sup>.</p>
</sec>
<sec><title>Read-Mapping Approach for the Detection of <italic>Bifidobacterium</italic> and <italic>Treponema</italic> Functions Involved in the Adaptation to the Gut Environment</title>
<p>High quality reads for each sample were aligned to <italic>Bifidobacterium-</italic> or <italic>Treponema</italic>-assigned genes encoding bile acid adaptation, host interaction, and polysaccharide catabolism using bowtie2 and setting the alignment parameters to &#x201C;&#x2013;sensitive-local&#x201D; (See Supplementary Table <xref ref-type="supplementary-material" rid="SM1">S1</xref> for the list of genes used for this analysis). As reference for the alignment, two different databases containing orthologous genes from the NCBI genomes of the previously detected <italic>Treponema</italic> or <italic>Bifidobacterium</italic> species were created. Specifically, the databases contain genes for alpha-amylase, beta-galactosidase, mannanase, cellulase, pectinase, and xylanase. Furthermore, the databases were implemented with the sequences of the bile e&#xFB04;ux pump, bile salt hydrolase, exopolysaccharide synthase, fimbrial subunit FimQ, sortase, galactosyl transferase, and undecaprenyl-phosphate phosphotransferase, since they were recently reported as genes that facilitate commensal-host cross-talk in <italic>Bifidobacterium</italic> (<xref ref-type="bibr" rid="B12">Gueimonde et al., 2009</xref>; <xref ref-type="bibr" rid="B10">Fanning et al., 2012</xref>; <xref ref-type="bibr" rid="B31">Turroni et al., 2013</xref>; <xref ref-type="bibr" rid="B11">Ferrario et al., 2016</xref>). In the event that the <italic>Bifidobacterium</italic> or <italic>Treponema</italic> NCBI genomes did not contain the above-mentioned genes, we supplemented the databases with genes belonging to the taxonomically closest annotated microorganism. The reads that aligned with a reference using bowtie2, were extracted and their taxonomy was further verified by blastn against the entire NCBI nucleotide database. Notably, in the case that the best hits of the blastn search were not assigned to <italic>Bifidobacterium</italic> or <italic>Treponema</italic>, we did not consider those reads for further analysis. The number of hits for each gene was normalized by the number of base pairs in the input file and in the correspondent reference in order to compare the results.</p>
</sec>
</sec>
<sec><title>Results</title>
<sec><title>Variation in the Gut Microbiome Carbohydrate-Degrading Repertoire as a Result of the Exclusive Presence of <italic>Bifidobacterium</italic> or <italic>Treponema</italic> in the Italian and Hadza Gut Microbial Ecosystems</title>
<p>By studying hunter-gatherers and rural populations, we previously depicted an approximation of the co-evolutional trajectory of the human GM structure in the recent times (<xref ref-type="bibr" rid="B25">Schnorr et al., 2014</xref>). In particular, we characterized the taxonomic and functional metabolic potential of the GM in the Hadza hunter-gatherers and urban Italian adults, reporting considerable differences between the two populations that map on to their respective life-ways (<xref ref-type="bibr" rid="B25">Schnorr et al., 2014</xref>; <xref ref-type="bibr" rid="B23">Rampelli et al., 2015</xref>). Among them, we observed a loss of <italic>Treponema</italic> and an acquisition of <italic>Bifidobacterium</italic> in the Italian GM. Considering the importance of <italic>Bifidobacterium</italic> as a health-promoting commensal of the GM in Western populations, it is important to understand how this beneficial role may have developed through the lens of co-evolution.</p>
<p>We first identified the diversity of the <italic>Bifidobacterium</italic> and <italic>Treponema</italic> species in urban Italian and Hadza GM by reconstructing the full 16S rDNA gene from assembled metagenomes. Italian samples contain sequences belonging to <italic>Bifidobacterium faecale</italic>, <italic>Bifidobacterium pseudocatenulatum</italic>, <italic>Bifidobacterium adolescentis</italic>, <italic>Bifidobacterium coryneforme</italic>, <italic>Bifidobacterium bifidum</italic>, <italic>Bifidobacterium longum</italic>, <italic>Bifidobacterium angulatum</italic>, and <italic>Bifidobacterium dentium</italic>. On the other hand, we detect the 16S rDNA sequences assigned to <italic>Treponema porcinum</italic>, <italic>Treponema bryantii</italic>, <italic>Treponema succinifaciens</italic>, <italic>Treponema parvum</italic>, and <italic>Treponema berlinense</italic> in the GM of the Hadza hunter-gatherers. However, we must acknowledge that these taxonomic assignments are limited by present whole genome for <italic>Treponema</italic> species, most of which have been characterized by work on human pathogens, rather than commensal members of the GM. Ongoing work with species-level variation should be crucial to resolving the disparities in our classification and resolution of these species and their functions in the future (<xref ref-type="bibr" rid="B21">Obregon-Tito et al., 2015</xref>). In order to compare the specific saccharolytic functions conferred by <italic>Bifidobacterium</italic> and <italic>Treponema</italic> in the Italian and Hadza microbiomes, we identified a total of 5.4 million of ORFs, of which 14,512 mapped to CAZymes for the Italian samples and 74,651 for the Hadza samples (<bold>Figure <xref ref-type="fig" rid="F1">1A</xref></bold>). Notably, the Hadza metagenomes contain significantly more CAZymes per subject, in terms of ORFs assigned to CAZymes per million of reads, respect to the Italian metagenomes (mean &#x00B1; SD, Hadza: 233 &#x00B1; 86, Italians: 137 &#x00B1; 78), as reported by <xref ref-type="bibr" rid="B23">Rampelli et al. (2015)</xref>. We then profiled the saccharolytic repertoire of <italic>Bifidobacterium</italic> and <italic>Treponema</italic> in the Italian and Hadza GM as relative abundance at the CAZyme category level based on the coverage of taxonomically assigned contigs (<bold>Figure <xref ref-type="fig" rid="F2">2A</xref></bold>). <italic>Bifidobacterium</italic> showed a higher presence of glycosyl transferase (GT) and carbohydrate esterase (CE), with respect to <italic>Treponema</italic> (GT relative abundance, rel. ab.: 54% for <italic>Bifidobacterium</italic> and 43% for <italic>Treponema</italic>; CE rel. ab.: 6% for <italic>Bifidobacterium</italic> and 5% for <italic>Treponema</italic>). On the other hand, <italic>Treponema</italic> were more enriched in glycoside hydrolase (GH) and carbohydrate binding module (CBM) (GH rel. ab.: 33% for <italic>Treponema</italic> and 31% for <italic>Bifidobacterium</italic>; CBM rel. ab.: 15% for <italic>Treponema</italic> and 8% for <italic>Bifidobacterium</italic>). At the CAZyme family level, we revealed the four families that constitute the core <italic>Bifidobacterium</italic> CAZyme repertoire: GH13 (GH family acting on substrates containing &#x03B1;-glycosidic linkages), GH3 (GH family that groups together exo-acting &#x03B2;-<sc>D</sc>-glucosidases, &#x03B1;-<sc>L</sc>-arabinofuranosidase, &#x03B2;-<sc>D</sc>-xylopyranosidase and <italic>N</italic>-acetyl-&#x03B2;-<sc>D</sc>-glucosaminidase), GT2 (GT family containing cellulose synthase, mannan synthase, and several monosaccharide-/oligosaccharide-transferases), and GT4 (GT family containing sucrose synthase, glucosyl transferase, and several phosphorylases). The sum of ORFs assigned to these four major families comprises 77% of the total detected CAZyme cohort. The relative abundance of the <italic>Bifidobacterium</italic> and <italic>Treponema</italic> CAZyme families detected in the Italian and Hadza samples reveals several differences in the potential saccharolytic functional contributions of these two microorganisms (<bold>Figure <xref ref-type="fig" rid="F2">2B</xref></bold>). In particular, <italic>Bifidobacterium</italic> have a greater abundance of genes involved in the degradation of lactate, which is produced from pyruvate in the fermentation of simple sugar and commonly found in sour milk as well as in other lacto-fermented foods (family GH2). This finding also corresponds to the exclusive contribution of dairy carbohydrates (&#x223C;5% of total carbohydrates), in the Italian cohort diet (Supplementary Table <xref ref-type="supplementary-material" rid="SM2">S2</xref>). Emphasis on monosaccharide catabolism is evidenced by enrichment in gene families containing enzymes that metabolize mannose, xylose and arabinose (GH2, GH31, and GH43), which are highly represented in plant and fruit glycans. As a variety of genes were also found that are involved in the degradation of &#x03B1;- and &#x03B2;-glucans (GH3 and GH31), this illustrates an ability of <italic>Bifidobacterium</italic> to retrieve energy also from more complex carbohydrates that are commonly present in the cellulosic biomass of plant foods in the Italian diet: salads, fruits, nuts, cereals and their product derivatives. In addition, <italic>Bifidobacterium</italic> are also enriched in genes involved in the catabolism of sucrose (GH31), which is widely distributed in nature, but robustly manifest in the industrial food products that are consumed daily by most urban populations. Further evidence of these functions comes from detection of a higher abundance of CBM families for lactose, galactose and &#x03B2;-glucans (CBM4, CBM13, and CBM32) in <italic>Bifidobacterium</italic>, with respect to <italic>Treponema</italic>. In contrast, the CAZyme profile of <italic>Treponema</italic> within the Hadza metagenome is mainly devoted to degradation of glucans, galactans, and fructans (GH16, GH32, and GH53), which are sugar polymers that comprise hemicellulose (galactans) and inulin (fructans). The monosaccharide of galactans, galactose, is also expressed in mucilages and glycoproteins that derive from the human host, as well as a number of vegetable-derived carbohydrates. Both sugar polymers are largely implicit in difficult-to-digest plant polysaccharides that escape small intestine absorption and are instead fermented by the colonic microbiota. The Hadza diet is replete with such unrefined plant foods that contain indigestible polysaccharides such as berries, baobab fruit, and particularly tubers (<xref ref-type="bibr" rid="B26">Schnorr et al., 2015</xref>). <italic>Treponema</italic> are also enriched in two CAZyme &#x03B1;-amylase families (GH57 and GH77), which are unlike the typical &#x03B1;-amylase GH13 family because they have a conserved trans-glycosylating region (<xref ref-type="bibr" rid="B14">MacGregor et al., 2001</xref>). Finally, <italic>Treponema</italic> are better equipped to metabolize peptidoglycans due to a wide range of acetylglucosaminases and peptidoglycan lyases (GH23, GH73, and GH109). These activities were confirmed by the detection of high levels of CBM families for peptidoglycans and &#x03B1;-glucans (CBM50 and CBM48). To specifically investigate the effective differences in polysaccharide-degrading repertoires between <italic>Bifidobacterium</italic> and <italic>Treponema</italic> metagenome functions, we aligned the reads of Italian and Hadza samples to a custom database containing the sequences for alpha-amylase, beta-galactosidase, mannanase, cellulase, pectinase, and xylanase. The taxonomic specificity of the aligned reads was verified by blasting the sequences to the nucleotide database of NCBI (<bold>Figure <xref ref-type="fig" rid="F1">1B</xref></bold>). Unlike <italic>Treponema</italic>, which appears to be functionally equipped to derive energy from a broad spectrum of polysaccharides through increased relative representation of pectinase, xylanase, and cellulase, the polysaccharide-degrading functions assigned to <italic>Bifidobacterium</italic> are enriched in genes coding for &#x03B2;-galactosidase and mannanase (<bold>Figure <xref ref-type="fig" rid="F3">3</xref></bold>). Interestingly, both genera are configured to hydrolyze 1&#x2013;4 &#x03B1;-glycosidic bonds (via &#x03B1;-amylase) to a similar degree. The catabolic configurations of <italic>Bifidobacterium</italic> and <italic>Treponema</italic> for polysaccharides are reflective of differences in dietary carbohydrates of Italians and Hadza: pasta, legumes, milk, and dairy products, versus berries, tubers, baobab, and honey, respectively, (Supplementary Table <xref ref-type="supplementary-material" rid="SM3">S3</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption><p><bold>Schematic representation of the analysis workflow. (A)</bold> Pipeline for the identification and assignment of <italic>Treponema</italic> and <italic>Bifidobacterium</italic> CAZymes on assembled metagenomes: (i) ORFs detection using FragGeneScan; (ii) detection of the CAZyme-coding ORFs by using hmmscan against the dbCAN CAZy database; (iii) taxonomy assignment to CAZyme-coding sequences by blastn against the NCBI nucleotide database. <bold>(B)</bold> Pipeline for the identification of <italic>Treponema</italic> and <italic>Bifidobacterium</italic> sequences coding for functions involved in the adaptation to the gut environment: (i) alignment of high quality reads to databases containing the selected <italic>Treponema</italic> or <italic>Bifidobacterium</italic> functions using bowtie2; (ii) blasting of the successfully aligned reads against the NCBI nucleotide database to confirm the taxonomy.</p></caption>
<graphic xlink:href="fmicb-07-01058-g001.tif"/>
</fig>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption><p><bold>Overview of <italic>Bifidobacterium</italic> and <italic>Treponema</italic> CAZyme repertoires in the Hadza and Italian samples. (A)</bold> Normalized relative abundance of the CAZyme category levels for <italic>Bifidobacterium</italic> and <italic>Treponema</italic>: auxiliary activity (AA), carbohydrate binding module (CBM), carbohydrate esterase (CE), glycoside hydrolase (GH), and glycosyl transferase (GT) categories. <bold>(B)</bold> Comparison between the <italic>Bifidobacterium</italic> and <italic>Treponema</italic> CAZyme family profiles. The relative abundance of each family is indicated by the color key.</p></caption>
<graphic xlink:href="fmicb-07-01058-g002.tif"/>
</fig>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption><p><bold>Profile of <italic>Bifidobacterium</italic> and <italic>Treponema</italic> functions involved in the adaptation to the host environment</bold>. Polysaccharide metabolism <bold>(A)</bold>; interaction with the host and immune modulation <bold>(B)</bold>; adaptation to the human gastrointestinal environment <bold>(C)</bold>. Color key represents reads mapped per giga-bases of the sample of origin. See Supplementary Table <xref ref-type="supplementary-material" rid="SM1">S1</xref> for the list of reference genes used.</p></caption>
<graphic xlink:href="fmicb-07-01058-g003.tif"/>
</fig>
</sec>
<sec><title>Differences in the Gut Metagenome Functions Involved in Bacteria/Host Interaction and Gut Adaptation Corresponding to the Presence of <italic>Bifidobacterium</italic> or <italic>Treponema</italic> in the Italian and Hadza Gut Metagenome</title>
<p>We further investigated the presence of the genes involved in host interaction and immune modulation in the ORFs attributed to <italic>Bifidobacterium</italic> and <italic>Treponema</italic>. In this scenario, the sortase-dependent pili are considered key surface molecules in establishing bacterial adherence to the host epithelium and they have been proposed as potential mediator of mucosal immune response (<xref ref-type="bibr" rid="B31">Turroni et al., 2013</xref>). Another recognized marker of interaction with the host that affects the immune response is the surface-exopolysaccharide (EPS). Specifically, bacteria producing EPS failed to elicit a strong immune B-cell response compared to EPS-deficient strains (<xref ref-type="bibr" rid="B10">Fanning et al., 2012</xref>). Our analyses confirm that the enzymes involved in the production of EPS and pili, namely EPS synthase, undecaprenyl-phosphate phosphotransferase, galactosyl transferase, sortase, and fimbrial subunit FimQ, are typical of the <italic>Bifidobacterium</italic> ORFs detected in the Italian metagenome, while virtually absent in the <italic>Treponema</italic> ORFs retrieved from the Hadza GM ecosystem (<bold>Figure <xref ref-type="fig" rid="F3">3</xref></bold>). Finally, we investigated the presence of genes involved in bile tolerance as mechanisms of bacterial adaptation to the human host. Bile salts are detergent-like compounds with strong antimicrobial activity (<xref ref-type="bibr" rid="B1">Begley et al., 2005</xref>), and intestinal bacteria have had to evolve strategies to tolerate physiological concentrations of bile salts to colonize the intestine (<xref ref-type="bibr" rid="B12">Gueimonde et al., 2009</xref>). Interestingly, two representative enzymes, which contribute to bile resistance and adaptation to gut environment, the bile-inducible e&#xFB04;ux transporters and the bile salt hydrolase, are present in within the <italic>Bifidobacterium</italic> gut metagenome functions, but are not detected in <italic>Treponema</italic> ORFs (<bold>Figure <xref ref-type="fig" rid="F3">3</xref></bold>).</p>
</sec>
</sec>
<sec><title>Discussion</title>
<p>In our work we explore how the apparent dichotomy between two instance-specific mutually exclusive gut inhabitants, <italic>Bifidobacterium</italic> and <italic>Treponema</italic>, reflects specific functional roles within the human GM ecosystem. We observe that the complete replacement of these microorganisms within the human gut in modern populations follows a curious pattern of lifestyle-associated differences, which stand at opposite ends of an entire subsistence spectrum, from hunting and gathering to post-industrial urban life. This may hint at a more ancient pattern of transition that occurred along more recent human evolutionary history as early settlements permitted the adoption of a fully agro-pastoral regime in place of mobile foraging. As seen in modern rural farmers, these earlier human groups may have harbored both communities of bifidobacteria and treponemes, suggesting an intermediary phase of co-habitation (<xref ref-type="bibr" rid="B8">De Filippo et al., 2010</xref>; <xref ref-type="bibr" rid="B21">Obregon-Tito et al., 2015</xref>). These patterns are instructive, and we can hypothesize that the stimulus for bifidobacterial acquisition and proliferation in the human GM followed much earlier changes to diet and lifestyle that occurred during the Neolithic transition. Specifically, these changes include permanent settlements, dense population structures, animal co-habitation, plant domestication, and shifts in the type and amount of carbohydrate-based foods (<xref ref-type="bibr" rid="B6">Childe, 1936</xref>; <xref ref-type="bibr" rid="B2">Binford, 1968</xref>).</p>
<p>While bifidobacteria are absent from the GM ecosystem of modern hunter-gatherers, and thus possibly acquired as a post-weaning symbionts among communities who adopted an agro-pastoral subsistence, <italic>Treponema</italic> has been lost from the human GM along the transition from small-scale rural communities to post-industrial urbanized society. This raises key questions such as &#x2018;what are the functional gains and losses in the human GM corresponding to this process?&#x2019;, and &#x2018;does an adaptive nature of these changes exist?&#x2019; To begin answering these questions, we compared metagenome functions assigned to <italic>Bifidobacterium</italic> and <italic>Treponema</italic> in the gut metagenome samples from urban Italians and Hadza hunter-gatherers, respectively. In particular, we focused on functional features key in the context of the GM-human host mutualism, such as the metabolism of complex polysaccharides, which is essential to provide the host with short-chain fatty acids (SCFA), and the microbe&#x2013;host interaction and immunomodulation processes.</p>
<p>As an initial approach, we described the <italic>Bifidobacterium</italic> and <italic>Treponema</italic> diversity in the urban Italian and Hadza GM. Eight different bifidobacterial species were detected in the Italian GM ecosystem, including the plant polysaccharides degrader <italic>B. adolescentis</italic>, and the milk fermenters <italic>B. bifidum</italic> and <italic>B. longum</italic>. Other bifidobacterial species detected in the Italian GM were <italic>B. faecale</italic>, <italic>B. pseudocatenulatum</italic>, <italic>B. coryneforme</italic>, <italic>B. angulatum</italic>, and <italic>B. dentium</italic>, highlighting considerable bifidobacterial diversity within the Italian GM ecosystem. Given our current limitation in resolving novel <italic>Treponema</italic> taxa, we were still able to identify five species: <italic>T. porcinum</italic>, <italic>T. bryantii</italic>, <italic>T. succinifaciens</italic>, <italic>T. parvum</italic>, and <italic>T. berlinense</italic>. These species were already described in association with the GM of Matses hunter-gatherers from Peru (<xref ref-type="bibr" rid="B21">Obregon-Tito et al., 2015</xref>).</p>
<p>According to our findings, <italic>Bifidobacterium</italic> and <italic>Treponema</italic> provide the host metagenome with a different configuration of CAZyme categories. Indeed, the <italic>Bifidobacterium</italic> genes encoding for CAZymes in the Italian gut metagenome are enriched in GT and CE, while the <italic>Treponema</italic> saccharolytic repertoire in the Hadza gut is principally devoted to GH and CBM functions, suggesting the higher propensity of the latter to act as a primary degrader of complex polysaccharides in the gut. Focusing on the CAZyme family level, we observed that <italic>Bifidobacterium</italic> provides the host metagenome with an extremely versatile panel of saccharolytic functions. Bifidobacterial CAZyme families enriched in the Italian samples range from the metabolism of simple sugars, such as lactate, mannose, xylose, and arabinose, to more complex carbohydrates from plant sources, such as &#x03B1;- and &#x03B2;-glucans, lactose, or galactose. These findings are in general agreement with the overall structure of the bifidobacterial glycobiome (<xref ref-type="bibr" rid="B17">Milani et al., 2015a</xref>). According to <xref ref-type="bibr" rid="B17">Milani et al. (2015a)</xref>, when comparing the bifidobacterial GH repertoire with the sequenced members of the human GM, <italic>Bifidobacterium</italic> appears enriched in GH3 and GH43 families for the degradation of plant carbohydrates. These findings are in agreement with our observations that the <italic>Bifidobacterium</italic>-assigned saccharolytic functions emphasize substrate-specific specializations that are not observed among <italic>Treponema</italic>-assigned genes. Instead, the saccharolytic functions supplied by <italic>Treponema</italic> to the Hadza gut metagenome are mainly devoted to the degradation of indigestible polysaccharides, such as galactans, fructans, and glucans, that are present in unrefined and wild plant foods. Other functional peculiarities characteristic of <italic>Treponema</italic> include a completely unique suite of &#x03B1;-amylases and a greater capacity to bind and metabolize peptidoglycan. When we specifically explored gut metagenome differences in polysaccharide-degrading genes, we confirmed the overall greater metabolic polysaccharolytic potential for <italic>Treponema</italic> while <italic>Bifidobacterium</italic> was enriched in beta-galactosidase and mannase.</p>
<p>Interestingly, the specific glycobiomes conferred by <italic>Bifidobacterium</italic> and <italic>Treponema</italic> match the respective dietary habits of the two populations. Hadza consume a heavily plant-based diet, particularly in the rainy season, in which &#x223C;70% of kcal are derived from fibrous wild plant foods and the 30% from wild game meat (<xref ref-type="bibr" rid="B15">Marlowe, 2010</xref>). Furthermore, Hadza entirely lack dairy in their diet. Correspondingly, the <italic>Treponema</italic>-assigned glycobiome in the Hadza is devoted to the metabolism of the vast array of refractory glucans, galactans, and fructans that are present in the wild Hadza plant foods. The Mediterranean diet, characteristic for the Italian cohort, is abundant in plant foods (salads, fruits, sauces), pasta, bread and olive oil with low to moderate inclusion of dairy, poultry, fish and read meat (Supplementary Table <xref ref-type="supplementary-material" rid="SM3">S3</xref>). Curiously, certain glycobiome functions conferred by bifidobacteria to the Italian GM &#x2013; e.g., GH2, GH31, and GH43, as well as CBM4, CBM13, and CBM32 &#x2013; are well suited to deal with dairy carbohydrates and some types of plant glycans, all of which are relatively abundant in the Italian diet.</p>
<p>When we looked for functions associated with host&#x2013;microbe interactions, we found matches only among <italic>Bifidobacterium</italic>. Conversely, <italic>Treponema</italic> does not provide the Hadza metagenome with known functionalities for direct host interaction. The <italic>Bifidobacterium</italic> functional potential encoded within strains retrieved from the Italian samples included enzymes for the biosynthesis of pili and EPS structures. Both cell components are essential in the context of a mutualistic strategy that permits close interaction and initiation of a tolerated immunological dialog between the bacterium and host (<xref ref-type="bibr" rid="B10">Fanning et al., 2012</xref>). Finally, as a functional marker of the bacterial adaptation to the gut environment (<xref ref-type="bibr" rid="B1">Begley et al., 2005</xref>), we sought functions involved in bile salts tolerance. Interestingly, only bifidobacteria seem to possess bile-inducible e&#xFB04;ux transporters and the bile salt hydrolase required for bile salt mitigation and detoxification. Even if biases in the assignment of <italic>Treponema</italic> functions as a result of the paucity of reference genomes cannot be excluded, our findings could indicate different ecological strategies for <italic>Bifidobacterium</italic> and <italic>Treponema</italic> in the human gut.</p>
</sec>
<sec><title>Conclusion</title>
<p>Our findings suggest possible co-evolutionary implications for the loss of <italic>Treponema</italic> and the acquisition of <italic>Bifidobacterium</italic> as a stable component of the post-weaning GM ecosystem from post-industrial urban populations. Capable of heterogeneous saccharolytic metabolism, which ranges from complex plant polysaccharides to simpler sugars such as lactose and sucrose, <italic>Bifidobacterium</italic> are well suited to handle the degradative demands imposed by a typical Western diet. Conversely, the progressive loss of more challenging microbiota accessible carbohydrates in the Western diet (<xref ref-type="bibr" rid="B28">Sonnenburg et al., 2016</xref>), such as hemicellulose and inulin, would help partially explain the extinction of a more specialized fiber degrader such as <italic>Treponema</italic> from the Western GM ecosystem. Furthermore, unlike <italic>Treponema</italic>, <italic>Bifidobacterium</italic> evolved the capacity to establish an intense microbe&#x2013;host connection, which may help support a continuous and abundant bifidobacterial presence in adults, allowing this commensal to outcompete other opportunistic, but functionally diverse, microbiota. The acquisition of <italic>Bifidobacterium</italic> as a stable component of the GM ecosystem in small-scale rural agriculturalists, reminiscent of early human farmers, and modern Westernized populations, may therefore engage the functionalities of the host immune system, providing new adaptive solutions in response to changing selective pressures during the restructuring of human diet and society.</p>
</sec>
<sec><title>Author Contributions</title>
<p>SR and MC conceived the study with the participation of MS in experimental design. MS and SR performed the bioinformatics analysis. MS, SR, SLS, ST, and MC carried out data interpretation and wrote the paper. SQ, AC, EB, and PB participated to the discussion of the results and in revising the final draft of the paper. All authors read and approved the final manuscript.</p>
</sec>
<sec><title>Conflict of Interest Statement</title>
<p>The 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.</p>
</sec>
</body>
<back>
<sec sec-type="supplementary material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="http://journal.frontiersin.org/article/10.3389/fmicb.2016.01058">http://journal.frontiersin.org/article/10.3389/fmicb.2016.01058</ext-link></p>
<supplementary-material xlink:href="Table_1.XLSX" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_2.PDF" id="SM2" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_3.PDF" id="SM3" mimetype="application/pdf" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Begley</surname> <given-names>M.</given-names></name> <name><surname>Gahan</surname> <given-names>C. G.</given-names></name> <name><surname>Hill</surname> <given-names>C.</given-names></name></person-group> (<year>2005</year>). <article-title>The interaction between bacteria and bile.</article-title> <source><italic>FEMS Microbiol. Rev.</italic></source> <volume>29</volume> <fpage>625</fpage>&#x2013;<lpage>651</lpage>. <pub-id pub-id-type="doi">10.1016/j.femsre.2004.09.003</pub-id></citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Binford</surname> <given-names>L. R.</given-names></name></person-group> (<year>1968</year>). <article-title>&#x201C;Post pleistocene adaptations,&#x201D; in</article-title> <source><italic>New Perspectives in Archaeology</italic>,</source> <role>eds</role> <person-group person-group-type="editor"><name><surname>Binford</surname> <given-names>S.</given-names></name> <name><surname>Binford</surname> <given-names>L. R.</given-names></name></person-group> (<publisher-loc>Chicago, IL</publisher-loc>: <publisher-name>Adeline</publisher-name>), <fpage>313</fpage>&#x2013;<lpage>342</lpage>.</citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Blaser</surname> <given-names>M. J.</given-names></name> <name><surname>Falkow</surname> <given-names>S.</given-names></name></person-group> (<year>2009</year>). <article-title>What are the consequences of the disappearing human microbiota?</article-title> <source><italic>Nat. Rev. Microbiol.</italic></source> <volume>7</volume> <fpage>887</fpage>&#x2013;<lpage>894</lpage>. <pub-id pub-id-type="doi">10.1038/nrmicro2245</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bottacini</surname> <given-names>F.</given-names></name> <name><surname>Ventura</surname> <given-names>M.</given-names></name> <name><surname>Sinderen</surname> <given-names>D.</given-names></name> <name><surname>Motherway</surname> <given-names>M.</given-names></name></person-group> (<year>2014</year>). <article-title>Diversity, ecology and intestinal function of bifidobacteria.</article-title> <source><italic>Microb. Cell Fact.</italic></source> <volume>13</volume>:<issue>S4</issue>. <pub-id pub-id-type="doi">10.1186/1475-2859-13-S1-S4</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Caporaso</surname> <given-names>J. G.</given-names></name> <name><surname>Kuczynski</surname> <given-names>J.</given-names></name> <name><surname>Stombaugh</surname> <given-names>J.</given-names></name> <name><surname>Bittinger</surname> <given-names>K.</given-names></name> <name><surname>Bushman</surname> <given-names>F. D.</given-names></name> <name><surname>Costello</surname> <given-names>E. K.</given-names></name><etal/></person-group> (<year>2010</year>). <article-title>QIIME allows analysis of high-throughput community sequencing data.</article-title> <source><italic>Nat. Methods</italic></source> <volume>7</volume> <fpage>335</fpage>&#x2013;<lpage>336</lpage>. <pub-id pub-id-type="doi">10.1038/nmeth.f.303</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Childe</surname> <given-names>V. G.</given-names></name></person-group> (<year>1936</year>). <source><italic>Man Makes Himself.</italic></source> <publisher-loc>London</publisher-loc>: <publisher-name>Watts</publisher-name>.</citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Culhane</surname> <given-names>A. C.</given-names></name> <name><surname>Thioulouse</surname> <given-names>J.</given-names></name> <name><surname>Perri&#x00E8;re</surname> <given-names>G.</given-names></name> <name><surname>Higgins</surname> <given-names>D. G.</given-names></name></person-group> (<year>2005</year>). <article-title>MADE4: an R package for multivariate analysis of gene expression data.</article-title> <source><italic>Bioinformatics</italic></source> <volume>21</volume> <fpage>2789</fpage>&#x2013;<lpage>2790</lpage>. <pub-id pub-id-type="doi">10.1093/bioinformatics/bti394</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>De Filippo</surname> <given-names>C.</given-names></name> <name><surname>Cavalieri</surname> <given-names>D.</given-names></name> <name><surname>Di Paola</surname> <given-names>M.</given-names></name> <name><surname>Ramazzotti</surname> <given-names>M.</given-names></name> <name><surname>Poullet</surname> <given-names>J. B.</given-names></name> <name><surname>Massart</surname> <given-names>S.</given-names></name><etal/></person-group> (<year>2010</year>). <article-title>Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>107</volume> <fpage>14691</fpage>&#x2013;<lpage>14696</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1005963107</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Eddy</surname> <given-names>S. R.</given-names></name></person-group> (<year>2011</year>). <article-title>Accelerated profile HMM searches.</article-title> <source><italic>PLoS Comput. Biol.</italic></source> <volume>7</volume>:<issue>e1002195</issue>. <pub-id pub-id-type="doi">10.1371/journal.pcbi.1002195</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fanning</surname> <given-names>S.</given-names></name> <name><surname>Hall</surname> <given-names>L. J.</given-names></name> <name><surname>Cronin</surname> <given-names>M.</given-names></name> <name><surname>Zomer</surname> <given-names>A.</given-names></name> <name><surname>MacSharry</surname> <given-names>J.</given-names></name> <name><surname>Goulding</surname> <given-names>D.</given-names></name><etal/></person-group> (<year>2012</year>). <article-title>Bifidobacterial surface-exopolysaccharide facilitates commensal-host interaction through immune modulation and pathogen protection.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>109</volume> <fpage>2108</fpage>&#x2013;<lpage>2113</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1115621109</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ferrario</surname> <given-names>C.</given-names></name> <name><surname>Milani</surname> <given-names>C.</given-names></name> <name><surname>Mancabelli</surname> <given-names>L.</given-names></name> <name><surname>Lugli</surname> <given-names>G. A.</given-names></name> <name><surname>Duranti</surname> <given-names>S.</given-names></name> <name><surname>Mangifesta</surname> <given-names>M.</given-names></name><etal/></person-group> (<year>2016</year>). <article-title>Modulation of the eps-ome transcription of bifidobacteria through simulation of human intestinal environment.</article-title> <source><italic>FEMS Microbiol. Ecol</italic></source> <volume>92</volume>:<issue>fiw056</issue>. <pub-id pub-id-type="doi">10.1093/femsec/fiw056</pub-id></citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gueimonde</surname> <given-names>M.</given-names></name> <name><surname>Garrigues</surname> <given-names>C.</given-names></name> <name><surname>van Sinderen</surname> <given-names>D.</given-names></name> <name><surname>de los Reyes-Gavilan</surname> <given-names>C. G.</given-names></name> <name><surname>Margolles</surname> <given-names>A.</given-names></name></person-group> (<year>2009</year>). <article-title>Bile-inducible e&#xFB04;ux transporter from <italic>Bifidobacterium longum</italic> NCC2705, conferring bile resistance.</article-title> <source><italic>Appl. Environ. Microbiol.</italic></source> <volume>75</volume> <fpage>3153</fpage>&#x2013;<lpage>3160</lpage>. <pub-id pub-id-type="doi">10.1128/AEM.00172-09</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kau</surname> <given-names>A. L.</given-names></name> <name><surname>Ahern</surname> <given-names>P. P.</given-names></name> <name><surname>Griffin</surname> <given-names>N. W.</given-names></name> <name><surname>Goodman</surname> <given-names>A. L.</given-names></name> <name><surname>Gordon</surname> <given-names>J. I.</given-names></name></person-group> (<year>2011</year>). <article-title>Human nutrition, the gut microbiome and the immune system.</article-title> <source><italic>Nature</italic></source> <volume>474</volume> <fpage>327</fpage>&#x2013;<lpage>336</lpage>. <pub-id pub-id-type="doi">10.1038/nature10213</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>MacGregor</surname> <given-names>E. A.</given-names></name> <name><surname>Janecek</surname> <given-names>S.</given-names></name> <name><surname>Svensson</surname> <given-names>B.</given-names></name></person-group> (<year>2001</year>). <article-title>Relationship of sequence and structure to specificity in the alpha-amylase family of enzymes.</article-title> <source><italic>Biochim. Biophys. Acta</italic></source> <volume>1546</volume> <fpage>1</fpage>&#x2013;<lpage>20</lpage>. <pub-id pub-id-type="doi">10.1016/S0167-4838(00)00302-2</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Marlowe</surname> <given-names>F. W.</given-names></name></person-group> (<year>2010</year>). <source><italic>The Hadza.</italic></source> <publisher-loc>Oakland, CA</publisher-loc>: <publisher-name>University of California Press</publisher-name>.</citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Mart&#x00ED;nez</surname> <given-names>I.</given-names></name> <name><surname>Stegen</surname> <given-names>J. C.</given-names></name> <name><surname>Maldonado-G&#x00F3;mez</surname> <given-names>M. X.</given-names></name> <name><surname>Eren</surname> <given-names>A. M.</given-names></name> <name><surname>Siba</surname> <given-names>P. M.</given-names></name> <name><surname>Greenhill</surname> <given-names>A. R.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>The gut microbiota of rural papua new guineans: composition, diversity patterns, and ecological processes.</article-title> <source><italic>Cell Rep.</italic></source> <volume>11</volume> <fpage>527</fpage>&#x2013;<lpage>538</lpage>. <pub-id pub-id-type="doi">10.1016/j.celrep.2015.03.049</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Milani</surname> <given-names>C.</given-names></name> <name><surname>Lugli</surname> <given-names>G.</given-names></name> <name><surname>Duranti</surname> <given-names>S.</given-names></name> <name><surname>Turroni</surname> <given-names>F.</given-names></name> <name><surname>Mancabelli</surname> <given-names>L.</given-names></name> <name><surname>Ferrario</surname> <given-names>C.</given-names></name><etal/></person-group> (<year>2015a</year>). <article-title>Bifidobacteria exhibit social behavior through carbohydrate resource sharing in the gut.</article-title> <source><italic>Sci. Rep.</italic></source> <volume>5</volume> <issue>15782</issue>. <pub-id pub-id-type="doi">10.1038/srep15782</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Milani</surname> <given-names>C.</given-names></name> <name><surname>Mancabelli</surname> <given-names>L.</given-names></name> <name><surname>Lugli</surname> <given-names>G. A.</given-names></name> <name><surname>Duranti</surname> <given-names>S.</given-names></name> <name><surname>Turroni</surname> <given-names>F.</given-names></name> <name><surname>Ferrario</surname> <given-names>C.</given-names></name><etal/></person-group> (<year>2015b</year>). <article-title>Exploring vertical transmission of bifidobacteria from mother to child.</article-title> <source><italic>Appl. Environ. Microbiol.</italic></source> <volume>81</volume> <fpage>7078</fpage>&#x2013;<lpage>7087</lpage>. <pub-id pub-id-type="doi">10.1128/AEM.02037-15</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Namiki</surname> <given-names>T.</given-names></name> <name><surname>Hachiya</surname> <given-names>T.</given-names></name> <name><surname>Tanaka</surname> <given-names>H.</given-names></name> <name><surname>Sakakibara</surname> <given-names>Y.</given-names></name></person-group> (<year>2012</year>). <article-title>MetaVelvet: an extension of Velvet assembler to de novo metagenome assembly from short sequence reads.</article-title> <source><italic>Nucleic Acids Res.</italic></source> <volume>40</volume>:<issue>e155</issue>. <pub-id pub-id-type="doi">10.1093/nar/gks678</pub-id></citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nicholson</surname> <given-names>J. K.</given-names></name> <name><surname>Holmes</surname> <given-names>E.</given-names></name> <name><surname>Kinross</surname> <given-names>J.</given-names></name> <name><surname>Burcelin</surname> <given-names>R.</given-names></name> <name><surname>Gibson</surname> <given-names>G.</given-names></name> <name><surname>Jia</surname> <given-names>W.</given-names></name><etal/></person-group> (<year>2012</year>). <article-title>Host-gut microbiota metabolic interactions.</article-title> <source><italic>Science</italic></source> <volume>336</volume> <fpage>1262</fpage>&#x2013;<lpage>1267</lpage>. <pub-id pub-id-type="doi">10.1126/science.1223813</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Obregon-Tito</surname> <given-names>A. J.</given-names></name> <name><surname>Tito</surname> <given-names>R. Y.</given-names></name> <name><surname>Metcalf</surname> <given-names>J.</given-names></name> <name><surname>Sankaranarayanan</surname> <given-names>K.</given-names></name> <name><surname>Clemente</surname> <given-names>J. C.</given-names></name> <name><surname>Ursell</surname> <given-names>L. K.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>Subsistence strategies in traditional societies distinguish gut microbiomes.</article-title> <source><italic>Nat. Commun.</italic></source> <volume>6</volume>:<issue>6505</issue>. <pub-id pub-id-type="doi">10.1038/ncomms7505</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Quercia</surname> <given-names>S.</given-names></name> <name><surname>Candela</surname> <given-names>M.</given-names></name> <name><surname>Giuliani</surname> <given-names>C.</given-names></name> <name><surname>Turroni</surname> <given-names>S.</given-names></name> <name><surname>Luiselli</surname> <given-names>D.</given-names></name> <name><surname>Rampelli</surname> <given-names>S.</given-names></name><etal/></person-group> (<year>2014</year>). <article-title>From lifetime to evolution: timescales of human gut microbiota adaptation.</article-title> <source><italic>Front. Microbiol.</italic></source> <volume>5</volume>:<issue>587</issue>. <pub-id pub-id-type="doi">10.3389/fmicb.2014.00587</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rampelli</surname> <given-names>S.</given-names></name> <name><surname>Schnorr</surname> <given-names>S. L.</given-names></name> <name><surname>Consolandi</surname> <given-names>C.</given-names></name> <name><surname>Turroni</surname> <given-names>S.</given-names></name> <name><surname>Severgnini</surname> <given-names>M.</given-names></name> <name><surname>Peano</surname> <given-names>C.</given-names></name><etal/></person-group> (<year>2015</year>). <article-title>Metagenome sequencing of the Hadza hunter-gatherer gut microbiota.</article-title> <source><italic>Curr. Biol.</italic></source> <volume>25</volume> <fpage>1682</fpage>&#x2013;<lpage>1693</lpage>. <pub-id pub-id-type="doi">10.1016/j.cub.2015.04.055</pub-id></citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Rho</surname> <given-names>M.</given-names></name> <name><surname>Tang</surname> <given-names>H.</given-names></name> <name><surname>Ye</surname> <given-names>Y.</given-names></name></person-group> (<year>2010</year>). <article-title>FragGeneScan: predicting genes in short and error-prone reads.</article-title> <source><italic>Nucleic Acids Res.</italic></source> <volume>38</volume>:<issue>e191</issue>. <pub-id pub-id-type="doi">10.1093/nar/gkq747</pub-id></citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schnorr</surname> <given-names>S. L.</given-names></name> <name><surname>Candela</surname> <given-names>M.</given-names></name> <name><surname>Rampelli</surname> <given-names>S.</given-names></name> <name><surname>Centanni</surname> <given-names>M.</given-names></name> <name><surname>Consolandi</surname> <given-names>C.</given-names></name> <name><surname>Basaglia</surname> <given-names>G.</given-names></name><etal/></person-group> (<year>2014</year>). <article-title>Gut microbiome of the Hadza hunter-gatherers.</article-title> <source><italic>Nat. Commun.</italic></source> <volume>5</volume>:<issue>3654</issue>. <pub-id pub-id-type="doi">10.1038/ncomms4654</pub-id></citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Schnorr</surname> <given-names>S. L.</given-names></name> <name><surname>Crittenden</surname> <given-names>A. N.</given-names></name> <name><surname>Venema</surname> <given-names>K.</given-names></name> <name><surname>Marlowe</surname> <given-names>F. W.</given-names></name> <name><surname>Henry</surname> <given-names>A. G.</given-names></name></person-group> (<year>2015</year>). <article-title>Assessing digestibility of Hadza tubers using a dynamic in-vitro model.</article-title> <source><italic>Am. J. Phys. Anthropol.</italic></source> <volume>158</volume> <fpage>371</fpage>&#x2013;<lpage>385</lpage>. <pub-id pub-id-type="doi">10.1002/ajpa.22805</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Segata</surname> <given-names>N.</given-names></name></person-group> (<year>2015</year>). <article-title>Gut microbiome: westernization and the disappearance of intestinal diversity.</article-title> <source><italic>Curr. Biol.</italic></source> <volume>25</volume> <fpage>R611</fpage>&#x2013;<lpage>R613</lpage>. <pub-id pub-id-type="doi">10.1016/j.cub.2015.05.040</pub-id></citation></ref>
<ref id="B28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sonnenburg</surname> <given-names>E. D.</given-names></name> <name><surname>Smits</surname> <given-names>S. A.</given-names></name> <name><surname>Tikhonov</surname> <given-names>M.</given-names></name> <name><surname>Higginbottom</surname> <given-names>S. K.</given-names></name> <name><surname>Wingreen</surname> <given-names>N. S.</given-names></name> <name><surname>Sonnenburg</surname> <given-names>J. L.</given-names></name></person-group> (<year>2016</year>). <article-title>Diet-induced extinctions in the gut microbiota compound over generations.</article-title> <source><italic>Nature</italic></source> <volume>529</volume> <fpage>212</fpage>&#x2013;<lpage>215</lpage>. <pub-id pub-id-type="doi">10.1038/nature16504</pub-id></citation></ref>
<ref id="B29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sonnenburg</surname> <given-names>E. D.</given-names></name> <name><surname>Sonnenburg</surname> <given-names>J. L.</given-names></name></person-group> (<year>2014</year>). <article-title>Starving our microbial self: the deleterious consequences of a diet deficient in microbiota-accessible carbohydrates.</article-title> <source><italic>Cell Metab.</italic></source> <volume>20</volume> <fpage>779</fpage>&#x2013;<lpage>786</lpage>. <pub-id pub-id-type="doi">10.1016/j.cmet.2014.07.003</pub-id></citation></ref>
<ref id="B30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Turroni</surname> <given-names>F.</given-names></name> <name><surname>Peano</surname> <given-names>C.</given-names></name> <name><surname>Pass</surname> <given-names>D. A.</given-names></name> <name><surname>Foroni</surname> <given-names>E.</given-names></name> <name><surname>Severgnini</surname> <given-names>M.</given-names></name> <name><surname>Claesson</surname> <given-names>M. J.</given-names></name><etal/></person-group> (<year>2012</year>). <article-title>Diversity of bifidobacteria within the infant gut microbiota.</article-title> <source><italic>PLoS ONE</italic></source> <volume>7</volume>:<issue>e36957</issue>. <pub-id pub-id-type="doi">10.1371/journal.pone.0036957</pub-id></citation></ref>
<ref id="B31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Turroni</surname> <given-names>F.</given-names></name> <name><surname>Serafini</surname> <given-names>F.</given-names></name> <name><surname>Foroni</surname> <given-names>E.</given-names></name> <name><surname>Duranti</surname> <given-names>S.</given-names></name> <name><surname>O&#x2019;Connell Motherway</surname> <given-names>M.</given-names></name> <name><surname>Taverniti</surname> <given-names>V.</given-names></name><etal/></person-group> (<year>2013</year>). <article-title>Role of sortase-dependent pili of Bifidobacterium bifidum PRL2010 in modulating bacterium-host interactions.</article-title> <source><italic>Proc. Natl. Acad. Sci. U.S.A.</italic></source> <volume>110</volume> <fpage>11151</fpage>&#x2013;<lpage>11156</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1303897110</pub-id></citation></ref>
<ref id="B32"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Warinner</surname> <given-names>C.</given-names></name> <name><surname>Speller</surname> <given-names>C.</given-names></name> <name><surname>Collins</surname> <given-names>M. J.</given-names></name> <name><surname>Lewis</surname> <given-names>C. M.</given-names></name></person-group> (<year>2015</year>). <article-title>Ancient human microbiomes.</article-title> <source><italic>J. Hum. Evol.</italic></source> <volume>79</volume> <fpage>125</fpage>&#x2013;<lpage>136</lpage>. <pub-id pub-id-type="doi">10.1016/j.jhevol.2014.10.016</pub-id></citation></ref>
<ref id="B33"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yatsunenko</surname> <given-names>T.</given-names></name> <name><surname>Rey</surname> <given-names>F. E.</given-names></name> <name><surname>Manary</surname> <given-names>M. J.</given-names></name> <name><surname>Trehan</surname> <given-names>I.</given-names></name> <name><surname>Dominguez-Bello</surname> <given-names>M. G.</given-names></name> <name><surname>Contreras</surname> <given-names>M.</given-names></name><etal/></person-group> (<year>2012</year>). <article-title>Human gut microbiome viewed across age and geography.</article-title> <source><italic>Nature</italic></source> <volume>486</volume> <fpage>222</fpage>&#x2013;<lpage>227</lpage>. <pub-id pub-id-type="doi">10.1038/nature11053</pub-id></citation></ref>
<ref id="B34"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yin</surname> <given-names>Y.</given-names></name> <name><surname>Mao</surname> <given-names>X.</given-names></name> <name><surname>Yang</surname> <given-names>J.</given-names></name> <name><surname>Chen</surname> <given-names>X.</given-names></name> <name><surname>Mao</surname> <given-names>F.</given-names></name> <name><surname>Xu</surname> <given-names>Y.</given-names></name></person-group> (<year>2012</year>). <article-title>dbCAN: a web resource for automated carbohydrate-active enzyme annotation.</article-title> <source><italic>Nucleic Acids Res.</italic></source> <volume>40</volume> <fpage>W445</fpage>&#x2013;<lpage>W451</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gks479</pub-id></citation></ref>
</ref-list>
<fn-group>
<fn id="fn01"><label>1</label><p><ext-link ext-link-type="uri" xlink:href="http://greengenes.secondgenome.com/downloads">http://greengenes.secondgenome.com/downloads</ext-link></p></fn>
<fn id="fn02"><label>2</label><p><ext-link ext-link-type="uri" xlink:href="https://github.com/carden24/Bioinformatics_scripts/blob/master/hmmscan-parser.sh">https://github.com/carden24/Bioinformatics_scripts/blob/master/hmmscan-parser.sh</ext-link></p></fn>
<fn id="fn03"><label>3</label><p><ext-link ext-link-type="uri" xlink:href="https://stat.ethz.ch/R-manual/R-devel/library/stats/html/00Index.html">https://stat.ethz.ch/R-manual/R-devel/library/stats/html/00Index.html</ext-link></p></fn>
</fn-group>
</back>
</article>