Microbial-enrichment method enables high-throughput metagenomic characterization from host-rich samples
Abstract
Host–microbe interactions have been linked to health and disease states through the use of microbial taxonomic profiling, mostly via 16S ribosomal RNA gene sequencing. However, many mechanistic insights remain elusive, in part because studying the genomes of microbes associated with mammalian tissue is difficult due to the high ratio of host to microbial DNA in such samples. Here we describe a microbial-enrichment method (MEM), which we demonstrate on a wide range of sample types, including saliva, stool, intestinal scrapings, and intestinal mucosal biopsies. MEM enabled high-throughput characterization of microbial metagenomes from human intestinal biopsies by reducing host DNA more than 1,000-fold with minimal microbial community changes (roughly 90% of taxa had no significant differences between MEM-treated and untreated control groups). Shotgun sequencing of MEM-treated human intestinal biopsies enabled characterization of both high- and low-abundance microbial taxa, pathways and genes longitudinally along the gastrointestinal tract. We report the construction of metagenome-assembled genomes directly from human intestinal biopsies for bacteria and archaea at relative abundances as low as 1%. Analysis of metagenome-assembled genomes reveals distinct subpopulation structures between the small and large intestine for some taxa. MEM opens a path for the microbiome field to acquire deeper insights into host–microbe interactions by enabling in-depth characterization of host-tissue-associated microbial communities.
Copyright and License
© The Author(s), under exclusive licence to Springer Nature America, Inc. 2023.
Acknowledgement
We acknowledge assistance with animal experiments from Caltech Office of Laboratory Animal Research. We thank M. Ratanapanichkich (California Institute of Technology) for assistance on manual refinement of metagenomic bins and feedback on figure design. We thank A. Carter (California Institute of Technology) for assistance with Quant-seq library preparation, ddPCR measurements and feedback during manuscript preparation. We thank M. Cooper (California Institute of Technology) for identifying appropriate statistical tests, guidance during Quant-seq analysis and feedback on figure design. We thank S. R. Bogatyrev for preliminary investigations, discussions and advice. We thank O. Pradhan (California Institute of Technology) and R. Akana (California Institute of Technology) for advice and feedback during manuscript preparation. We thank B. McDonald (University of Chicago) for providing his expertise and advice on clinical sample collection and processing. We thank A. Wang (University of Chicago) for her assistance in the processing of the human tissue for Figs. 3–6. We thank N. Shelby (California Institute of Technology) for contributions to writing and editing this manuscript. This work was funded in part by a grant from the Kenneth Rainin Foundation (grant no. 2018-1207 to R.F.I.), the Army Research Office Multidisciplinary University Research Initiative (grant no. W911NF-17-1-0402 to R.F.I.), the Jacobs Institute for Molecular Engineering for Medicine, a NIH NIDDK grant (no. RC2 DK133947 to R.F.I. and B.J.), a National Science Foundation Graduate Research Fellowship (grant no. DGE‐1745301 to N.J.W.-W.), and a National Institutes of Health Biotechnology Leadership Pre-doctoral Training Program fellowship from Caltech's Donna and Benjamin M. Rosen Bioengineering Center (grant no. T32GM112592, to J.T.B.), a Helmsley Foundation grant (to F.T.), a NIH NIDDK grant (no. RC2 DK122394, to F.T.), a F30 (grant no. 5F30DK121470, to D.G.S.), a R01 (grant no. DK067180, to B.J.) and the Digestive Diseases Research Core Center grant no. P30 DK42086 at the University of Chicago (to B.J.). The funders had no role in the design of the study, the collection, analysis and interpretation of data, nor in writing the manuscript.
Contributions
These authors contributed equally: Natalie J. Wu-Woods, Jacob T. Barlow.
N.J.W.-W. and J.T.B. conceived and optimized MEM. J.T.B. designed sample collection and analyzed 16S sequencing. D.G.S. codesigned and performed human biopsy collection. N.J.W.-W. and F.T. analyzed shotgun sequencing. A.E.R. performed library preparation. R.F.I. contributed to the design and implementation of the study and to obtaining funding. A.M.E. oversaw the bioinformatic analysis, contributed to the design and implementation of the study and to obtaining funding. B.J. supervised the clinical work, contributed to the design and implementation of the study and to obtaining funding. All authors edited the manuscript. A detailed author contribution statement is available in the Supplementary Information.
Data Availability
The datasets generated and analyzed during the current study are available at CaltechDATA, https://doi.org/10.22002/gx69z-wec80. Microbial sequencing data are available at NCBI Accession no. PRJNA991155. Sequencing data from human samples have been host scrubbed using STAT78 sra-human-scrubber (https://github.com/ncbi/sra-human-scrubber) followed by alignment to CHM13 (ref. 79). Source data are provided with this paper.
Code Availability
The code used in data processing and analysis is available at CaltechDATA, https://doi.org/10.22002/gx69z-wec80.
Conflict of Interest
The work in this paper is the subject of a patent application filed by Caltech (R.F.I., N.J.W.-W., J.T.B. and A.E.R.). The other authors declare no competing interests.
Files
Name | Size | Download all |
---|---|---|
md5:bdae984b2abded7c9aa1b84f21ced6ed
|
36.5 kB | Download |
md5:828cd0cf6104eb3dd2df36b4db7f549f
|
17.4 kB | Download |
md5:86446fa38a615ca764c6be83cc64a1af
|
1.4 MB | Download |
md5:1343634706175c3d884b52a1ba450d83
|
728.8 kB | Download |
md5:9b19936c422a854208eb6810d2ea7173
|
52.2 kB | Preview Download |
md5:8f8d974d819faf359cdba8948fb0c8d6
|
52.8 kB | Preview Download |
md5:69fcd4d040fc71a9ef6f2ded3fb049d9
|
149.2 kB | Download |
md5:22d899d988e82dbfd9f9e88ecd4fc8e8
|
8.2 MB | Download |
md5:82e09ef3b6aca0cb4ce9bc3d3194c402
|
22.4 kB | Download |
md5:ff287ab65495a2e9a686cf19876e8efd
|
60.9 kB | Download |
md5:1a8627ad340ae6460cd5778f0ea219c3
|
16.6 kB | Download |
md5:37bdfdd3c5025f0d5b2085f9aa13833d
|
233.4 kB | Preview Download |
md5:d2ff5b802ad048e83a14e33601ad19eb
|
1.6 MB | Preview Download |
md5:409eed7e1b1009499c5ecc86b75e1b37
|
25.8 kB | Download |
md5:ad12786ad0abc5ddb6316e38c464f991
|
251.5 kB | Download |
md5:3d7941add7a46dda03155dcd8082a9a5
|
43.3 kB | Download |
md5:812ad840369fc25bdeb06ab8a06b2b7c
|
95.0 kB | Preview Download |
md5:3aa0e46bec0b7375d1cede20fbe9bf8c
|
180.3 kB | Preview Download |
Additional details
- ISSN
- 1548-7105
- URL
- https://rdcu.be/dowgH
- Kenneth Rainin Foundation
- 2018-1207
- United States Army Research Office
- W911NF-17-1-0402
- National Institutes of Health
- RC2 DK133947
- National Science Foundation
- DGE‐1745301
- National Institutes of Health
- T32GM112592
- Leona M. and Harry B. Helmsley Charitable Trust
- National Institutes of Health
- RC2 DK122394
- National Institutes of Health
- 5F30DK121470
- National Institutes of Health
- DK067180
- National Institutes of Health
- P30 DK42086
- Caltech groups
- Jacobs Institute for Molecular Engineering for Medicine, Rosen Bioengineering Center, Division of Biology and Biological Engineering