Uncovering Endolysins against Methicillin-Resistant Staphylococcus aureus Using a Microbial Single-Cell Genome Database.

Endolysins, peptidoglycan hydrolases derived from bacteriophages (phages), are being developed as a promising alternative to conventional antibiotics. To obtain highly active endolysins, a diverse library of these endolysins is vital. We propose here microbial single-cell genome sequencing as an efficient tool to discover dozens of previously unknown endolysins, owing to its culture-independent sequencing method. As a proof of concept, we analyzed and recovered endolysin genes within prophage regions of Staphylococcus single-amplified genomes in human skin microbiome samples. We constructed a library of chimeric endolysins by shuffling domains of the natural endolysins and performed high-throughput screening against Staphylococcus aureus. One of the lead endolysins, bbst1027, exhibited desirable antimicrobial properties, such as rapid bactericidal activity, no detectable resistance development, and in vivo efficacy. We foresee that this endolysin discovery pipeline is in principle applicable to any bacterial target and boost the development of novel antimicrobial agents.

Analysis of microbial dynamics in the soybean root-associated environments from community to single-cell levels.

Plant root-associated environments such as the rhizosphere, rhizoplane, and endosphere, are notably different from non-root-associated soil environments. However, the microbial dynamics in these spatially divided compartments remain unexplored. In this study, we propose a combinational analysis of single-cell genomics with 16S rRNA gene sequencing. This method enabled us to understand the entire soil microbiome and individual root-associated microorganisms. We applied this method to soybean microbiomes and revealed that their composition was different between the rhizoplane and rhizosphere in the early growth stages, but became more similar as growth progressed. In addition, a total of 610 medium- to high-quality single-amplified genomes (SAGs) were acquired, including plant growth-promoting rhizobacteria (PGPR) candidates while genomes with high GC content tended to be missed by SAGs. The whole-genome analyses of the SAGs suggested that rhizoplane-enriched Flavobacterium solubilizes organophosphate actively and Bacillus colonizes roots more efficiently. Single-cell genomics, together with 16S rRNA gene sequencing, enabled us to connect microbial taxonomy and function, and assess microorganisms at a strain resolution even in the complex soil microbiome.

Tools for microbial single-cell genomics for obtaining uncultured microbial genomes.

The advent of next-generation sequencing technologies has facilitated the acquisition of large amounts of DNA sequence data at a relatively low cost, leading to numerous breakthroughs in decoding microbial genomes. Among the various genome sequencing activities, metagenomic analysis, which entails the direct analysis of uncultured microbial DNA, has had a profound impact on microbiome research and has emerged as an indispensable technology in this field. Despite its valuable contributions, metagenomic analysis is a "bulk analysis" technique that analyzes samples containing a wide diversity of microbes, such as bacteria, yielding information that is averaged across the entire microbial population. In order to gain a deeper understanding of the heterogeneous nature of the microbial world, there is a growing need for single-cell analysis, similar to its use in human cell biology. With this paradigm shift in mind, comprehensive single-cell genomics technology has become a much-anticipated innovation that is now poised to revolutionize microbiome research. It has the potential to enable the discovery of differences at the strain level and to facilitate a more comprehensive examination of microbial ecosystems. In this review, we summarize the current state-of-the-art in microbial single-cell genomics, highlighting the potential impact of this technology on our understanding of the microbial world. The successful implementation of this technology is expected to have a profound impact in the field, leading to new discoveries and insights into the diversity and evolution of microbes.

Uncultured prokaryotic genomes in the spotlight: An examination of publicly available data from metagenomics and single-cell genomics.

