1. Genetics and Genomics
  2. Microbiology and Infectious Disease

A phylogenetic transform enhances analysis of compositional microbiota data

  1. Justin D Silverman
  2. Alex D Washburne
  3. Sayan Mukherjee
  4. Lawrence A David  Is a corresponding author
  1. Duke University, United States
https://doi.org/10.7554/eLife.21887
Share this article
Cite this article
  1. Justin D Silverman
  2. Alex D Washburne
  3. Sayan Mukherjee
  4. Lawrence A David
(2017)
A phylogenetic transform enhances analysis of compositional microbiota data
eLife 6:e21887.
https://doi.org/10.7554/eLife.21887
  2. Download BibTeX
  3. Download .RIS

Abstract

Surveys of microbial communities (microbiota), typically measured as relative abundance of species, have illustrated the importance of these communities in human health and disease. Yet, statistical artifacts commonly plague the analysis of relative abundance data. Here, we introduce the PhILR transform, which incorporates microbial evolutionary models with the isometric log-ratio transform to allow off-the-shelf statistical tools to be safely applied to microbiota surveys. We demonstrate that analyses of community-level structure can be applied to PhILR transformed data with performance on benchmarks rivaling or surpassing standard tools. Additionally, By decomposing distance in the PhILR transformed space, we identified neighboring clades that may have adapted to distinct human body sites. Decomposing variance revealed that covariation of bacterial clades within human body sites increases with phylogenetic relatedness. Together, these findings illustrate how the PhILR transform combines statistical and phylogenetic models to overcome compositional data challenges and enable evolutionary insights relevant to microbial communities.

Data availability

The following previously published data sets were used
    1. Human Microbiome Project Consortium
    (2010) Human Microbiome Project
    Publicly available at HMPDACC (v35 download of files 6, 9, and 10).
    http://hmpdacc.org/HMQCP/
    1. Costello EK
    2. Lauber CL
    3. Hamady M
    4. Fierer N
    5. Gordon JI
    6. Knight R
    (2009) Costello Skin Sites
    Publicly available as part of the FEMS Benchmark dataset (2011) provided Dan Knights.
    http://www.knightslab.org/data

Article and author information

Author details

  1. Justin D Silverman

    Program in Computational Biology and Bioinformatics, Duke University, Durham, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3063-2098

  2. Alex D Washburne

    Nicholas School of the Environment, Duke University, Durham, United States
    Competing interests
    The authors declare that no competing interests exist.

  3. Sayan Mukherjee

    Program in Computational Biology and Bioinformatics, Duke University, Durham, United States
    Competing interests
    The authors declare that no competing interests exist.

  4. Lawrence A David

    Program in Computational Biology and Bioinformatics, Duke University, Durham, United States
    For correspondence
    lawrence.david@duke.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3570-4767

Funding

Global Probiotics Council (Young Investigator Grant for Probiotics Research)

  • Lawrence A David

Searle Scholars Program (15-SSP-184 Research Agreement)

  • Lawrence A David

Alfred P. Sloan Foundation (BR2014-003)

  • Lawrence A David

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Reviewing Editor

  1. Anthony Fodor, University of North Carolina at Charlotte

Version history

  1. Received: September 27, 2016
  2. Accepted: February 13, 2017
  3. Accepted Manuscript published: February 15, 2017 (version 1)
  4. Version of Record published: February 27, 2017 (version 2)

Copyright

© 2017, Silverman et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 11,393
    views
  • 1,717
    downloads
  • 236
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Justin D Silverman
  2. Alex D Washburne
  3. Sayan Mukherjee
  4. Lawrence A David
(2017)
A phylogenetic transform enhances analysis of compositional microbiota data
eLife 6:e21887.
https://doi.org/10.7554/eLife.21887

Share this article

https://doi.org/10.7554/eLife.21887

Categories and tags

Research organism

Further reading

    1. Computational and Systems Biology
    2. Genetics and Genomics

    Imputation of 3D genome structure by genetic–epigenetic interaction modeling in mice

    Lauren Kuffler, Daniel A Skelly ... Gregory W Carter
    Research Article

    Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic–epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic–epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.

    1. Genetics and Genomics
    2. Neuroscience

    Genetic Chimerism: Marmosets contain multitudes

    Kenneth Chiou, Noah Snyder-Mackler
    Insight

    Single-cell RNA sequencing reveals the extent to which marmosets carry genetically distinct cells from their siblings.

    1. Genetics and Genomics

    Elimination of subtelomeric repeat sequences exerts little effect on telomere essential functions in Saccharomyces cerevisiae

    Can Hu, Xue-Ting Zhu ... Jin-Qiu Zhou
    Research Article

    Telomeres, which are chromosomal end structures, play a crucial role in maintaining genome stability and integrity in eukaryotes. In the baker’s yeast Saccharomyces cerevisiae, the X- and Y’-elements are subtelomeric repetitive sequences found in all 32 and 17 telomeres, respectively. While the Y’-elements serve as a backup for telomere functions in cells lacking telomerase, the function of the X-elements remains unclear. This study utilized the S. cerevisiae strain SY12, which has three chromosomes and six telomeres, to investigate the role of X-elements (as well as Y’-elements) in telomere maintenance. Deletion of Y’-elements (SY12YΔ), X-elements (SY12XYΔ+Y), or both X- and Y’-elements (SY12XYΔ) did not impact the length of the terminal TG1-3 tracks or telomere silencing. However, inactivation of telomerase in SY12YΔ, SY12XYΔ+Y, and SY12XYΔ cells resulted in cellular senescence and the generation of survivors. These survivors either maintained their telomeres through homologous recombination-dependent TG1-3 track elongation or underwent microhomology-mediated intra-chromosomal end-to-end joining. Our findings indicate the non-essential role of subtelomeric X- and Y’-elements in telomere regulation in both telomerase-proficient and telomerase-null cells and suggest that these elements may represent remnants of S. cerevisiae genome evolution. Furthermore, strains with fewer or no subtelomeric elements exhibit more concise telomere structures and offer potential models for future studies in telomere biology.