Single-cell RNA sequencing offers opportunities to explore the depth of physiology, adaptation, and biochemistry in non-model organisms exposed to pollution

Single-cell Sequencing technology (scSeq) has revolutionized our understanding of individual cells, uncovering unprecedented heterogeneity within tissues and cell populations, principality through single-cell RNA Sequencing (scRNA-Seq). This short review highlights the pivotal role of scRNA-Seq in elucidating genotype-phenotype relationships, particularly in biological systems. Based on published articles, our analysis involved manual curation and automated Scopus tools to illustrate recent advances in the application of scRNA-Seq. The results reveal that scRNA-Seq has been extensively utilized in various biological areas, including biochemistry, genetics, molecular biology, immunology, and microbiology, followed by health sciences covering studies related to the nervous system, immune system, human health, development, and diseases, with a particular focus on cancer research. However, the potential of scRNA-Seq extends beyond disease research, offering insights into non-model organisms' responses to environmental contaminants. By enabling the study of cellular reactions at a molecular level, scRNA-Seq provides a comprehensive understanding of intracellular heterogeneity that enhances our comprehension of physiological, biochemical, and pathological environmental impacts on non-model organisms exposed to pollution. This understanding has many practical benefits, as it can aid in regulation and conservation efforts that benefit the environment and the use of economically essential and ecologically relevant organisms.

Enhanced feature matching in single-cell proteomics characterizes IFN-γ response and co-existence of cell states

Proteome analysis by data-independent acquisition (DIA) has become a powerful approach to obtain deep proteome coverage, and has gained recent traction for label-free analysis of single cells. However, optimal experimental design for DIA-based single-cell proteomics has not been fully explored, and performance metrics of subsequent data analysis tools remain to be evaluated. Therefore, we here formalize and comprehensively evaluate a DIA data analysis strategy that exploits the co-analysis of low-input samples with a so-called matching enhancer (ME) of higher input, to increase sensitivity, proteome coverage, and data completeness. We assess the matching specificity of DIA-ME by a two-proteome model, and demonstrate that false discovery and false transfer are maintained at low levels when using DIA-NN software, while preserving quantification accuracy. We apply DIA-ME to investigate the proteome response of U-2 OS cells to interferon gamma (IFN-γ) in single cells, and recapitulate the time-resolved induction of IFN-γ response proteins as observed in bulk material. Moreover, we uncover co- and anti-correlating patterns of protein expression within the same cell, indicating mutually exclusive protein modules and the co-existence of different cell states. Collectively our data show that DIA-ME is a powerful, scalable, and easy-to-implement strategy for single-cell proteomics.

Deciphering deep-sea chemosynthetic symbiosis by single-nucleus RNA-sequencing

Bathymodioline mussels dominate deep-sea methane seep and hydrothermal vent habitats and obtain nutrients and energy primarily through chemosynthetic endosymbiotic bacteria in the bacteriocytes of their gill. However, the molecular mechanisms that orchestrate mussel host-symbiont interactions remain unclear. Here, we constructed a comprehensive cell atlas of the gill in the mussel Gigantidas platifrons from the South China Sea methane seeps (1100 m depth) using single-nucleus RNA-sequencing (snRNA-seq) and whole-mount in situ hybridisation. We identified 13 types of cells, including three previously unknown ones, and uncovered unknown tissue heterogeneity. Every cell type has a designated function in supporting the gill's structure and function, creating an optimal environment for chemosynthesis, and effectively acquiring nutrients from the endosymbiotic bacteria. Analysis of snRNA-seq of in situ transplanted mussels clearly showed the shifts in cell state in response to environmental oscillations. Our findings provide insight into the principles of host-symbiont interaction and the bivalves' environmental adaption mechanisms.

Scalable integration of multiomic single-cell data using generative adversarial networks

MOTIVATION

Single-cell profiling has become a common practice to investigate the complexity of tissues, organs, and organisms. Recent technological advances are expanding our capabilities to profile various molecular layers beyond the transcriptome such as, but not limited to, the genome, the epigenome, and the proteome. Depending on the experimental procedure, these data can be obtained from separate assays or the very same cells. Yet, integration of more than two assays is currently not supported by the majority of the computational frameworks avaiable.

RESULTS

We here propose a Multi-Omic data integration framework based on Wasserstein Generative Adversarial Networks suitable for the analysis of paired or unpaired data with a high number of modalities (>2). At the core of our strategy is a single network trained on all modalities together, limiting the computational burden when many molecular layers are evaluated.

AVAILABILITY AND IMPLEMENTATION

Source code of our framework is available at https://github.com/vgiansanti/MOWGAN.

The chromatin accessibility dynamics during cell fate specifications in zebrafish early embryogenesis

Chromatin accessibility plays a critical role in the regulation of cell fate decisions. Although gene expression changes have been extensively profiled at the single-cell level during early embryogenesis, the dynamics of chromatin accessibility at cis-regulatory elements remain poorly studied. Here, we used a plate-based single-cell ATAC-seq method to profile the chromatin accessibility dynamics of over 10 000 nuclei from zebrafish embryos. We investigated several important time points immediately after zygotic genome activation (ZGA), covering key developmental stages up to dome. The results revealed key chromatin signatures in the first cell fate specifications when cells start to differentiate into enveloping layer (EVL) and yolk syncytial layer (YSL) cells. Finally, we uncovered many potential cell-type specific enhancers and transcription factor motifs that are important for the cell fate specifications.

Viral vector-mediated transgene delivery with novel recombinase systems for targeting neuronal populations defined by multiple features

A comprehensive understanding of neuronal diversity and connectivity is essential for understanding the anatomical and cellular mechanisms that underlie functional contributions. With the advent of single-cell analysis, growing information regarding molecular profiles leads to the identification of more heterogeneous cell types. Therefore, the need for additional orthogonal recombinase systems is increasingly apparent, as heterogeneous tissues can be further partitioned into increasing numbers of specific cell types defined by multiple features. Critically, new recombinase systems should work together with pre-existing systems without cross-reactivity in vivo. Here, we introduce novel site-specific recombinase systems based on ΦC31 bacteriophage recombinase for labeling multiple cell types simultaneously and a novel viral strategy for versatile and robust intersectional expression of any transgene. Together, our system will help researchers specifically target different cell types with multiple features in the same animal.

Leading the pack: Best practices in comparative canine cancer genomics to inform human oncology