Owing to the ineffectiveness of traditional culture techniques for the vast majority of microbial species, culture-independent analyses utilizing next-generation sequencing and bioinformatics have become essential for gaining insight into microbial ecology and function. This mini-review focuses on two essential methods for obtaining genetic information from uncultured prokaryotes, metagenomics and single-cell genomics. We analyzed the registration status of uncultured prokaryotic genome data from major public databases and assessed the advantages and limitations of both the methods. Metagenomics generates a significant quantity of sequence data and multiple prokaryotic genomes using straightforward experimental procedures. However, in ecosystems with high microbial diversity, such as soil, most genes are presented as brief, disconnected contigs, and lack association of highly conserved genes and mobile genetic elements with individual species genomes. Although technically more challenging, single-cell genomics offers valuable insights into complex ecosystems by providing strain-resolved genomes, addressing issues in metagenomics. Recent technological advancements, such as long-read sequencing, machine learning algorithms, and in silico protein structure prediction, in combination with vast genomic data, have the potential to overcome the current technical challenges and facilitate a deeper understanding of uncultured microbial ecosystems and microbial dark matter genes and proteins. In light of this, it is imperative that continued innovation in both methods and technologies take place to create high-quality reference genome databases that will support future microbial research and industrial applications.

Enhancing the sensitivity of bacterial single-cell RNA sequencing using RamDA-seq and Cas9-based rRNA depletion.

Bacterial populations exhibit heterogeneity in gene expression, which facilitates their survival and adaptation to unstable and unpredictable environments through the bet-hedging strategy. However, unraveling the rare subpopulations and heterogeneity in gene expression using population-level gene expression analysis remains a challenging task. Single-cell RNA sequencing (scRNA-seq) has the potential to identify rare subpopulations and capture heterogeneity in bacterial populations, but standard methods for scRNA-seq in bacteria are still under development, mainly due to differences in mRNA abundance and structure between eukaryotic and prokaryotic organisms. In this study, we present a hybrid approach that combines random displacement amplification sequencing (RamDA-seq) with Cas9-based rRNA depletion for scRNA-seq in bacteria. This approach allows cDNA amplification and subsequent sequencing library preparation from low-abundance bacterial RNAs. We evaluated its sequenced read proportion, gene detection sensitivity, and gene expression patterns from the dilution series of total RNA or the sorted single Escherichia coli cells. Our results demonstrated the detection of more than 1000 genes, about 24% of the genes in the E. coli genome, from single cells with less sequencing effort compared to conventional methods. We observed gene expression clusters between different cellular proliferation states or heat shock treatment. The approach demonstrated high detection sensitivity in gene expression analysis compared to current bacterial scRNA-seq methods and proved to be an invaluable tool for understanding the ecology of bacterial populations and capturing the heterogeneity of bacterial gene expression.

Target enrichment of uncultured human oral bacteria with phage-derived molecules found by single-cell genomics.

Advances in culture-independent microbial analysis, such as metagenomics and single-cell genomics, have significantly increased our understanding of microbial lineages. While these methods have uncovered a large number of novel microbial taxa, many remain uncultured, and their function and mode of existence in the environment are still unknown. This study aims to explore the use of bacteriophage-derived molecules as probes for detecting and isolating uncultured bacteria. Here, we proposed multiplex single-cell sequencing to obtain massive uncultured oral bacterial genomes and searched prophage sequences from over 450 obtained human oral bacterial single-amplified genomes (SAGs). The focus was on the cell wall binding domain (CBD) in phage endolysin, and fluorescent protein-fused CBDs were generated based on several CBD gene sequences predicted from Streptococcus SAGs. The ability of the Streptococcus prophage-derived CBDs to detect and enrich specific Streptococcus species from human saliva while maintaining cell viability was confirmed by magnetic separation and flow cytometry. The approach to phage-derived molecule generation based on uncultured bacterial SAG is expected to improve the process of designing molecules that selectively capture or detect specific bacteria, notably from uncultured gram-positive bacteria, and will have applications in isolation and in situ detection of beneficial or pathogenic bacteria.

Combined actions of bacteriophage-encoded genes in Wolbachia-induced male lethality.

