Search

Abstract Background The increase in metagenome-assembled genomes (MAGs) has advanced our understanding of the functional characterization and taxonomic assignment within the human microbiome. However, MAGs, as population consensus genomes, often aggregate heterogeneity among species and strains, thereby obfuscating the precise relationships between microbial hosts and mobile genetic elements (MGEs). In contrast, single amplified genomes (SAGs) derived via single-cell genome sequencing can capture individual genomic content, including MGEs. Results We introduce the first substantial SAG dataset (bbsag20) from the human oral and gut microbiome, comprising 17,202 SAGs above medium-quality without co-assembly. This collection unveils a diversity of bacterial lineages across 312 oral and 647 gut species, demonstrating different taxonomic compositions from MAGs. Moreover, the SAGs showed cellular-level evidence of the translocation of oral bacteria to the gut. We also identified broad-host-range MGEs harboring antibiotic resistance genes (ARGs), which were not detected in the MAGs. Conclusions The difference in taxonomic composition between SAGs and MAGs indicates that combining both methods would be effective in expanding the genome catalog. By connecting mobilomes and resistomes in individual samples, SAGs could meticulously chart a dynamic network of ARGs on MGEs, pinpointing potential ARG reservoirs and their spreading patterns in the microbial community. Video Abstract

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.

Abstract 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

Abstract 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.

Abstract 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.

Summary: 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.

Abstract 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.

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.

Abstract Background Obtaining high-quality (HQ) reference genomes from microbial communities is crucial for understanding the phylogeny and function of uncultured microbes in complex microbial ecosystems. Despite improvements in bioinformatic approaches to generate curated metagenome-assembled genomes (MAGs), existing metagenome binners obtain population consensus genomes but they are nowhere comparable to genomes sequenced from isolates in terms of strain level resolution. Here, we present a framework for the integration of single-cell genomics and metagenomics, referred to as single-cell (sc) metagenomics, to reconstruct strain-resolved genomes from microbial communities at once. Results Our sc-metagenomics integration framework, termed SMAGLinker, uses single-cell amplified genomes (SAGs) generated using microfluidic technology as binning guides and integrates them with metagenome-assembled genomes (MAGs) to recover improved draft genomes. We compared sc-metagenomics with the metagenomics-alone approach using conventional metagenome binners. The sc-metagenomics approach showed precise contig binning and higher recovery rates (>97%) of rRNA and plasmids than conventional metagenomics in genome reconstruction from the cell mock community. In human microbiota samples, sc-metagenomics recovered the largest number of genomes with a total of 103 gut microbial genomes (21 HQ, with 65 showing >90% completeness) and 45 skin microbial genomes (10 HQ, with 40 showing >90% completeness), respectively. Conventional metagenomics recovered one Staphylococcus hominis genome, whereas sc-metagenomics recovered two S. hominis genomes from identical skin microbiota sample. Single-cell sequencing revealed that these S. hominis genomes were derived from two distinct strains harboring specifically different plasmids. We found that all conventional S. hominis MAGs had a substantial lack or excess of genome sequences and contamination from other Staphylococcus species (S. epidermidis). Conclusions SMAGLinker enabled us to obtain strain-resolved genomes in the mock community and human microbiota samples by assigning metagenomic sequences correctly and covering both highly conserved genes such as rRNA genes and unique extrachromosomal elements, including plasmids. SMAGLinker will provide HQ genomes that are difficult to obtain using metagenomics alone and will facilitate the understanding of microbial ecosystems by elucidating detailed metabolic pathways and horizontal gene transfer networks. SMAGLinker is available at https://github.com/kojiari/smaglinker . Video abstract

IntroductionA 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.MethodTwo 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.ResultsUsing 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.ConclusionOur integrated spatial analysis method would be applied to various FFPE pathology specimens providing area-specific gene expression and gene mutation information.

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.

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.

