Emily R. Davenport

Taxonomic classification of bacterial 16S rRNA genes using short sequencing reads: evaluation of effective study designs.

Massively parallel high throughput sequencing technologies allow us to interrogate the microbial composition of biological samples at unprecedented resolution. The typical approach is to perform high-throughout sequencing of 16S rRNA genes, which are then taxonomically classified based on similarity to known sequences in existing databases. Current technologies cause a predicament though, because although they enable deep coverage of samples, they are limited in the length of sequence they can produce. As a result, high-throughout studies of microbial communities often do not sequence the entire 16S rRNA gene. The challenge is to obtain reliable representation of bacterial communities through taxonomic classification of short 16S rRNA gene sequences. In this study we explored properties of different study designs and developed specific recommendations for effective use of short-read sequencing technologies for the purpose of interrogating bacterial communities, with a focus on classification using naïve Bayesian classifiers. To assess precision and coverage of each design, we used a collection of ∼8,500 manually curated 16S rRNA gene sequences from cultured bacteria and a set of over one million bacterial 16S rRNA gene sequences retrieved from environmental samples, respectively. We also tested different configurations of taxonomic classification approaches using short read sequencing data, and provide recommendations for optimal choice of the relevant parameters. We conclude that with a judicious selection of the sequenced region and the corresponding choice of a suitable training set for taxonomic classification, it is possible to explore bacterial communities at great depth using current technologies, with only a minimal loss of taxonomic resolution.

Seasonal Variation in Human Gut Microbiome Composition

The composition of the human gut microbiome is influenced by many environmental factors. Diet is thought to be one of the most important determinants, though we have limited understanding of the extent to which dietary fluctuations alter variation in the gut microbiome between individuals. In this study, we examined variation in gut microbiome composition between winter and summer over the course of one year in 60 members of a founder population, the Hutterites. Because of their communal lifestyle, Hutterite diets are similar across individuals and remarkably stable throughout the year, with the exception that fresh produce is primarily served during the summer and autumn months. Our data indicate that despite overall gut microbiome stability within individuals over time, there are consistent and significant population-wide shifts in microbiome composition across seasons. We found seasonal differences in both (i) the abundance of particular taxa (false discovery rate <0.05), including highly abundant phyla Bacteroidetes and Firmicutes, and (ii) overall gut microbiome diversity (by Shannon diversity; P = 0.001). It is likely that the dietary fluctuations between seasons with respect to produce availability explain, at least in part, these differences in microbiome composition. For example, high levels of produce containing complex carbohydrates consumed during the summer months might explain increased abundance of Bacteroidetes, which contain complex carbohydrate digesters, and decreased levels of Actinobacteria, which have been negatively correlated to fiber content in food questionnaires. Our observations demonstrate the plastic nature of the human gut microbiome in response to variation in diet.

Cross-species comparisons of host genetic associations with the microbiome.

Recent studies in human populations and mouse models reveal notable congruences in gut microbial taxa whose abundances are partly regulated by host genotype. Host genes associating with these taxa are related to diet sensing, metabolism, and immunity. These broad patterns are further validated in similar studies of nonmammalian microbiomes. The next generation of genome-wide association studies will expand the size of the data sets and refine the microbial phenotypes to fully capture these intriguing signatures of host-microbiome coevolution.

Genome-Wide Association Studies of the Human Gut Microbiota

The bacterial composition of the human fecal microbiome is influenced by many lifestyle factors, notably diet. It is less clear, however, what role host genetics plays in dictating the composition of bacteria living in the gut. In this study, we examined the association of ~200K host genotypes with the relative abundance of fecal bacterial taxa in a founder population, the Hutterites, during two seasons (n = 91 summer, n = 93 winter, n = 57 individuals collected in both). These individuals live and eat communally, minimizing variation due to environmental exposures, including diet, which could potentially mask small genetic effects. Using a GWAS approach that takes into account the relatedness between subjects, we identified at least 8 bacterial taxa whose abundances were associated with single nucleotide polymorphisms in the host genome in each season (at genome-wide FDR of 20%). For example, we identified an association between a taxon known to affect obesity (genus Akkermansia) and a variant near PLD1, a gene previously associated with body mass index. Moreover, we replicate a previously reported association from a quantitative trait locus (QTL) mapping study of fecal microbiome abundance in mice (genus Lactococcus, rs3747113, P = 3.13 x 10−7). Finally, based on the significance distribution of the associated microbiome QTLs in our study with respect to chromatin accessibility profiles, we identified tissues in which host genetic variation may be acting to influence bacterial abundance in the gut.

