BACKGROUND
our body hosts a vast and mostly unexplored microbial world known as the human microbiota  <cit> .
the microbiota is likely our most intimate connection with the environment.
recent investigations have highlighted the integral role these microorganisms play in human physiology, health and disease  <cit> .
in the gastrointestinal tract, the gut microbiota is mostly composed of anaerobic bacteria of the firmicutes and bacteroidetes phyla  <cit> .
these microorganisms are beneficial to the host since they confer resistance to pathogens  <cit> , stimulate the proliferation of the gut epithelium  <cit> , synthesize essential vitamins and regulate fat storage  <cit> .
however, dysbiosis is associated with clinical conditions, such as obesity  <cit> , diabetes  <cit>  and cancer  <cit> .
regarding obesity, it has been shown that dietary changes alter the gut microbiota in a way that causes its metabolic activity to favor energy acquisition from ingested food, contribute with nutrient absorption and facilitate being stocked in adipose tissue  <cit>  through a diversity of mechanisms  <cit> .
it has been demonstrated that an increase caloric intake, either produced by a high-fat diet  <cit>  or by overfeeding in genetically obese mice  <cit> , selects an obesogenic microbiota.
studies in animal models have shown that shifts in the gut microbiota following weight gain occur in a way that causes obese animals to have proportionally less bacteroidetes and more firmicutes than lean animals  <cit> .
in humans, however, evidence is less clear.
ley et al.
<cit>  studied  <dig> obese individuals following different low-calorie diets and found that weight reduction increased the proportion of bacteroidetes and reduced that of firmicutes, eventually reaching the composition of lean subjects.
in contrast, other authors have described modifications in the composition of the gut microbiota with weight gain in different directions  <cit> .
one aspect that must be noted in the human studies is that most of them have focused on americans or europeans  <cit> , giving a limited picture of the human gut microbiome.
it has been established that the composition of the gut microbiota dramatically varies among individuals  <cit>  and populations  <cit>  according to the geographic  <cit>  and ethnic origin  <cit> , diet  <cit> , host genetics  <cit> , age  <cit>  and several other factors  <cit> .
an open question is how these factors interact with bmi and explain discordant results about the composition of the gut microbiota in lean and obese subjects.
to expand our knowledge of the human microbiome and determine the generality of previous observations concerning shifts in the composition of the gut microbiota following weight gain, we describe, for the first time, the gut microbiota of a group of colombian adults using high throughput dna sequencing and compare it with data previously obtained in other populations .
note that colombians differ from europeans, americans and asians in genetic terms, since they constitute an admixed population involving native american, european and african ancestry in variable proportions  <cit>  and have likely been exposed to different environmental conditions, including dietary habits and lifestyle  <cit> .
we first asked whether the composition of the gut microbiota differs with the geographic origin of the host population.
next, we explored how bmi affects the taxonomic composition of the gut microbiota and determined whether shifts in the composition of this bacterial community following weight gain operated at broad phylogenetic scales  or if they were produced by a reduced number of bacterial phylotypes that, eventually, might become targets to modulate or control obesity.
methods
we analyzed five datasets : original data contributed by us from a group of  <dig> colombian volunteers and four publicly available datasets from the usa, europe, japan and south korea .
the latter datasets were chosen because they represent comprehensive data from populations with distinct geographic origins, were directly comparable with the colombian dataset in terms of the target population , used similar methods to characterize the gut microbiota  and sequenced overlapping regions of the 16s gene .
we first describe how the new data from colombians were obtained and then how we retrieved other data.
colombian dataset
we performed a cross-sectional study with apparently healthy adults of both genders from the general population living in medellin, colombia south america.
volunteers fulfilled the following inclusion criteria: bmi ≥ <dig>  kg/m <dig>  were non smokers, had not been diagnosed with gastrointestinal disease, had not consumed antibiotics or antiparasitics in the last four months, had not consumed laxatives in the last two months, were not enrolled in any weight-reduction program, were not consuming weight-loss supplements, consumed less than  <dig>  or  <dig>  drinks of alcohol per week, and did not exercise for more than 10 hours per week.
we enrolled  <dig> volunteers  who fulfilled these criteria.
note that this sample size did not target any statistical power, since there are no previous data on colombians and the results of studies performed on other populations are highly variable and, in many cases, contradictory  <cit> .
even among studies showing the same pattern, the magnitude of differences between lean and obese individuals is very different.
therefore, the choice of one study or another to calculate a sample size would have been totally arbitrary.
since this study constitutes a first attempt to evaluate the statistical variability of the gut microbiota among colombians, it should be considered a pilot study.
nonetheless, its sample size is comparable to that of previous influential studies  <cit> .
ethical approval
the present study was conducted according to the guidelines laid down in the declaration of helsinki.
in addition, it was considered to have minimal risk according to the colombian ministry of health .
all the volunteers were thoroughly informed about the study and procedures by a member of the team.
participants were assured of anonymity and confidentiality.
written informed consent was obtained from all the volunteers before beginning the study.
the institutional ethics committee of the university of la sabana  reviewed the protocol and the consent forms and approved all the procedures described here.
anthropometric evaluation
weight, height and waist circumference were measured with international techniques after training and standardizing evaluators  <cit> .
weight was measured with a tanita digital scale , height with seca mechanical measuring rods  and waist circumference with mabis tape measures .
each measure was evaluated twice and the average of the two measures was reported.
we then calculated bmi  of participants to classify and select them according to three categories: lean , overweight  or obese .
stool sample preparation
each participant collected a fecal sample in a hermetic, sterile recipient provided by the research team.
samples were immediately refrigerated in household freezers and brought to the laboratory within 12 hours, where a homogenized fraction was lyophilized in a labconco  <dig> freeze dry system  at −50°c and 25 × 10−3 μbar during 48 hours, or until complete desiccation.
dna extraction, sequencing and taxonomic identification of bacteria
one gram of each lyophilized sample was diluted in a sterile saline solution for dna extraction.
dna extraction was performed using the qiaamp dna stool mini kit  according to the manufacturer’s instructions, using 200 μl of diluted samples.
dna was eluted from the column with 50 μl of water and diluted according to a final concentration of 20 ng/μl.
dna was quantified using a nanodrop spectrophotometer  and sent to the research & testing laboratory  for sequencing.
dna sequencing of the 16s rdna was performed with the bacterial tag-encoded flx amplicon pyrosequencing  using 28f 5′tttgatcntggctcag and 519r 5′gtnttacngcggckgctg primers to survey the v <dig>  v <dig> and v <dig> variable regions.
initial generation of the sequencing library utilized a one-step pcr with a total of 30 cycles, a mixture of hot start and hot star high fidelity taq polymerases, and amplicons originating and extending from the 28f primer for bacterial diversity.
the btefap utilized the roche  <dig> flx instrument with titanium reagents and titanium procedures.
the average sequencing depth was 10 k reads per assay.
following dna sequencing, all failed sequence reads ; briefly, the signal processing performs a series of normalization, correction and quality filtering steps and outputs the remaining  signals into flowgrams for each read), low quality sequence ends , barcodes and primers were removed, and sequence collections depleted of any non-bacterial rdna sequence and chimeras using b2c <dig>  <cit> .
to determine the identity of bacteria in the remaining reads, dna sequences were filtered , assigned to samples based on their nucleotide barcode, assembled into clusters of operational taxonomic units  based on their sequence similarity using uclust  <cit>  and pynast  <cit> , and queried against the greengenes database  <cit> , 12_ <dig> release, using the rdp classifier  <cit>  implemented in qiime  <dig> .0-dev  <cit> .
sequence identity ≥80%, ≥95% and ≥97% delimited taxonomy at the phylum, genus and species levels, respectively.
although determining exactly how otus should be defined is an active area of research  <cit> , we adhered to these commonly used values for the sake of comparability with previous studies  <cit> .
phylogenetic trees and otu tables were constructed for each dataset with qiime.
the analysis pipeline is provided as additional file 1: figure s <dig>  raw sequences were deposited at the european nucleotide archive .
assembled sequences are available as additional file  <dig>
american, european and asian datasets
we retrieved and analyzed 16s rdna sequences from some previous studies: usa  <cit> , europe  <cit> , south korea  <cit>  and japan  <cit> .
although in all these studies the bmi of volunteers was recorded, in the usa and european datasets only lean, overweight and obese volunteers were recruited; the japanese and korean datasets focused almost exclusively on lean individuals.
for the usa dataset, where the gut microbiota of obese and lean female twins and their mothers was characterized  <cit> , we downloaded final 454-generated v <dig> 16s rdna sequences  and extracted reads from  <dig> twins .
we refrained from analyzing the two twins or the twin-mother couple because relatedness is a source of within-population community similarity  that might exacerbate statistical differences among populations.
in addition, by restricting analyses to unrelated individuals we made all datasets directly comparable.
also, for the sake of comparability, we only analyzed the first fecal sample  of each subject.figure  <dig>
taxonomic profiles of the gut microbiota of colombians and americans.
relative abundance of phylum-level otus.
relative abundance of the most frequent genus-level otus , colored by their respective phylum .
unclassified phylotypes are marked with asterisk.
upper bars = colombians; lower bars = americans.
the european dataset consisted of the subset of  <dig> healthy volunteers  among spanish, french and danish inhabitants whose microbiomes were published by the metahit consortium  <cit> .
for the sake of comparability with the other studies, italians were not analyzed because they were elders.
raw sequences were downloaded from the ncbi trace archive  and complete 16s gene fragments extracted using blastn searches against the greengenes 12_ <dig> database .
the korean dataset consisted of  <dig> lean and four overweight individuals  in which the v <dig>  v <dig> and v <dig> 16s rdna regions were sequenced  <cit> .
originally denoised and filtered 454-generated sequences were kindly provided by dr. young-do nam.
for comparability with the other studies, we analyzed only the first stool sample  of individuals a–f , and the only stool sample of individuals i–t.
for comparability with the other datasets, we did not analyze individuals g and h since they were children six and four years old, respectively.
finally, the japanese dataset consisted of 454-generated v <dig> and v <dig> 16s rdna sequences of  <dig> lean and one overweight adults  that participated in an intervention with probiotics  <cit> , available at the ncbi’s sra database .
for comparability with the other datasets, we only analyzed sequences obtained before the probiotic intervention.
to compare datasets, we extracted the only common 16s rdna region to the five studies  using the v-xtractor  <dig>   <cit> .
v <dig> sequences were assembled into otus, aligned and queried against greengenes 12_ <dig> using the same procedures described above .
statistical analysis
the gut microbiota of each individual in each dataset was first summarized by taxonomic composition to obtain α-diversity estimates.
rarefaction curves were constructed using chao- <dig>  the number of species-level otus and phylogenetic distance using qiime.
we then assessed the β-diversity of the gut microbiota using multivariate unifrac analyses.
unifrac measures differences between microbial communities based on phylogenetic information; its premise is that two microbial communities with a shared evolutionary history share branches on a phylogenetic tree and that the fraction of branch length shared can be quantified and interpreted as the degree of community similarity.
we restricted analyses to unweighted unifrac distances because heterogeneity in sequencing depth between studies.
unweighted distances consider only changes in species composition   <cit> .
unifrac distances were obtained with fast unifrac  <cit>  using rarefied data .
comparisons among populations , bmi categories  and gender  utilized the analysis of similarity  and the adonis function for permutational multivariate analysis of variance implemented in qiime.
next, we tested hypotheses put forward in previous studies concerning shifts in the taxonomic composition of the gut microbiota between lean and obese subjects in more detail.
for this, we performed linear regressions on the proportions  of phylum-level otus using population, bmi, age and gender as independent variables.
in addition, since it has recently been suggested that latitude would be the main underlying factor explaining between-population differences in firmicutes and bacteroidetes  <cit> , we correlated latitude with the proportions of these two phyla using pearson’s r. when comparing populations, analyses were performed on bacterial proportions because total bacterial counts were significantly different among datasets .
since the colombian, usa and european datasets contained lean, overweight and obese individuals, we analyzed them separately to test the effect of bmi on the composition of the gut microbiota in each population independently.
in these cases, we analyzed the proportions as well as the counts of phylum-level otus and controlled for possible confounding factors .
in addition, we performed univariate f-tests and correlation analysis  in these three datasets to investigate the correlations between genus-level otus and bmi.
where necessary, p-values were adjusted for multiple comparisons  <cit> .
in all analyses, bacterial counts were log-transformed and proportions were arcsin-square-root transformed to guarantee the normal distribution of residuals and homoscedasticity, tested using the shapiro-wilk and fligner-killeen tests, respectively.
note that in genus-level analyses, some individuals had no bacterium of a given genus  and logarithmic transformation was impossible.
however, these data were important because they represented extreme values.
rather than removing them, in these analyses we used the transformation log.
general statistical analyses were performed with r  <dig> .
<dig>  <cit> .
RESULTS
some characteristics of the different datasets are shown in table  <dig>  this table indicated that individuals with excess weight tended to be older than lean individuals; although the tendency was not significant, except in the japanese dataset, it justified controlling for age in statistical models.
table  <dig> also showed that, in the colombian dataset, waist circumference increased significantly at higher bmi, indicating central obesity and justifying taking this variable into account in analyses.table  <dig>
general characteristics of the different datasets
dataset – variable
lean
overweight
obese
p
colombia
europe
japan
south korea
usa
data presented as average ± standard deviation; p-values from anova testing differences among lean, overweight and obese subjects.
wc = waist circumference; na = not available.
geographic variability of the gut microbiota
in the new dataset contributed here on colombians, we obtained  <dig>  16s rdna sequences from the stool samples of the  <dig> volunteers.
of these,  <dig>  sequences passed qiime quality filters and were subsequently analyzed.
the minimum/average/maximum sequence counts per individual were  <dig> / <dig> / <dig> , respectively, and the minimum/average/maximum sequence length was 173/318/529 bp, respectively.
these sequences clustered into  <dig>  different species-level otus , of which  <dig>  could be assigned a phylum name and  <dig> a genus name.
note that the remaining sequences were correctly clustered by percentage of identity but were assigned a higher taxonomic rank.
the comparison between the numbers of observed species-level otus and the chao- <dig> estimator in the colombian dataset suggested that, at the depth of our sequencing, we sampled about half the bacterial diversity hosted in the gut of these volunteers.
additional sampling would be needed to capture the remaining diversity, made of species present at very low abundance .
the tendency was rather similar in the other datasets .
this is a common limitation to most bacterial diversity studies and indicates that rare components of the gut microbiota are difficult to detect at the depth of common sequencing.
however, even if rare species make an important contribution to the total gut diversity, dominant species  are expected to be the main determinants of ecosystem processes  <cit> .
it is, therefore, reasonable to focus on dominant species to investigate differences between populations and test previous observations concerning shifts in the gut microbiota following weight gain.
we found that the gut microbiota of colombians was mostly composed of firmicutes  and bacteroidetes , followed by other phyla present in minor frequencies .
however, variation among volunteers in the proportion of these phyla was notorious, with some individuals having up to 97% of their gut microbiota composed of firmicutes and less than 2% of bacteroidetes, and others having 40% of firmicutes and 53% of bacteroidetes .
the remaining datasets had lower proportions of firmicutes and higher proportions of bacteroidetes , but dispersion of data among individuals was equally notorious than in the colombian dataset .
in the japanese, there was a higher proportion of actinobacteria than in the other datasets .table  <dig>
taxonomic composition of the gut microbiota in the different datasets
phylum
colombia
europe
japan
south korea
usa
p
data presented as average ± standard deviation; p-values from anova testing differences among lean, overweight and obese subjects.
wc = waist circumference; na = not available.
the unifrac analysis indicated that the gut microbiota of colombians was significantly different from that of americans, europeans and asians .
indeed, the geographic origin of the population was the most relevant grouping factor in the analysis of the gut microbiota of colombians, americans, europeans, japanese and koreans, above bmi  or gender  .
on the other hand, our results did not support the hypothesis that firmicutes increases and bacteroidetes decreases with latitude  <cit> .
in contrast, in the five datasets analyzed here, we found that the relative abundance of firmicutes decreased with latitude  and that of bacteroidetes increased with latitude  .figure  <dig>
principal correspondence analysis of unifrac distances.
differences in the composition of the gut microbiota according to the geographic origin of the sampled population , nutritional status  and gender .
r <dig> and p-value from permutational multivariate analysis of variance .
composition of the gut microbiota following weight gain
we found that firmicutes tended to be less abundant at a higher bmi in the colombian dataset when controlling for gender, age and waist circumference .
no change was observed for bacteroidetes though  .
a similar result was found in the usa dataset .
in the european dataset there was no change in firmicutes or bacteroidetes with bmi .figure  <dig>
changes in the abundance of phylum-level and genus-level otus with bmi in the colombian dataset.
a-b: phylum-level otus; c-g: genus-level otus.
background color: green = lean; yellow = overweight; red = obese.
pearson’s r from correlation analyses and p-value from linear models.
we then looked in more detail to see which of the most representative phylotypes, binned at 95% sequence identity , changed their abundance with an increasing bmi.
in the colombian dataset,  <dig> different genus-level otus were identified;  <dig> of them occurred at frequencies greater than  <dig> % and, together, represented  <dig> % of the total diversity of the gut bacteria .
in this dataset, an undetermined lachnospiraceae, faecalibacterium and roseburia were predominant among firmicutes, whereas bacteroides and prevotella were the most abundant bacteroidetes .
we detected that five out of the  <dig> most abundant genus-level phylotypes present in this dataset suffered reductions with an increasing bmi: four firmicutes  and one verrucomicrobia  .
the other datasets had lower species richness but similar numbers of the most prevalent phylotypes than the colombian dataset.
in the usa dataset, among the most prevalent genera bacteroides, coprococcus, oscillospira, parabacteroides, undetermined clostridia, clostridiales, rikenellaceae and ruminococcaceae diminished with bmi whereas catenibacterium became more abundant at a higher bmi.
in europeans, bacteroides became less abundant and undetermined veillonellaceae, which was more abundant at a higher bmi .
discussion
geographic variation of the gut microbiota
the gut microbiota is currently recognized as an organ that interacts in a complex way with the human body.
these bacteria play a fundamental role both in maintaining gut health and contributing to several pathologies  <cit> .
recent research has emphasized the relationship between bacterial composition and obesity  <cit> .
however, there is no consensus about what the “typical” gut microbiota of obese and lean subjects would be.
one of the reasons for this is that we have a limited understanding of the extent to which this relationship is affected by factors such as the geographic origin of the surveyed population.
most studies in humans have focused on americans or europeans  <cit>  and few have done tests in populations with different geographic and genetic origins  <cit> .
yet, it has been established that genetic background and geography are some of the most important determinants of the gut bacterial composition  <cit> .
for instance, a study comparing the gut microbiota of subjects from the amazonas of venezuela, rural malawi and usa metropolitan areas found that the origin of the population primarily explains the variation in the composition of this bacterial community  <cit> .
likewise, another recent study, in which the gut microbiota of hazdas, burkinabes, malawians, italians and americans was compared, found that geography was clearly the most important grouping factor  <cit> .
in agreement with this, we here show, using new data in overlooked colombians, that the origin of the population explains more variability in the composition of the gut microbiota than factors such as bmi or gender.
a recent study suggested a link between latitude and the prevalence of firmicutes and bacteroidetes in a sort of bergman’s rule, where populations living in higher latitudes tend to have a larger body mass and relatively more firmicutes and less bacteroidetes than in populations at lower latitudes  <cit> .
in contrast with such a hypothesis, our results on colombians suggest that individuals from this population have a higher proportion of firmicutes and a lower proportion of bacteroidetes than expected according to colombia’s latitude.
an interesting result obtained with the unifrac analysis was that the taxonomic composition of the gut microbiota of colombians and koreans, and europeans and japanese were partially overlapping.
whereas it is hard to impute such resemblance to host genetic similarities, it is tempting to assign it to shared environmental factors, such as macronutrient intake.
according to national health and nutrition surveys, the energy intake of colombians   <cit>  is closer to that of koreans   <cit>  than to intakes of japanese   <cit> , americans   <cit>  or europeans   <cit> .
such lower energy intake in colombians and koreans is due to an average diet lower in total fat , lower protein content  and higher carbohydrate intake .
likewise, fiber intake seems to be higher in koreans   <cit>  and colombians   <cit>  than americans  or japanese   <cit> .
although this is mere speculation and we do not pretend to claim causality with such rough values, it would be interesting to tease apart the effect of diet and geography on the composition of the gut microbiota.
composition of the gut microbiota in lean and obese individuals
several authors have given support to the observation that firmicutes increases and bacteroidetes decreases in obese compared to lean subjects  <cit> .
in one of the most influential studies to date analyzing the gut microbiota of  <dig> individuals  with different bmi by means of  <dig> pyrosequencing of the v <dig>  v <dig>  complete 16s and whole metagenome, turnbaugh et al.
<cit>  found less bacteroidetes in obese subjects than in those who were lean.
however, they did not detect any difference among firmicutes.
the difference between the original publication and our results in the usa dataset  is likely that the two studies performed different analyses  and that turnbaugh et al.
<cit>  analyzed much more data than we did, which gave them greater statistical power.
using a smaller sample  and a different bacterial identification technique , armougom et al.
<cit>  found the same results of turnbaugh et al.
<cit> .
in contrast with these studies, other authors have described shifts in the gut microbiota with bmi in the opposite sense: a higher proportion of bacteroidetes  <cit>  and a lower proportion of firmicutes  <cit>  in individuals with excess weight compared to lean subjects.
furthermore, other studies have detected increases in both phyla  <cit>  or, more commonly, no difference in their abundance with increasing bmi  <cit> .
our results indicated that, similar to previous studies  <cit> , firmicutes and bacteroidetes were the dominant bacterial phyla colonizing the gut of colombians.
these two phyla constituted >95% of the phylotypes detected in this dataset.
nonetheless, differences in their abundance between individuals, which also occur in the other datasets analyzed here, suggest that there are complex genotype-by-environment interactions that contribute to maintain the bacterial community structure in the face of immune, environmental and lifestyle/dietary exposures.
the uniqueness of each individual’s microbial community is a universal feature of the human microbiome  <cit> .
however, results in the colombian dataset did not agree with the observation of increased firmicutes and reduced bacteroidetes in individuals with a higher bmi.
we found less firmicutes in volunteers with a higher bmi, as observed by others  <cit> , and no shift in bacteroidetes  <cit> .
contradictory results between studies on obesity and phylum-level changes on the gut microbiota are common and have deserved explanations.
inspection of studies revealed they are heterogeneous in several aspects.
whereas some of them, including the new data contributed by us in colombians, assessed bacterial diversity using broad rdna surveys and high throughput sequencing  <cit> , others performed analyses based on taxon-specific oligonucleotide probes  <cit> .
the latter techniques are limited by the specificity of the selected probes, which is uncertain in the absence of large rdna surveys that assess the overall diversity within a sample  <cit> .
another methodological issue that affects comparability between studies is the use of different taxonomic databases to classify 16s rdna sequences.
our choice of greengenes was based on the fact that this is a curated, quality-checked database with millions of sequences that has been proved to improve the classification of 75% of the sequences by one or more taxonomic ranks related to the ncbi  <cit> .
sample size is another issue that can contribute to disagreement among studies.
while some of them analyzed as few as nine or  <dig> individuals  <cit>  others sampled  <dig> subjects or more  <cit> .
a higher sample size reduces sampling stochasticity and increases statistical power.
other factors, such as the duration of the fasting period at the moment of sampling  <cit>  or the storage conditions of stool samples prior to dna extraction  <cit> , could also contribute to differences among studies.
however, as suggested above, a more fundamental aspect that profoundly affects comparability among studies is the geographic origin of the sampled population.
populations differ in two domains: genetic  and environmental .
studies in laboratories with animal models usually lack genetic variation and control macro-environmental variables, which might explain why results in obese and lean animals are more consistent than in humans  <cit> .
since in human studies such controls are not possible, it is important to split apart the contributions of geography and bmi  to changes in this bacterial community.
although pioneering studies associated obesity with phylum-level changes in the gut microbiota, studies finding correlations at lower taxonomic levels are becoming more abundant.
ley et al.
<cit>  did not find differences in any particular subgroup of firmicutes or bacteroidetes with obesity, which made them speculate that factors driving shifts in the gut microbiota composition must operate on highly conserved traits shared by a variety of bacteria within these phyla  <cit> .
however, more recent evidence suggested that specific bacteria might play determinant roles in the maintenance of normal weight  <cit> , in the development of obesity  <cit>  or in disease  <cit> .
in this study, we found that a reduced set of genus-level phylotypes was responsible for the reductions at the phylum level with an increasing bmi.
in colombians, the phylotypes that became less abundant in obese subjects were related to degradation of complex carbohydrates  <cit>  and had been found to correlate with normal weight  <cit> .
results in this population suggest that a lower bmi associates with the presence of primary-fiber degraders and that these bacteria impact the energy balance of the host.
they might represent promising avenues to modulate or control obesity in this population.
CONCLUSIONS
studies examining the gut microbiota outside the usa and europe are beginning to be accumulated.
they expand our knowledge of the human microbiome.
this study contributed to this aim by describing, for the first time, the gut microbiota of unstudied colombians.
we showed that the geographic origin of the studied population was a more important factor driving the taxonomic composition of the gut microbiota than bmi or gender.
strategies to modulate or control obesity via intervention of the gut microbiota should take this effect into account.
availability of supporting data
raw sequences of colombians supporting the results of this article are available in the european nucleotide archive , of europeans in the ncbi trace archive , and of americans, japanese and koreans in the ncbi short read archive .
additional files
additional file 1: table s <dig>  some characteristics of the different datasets analyzed in this study.
figure s <dig> - analysis pipeline.
figure s <dig> - rarefaction curves in the different datasets.
figure s <dig> - inter-individual variability of the gut microbiota among colombians.
figure s <dig> - correlations between the relative abundance of firmicutes and bacteroidetes with latitude.
additional file 2:
assembled sequences of the colombian dataset .
additional file 3:
correlation analyses between genus-level otu abundance and bmi for the colombian, american and european datasets.
abbreviations
anosimanalysis of similarity
bmibody mass index
btefapbacterial tag-encoded flx amplicon pyrosequencing
otuoperational taxonomic unit
rdnaribosomal dna
competing interests
at the moment when the research took place, jse, bev and gma were employed by a food company interested in pre- and probiotics.
bk had no conflict of interests.
authors’ contributions
gma, bk and bev conceived and designed the study of colombians, participated in their recruitment and sample collection; jse obtained data from previous studies, performed analyses and drafted the manuscript.
all authors read and approved the final manuscript.
authors’ information
jse is a researcher at vidarium; bk is currently employed by a food company that does research and market probiotic dairy products; bev and gma are currently employed by universities.
we are especially indebted to all the colombian volunteers who agreed to participate in this study.
we thank p. turnbaugh, m. arumugam and y-d. nam for providing data from previous studies, the qiime staff for their valuable help with sequence submission to embl-ebi, and t. yatsunenko and one anonymous reviewer for valuable comments on previous versions of the manuscript.
this study was funded by the vidarium nutrition, health and wellness research center and the university of la sabana.
funders did not have any role in the design, collection, analysis or interpretation of data.

 our body hosts a vast and mostly unexplored microbial world known as the human microbiota. the microbiota is likely our most intimate connection with the environment. recent investigations have highlighted the integral role these microorganisms play in human physiology, health and disease. in the gastrointestinal tract, the gut microbiota is mostly composed of anaerobic bacteria of the firmicutes and bacteroidetes phyla. these microorganisms are beneficial to the host since they confer resistance to pathogens, stimulate the proliferation of the gut epithelium, synthesize essential vitamins and regulate fat storage. however, dysbiosis is associated with clinical conditions, such as obesity, diabetes and cancer. regarding obesity, it has been shown that dietary changes alter the gut microbiota in a way that causes its metabolic activity to favor energy acquisition from ingested food, contribute with nutrient absorption and facilitate being stocked in adipose tissue through a diversity of mechanisms. it has been demonstrated that an increase caloric intake, either produced by a high-fat diet or by overfeeding in genetically obese mice, selects an obesogenic microbiota. studies in animal models have shown that shifts in the gut microbiota following weight gain occur in a way that causes obese animals to have proportionally less bacteroidetes and more firmicutes than lean animals. in humans, however, evidence is less clear. ley et al. studied obese individuals following different low-calorie diets and found that weight reduction increased the proportion of bacteroidetes and reduced that of firmicutes, eventually reaching the composition of lean subjects. in contrast, other authors have described modifications in the composition of the gut microbiota with weight gain in different directions. one aspect that must be noted in the human studies is that most of them have focused on americans or europeans, giving a limited picture of the human gut microbiome. it has been established that the composition of the gut microbiota dramatically varies among individuals and populations according to the geographic and ethnic origin, diet, host genetics, age and several other factors. an open question is how these factors interact with bmi and explain discordant results about the composition of the gut microbiota in lean and obese subjects. to expand our knowledge of the human microbiome and determine the generality of previous observations concerning shifts in the composition of the gut microbiota following weight gain, we describe, for the first time, the gut microbiota of a group of colombian adults using high throughput dna sequencing and compare it with data previously obtained in other populations. note that colombians differ from europeans, americans and asians in genetic terms, since they constitute an admixed population involving native american, european and african ancestry in variable proportions and have likely been exposed to different environmental conditions, including dietary habits and lifestyle. we first asked whether the composition of the gut microbiota differs with the geographic origin of the host population. next, we explored how bmi affects the taxonomic composition of the gut microbiota and determined whether shifts in the composition of this bacterial community following weight gain operated at broad phylogenetic scales or if they were produced by a reduced number of bacterial phylotypes that, eventually, might become targets to modulate or control obesity. methods we analyzed five datasets: original data contributed by us from a group of colombian volunteers and four publicly available datasets from the usa, europe, japan and south korea. the latter datasets were chosen because they represent comprehensive data from populations with distinct geographic origins, were directly comparable with the colombian dataset in terms of the target population, used similar methods to characterize the gut microbiota and sequenced overlapping regions of the 16s gene. we first describe how the new data from colombians were obtained and then how we retrieved other data. colombian dataset we performed a cross-sectional study with apparently healthy adults of both genders from the general population living in medellin, colombia south america. volunteers fulfilled the following inclusion criteria: bmi ≥  kg/m were non smokers, had not been diagnosed with gastrointestinal disease, had not consumed antibiotics or antiparasitics in the last four months, had not consumed laxatives in the last two months, were not enrolled in any weight-reduction program, were not consuming weight-loss supplements, consumed less than or drinks of alcohol per week, and did not exercise for more than 10 hours per week. we enrolled volunteers who fulfilled these criteria. note that this sample size did not target any statistical power, since there are no previous data on colombians and the results of studies performed on other populations are highly variable and, in many cases, contradictory. even among studies showing the same pattern, the magnitude of differences between lean and obese individuals is very different. therefore, the choice of one study or another to calculate a sample size would have been totally arbitrary. since this study constitutes a first attempt to evaluate the statistical variability of the gut microbiota among colombians, it should be considered a pilot study. nonetheless, its sample size is comparable to that of previous influential studies. ethical approval the present study was conducted according to the guidelines laid down in the declaration of helsinki. in addition, it was considered to have minimal risk according to the colombian ministry of health. all the volunteers were thoroughly informed about the study and procedures by a member of the team. participants were assured of anonymity and confidentiality. written informed consent was obtained from all the volunteers before beginning the study. the institutional ethics committee of the university of la sabana reviewed the protocol and the consent forms and approved all the procedures described here. anthropometric evaluation weight, height and waist circumference were measured with international techniques after training and standardizing evaluators. weight was measured with a tanita digital scale, height with seca mechanical measuring rods and waist circumference with mabis tape measures. each measure was evaluated twice and the average of the two measures was reported. we then calculated bmi of participants to classify and select them according to three categories: lean, overweight or obese. stool sample preparation each participant collected a fecal sample in a hermetic, sterile recipient provided by the research team. samples were immediately refrigerated in household freezers and brought to the laboratory within 12 hours, where a homogenized fraction was lyophilized in a labconco freeze dry system at −50°c and 25 × 10−3 μbar during 48 hours, or until complete desiccation. dna extraction, sequencing and taxonomic identification of bacteria one gram of each lyophilized sample was diluted in a sterile saline solution for dna extraction. dna extraction was performed using the qiaamp dna stool mini kit according to the manufacturer’s instructions, using 200 μl of diluted samples. dna was eluted from the column with 50 μl of water and diluted according to a final concentration of 20 ng/μl. dna was quantified using a nanodrop spectrophotometer and sent to the research & testing laboratory for sequencing. dna sequencing of the 16s rdna was performed with the bacterial tag-encoded flx amplicon pyrosequencing using 28f 5′tttgatcntggctcag and 519r 5′gtnttacngcggckgctg primers to survey the v v and v variable regions. initial generation of the sequencing library utilized a one-step pcr with a total of 30 cycles, a mixture of hot start and hot star high fidelity taq polymerases, and amplicons originating and extending from the 28f primer for bacterial diversity. the btefap utilized the roche flx instrument with titanium reagents and titanium procedures. the average sequencing depth was 10 k reads per assay. following dna sequencing, all failed sequence reads; briefly, the signal processing performs a series of normalization, correction and quality filtering steps and outputs the remaining signals into flowgrams for each read), low quality sequence ends, barcodes and primers were removed, and sequence collections depleted of any non-bacterial rdna sequence and chimeras using b2c. to determine the identity of bacteria in the remaining reads, dna sequences were filtered, assigned to samples based on their nucleotide barcode, assembled into clusters of operational taxonomic units based on their sequence similarity using uclust and pynast, and queried against the greengenes database, 12_ release, using the rdp classifier implemented in qiime.0-dev. sequence identity ≥80%, ≥95% and ≥97% delimited taxonomy at the phylum, genus and species levels, respectively. although determining exactly how otus should be defined is an active area of research, we adhered to these commonly used values for the sake of comparability with previous studies. phylogenetic trees and otu tables were constructed for each dataset with qiime. the analysis pipeline is provided as additional file 1: figure s raw sequences were deposited at the european nucleotide archive. assembled sequences are available as additional file american, european and asian datasets we retrieved and analyzed 16s rdna sequences from some previous studies: usa, europe, south korea and japan. although in all these studies the bmi of volunteers was recorded, in the usa and european datasets only lean, overweight and obese volunteers were recruited; the japanese and korean datasets focused almost exclusively on lean individuals. for the usa dataset, where the gut microbiota of obese and lean female twins and their mothers was characterized, we downloaded final 454-generated v 16s rdna sequences and extracted reads from twins. we refrained from analyzing the two twins or the twin-mother couple because relatedness is a source of within-population community similarity that might exacerbate statistical differences among populations. in addition, by restricting analyses to unrelated individuals we made all datasets directly comparable. also, for the sake of comparability, we only analyzed the first fecal sample of each subject.figure taxonomic profiles of the gut microbiota of colombians and americans. relative abundance of phylum-level otus. relative abundance of the most frequent genus-level otus, colored by their respective phylum. unclassified phylotypes are marked with asterisk. upper bars = colombians; lower bars = americans. the european dataset consisted of the subset of healthy volunteers among spanish, french and danish inhabitants whose microbiomes were published by the metahit consortium. for the sake of comparability with the other studies, italians were not analyzed because they were elders. raw sequences were downloaded from the ncbi trace archive and complete 16s gene fragments extracted using blastn searches against the greengenes 12_ database. the korean dataset consisted of lean and four overweight individuals in which the v v and v 16s rdna regions were sequenced. originally denoised and filtered 454-generated sequences were kindly provided by dr. young-do nam. for comparability with the other studies, we analyzed only the first stool sample of individuals a–f, and the only stool sample of individuals i–t. for comparability with the other datasets, we did not analyze individuals g and h since they were children six and four years old, respectively. finally, the japanese dataset consisted of 454-generated v and v 16s rdna sequences of lean and one overweight adults that participated in an intervention with probiotics, available at the ncbi’s sra database. for comparability with the other datasets, we only analyzed sequences obtained before the probiotic intervention. to compare datasets, we extracted the only common 16s rdna region to the five studies using the v-xtractor. v sequences were assembled into otus, aligned and queried against greengenes 12_ using the same procedures described above. statistical analysis the gut microbiota of each individual in each dataset was first summarized by taxonomic composition to obtain α-diversity estimates. rarefaction curves were constructed using chao- the number of species-level otus and phylogenetic distance using qiime. we then assessed the β-diversity of the gut microbiota using multivariate unifrac analyses. unifrac measures differences between microbial communities based on phylogenetic information; its premise is that two microbial communities with a shared evolutionary history share branches on a phylogenetic tree and that the fraction of branch length shared can be quantified and interpreted as the degree of community similarity. we restricted analyses to unweighted unifrac distances because heterogeneity in sequencing depth between studies. unweighted distances consider only changes in species composition. unifrac distances were obtained with fast unifrac using rarefied data. comparisons among populations, bmi categories and gender utilized the analysis of similarity and the adonis function for permutational multivariate analysis of variance implemented in qiime. next, we tested hypotheses put forward in previous studies concerning shifts in the taxonomic composition of the gut microbiota between lean and obese subjects in more detail. for this, we performed linear regressions on the proportions of phylum-level otus using population, bmi, age and gender as independent variables. in addition, since it has recently been suggested that latitude would be the main underlying factor explaining between-population differences in firmicutes and bacteroidetes, we correlated latitude with the proportions of these two phyla using pearson’s r. when comparing populations, analyses were performed on bacterial proportions because total bacterial counts were significantly different among datasets. since the colombian, usa and european datasets contained lean, overweight and obese individuals, we analyzed them separately to test the effect of bmi on the composition of the gut microbiota in each population independently. in these cases, we analyzed the proportions as well as the counts of phylum-level otus and controlled for possible confounding factors. in addition, we performed univariate f-tests and correlation analysis in these three datasets to investigate the correlations between genus-level otus and bmi. where necessary, p-values were adjusted for multiple comparisons. in all analyses, bacterial counts were log-transformed and proportions were arcsin-square-root transformed to guarantee the normal distribution of residuals and homoscedasticity, tested using the shapiro-wilk and fligner-killeen tests, respectively. note that in genus-level analyses, some individuals had no bacterium of a given genus and logarithmic transformation was impossible. however, these data were important because they represented extreme values. rather than removing them, in these analyses we used the transformation log. general statistical analyses were performed with r..  some characteristics of the different datasets are shown in table  this table indicated that individuals with excess weight tended to be older than lean individuals; although the tendency was not significant, except in the japanese dataset, it justified controlling for age in statistical models. table  also showed that, in the colombian dataset, waist circumference increased significantly at higher bmi, indicating central obesity and justifying taking this variable into account in analyses.table general characteristics of the different datasets dataset – variable lean overweight obese p colombia europe japan south korea usa data presented as average ± standard deviation; p-values from anova testing differences among lean, overweight and obese subjects. wc = waist circumference; na = not available. geographic variability of the gut microbiota in the new dataset contributed here on colombians, we obtained 16s rdna sequences from the stool samples of the volunteers. of these, sequences passed qiime quality filters and were subsequently analyzed. the minimum/average/maximum sequence counts per individual were / /, respectively, and the minimum/average/maximum sequence length was 173/318/529 bp, respectively. these sequences clustered into different species-level otus, of which could be assigned a phylum name and a genus name. note that the remaining sequences were correctly clustered by percentage of identity but were assigned a higher taxonomic rank. the comparison between the numbers of observed species-level otus and the chao- estimator in the colombian dataset suggested that, at the depth of our sequencing, we sampled about half the bacterial diversity hosted in the gut of these volunteers. additional sampling would be needed to capture the remaining diversity, made of species present at very low abundance. the tendency was rather similar in the other datasets. this is a common limitation to most bacterial diversity studies and indicates that rare components of the gut microbiota are difficult to detect at the depth of common sequencing. however, even if rare species make an important contribution to the total gut diversity, dominant species are expected to be the main determinants of ecosystem processes. it is, therefore, reasonable to focus on dominant species to investigate differences between populations and test previous observations concerning shifts in the gut microbiota following weight gain. we found that the gut microbiota of colombians was mostly composed of firmicutes and bacteroidetes, followed by other phyla present in minor frequencies. however, variation among volunteers in the proportion of these phyla was notorious, with some individuals having up to 97% of their gut microbiota composed of firmicutes and less than 2% of bacteroidetes, and others having 40% of firmicutes and 53% of bacteroidetes. the remaining datasets had lower proportions of firmicutes and higher proportions of bacteroidetes, but dispersion of data among individuals was equally notorious than in the colombian dataset. in the japanese, there was a higher proportion of actinobacteria than in the other datasets.table taxonomic composition of the gut microbiota in the different datasets phylum colombia europe japan south korea usa p data presented as average ± standard deviation; p-values from anova testing differences among lean, overweight and obese subjects. wc = waist circumference; na = not available. the unifrac analysis indicated that the gut microbiota of colombians was significantly different from that of americans, europeans and asians. indeed, the geographic origin of the population was the most relevant grouping factor in the analysis of the gut microbiota of colombians, americans, europeans, japanese and koreans, above bmi or gender. on the other hand, our results did not support the hypothesis that firmicutes increases and bacteroidetes decreases with latitude. in contrast, in the five datasets analyzed here, we found that the relative abundance of firmicutes decreased with latitude and that of bacteroidetes increased with latitude.figure principal correspondence analysis of unifrac distances. differences in the composition of the gut microbiota according to the geographic origin of the sampled population, nutritional status and gender. r and p-value from permutational multivariate analysis of variance. composition of the gut microbiota following weight gain we found that firmicutes tended to be less abundant at a higher bmi in the colombian dataset when controlling for gender, age and waist circumference. no change was observed for bacteroidetes though. a similar result was found in the usa dataset. in the european dataset there was no change in firmicutes or bacteroidetes with bmi.figure changes in the abundance of phylum-level and genus-level otus with bmi in the colombian dataset. a-b: phylum-level otus; c-g: genus-level otus. background color: green = lean; yellow = overweight; red = obese. pearson’s r from correlation analyses and p-value from linear models. we then looked in more detail to see which of the most representative phylotypes, binned at 95% sequence identity, changed their abundance with an increasing bmi. in the colombian dataset, different genus-level otus were identified; of them occurred at frequencies greater than % and, together, represented % of the total diversity of the gut bacteria. in this dataset, an undetermined lachnospiraceae, faecalibacterium and roseburia were predominant among firmicutes, whereas bacteroides and prevotella were the most abundant bacteroidetes. we detected that five out of the most abundant genus-level phylotypes present in this dataset suffered reductions with an increasing bmi: four firmicutes and one verrucomicrobia. the other datasets had lower species richness but similar numbers of the most prevalent phylotypes than the colombian dataset. in the usa dataset, among the most prevalent genera bacteroides, coprococcus, oscillospira, parabacteroides, undetermined clostridia, clostridiales, rikenellaceae and ruminococcaceae diminished with bmi whereas catenibacterium became more abundant at a higher bmi. in europeans, bacteroides became less abundant and undetermined veillonellaceae, which was more abundant at a higher bmi. discussion geographic variation of the gut microbiota the gut microbiota is currently recognized as an organ that interacts in a complex way with the human body. these bacteria play a fundamental role both in maintaining gut health and contributing to several pathologies. recent research has emphasized the relationship between bacterial composition and obesity. however, there is no consensus about what the “typical” gut microbiota of obese and lean subjects would be. one of the reasons for this is that we have a limited understanding of the extent to which this relationship is affected by factors such as the geographic origin of the surveyed population. most studies in humans have focused on americans or europeans and few have done tests in populations with different geographic and genetic origins. yet, it has been established that genetic background and geography are some of the most important determinants of the gut bacterial composition. for instance, a study comparing the gut microbiota of subjects from the amazonas of venezuela, rural malawi and usa metropolitan areas found that the origin of the population primarily explains the variation in the composition of this bacterial community. likewise, another recent study, in which the gut microbiota of hazdas, burkinabes, malawians, italians and americans was compared, found that geography was clearly the most important grouping factor. in agreement with this, we here show, using new data in overlooked colombians, that the origin of the population explains more variability in the composition of the gut microbiota than factors such as bmi or gender. a recent study suggested a link between latitude and the prevalence of firmicutes and bacteroidetes in a sort of bergman’s rule, where populations living in higher latitudes tend to have a larger body mass and relatively more firmicutes and less bacteroidetes than in populations at lower latitudes. in contrast with such a hypothesis, our results on colombians suggest that individuals from this population have a higher proportion of firmicutes and a lower proportion of bacteroidetes than expected according to colombia’s latitude. an interesting result obtained with the unifrac analysis was that the taxonomic composition of the gut microbiota of colombians and koreans, and europeans and japanese were partially overlapping. whereas it is hard to impute such resemblance to host genetic similarities, it is tempting to assign it to shared environmental factors, such as macronutrient intake. according to national health and nutrition surveys, the energy intake of colombians is closer to that of koreans than to intakes of japanese, americans or europeans. such lower energy intake in colombians and koreans is due to an average diet lower in total fat, lower protein content and higher carbohydrate intake. likewise, fiber intake seems to be higher in koreans and colombians than americans or japanese. although this is mere speculation and we do not pretend to claim causality with such rough values, it would be interesting to tease apart the effect of diet and geography on the composition of the gut microbiota. composition of the gut microbiota in lean and obese individuals several authors have given support to the observation that firmicutes increases and bacteroidetes decreases in obese compared to lean subjects. in one of the most influential studies to date analyzing the gut microbiota of individuals with different bmi by means of pyrosequencing of the v v complete 16s and whole metagenome, turnbaugh et al. found less bacteroidetes in obese subjects than in those who were lean. however, they did not detect any difference among firmicutes. the difference between the original publication and our results in the usa dataset is likely that the two studies performed different analyses and that turnbaugh et al. analyzed much more data than we did, which gave them greater statistical power. using a smaller sample and a different bacterial identification technique, armougom et al. found the same results of turnbaugh et al.. in contrast with these studies, other authors have described shifts in the gut microbiota with bmi in the opposite sense: a higher proportion of bacteroidetes and a lower proportion of firmicutes in individuals with excess weight compared to lean subjects. furthermore, other studies have detected increases in both phyla or, more commonly, no difference in their abundance with increasing bmi. our results indicated that, similar to previous studies, firmicutes and bacteroidetes were the dominant bacterial phyla colonizing the gut of colombians. these two phyla constituted >95% of the phylotypes detected in this dataset. nonetheless, differences in their abundance between individuals, which also occur in the other datasets analyzed here, suggest that there are complex genotype-by-environment interactions that contribute to maintain the bacterial community structure in the face of immune, environmental and lifestyle/dietary exposures. the uniqueness of each individual’s microbial community is a universal feature of the human microbiome. however, results in the colombian dataset did not agree with the observation of increased firmicutes and reduced bacteroidetes in individuals with a higher bmi. we found less firmicutes in volunteers with a higher bmi, as observed by others, and no shift in bacteroidetes. contradictory results between studies on obesity and phylum-level changes on the gut microbiota are common and have deserved explanations. inspection of studies revealed they are heterogeneous in several aspects. whereas some of them, including the new data contributed by us in colombians, assessed bacterial diversity using broad rdna surveys and high throughput sequencing, others performed analyses based on taxon-specific oligonucleotide probes. the latter techniques are limited by the specificity of the selected probes, which is uncertain in the absence of large rdna surveys that assess the overall diversity within a sample. another methodological issue that affects comparability between studies is the use of different taxonomic databases to classify 16s rdna sequences. our choice of greengenes was based on the fact that this is a curated, quality-checked database with millions of sequences that has been proved to improve the classification of 75% of the sequences by one or more taxonomic ranks related to the ncbi. sample size is another issue that can contribute to disagreement among studies. while some of them analyzed as few as nine or individuals others sampled subjects or more. a higher sample size reduces sampling stochasticity and increases statistical power. other factors, such as the duration of the fasting period at the moment of sampling or the storage conditions of stool samples prior to dna extraction, could also contribute to differences among studies. however, as suggested above, a more fundamental aspect that profoundly affects comparability among studies is the geographic origin of the sampled population. populations differ in two domains: genetic and environmental. studies in laboratories with animal models usually lack genetic variation and control macro-environmental variables, which might explain why results in obese and lean animals are more consistent than in humans. since in human studies such controls are not possible, it is important to split apart the contributions of geography and bmi to changes in this bacterial community. although pioneering studies associated obesity with phylum-level changes in the gut microbiota, studies finding correlations at lower taxonomic levels are becoming more abundant. ley et al. did not find differences in any particular subgroup of firmicutes or bacteroidetes with obesity, which made them speculate that factors driving shifts in the gut microbiota composition must operate on highly conserved traits shared by a variety of bacteria within these phyla. however, more recent evidence suggested that specific bacteria might play determinant roles in the maintenance of normal weight, in the development of obesity or in disease. in this study, we found that a reduced set of genus-level phylotypes was responsible for the reductions at the phylum level with an increasing bmi. in colombians, the phylotypes that became less abundant in obese subjects were related to degradation of complex carbohydrates and had been found to correlate with normal weight. results in this population suggest that a lower bmi associates with the presence of primary-fiber degraders and that these bacteria impact the energy balance of the host. they might represent promising avenues to modulate or control obesity in this population.  studies examining the gut microbiota outside the usa and europe are beginning to be accumulated. they expand our knowledge of the human microbiome. this study contributed to this aim by describing, for the first time, the gut microbiota of unstudied colombians. we showed that the geographic origin of the studied population was a more important factor driving the taxonomic composition of the gut microbiota than bmi or gender. strategies to modulate or control obesity via intervention of the gut microbiota should take this effect into account. availability of supporting data raw sequences of colombians supporting the results of this article are available in the european nucleotide archive, of europeans in the ncbi trace archive, and of americans, japanese and koreans in the ncbi short read archive. additional files additional file 1: table s some characteristics of the different datasets analyzed in this study. figure s - analysis pipeline. figure s - rarefaction curves in the different datasets. figure s - inter-individual variability of the gut microbiota among colombians. figure s - correlations between the relative abundance of firmicutes and bacteroidetes with latitude. additional file 2: assembled sequences of the colombian dataset. additional file 3: correlation analyses between genus-level otu abundance and bmi for the colombian, american and european datasets. abbreviations anosimanalysis of similarity bmibody mass index btefapbacterial tag-encoded flx amplicon pyrosequencing otuoperational taxonomic unit rdnaribosomal dna competing interests at the moment when the research took place, jse, bev and gma were employed by a food company interested in pre- and probiotics. bk had no conflict of interests. authors’ contributions gma, bk and bev conceived and designed the study of colombians, participated in their recruitment and sample collection; jse obtained data from previous studies, performed analyses and drafted the manuscript. all authors read and approved the final manuscript. authors’ information jse is a researcher at vidarium; bk is currently employed by a food company that does research and market probiotic dairy products; bev and gma are currently employed by universities. we are especially indebted to all the colombian volunteers who agreed to participate in this study. we thank p. turnbaugh, m. arumugam and y-d. nam for providing data from previous studies, the qiime staff for their valuable help with sequence submission to embl-ebi, and t. yatsunenko and one anonymous reviewer for valuable comments on previous versions of the manuscript. this study was funded by the vidarium nutrition, health and wellness research center and the university of la sabana. funders did not have any role in the design, collection, analysis or interpretation of data.