Pet dogs develop spontaneous cancers at a rate estimated to be five times higher than that of humans, providing a unique opportunity to study disease biology and evaluate novel therapeutic strategies in a model system that possesses an intact immune system and mirrors key aspects of human cancer biology. Despite decades of interest, effective utilization of pet dog cancers has been hindered by a limited repertoire of necessary cellular and molecular reagents for both in vitro and in vivo studies, as well as a dearth of information regarding the genomic landscape of these cancers. Recently, many of these critical gaps have been addressed through the generation of a highly annotated canine reference genome, the creation of several tools necessary for multi-omic analysis of canine tumours, and the development of a centralized repository for key genomic and associated clinical information from canine cancer patients, the Integrated Canine Data Commons. Together, these advances have catalysed multidisciplinary efforts designed to integrate the study of pet dog cancers more effectively into the translational continuum, with the ultimate goal of improving human outcomes. The current review summarizes this recent progress and provides a guide to resources and tools available for comparative study of pet dog cancers.

Mapping the castor bean endosperm proteome revealed a metabolic interaction between plastid, mitochondria, and peroxisomes to optimize seedling growth

In this work, we studied castor-oil plant Ricinus communis as a classical system for endosperm reserve breakdown. The seeds of castor beans consist of a centrally located embryo with the two thin cotyledons surrounded by the endosperm. The endosperm functions as major storage tissue and is packed with nutritional reserves, such as oil, proteins, and starch. Upon germination, mobilization of the storage reserves requires inter-organellar interplay of plastids, mitochondria, and peroxisomes to optimize growth for the developing seedling. To understand their metabolic interactions, we performed a large-scale organellar proteomic study on castor bean endosperm. Organelles from endosperm of etiolated seedlings were isolated and subjected to liquid chromatography-tandem mass spectrometry (LC-MS/MS). Computer-assisted deconvolution algorithms were applied to reliably assign the identified proteins to their correct subcellular localization and to determine the abundance of the different organelles in the heterogeneous protein samples. The data obtained were used to build a comprehensive metabolic model for plastids, mitochondria, and peroxisomes during storage reserve mobilization in castor bean endosperm.

Modeling group heteroscedasticity in single-cell RNA-seq pseudo-bulk data

Group heteroscedasticity is commonly observed in pseudo-bulk single-cell RNA-seq datasets and its presence can hamper the detection of differentially expressed genes. Since most bulk RNA-seq methods assume equal group variances, we introduce two new approaches that account for heteroscedastic groups, namely voomByGroup and voomWithQualityWeights using a blocked design (voomQWB). Compared to current gold-standard methods that do not account for group heteroscedasticity, we show results from simulations and various experiments that demonstrate the superior performance of voomByGroup and voomQWB in terms of error control and power when group variances in pseudo-bulk single-cell RNA-seq data are unequal.

Ultrasound frequency sonication facilitates high-throughput and uniform dissociation of cellular aggregates and tissues

A sample preparation step involving dissociation of tissues into their component cells is often required to conduct analysis of nucleic acids and other constituents from tissue samples. Frequently, the extracellular matrix and cell-cell adhesions are disrupted via treatment with a chemical dissociating reagent or various mechanical forces. In this work, a new, high-throughput, multiplexed method of dissociating tissues and cellular aggregates into single cells using ultrasound frequency bath sonication is explored and characterized. Different operating parameters are evaluated, and a treatment protocol with potential for uniform, high-throughput tissue dissociation is compared to the existing best chemical and orbital plate shaking protocol. Metrics such as percent dissociation, cellular recovery, average aggregate size, proportion of various aggregate sizes, membrane circularity, and cellular viability are subsequently assessed and found to be favorable. In optimized conditions, 53 ± 8% of 1 mm biopsy cores are dissociated within 30 min using sonication alone, surpassing leading high-throughput orbital plate shaking techniques five-fold. Chemical digestion is also 2 times more effective when complexed with sonication rather than orbital plate shaking. RNA content, quality, and expression are found to be superior to the standard protocol in terms of transcriptional preservation.

SPICEMIX enables integrative single-cell spatial modeling of cell identity

Spatial transcriptomics can reveal spatially resolved gene expression of diverse cells in complex tissues. However, the development of computational methods that can use the unique properties of spatial transcriptome data to unveil cell identities remains a challenge. Here we introduce SPICEMIX, an interpretable method based on probabilistic, latent variable modeling for joint analysis of spatial information and gene expression from spatial transcriptome data. Both simulation and real data evaluations demonstrate that SPICEMIX markedly improves on the inference of cell types and their spatial patterns compared with existing approaches. By applying to spatial transcriptome data of brain regions in human and mouse acquired by seqFISH+, STARmap and Visium, we show that SPICEMIX can enhance the inference of complex cell identities, reveal interpretable spatial metagenes and uncover differentiation trajectories. SPICEMIX is a generalizable analysis framework for spatial transcriptome data to investigate cell-type composition and spatial organization of cells in complex tissues.

HT-smFISH: a cost-effective and flexible workflow for high-throughput single-molecule RNA imaging

The ability to visualize RNA in its native subcellular environment by using single-molecule fluorescence in situ hybridization (smFISH) has reshaped our understanding of gene expression and cellular functions. A major hindrance of smFISH is the difficulty to perform systematic experiments in medium- or high-throughput formats, principally because of the high cost of generating the individual fluorescent probe sets. Here, we present high-throughput smFISH (HT-smFISH), a simple and cost-efficient method for imaging hundreds to thousands of single endogenous RNA molecules in 96-well plates. HT-smFISH uses RNA probes transcribed in vitro from a large pool of unlabeled oligonucleotides. This allows the generation of individual probes for many RNA species, replacing commercial DNA probe sets. HT-smFISH thus reduces costs per targeted RNA compared with many smFISH methods and is easily scalable and flexible in design. We provide a protocol that combines oligo pool design, probe set generation, optimized hybridization conditions and guidelines for image acquisition and analysis. The pipeline requires knowledge of standard molecular biology tools, cell culture and fluorescence microscopy. It is achievable in ~20 d. In brief, HT-smFISH is tailored for medium- to high-throughput screens that image RNAs at single-molecule sensitivity.

Single-cell RNA-seq analysis to identify potential biomarkers for diagnosis, and prognosis of non-small cell lung cancer by using comprehensive bioinformatics approaches