Some Wolbachia endosymbionts induce male killing, whereby male offspring of infected females are killed during development; however, the origin and diversity of the underlying mechanisms remain unclear. In this study, we identified a 76 kbp prophage region specific to male-killing Wolbachia hosted by the moth Homona magnanima. The prophage encoded a homolog of the male-killing gene oscar in Ostrinia moths and the wmk gene that induces various toxicities in Drosophila melanogaster. Upon overexpressing these genes in D. melanogaster, wmk-1 and wmk-3 killed all males and most females, whereas Hm-oscar, wmk-2, and wmk-4 had no impact on insect survival. Strikingly, co-expression of tandemly arrayed wmk-3 and wmk-4 killed 90% of males and restored 70% of females, suggesting their conjugated functions for male-specific lethality. While the male-killing gene in the native host remains unknown, our findings highlight the role of bacteriophages in male-killing evolution and differences in male-killing mechanisms among insects.

Revealing within-species diversity in uncultured human gut bacteria with single-cell long-read sequencing.

Obtaining complete and accurate bacterial genomes is vital for studying the characteristics of uncultured bacteria. Single-cell genomics is a promising approach for the culture-independent recovery of bacterial genomes from individual cells. However, single-amplified genomes (SAGs) often have fragmented and incomplete sequences due to chimeric and biased sequences introduced during the genome amplification process. To address this, we developed a single-cell amplified genome long-read assembly (scALA) workflow to construct complete circular SAGs (cSAGs) from long-read single-cell sequencing data of uncultured bacteria. We used the SAG-gel platform, which is both cost-effective and high-throughput, to obtain hundreds of short-read and long-read sequencing data for specific bacterial strains. The scALA workflow generated cSAGs by repeated in silico processing for sequence bias reduction and contig assembly. From 12 human fecal samples, including two cohabitant groups, scALA generated 16 cSAGs of three specifically targeted bacterial species: Anaerostipes hadrus, Agathobacter rectalis, and Ruminococcus gnavus. We discovered strain-specific structural variations shared among cohabiting hosts, while all cSAGs of the same species showed high homology in aligned genomic regions. A. hadrus cSAGs exhibited 10 kbp-long phage insertions, various saccharide metabolic capabilities, and different CRISPR-Cas systems in each strain. The sequence similarity of A. hadrus genomes did not necessarily correspond with orthologous functional genes, while host geographical regionality seemed to be highly related to gene possession. scALA allowed us to obtain closed circular genomes of specifically targeted bacteria from human microbiota samples, leading to an understanding of within-species diversities, including structural variations and linking mobile genetic elements, such as phages, to hosts. These analyses provide insight into microbial evolution, the adaptation of the community to environmental changes, and interactions with hosts. cSAGs constructed using this method can expand bacterial genome databases and our understanding of within-species diversities in uncultured bacteria.

Linking antigen specific T-cell dynamics in a microfluidic chip to single cell transcription patterns.

A new non-invasive screening profile has been realized that can aid in determining T-cell activation state at single-cell level. Production of activated T-cells with good specificity and stable proliferation is greatly beneficial for advancing adoptive immunotherapy as innate immunological cells are not effective in recognizing and eliminating cancer as expected. The screening method is realized by relating intracellular Ca2+ intensity and motility of T-cells interacting with APC (Antigen Presenting Cells) in a microfluidic chip. The system is tested using APC pulsed with OVA257-264 peptide and its modified affinities (N4, Q4, T4 and V4), and the T-cells from OT-1 mice. In addition, single cell RNA sequencing reveals the activation states of the cells and the clusters from the derived profiles can be indicative of the T-cell activation state. The presented system here can be versatile for a comprehensive application to proceed with T-cell-based immunotherapy and screen the antigen-specific T-cells with excellent efficiency and high proliferation.

Targeted single-cell genomics reveals novel host adaptation strategies of the symbiotic bacteria Endozoicomonas in Acropora tenuis coral.