any strategy aiming to modulate or control obesity via manipulation of this bacterial community should consider this effect. a more detailed inspection of the colombian dataset revealed that five fiber-degrading bacteria, including akkermansia, dialister, oscillospira, ruminococcaceae and clostridiales, became less abundant in obese subjects. thus, we describe here, for the first time, the gut microbiota of colombian adults via the pyrosequencing of the 16s ribosomal dna, comparing it with results obtained in americans, europeans, japanese and south koreans, and testing the generality of previous observations concerning changes in firmicutes and bacteroidetes with increasing body mass index. a similar result was found in americans. to determine the generality of previous observations and expand our knowledge of the human gut microbiota, it is important to replicate studies in overlooked populations. the geographic origin of the population explained more variance in the composition of this bacterial community than bmi or gender. we contributed data from unstudied colombians that showed that the geographic origin of the studied population had a greater impact on the composition of the gut microbiota than bmi or gender. the composition of the gut microbiota has recently been associated with health and disease, particularly with obesity. we found that the composition of the gut microbiota of colombians was significantly different from that of americans, europeans and asians.
 the composition of the gut microbiota has recently been associated with health and disease, particularly with obesity. some studies suggested a higher proportion of firmicutes and a lower proportion of bacteroidetes in obese compared to lean people; others found discordant patterns. most studies, however, focused on americans or europeans, giving a limited picture of the gut microbiome. to determine the generality of previous observations and expand our knowledge of the human gut microbiota, it is important to replicate studies in overlooked populations. thus, we describe here, for the first time, the gut microbiota of colombian adults via the pyrosequencing of the 16s ribosomal dna, comparing it with results obtained in americans, europeans, japanese and south koreans, and testing the generality of previous observations concerning changes in firmicutes and bacteroidetes with increasing body mass index.  we found that the composition of the gut microbiota of colombians was significantly different from that of americans, europeans and asians. the geographic origin of the population explained more variance in the composition of this bacterial community than bmi or gender. concerning changes in firmicutes and bacteroidetes with obesity, in colombians we found a tendency in firmicutes to diminish with increasing bmi, whereas no change was observed in bacteroidetes. a similar result was found in americans. a more detailed inspection of the colombian dataset revealed that five fiber-degrading bacteria, including akkermansia, dialister, oscillospira, ruminococcaceae and clostridiales, became less abundant in obese subjects.  we contributed data from unstudied colombians that showed that the geographic origin of the studied population had a greater impact on the composition of the gut microbiota than bmi or gender. any strategy aiming to modulate or control obesity via manipulation of this bacterial community should consider this effect. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords bacterial diversitymicrobiomegeographyfirmicutesbacteroidetesobesitypyrosequencinglatin americacolombiaissue-copyright-statement© the author 2014 