Skip to main content
  • Research article
  • Open access
  • Published:

mtDNA haplogroup and single nucleotide polymorphisms structure human microbiome communities

Abstract

Background

Although our microbial community and genomes (the human microbiome) outnumber our genome by several orders of magnitude, to what extent the human host genetic complement informs the microbiota composition is not clear. The Human Microbiome Project (HMP) Consortium established a unique population-scale framework with which to characterize the relationship of microbial community structure with their human hosts. A wide variety of taxa and metabolic pathways have been shown to be differentially distributed by virtue of race/ethnicity in the HMP. Given that mtDNA haplogroups are the maternally derived ancestral genomic markers and mitochondria’s role as the generator for cellular ATP, characterizing the relationship between human mtDNA genomic variants and microbiome profiles becomes of potential marked biologic and clinical interest.

Results

We leveraged sequencing data from the HMP to investigate the association between microbiome community structures with its own host mtDNA variants. 15 haplogroups and 631 mtDNA nucleotide polymorphisms (mean sequencing depth of 280X on the mitochondria genome) from 89 individuals participating in the HMP were accurately identified. 16S rRNA (V3-V5 region) sequencing generated microbiome taxonomy profiles and whole genome shotgun sequencing generated metabolic profiles from various body sites were treated as traits to conduct association analysis between haplogroups and host clinical metadata through linear regression. The mtSNPs of individuals with European haplogroups were associated with microbiome profiles using PLINK quantitative trait associations with permutation and adjusted for multiple comparisons. We observe that among 139 stool and 59 vaginal posterior fornix samples, several haplogroups show significant association with specific microbiota (q-value < 0.05) as well as their aggregate community structure (Chi-square with Monte Carlo, p < 0.005), which confirmed and expanded previous research on the association of race and ethnicity with microbiome profile. Our results further indicate that mtDNA variations may render different microbiome profiles, possibly through an inflammatory response to different levels of reactive oxygen species activity.

Conclusions

These data provide initial evidence for the association between host ancestral genome with the structure of its microbiome.

Background

Humans are remarkable hosts to microbes, and we have in fact co-evolved as highly plethoric communities. The NIH Roadmap initiative, known as the Human Microbiome Project (HMP) [1], enabled sequence-based comprehensive characterization of the adult human microbiota. Human-associated microorganisms are present in numbers exceeding the quantities of human cells by at least 10-fold beginning in the neonatal period, and the collective microbiome metagenome exceeds our human genome in terms of gene content by orders of magnitude (>150 fold) [24]. We now appreciate that the microbiota are a metabolically and antigenically vibrant diverse community, which may function as mutualists (symbiotically beneficial), commensalists (of neither harm nor benefit), or pathogens (detrimental to the host). However, we do not yet understand how the host genomic content influences its establishment.

The HMP Consortium [25] established just such a population-scale framework with which to characterize the relationship of microbial communities with their human hosts. The signature frameworks of as many as 18 body sites of 242 screened and phenotyped adults from the target population of 300 subjects have been described. To minimize exogenous and environmental exposure that may influence taxonomy abundance, rigorous clinical standards were applied to screen subjects to assure that the cohort were similar in baseline health status [4]. The majority of the subjects were non-vegetarian, non-smokers, non-obese and yet diverged with respect to race/ethnicity, and parental country of origin. Of interest, although no taxa were universally present among all body habitats and individuals, the carriage of metabolic pathways was surprisingly alike, with a greater degree of similarity observed among related race or ethnic groups [2, 3]. These carriage patterns were functionally relevant, and genomic variation in microbial strains (gains, losses, and polymorphisms) underscored inter-individual variation in the microbiome. Taxonomic profiling associating both clades and metabolism with host covariates (namely age, gender, BMI, blood pressure, race and ethnicity, etc) demonstrated that most microbial variations are not well explained by examined clinical covariates other than race/ethnicity [2].

Race and ethnicity exert their effects through innate or genetically determined biologic mechanisms, and have broader implications with relation to socioeconomic status, diet habit, life style, etc. Therefore, it is not surprising to see the strong association of race/ethnicity with the microbiome. However, self-defined race/ethnicity is not always accurate and further complicated by secondary associations of race and ethnicity with diet, birth country, etc. Thus, further investigation down to the molecular level is essential to gain more knowledge on the underlying mechanism of association and to prevent potential misclassification bias.

One recent study has reported on both the commonality and the distinctions in the gut metagenome when compared among children and adults from rural Venezuela, Malawi, and the urban U.S. Of note, the study cohort was comprised of 531 subjects from a limited number of families (151). While relatively few distinctions in the gut microbiome were observed across all cohorts through the first 3 years of life, pronounced differences in the gut microbiome and functional gene repertoires were noted among geographic locations. However, there was no clustering observed among adult Malawians, Amerindians, nor regional U.S. populations. The investigators concluded that the host age/stage of development and geography served as primary determinants of the gut microbiome [6]. However, it bears mention that the potential for shared or divergent host ancestral genomic variation (mitochondrial or nuclear DNA) among their study cohorts was not investigated, and thus, this study may have failed to recognize an independent contribution of the host genome by examining likely homogenous populations (e.g., 34 Malawian and 19 Venezuelan families with a significant number of mono and dizygotic twin pairs). Moreover, when analyzed by the degree of familial relatedness (mono and dizygotic pairs, siblings, and unrelated children and adults) and cohabitation, significant UniFrac distance metrics was observed with greater heterogeneity. Similar studies have not been duplicated in larger genetically divergent population-based cohorts living in the same region with relatively common diets and exogenous exposures [7].

In contrast, there are several lines of evidence suggesting that the host genomic ancestry may structure the microbiome. A large murine advanced intercross line detected a core measurable microbiota (CMM) consisting of 64 conserved taxonomic groups. In this cohort, 13 murine genomic regions and five quantitative trait loci significantly associated with 26 CMM taxa at both the genus and species level of operational taxonomic unit (OTU) projections [8]. With respect to humans, previous studies have largely utilized twin and unrelated sibling pairs as surrogates for host genomic identity [7, 9, 10]. Turnbaugh et al [7] employed 16S rRNA based analysis and observed that the gut microbial community structures of adult monozygotic (haplotype identical) twin pairs had a degree of similarity comparable to dizygotic pairs, and only slightly more similar to their mothers. Of note, there was a roughly comparable degree of covariation between adult monozygotic and dizygotic twin pairs and deviations from this “core” gut microbiome were associated with adult obesity [7, 9, 10]. Total community DNA and RNA was relatively deep sequenced from one set of obese monozygotic twins, and comparative analysis indicated that the majority of species-level projections were shared with relatively significant variations in abundance [10]. These findings expanded upon earlier unrelated infant and dizygotic twin pair studies across the first year of life, which similarly demonstrated that a higher degree of host relatedness was associated with akin gut microbiome profiles [11]. However, these studies were limited to associations with gut microbiota, failed to consider whether the host nuclear or mitochondrial genome (which is solely maternally inherited) might be driving these associations, nor controlled for the potential impact of perinatal co-morbidities known to accompany human twin gestations [12]. This is of potential importance, as twinning in human is an independent risk factor for the development of common perinatal morbidities and mortalities as prematurity, growth discordance, and twin-twin transfusion syndrome, all of which are recognized to render risk of adult obesity [1315]. Despite these potential limitations, such early studies at relatively lower microbial fingerprinting resolution collectively suggest that that the structure of microbial communities is in some part influenced by their human host’s genome.