Recent studies have shown that mutation at Ser522 causes inhibition of collapsin response mediator protein 2 (CRMP2) phosphorylation and induces axon elongation and partial recovery of the lost sensorimotor function after spinal cord injury (SCI). We aimed to reveal the intracellular mechanism in axotomized neurons in the CRMP2 knock-in (CRMP2KI) mouse model by performing transcriptome analysis in mouse sensorimotor cortex using micro-dissection punching system. Prior to that, we analyzed the structural pathophysiology in axotomized or neighboring neurons after SCI and found that somatic atrophy and dendritic spine reduction in sensorimotor cortex were suppressed in CRMP2KI mice. Further analysis of the transcriptome has aided in the identification of four hemoglobin genes Hba-a1, Hba-a2, Hbb-bs, and Hbb-bt that are significantly upregulated in wild-type mice with concomitant upregulation of genes involved in the oxidative phosphorylation and ribosomal pathways after SCI. However, we observed substantial upregulation in channel activity genes and downregulation of genes regulating vesicles, synaptic function, glial cell differentiation in CRMP2KI mice. Moreover, the transcriptome profile of CRMP2KI mice has been discussed wherein energy metabolism and neuronal pathways were found to be differentially regulated. Our results showed that CRMP2KI mice displayed improved SCI pathophysiology not only via microtubule stabilization in neurons, but also possibly via the whole metabolic system in the central nervous system, response changes in glial cells, and synapses. Taken together, we reveal new insights on SCI pathophysiology and the regenerative mechanism of central nervous system by the inhibition of CRMP2 phosphorylation at Ser522. All these experiments were performed in accordance with the guidelines of the Institutional Animal Care and Use Committee at Waseda University, Japan (2017-A027 approved on March 21, 2017; 2018-A003 approved on March 25, 2018; 2019-A026 approved on March 25, 2019).

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.

Asynchronicity in fungal spore germination makes transcriptomic analysis of the process challenging. Here, the authors assay single cell transcriptomes of germinating yeast cells and find that one of the histone H3 genes shows fluctuating expression, disruption of which causes germination defects.

Abstract Background The gut microbiota can have dramatic effects on host metabolism; however, current genomic strategies for uncultured bacteria have several limitations that hinder their ability to identify responders to metabolic changes in the microbiota. In this study, we describe a novel single-cell genomic sequencing technique that can identify metabolic responders at the species level without the need for reference genomes, and apply this method to identify bacterial responders to an inulin-based diet in the mouse gut microbiota. Results Inulin-feeding changed the mouse fecal microbiome composition to increase Bacteroides spp., resulting in the production of abundant succinate in the mouse intestine. Using our massively parallel single-cell genome sequencing technique, named SAG-gel platform, we obtained 346 single-amplified genomes (SAGs) from mouse gut microbes before and after dietary inulin supplementation. After quality control, the SAGs were classified as 267 bacteria, spanning 2 phyla, 4 classes, 7 orders, and 14 families, and 31 different strains of SAGs were graded as high- and medium-quality draft genomes. From these, we have successfully obtained the genomes of the dominant inulin-responders, Bacteroides spp., and identified their polysaccharide utilization loci and their specific metabolic pathways for succinate production. Conclusions Our single-cell genomics approach generated a massive amount of SAGs, enabling a functional analysis of uncultured bacteria in the intestinal microbiome. This enabled us to estimate metabolic lineages involved in the bacterial fermentation of dietary fiber and metabolic outcomes such as short-chain fatty acid production in the intestinal environment based on the fibers ingested. The technique allows the in-depth isolation and characterization of uncultured bacteria with specific functions in the microbiota and could be exploited to improve human and animal health. Video abstract.

The central nucleus of the amygdala (CeA) and the lateral division of the bed nucleus of the stria terminalis (BNST) are the two major nuclei of the central extended amygdala that plays essential roles in threat processing, responsible for emotional states such as fear and anxiety. While some studies suggested functional differences between these nuclei, others showed anatomical and neurochemical similarities. Despite their complex subnuclear organization, subnuclei-specific functional impact on behavior and their underlying molecular profiles remain obscure. We here constitutively inhibited neurotransmission of protein kinase C-δ-positive (PKCδ+) neurons—a major cell type of the lateral subdivision of the CeA (CeL) and the oval nucleus of the BNST (BNSTov)—and found striking subnuclei-specific effects on fear- and anxiety-related behaviors, respectively. To obtain molecular clues for this dissociation, we conducted RNA sequencing in subnuclei-targeted micropunch samples. The CeL and the BNSTov displayed similar gene expression profiles at the basal level; however, both displayed differential gene expression when animals were exposed to fear-related stimuli, with a more robust expression change in the CeL. These findings provide novel insights into the molecular makeup and differential engagement of distinct subnuclei of the extended amygdala, critical for regulation of threat processing.

