Heritability and genomics of gene expression in peripheral blood.

We assessed gene expression profiles in 2,752 twins, using a classic twin design to quantify expression heritability and quantitative trait loci (eQTLs) in peripheral blood. The most highly heritable genes (∼777) were grouped into distinct expression clusters, enriched in gene-poor regions, associated with specific gene function or ontology classes, and strongly associated with disease designation. The design enabled a comparison of twin-based heritability to estimates based on dizygotic identity-by-descent sharing and distant genetic relatedness. Consideration of sampling variation suggests that previous heritability estimates have been upwardly biased. Genotyping of 2,494 twins enabled powerful identification of eQTLs, which we further examined in a replication set of 1,895 unrelated subjects. A large number of non-redundant local eQTLs (6,756) met replication criteria, whereas a relatively small number of distant eQTLs (165) met quality control and replication standards. Our results provide a new resource toward understanding the genetic control of transcription.

Quantification of the relative roles of niche and neutral processes in structuring gastrointestinal microbiomes.

The theoretical description of the forces that shape ecological communities focuses around two classes of models. In niche theory, deterministic interactions between species, individuals, and the environment are considered the dominant factor, whereas in neutral theory, stochastic forces, such as demographic noise, speciation, and immigration, are dominant. Species abundance distributions predicted by the two classes of theory are difficult to distinguish empirically, making it problematic to deduce ecological dynamics from typical measures of diversity and community structure. Here, we show that the fusion of species abundance data with genome-derived measures of evolutionary distance can provide a clear indication of ecological dynamics, capable of quantifying the relative roles played by niche and neutral forces. We apply this technique to six gastrointestinal microbiomes drawn from three different domesticated vertebrates, using high-resolution surveys of microbial species abundance obtained from carefully curated deep 16S rRNA hypervariable tag sequencing data. Although the species abundance patterns are seemingly well fit by the neutral theory of metacommunity assembly, we show that this theory cannot account for the evolutionary patterns in the genomic data; moreover, our analyses strongly suggest that these microbiomes have, in fact, been assembled through processes that involve a significant nonneutral (niche) contribution. Our results demonstrate that high-resolution genomics can remove the ambiguities of process inference inherent in classic ecological measures and permits quantification of the forces shaping complex microbial communities.

Phylogenetic characterization of fecal microbial communities of dogs fed diets with or without supplemental dietary fiber using 454 pyrosequencing.

BACKGROUND: Dogs suffer from many of the same maladies as humans that may be affected by the gut microbiome, but knowledge of the canine microbiome is incomplete. This work aimed to use 16S rDNA tag pyrosequencing to phylogenetically characterize hindgut microbiome in dogs and determine how consumption of dietary fiber affects community structure. PRINCIPAL FINDINGS: Six healthy adult dogs were used in a crossover design. A control diet without supplemental fiber and a beet pulp-supplemented (7.5%) diet were fed. Fecal DNA was extracted and the V3 hypervariable region of the microbial 16S rDNA gene amplified using primers suitable for 454-pyrosequencing. Microbial diversity was assessed on random 2000-sequence subsamples of individual and pooled DNA samples by diet. Our dataset comprised 77,771 reads with an average length of 141 nt. Individual samples contained approximately 129 OTU, with Fusobacteria (23-40% of reads), Firmicutes (14-28% of reads) and Bacteroidetes (31-34% of reads) being co-dominant phyla. Feeding dietary fiber generally decreased Fusobacteria and increased Firmicutes, but these changes were not equally apparent in all dogs. UniFrac analysis revealed that structure of the gut microbiome was affected by diet and Firmicutes appeared to play a strong role in by-diet clustering. CONCLUSIONS: Our data suggest three co-dominant bacterial phyla in the canine hindgut. Furthermore, a relatively small amount of dietary fiber changed the structure of the gut microbiome detectably. Our data are among the first to characterize the healthy canine gut microbiome using pyrosequencing and provide a basis for studies focused on devising dietary interventions for microbiome-associated diseases.

Robust computational analysis of rRNA hypervariable tag datasets.

Next-generation DNA sequencing is increasingly being utilized to probe microbial communities, such as gastrointestinal microbiomes, where it is important to be able to quantify measures of abundance and diversity. The fragmented nature of the 16S rRNA datasets obtained, coupled with their unprecedented size, has led to the recognition that the results of such analyses are potentially contaminated by a variety of artifacts, both experimental and computational. Here we quantify how multiple alignment and clustering errors contribute to overestimates of abundance and diversity, reflected by incorrect OTU assignment, corrupted phylogenies, inaccurate species diversity estimators, and rank abundance distribution functions. We show that straightforward procedural optimizations, combining preexisting tools, are effective in handling large (10(5)-10(6)) 16S rRNA datasets, and we describe metrics to measure the effectiveness and quality of the estimators obtained. We introduce two metrics to ascertain the quality of clustering of pyrosequenced rRNA data, and show that complete linkage clustering greatly outperforms other widely used methods.

Comparative metagenomics reveals host specific metavirulomes and horizontal gene transfer elements in the chicken cecum microbiome.

BACKGROUND: The complex microbiome of the ceca of chickens plays an important role in nutrient utilization, growth and well-being of these animals. Since we have a very limited understanding of the capabilities of most species present in the cecum, we investigated the role of the microbiome by comparative analyses of both the microbial community structure and functional gene content using random sample pyrosequencing. The overall goal of this study was to characterize the chicken cecal microbiome using a pathogen-free chicken and one that had been challenged with Campylobacter jejuni. METHODOLOGY/PRINCIPAL FINDINGS: Comparative metagenomic pyrosequencing was used to generate 55,364,266 bases of random sampled pyrosequence data from two chicken cecal samples. SSU rDNA gene tags and environmental gene tags (EGTs) were identified using SEED subsystems-based annotations. The distribution of phylotypes and EGTs detected within each cecal sample were primarily from the Firmicutes, Bacteroidetes and Proteobacteria, consistent with previous SSU rDNA libraries of the chicken cecum. Carbohydrate metabolism and virulence genes are major components of the EGT content of both of these microbiomes. A comparison of the twelve major pathways in the SEED Virulence Subsystem (metavirulome) represented in the chicken cecum, mouse cecum and human fecal microbiomes showed that the metavirulomes differed between these microbiomes and the metavirulomes clustered by host environment. The chicken cecum microbiomes had the broadest range of EGTs within the SEED Conjugative Transposon Subsystem, however the mouse cecum microbiomes showed a greater abundance of EGTs in this subsystem. Gene assemblies (32 contigs) from one microbiome sample were predominately from the Bacteroidetes, and seven of these showed sequence similarity to transposases, whereas the remaining sequences were most similar to those from catabolic gene families. CONCLUSION/SIGNIFICANCE: This analysis has demonstrated that mobile DNA elements are a major functional component of cecal microbiomes, thus contributing to horizontal gene transfer and functional microbiome evolution. Moreover, the metavirulomes of these microbiomes appear to associate by host environment. These data have implications for defining core and variable microbiome content in a host species. Furthermore, this suggests that the evolution of host specific metavirulomes is a contributing factor in disease resistance to zoonotic pathogens.