Non-small cell lung cancer (NSCLC) is the most common type of lung cancer and the leading cause of cancer-related deaths worldwide. Identification of gene biomarkers and their regulatory factors and signaling pathways is very essential to reveal the molecular mechanisms of NSCLC initiation and progression. Thus, the goal of this study is to identify gene biomarkers for NSCLC diagnosis and prognosis by using scRNA-seq data through bioinformatics techniques. scRNA-seq data were obtained from the GEO database to identify DEGs. A total of 158 DEGs (including 48 upregulated and 110 downregulated) were detected after gene integration. Gene Ontology enrichment and KEGG pathway analysis of DEGs were performed by FunRich software. A PPI network of DEGs was then constructed using the STRING database and visualized by Cytoscape software. We identified 12 key genes (KGs) including MS4A1, CCL5, and GZMB, by using two topological methods based on the PPI networking results. The diagnostic, expression, and prognostic potentials of the identified 12 key genes were assessed using the receiver operating characteristics (ROC) curve and a web-based tool, SurvExpress. From the regulatory network analysis, we extracted the 7 key transcription factors (TFs) (FOXC1, YY1, CEBPB, TFAP2A, SREBF2, RELA, and GATA2), and 8 key miRNAs (hsa-miR-124-3p, hsa-miR-34a-5p, hsa-miR-21-5p, hsa-miR-155-5p, hsa-miR-449a, hsa-miR-24-3p, hsa-let-7b-5p, and hsa-miR-7-5p) associated with the KGs were evaluated. Functional enrichment and pathway analysis, survival analysis, ROC analysis, and regulatory network analysis highlighted crucial roles of the key genes. Our findings might play a significant role as candidate biomarkers in NSCLC diagnosis and prognosis.

Aquila: a spatial omics database and analysis platform

Spatial omics is a rapidly evolving approach for exploring tissue microenvironment and cellular networks by integrating spatial knowledge with transcript or protein expression information. However, there is a lack of databases for users to access and analyze spatial omics data. To address this limitation, we developed Aquila, a comprehensive platform for managing and analyzing spatial omics data. Aquila contains 107 datasets from 30 diseases, including 6500+ regions of interest, and 15.7 million cells. The database covers studies from spatial transcriptome and proteome analyses, 2D and 3D experiments, and different technologies. Aquila provides visualization of spatial omics data in multiple formats such as spatial cell distribution, spatial expression and co-localization of markers. Aquila also lets users perform many basic and advanced spatial analyses on any dataset. In addition, users can submit their own spatial omics data for visualization and analysis in a safe and secure environment. Finally, Aquila can be installed as an individual app on a desktop and offers the RESTful API service for power users to access the database. Overall, Aquila provides a detailed insight into transcript and protein expression in tissues from a spatial perspective. Aquila is available at https://aquila.cheunglab.org.

Omnibus and robust deconvolution scheme for bulk RNA sequencing data integrating multiple single-cell reference sets and prior biological knowledge

MOTIVATION

Cell-type deconvolution of bulk tissue RNA sequencing (RNA-seq) data is an important step toward understanding the variations in cell-type composition among disease conditions. Owing to recent advances in single-cell RNA sequencing (scRNA-seq) and the availability of large amounts of bulk RNA-seq data in disease-relevant tissues, various deconvolution methods have been developed. However, the performance of existing methods heavily relies on the quality of information provided by external data sources, such as the selection of scRNA-seq data as a reference and prior biological information.

RESULTS

We present the Integrated and Robust Deconvolution (InteRD) algorithm to infer cell-type proportions from target bulk RNA-seq data. Owing to the innovative use of penalized regression with a new evaluation criterion for deconvolution, InteRD has three primary advantages. First, it is able to effectively integrate deconvolution results from multiple scRNA-seq datasets. Second, InteRD calibrates estimates from reference-based deconvolution by taking into account extra biological information as priors. Third, the proposed algorithm is robust to inaccurate external information imposed in the deconvolution system. Extensive numerical evaluations and real-data applications demonstrate that InteRD yields more accurate and robust cell-type proportion estimates that agree well with known biology.

AVAILABILITY AND IMPLEMENTATION

The proposed InteRD framework is implemented in R and the package is available at https://cran.r-project.org/web/packages/InteRD/index.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Investigation of Brassica and its relative genomes in the post-genomics era

The Brassicaceae family includes many economically important crop species, as well as cosmopolitan agricultural weed species. In addition, Arabidopsis thaliana, a member of this family, is used as a molecular model plant species. The genus Brassica is mesopolyploid, and the genus comprises comparatively recently originated tetrapolyploid species. With these characteristics, Brassicas have achieved the commonly accepted status of model organisms for genomic studies. This paper reviews the rapid research progress in the Brassicaceae family from diverse omics studies, including genomics, transcriptomics, epigenomics, and three-dimensional (3D) genomics, with a focus on cultivated crops. The morphological plasticity of Brassicaceae crops is largely due to their highly variable genomes. The origin of several important Brassicaceae crops has been established. Genes or loci domesticated or contributing to important traits are summarized. Epigenetic alterations and 3D structures have been found to play roles in subgenome dominance, either in tetraploid Brassica species or their diploid ancestors. Based on this progress, we propose future directions and prospects for the genomic investigation of Brassicaceae crops.

Well-Paired-Seq: A Size-Exclusion and Locally Quasi-Static Hydrodynamic Microwell Chip for Single-Cell RNA-Seq

Single-cell RNA sequencing (scRNA-seq) is a powerful technology for revealing the heterogeneity of cellular states. However, existing scRNA-seq platforms that utilize bead-based technologies suffer from a large number of empty microreactors and a low cell/bead capture efficiency. Here, Well-paired-seq is presented, which consists of thousands of size exclusion and quasi-static hydrodynamic dual wells to address these limitations. The size-exclusion principle allows one cell and one bead to be trapped in the bottom well (cell-capture-well) and the top well (bead-capture-well), respectively, while the quasi-static hydrodynamic principle ensures that the trapped cells are difficult to escape from cell-capture-wells, achieving cumulative capture of cells and effective buffer exchange. By the integration of quasi-static hydrodynamic and size-exclusion principles, the dual wells ensure single cells/beads pairing with high density, achieving excellent efficiency of cell capture (≈91%), cell/bead pairing (≈82%), and cell-free RNA removal. The high utilization of microreactors and single cells/beads enable to achieve a high throughput (≈10⁵ cells) with low collision rates. The technical performance of Well-paired-seq is demonstrated by collecting transcriptome data from around 200 000 cells across 21 samples, successfully revealing the heterogeneity of single cells and showing the wide applicability of Well-paired-seq for basic and clinical research.

Phenotype Identification of HeLa Cells Knockout CDK6 Gene Based on Label-Free Raman Imaging