Abstract Since G-protein coupled receptors (GPCRs) are linked to various diseases, screening of functional ligands against GPCRs is vital for drug discovery. In the present study, we developed a high-throughput functional cell-based assay by combining human culture cells producing a GPCR, yeast cells secreting randomized peptide ligands, and a droplet microfluidic device. We constructed a reporter human cell line that emits fluorescence in response to the activation of human glucagon-like peptide-1 receptor (hGLP1R). We then constructed a yeast library secreting an agonist of hGLP1R or randomized peptide ligands. We demonstrated that high-throughput identification of functional ligands against hGLP1R could be performed by co-culturing the reporter cells and the yeast cells in droplets. We identified functional ligands, one of which had higher activity than that of an original sequence. The result suggests that our system could facilitate the discovery of functional peptide ligands of GPCRs.

Abstract Objective Resting CD4+ T cells are major reservoirs of latent HIV-1 infection, and may be formed during the early phase of the infection. Although CCR5-tropic (R5) HIV-1 is highly transmissible during the early phase, newly infected individuals have usually been exposed to a mixture of R5 and CXCR4-tropic (X4) viruses, and X4 viral DNA is also detectable in the host. Our aim was to identify which subsets of resting CD4+ T cells contribute to forming the latent reservoir in the presence of both X4 and R5 viruses. Results Primary resting CD4+ naïve T (TN) cells, CCR5− memory T (TM) cells, and CCR5+ TM cells isolated by flow cytometry were infected simultaneously with X4 and R5 HIV-1, which harbored different reporter genes, and were cultured in the resting condition. Flow cytometry at 3 days post-infection demonstrated that X4 HIV-1+ cells were present in all three subsets of cells, whereas R5 HIV-1+ cells were present preferentially in CCR5+ TM cells, but not in TN cells. Following CD3/CD28-mediated activation at 3 days post-infection, numbers of R5 HIV-1+ cells and X4 HIV-1+ cells increased significantly only in the CCR5+ TM subset, suggesting that it provides a major reservoir of replication-competent, latently infected viruses.

Abstract Single-cell RNA-seq is a powerful tool for revealing heterogeneity in cancer cells. However, each of the current single-cell RNA-seq platforms has inherent advantages and disadvantages. Here, we show that combining the different single-cell RNA-seq platforms can be an effective approach to obtaining complete information about expression differences and a sufficient cellular population to understand transcriptional heterogeneity in cancers. We demonstrate that it is possible to estimate missing expression information. We further demonstrate that even in the cases where precise information for an individual gene cannot be inferred, the activity of given transcriptional modules can be analyzed. Interestingly, we found that two distinct transcriptional modules, one associated with the Aurora kinase gene and the other with the DUSP gene, are aberrantly regulated in a minor population of cells and may thus contribute to the possible emergence of dormancy or eventual drug resistance within the population.

After publication of the original article [1], the authors became aware of a miscalculation in the original Fig. 2d.

Abstract Massively parallel single-cell genome sequencing is required to further understand genetic diversities in complex biological systems. Whole genome amplification (WGA) is the first step for single-cell sequencing, but its throughput and accuracy are insufficient in conventional reaction platforms. Here, we introduce single droplet multiple displacement amplification (sd-MDA), a method that enables massively parallel amplification of single cell genomes while maintaining sequence accuracy and specificity. Tens of thousands of single cells are compartmentalized in millions of picoliter droplets and then subjected to lysis and WGA by passive droplet fusion in microfluidic channels. Because single cells are isolated in compartments, their genomes are amplified to saturation without contamination. This enables the high-throughput acquisition of contamination-free and cell specific sequence reads from single cells (21,000 single-cells/h), resulting in enhancement of the sequence data quality compared to conventional methods. This method allowed WGA of both single bacterial cells and human cancer cells. The obtained sequencing coverage rivals those of conventional techniques with superior sequence quality. In addition, we also demonstrate de novo assembly of uncultured soil bacteria and obtain draft genomes from single cell sequencing. This sd-MDA is promising for flexible and scalable use in single-cell sequencing.