Heritable components of the human fecal microbiome are associated with visceral fat

BackgroundVariation in the human fecal microbiota has previously been associated with body mass index (BMI). Although obesity is a global health burden, the accumulation of abdominal visceral fat is the specific cardio-metabolic disease risk factor. Here, we explore links between the fecal microbiota and abdominal adiposity using body composition as measured by dual-energy X-ray absorptiometry in a large sample of twins from the TwinsUK cohort, comparing fecal 16S rRNA diversity profiles with six adiposity measures.ResultsWe profile six adiposity measures in 3666 twins and estimate their heritability, finding novel evidence for strong genetic effects underlying visceral fat and android/gynoid ratio. We confirm the association of lower diversity of the fecal microbiome with obesity and adiposity measures, and then compare the association between fecal microbial composition and the adiposity phenotypes in a discovery subsample of twins. We identify associations between the relative abundances of fecal microbial operational taxonomic units (OTUs) and abdominal adiposity measures. Most of these results involve visceral fat associations, with the strongest associations between visceral fat and Oscillospira members. Using BMI as a surrogate phenotype, we pursue replication in independent samples from three population-based cohorts including American Gut, Flemish Gut Flora Project and the extended TwinsUK cohort. Meta-analyses across the replication samples indicate that 8 OTUs replicate at a stringent threshold across all cohorts, while 49 OTUs achieve nominal significance in at least one replication sample. Heritability analysis of the adiposity-associated microbial OTUs prompted us to assess host genetic-microbe interactions at obesity-associated human candidate loci. We observe significant associations of adiposity-OTU abundances with host genetic variants in the FHIT, TDRG1 and ELAVL4 genes, suggesting a potential role for host genes to mediate the link between the fecal microbiome and obesity.ConclusionsOur results provide novel insights into the role of the fecal microbiota in cardio-metabolic disease with clear potential for prevention and novel therapies.

Epigenetic modifications are associated with inter-species gene expression variation in primates

BackgroundChanges in gene regulation have long been thought to play an important role in evolution and speciation, especially in primates. Over the past decade, comparative genomic studies have revealed extensive inter-species differences in gene expression levels, yet we know much less about the extent to which regulatory mechanisms differ between species.ResultsTo begin addressing this gap, we perform a comparative epigenetic study in primate lymphoblastoid cell lines, to query the contribution of RNA polymerase II and four histone modifications, H3K4me1, H3K4me3, H3K27ac, and H3K27me3, to inter-species variation in gene expression levels. We find that inter-species differences in mark enrichment near transcription start sites are significantly more often associated with inter-species differences in the corresponding gene expression level than expected by chance alone. Interestingly, we also find that first-order interactions among the five marks, as well as chromatin states, do not markedly contribute to the degree of association between the marks and inter-species variation in gene expression levels, suggesting that the marginal effects of the five marks dominate this contribution.ConclusionsOur observations suggest that epigenetic modifications are substantially associated with changes in gene expression levels among primates and may represent important molecular mechanisms in primate evolution.

A Quick Introduction to Version Control with Git and GitHub