Identifying cell phenotypes is essential for understanding the function of biological macromolecules and molecular biology. We developed a noninvasive, label-free, single-cell Raman imaging analysis platform to distinguish between the cell phenotypes of the HeLa cell wild type (WT) and cyclin-dependent kinase 6 (CDK6) gene knockout (KO) type. Via large-scale Raman spectral and imaging analysis, two phenotypes of the HeLa cells were distinguished by their intrinsic biochemical profiles. A significant difference was found between the two cell lines: large lipid droplets formed in the knockout HeLa cells but were not observed in the WT cells, which was confirmed by Oil Red O staining. The band ratio of the Raman spectrum of saturated/unsaturated fatty acids was identified as the Raman spectral marker for HeLa cell WT or gene knockout type differentiation. The interaction between organelles involved in lipid metabolism was revealed by Raman imaging and Lorentz fitting, where the distribution intensity of the mitochondria and the endoplasmic reticulum membrane decreased. At the same time, lysosomes increased after the CDK6 gene knockout. The parameters obtained from Raman spectroscopy are based on hierarchical cluster analysis and one-way ANOVA, enabling highly accurate cell classification.

Comparison of methods and resources for cell-cell communication inference from single-cell RNA-Seq data

The growing availability of single-cell data, especially transcriptomics, has sparked an increased interest in the inference of cell-cell communication. Many computational tools were developed for this purpose. Each of them consists of a resource of intercellular interactions prior knowledge and a method to predict potential cell-cell communication events. Yet the impact of the choice of resource and method on the resulting predictions is largely unknown. To shed light on this, we systematically compare 16 cell-cell communication inference resources and 7 methods, plus the consensus between the methods' predictions. Among the resources, we find few unique interactions, a varying degree of overlap, and an uneven coverage of specific pathways and tissue-enriched proteins. We then examine all possible combinations of methods and resources and show that both strongly influence the predicted intercellular interactions. Finally, we assess the agreement of cell-cell communication methods with spatial colocalisation, cytokine activities, and receptor protein abundance and find that predictions are generally coherent with those data modalities. To facilitate the use of the methods and resources described in this work, we provide LIANA, a LIgand-receptor ANalysis frAmework as an open-source interface to all the resources and methods.

Explainable multiview framework for dissecting spatial relationships from highly multiplexed data

The advancement of highly multiplexed spatial technologies requires scalable methods that can leverage spatial information. We present MISTy, a flexible, scalable, and explainable machine learning framework for extracting relationships from any spatial omics data, from dozens to thousands of measured markers. MISTy builds multiple views focusing on different spatial or functional contexts to dissect different effects. We evaluated MISTy on in silico and breast cancer datasets measured by imaging mass cytometry and spatial transcriptomics. We estimated structural and functional interactions coming from different spatial contexts in breast cancer and demonstrated how to relate MISTy's results to clinical features.

Behavioral Neuroscience in the Era of Genomics: Tools and Lessons for Analyzing High-Dimensional Datasets

Behavioral neuroscience underwent a technology-driven revolution with the emergence of machine-vision and machine-learning technologies. These technological advances facilitated the generation of high-resolution, high-throughput capture and analysis of complex behaviors. Therefore, behavioral neuroscience is becoming a data-rich field. While behavioral researchers use advanced computational tools to analyze the resulting datasets, the search for robust and standardized analysis tools is still ongoing. At the same time, the field of genomics exploded with a plethora of technologies which enabled the generation of massive datasets. This growth of genomics data drove the emergence of powerful computational approaches to analyze these data. Here, we discuss the composition of a large behavioral dataset, and the differences and similarities between behavioral and genomics data. We then give examples of genomics-related tools that might be of use for behavioral analysis and discuss concepts that might emerge when considering the two fields together.

Question-driven stepwise experimental discoveries in biochemistry: two case studies

Philosophers of science diverge on the question what drives the growth of scientific knowledge. Most of the twentieth century was dominated by the notion that theories propel that growth whereas experiments play secondary roles of operating within the theoretical framework or testing theoretical predictions. New experimentalism, a school of thought pioneered by Ian Hacking in the early 1980s, challenged this view by arguing that theory-free exploratory experimentation may in many cases effectively probe nature and potentially spawn higher evidence-based theories. Because theories are often powerless to envisage workings of complex biological systems, theory-independent experimentation is common in the life sciences. Some such experiments are triggered by compelling observation, others are prompted by innovative techniques or instruments, whereas different investigations query big data to identify regularities and underlying organizing principles. A distinct fourth type of experiments is motivated by a major question. Here I describe two question-guided experimental discoveries in biochemistry: the cyclic adenosine monophosphate mediator of hormone action and the ubiquitin-mediated system of protein degradation. Lacking underlying theories, antecedent data bases, or new techniques, the sole guides of the two discoveries were respective substantial questions. Both research projects were similarly instigated by theory-free exploratory experimentation and continued in alternating phases of results-based interim working hypotheses, their examination by experiment, provisional hypotheses again, and so on. These two cases designate theory-free, question-guided, stepwise biochemical investigations as a distinct subtype of the new experimentalism mode of scientific enquiry.

Characterization of transcript enrichment and detection bias in single-nucleus RNA-seq for mapping of distinct human adipocyte lineages

Single-cell RNA sequencing (scRNA-seq) enables molecular characterization of complex biological tissues at high resolution. The requirement of single-cell extraction, however, makes it challenging for profiling tissues such as adipose tissue, for which collection of intact single adipocytes is complicated by their fragile nature. For such tissues, single-nucleus extraction is often much more efficient and therefore single-nucleus RNA sequencing (snRNA-seq) presents an alternative to scRNA-seq. However, nuclear transcripts represent only a fraction of the transcriptome in a single cell, with snRNA-seq marked with inherent transcript enrichment and detection biases. Therefore, snRNA-seq may be inadequate for mapping important transcriptional signatures in adipose tissue. In this study, we compare the transcriptomic landscape of single nuclei isolated from preadipocytes and mature adipocytes across human white and brown adipocyte lineages, with whole-cell transcriptome. We show that snRNA-seq is capable of identifying the broad cell types present in scRNA-seq at all states of adipogenesis. However, we also explore how and why the nuclear transcriptome is biased and limited, as well as how it can be advantageous. We robustly characterize the enrichment of nuclear-localized transcripts and adipogenic regulatory lncRNAs in snRNA-seq, while also providing a detailed understanding for the preferential detection of long genes upon using this technique. To remove such technical detection biases, we propose a normalization strategy for a more accurate comparison of nuclear and cellular data. Finally, we show successful integration of scRNA-seq and snRNA-seq data sets with existing bioinformatic tools. Overall, our results illustrate the applicability of snRNA-seq for the characterization of cellular diversity in the adipose tissue.

Capturing the third dimension in drug discovery: Spatially-resolved tools for interrogation of complex 3D cell models