BACKGROUND: Endozoicomonas bacteria symbiosis with various marine organisms is hypothesized as a potential indicator of health in corals. Although many amplicon analyses using 16S rRNA gene have suggested the diversity of Endozoicomonas species, genome analysis has been limited due to contamination of host-derived sequences and difficulties in culture and metagenomic analysis. Therefore, the evolutionary and functional potential of individual Endozoicomonas species symbiotic with the same coral species remains unresolved. RESULTS: In this study, we applied a novel single-cell genomics technique using droplet microfluidics to obtain single-cell amplified genomes (SAGs) for uncultured coral-associated Endozoicomonas spp. We obtained seven novel Endozoicomonas genomes and quantitative bacterial composition from Acropora tenuis corals at four sites in Japan. Our quantitative 16S rRNA gene and comparative genomic analysis revealed that these Endozoicomonas spp. belong to different lineages (Clade A and Clade B), with widely varying abundance among individual corals. Furthermore, each Endozoicomonas species possessed various eukaryotic-like genes in clade-specific genes. It was suggested that these eukaryotic-like genes might have a potential ability of different functions in each clade, such as infection of the host coral or suppression of host immune pathways. These Endozoicomonas species may have adopted different host adaptation strategies despite living symbiotically on the same coral. CONCLUSIONS: This study suggests that coral-associated Endozoicomonas spp. on the same species of coral have different evolutional strategies and functional potentials in each species and emphasizes the need to analyze the genome of each uncultured strain in future coral-Endozoicomonas relationships studies. Video Abstract.

Integrated spatial analysis of gene mutation and gene expression for understanding tumor diversity in formalin-fixed paraffin-embedded lung adenocarcinoma.

Introduction: A deeper understanding of intratumoral heterogeneity is essential for prognosis prediction or accurate treatment plan decisions in clinical practice. However, due to the cross-links and degradation of biomolecules within formalin-fixed paraffin-embedded (FFPE) specimens, it is challenging to analyze them. In this study, we aimed to optimize the simultaneous extraction of mRNA and DNA from microdissected FFPE tissues (φ = 100 µm) and apply the method to analyze tumor diversity in lung adenocarcinoma before and after erlotinib administration. Method: Two magnetic beads were used for the simultaneous extraction of mRNA and DNA. The decross-linking conditions were evaluated for gene mutation and gene expression analyses of microdissected FFPE tissues. Lung lymph nodes before treatment and lung adenocarcinoma after erlotinib administration were collected from the same patient and were preserved as FFPE specimens for 4 years. Gene expression and gene mutations between histologically classified regions of lung adenocarcinoma (pre-treatment tumor in lung lymph node biopsies and post-treatment tumor, normal lung, tumor stroma, and remission stroma, in resected lung tissue) were compared in a microdissection-based approach. Results: Using the optimized simultaneous extraction of DNA and mRNA and whole-genome amplification, we detected approximately 4,000-10,000 expressed genes and the epidermal growth factor receptor (EGFR) driver gene mutations from microdissected FFPE tissues. We found the differences in the highly expressed cancer-associated genes and the positive rate of EGFR exon 19 deletions among the tumor before and after treatment and tumor stroma, even though they were collected from tumors of the same patient or close regions of the same specimen. Conclusion: Our integrated spatial analysis method would be applied to various FFPE pathology specimens providing area-specific gene expression and gene mutation information.

Reproducible and sensitive micro-tissue RNA sequencing from formalin-fixed paraffin-embedded tissues for spatial gene expression analysis.