Abstract Site-specific gene expression analyses are important for understanding tissue functions. Despite rapid developments in DNA-related technologies, the site-specific analysis of whole genome expression for a tissue remains challenging. Thus, a new tool is required for capturing multiple tissue micro-dissections or single cells while retaining the positional information. Here, we describe the development of such a system, which can pick up micro-dissections by punching a tissue repeatedly in a very short period, e.g., 5 s/sampling cycle. A photo of the punched tissue provides information on the dissected positions, allowing site-specific gene expression analysis. We demonstrate the site-specific analysis of a frozen tissue slice of mouse brain by analyzing many micro-dissections produced from the tissue at a 300-μm pitch. The site-specific analysis provided new insights into the gene expression profiles in a tissue and on tissue functions. The analysis of site-specific whole genome expression may therefore, open new avenues in life science.

Summary: Severe invasive group A Streptococcus (GAS) infection evades anti-bacterial immunity by attenuating the cellular components of innate immune responses. However, this loss of protection is compensated for by interferon (IFN)-γ-producing immature myeloid cells (γIMCs), which are selectively recruited upon severe invasive GAS infection in mice. Here, we demonstrate that γIMCs provide this IFN-γ-mediated protection by sequentially sensing GAS through two distinct pattern recognition receptors. In a mouse model, GAS is initially recognized by Toll-like receptor 2 (TLR2), which promptly induces interleukin (IL)-6 production in γIMCs. γIMC-derived IL-6 promotes the upregulation of a recently identified GAS-sensing receptor, macrophage-inducible C-type lectin (Mincle), in an autocrine or paracrine manner. Notably, blockade of γIMC-derived IL-6 abrogates Mincle expression, downstream IFN-γ production, and γIMC-mediated protection against severe invasive GAS infection. Thus, γIMCs regulate host protective immunity against severe invasive GAS infection via a TLR2–IL-6–Mincle axis. : Matsumura et al. show that γIMCs sequentially sense group A Streptococcus (GAS) through TLR2 and Mincle. Specifically, TLR2-triggered production of IL-6 functions as an intermediate that amplifies Mincle expression to maximize host protection through IFN-γ production. The sequential sensing is a distinct feature in γIMCs following severe invasive GAS infection. Keywords: group A Streptococcus infections, immature myeloid cells, Toll-like receptors, C-type lectin receptors, cytokines

The study of spatial distribution of secondary metabolites within microbial cells facilitates the screening of candidate strains from marine environments for functional metabolites and allows for the subsequent assessment of the production of metabolites, such as antibiotics. This paper demonstrates the first application of Raman microspectroscopy for in situ detection of the antifungal antibiotic amphotericin B (AmB) produced by actinomycetes—Streptomyces nodosus. Raman spectra measured from hyphae of S. nodosus show the specific Raman bands, caused by resonance enhancement, corresponding to the polyene chain of AmB. In addition, Raman microspectroscopy enabled us to monitor the time-dependent change of AmB production corresponding to the growth of mycelia. The Raman images of S. nodosus reveal the heterogeneous distribution of AmB within the mycelia and individual hyphae. Moreover, the molecular association state of AmB in the mycelia was directly identified by observed Raman spectral shifts. These findings suggest that Raman microspectroscopy could be used for in situ monitoring of antibiotic production directly in marine microorganisms with a method that is non-destructive and does not require labeling.

Whole genome amplification (WGA) is essential for obtaining genome sequences from single bacterial cells because the quantity of template DNA contained in a single cell is very low. Multiple displacement amplification (MDA), using Phi29 DNA polymerase and random primers, is the most widely used method for single-cell WGA. However, single-cell MDA usually results in uneven genome coverage because of amplification bias, background amplification of contaminating DNA, and formation of chimeras by linking of non-contiguous chromosomal regions. Here, we present a novel MDA method, termed droplet MDA, that minimizes amplification bias and amplification of contaminants by using picoliter-sized droplets for compartmentalized WGA reactions. Extracted DNA fragments from a lysed cell in MDA mixture are divided into 105 droplets (67 pL) within minutes via flow through simple microfluidic channels. Compartmentalized genome fragments can be individually amplified in these droplets without the risk of encounter with reagent-borne or environmental contaminants. Following quality assessment of WGA products from single Escherichia coli cells, we showed that droplet MDA minimized unexpected amplification and improved the percentage of genome recovery from 59% to 89%. Our results demonstrate that microfluidic-generated droplets show potential as an efficient tool for effective amplification of low-input DNA for single-cell genomics and greatly reduce the cost and labor investment required for determination of nearly complete genome sequences of uncultured bacteria from environmental samples.