Advanced three-dimensional (3D) cell models have proven to be capable of depicting architectural and microenvironmental features of several tissues. By providing data of higher physiological and pathophysiological relevance, 3D cell models have been contributing to a better understanding of human development, pathology onset and progression mechanisms, as well as for 3D cell-based assays for drug discovery. Nonetheless, the characterization and interrogation of these tissue-like structures pose major challenges on the conventional analytical methods, pushing the development of spatially-resolved technologies. Herein, we review recent advances and pioneering technologies suitable for the interrogation of multicellular 3D models, while capable of retaining biological spatial information. We focused on imaging technologies and omics tools, namely transcriptomics, proteomics and metabolomics. The advantages and shortcomings of these novel methodologies are discussed, alongside the opportunities to intertwine data from the different tools.

A Modified SMART-Seq Method for Single-Cell Transcriptomic Analysis of Embryoid Body Differentiation

Embryoid bodies (EBs) are aggregate of cells that contain three embryonic germ layers. They can be formed by direct differentiation from pluripotent embryonic stem cells (ESCs), which serves as a useful model for understanding early embryo development. Due to the mixture of different cell types, it is necessary to investigate EBs at the single-cell level. Here, we describe a robust and straightforward method for single-cell gene expression profiling during mouse EB differentiation from mouse ESCs (mESCs). The protocol is modified from a widely used method in the SMART-seq family, which only requires standard molecular biology techniques and lab equipment. It allows for accurate 3' counting of transcript at the single-cell level, which helps reveal cellular identities during EB formation. Combined with perturbation experiments, the method provides an opportunity for mechanistic studies of embryo development at the single-cell level.

The Stochastic Genome and Its Role in Gene Expression

Mammalian genomes have distinct levels of spatial organization and structure that have been hypothesized to play important roles in transcription regulation. Although much has been learned about these architectural features with ensemble techniques, single-cell studies are showing a new universal trend: Genomes are stochastic and dynamic at every level of organization. Stochastic gene expression, on the other hand, has been studied for years. In this review, we probe whether there is a causative link between the two phenomena. We specifically discuss the functionality of chromatin state, topologically associating domains (TADs), and enhancer biology in light of their stochastic nature and their specific roles in stochastic gene expression. We highlight persistent fundamental questions in this area of research.

Confronting false discoveries in single-cell differential expression

Differential expression analysis in single-cell transcriptomics enables the dissection of cell-type-specific responses to perturbations such as disease, trauma, or experimental manipulations. While many statistical methods are available to identify differentially expressed genes, the principles that distinguish these methods and their performance remain unclear. Here, we show that the relative performance of these methods is contingent on their ability to account for variation between biological replicates. Methods that ignore this inevitable variation are biased and prone to false discoveries. Indeed, the most widely used methods can discover hundreds of differentially expressed genes in the absence of biological differences. To exemplify these principles, we exposed true and false discoveries of differentially expressed genes in the injured mouse spinal cord.

Cross-platform transcriptional profiling identifies common and distinct molecular pathologies in Lewy body diseases

Parkinson's disease (PD), Parkinson's disease with dementia (PDD) and dementia with Lewy bodies (DLB) are three clinically, genetically and neuropathologically overlapping neurodegenerative diseases collectively known as the Lewy body diseases (LBDs). A variety of molecular mechanisms have been implicated in PD pathogenesis, but the mechanisms underlying PDD and DLB remain largely unknown, a knowledge gap that presents an impediment to the discovery of disease-modifying therapies. Transcriptomic profiling can contribute to addressing this gap, but remains limited in the LBDs. Here, we applied paired bulk-tissue and single-nucleus RNA-sequencing to anterior cingulate cortex samples derived from 28 individuals, including healthy controls, PD, PDD and DLB cases (n = 7 per group), to transcriptomically profile the LBDs. Using this approach, we (i) found transcriptional alterations in multiple cell types across the LBDs; (ii) discovered evidence for widespread dysregulation of RNA splicing, particularly in PDD and DLB; (iii) identified potential splicing factors, with links to other dementia-related neurodegenerative diseases, coordinating this dysregulation; and (iv) identified transcriptomic commonalities and distinctions between the LBDs that inform understanding of the relationships between these three clinical disorders. Together, these findings have important implications for the design of RNA-targeted therapies for these diseases and highlight a potential molecular "window" of therapeutic opportunity between the initial onset of PD and subsequent development of Lewy body dementia.

Cell Type-Specific Expression Analysis of the Inner Ear: A Technical Report

OBJECTIVE

The cellular diversity of the inner ear has presented a technical challenge in obtaining molecular insight into its development and function. The application of technological advancements in cell type-specific expression enable clinicians and researchers to leap forward from classic genetics to obtaining mechanistic understanding of congenital and acquired hearing loss. This understanding is essential for development of therapeutics to prevent and reverse diseases of the inner ear, including hearing loss. The objective of this study is to describe and compare the available tools for cell type-specific analysis of the ear, as a means to support decision making in study design.

STUDY DESIGN

Three major approaches for cell type-specific analysis of the ear including fluorescence-activated cell sorting (FACS), ribosomal and RNA pulldown techniques, and single cell RNA-seq (scRNA-seq) are compared and contrasted using both published and original data.

RESULTS

We demonstrate the strength and weaknesses of these approaches leading to the inevitable conclusion that to maximize the utility of these approaches, it is important to match the experimental approach with the tissue of origin, cell type of interest, and the biological question. Often, a combined approach (eg, cell sorting and scRNA-seq or expression analysis using 2 separate approaches) is required. Finally, new tools for visualization and analysis of complex expression data, such as the gEAR platform (umgear.org), collate cell type-specific gene expression from the ear field and provide unprecedented access to both clinicians and researchers.

LEVEL OF EVIDENCE

N/A Laryngoscope, 131:S1-S16, 2021.

RNA imaging in bacteria

The spatiotemporal regulation of gene expression plays an essential role in many biological processes. Recently, several imaging-based RNA labeling and detection methods, both in fixed and live cells, were developed and now enable the study of transcript abundance, localization and dynamics. Here, we review the main single-cell techniques for RNA visualization with fluorescence microscopy and describe their applications in bacteria.

Diversity of Reactive Astrogliosis in CNS Pathology: Heterogeneity or Plasticity?

