High-throughput single-microbe RNA sequencing reveals adaptive state heterogeneity and host-phage activity associations in human gut microbiome

Abstract

Microbial communities such as those residing in the human gut are highly diverse and complex, and many with important implications for health and diseases. The effects and functions of these microbial communities are determined not only by their species compositions and diversities but also by the dynamic intra- and inter-cellular states at the transcriptional level. Powerful and scalable technologies capable of acquiring single-microbe-resolution RNA sequencing information in order to achieve a comprehensive understanding of complex microbial communities together with their hosts are therefore utterly needed. Here we report the development and utilization of a droplet-based smRNA-seq (single-microbe RNA sequencing) method capable of identifying large species varieties in human samples, which we name smRandom-seq2. Together with a triple-module computational pipeline designed for the bacteria and bacteriophage sequencing data by smRandom-seq2 in four human gut samples, we established a single-cell level bacterial transcriptional landscape of human gut microbiome, which included 29,742 single microbes and 329 unique species. Distinct adaptive response states among species in Prevotella and Roseburia genera and intrinsic adaptive strategy heterogeneity in Phascolarctobacterium succinatutens were uncovered. Additionally, we identified hundreds of novel host-phage transcriptional activity associations in the human gut microbiome. Our results indicated that smRandom-seq2 is a high-throughput and high-resolution smRNA-seq technique that is highly adaptable to complex microbial communities in real-world situations and promises new perspectives in the understanding of human microbiomes.

Keywords

single-microbe RNA sequencing (smRNA-seq) / droplet microfluidics / microbiome / host-phage association / smRandom-seq2

References

[1] Almeida A , Nayfach S , Boland M . A unified catalog of 204,938 reference genomes from the human gut microbiome. Nat Biotechnol 2021; 39: 105- 114.

[2] Blattman SB , Jiang W , Oikonomou P . Prokaryotic single-cell RNA sequencing by in situ combinatorial indexing. Nat Microbiol 2020; 5: 1192- 1201.

[3] Camarillo-Guerrero LF , Almeida A , Rangel-Pineros G . Massive expansion of human gut bacteriophage diversity. Cell 2021; 184: 1098- 1109.e9.

[4] Chong S , Chen C , Ge H . Mechanism of transcriptional bursting in bacteria. Cell 2014; 158: 314- 326.

[5] Dar D , Dar N , Cai L . Spatial transcriptomics of planktonic and sessile bacterial populations at single-cell resolution. Science 2021; 373: eabi4882.

[6] Dobin A , Davis CA , Schlesinger F . STAR: ultrafast universal RNA-seq aligner. Bioinformatics 2013; 29: 15- 21.

[7] Fujimoto K , Kimura Y , Allegretti JR . Functional restoration of bacteriomes and viromes by fecal microbiota transplantation. Gastroenterology 2021; 160: 2089- 2102.e12.

[8] Han J , Depinho RA , Maitra A . Single-cell RNA sequencing in pancreatic cancer. Nat Rev Gastroenterol Hepatol 2021; 18: 451- 452.

[9] Hao Y , Hao S , Andersen-Nissen E . Integrated analysis of multimodal single-cell data. Cell 2021; 184: 3573- 3587.e29.

[10] Ikeyama N , Murakami T , Toyoda A . Microbial interaction between the succinate-utilizing bacterium Phascolarctobacterium faecium and the gut commensal Bacteroides thetaiotaomicron. Microbiologyopen 2020; 9: e1111.

[11] Imdahl F , Vafadarnejad E , Homberger C . Single-cell RNAsequencing reports growth-condition-specific global transcriptomes of individual bacteria. Nat Microbiol 2020; 5: 1202- 1206.

[12] Klein AM , Mazutis L , Akartuna I . Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells. Cell 2015; 161: 1187- 1201.

[13] Ko J , Wang Y , Sheng K . Sequencing-based protein analysis of single extracellular vesicles. ACS Nano 2021; 15: 5631- 5638.

[14] Kopylova E , Noé L , Touzet H . SortMeRNA: fast and accurate filtering of ribosomal RNAs in metatranscriptomic data. Bioinformatics 2012; 28: 3211- 3217.

