A restriction enzyme reduced representation sequencing approach for low-cost, high-throughput metagenome profiling
Open Access
- 3 April 2020
- journal article
- research article
- Published by Public Library of Science (PLoS) in PLOS ONE
- Vol. 15 (4) , e0219882
Abstract
Microbial community profiles have been associated with a variety of traits, including methane emissions in livestock. These profiles can be difficult and expensive to obtain for thousands of samples (e.g. for accurate association of microbial profiles with traits), therefore 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 ApeKI or PstI, and two bioinformatic pipelines (reference-based and reference-free) were compared to bacterial 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 RE-RRS. The resulting sequences were used to investigate the repeatability of the rumen microbial community profiles, the effect of laboratory and analytical method, and the relationship with methane production. The results suggested that the best method was PstI RE-RRS analyzed with the reference-free approach, which accounted for 53.3±5.9% of reads, and had repeatabilities of 0.49±0.07 and 0.50±0.07 for the first two principal components (PC1 and PC2), phenotypic correlations with methane yield of 0.43±0.06 and 0.46±0.06 for PC1 and PC2, and explained 41±8% of the variation in methane yield. These results were significantly better than for bacterial 16S rRNA gene sequencing of the same samples (p<0.05) except for the correlation between PC2 and methane yield. A Sensitivity study suggested approximately 2000 samples could be sequenced in a single lane on an Illumina HiSeq 2500, meaning the current work using 118 samples/lane and future proposed 384 samples/lane are well within that threshold. With minor adaptations, our approach could be used to obtain microbial profiles from other metagenomic samples.Keywords
All Related Versions
This publication has 36 references indexed in Scilit:
- Metagenomic Predictions: From Microbiome to Complex Health and Environmental Phenotypes in Humans and CattlePLOS ONE, 2013
- Heritability estimates of methane emissions from sheepAnimal, 2013
- The “Most Wanted” Taxa from the Human Microbiome for Whole Genome SequencingPLOS ONE, 2012
- A framework for human microbiome researchNature, 2012
- 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
- QIIME allows analysis of high-throughput community sequencing dataNature Methods, 2010
- BLAST+: architecture and applicationsBMC Bioinformatics, 2009
- MEGAN analysis of metagenomic dataGenome Research, 2007
- Greengenes, a Chimera-Checked 16S rRNA Gene Database and Workbench Compatible with ARBApplied and Environmental Microbiology, 2006