Astrocytes are essential for the development and homeostatic maintenance of the central nervous system (CNS). They are also critical players in the CNS injury response during which they undergo a process referred to as "reactive astrogliosis." Diversity in astrocyte morphology and gene expression, as revealed by transcriptional analysis, is well-recognized and has been reported in several CNS pathologies, including ischemic stroke, CNS demyelination, and traumatic injury. This diversity appears unique to the specific pathology, with significant variance across temporal, topographical, age, and sex-specific variables. Despite this, there is limited functional data corroborating this diversity. Furthermore, as reactive astrocytes display significant environmental-dependent plasticity and fate-mapping data on astrocyte subsets in the adult CNS is limited, it remains unclear whether this diversity represents heterogeneity or plasticity. As astrocytes are important for neuronal survival and CNS function post-injury, establishing to what extent this diversity reflects distinct established heterogeneous astrocyte subpopulations vs. environmentally dependent plasticity within established astrocyte subsets will be critical for guiding therapeutic development. To that end, we review the current state of knowledge on astrocyte diversity in the context of three representative CNS pathologies: ischemic stroke, demyelination, and traumatic injury, with the goal of identifying key limitations in our current knowledge and suggesting future areas of research needed to address them. We suggest that the majority of identified astrocyte diversity in CNS pathologies to date represents plasticity in response to dynamically changing post-injury environments as opposed to heterogeneity, an important consideration for the understanding of disease pathogenesis and the development of therapeutic interventions.

Seq-ing out cell types across the isocortex and hippocampal formation

A central quest in neuroscience is to gain a holistic understanding of all cell types in the brain. In this issue of Cell, Yao et al. establish a molecular architectural view of cell types across the entire adult mouse isocortex and hippocampal formation and reveal surprising similarities of cell types in these two brain regions.

Cryopreservation Preserves Cell-Type Composition and Gene Expression Profiles in Bone Marrow Aspirates From Multiple Myeloma Patients

Single-cell RNA sequencing reveals gene expression differences between individual cells and also identifies different cell populations that are present in the bulk starting material. To obtain an accurate assessment of patient samples, single-cell suspensions need to be generated as soon as possible once the tissue or sample has been collected. However, this requirement poses logistical challenges for experimental designs involving multiple samples from the same subject since these samples would ideally be processed at the same time to minimize technical variation in data analysis. Although cryopreservation has been shown to largely preserve the transcriptome, it is unclear whether the freeze-thaw process might alter gene expression profiles in a cell-type specific manner or whether changes in cell-type proportions might also occur. To address these questions in the context of multiple myeloma clinical studies, we performed single-cell RNA sequencing (scRNA-seq) to compare fresh and frozen cells isolated from bone marrow aspirates of six multiple myeloma patients, analyzing both myeloma cells (CD138⁺) and cells constituting the microenvironment (CD138−). We found that cryopreservation using 90% fetal calf serum and 10% dimethyl sulfoxide resulted in highly consistent gene expression profiles when comparing fresh and frozen samples from the same patient for both CD138⁺ myeloma cells (R ≥ 0.96) and for CD138^– cells (R ≥ 0.9). We also demonstrate that CD138^– cell-type proportions showed minimal alterations, which were mainly related to small differences in immune cell subtype sensitivity to the freeze-thaw procedures. Therefore, when processing fresh multiple myeloma samples is not feasible, cryopreservation is a useful option in single-cell profiling studies.

Two reactive behaviors of chondrocytes in an IL-1β-induced inflammatory environment revealed by the single-cell RNA sequencing

Objective: To investigate the heterogeneous responses of in vitro expanded chondrocytes, which were cultured in an interleukin (IL)-1β -induced inflammatory environment.

Method: Human articular chondrocytes were expanded, in vitro, for 13 days and treated with IL-1β for 0, 24, and 48 h. Cells were collected and subjected to single-cell RNA sequencing. Multiple bioinformatics tools were used to determine the signatures that define chondrocyte physiology.

Results: Two major cell clusters with distinct expression patterns were identified at the initial phase and were with heterogeneous variation that coincides with inflammation progress. They transformed into two terminal cell clusters one of which exhibited OA-phenotype and proinflammatory characteristics through two paths, “response-to-inflammation” and “atypical response-to-inflammation”, respectively. The involved cell clusters exhibited intrinsic relationship with cell types within native cartilage from OA patients. Genes controlling cell transformation to OA-phenotype were relating to the tumor necrosis factor (TNF) signaling pathway via NFKB, up-regulated KRAS signaling and the IL2/STAT5 signaling pathway and pathways relating to apoptosis and reactive oxygen species.

Conclusion: The in vitro expanded chondrocytes under IL-1β-induced inflammatory progression behave heterogeneously. One of the initial cell clusters could transform into a proinflammatory subpopulation through a termed response-to-inflammation path, which may serve as the core target to alleviate OA progression.

Novel Technologies in Studying Brain Immune Response

Over the past few decades, the immune system, including both the adaptive and innate immune systems, proved to be essential and critical to brain damage and recovery in the pathogenesis of several diseases, opening a new avenue for developing new immunomodulatory therapies and novel treatments for many neurological diseases. However, due to the specificity and structural complexity of the central nervous system (CNS), and the limit of the related technologies, the biology of the immune response in the brain is still poorly understood. Here, we discuss the application of novel technologies in studying the brain immune response, including single-cell RNA analysis, cytometry by time-of-flight, and whole-genome transcriptomic and proteomic analysis. We believe that advancements in technology related to immune research will provide an optimistic future for brain repair.

FlsnRNA-seq: protoplasting-free full-length single-nucleus RNA profiling in plants

The broad application of single-cell RNA profiling in plants has been hindered by the prerequisite of protoplasting that requires digesting the cell walls from different types of plant tissues. Here, we present a protoplasting-free approach, flsnRNA-seq, for large-scale full-length RNA profiling at a single-nucleus level in plants using isolated nuclei. Combined with 10x Genomics and Nanopore long-read sequencing, we validate the robustness of this approach in Arabidopsis root cells and the developing endosperm. Sequencing results demonstrate that it allows for uncovering alternative splicing and polyadenylation-related RNA isoform information at the single-cell level, which facilitates characterizing cell identities.

A hitchhiker's guide to single-cell transcriptomics and data analysis pipelines

Single-cell transcriptomics (SCT) is a tour de force in the era of big omics data that has led to the accumulation of massive cellular transcription data at an astounding resolution of single cells. It provides valuable insights into cells previously unachieved by bulk cell analysis and is proving crucial in uncovering cellular heterogeneity, identifying rare cell populations, distinct cell-lineage trajectories, and mechanisms involved in complex cellular processes. SCT data is highly complex and necessitates advanced statistical and computational methods for analysis. This review provides a comprehensive overview of the steps in a typical SCT workflow, starting from experimental protocol to data analysis, deliberating various pipelines used. We discuss recent trends, challenges, machine learning methods for data analysis, and future prospects. We conclude by listing the multitude of scRNA-seq data applications and how it shall revolutionize our understanding of cellular biology and diseases.