Spatial transcriptome analysis of formalin-fixed paraffin-embedded (FFPE) tissues using RNA-sequencing (RNA-seq) provides interactive information on morphology and gene expression, which is useful for clinical applications. However, despite the advantages of long-term storage at room temperature, FFPE tissues may be severely damaged by methylene crosslinking and provide less gene information than fresh-frozen tissues. In this study, we proposed a sensitive FFPE micro-tissue RNA-seq method that combines the punching of tissue sections (diameter: 100 µm) and the direct construction of RNA-seq libraries. We evaluated a method using mouse liver tissues at two years after fixation and embedding and detected approximately 7000 genes in micro-punched tissue-spots (thickness: 10 µm), similar to that detected with purified total RNA (2.5 ng) equivalent to the several dozen cells in the spot. We applied this method to clinical FFPE specimens of lung cancer that had been fixed and embedded 6 years prior, and found that it was possible to determine characteristic gene expression in the microenvironment containing tumor and non-tumor cells of different morphologies. This result indicates that spatial gene expression analysis of the tumor microenvironment is feasible using FFPE tissue sections stored for extensive periods in medical facilities.

Massively parallel single-cell genomics of microbiomes in rice paddies.

Plant growth-promoting microbes (PGPMs) have attracted increasing attention because they may be useful in increasing crop yield in a low-input and sustainable manner to ensure food security. Previous studies have attempted to understand the principles underlying the rhizosphere ecology and interactions between plants and PGPMs using ribosomal RNA sequencing, metagenomic sequencing, and genome-resolved metagenomics; however, these approaches do not provide comprehensive genomic information for individual species and do not facilitate detailed analyses of plant-microbe interactions. In the present study, we developed a pipeline to analyze the genomic diversity of the rice rhizosphere microbiome at single-cell resolution. We isolated microbial cells from paddy soil and determined their genomic sequences by using massively parallel whole-genome amplification in microfluidic-generated gel capsules. We successfully obtained 3,237 single-amplified genomes in a single experiment, and these genomic sequences provided insights into microbial functions in the paddy ecosystem. Our approach offers a promising platform for gaining novel insights into the roles of microbes in the rice rhizomicrobiome and to develop microbial technologies for improved and sustainable rice production.

Cancer Cachexia among Patients with Advanced Non-Small-Cell Lung Cancer on Immunotherapy: An Observational Study with Exploratory Gut Microbiota Analysis.

Cancer cachexia exerts a negative clinical influence on patients with advanced non-small-cell lung cancer (NSCLC) treated with immune checkpoint inhibitors (ICI). The prognostic impact of body weight change during ICI treatment remains unknown. The gut microbiota (GM) is a key contributor to the response to ICI therapy in cancer patients. However, the association between cancer cachexia and GM and their association with the response to ICIs remains unexplored. This study examined the association of cancer cachexia with GM composition and assessed the impact of GM on clinical outcomes in patients with NSCLC treated with ICIs. In this observational, prospective study, which included 113 Japanese patients with advanced NSCLC treated with ICIs, the prevalence of cachexia was 50.4% (57/113). The median progression-free survival (PFS) and overall survival (OS) were significantly shorter in the cachexia group than in the non-cachexia group (4.3 vs. 11.6 months (p = 0.003) and 12.0 months vs. not reached (p = 0.02), respectively). A multivariable analysis revealed that baseline cachexia was independently associated with a shorter PFS. Moreover, a gain in body weight from the baseline (reversible cachexia) was associated with a significantly longer PFS and OS compared to irreversible cachexia. Microbiome profiling with 16S rRNA analysis revealed that the cachexia group presented an overrepresentation of the commensal bacteria, Escherichia-Shigella and Hungatella, while the non-cachexia group had a preponderance of Anaerostipes, Blautia, and Eubacterium ventriosum. Anaerostipes and E. ventriosum were associated with longer PFS and OS. Moreover, a cachexia status correlated with the systemic inflammatory marker-derived-neutrophil-to-lymphocytes ratio (dNLR) and Lung Immune Prognostic Index (LIPI) indexes. Our study demonstrates that cachexia and longitudinal bodyweight change have a prognostic impact on patients with advanced NSCLC treated with ICI therapy. Moreover, our study demonstrates that bacteria associated with ICI resistance are also linked to cachexia. Targeted microbiota interventions may represent a new type of treatment to overcome cachexia in patients with NSCLC.