In this paper, we present a novel cell counting method accomplished using a single-cell array fabricated on an image sensor, complementary metal oxide semiconductor sensor. The single-cell array was constructed using a microcavity array, which can trap up to 7,500 single cells on microcavities periodically arranged on a plane metallic substrate via the application of a negative pressure. The proposed method for cell counting is based on shadow imaging, which uses a light diffraction pattern generated by the microcavity array and trapped cells. Under illumination, the cell-occupied microcavities are visualized as shadow patterns in an image recorded by the complementary metal oxide semiconductor sensor due to light attenuation. The cell count is determined by enumerating the uniform shadow patterns created from one-on-one relationships with single cells trapped on the microcavities in digital format. In the experiment, all cell counting processes including entrapment of non-labeled HeLa cells from suspensions on the array and image acquisition of a wide-field-of-view of 30 mm(2) in 1/60 seconds were implemented in a single integrated device. As a result, the results from the digital cell counting had a linear relationship with those obtained from microscopic observation (r(2) = 0.99). This platform could be used at extremely low cell concentrations, i.e., 25-15,000 cells/mL. Our proposed system provides a simple and rapid miniaturized cell counting device for routine laboratory use.

Epithelial cell adhesion molecule (EpCAM)-based enumeration of circulating tumor cells (CTC) has prognostic value in patients with solid tumors, such as advanced breast, colon, and prostate cancer. However, poor sensitivity has been reported for non-small cell lung cancer (NSCLC). To address this problem, we developed a microcavity array (MCA) system integrated with a miniaturized device for CTC isolation without relying on EpCAM expression. Here, we report the results of a clinical study on CTCs of advanced lung cancer patients in which we compared the MCA system with the CellSearch system, which employs the conventional EpCAM-based method.Paired peripheral blood samples were collected from 43 metastatic lung cancer patients to enumerate CTCs using the CellSearch system according to the manufacturer's protocol and the MCA system by immunolabeling and cytomorphological analysis. The presence of CTCs was assessed blindly and independently by both systems.CTCs were detected in 17 of 22 NSCLC patients using the MCA system versus 7 of 22 patients using the CellSearch system. On the other hand, CTCs were detected in 20 of 21 small cell lung cancer (SCLC) patients using the MCA system versus 12 of 21 patients using the CellSearch system. Significantly more CTCs in NSCLC patients were detected by the MCA system (median 13, range 0-291 cells/7.5 mL) than by the CellSearch system (median 0, range 0-37 cells/7.5 ml) demonstrating statistical superiority (p = 0.0015). Statistical significance was not reached in SCLC though the trend favoring the MCA system over the CellSearch system was observed (p = 0.2888). The MCA system also isolated CTC clusters from patients who had been identified as CTC negative using the CellSearch system.The MCA system has a potential to isolate significantly more CTCs and CTC clusters in advanced lung cancer patients compared to the CellSearch system.

Despite recent advancements, it is still difficult to evaluate in vivo responses to toxicants in humans. Development of a system that can mimic the in vivo responses of human cells will enable more accurate health risk assessments. A surrogate human hematopoietic lineage can be established in NOD/Shi-scid/IL-2Rγ(null) (NOG) mice by transplanting human hematopoietic stem/progenitor cells (Hu-NOG mice). Here, we first evaluated the toxic response of human-like hematopoietic lineage in NOG mice to a representative toxic agent, benzene. Flow cytometric analysis showed that benzene caused a significant decrease in the number of human hematopoietic stem/progenitor cells in the bone marrow and the number of human leukocytes in the peripheral blood and hematopoietic organs. Next, we established chimeric mice by transplanting C57BL/6 mouse-derived bone marrow cells into NOG mice (Mo-NOG mice). A comparison of the degree of benzene-induced hematotoxicity in donor-derived hematopoietic lineage cells within Mo-NOG mice indicated that the toxic response of Hu-NOG mice reflected interspecies differences in susceptibilities to benzene. Responses to the toxic effects of benzene were greater in lymphoid cells than in myeloid cells in Mo-NOG and Hu-NOG mice. These findings suggested that Hu-NOG mice may be a powerful in vivo tool for assessing hematotoxicity in humans, while accounting for interspecies differences.