Meta-analysis of Caenorhabditis elegans single-cell developmental data reveals multi-frequency oscillation in gene activation

MOTIVATION

The advent of in vivo automated techniques for single-cell lineaging, sequencing and analysis of gene expression has begun to dramatically increase our understanding of organismal development. We applied novel meta-analysis and visualization techniques to the EPIC single-cell-resolution developmental gene expression dataset for Caenorhabditis elegans from Bao, Murray, Waterston et al. to gain insights into regulatory mechanisms governing the timing of development.

RESULTS

Our meta-analysis of the EPIC dataset revealed that a simple linear combination of the expression levels of the developmental genes is strongly correlated with the developmental age of the organism, irrespective of the cell division rate of different cell lineages. We uncovered a pattern of collective sinusoidal oscillation in gene activation, in multiple dominant frequencies and in multiple orthogonal axes of gene expression, pointing to the existence of a coordinated, multi-frequency global timing mechanism. We developed a novel method based on Fisher's Discriminant Analysis to identify gene expression weightings that maximally separate traits of interest, and found that remarkably, simple linear gene expression weightings are capable of producing sinusoidal oscillations of any frequency and phase, adding to the growing body of evidence that oscillatory mechanisms likely play an important role in the timing of development. We cross-linked EPIC with gene ontology and anatomy ontology terms, employing Fisher's Discriminant Analysis methods to identify previously unknown positive and negative genetic contributions to developmental processes and cell phenotypes. This meta-analysis demonstrates new evidence for direct linear and/or sinusoidal mechanisms regulating the timing of development. We uncovered a number of previously unknown positive and negative correlations between developmental genes and developmental processes or cell phenotypes. Our results highlight both the continued relevance of the EPIC technique, and the value of meta-analysis of previously published results. The presented analysis and visualization techniques are broadly applicable across developmental and systems biology.

AVAILABILITY AND IMPLEMENTATION

Analysis software available upon request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Systematic Evaluation of Interneuron Morphology Representations for Cell Type Discrimination

Quantitative analysis of neuronal morphologies usually begins with choosing a particular feature representation in order to make individual morphologies amenable to standard statistics tools and machine learning algorithms. Many different feature representations have been suggested in the literature, ranging from density maps to intersection profiles, but they have never been compared side by side. Here we performed a systematic comparison of various representations, measuring how well they were able to capture the difference between known morphological cell types. For our benchmarking effort, we used several curated data sets consisting of mouse retinal bipolar cells and cortical inhibitory neurons. We found that the best performing feature representations were two-dimensional density maps, two-dimensional persistence images and morphometric statistics, which continued to perform well even when neurons were only partially traced. Combining these feature representations together led to further performance increases suggesting that they captured non-redundant information. The same representations performed well in an unsupervised setting, implying that they can be suitable for dimensionality reduction or clustering.

Drosophila PDGF/VEGF signaling from muscles to hepatocyte-like cells protects against obesity

PDGF/VEGF ligands regulate a plethora of biological processes in multicellular organisms via autocrine, paracrine, and endocrine mechanisms. We investigated organ-specific metabolic roles of Drosophila PDGF/VEGF-like factors (Pvfs). We combine genetic approaches and single-nuclei sequencing to demonstrate that muscle-derived Pvf1 signals to the Drosophila hepatocyte-like cells/oenocytes to suppress lipid synthesis by activating the Pi3K/Akt1/TOR signaling cascade in the oenocytes. Functionally, this signaling axis regulates expansion of adipose tissue lipid stores in newly eclosed flies. Flies emerge after pupation with limited adipose tissue lipid stores and lipid level is progressively accumulated via lipid synthesis. We find that adult muscle-specific expression of pvf1 increases rapidly during this stage and that muscle-to-oenocyte Pvf1 signaling inhibits expansion of adipose tissue lipid stores as the process reaches completion. Our findings provide the first evidence in a metazoan of a PDGF/VEGF ligand acting as a myokine that regulates systemic lipid homeostasis by activating TOR in hepatocyte-like cells.

Probing infectious disease by single-cell RNA sequencing: Progresses and perspectives

The increasing application of single-cell RNA sequencing (scRNA-seq) technology in life science and biomedical research has significantly increased our understanding of the cellular heterogeneities in immunology, oncology and developmental biology. This review will summarize the development of various scRNA-seq technologies; primarily discussing the application of scRNA-seq on infectious diseases, and exploring the current development, challenges, and potential applications of scRNA-seq technology in the future.

bayNorm: Bayesian gene expression recovery, imputation and normalization for single-cell RNA-sequencing data

Motivation

Normalization of single-cell RNA-sequencing (scRNA-seq) data is a prerequisite to their interpretation. The marked technical variability, high amounts of missing observations and batch effect typical of scRNA-seq datasets make this task particularly challenging. There is a need for an efficient and unified approach for normalization, imputation and batch effect correction.

Results

Here, we introduce bayNorm, a novel Bayesian approach for scaling and inference of scRNA-seq counts. The method’s likelihood function follows a binomial model of mRNA capture, while priors are estimated from expression values across cells using an empirical Bayes approach. We first validate our assumptions by showing this model can reproduce different statistics observed in real scRNA-seq data. We demonstrate using publicly available scRNA-seq datasets and simulated expression data that bayNorm allows robust imputation of missing values generating realistic transcript distributions that match single molecule fluorescence in situ hybridization measurements. Moreover, by using priors informed by dataset structures, bayNorm improves accuracy and sensitivity of differential expression analysis and reduces batch effect compared with other existing methods. Altogether, bayNorm provides an efficient, integrated solution for global scaling normalization, imputation and true count recovery of gene expression measurements from scRNA-seq data.

Availability and implementation

The R package ‘bayNorm’ is publishd on bioconductor at https://bioconductor.org/packages/release/bioc/html/bayNorm.html. The code for analyzing data in this article is available at https://github.com/WT215/bayNorm_papercode.

Supplementary information

Supplementary data are available at Bioinformatics online.

Improved haplotype inference by exploiting long-range linking and allelic imbalance in RNA-seq datasets