Many scientists write code as part of their research. Just as experiments are logged in laboratory notebooks, it is important to document the code you use for analysis. However, a few key problems can arise when iteratively developing code that make it difficult to document and track which code version was used to create each result. First, you often need to experiment with new ideas, such as adding new features to a script or increasing the speed of a slow step, but you do not want to risk breaking the currently working code. One often-utilized solution is to make a copy of the script before making new edits. However, this can quickly become a problem because it clutters your file system with uninformative filenames, e.g., analysis.sh, analysis_02.sh, analysis_03.sh, etc. It is difficult to remember the differences between the versions of the files and, more importantly, which version you used to produce specific results, especially if you return to the code months later. Second, you will likely share your code with multiple lab mates or collaborators, and they may have suggestions on how to improve it. If you email the code to multiple people, you will have to manually incorporate all the changes each of them sends. Fortunately, software engineers have already developed software to manage these issues: version control. A version control system (VCS) allows you to track the iterative changes you make to your code. Thus, you can experiment with new ideas but always have the option to revert to a specific past version of the code you used to generate particular results. Furthermore, you can record messages as you save each successive version so that you (or anyone else) reviewing the development history of the code is able to understand the rationale for the given edits. It also facilitates collaboration. Using a VCS, your collaborators can make and save changes to the code, and you can automatically incorporate these changes to the main code base. The collaborative aspect is enhanced with the emergence of websites that host version-controlled code. In this quick guide, we introduce you to one VCS, Git (https://git-scm.com), and one online hosting site, GitHub (https://github.com), both of which are currently popular among scientists and programmers in general. More importantly, we hope to convince you that although mastering a given VCS takes time, you can already achieve great benefits by getting started using a few simple commands. Furthermore, not only does using a VCS solve many common problems when writing code, it can also improve the scientific process. By tracking your code development with a VCS and hosting it online, you are performing science that is more transparent, reproducible, and open to collaboration [1,2]. There is no reason this framework needs to be limited only to code; a VCS is well-suited for tracking any plain-text files: manuscripts, electronic lab notebooks, protocols, etc.

Elucidating the role of the host genome in shaping microbiome composition

ABSTRACT A major goal of microbiome research is to identify the factors that determine bacterial composition within and upon a host. Environmental factors are thought to play a large role, such as diet in determining gut microbiome composition and moisture in determining skin microbiome composition. The role of host genetics, however, has been a source of debate in the literature. Recently, we examined the association of host genetics with human gut microbiome composition in the Hutterites, a population that lives and eats communally. We identified heritable bacterial taxa and host genetic loci associated with their abundances. In this addendum, I put these results into a broader context along with other recent studies of microbiome heritability, and synthesize common themes that appear across organisms and tissues, such as the relatively small extent genetics plays compared to environment and the role of host genetic variation in immune response and barrier integrity.

The Relationship Between the Human Genome and Microbiome Comes into View

The bodys microbiome, composed of microbial cells that number in the trillions, is involved in human health and disease in ways that are just starting to emerge. The microbiome is assembled at birth, develops with its host, and is greatly influenced by environmental factors such as diet and other exposures. Recently, a role for human genetic variation has emerged as also influential in accounting for interpersonal differences in microbiomes. Thus, human genes may influence health directly or by promoting a beneficial microbiome. Studies of the heritability of gut microbiotas reveal a subset of microbes whose abundances are partly genetically determined by the host. However, the use of genome-wide association studies (GWASs) to identify human genetic variants associated with microbiome phenotypes has proven challenging. Studies to date are small by GWAS standards, and cross-study comparisons are hampered by differences in analytical approaches. Nevertheless, associations between microbes or microbial genes and human genes have emerged that are consistent between human populations. Most notably, higher levels of beneficial gut bacteria called Bifidobacteria are associated with the human lactase nonpersister genotype, which typically confers lactose intolerance, in several different human populations. It is time for the microbiome to be incorporated into studies that quantify interactions among genotype, environment, and the microbiome in order to predict human disease susceptibility.

ABO antigen and secretor statuses are not associated with gut microbiota composition in 1,500 twins.

BackgroundHost genetics is one of several factors known to shape human gut microbiome composition, however, the physiological processes underlying the heritability are largely unknown. Inter-individual differences in host factors secreted into the gut lumen may lead to variation in microbiome composition. One such factor is the ABO antigen. This molecule is not only expressed on the surface of red blood cells, but is also secreted from mucosal surfaces in individuals containing an intact FUT2 gene (secretors). Previous studies report differences in microbiome composition across ABO and secretor genotypes. However, due to methodological limitations, the specific bacterial taxa involved remain unknown.ResultsHere, we sought to determine the relationship of the microbiota to ABO blood group and secretor status in a large panel of 1503 individuals from a cohort of twins from the United Kingdom. Contrary to previous reports, robust associations between either ABO or secretor phenotypes and gut microbiome composition were not detected. Overall community structure, diversity, and the relative abundances of individual taxa were not significantly associated with ABO or secretor status. Additionally, joint-modeling approaches were unsuccessful in identifying combinations of taxa that were predictive of ABO or secretor status.ConclusionsDespite previous reports, the taxonomic composition of the microbiota does not appear to be strongly associated with ABO or secretor status in 1503 individuals from the United Kingdom. These results highlight the importance of replicating microbiome-associated traits in large, well-powered cohorts to ensure results are robust.