A restriction enzyme reduced representation sequencing approach for low-cost, high-throughput metagenome profiling
Authors:
- Melanie K. HessMHMelanie K. Hess
- Suzanne J. RoweSuzanne J. Rowe
- Tracey C. Van StijnTSTracey C. Van Stijn
- Hannah M. HenryHHHannah M. Henry
- Sharon M. HickeySHSharon M. Hickey
- Rudiger BrauningRBRudiger Brauning
- Alan F. McCullochAMAlan F. McCulloch
- Andrew S. HessAHAndrew S. Hess
- Michelle R. KirkMKMichelle R. Kirk
- Sandra KittelmannSKSandra Kittelmann
- Graham R. WoodGWGraham R. Wood
- Peter H. JanssenPJPeter H. Janssen
Preprint
Abstract
Microbial community profiles have been associated with a variety of traits, including methane emissions in livestock, however, these profiles can be difficult and expensive to obtain for thousands of samples. The objective of this work was to develop a low-cost, high-throughput approach to capture the diversity of the rumen microbiome. Restriction enzyme reduced representation sequencing (RE-RRS) using ApeK I or Pst I, and two bioinformatic pipelines (reference-based and reference-free) were compared to 16S rRNA gene sequencing using repeated samples collected two weeks apart from 118 sheep that were phenotypically extreme (60 high and 58 low) for methane emitted per kg dry matter intake (n=236). DNA was extracted from freeze-dried rumen samples using a phenol chloroform and bead-beating protocol prior to sequencing. The resulting sequences were used to investigate the repeatability of the rumen microbial community profiles, the effect of host genetics, laboratory and analytical method, and the genetic and phenotypic correlations with methane production. The results suggested that the best method was Pst I RE-RRS analyzed with the reference-free approach via a correspondence analysis, with estimates for repeatability of 0.62±0.06, heritability 0.31±0.29, and genetic and phenotypic correlation with methane emissions of 0.88±0.25 and 0.64±0.05 respectively for the first component of correspondence analysis. The reference-free approach assigned 62.0±5.7% of reads to common 65 bp tags, much higher than the reference-based approach of 6.8±1.8% of reads assigned. Sensitivity studies suggested approximately 2000 samples could be sequenced in a single lane on an Illumina HiSeq 2500, therefore the current work of 118 samples/lane and future proposed 384 samples/lane are well within that threshold. Our approach is now being used to investigate host factors affecting the rumen and its association with a variety of production and environmental traits. With minor adaptations, our approach could be used to obtain microbial profiles from other metagenomic samples.Keywords
All Related Versions
- Published version: PLOS ONE, 15 (4), e0219882.
This publication has 26 references indexed in Scilit:
- Improved taxonomic assignment of rumen bacterial 16S rRNA sequences using a revised SILVA taxonomic frameworkPeerJ, 2019
- Host genetics and the rumen microbiome jointly associate with methane emissions in dairy cowsPLOS Genetics, 2018
- Inflammatory Bowel Disease2018
- Coeliac diseaseThe Lancet, 2018
- Specific microbiome-dependent mechanisms underlie the energy harvest efficiency of ruminantsThe ISME Journal, 2016
- Improved OTU-picking using long-read 16S rRNA gene amplicon sequencing and generic hierarchical clusteringMicrobiome, 2015
- Two Different Bacterial Community Types Are Linked with the Low-Methane Emission Trait in SheepPLOS ONE, 2014
- Methane yield phenotypes linked to differential gene expression in the sheep rumen microbiomeGenome Research, 2014
- Improved Detection of Bifidobacteria with Optimised 16S rRNA-Gene Based PyrosequencingPLOS ONE, 2012
- A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity SpeciesPLOS ONE, 2011