Haplotype reconstruction of distant genetic variants remains an unsolved problem due to the short-read length of common sequencing data. Here, we introduce HapTree-X, a probabilistic framework that utilizes latent long-range information to reconstruct unspecified haplotypes in diploid and polyploid organisms. It introduces the observation that differential allele-specific expression can link genetic variants from the same physical chromosome, thus even enabling using reads that cover only individual variants. We demonstrate HapTree-X's feasibility on in-house sequenced Genome in a Bottle RNA-seq and various whole exome, genome, and 10X Genomics datasets. HapTree-X produces more complete phases (up to 25%), even in clinically important genes, and phases more variants than other methods while maintaining similar or higher accuracy and being up to 10×  faster than other tools. The advantage of HapTree-X's ability to use multiple lines of evidence, as well as to phase polyploid genomes in a single integrative framework, substantially grows as the amount of diverse data increases.

Computational Methods for Single-Cell RNA Sequencing

Single-cell RNA sequencing (scRNA-seq) has provided a high-dimensional catalog of millions of cells across species and diseases. These data have spurred the development of hundreds of computational tools to derive novel biological insights. Here, we outline the components of scRNA-seq analytical pipelines and the computational methods that underlie these steps. We describe available methods, highlight well-executed benchmarking studies, and identify opportunities for additional benchmarking studies and computational methods. As the biochemical approaches for single-cell omics advance, we propose coupled development of robust analytical pipelines suited for the challenges that new data present and principled selection of analytical methods that are suited for the biological questions to be addressed.

The barcode, UMI, set format and BUStools

SUMMARY

We introduce the Barcode-UMI-Set format (BUS) for representing pseudoalignments of reads from single-cell RNA-seq experiments. The format can be used with all single-cell RNA-seq technologies, and we show that BUS files can be efficiently generated. BUStools is a suite of tools for working with BUS files and facilitates rapid quantification and analysis of single-cell RNA-seq data. The BUS format therefore makes possible the development of modular, technology-specific and robust workflows for single-cell RNA-seq analysis.

AVAILABILITY AND IMPLEMENTATION

http://BUStools.github.io/ and http://pachterlab.github.io/kallisto/singlecell.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Solo: Doublet Identification in Single-Cell RNA-Seq via Semi-Supervised Deep Learning

Single-cell RNA sequencing (scRNA-seq) measurements of gene expression enable an unprecedented high-resolution view into cellular state. However, current methods often result in two or more cells that share the same cell-identifying barcode; these "doublets" violate the fundamental premise of single-cell technology and can lead to incorrect inferences. Here, we describe Solo, a semi-supervised deep learning approach that identifies doublets with greater accuracy than existing methods. Solo embeds cells unsupervised using a variational autoencoder and then appends a feed-forward neural network layer to the encoder to form a supervised classifier. We train this classifier to distinguish simulated doublets from the observed data. Solo can be applied in combination with experimental doublet detection methods to further purify scRNA-seq data to true single cells. It is freely available from https://github.com/calico/solo. A record of this paper's transparent peer review process is included in the Supplemental Information.

A scalable SCENIC workflow for single-cell gene regulatory network analysis

This protocol explains how to perform a fast SCENIC analysis alongside standard best practices steps on single-cell RNA-sequencing data using software containers and Nextflow pipelines. SCENIC reconstructs regulons (i.e., transcription factors and their target genes) assesses the activity of these discovered regulons in individual cells and uses these cellular activity patterns to find meaningful clusters of cells. Here we present an improved version of SCENIC with several advances. SCENIC has been refactored and reimplemented in Python (pySCENIC), resulting in a tenfold increase in speed, and has been packaged into containers for ease of use. It is now also possible to use epigenomic track databases, as well as motifs, to refine regulons. In this protocol, we explain the different steps of SCENIC: the workflow starts from the count matrix depicting the gene abundances for all cells and consists of three stages. First, coexpression modules are inferred using a regression per-target approach (GRNBoost2). Next, the indirect targets are pruned from these modules using cis-regulatory motif discovery (cisTarget). Lastly, the activity of these regulons is quantified via an enrichment score for the regulon's target genes (AUCell). Nonlinear projection methods can be used to display visual groupings of cells based on the cellular activity patterns of these regulons. The results can be exported as a loom file and visualized in the SCope web application. This protocol is illustrated on two use cases: a peripheral blood mononuclear cell data set and a panel of single-cell RNA-sequencing cancer experiments. For a data set of 10,000 genes and 50,000 cells, the pipeline runs in <2 h.

Advances of single-cell genomics and epigenomics in human disease: where are we now?

Cellular heterogeneity is revolutionizing the way to study, monitor and dissect complex diseases. This has been possible with the technological and computational advances associated to single-cell genomics and epigenomics. Deeper understanding of cell-to-cell variation and its impact on tissue function will open new avenues for early disease detection, accurate diagnosis and personalized treatments, all together leading to the next generation of health care. This review focuses on the recent discoveries that single-cell genomics and epigenomics have facilitated in the context of human health. It highlights the potential of single-cell omics to further advance the development of personalized treatments and precision medicine in cancer, diabetes and chronic age-related diseases. The promise of single-cell technologies to generate new insights about the differences in function between individual cells is just emerging, and it is paving the way for identifying biomarkers and novel therapeutic targets to tackle age, complex diseases and understand the effect of life style interventions and environmental factors.

Mapping human cell phenotypes to genotypes with single-cell genomics

The cumulative activity of all of the body's cells, with their myriad interactions, life histories, and environmental experiences, gives rise to a condition that is distinctly human and specific to each individual. It is an enduring goal to catalog our human cell types, to understand how they develop, how they vary between individuals, and how they fail in disease. Single-cell genomics has revolutionized this endeavor because sequencing-based methods provide a means to quantitatively annotate cell states on the basis of high-information content and high-throughput measurements. Together with advances in stem cell biology and gene editing, we are in the midst of a fascinating journey to understand the cellular phenotypes that compose human bodies and how the human genome is used to build and maintain each cell. Here, we will review recent advances into how single-cell genomics is being used to develop personalized phenotyping strategies that cross subcellular, cellular, and tissue scales to link our genome to our cumulative cellular phenotypes.

High-dimensional analysis of intestinal immune cells during helminth infection

Single cell isolation from helminth-infected murine intestines has been notoriously difficult, due to the strong anti-parasite type 2 immune responses that drive mucus production, tissue remodeling and immune cell infiltration. Through the systematic optimization of a standard intestinal digestion protocol, we were able to successfully isolate millions of immune cells from the heavily infected duodenum. To validate that these cells gave an accurate representation of intestinal immune responses, we analyzed them using a high-dimensional spectral flow cytometry panel and confirmed our findings by confocal microscopy. Our cell isolation protocol and high-dimensional analysis allowed us to identify many known hallmarks of anti-parasite immune responses throughout the entire course of helminth infection and has the potential to accelerate single-cell discoveries of local helminth immune responses that have previously been unfeasible.