[15] Kuchina A , Brettner LM , Paleologu L . Microbial single-cell RNA sequencing by split-pool barcoding. Science 2021; 371: eaba5257.

[16] Liao Y , Smyth GK , Shi W . featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 2014; 30: 923- 930.

[17] Lu J , Rincon N , Wood DE . Metagenome analysis using the Kraken software suite. Nat Protoc 2022; 17: 2815- 2839.

[18] Ma P , Amemiya HM , He LL . Bacterial droplet-based single-cell RNA-seq reveals antibiotic-associated heterogeneous cellular states. Cell 2023; 186: 877- 891.e14.

[19] Manrique P , Bolduc B , Walk S T . Healthy human gut phageome. Proc Natl Acad Sci U S A 2016; 113: 10400- 10405.

[20] Mayneris-Perxachs J , Castells-Nobau A , Arnoriaga-Rodríguez M . Microbiota alterations in proline metabolism impact depression. Cell Metab 2022; 34: 681- 701.e10.

[21] Mcnulty R , Sritharan D , Pahng SH . Probe-based bacterial single-cell RNA sequencing predicts toxin regulation. Nat Microbiol 2023; 8: 934- 945.

[22] Natarajan A , Bhatt AS . Microbes and microbiomes in 2020 and beyond. Nat Commun 2020; 11: 4988.

[23] Paik DT , Cho S , Tian L . Single-cell RNA sequencing in cardiovascular development, disease and medicine. Nat Rev Cardiol 2020; 17: 457- 473.

[24] Papalexi E , Satija R . Single-cell RNA sequencing to explore immune cell heterogeneity. Nat Rev Immunol 2018; 18: 35- 45.

[25] Patterson AM , Mulder IE , Travis AJ . Human gut symbiont roseburia hominis promotes and regulates innate immunity. Front Immunol 2017; 8: 1166.

[26] Prasoodanan PKV , Sharma AK , Mahajan S . Western and non-western gut microbiomes reveal new roles of Prevotella in carbohydrate metabolism and mouth-gut axis. npj Biofilms Microbiomes 2021; 7: 77.

[27] Roemhild R , Bollenbach T , Andersson DI . The physiology and genetics of bacterial responses to antibiotic combinations. Nat Rev Microbiol 2022; 20: 478- 490.

[28] Sharma A , Gilbert JA . Microbial exposure and human health. Curr Opin Microbiol 2018; 44: 79- 87.

[29] Sheng K , Cao W , Niu Y . Effective detection of variation in single-cell transcriptomes using MATQ-seq. Nat Methods 2017; 14: 267- 270.

[30] Shkoporov AN , Turkington CJ , Hill C . Mutualistic interplay between bacteriophages and bacteria in the human gut. Nat Rev Microbiol 2022; 20: 737- 749.

[31] Smith T , Heger A , Sudbery I . UMI-tools: modeling sequencing errors in Unique Molecular Identifiers to improve quantification accuracy. Genome Res 2017; 27: 491- 499.

[32] Smith, WPJ , Wucher BR , Nadell CD . Bacterial defences:mechanisms, evolution and antimicrobial resistance. Nat Rev Microbiol 2023; 21: 519- 534.

[33] Van De Sande B , Lee JS , Mutasa-Gottgens E . applications of single-cell RNA sequencing in drug discovery and development. Nat Rev Drug Discov 2023; 22: 496- 520.

[34] Van Der Leun AM , Thommen DS , Schumacher TN . CD8⁺ T cell states in human cancer: insights from single-cell analysis. Nat Rev Cancer 2020; 20: 218- 232.

[35] Wang Y , Cao T , Ko J . Dissolvable polyacrylamide beads for high-throughput droplet DNA barcoding. Adv Sci(Weinh) 2020; 7: 1903463.

[36] Xu Z , Wang Y , Sheng K . Droplet-based high-throughput single microbe RNA sequencing by smRandom-seq. Nat Commun 2023a; 14: 5130.

[37] Xu Z , Zhang T , Chen H . High-throughput single nucleus total RNA sequencing of formalin-fixed paraffin-embedded tissues by snRandom-seq. Nat Commun 2023b; 14: 2734.

RIGHTS & PERMISSIONS

The Author(s). Published by Oxford University Press on behalf of Higher Education Press.