Exploring strain diversity of dominant human skin bacterial species using single-cell genome sequencing.

To understand the role of the skin commensal bacterial community in skin health and the spread of pathogens, it is crucial to identify genetic differences in the bacterial strains corresponding to human individuals. A culture-independent genomics approach is an effective tool for obtaining massive high-quality bacterial genomes. Here we present a single-cell genome sequencing to obtain comprehensive whole-genome sequences of uncultured skin bacteria from skin swabs. We recovered 281 high-quality (HQ) and 244 medium-quality single-amplified genomes (SAGs) of multiple skin bacterial species from eight individuals, including cohabiting group. Single-cell sequencing outperformed in the genome recovery from the same skin swabs, showing 10-fold non-redundant strain genomes compared to the shotgun metagenomic sequencing and binning approach. We then focused on the abundant skin bacteria and identified intra-species diversity, especially in 47 Moraxella osloensis derived HQ SAGs, characterizing the strain-level heterogeneity at mobile genetic element profiles, including plasmids and prophages. Even between the cohabiting individual hosts, they have unique skin bacterial strains in the same species, which shows microdiversity in each host. Genetic and functional differences between skin bacterial strains are predictive of in vivo competition to adapt bacterial genome to utilize the sparse nutrients available on the skin or produce molecules that inhibit the colonization of other microbes or alter their behavior. Thus, single-cell sequencing provides a large number of genomes of higher resolution and quality than conventional metagenomic analysis and helps explore the skin commensal bacteria at the strain level, linking taxonomic and functional information.

Identification of two cancer stem cell-like populations in triple-negative breast cancer xenografts.

Gene expression analysis at the single-cell level by next-generation sequencing has revealed the existence of clonal dissemination and microheterogeneity in cancer metastasis. The current spatial analysis technologies can elucidate the heterogeneity of cell-cell interactions in situ. To reveal the regional and expressional heterogeneity in primary tumors and metastases, we performed transcriptomic analysis of microtissues dissected from a triple-negative breast cancer (TNBC) cell line MDA-MB-231 xenograft model with our automated tissue microdissection punching technology. This multiple-microtissue transcriptome analysis revealed three cancer cell-type clusters in the primary tumor and axillary lymph node metastasis, two of which were cancer stem cell (CSC)-like clusters (CD44/MYC-high, HMGA1-high). Reanalysis of public single-cell RNA-sequencing datasets confirmed that the two CSC-like populations existed in TNBC xenograft models and in TNBC patients. The diversity of these multiple CSC-like populations could cause differential anticancer drug resistance, increasing the difficulty of curing this cancer.

Identification of lipolytic enzymes using high-throughput single-cell screening and sorting of a metagenomic library.

The demand for novel, robust microbial biocatalysts for use in industrial and pharmaceutical applications continues to increase rapidly. As a result, there is a need to develop advanced tools and technologies to exploit the vast metabolic potential of unculturable microorganisms found in various environments. Single-cell and functional metagenomics studies can explore the enzymatic potential of entire microbial communities in a given environment without the need to culture the microorganisms. This approach has contributed substantially to the discovery of unique microbial genes for industrial and medical applications. Functional metagenomics involves the extraction of microbial DNA directly from environmental samples, constructing expression libraries comprising the entire microbial genome, and screening of the libraries for the presence of desired phenotypes. In this study, lipolytic enzymes from the Red Sea were targeted. A high-throughput single-cell microfluidic platform combined with a laser-based fluorescent screening bioassay was employed to discover new genes encoding lipolytic enzymes. Analysis of the metagenomic library led to the identification of three microbial genes encoding lipases based on their functional similarity and sequence homology to known lipases. The results demonstrated that microfluidics is a robust technology that can be used for screening in functional metagenomics. The results also indicate that the Red Sea is a promising, under-investigated source of new genes and gene products.