Human populations can be divided into mtDNA haplogroups based on SNPs scattered throughout the mitochondrial genome, reflecting mutations accumulated by the maternal lineage. Although there is association between ethnic groups and haplogroups, mtDNA variants are most broadly representative of ancient ancestral roots from tens of thousands of years ago. As people migrated to form isolated groups, small changes accumulated over generations to set genetic diversity. Ergo, haplogroups may be regarded as broadly representative of genomic ancestry but are neither the sole nor definitive ancestral tags and may exist in linkage with other nuclear DNA markers [16, 17]. Given that race and ethnicity are delineated by haplogroup, and that mtDNA is more susceptible to DNA damage and acquires mutations across an individual’s lifetime at a higher rate than nuclear DNA, we hypothesized that mtDNA haplogroup and polymorphisms may be associated with variations in the human microbiome [2].

There are additional inherent characteristics of mtDNA, which make it an attractive candidate for host-microbiome association studies. Mitochondria derive from ancestral endosymbiont bacteria and have 16.5 kb circular double-stranded DNA molecules in multiple copies per cell (heteroplasmy). Consistent with their role as the generators of cellular ATP by oxidative phosphorylation, they play crucial roles in energy metabolism and apoptosis [18], and serve as the primary cellular source of reactive oxygen species [19]. The central role of mitochondrial proteins in cellular energy also makes mtDNA an ideal system for rapid human adaptation to new climate and dietary conditions [20]. The mtDNA genome encodes 13 protein-coding genes involved in respiration and oxidative phosphorylation, alongside two rRNAs and a complete set of 22 tRNAs that are important for protein synthesis (http://www.mitomap.org). Likely due to high levels of reactive oxygen species, lack of protective histones, and a limited DNA repair capacity, somatic mtDNA mutations (primarily in the D-Loop) are known to accumulate in individuals over time [21]. However, the majority of mtDNA variants are benign polymorphisms related to SNPs in the nuclear genome. The frequencies of these variants differ among populations and serve as the foundation of haplogrouping. The mutations in mtDNA have been related to a series of diseases related to neurological, muscular or metabolic disorder [22, 23]. Recent population genetics studies also indicate the association of mtDNA variants with complex human diseases, such as Alzheimer’s, Parkinson disease and cancer [24, 25]. Of interest to our studies, mitochondria have been recently suggested to play a pivotal role in the innate immune response [19, 21, 2629].

In this study, we sought to provide initial evidence regarding the contribution of the host genome to its microbiome through robust interrogation of sequencing data and clinical metadata from the HMP [24]. Our aim was to determine whether there exists a significant association between haplogroups and mtDNA variants with microbiome taxonomic abundance and functional profiling.

Results

16S and WGS profile construction

The Human Microbiome Project screened 554 individuals to enroll 300 subjects (149 males, 151 females, mean age 26, mean BMI 24, 20.0% racial minority and 10.7% Hispanic) [4]. A longitudinal sampling strategy yielded 11,174 primary specimens, from which 12,479 DNA samples were submitted to four centers for metagenomic sequencing. While the majority of the samples were targeted for 16S (V3-V5 region) profiling on 454 FLX Titanium platform [2, 3], whole genome shotgun sequencing (WGS) data were generated for a subset of 681 samples. In this analysis, we utilized a robust cohort of samples from 89 individuals, which were sequenced at a single center and retained both 16S and WGS data. The Operational Taxonomic Unit (OTU) table [30] was generated using the HMP 16S pipeline with high stringency approach on V3 to V5 (v35) variable regions [24]. For each body site, OTU counts from samples belonging to individuals in the cohort were extracted to construct the OTU table (Additional file 1). The abundance of each taxonomic level (Phylum, Order, Family, Genus, Species) was binned and calculated based on the OTU table [2].

WGS was performed on the Illumina GAIIx platform. After identification and removal of human reads, microbial sequences were quality filtered and trimmed. The remaining sequences were aligned to protein families, and the abundance of each KEGG pathway and module [31] was generated to reveal the functional activities of each sample.

Host SNP and haplogroup identification

Mitochondrial genomic sequences were extracted from whole genome sequencing of human samples by aligning with Cambridge Reference Sequence (NC_012920) through BWA (Burrows-Wheeler Aligner) [32]. High sequencing depth was achieved on all samples (mean 280X), thereby enabling a high degree of accuracy and precision of variant call. In total, 631 mitochondrial single nucleotide polymorphisms (mtSNP) from across the whole mitochondrial genome were identified by both GATK (Genome Analysis Tool Kit) [33] and SAMtools (Sequences Alignment/Map tool set) [34] as detailed in Methods.

Human mtDNA is characterized by variants, which in turn define haplogroups and polymorphisms. Mitochondria haplogroups are defined on the basis of haplogroup associating mtSNPs corresponding to the original RFLP (Restriction fragment length polymorphism)-defining loci [35]. All common European haplogroups (H, J, K, T, U, V, W, X, I), Asian haplogroups (B, F), African American haplogroups (L2, L3) and Mexican American haplogroups (A,C) were observed at high confidence with correlation to subject self-identified race and ethnicity [36]. Since the majority of enrolled subjects in HMP are Caucasian, the dominant haplogroups in our study are European haplogroups, especially HV and UK, and are in accordance with prior population-based cohorts [37].

Association with haplogroup

In our first pass analysis, we tried to observe the difference of 16s-based microbiome community profiles among haplogroups at various body sites. OTUs from four body sites (stool, posterior fornix, anterior nares and tongue dorsum) representing four distinct body areas (gastrointestinal tract, vagina, skin and oral) were extracted separately to compare the distribution of top family level abundance in relation to haplogroup. To reduce the effect of low counts in less common haplogroups, based on the distance in the accepted haplogroup evolutionary tree [36], individuals were further binned into HV, JT, IWX, UK, A, BF, C and L2L3 groups [38]. As shown in Figure 1, although the overall family abundance at each body site is consistent with anticipated microbial community profiles [2], the abundance varies by virtue of host mtDNA haplogroup. Lactobacillaceae is decidedly the dominant taxa in the vaginal posterior fornix among most of subjects, yet the relative abundance of Lactobacillaceae is lower in haplogroup BF with high variance (mean 0.63, s.d. 0.42) compared with other haplogroups, such as haplogroup UK, which has a high and relatively stable Lactobacillaceae abundance (mean 0.96, s.d. 0.03) (Figure 1). Similarly, stool has a diverse taxonomic profile with familial taxa-projections dominated by the high abundance of Bacteroidaceae, Ruminococcaceae and Lachnospiraceae; however, a clear difference in Bacteroidaceae abundance between BF (mean 0.67, s.d. 0.15) and IWX (mean 0.29, s.d. 0.16) haplogroups is observed (Figure 1). The anterior nares has the highest number of families identified and the tongue dorsum samples are the most ecologically rich, the difference among haplogroups is also observable. There is marked difference of taxonomy abundance among individuals from various haplogroups, which may reflect the influence of maternal lineage on the microbiome.

Figure 1
figure 1

Family level abundance distribution over haplogroups at various body sites. The abundance of top five most abundant families at each body site (posterior fornix, stool, anterior nares, tongue dorsum) was selected to represent the community. Top families represent more than 90% of the total abundance for each body site.

Given our initial observation described above, we next sought to interrogate genus-level association with haplogroups using stool samples. Recent research shows that niche selection, rather than neutral process, drives the assembly of the gut microbiome [39]. We therefore conducted multiple linear regression analyses to derive stool genus level abundance association with haplogroup together with other metadata. Assigning genus abundance as the dependent variable, we modeled host haplogroup and clinical metadata (including age, gender, race, BMI, blood pressure) as independent variables [2]. The linear regression modeling revealed a persistent association of genus abundance by virtue of haplogroups. For example, the IWX haplogroup has the highest abundance of Coprococcus (p value 0.016), while the BF haplogroup has the lowest abundance of Coprococcus. Similarly, the UK haplogroup demonstrated a preferential high abundance of Roseburia (p value 0.00025) and Streptococcus (p value 0.0086); JT haplogroup demonstrated a low abundance of Akkermansia (p value 0.0038) (Figure 2). Gender, BMI, age and blood pressure were neither modifiers of the haplogroup association, nor were they associated with projections of above genus abundance.

Figure 2
figure 2

The distribution of four genus level abundance over haplogroups in stool samples. The Y-axis shows the relative abundance of the genus in corresponding haplogroups. The haplogroup box marked in red is the haplogroup associated with corresponding genus, which is discovered by multiple linear regression using genus level relative abundance as dependent variable and using haplogroup assignment, clinical metadata as independent variables. The number of samples in each haplogroup is parenthetically annotated, A(6), BF(11), C(10), HV(56), IWX(11), JT(18), L2L3(5), UK(22).

Association with microbiome traits

Despite the evident complexity of the microbiome, two studies [40, 41] have described microbial clusters dominated by Bacteroides, Prevotella and Ruminococcus in stool [40] and the relative dominance, absence, or presence of various species of Lactobacillus in vaginal samples [41]. However, the existence of discrete clusters is challenged by recent research [42], which demonstrates community gradients based on taxonomic abundance. As it was not a primary aim of our study to define or refute the existence of discrete clusters, we alternately tested the association of haplogroups with microbiome traits generated from gradients-based taxonomic abundance. We therefore investigated if individuals in our study of 89 subjects displayed identifiable traits based on taxonomic profile, then tested if these microbiome traits associated with the host haplogroup in the cohort.

Since the vaginal microbiome is dominated by Lactobacillus in most individuals, variance is primarily determined by the relative abundance of Lactobacillus spp [41]. Based on recent work on structure, function and diversity of the microbiome [2, 3, 43], posterior fornix taxonomy profiles were constructed using species abundance of Lactobacillus and genus abundance of other microbiota. All subjects could be assigned to one of groups primarily based on the abundance of specific Lactobacillus species, subtype I dominated by L. crispatus, subtype II dominated by L. gasseri, subtype III dominated by L. iners, subtype V dominated by L. jensenii, subtype IV not dominated by Lactobacillus (Figure 3A). These subtypes could then serve as traits for subsequent tests of association. As shown in Figure 3B, the distribution of subtype varies by virtue of host haplogroups (Chi-square test with Monte Carlo test of 2000 replicates, p value 0.001). 60% of individuals in Subtype IV, which interestingly is not dominated by Lactobacillus spp, have non-European haplogroups, while 90% of subtype I and all subtype V individuals are in European haplogroups (Figure 3B).

Figure 3
figure 3

Haplogroup associations with vaginal microbial community profiles. (A) Heatmap of log10-transformed proportion of microbial taxa found in the posterior fornix samples (color key is indicated in the upper left corner). Subtype (trait) I-V assigned based on the species composition and abundance of vaginal microbial are shown on the top of the heatmap. (B) Vertical bar represent the distribution of each subtype I-V across haplogroups. The number of samples in each haplogroup is parenthetically annotated, A(2), BF(6), C(2), HV(26), IWX(4), JT(6), L2L3(4), UK(9).

Although PAM (Partitioning Around Medoids) clustering failed to identify enterotype-like clustering in our cohort, we observed that the variance of the relative abundance of Bacteroides and Prevotella are the highest among all genus level clades, with an observed trade-off between Prevotella and Bacteroides as previously described [6, 42]. Samples were therefore assigned to one of two gut microbiome traits by virtue of their gradient Bacteroides/Prevotella ratio: group 1 with higher Prevotella abundance, and group 2 with higher Bacteroides abundance (Figure 4A). Akin to the vaginal microbiome trait projections, we again observed significant variation by virtue of host haplogroup (Chi-square test with Monte Carlo test of 2000 replicates, p value 0.002). The proportion of group 1 is significantly higher among the IWX haplogroup, and lower among the haplogroup BF (Figure 4B). The proportion of the two trait groups varies by large region-defined haplogroups, with bias by virtue of European haplogroups (HV, IWX, JT, UK) relative to Asian, Mexican American (BF, C) haplogroups. Of interest, we failed to observe a strong association of non-haplogroup defining mtSNPs to any microbiome trait.

Figure 4
figure 4

Association of haplogroups with gut microbial community profiles. (A) Relative abundance contributing to microbiome traits. The number on the X-axis indicates group 1 and 2. (B) Vertical bar represent the distribution of each group across haplogroup. The number of samples in each haplogroup is parenthetically annotated, A(6), BF(11), C(10), HV(56), IWX(11), JT(18), L2L3(5), UK(22).

Association with SNP

To avoid potential confounding with subject-defined population stratification, we alternately stratified the cohort by haplogroup and performed SNP association exclusively on individuals classified into European haplogroups (namely H, J, K, T, U, V, W, X, I, 57 stool samples, 25 posterior fornix samples) [36]; this is as anticipated given that the majority of our subjects self-identified their ethnicity as Caucasian. Adjacent mtSNPs with similar profile could be represented by tagging SNPs instead of genotyping every individual mtSNP [44, 45]. 158 tagging SNPs were identified using Tagger through Haploview [46] on mtSNPs identified in European haplogroups. These tagging mtSNPs are located in both regulatory and coding regions of the mitochondrial genome.

The relative abundance of each clade was treated as a quantitative trait and tested for association with 158 tagging mtSNPs observed in the European cohort. In order to overcome the potential normality of continuous variable and rare allele testing, we adopted PLINK [47] quantitative trait associations with permutation to derive mtSNP associations. In gut samples, A13434G (q value 0.045) and T15784C (q value 0.045) are significantly associated with Eubacterium and Roseburia, which belong to Clostridiale. G16390A (q value 0.003) was observed in robust association with Deltaproteobacteria and Desulfovibrionaceae, which are both in the phylum of Proteobacteria (Table 1) (Figure 5A). In posterior fornix samples, a non-synonymous point mutation in cytochrome b, T14798C (q value 0.065), is associated with the Veillonellaceae family and one of its genus Dialister (Figure 5B). Overall, SNPs in the mtDNA genomic regions of ND5, CYTB and HV1 were in modest association with stool taxonomies. SNPs associated with microbiome profiles derived from the posterior fornix samples were distributed across the mtDNA genomic 12S region, ND3, CYTB, HV1 and HV2 regions.

Table 1 mtSNPs association with microbiome taxonomy profile in stool and posterior fornix samples
Figure 5
figure 5

Specific regions on mtDNA were shown to associate with stool (panel A) and posterior fornix (panel B) taxonomy profile. -log10(pvalue) is used to represent the strength of association between specific mtSNP and clade. The –log10(pvalue) for each mtSNP on each clade were displayed in one Manhattan plot. The X-axis is labeled with mitochondrial gene names corresponding to their relative position on the chromosome. Significant association between mtSNPs and clade were marked in the plot.

The results from whole genome sequencing of metagenomic samples were represented by KEGG pathway and module relative abundance for carriage patterns. Interestingly, we again observed a strong association between host mtSNP and gut carriage patterns (Table 2). Notably, the C3333T (q value 0.026) found in the ND1 coding region (synonymous mutation) exhibited a strong association with microbial metabolic pathways, including Phenylalanine metabolism and Styrene degradation (Table 2). In posterior fornix samples, a second synonymous point mutation T6776C (q value 0.049) in COI was observed in strong association with fatty acid metabolism.

Table 2 mtSNPs association with KEGG pathways and modules in stool and posterior fornix samples

Discussion

By enabling concomitant analysis of host genomic variants with their respective microbiome community profiles, we can begin to address a number of crucial gaps in our understanding. In this study, we observed significant associations among both mtDNA haplogroups and mtSNPs with microbiota taxa alongside their carriage patterns. These studies serve as robust initial population-based interrogations into host genome-microbiome associations, and further clarify previously described microbiome-ethnicity associations at the molecular level.

First, mtSNPs were identified for individuals in the HMP cohort with high confidence. Haplogroups were identified for each individual based on mtSNPs. Given the fundamental role of the mitochondrial genome in cellular metabolism, evidence has already accumulated that different human mtDNA lineages are functionally different [48, 49]. mtDNA haplogroups are also thought to modify mitochondrial function because there is evidence that they are formed by adaptation to thermal environments [50] and subsequent natural selection [20]. It is clear that mitochondria haplogroup is more than a genomic marker. It reflects the ancestral difference in the human genome among the population. In publications from HMP consortium, association of clades and metabolic pathways with host properties, such as age, gender, BMI and race, have been preliminarily examined and reported [2, 3]. A wide variety of taxa and metabolic pathways were differentially distributed by virtue of subject ethnicity at limited body sites. Although age was shown to be associated with several metabolic pathways and one skin clade, most other metadata (BMI, gender, etc) are generally modest and non-significant [2]. Ergo, other factors, such as host genetics, may play roles in shaping the host microbiome profile [2]. We acknowledge as a potential limitation to our study that while we observed an association between certain mtDNA SNPs and haplogroups with community microbial profiles, these associations may be more broadly representative of ancestral differences, which are in linkage with the mitochondrial genome.

Of interest to our data, Ravel et. al.[41] previously reported on the relationship of racial or ethnic background to vaginal bacterial community composition. Using a cohort of self-sampled subjects, they have previously reported on the proportion of vaginal subtypes, which significantly varied among Caucasian, Asian, African American and Latino subjects. Similarly, other investigators with the MetaHIT have shown that stool enterotype 1 is observed in strong association with self-identifying Japanese subjects [40]. Given that mtDNA haplogroups are the maternally derived genomic markers of race/ethnicity, our results have confirmed and extended these observations of others to molecular level.

Nevertheless, having defined mtDNA haplogroups existing in significant association with specific taxa, mtSNPs across mitochondria genome were further tested to identify regions with higher association. We observed several findings of potential relevance.

T14798C is a non-synonymous SNP that encodes an amino acid substitution of phenylalanine (Phe) to leucine (Leu). This mtSNP maps to cytochrome b and was observed in strong association with differential abundance of Dialister in the Veillonellaceae family specifically in the vaginal posterior fornix. This is of potential biological importance as the posterior fornix of the vagina is most proximal to the cervix and uterus and Dialister has been found in both amniotic fluid and placental tissue of women with preterm premature rupture of the amniotic sac [51]. There is long-standing racial and ethnic disparity in the risk of preterm birth [52].

Turning to the gut microbiome, A13434G is a synonymous mtSNP on ND5 and T15784C is a synonymous mutation also on cytochrome b. Although synonymous mutations are usually referred to as “silent”, increasing evidence demonstrates significant effects on transcription, splicing, and mRNA transport or translation, all of which would alter phenotype [53]. Both of the above synonymous mutations were observed in strong association with Eubacterium (belonging to Clostridium cluster IV) and Roseburia (Clostridium cluster XIVa). Clostridial clusters IV and XIVa are highly oxygen-sensitive anaerobes and produce butyrate along the GI tract, which is excreted in feces [54]. Notably, decreased abundance in members of the butyrate-producing Clostridial clusters IV and XIVa have been reported in the gut of inflammatory bowel disease (IBD) patients, another prevalent human disease which displays variable risk by virtue of race and ethnicity [55]. Functionally extending these observations to the microbiota metabolic profile, several mtSNPs mapping to the mitochondrial encoded NADH dehydrogenase genes were observed in significant association with pathways in amino acid metabolism and protein secretion system from stool samples. It has long been debated whether amino acid derived molecules produced by intestinal bacteria affect host health by regulating host immunity and cell function, or alternately by varying microbial composition and community metabolism [56]. Our findings suggest that the host mtDNA variants in key redox pathways inherently define the gut microbiome, which in turn will structure their community and carriage patterns. This is consistent with long standing demonstrations that in microbiota, protein secretion transport effector molecules function from the interior to exterior. Moreover, protein secretion is known to play a central role in modulating the interactions of bacteria with their environments, particularly when symbiotic bacteria interact with host cellular constituents [57].

The primary role of the mitochondria is to produce energy for the cell. Overall, we observed more clades and metabolic pathways to occur in significant association with mtSNPs mapping to genomic loci encoding for electron transport. The mitochondrial respiratory chain, which is comprised of five multi-subunit protein complexes, carries out oxidative phosphorylation. During oxidative phosphorylation, mitochondrial enzyme complexes generate an electrical charge on either side of the inner mitochondrial membrane through transferring of electrons (membrane potential). This potential difference in electrical charge provides the energy for ATP production. NADH dehydrogenase genes encoded by mtDNA comprise complex I, which is responsible for the first step in the electron transport process. Mitochondrial cytochrome b is a subunit of complex III, which passes electrons to cytochrome c, which relays them to complex IV; three of protein subunits of complex IV are encoded by mitochondrial originated cytochrome c oxidase.

Our observed mutations on mitochondrial genes are known to disrupt normal activity of the electron transport chain, which will affect the production of ATP and increase the production of reactive oxygen species (ROS)[58]. ROS have important roles in cell signaling and homeostasis, but excessive amount of ROS could also cause significant damage to cell structures [59]. In human disease studies, oxidative stress is involved in the development of many diseases, including cancer [60], Parkinson’s disease, Alzheimer’s disease [61], and inflammatory bowel disease [62]. In addition to generating cellular energy, mitochondria are also involved in cell differentiation, cell death, as well as the control of cell cycle and growth [63]. It is also well established that the mitochondrial genomes are dynamic structures whose quantity and quality alter in response to cellular oxidative and metabolic demands. Although the mitochondrial genome is small, mtDNA encodes genes essential to perform the above functions and maintain symbiotic host-microbiome relationships.

In sum, the regulation of metabolic function is highly conserved in higher eukaryotes and is essential to the establishment of microbial communities. Although pathogenic mutations on mtDNA result in severe disease, milder mutations may have subtle phenotypic consequences. Ergo, the mtDNA variations among individuals in the HMP cohort may still render different microbiome profiles, possibly through inflammation response to different levels of ROS activity [58, 64]. Although the detailed mechanism awaits further analysis, individuals who have mtSNPs associated with haplogroups discovered in our study may show differences in their capability of energy metabolism and thus affect their own microbiome profile. We speculate that larger population-based and prospective analyses, which would enroll both healthy and disease-afflicted subjects, will provide further robust evidence for our initial associations proposed herein.

Conclusions

We have described our approach to leverage sequencing data from the HMP to demonstrate significant association between the human microbiome with host mtDNA variants (haplogroup and nucleotide polymorphisms). We observe that among stool and vaginal posterior fornix samples, several haplogroups and mtDNA variants show significant association with specific microbiota and their community structure traits. Delving deeply into both mtSNPs and metagenomic-derived carriage patterns, we have further shown that our observed associations between host and microbe are of likely functional relevance. Given the long-standing described role of mitochondria in cellular metabolism and oxidative stress, and emerging data describing its role in innate immunity, our findings may be of likely high significance. In sum, these data provide initial evidence for the host mitochondrial genome influencing the structure of its microbiome, and underscore the capacity for metagenomics to explore host-microbe interactions.

Methods

Subjects

Our initial study consists of 89 subjects with self identified NIH-defined race (74 Caucasians, 7 Asians, 4 African American and 4 Mexican American, 43 female, 46 male) and ethnicity (Hispanic versus non-Hispanic) in the Human Microbiome Project recruited at Houston, Texas and St. Louis, Missouri. Subjects were all healthy individuals and were sampled one to three times at 15 (male) or 18 (female) body habitats following a common sampling protocol [4]. The subjects were between 18-40 (mean 26 s.d. 5) years old to minimize the variability due to growth and aging. The mean of BMI is 24 (s.d. 4). For 16S, there are 139 stool samples from 81 individuals and 59 posterior fornix samples from 37 individuals in the cohort. For WGS, there are 116 stool samples from 79 individuals and 50 posterior fornix samples from 37 individuals available. The HMP study shows that within-subject variation over time is consistently lower than between-subject variation at both taxonomy and metabolic level [2]. Since each individual’s microbial community is stable over time relative to the population as a whole, the microbial clades association of clinical metadata analysis in HMP includes data from all visits. Only the data from first visit was used for mtSNPs association analysis in our study.

The 16S and WGS data used in this study were downloaded from hmpdacc.org. The detailed procedure for data processing is described in the HMP publication. Briefly, raw V3-V5 region 16S sequences were demultiplexed using QIIME [65]. OTU picking (including error correction, chimera checking through QIIME and clustering via UCLUST [66]) was performed by OTUPipe on V3-5 region. Taxonomy was assigned using the RDP classifier 2.2 [67]. For metabolic reconstruction, the HMP Unified Metabolic Analysis Network (HUMAnN) [68] was used to infer KEGG based community functions from WGS reads. The OTU counts, pathway/module abundance for subjects in our study were extracted from summary OTU table and abundance files (HMQCP for 16S, HMMRC for WGS metabolic reconstruction).

Variants identification

Sequence alignment, quality control, and variant calling were performed with BWA (Burrows-Wheeler Aligner) [32], SAMtools [34], Picard and the Genomic Analysis Toolkit (GATK) [33]. In this study, we utilized the small percentage of human reads from whom the microbiome samples were collected. The whole genome sequencing data without removal of human sequences could be obtained from dbGaP (phs000228.v2.p1). Prior to variant calling, mitochondrial genomic sequences were extracted from whole genome sequencing of human microbiome samples by aligning with Cambridge Reference Sequence (NC_012920) through BWA. SAMtools were used to convert, sort, and index the aligned data files. Picard was then used to identify and remove duplicate reads from each lane. Variants were identified with GATK’s variant detection tools. Following base quality recalibration, indel realignment and unified genotyper, the single-nucleotide variants were filtered for the removal of low-quality variant calls with GATK’s VariantFiltrationWalker tool. In addition to GATK, SAMtools was also used to call targeted bases, and any base call that deviates from reference base was regarded as a potential variation.

Haplogroup analysis

Based on published references [23], PhyloTree [69] and the MITOMAP [36] database, genotypes of several mtSNPs were combined to construct the haplogroups. The most common European mtDNA haplogroups include H, HV, I, J, K, T, U, V, W and X. 70 individuals were identified as European haplogroups. Typical haplogroups for Asian, African American, Mexican American and the rest of Caucasian were identified as A, B, C, F and L2. 19 individuals were classified into non-European haplogroups.

Microbiome trait identification

Since the vaginal microbiome is dominated by one or more species of Lactobacillus, the representative sequences of each OTU in the vagina microbiome was aligned to Greengenes (4feb2011 version) [70] through blast to achieve species level assignment. Greengenes maintains a consistent multiple-sequence alignment of both archaeal and bacterial 16S sequences to facilitate this process. Major Lactobacillus species in the vaginal microbiota (L.crispatus, L.iners, Ljensenii and L.gasseri) were further checked for alignments. OTUs with assignment to Lactobacillus species and other genera based on RDP classifier were normalized (OTU relative abundance were combined based on taxonomy assignment. The abundance of Lactobacillus species and other genera abundance sum to 1.) and then clustered through complete linkage hierarchical clustering using R packages. Subjects were assigned to vaginal community memberships according to dendrogram in the heat map (Figure 3A).

The reconstructed OTU table of stool samples was normalized (The OTU table column sum to 1.) and clustered using partitioning around medoids (PAM) clustering over Euclidian distance. The optimal number of clusters is two, which is validated by silhouette validation technique (silhouette score > 0.5). Group 1 has higher abundance of Prevotella, group 2 has higher abundance of Bacteroides. Our assignment of groups is generally consistent with Wu et al[71].

Haplogroup associated mtSNPs

Since groups of mtSNPs were used for identifying haplogroups, it remained a formal possibility that some other non-haplogroup identifying mtSNPs may be also linked to certain specific haplogroup. These haplogroup-associated mtSNPs may be not as informative as non-haplogroup associated SNPs. To identify haplogroup associated SNPs among European, information from two sources was used. Saxena et al. [37] cataloged all common sites in mtDNA (excluding the control region) from 928 Europeans. 144 sites has a frequency of >1% in these individuals. 64 tagging SNPs were identified as being necessary to tag common mitochondrial variation as well as nine haplogroups (H, V, J, T, U, K, I, W, X) with an r2 (squared correlation of the alleles at two SNPs) of 0.8 using Tagger. Based on Table four in Soxena et al., 69 SNPs out of 144 SNPs (19 tagging SNPs out of 64 tagging SNPs) are identified as potential haplogroup associated SNPs. The other source used is phylotree.org, which includes updated comprehensive phylogeny of global human mtDNA variation and haplogroups. SNPs used for assigning haplogroups down to one level lower than the nine haplogroups mentioned above were identified as potential haplogroup associated SNPs. Only SNPs identified from both sources are marked as haplogroup associated SNPs.

Association analysis

Metadata for HMP individuals are available from dbGaP. All OTU counts in the cohort stool or posterior fornix communities were normalized to relative abundance, where the OTU table columns sum to 1, and relative abundance for each clade was inferred based on the RDP annotation (The abundance of each genus in the same sample sum to 1). There are 109 genera identified from HMP stool samples and 51 genera identified from HMP posterior fornix samples. To avoid bias caused by rare genus, genus appears in less 5% of samples were not included in association analysis. Haplogroup association with quantitative trait was done through multiple linear regression, with clade abundance as dependent variable, and each identified haplogroup (with other haplogroups as one group), BMI, gender, race, blood pressure and age as independent variables. The linear regression was conducted in R. The significance of each independent variable was adjusted for multiple comparisons using R package qvalue [72] with a list of p-values resulting from testing on each clade as input and estimated q-value as output. The q-value measures the false discovery rate of the test.

Since several SNPs were correlated with other SNPs, we used Tagger implemented through haploview [46] to find tagging SNPs which have pairwise values r2 > =0.8. Those mtSNPs tagged by other mtSNPs were not included in the association analysis. mtSNP association analysis with each quantitative trait (either arcsin square root transformed clad abundance or pathway/module abundance) is carried through PLINK [47]. Quantitative trait association with adaptive permutation procedure using default parameters. Given the property of our data, permutation test could relax the assumption about normality of continuous phenotypes and dealing well with rare alleles, as well as small sample size. To control for multiple testing, empirical p-value for each SNP from PLINK was adjusted for multiple comparisons through R package qvalue. SNPs with a q-value less than 0.05 are considered significant.

Data availability

The 16S and WGS data used in this study were downloaded from hmpdacc.org (HMQCP for 16S, HMMRC for WGS metabolic reconstruction). Since there are large amount of human sequences in the whole genome sequencing data, the sequencing data without removal of human sequences could be only obtained from dbGaP (phs000228.v2.p1). Metadata for HMP individuals are also available from dbGaP.

Abbreviations

BMI:

Body mass index

CMM:

Core measurable microbiota

FDR:

False discovery rate

GATK:

Genomic analysis toolkit

HMP:

Human microbiome project

mtDNA:

mitochondrial DNA

mtSNP:

mitochondrial single nucleotide polymorphisms

OTU:

Operational taxonomic unit

PAM:

Partitioning around medoids

RFLP:

Restriction fragment length polymorphism

WGS:

Whole-genome shotgun.

References

  1. Peterson J, Garges S, Giovanni M, McInnes P, Wang L, Schloss JA, Bonazzi V, McEwen JE, Wetterstrand KA, Deal C, Baker CC, Di Francesco V, Howcroft TK, Karp RW, Lunsford RD, Wellington CR, Belachew T, Wright M, Giblin C, David H, Mills M, Salomon R, Mullins C, Akolkar B, Begg L, Davis C, Grandison L, Humble M, Khalsa J, Little AR, et al: The NIH Human Microbiome Project. Genome Res. 2009, 19: 2317-2323.

    Article  PubMed Central  PubMed  Google Scholar 

  2. Human Microbiome Project Consortium: Structure, function and diversity of the healthy human microbiome. Nature. 2012, 486: 207-214. 10.1038/nature11234.

    Article  Google Scholar 

  3. Human Microbiome Project Consortium: A framework for human microbiome research. Nature. 2012, 486: 215-221. 10.1038/nature11209.

    Article  Google Scholar 

  4. Aagaard K, Petrosino J, Keitel W, Watson M, Katancik J, Garcia N, Patel S, Cutting M, Madden T, Hamilton H, Harris E, Gevers D, Simone G, McInnes P, Versalovic J: The Human Microbiome Project strategy for comprehensive sampling of the human microbiome and why it matters. FASEB J. 2013, 27: 1012-1022. 10.1096/fj.12-220806.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  5. Faust K, Sathirapongsasuti JF, Izard J, Segata N, Gevers D, Raes J, Huttenhower C: Microbial co-occurrence relationships in the human microbiome. PLoS Comput Biol. 2012, 8: e1002606-10.1371/journal.pcbi.1002606.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  6. Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, Magris M, Hidalgo G, Baldassano RN, Anokhin AP, Heath AC, Warner B, Reeder J, Kuczynski J, Caporaso JG, Lozupone CA, Lauber C, Clemente JC, Knights D, Knight R, Gordon JI: Human gut microbiome viewed across age and geography. Nature. 2012, 486: 222-227.

    CAS  PubMed Central  PubMed  Google Scholar 

  7. Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, Egholm M, Henrissat B, Heath AC, Knight R, Gordon JI: A core gut microbiome in obese and lean twins. Nature. 2009, 457: 480-484. 10.1038/nature07540.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  8. Benson AK, Kelly SA, Legge R, Ma F, Low SJ, Kim J, Zhang M, Oh PL, Nehrenberg D, Hua K, Kachman SD, Moriyama EN, Walter J, Peterson DA, Pomp D: Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors. Proc Natl Acad Sci U S A. 2010, 107: 18933-18938. 10.1073/pnas.1007028107.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  9. Turnbaugh PJ, Gordon JI: The core gut microbiome, energy balance and obesity. J Physiol. 2009, 587: 4153-4158. 10.1113/jphysiol.2009.174136.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  10. Turnbaugh PJ, Quince C, Faith JJ, McHardy AC, Yatsunenko T, Niazi F, Affourtit J, Egholm M, Henrissat B, Knight R, Gordon JI: Organismal, genetic, and transcriptional variation in the deeply sequenced gut microbiomes of identical twins. Proc Natl Acad Sci U S A. 2010, 107: 7503-7508. 10.1073/pnas.1002355107.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO: Development of the human infant intestinal microbiota. PLoS Biol. 2007, 5: e177-10.1371/journal.pbio.0050177.

    Article  PubMed Central  PubMed  Google Scholar 

  12. Bruder CE, Piotrowski A, Gijsbers AA, Andersson R, Erickson S, Diaz de Stahl T, Menzel U, Sandgren J, von Tell D, Poplawski A, Crowley M, Crasto C, Partridge EC, Tiwari H, Allison DB, Komorowski J, van Ommen GJ, Boomsma DI, Pedersen NL, den Dunnen JT, Wirdefeldt K, Dumanski JP: Phenotypically concordant and discordant monozygotic twins display different DNA copy-number-variation profiles. Am J Hum Genet. 2008, 82: 763-771. 10.1016/j.ajhg.2007.12.011.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  13. Lopriore E, Oepkes D, Walther FJ: Neonatal morbidity in twin-twin transfusion syndrome. Early Hum Dev. 2011, 87: 595-599. 10.1016/j.earlhumdev.2011.07.006.

    Article  PubMed  Google Scholar 

  14. Rossi AC, Mullin PM, Chmait RH: Neonatal outcomes of twins according to birth order, presentation and mode of delivery: a systematic review and meta-analysis. BJOG. 2011, 118: 523-532. 10.1111/j.1471-0528.2010.02836.x.

    Article  CAS  PubMed  Google Scholar 

  15. Chauhan SP, Scardo JA, Hayes E, Abuhamad AZ, Berghella V: Twins: prevalence, problems, and preterm births. Am J Obstet Gynecol. 2010, 203: 305-315. 10.1016/j.ajog.2010.04.031.

    Article  PubMed  Google Scholar 

  16. Chen A, Raule N, Chomyn A, Attardi G: Decreased reactive oxygen species production in cells with mitochondrial haplogroups associated with longevity. PLoS ONE. 2012, 7: e46473-10.1371/journal.pone.0046473.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  17. Fernandez-Moreno M, Soto-Hermida A, Oreiro N, Pertega S, Fenandez-Lopez C, Rego-Perez I, Blanco FJ: Mitochondrial haplogroups define two phenotypes of osteoarthritis. Front Physiol. 2012, 3: 129-

    Article  PubMed Central  PubMed  Google Scholar 

  18. Phillips D, Covian R, Aponte AM, Glancy B, Taylor JF, Chess D, Balaban RS: Regulation of oxidative phosphorylation complex activity: effects of tissue-specific metabolic stress within an allometric series and acute changes in workload. Am J Physiol Regul Integr Comp Physiol. 2012, 302: R1034-R1048. 10.1152/ajpregu.00596.2011.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  19. Cook CC, Kim A, Terao S, Gotoh A, Higuchi M: Consumption of oxygen: a mitochondrial-generated progression signal of advanced cancer. Cell Death Dis. 2012, 3: e258-10.1038/cddis.2011.141.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  20. Mishmar D, Ruiz-Pesini E, Golik P, Macaulay V, Clark AG, Hosseini S, Brandon M, Easley K, Chen E, Brown MD, Sukernik RI, Olckers A, Wallace DC: Natural selection shaped regional mtDNA variation in humans. Proc Natl Acad Sci U S A. 2003, 100: 171-176. 10.1073/pnas.0136972100.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  21. Mambo E, Gao X, Cohen Y, Guo Z, Talalay P, Sidransky D: Electrophile and oxidant damage of mitochondrial DNA leading to rapid evolution of homoplasmic mutations. Proc Natl Acad Sci U S A. 2003, 100: 1838-1843. 10.1073/pnas.0437910100.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Escames G, Lopez LC, Garcia JA, Garcia-Corzo L, Ortiz F, Acuna-Castroviejo D: Mitochondrial DNA and inflammatory diseases. Hum Genet. 2012, 131: 161-173. 10.1007/s00439-011-1057-y.

    Article  CAS  PubMed  Google Scholar 

  23. Guo Y, Yang TL, Liu YZ, Shen H, Lei SF, Yu N, Chen J, Xu T, Cheng Y, Tian Q, Yu P, Deng HW: Mitochondria-wide association study of common variants in osteoporosis. Ann Hum Genet. 2011, 75: 569-574. 10.1111/j.1469-1809.2011.00663.x.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  24. Yakes FM, Van Houten B: Mitochondrial DNA damage is more extensive and persists longer than nuclear DNA damage in human cells following oxidative stress. Proc Natl Acad Sci U S A. 1997, 94: 514-519. 10.1073/pnas.94.2.514.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  25. Yang Y, Shou Z, Zhang P, He Q, Xiao H, Xu Y, Li C, Chen J: Mitochondrial DNA haplogroup R predicts survival advantage in severe sepsis in the Han population. Genet Med. 2008, 10: 187-192. 10.1097/GIM.0b013e318163c343.

    Article  PubMed  Google Scholar 

  26. Seth RB, Sun L, Ea CK, Chen ZJ: Identification and characterization of MAVS, a mitochondrial antiviral signaling protein that activates NF-kappaB and IRF 3. Cell. 2005, 122: 669-682. 10.1016/j.cell.2005.08.012.

    Article  CAS  PubMed  Google Scholar 

  27. Sims GP, Rowe DC, Rietdijk ST, Herbst R, Coyle AJ: HMGB1 and RAGE in inflammation and cancer. Annu Rev Immunol. 2010, 28: 367-388. 10.1146/annurev.immunol.021908.132603.

    Article  CAS  PubMed  Google Scholar 

  28. Zhang Q, Raoof M, Chen Y, Sumi Y, Sursal T, Junger W, Brohi K, Itagaki K, Hauser CJ: Circulating mitochondrial DAMPs cause inflammatory responses to injury. Nature. 2010, 464: 104-107. 10.1038/nature08780.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  29. Zhang Q, Itagaki K, Hauser CJ: Mitochondrial DNA is released by shock and activates neutrophils via p38 map kinase. Shock. 2010, 34: 55-59. 10.1097/SHK.0b013e3181cd8c08.

    Article  PubMed  Google Scholar 

  30. Knights D, Kuczynski J, Charlson ES, Zaneveld J, Mozer MC, Collman RG, Bushman FD, Knight R, Kelley ST: Bayesian community-wide culture-independent microbial source tracking. Nat Methods. 2011, 8: 761-763. 10.1038/nmeth.1650.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  31. Kanehisa M, Goto S, Sato Y, Furumichi M, Tanabe M: KEGG for integration and interpretation of large-scale molecular data sets. Nucleic Acids Res. 2012, 40: D109-D114. 10.1093/nar/gkr988.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo MA: The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20: 1297-1303. 10.1101/gr.107524.110.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  34. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The sequence alignment/map format and SAMtools. Bioinformatics. 2009, 25: 2078-2079. 10.1093/bioinformatics/btp352.

    Article  PubMed Central  PubMed  Google Scholar 

  35. Kim W, Yoo TK, Shin DJ, Rho HW, Jin HJ, Kim ET, Bae YS: Mitochondrial DNA haplogroup analysis reveals no association between the common genetic lineages and prostate cancer in the Korean population. PLoS ONE. 2008, 3: e2211-10.1371/journal.pone.0002211.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Ruiz-Pesini E, Lott MT, Procaccio V, Poole JC, Brandon MC, Mishmar D, Yi C, Kreuziger J, Baldi P, Wallace DC: An enhanced MITOMAP with a global mtDNA mutational phylogeny. Nucleic Acids Res. 2007, 35: D823-D828. 10.1093/nar/gkl927.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  37. Saxena R, de Bakker PI, Singer K, Mootha V, Burtt N, Hirschhorn JN, Gaudet D, Isomaa B, Daly MJ, Groop L, Ardlie KG, Altshuler D: Comprehensive association testing of common mitochondrial DNA variation in metabolic disease. Am J Hum Genet. 2006, 79: 54-61. 10.1086/504926.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Autere J, Moilanen JS, Finnila S, Soininen H, Mannermaa A, Hartikainen P, Hallikainen M, Majamaa K: Mitochondrial DNA polymorphisms as risk factors for Parkinson’s disease and Parkinson’s disease dementia. Hum Genet. 2004, 115: 29-35. 10.1007/s00439-004-1123-9.

    Article  CAS  PubMed  Google Scholar 

  39. Jeraldo P, Sipos M, Chia N, Brulc JM, Dhillon AS, Konkel ME, Larson CL, Nelson KE, Qu A, Schook LB, Yang F, White BA, Goldenfeld N: Quantification of the relative roles of niche and neutral processes in structuring gastrointestinal microbiomes. Proc Natl Acad Sci U S A. 2012, 109: 9692-9698. 10.1073/pnas.1206721109.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  40. Arumugam M, Raes J, Pelletier E, Le Paslier D, Yamada T, Mende DR, Fernandes GR, Tap J, Bruls T, Batto JM, Bertalan M, Borruel N, Casellas F, Fernandez L, Gautier L, Hansen T, Hattori M, Hayashi T, Kleerebezem M, Kurokawa K, Leclerc M, Levenez F, Manichanh C, Nielsen HB, Nielsen T, Pons N, Poulain J, Qin J, Sicheritz-Ponten T, Tims S, et al: Enterotypes of the human gut microbiome. Nature. 2011, 473: 174-180. 10.1038/nature09944.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  41. Ravel J, Gajer P, Abdo Z, Schneider GM, Koenig SS, McCulle SL, Karlebach S, Gorle R, Russell J, Tacket CO, Brotman RM, Davis CC, Ault K, Peralta L, Forney LJ: Vaginal microbiome of reproductive-age women. Proc Natl Acad Sci U S A. 2011, 108 (Suppl 1): 4680-4687.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Koren O, Knights D, Gonzalez A, Waldron L, Segata N, Knight R, Huttenhower C, Ley RE: A guide to enterotypes across the human body: meta-analysis of microbial community structures in human microbiome datasets. PLoS Comput Biol. 2013, 9: e1002863-10.1371/journal.pcbi.1002863.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  43. Aagaard K, Riehle K, Ma J, Segata N, Mistretta TA, Coarfa C, Raza S, Rosenbaum S, Van den Veyver I, Milosavljevic A, Gevers D, Huttenhower C, Petrosino J, Versalovic J: A metagenomic approach to characterization of the vaginal microbiome signature in pregnancy. PLoS ONE. 2012, 7: e36466-10.1371/journal.pone.0036466.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  44. McRae AF, Byrne EM, Zhao ZZ, Montgomery GW, Visscher PM: Power and SNP tagging in whole mitochondrial genome association studies. Genome Res. 2008, 18: 911-917. 10.1101/gr.074872.107.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  45. Bai RK, Leal SM, Covarrubias D, Liu A, Wong LJ: Mitochondrial genetic background modifies breast cancer risk. Cancer Res. 2007, 67: 4687-4694. 10.1158/0008-5472.CAN-06-3554.

    Article  CAS  PubMed  Google Scholar 

  46. Barrett JC: Haploview: wisualization and analysis of SNP genotype data. Cold Spring Harb Protoc. 2009, 2009: pdb ip71-

    Article  PubMed  Google Scholar 

  47. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007, 81: 559-575. 10.1086/519795.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  48. Ruiz-Pesini E, Lapena AC, Diez-Sanchez C, Perez-Martos A, Montoya J, Alvarez E, Diaz M, Urries A, Montoro L, Lopez-Perez MJ, Enriquez JA: Human mtDNA haplogroups associated with high or reduced spermatozoa motility. Am J Hum Genet. 2000, 67: 682-696. 10.1086/303040.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  49. Niemi AK, Hervonen A, Hurme M, Karhunen PJ, Jylha M, Majamaa K: Mitochondrial DNA polymorphisms associated with longevity in a Finnish population. Hum Genet. 2003, 112: 29-33. 10.1007/s00439-002-0843-y.

    Article  CAS  PubMed  Google Scholar 

  50. Ruiz-Pesini E, Mishmar D, Brandon M, Procaccio V, Wallace DC: Effects of purifying and adaptive selection on regional variation in human mtDNA. Science. 2004, 303: 223-226. 10.1126/science.1088434.

    Article  CAS  PubMed  Google Scholar 

  51. Evaldson G, Carlstrom G, Lagrelius A, Malmborg AS, Nord CE: Microbiological findings in pregnant women with premature rupture of the membranes. Med Microbiol Immunol. 1980, 168: 283-297. 10.1007/BF02121812.

    Article  CAS  PubMed  Google Scholar 

  52. Ganu RS, Ma J, Aagaard KM: The role of microbial communities in parturition: is there evidence of association with preterm birth and perinatal morbidity and mortality?. Am J Perinatol. 2013, 30: 613-624.

    Article  PubMed  Google Scholar 

  53. Goymer P: Synonymous mutations break their silence. Nat Rev Genet. 2007, 8: 92-

    CAS  Google Scholar 

  54. Louis P, Young P, Holtrop G, Flint HJ: Diversity of human colonic butyrate-producing bacteria revealed by analysis of the butyryl-CoA:acetate CoA-transferase gene. Environ Microbiol. 2010, 12: 304-314. 10.1111/j.1462-2920.2009.02066.x.

    Article  CAS  PubMed  Google Scholar 

  55. Van Immerseel F, Ducatelle R, De Vos M, Boon N, Van De Wiele T, Verbeke K, Rutgeerts P, Sas B, Louis P, Flint HJ: Butyric acid-producing anaerobic bacteria as a novel probiotic treatment approach for inflammatory bowel disease. J Med Microbiol. 2010, 59: 141-143. 10.1099/jmm.0.017541-0.

    Article  PubMed  Google Scholar 

  56. Dai ZL, Wu G, Zhu WY: Amino acid metabolism in intestinal bacteria: links between gut ecology and host health. Front Biosci. 2011, 16: 1768-1786. 10.2741/3820.

    Article  CAS  Google Scholar 

  57. Gerlach RG, Hensel M: Protein secretion systems and adhesins: the molecular armory of Gram-negative pathogens. Int J Med Microbiol. 2007, 297: 401-415. 10.1016/j.ijmm.2007.03.017.

    Article  CAS  PubMed  Google Scholar 

  58. West AP, Shadel GS, Ghosh S: Mitochondria in innate immune responses. Nat Rev Immunol. 2011, 11: 389-402. 10.1038/nri2975.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  59. Kirkinezos IG, Moraes CT: Reactive oxygen species and mitochondrial diseases. Semin Cell Dev Biol. 2001, 12: 449-457. 10.1006/scdb.2001.0282.

    Article  CAS  PubMed  Google Scholar 

  60. Halliwell B: Oxidative stress and cancer: have we moved forward?. Biochem J. 2007, 401: 1-11.

    Article  CAS  PubMed  Google Scholar 

  61. Valko M, Leibfritz D, Moncol J, Cronin MT, Mazur M, Telser J: Free radicals and antioxidants in normal physiological functions and human disease. Int J Biochem Cell Biol. 2007, 39: 44-84. 10.1016/j.biocel.2006.07.001.

    Article  CAS  PubMed  Google Scholar 

  62. Beltran B, Nos P, Dasi F, Iborra M, Bastida G, Martinez M, O’Connor JE, Saez G, Moret I, Ponce J: Mitochondrial dysfunction, persistent oxidative damage, and catalase inhibition in immune cells of naive and treated Crohn’s disease. Inflamm Bowel Dis. 2010, 16: 76-86. 10.1002/ibd.21027.

    Article  PubMed  Google Scholar 

  63. McBride HM, Neuspiel M, Wasiak S: Mitochondria: more than just a powerhouse. Curr Biol. 2006, 16: R551-R560. 10.1016/j.cub.2006.06.054.

    Article  CAS  PubMed  Google Scholar 

  64. Honda K, Littman DR: The microbiome in infectious disease and inflammation. Annu Rev Immunol. 2012, 30: 759-795. 10.1146/annurev-immunol-020711-074937.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  65. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, Fierer N, Pena AG, Goodrich JK, Gordon JI, Huttley GA, Kelley ST, Knights D, Koenig JE, Ley RE, Lozupone CA, McDonald D, Muegge BD, Pirrung M, Reeder J, Sevinsky JR, Turnbaugh PJ, Walters WA, Widmann J, Yatsunenko T, Zaneveld J, Knight R: QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010, 7: 335-336. 10.1038/nmeth.f.303.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  66. Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26: 2460-2461. 10.1093/bioinformatics/btq461.

    Article  CAS  PubMed  Google Scholar 

  67. Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, Kulam-Syed-Mohideen AS, McGarrell DM, Marsh T, Garrity GM, Tiedje JM: The ribosomal database project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009, 37: D141-D145. 10.1093/nar/gkn879.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  68. Abubucker S, Segata N, Goll J, Schubert AM, Izard J, Cantarel BL, Rodriguez-Mueller B, Zucker J, Thiagarajan M, Henrissat B, White O, Kelley ST, Methe B, Schloss PD, Gevers D, Mitreva M, Huttenhower C: Metabolic reconstruction for metagenomic data and its application to the human microbiome. PLoS Comput Biol. 2012, 8: e1002358-10.1371/journal.pcbi.1002358.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  69. van Oven M, Kayser M: Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum Mutat. 2009, 30: E386-E394. 10.1002/humu.20921.

    Article  PubMed  Google Scholar 

  70. DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, Huber T, Dalevi D, Hu P, Andersen GL: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006, 72: 5069-5072. 10.1128/AEM.03006-05.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  71. Wu GD, Chen J, Hoffmann C, Bittinger K, Chen YY, Keilbaugh SA, Bewtra M, Knights D, Walters WA, Knight R, Sinha R, Gilroy E, Gupta K, Baldassano R, Nessel L, Li H, Bushman FD, Lewis JD: Linking long-term dietary patterns with gut microbial enterotypes. Science. 2011, 334: 105-108. 10.1126/science.1208344.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  72. Alan Dabney JDS, Warnes GR: qvalue: Q-value estimation for false discovery rate control. R package version 1.30.0 edition. http://CRAN.R-project.org/package=qvalue

Download references

Acknowledgements

We wish to acknowledge Drs. Alan Harris and Melissa Suter (Baylor College of Medicine), Jianzhong Ma (MD Anderson Cancer Center), and Rob Knight (University of Colorado-Boulder) for their critical review of the manuscript and insightful discussions.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Kjersti Aagaard.

Additional information

Competing interests

The authors declare no competing financial interests.

Authors’ contributions

KA and JV are investigators with the Human Microbiome Project consortium. KA, JM, JV are currently funded through NICHD and/or the Burroughs Welcome Fund Preterm Birth Initiative on microbiome profiling in maternal health, pregnancy, and parturition. KA is a Maternal-Fetal Medicine specialist, and PB is an expert in mitochondrial DNA variation and mtDNA genomic sequencing and assembly with a focus on mitochondrial disorders and metabolism. JM, CC, XQ and AM are bionformaticians with the Bioinformatics Research Lab and the Human Genome Sequencing Center at Baylor College of Medicine. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2013_6011_MOESM1_ESM.xlsx

Additional file 1: OTU tables and KEGG module/pathway tables. This spreadsheet has 8 tabs. The OTU tables of posterior fornix, stool, anterior nares, tongue dorsum and HUMAnN generated KEGG module/pathway activity files for posterior fornix, stool samples used in the study were provided. (XLSX 16 MB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Ma, J., Coarfa, C., Qin, X. et al. mtDNA haplogroup and single nucleotide polymorphisms structure human microbiome communities. BMC Genomics 15, 257 (2014). https://doi.org/10.1186/1471-2164-15-257

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-257

Keywords