Strain-level profiling of viable microbial community by selective single-cell genome sequencing.

Culture-independent analysis with high-throughput sequencing has been widely used to characterize bacterial communities. However, signals derived from non-viable bacteria and non-cell DNA may inhibit its characterization. Here, we present a method for viable bacteria-targeted single-cell genome sequencing, called PMA-SAG-gel, to obtain comprehensive whole-genome sequences of surviving uncultured bacteria from microbial communities. PMA-SAG-gel uses gel matrixes that enable sequential enzymatic reactions for cell lysis and genome amplification of viable single cells from the microbial communities. PMA-SAG-gel removed the single-amplified genomes (SAGs) derived from dead bacteria and enabled selective sequencing of viable bacteria in the model samples of Escherichia coli and Bacillus subtilis. Next, we demonstrated the recovery of near-complete SAGs of eight oxygen-tolerant bacteria, including Bacteroides spp. and Phocaeicola spp., from 1331 human feces SAGs. We found the presence of two different strains in each species and identified their specific genes to investigate the metabolic functions. The survival profile of an entire population at the strain level will provide the information for understanding the characteristics of the surviving bacteria under the specific environments or sample processing and insights for quality assessment of live bacterial products or fecal microbiota transplantation and for understanding the effect of antimicrobial treatments.

Validation of the application of gel beads-based single-cell genome sequencing platform to soil and seawater.

Single-cell genomics is applied to environmental samples as a method to solve the problems of current metagenomics. However, in the fluorescence-activated cell sorting-based cell isolation and subsequent whole genome amplification, the sorting efficiency and the sequence quality are greatly affected by the type of target environment, limiting its adaptability. Here, we developed an improved single-cell genomics platform, named SAG-gel, which utilizes gel beads for single-cell isolation, lysis, and whole genome amplification. To validate the versatility of SAG-gel, single-cell genome sequencing was performed with model bacteria and microbial samples collected from eight environmental sites, including soil and seawater. Gel beads enabled multiple lysis treatments. The genome coverage with model bacteria was improved by 9.1-25%. A total of 734 single amplified genomes were collected from the diverse environmental samples, and almost full-length 16S rRNA genes were recovered from 57.8% of them. We also revealed two marine Rhodobacter strains harboring nearly identical 16S rRNA genes but having different genome contents. In addition, searching for viral sequences elucidated the virus-host linkage over the sampling sites, revealing the geographic distribution and diverse host range of viruses.

Single-cell metabolite detection and genomics reveals uncultivated talented producer.

The production of bioactive metabolites is increasingly recognized as an important function of host-associated bacteria. An example is defensive symbiosis that might account for much of the chemical richness of marine invertebrates including sponges (Porifera), 1 of the oldest metazoans. However, most bacterial members of sponge microbiomes have not been cultivated or sequenced, and therefore, remain unrecognized. Unequivocally linking metabolic functions to a cellular source in sponge microbiomes is, therefore, a challenge. Here, we report an analysis pipeline of microfluidic encapsulation, Raman microscopy, and integrated digital genomics (MERMAID) for an efficient identification of uncultivated producers. We applied this method to the chemically rich bacteriosponge (sponge that hosts a rich bacterial community) Theonella swinhoei, previously shown to contain 'Entotheonella' symbionts that produce most of the bioactive substances isolated from the sponge. As an exception, the antifungal aurantosides had remained unassigned to a source. Raman-guided single-bacterial analysis and sequencing revealed a cryptic, distinct multiproducer, 'Candidatus Poriflexus aureus' from a new Chloroflexi lineage as the aurantoside producer. Its exceptionally large genome contains numerous biosynthetic loci and suggested an even higher chemical richness of this sponge than previously appreciated. This study highlights the importance of complementary technologies to uncover microbiome functions, reveals remarkable parallels between distantly related symbionts of the same host, and adds functional support for diverse chemically prolific lineages being present in microbial dark matter.