Single-cell transcriptome sequencing: recent advances and remaining challenges

Single-cell transcriptome landscape of the kidney reveals potential innate immune regulation mechanisms in hybrid yellow catfish after Aeromonas hydrophila infection

Aeromonas hydrophila, the pathogen that is the causative agent of motile Aeromonas septicemia (MAS) disease, commonly attacks freshwater fishes, including yellow catfish (Pelteobagrus fulvidraco). Although the kidney is one of the most important organs involved in immunity in fish, its role in disease progression has not been fully elucidated. Understanding the cellular composition and innate immune regulation mechanisms of the kidney of yellow catfish is important for the treatment of MAS. In this study, single-cell RNA sequencing (scRNA-seq) was performed on the kidney of hybrid yellow catfish (Pelteobagrus fulvidraco ♀ × Pelteobagrus vachelli ♂) after A. hydrophila infection. Nine types of kidney cells were identified using marker genes, and a transcription module of marker genes in the main immune cells of hybrid yellow catfish kidney tissue was constructed using in-situ hybridization. In addition, the single-cell transcriptome data showed that the differentially expressed genes of macrophages were primarily enriched in the Toll-like receptor and Nod-like receptor signaling pathways. The expression levels of genes involved in these pathways were upregulated in macrophages following A. hydrophila infection. Transmission electron microscopy and TUNEL analysis revealed the cellular characteristics of macrophages before and after A. hydrophila infection. These data provide empirical support for in-depth research on the role of the kidney in the innate immune response of hybrid yellow catfish.

Performance comparison of high throughput single-cell RNA-Seq platforms in complex tissues

Single-cell transcriptomics has emerged as the preferred tool to define cell identity through the analysis of gene expression signatures. However, there are limited studies that have comprehensively compared the performance of different scRNAseq systems in complex tissues. Here, we present a systematic comparison of two well-established high throughput 3'-scRNAseq platforms: 10× Chromium and BD Rhapsody, using tumours that present high cell diversity. Our experimental design includes both fresh and artificially damaged samples from the same tumours, which also provides a comparable dataset to examine their performance under challenging conditions. The performance metrics used in this study consist of gene sensitivity, mitochondrial content, reproducibility, clustering capabilities, cell type representation and ambient RNA contamination. These analyses showed that BD Rhapsody and 10× Chromium have similar gene sensitivity, while BD Rhapsody has the highest mitochondrial content. Interestingly, we found cell type detection biases between platforms, including a lower proportion of endothelial and myofibroblast cells in BD Rhapsody and lower gene sensitivity in granulocytes for 10× Chromium. Moreover, the source of the ambient noise was different between plate-based and droplet-based platforms. In conclusion, our reported platform differential performance should be considered for the selection of the scRNAseq method during the study experimental designs.

Identifying prognostic biomarkers in oral squamous cell carcinoma: an integrated single-cell and bulk RNA sequencing study on mitophagy-related genes

Oral squamous cell carcinoma (OSCC) has an extremely poor prognosis. Recent studies have suggested that mitophagy-related genes (MRGs) are closely correlated with the development and occurrence of cancer, but the role they play in oral cancer has not yet been explained.We conducted a comprehensive analysis of integrated single-cell and bulk RNA sequencing (RNA-seq) data retrieved from Gene Expression Omnibus (GEO) datasets and The Cancer Genome Atlas (TCGA) database. Multiple methods were combined to provide a comprehensive understanding of the genetic expression patterns and biology of OSCC, such as analysis of pseudotime series, CellChat cell communication, immune infiltration, Gene Ontology (GO), LASSO Cox regression, gene set variation analysis (GSVA), Kyoto Encyclopedia of Genes and Genomes (KEGG), gene set enrichment analysis (GSEA), Tumor Mutation Burden (TMB) and drug sensitivity assessments. The findings of this study demonstrated significantly greater activity of MRGs in NK cells than in other cells in OSCC. A reliable prognostic model was developed using 12 candidate genes strongly associated with mitochondrial autophagy. T stage, N stage and risk score were revealed as independent prognostic factors. Distinctively enriched pathways and immune cells were observed in different risk groups. Notably, low-risk patients were more responsive to chemotherapy. In addition, a nomogram model with excellent predictive ability was established by combining the risk scores and clinical features. The activity of MRGs suggest the potential for the development of new targeted therapies. The construction of a robust prognostic model also provides reference value for individualized prediction and clinical decision-making in patients with OSCC.

Decoding cellular plasticity and niche regulation of limbal stem cells during corneal wound healing

BACKGROUND

Dysfunction or deficiency of corneal epithelium results in vision impairment or blindness in severe cases. The rapid and effective regeneration of corneal epithelial cells relies on the limbal stem cells (LSCs). However, the molecular and functional responses of LSCs and their niche cells to injury remain elusive.

METHODS

Single-cell RNA sequencing was performed on corneal tissues from normal mice and corneal epithelium defect models. Bioinformatics analysis was performed to confirm the distinct characteristics and cell fates of LSCs. Knockdown of Creb5 and OSM treatment experiment were performed to determine their roles of in corneal epithelial wound healing.

RESULTS

Our data defined the molecular signatures of LSCs and reconstructed the pseudotime trajectory of corneal epithelial cells. Gene network analyses characterized transcriptional landmarks that potentially regulate LSC dynamics, and identified a transcription factor Creb5, that was expressed in LSCs and significantly upregulated after injury. Loss-of-function experiments revealed that silencing Creb5 delayed the corneal epithelial healing and LSC mobilization. Through cell-cell communication analysis, we identified 609 candidate regeneration-associated ligand-receptor interaction pairs between LSCs and distinct niche cells, and discovered a unique subset of Arg1+ macrophages infiltrated after injury, which were present as the source of Oncostatin M (OSM), an IL-6 family cytokine, that were demonstrated to effectively accelerate the corneal epithelial wound healing.

CONCLUSIONS

This research provides a valuable single-cell resource and reference for the discovery of mechanisms and potential clinical interventions aimed at ocular surface reconstruction.

Twin Scholarships of Glycomedicine and Precision Medicine in Times of Single-Cell Multiomics

Systems biology and multiomics research expand the prospects of planetary health innovations. In this context, this mini-review unpacks the twin scholarships of glycomedicine and precision medicine in the current era of single-cell multiomics. A significant growth in glycan research has been observed over the past decade, unveiling and establishing co- and post-translational modifications as dynamic indicators of both pathological and physiological conditions. Systems biology technologies have enabled large-scale and high-throughput glycoprofiling and access to data-intensive biological repositories for global research. These advancements have established glycans as a pivotal third code of life, alongside nucleic acids and amino acids. However, challenges persist, particularly in the simultaneous analysis of the glycome and transcriptome in single cells owing to technical limitations. In addition, holistic views of the complex molecular interactions between glycomics and other omics types remain elusive. We underscore and call for a paradigm shift toward the exploration of integrative glycan platforms and analysis methods for single-cell multiomics research and precision medicine biomarker discovery. The integration of multiple datasets from various single-cell omics levels represents a crucial application of systems biology in understanding complex cellular processes and is essential for advancing the twin scholarships of glycomedicine and precision medicine.

Integrating phylogenies into single-cell RNA sequencing analysis allows comparisons across species, genes, and cells

Comparisons of single-cell RNA sequencing (scRNA-seq) data across species can reveal links between cellular gene expression and the evolution of cell functions, features, and phenotypes. These comparisons evoke evolutionary histories, as depicted by phylogenetic trees, that define relationships between species, genes, and cells. This Essay considers each of these in turn, laying out challenges and solutions derived from a phylogenetic comparative approach and relating these solutions to previously proposed methods for the pairwise alignment of cellular dimensional maps. This Essay contends that species trees, gene trees, cell phylogenies, and cell lineages can all be reconciled as descriptions of the same concept-the tree of cellular life. By integrating phylogenetic approaches into scRNA-seq analyses, challenges for building informed comparisons across species can be overcome, and hypotheses about gene and cell evolution can be robustly tested.

Techniques for investigating lncRNA transcript functions in neurodevelopment

Long noncoding RNAs (lncRNAs) are sequences of 200 nucleotides or more that are transcribed from a large portion of the mammalian genome. While hypothesized to have a variety of biological roles, many lncRNAs remain largely functionally uncharacterized due to unique challenges associated with their investigation. For example, some lncRNAs overlap with other genomic loci, are expressed in a cell-type-specific manner, and/or are differentially processed at the post-transcriptional level. The mammalian CNS contains a vast diversity of lncRNAs, and lncRNAs are highly abundant in the mammalian brain. However, interrogating lncRNA function in models of the CNS, particularly in vivo, can be complex and challenging. Here we review the breadth of methods used to investigate lncRNAs in the CNS, their merits, and the understanding they can provide with respect to neurodevelopment and pathophysiology. We discuss remaining challenges in the field and provide recommendations to assay lncRNAs based on current methods.

An ultra high-throughput, massively multiplexable, single-cell RNA-seq platform in yeasts

Yeasts are naturally diverse, genetically tractable, and easy to grow such that researchers can investigate any number of genotypes, environments, or interactions thereof. However, studies of yeast transcriptomes have been limited by the processing capabilities of traditional RNA sequencing techniques. Here we optimize a powerful, high-throughput single-cell RNA sequencing (scRNAseq) platform, SPLiT-seq (Split Pool Ligation-based Transcriptome sequencing), for yeasts and apply it to 43,388 cells of multiple species and ploidies. This platform utilizes a combinatorial barcoding strategy to enable massively parallel RNA sequencing of hundreds of yeast genotypes or growth conditions at once. This method can be applied to most species or strains of yeast for a fraction of the cost of traditional scRNAseq approaches. Thus, our technology permits researchers to leverage "the awesome power of yeast" by allowing us to survey the transcriptome of hundreds of strains and environments in a short period of time and with no specialized equipment. The key to this method is that sequential barcodes are probabilistically appended to cDNA copies of RNA while the molecules remain trapped inside of each cell. Thus, the transcriptome of each cell is labeled with a unique combination of barcodes. Since SPLiT-seq uses the cell membrane as a container for this reaction, many cells can be processed together without the need to physically isolate them from one another in separate wells or droplets. Further, the first barcode in the sequence can be chosen intentionally to identify samples from different environments or genetic backgrounds, enabling multiplexing of hundreds of unique perturbations in a single experiment. In addition to greater multiplexing capabilities, our method also facilitates a deeper investigation of biological heterogeneity, given its single-cell nature. For example, in the data presented here, we detect transcriptionally distinct cell states related to cell cycle, ploidy, metabolic strategies, and so forth, all within clonal yeast populations grown in the same environment. Hence, our technology has two obvious and impactful applications for yeast research: the first is the general study of transcriptional phenotypes across many strains and environments, and the second is investigating cell-to-cell heterogeneity across the entire transcriptome.

Single-cell transcriptomic analysis reveals transcriptional and cell subpopulation differences between human and pig immune cells

BACKGROUND

The pig is a promising donor candidate for xenotransplantation. Understanding the differences between human and swine immune systems is critical for addressing xenotransplant rejection and hematopoietic reconstitution. The gene transcriptional profile differences between human and pig immune cell subpopulations have not been studied. To assess the similarities and differences between pigs and humans at the levels of gene transcriptional profiles or cell subpopulations are important for better understanding the cross-species similarity of humans and pigs, and it would help establish the fundamental principles necessary to genetically engineer donor pigs and improve xenotransplantation.

OBJECTIVE

To assess the gene transcriptional similarities and differences between pigs and humans.

METHODS

Two pigs and two healthy humans' PBMCs were sorted for 10 × genomics single-cell sequence. We generated integrated human-pig scRNA-seq data from human and pig PBMCs and defined the overall gene expression landscape of pig peripheral blood immune cell subpopulations by updating the set of human-porcine homologous genes. The subsets of immune cells were detected by flow cytometry.

RESULTS

There were significantly less T cells, NK cells and monocytes but more B cells in pig peripheral blood than those in human peripheral blood. High oxidative phosphorylation, HIF-1, glycolysis, and lysosome-related gene expressions in pig CD14+ monocytes were observed, whereas pig CD14+ monocytes exhibited lower levels of cytokine receptors and JAK-STAT-related genes. Pig activated CD4+T cells decreased cell adhesion and inflammation, while enriched for migration and activation processes. Porcine GNLY+CD8+T cells reduced cytotoxicity and increased proliferation compared with human GNLY+CD8+T cells. Pig CD2+CD8+γδT cells were functionally homologous to human CD2+CD4+ γδT cells. Pig CD2-CD8-γδT cells expressed genes with quiescent and precursor characteristics, while CD2-CD8+γδT cells expressed migration and memory-related molecules. Pig CD24+ and CD5+B cells are associated with inflammatory responses.

CONCLUSION

Our research with integrated scRNA-seq assays identified the different distribution of pig immune cell subpopulations and the different transcriptional profiles of human and pig immune cells. This study enables a deeper understanding of the development and function of porcine immune cells.

A Data-Distribution and Successive Spline Points based discretization approach for evolving gene regulatory networks from scRNA-Seq time-series data using Cartesian Genetic Programming

The inference of gene regulatory networks (GRNs) is a widely addressed problem in Systems Biology. GRNs can be modeled as Boolean networks, which is the simplest approach for this task. However, Boolean models need binarized data. Several approaches have been developed for the discretization of gene expression data (GED). Also, the advance of data extraction technologies, such as single-cell RNA-Sequencing (scRNA-Seq), provides a new vision of gene expression and brings new challenges for dealing with its specificities, such as a large occurrence of zero data. This work proposes a new discretization approach for dealing with scRNA-Seq time-series data, named Distribution and Successive Spline Points Discretization (DSSPD), which considers the data distribution and a proper preprocessing step. Here, Cartesian Genetic Programming (CGP) is used to infer GRNs using the results of DSSPD. The proposal is compared with CGP with the standard data handling and five state-of-the-art algorithms on curated models and experimental data. The results show that the proposal improves the results of CGP in all tested cases and outperforms the state-of-the-art algorithms in most cases.

Normalizing need not be the norm: count-based math for analyzing single-cell data

Counting transcripts of mRNA are a key method of observation in modern biology. With advances in counting transcripts in single cells (single-cell RNA sequencing or scRNA-seq), these data are routinely used to identify cells by their transcriptional profile, and to identify genes with differential cellular expression. Because the total number of transcripts counted per cell can vary for technical reasons, the first step of many commonly used scRNA-seq workflows is to normalize by sequencing depth, transforming counts into proportional abundances. The primary objective of this step is to reshape the data such that cells with similar biological proportions of transcripts end up with similar transformed measurements. But there is growing concern that normalization and other transformations result in unintended distortions that hinder both analyses and the interpretation of results. This has led to an intense focus on optimizing methods for normalization and transformation of scRNA-seq data. Here, we take an alternative approach, by avoiding normalization and transformation altogether. We abandon the use of distances to compare cells, and instead use a restricted algebra, motivated by measurement theory and abstract algebra, that preserves the count nature of the data. We demonstrate that this restricted algebra is sufficient to draw meaningful and practical comparisons of gene expression through the use of the dot product and other elementary operations. This approach sidesteps many of the problems with common transformations, and has the added benefit of being simpler and more intuitive. We implement our approach in the package countland, available in python and R.

The Use of Non-Apoptotic Sperm Selected by Magnetic Activated Cell Sorting (MACS) to Enhance Reproductive Outcomes: What the Evidence Says

Sperm selection of the most competent sperm is a promising way to enhance reproductive outcomes. Apoptosis is the programmed cell death process to maintain tissue homeostasis, and MACS sperm selection of non-apoptotic cells enables the removal of apoptotic sperm from an ejaculate, thus leaving the non-apoptotic available to be microinjected, but given the associated costs of adding these sperm selection steps to the routine practice, there is a need for a careful examination of the literature available to answer questions such as who can benefit from this MACS, how significant this improvement is, and how robust the evidence and data available supporting this choice are. Thus, the aim of this narrative review was to objectively evaluate the available evidence regarding the potential benefits of the use of MACS. From the literature, there are controversial results since its implementation as an in vitro fertilization add-on, and this may be explained in part by the low quality of the evidence available, wrong designs, or even inadequate statistical analyses. We concluded that the benefits of adding MACS are unclear, and further methodologically sound research on specific populations is much needed before offering it clinically.

Fibroblasts - the cellular choreographers of wound healing

Injuries to our skin trigger a cascade of spatially- and temporally-synchronized healing processes. During such endogenous wound repair, the role of fibroblasts is multifaceted, ranging from the activation and recruitment of innate immune cells through the synthesis and deposition of scar tissue to the conveyor belt-like transport of fascial connective tissue into wounds. A comprehensive understanding of fibroblast diversity and versatility in the healing machinery may help to decipher wound pathologies whilst laying the foundation for novel treatment modalities. In this review, we portray the diversity of fibroblasts and delineate their unique wound healing functions. In addition, we discuss future directions through a clinical-translational lens.

Roles of Non-Coding RNA in Alzheimer's Disease Pathophysiology

Alzheimer's disease (AD) is a chronic neurodegenerative disorder that is accompanied by deficits in memory and cognitive functions. The disease is pathologically characterised by the accumulation and aggregation of an extracellular peptide referred to as amyloid-β (Aβ) in the form of amyloid plaques and the intracellular aggregation of a hyperphosphorelated protein tau in the form of neurofibrillary tangles (NFTs) that cause neuroinflammation, synaptic dysfunction, and oxidative stress. The search for pathomechanisms leading to disease onset and progression has identified many key players that include genetic, epigenetic, behavioural, and environmental factors, which lend support to the fact that this is a multi-faceted disease where failure in various systems contributes to disease onset and progression. Although the vast majority of individuals present with the sporadic (non-genetic) form of the disease, dysfunctions in numerous protein-coding and non-coding genes have been implicated in mechanisms contributing to the disease. Recent studies have provided strong evidence for the association of non-coding RNAs (ncRNAs) with AD. In this review, we highlight the current findings on changes observed in circular RNA (circRNA), microRNA (miRNA), short interfering RNA (siRNA), piwi-interacting RNA (piRNA), and long non-coding RNA (lncRNA) in AD. Variations in these ncRNAs could potentially serve as biomarkers or therapeutic targets for the diagnosis and treatment of Alzheimer's disease. We also discuss the results of studies that have targeted these ncRNAs in cellular and animal models of AD with a view for translating these findings into therapies for Alzheimer's disease.

Integrated analysis of plasma proteome and cortex single-cell transcriptome reveals the novel biomarkers during cortical aging

Background

With the increase of age, multiple physiological functions of people begin gradually degenerating. Regardless of natural aging or pathological aging, the decline in cognitive function is one of the most obvious features in the process of brain aging. Brain aging is a key factor for several neuropsychiatric disorders and for most neurodegenerative diseases characterized by onset typically occurring late in life and with worsening of symptoms over time. Therefore, the early prevention and intervention of aging progression are particularly important. Since there is no unified conclusion about the plasma diagnostic biomarkers of brain aging, this paper innovatively employed the combined multi-omics analysis to delineate the plasma markers of brain aging.

Methods

In order to search for specific aging markers in plasma during cerebral cortex aging, we used multi-omics analysis to screen out differential genes/proteins by integrating two prefrontal cortex (PFC) single-nucleus transcriptome sequencing (snRNA-seq) datasets and one plasma proteome sequencing datasets. Then plasma samples were collected from 20 young people and 20 elder people to verify the selected differential genes/proteins with ELISA assay.

Results

We first integrated snRNA-seq data of the post-mortem human PFC and generated profiles of 65,064 nuclei from 14 subjects across adult (44-58 years), early-aging (69-79 years), and late-aging (85-94 years) stages. Seven major cell types were classified based on established markers, including oligodendrocyte, excitatory neurons, oligodendrocyte progenitor cells, astrocytes, microglia, inhibitory neurons, and endotheliocytes. A total of 93 cell-specific genes were identified to be significantly associated with age. Afterward, plasma proteomics data from 2,925 plasma proteins across 4,263 young adults to nonagenarians (18-95 years old) were combined with the outcomes from snRNA-seq data to obtain 12 differential genes/proteins (GPC5, CA10, DGKB, ST6GALNAC5, DSCAM, IL1RAPL2, TMEM132C, VCAN, APOE, PYH1R, CNTN2, SPOCK3). Finally, we verified the 12 differential genes by ELISA and found that the expression trends of five biomarkers (DSCAM, CNTN2, IL1RAPL2, CA10, GPC5) were correlated with brain aging.

Conclusion

Five differentially expressed proteins (DSCAM, CNTN2, IL1RAPL2, CA10, GPC5) can be considered as one of the screening indicators of brain aging, and provide a scientific basis for clinical diagnosis and intervention.

Single-cell RNA sequencing in cornea research: Insights into limbal stem cells and their niche regulation

The corneal epithelium is composed of stratified squamous epithelial cells on the outer surface of the eye, which acts as a protective barrier and is critical for clear and stable vision. Its continuous renewal or wound healing depends on the proliferation and differentiation of limbal stem cells (LSCs), a cell population that resides at the limbus in a highly regulated niche. Dysfunction of LSCs or their niche can cause limbal stem cell deficiency, a disease that is manifested by failed epithelial wound healing or even blindness. Nevertheless, compared to stem cells in other tissues, little is known about the LSCs and their niche. With the advent of single-cell RNA sequencing, our understanding of LSC characteristics and their microenvironment has grown considerably. In this review, we summarized the current findings from single-cell studies in the field of cornea research and focused on important advancements driven by this technology, including the heterogeneity of the LSC population, novel LSC markers and regulation of the LSC niche, which will provide a reference for clinical issues such as corneal epithelial wound healing, ocular surface reconstruction and interventions for related diseases.

Single-Cell Analysis in Immuno-Oncology

The complexity of the cellular and non-cellular milieu surrounding human tumors plays a decisive role in the course and outcome of disease. The high variability in the distribution of the immune and non-immune compartments within the tumor microenvironments (TME) among different patients governs the mode of their response or resistance to current immunotherapeutic approaches. Through deciphering this diversity, one can tailor patients' management to meet an individual's needs. Single-cell (sc) omics technologies have given a great boost towards this direction. This review gathers recent data about how multi-omics profiling, including the utilization of single-cell RNA sequencing (scRNA-seq), assay for transposase-accessible chromatin with sequencing (scATAC-seq), T-cell receptor sequencing (scTCR-seq), mass, tissue-based, or microfluidics cytometry, and related bioinformatics tools, contributes to the high-throughput assessment of a large number of analytes at single-cell resolution. Unravelling the exact TCR clonotype of the infiltrating T cells or pinpointing the classical or novel immune checkpoints across various cell subsets of the TME provide a boost to our comprehension of adaptive immune responses, their antigen specificity and dynamics, and grant suggestions for possible therapeutic targets. Future steps are expected to merge high-dimensional data with tissue localization data, which can serve the investigation of novel multi-modal biomarkers for the selection and/or monitoring of the optimal treatment from the current anti-cancer immunotherapeutic armamentarium.

Genetic and Molecular Contributors in Eosinophilic Esophagitis

PURPOSE OF REVIEW

Eosinophilic esophagitis (EoE) is an allergic inflammatory esophageal disorder with a complex underlying genetic and molecular etiology. The interest of the scientific community in EoE has grown considerably over the past three decades, and the understanding of the genetic and molecular mechanisms involved in this disease has greatly increased.

RECENT FINDINGS

In this article, we aim to provide both historic aspects and updates on the recent genetic and molecular advances in the understanding of EoE. Although EoE is a relatively newly described disorder, much progress has been made toward identifying the genetic and molecular factors contributing to the disease pathogenesis by a variety of approaches with next-generation sequencing technologies, including genome-wide association study, whole exome sequencing, and bulk and single-cell RNA sequencing. This review highlights the multifaceted impacts of various findings that have shaped the current molecular and genetic landscape of EoE, providing insights that facilitate further understanding of the disease process.

SMURF: embedding single-cell RNA-seq data with matrix factorization preserving self-consistency

The advance in single-cell RNA-sequencing (scRNA-seq) sheds light on cell-specific transcriptomic studies of cell developments, complex diseases and cancers. Nevertheless, scRNA-seq techniques suffer from 'dropout' events, and imputation tools are proposed to address the sparsity. Here, rather than imputation, we propose a tool, SMURF, to extract the low-dimensional embeddings from cells and genes utilizing matrix factorization with a mixture of Poisson-Gamma divergent as objective while preserving self-consistency. SMURF exhibits feasible cell subpopulation discovery efficacy with obtained cell embeddings on replicated in silico and eight web lab scRNA datasets with ground truth cell types. Furthermore, SMURF can reduce the cell embedding to a 1D-oval space to recover the time course of cell cycle. SMURF can also serve as an imputation tool; the in silico data assessment shows that SMURF parades the most robust gene expression recovery power with low root mean square error and high Pearson correlation. Moreover, SMURF recovers the gene distribution for the WM989 Drop-seq data. SMURF is available at https://github.com/deepomicslab/SMURF.

BLTSA: pseudotime prediction for single cells by branched local tangent space alignment

MOTIVATION

The development of single-cell RNA sequencing (scRNA-seq) technology makes it possible to study the cellular dynamic processes such as cell cycle and cell differentiation. Due to the difficulties in generating genuine time-series scRNA-seq data, it is of great importance to computationally infer the pseudotime of the cells along differentiation trajectory based on their gene expression patterns. The existing pseudotime prediction methods often suffer from the high level noise of single-cell data, thus it is still necessary to study the single-cell trajectory inference methods.

RESULTS

In this study, we propose a branched local tangent space alignment (BLTSA) method to infer single-cell pseudotime for multi-furcation trajectories. By assuming that single cells are sampled from a low-dimensional self-intersecting manifold, BLTSA first identifies the tip and branching cells in the trajectory based on cells' local Euclidean neighborhoods. Local coordinates within the tangent spaces are then determined by each cell's local neighborhood after clustering all the cells to different branches iteratively. The global coordinates for all the single cells are finally obtained by aligning the local coordinates based on the tangent spaces. We evaluate the performance of BLTSA on four simulation datasets and five real datasets. The experimental results show that BLTSA has obvious advantages over other comparison methods.

AVAILABILITY AND IMPLEMENTATION

R codes are available at https://github.com/LiminLi-xjtu/BLTSA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

scRNA-seq data analysis method to improve analysis performance

With the development of single-cell RNA sequencing technology (scRNA-seq), we have the ability to study biological questions at the level of the individual cell transcriptome. Nowadays, many analysis tools, specifically suitable for single-cell RNA sequencing data, have been developed. In this review, the currently commonly used scRNA-seq protocols are discussed. The upstream processing flow pipeline of scRNA-seq data, including goals and popular tools for reads mapping and expression quantification, quality control, normalization, imputation, and batch effect removal is also introduced. Finally, methods to evaluate these tools in both cellular and genetic dimensions, clustering and differential expression analysis are presented.

Absolute scaling of single-cell transcriptomes identifies pervasive hypertranscription in adult stem and progenitor cells

Hypertranscription supports biosynthetically demanding cellular states through global transcriptome upregulation. Despite its potential widespread relevance, documented examples of hypertranscription remain few and limited to early development. Here, we demonstrate that absolute scaling of single-cell RNA-sequencing data enables the estimation of total transcript abundances per cell. We validate absolute scaling in known cases of developmental hypertranscription and apply it to adult cell types, revealing a remarkable dynamic range in transcriptional output. In adult organs, hypertranscription marks activated stem/progenitor cells with multilineage potential and is redeployed in conditions of tissue injury, where it precedes bursts of proliferation during regeneration. Our analyses identify a common set of molecular pathways associated with both adult and embryonic hypertranscription, including chromatin remodeling, DNA repair, ribosome biogenesis, and translation. These shared features across diverse cell contexts support hypertranscription as a general and dynamic cellular program that is pervasively employed during development, organ maintenance, and regeneration.

Improved downstream functional analysis of single-cell RNA-sequence data using DGAN

The dramatic increase in the number of single-cell RNA-sequence (scRNA-seq) investigations is indeed an endorsement of the new-fangled proficiencies of next generation sequencing technologies that facilitate the accurate measurement of tens of thousands of RNA expression levels at the cellular resolution. Nevertheless, missing values of RNA amplification persist and remain as a significant computational challenge, as these data omission induce further noise in their respective cellular data and ultimately impede downstream functional analysis of scRNA-seq data. Consequently, it turns imperative to develop robust and efficient scRNA-seq data imputation methods for improved downstream functional analysis outcomes. To overcome this adversity, we have designed an imputation framework namely deep generative autoencoder network [DGAN]. In essence, DGAN is an evolved variational autoencoder designed to robustly impute data dropouts in scRNA-seq data manifested as a sparse gene expression matrix. DGAN principally reckons count distribution, besides data sparsity utilizing a gaussian model whereby, cell dependencies are capitalized to detect and exclude outlier cells via imputation. When tested on five publicly available scRNA-seq data, DGAN outperformed every single baseline method paralleled, with respect to downstream functional analysis including cell data visualization, clustering, classification and differential expression analysis. DGAN is executed in Python and is accessible at https://github.com/dikshap11/DGAN .

Single-cell transcriptome reveals core cell populations and androgen-RXFP2 axis involved in deer antler full regeneration

Deer antlers constitute a unique mammalian model for the study of both organ formation in postnatal life and annual full regeneration. Previous studies revealed that these events are achieved through the proliferation and differentiation of antlerogenic periosteum (AP) cells and pedicle periosteum (PP) cells, respectively. As the cells resident in the AP and the PP possess stem cell attributes, both antler generation and regeneration are stem cell-based processes. However, the cell composition of each tissue type and molecular events underlying antler development remain poorly characterized. Here, we took the approach of single-cell RNA sequencing (scRNA-Seq) and identified eight cell types (mainly THY1⁺ cells, progenitor cells, and osteochondroblasts) and three core subclusters of the THY1⁺ cells (SC2, SC3, and SC4). Endothelial and mural cells each are heterogeneous at transcriptional level. It was the proliferation of progenitor, mural, and endothelial cells in the activated antler-lineage-specific tissues that drove the rapid formation of the antler. We detected the differences in the initial differentiation process between antler generation and regeneration using pseudotime trajectory analysis. These may be due to the difference in the degree of stemness of the AP-THY1⁺ and PP-THY1⁺ cells. We further found that androgen-RXFP2 axis may be involved in triggering initial antler full regeneration. Fully deciphering the cell composition for these antler tissue types will open up new avenues for elucidating the mechanism underlying antler full renewal in specific and regenerative medicine in general.

Integrated analyses reveal evolutionarily conserved and specific injury response genes in dorsal root ganglion

Rodent dorsal root ganglion (DRG) is widely used for studying axonal injury. Extensive studies have explored genome-wide profiles on rodent DRGs under peripheral nerve insults. However, systematic integration and exploration of these data still be limited. Herein, we re-analyzed 21 RNA-seq datasets and presented a web-based resource (DRGProfile). We identified 53 evolutionarily conserved injury response genes, including well-known injury genes (Atf3, Npy and Gal) and less-studied transcriptional factors (Arid5a, Csrnp1, Zfp367). Notably, we identified species-preference injury response candidates (e.g. Gpr151, Lipn, Anxa10 in mice; Crisp3, Csrp3, Vip, Hamp in rats). Temporal profile analysis reveals expression patterns of genes related to pre-regenerative and regenerating states. Finally, we found a large sex difference in response to sciatic nerve injury, and identified four male-specific markers (Uty, Eif2s3y, Kdm5d, Ddx3y) expressed in DRG. Our study provides a comprehensive integrated landscape for expression change in DRG upon injury which will greatly contribute to the neuroscience community.

Comprehensive multi-omics single-cell data integration reveals greater heterogeneity in the human immune system

Single-cell transcriptomics enables the definition of diverse human immune cell types across multiple tissues and disease contexts. Further deeper biological understanding requires comprehensive integration of multiple single-cell omics (transcriptomic, proteomic, and cell-receptor repertoire). To improve the identification of diverse cell types and the accuracy of cell-type classification in multi-omics single-cell datasets, we developed SuPERR, a novel analysis workflow to increase the resolution and accuracy of clustering and allow for the discovery of previously hidden cell subsets. In addition, SuPERR accurately removes cell doublets and prevents widespread cell-type misclassification by incorporating information from cell-surface proteins and immunoglobulin transcript counts. This approach uniquely improves the identification of heterogeneous cell types and states in the human immune system, including rare subsets of antibody-secreting cells in the bone marrow.

•

SuPERR removes heterotypic doublets and cell-type misclassifications in scRNA-seq

•

Sequential gating on cell-surface proteins resolves major cell lineages in scRNA-seq

•

Defining major cell lineages before clustering reduces cell-type misclassifications

•

Antibody counts from single-cell V(D)J matrix accurately identify plasma cells

Executable models of immune signaling pathways in HIV-associated atherosclerosis

Atherosclerosis (AS)-associated cardiovascular disease is an important cause of mortality in an aging population of people living with HIV (PLWH). This elevated risk has been attributed to viral infection, anti-retroviral therapy, chronic inflammation, and lifestyle factors. However, the rates at which PLWH develop AS vary even after controlling for length of infection, treatment duration, and for lifestyle factors. To investigate the molecular signaling underlying this variation, we sequenced 9368 peripheral blood mononuclear cells (PBMCs) from eight PLWH, four of whom have atherosclerosis (AS+). Additionally, a publicly available dataset of PBMCs from persons before and after HIV infection was used to investigate the effect of acute HIV infection. To characterize dysregulation of pathways rather than just measuring enrichment, we developed the single-cell Boolean Omics Network Invariant Time Analysis (scBONITA) algorithm. scBONITA infers executable dynamic pathway models and performs a perturbation analysis to identify high impact genes. These dynamic models are used for pathway analysis and to map sequenced cells to characteristic signaling states (attractor analysis). scBONITA revealed that lipid signaling regulates cell migration into the vascular endothelium in AS+ PLWH. Pathways implicated included AGE-RAGE and PI3K-AKT signaling in CD8+ T cells, and glucagon and cAMP signaling pathways in monocytes. Attractor analysis with scBONITA facilitated the pathway-based characterization of cellular states in CD8+ T cells and monocytes. In this manner, we identify critical cell-type specific molecular mechanisms underlying HIV-associated atherosclerosis using a novel computational method.

Wound healing, fibroblast heterogeneity, and fibrosis

Fibroblasts are highly dynamic cells that play a central role in tissue repair and fibrosis. However, the mechanisms by which they contribute to both physiologic and pathologic states of extracellular matrix deposition and remodeling are just starting to be understood. In this review article, we discuss the current state of knowledge in fibroblast biology and heterogeneity, with a primary focus on the role of fibroblasts in skin wound repair. We also consider emerging techniques in the field, which enable an increasingly nuanced and contextualized understanding of these complex systems, and evaluate limitations of existing methodologies and knowledge. Collectively, this review spotlights a diverse body of research examining an often-overlooked cell type-the fibroblast-and its critical functions in wound repair and beyond.

Differential Expression Analysis of Single-Cell RNA-Seq Data: Current Statistical Approaches and Outstanding Challenges

With the advent of single-cell RNA-sequencing (scRNA-seq), it is possible to measure the expression dynamics of genes at the single-cell level. Through scRNA-seq, a huge amount of expression data for several thousand(s) of genes over million(s) of cells are generated in a single experiment. Differential expression analysis is the primary downstream analysis of such data to identify gene markers for cell type detection and also provide inputs to other secondary analyses. Many statistical approaches for differential expression analysis have been reported in the literature. Therefore, we critically discuss the underlying statistical principles of the approaches and distinctly divide them into six major classes, i.e., generalized linear, generalized additive, Hurdle, mixture models, two-class parametric, and non-parametric approaches. We also succinctly discuss the limitations that are specific to each class of approaches, and how they are addressed by other subsequent classes of approach. A number of challenges are identified in this study that must be addressed to develop the next class of innovative approaches. Furthermore, we also emphasize the methodological challenges involved in differential expression analysis of scRNA-seq data that researchers must address to draw maximum benefit from this recent single-cell technology. This study will serve as a guide to genome researchers and experimental biologists to objectively select options for their analysis.

Identification of new marker genes from plant single-cell RNA-seq data using interpretable machine learning methods

An essential step in the analysis of single-cell RNA sequencing data is to classify cells into specific cell types using marker genes. In this study, we have developed a machine learning pipeline called single-cell predictive marker (SPmarker) to identify novel cell-type marker genes in the Arabidopsis root. Unlike traditional approaches, our method uses interpretable machine learning models to select marker genes. We have demonstrated that our method can: assign cell types based on cells that were labelled using published methods; project cell types identified by trajectory analysis from one data set to other data sets; and assign cell types based on internal GFP markers. Using SPmarker, we have identified hundreds of new marker genes that were not identified before. As compared to known marker genes, the new marker genes have more orthologous genes identifiable in the corresponding rice single-cell clusters. The new root hair marker genes also include 172 genes with orthologs expressed in root hair cells in five non-Arabidopsis species, which expands the number of marker genes for this cell type by 35-154%. Our results represent a new approach to identifying cell-type marker genes from scRNA-seq data and pave the way for cross-species mapping of scRNA-seq data in plants.

Novel insights into adipose tissue heterogeneity

When normalized to volume, adipose tissue is comprised mainly of large lipid metabolizing and storing cells called adipocytes. Strikingly, the numerical representation of non-adipocytes, composed of a wide variety of cell types found in the so-called stromal vascular fraction (SVF), outnumber adipocytes by far. Besides its function in energy storage, adipose tissue has emerged as a versatile organ that regulates systemic metabolism and has therefore constituted an attractive target for the treatment of metabolic diseases. Recent high-resolution single cells/nucleus RNA seq data exemplify an intriguingly profound diversity of both adipocytes and SVF cells in all adipose depots, and the current data, while limited, demonstrate the significance of the intra-tissue cell composition in shaping the overall functionality of this tissue. Due to the complexity of adipose tissue, our understanding of the biological relevance of this heterogeneity and plasticity is fractional. Therefore, establishing atlases of adipose tissue cell heterogeneity is the first step towards generating an understanding of these functionalities. In this review, we will describe the current knowledge on adipose tissue cell composition and the heterogeneity of single-cell RNA sequencing, including the technical limitations.

Single-cell RNA sequencing to track novel perspectives in HSC heterogeneity

As the importance of cell heterogeneity has begun to be emphasized, single-cell sequencing approaches are rapidly adopted to study cell heterogeneity and cellular evolutionary relationships of various cells, including stem cell populations. The hematopoietic stem and progenitor cell (HSPC) compartment contains HSC hematopoietic stem cells (HSCs) and distinct hematopoietic cells with different abilities to self-renew. These cells perform their own functions to maintain different hematopoietic lineages. Undeniably, single-cell sequencing approaches, including single-cell RNA sequencing (scRNA-seq) technologies, empower more opportunities to study the heterogeneity of normal and pathological HSCs. In this review, we discuss how these scRNA-seq technologies contribute to tracing origin and lineage commitment of HSCs, profiling the bone marrow microenvironment and providing high-resolution dissection of malignant hematopoiesis, leading to exciting new findings in HSC biology.

Single Cell RNA Sequencing in Autoimmune Inflammatory Rheumatic Diseases: Current Applications, Challenges and a Step Toward Precision Medicine

Single cell RNA sequencing (scRNA-seq) represents a new large scale and high throughput technique allowing analysis of the whole transcriptome at the resolution of an individual cell. It has emerged as an imperative method in life science research, uncovering complex cellular networks and providing indices that will eventually lead to the development of more targeted and personalized therapies. The importance of scRNA-seq has been particularly highlighted through the analysis of complex biological systems, in which cellular heterogeneity is a key aspect, such as the immune system. Autoimmune inflammatory rheumatic diseases represent a group of disorders, associated with a dysregulated immune system and high patient heterogeneity in both pathophysiological and clinical aspects. This complicates the complete understanding of underlying pathological mechanisms, associated with limited therapeutic options available and their long-term inefficiency and even toxicity. There is an unmet need to investigate, in depth, the cellular and molecular mechanisms driving the pathogenesis of rheumatic diseases and drug resistance, identify novel therapeutic targets, as well as make a step forward in using stratified and informed therapeutic decisions, which could now be achieved with the use of single cell approaches. This review summarizes the current use of scRNA-seq in studying different rheumatic diseases, based on recent findings from published in vitro, in vivo, and clinical studies, as well as discusses the potential implementation of scRNA-seq in the development of precision medicine in rheumatology.

Advances in genomics and genome editing for breeding next generation of fruit and nut crops

Fruit tree crops are an essential part of the food production systems and are key to achieve food and nutrition security. Genetic improvement of fruit trees by conventional breeding has been slow due to the long juvenile phase. Advancements in genomics and molecular biology have paved the way for devising novel genetic improvement tools like genome editing, which can accelerate the breeding of these perennial crops to a great extent. In this article, advancements in genomics of fruit trees covering genome sequencing, transcriptome sequencing, genome editing technologies (GET), CRISPR-Cas system based genome editing, potential applications of CRISPR-Cas9 in fruit tree crops improvement, the factors influencing the CRISPR-Cas editing efficiency and the challenges for CRISPR-Cas9 applications in fruit tree crops improvement are reviewed. Besides, base editing, a recently emerging more precise editing system, and the future perspectives of genome editing in the improvement of fruit and nut crops are covered.

Single-cell RNA sequencing in Drosophila: Technologies and applications

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful tool for investigating cell states and functions at the single-cell level. It has greatly revolutionized transcriptomic studies in many life science research fields, such as neurobiology, immunology, and developmental biology. With the fast development of both experimental platforms and bioinformatics approaches over the past decade, scRNA-seq is becoming economically feasible and experimentally practical for many biomedical laboratories. Drosophila has served as an excellent model organism for dissecting cellular and molecular mechanisms that underlie tissue development, adult cell function, disease, and aging. The recent application of scRNA-seq methods to Drosophila tissues has led to a number of exciting discoveries. In this review, I will provide a summary of recent scRNA-seq studies in Drosophila, focusing on technical approaches and biological applications. I will also discuss current challenges and future opportunities of making new discoveries using scRNA-seq in Drosophila. This article is categorized under:   Technologies > Analysis of the Transcriptome.

CALLR: a semi-supervised cell-type annotation method for single-cell RNA sequencing data

Motivation

Single-cell RNA sequencing (scRNA-seq) technology has been widely applied to capture the heterogeneity of different cell types within complex tissues. An essential step in scRNA-seq data analysis is the annotation of cell types. Traditional cell-type annotation is mainly clustering the cells first, and then using the aggregated cluster-level expression profiles and the marker genes to label each cluster. Such methods are greatly dependent on the clustering results, which are insufficient for accurate annotation.

Results

In this article, we propose a semi-supervised learning method for cell-type annotation called CALLR. It combines unsupervised learning represented by the graph Laplacian matrix constructed from all the cells and supervised learning using sparse logistic regression. By alternately updating the cell clusters and annotation labels, high annotation accuracy can be achieved. The model is formulated as an optimization problem, and a computationally efficient algorithm is developed to solve it. Experiments on 10 real datasets show that CALLR outperforms the compared (semi-)supervised learning methods, and the popular clustering methods.

Availability and implementation

The implementation of CALLR is available at https://github.com/MathSZhang/CALLR.

Supplementary information

Supplementary data are available at Bioinformatics online.

MP3RNA-seq: Massively parallel 3' end RNA sequencing for high-throughput gene expression profiling and genotyping

Transcriptome deep sequencing (RNA-seq) has become a routine method for global gene expression profiling. However, its application to large-scale experiments remains limited by cost and labor constraints. Here we describe a massively parallel 3' end RNA-seq (MP3RNA-seq) method that introduces unique sample barcodes during reverse transcription to permit sample pooling immediately following this initial step. MP3RNA-seq allows for handling of hundreds of samples in a single experiment, at a cost of about $6 per sample for library construction and sequencing. MP3RNA-seq is effective for not only high-throughput gene expression profiling, but also genotyping. To demonstrate its utility, we applied MP3RNA-seq to 477 double haploid lines of maize. We identified 19,429 genes expressed in at least 50% of the lines and 35,836 high-quality single nucleotide polymorphisms for genotyping analysis. Armed with these data, we performed expression and agronomic trait quantitative trait locus (QTL) mapping and identified 25,797 expression QTLs for 15,335 genes and 21 QTLs for plant height, ear height, and relative ear height. We conclude that MP3RNA-seq is highly reproducible, accurate, and sensitive for high-throughput gene expression profiling and genotyping, and should be generally applicable to most eukaryotic species.

Reference Transcriptomes of Porcine Peripheral Immune Cells Created Through Bulk and Single-Cell RNA Sequencing

Pigs are a valuable human biomedical model and an important protein source supporting global food security. The transcriptomes of peripheral blood immune cells in pigs were defined at the bulk cell-type and single cell levels. First, eight cell types were isolated in bulk from peripheral blood mononuclear cells (PBMCs) by cell sorting, representing Myeloid, NK cells and specific populations of T and B-cells. Transcriptomes for each bulk population of cells were generated by RNA-seq with 10,974 expressed genes detected. Pairwise comparisons between cell types revealed specific expression, while enrichment analysis identified 1,885 to 3,591 significantly enriched genes across all 8 cell types. Gene Ontology analysis for the top 25% of significantly enriched genes (SEG) showed high enrichment of biological processes related to the nature of each cell type. Comparison of gene expression indicated highly significant correlations between pig cells and corresponding human PBMC bulk RNA-seq data available in Haemopedia. Second, higher resolution of distinct cell populations was obtained by single-cell RNA-sequencing (scRNA-seq) of PBMC. Seven PBMC samples were partitioned and sequenced that produced 28,810 single cell transcriptomes distributed across 36 clusters and classified into 13 general cell types including plasmacytoid dendritic cells (DC), conventional DCs, monocytes, B-cell, conventional CD4 and CD8 αβ T-cells, NK cells, and γδ T-cells. Signature gene sets from the human Haemopedia data were assessed for relative enrichment in genes expressed in pig cells and integration of pig scRNA-seq with a public human scRNA-seq dataset provided further validation for similarity between human and pig data. The sorted porcine bulk RNAseq dataset informed classification of scRNA-seq PBMC populations; specifically, an integration of the datasets showed that the pig bulk RNAseq data helped define the CD4CD8 double-positive T-cell populations in the scRNA-seq data. Overall, the data provides deep and well-validated transcriptomic data from sorted PBMC populations and the first single-cell transcriptomic data for porcine PBMCs. This resource will be invaluable for annotation of pig genes controlling immunogenetic traits as part of the porcine Functional Annotation of Animal Genomes (FAANG) project, as well as further study of, and development of new reagents for, porcine immunology.

Tutorial: guidelines for annotating single-cell transcriptomic maps using automated and manual methods

Single-cell transcriptomics can profile thousands of cells in a single experiment and identify novel cell types, states and dynamics in a wide variety of tissues and organisms. Standard experimental protocols and analysis workflows have been developed to create single-cell transcriptomic maps from tissues. This tutorial focuses on how to interpret these data to identify cell types, states and other biologically relevant patterns with the objective of creating an annotated map of cells. We recommend a three-step workflow including automatic cell annotation (wherever possible), manual cell annotation and verification. Frequently encountered challenges are discussed, as well as strategies to address them. Guiding principles and specific recommendations for software tools and resources that can be used for each step are covered, and an R notebook is included to help run the recommended workflow. Basic familiarity with computer software is assumed, and basic knowledge of programming (e.g., in the R language) is recommended.

Alternative Polyadenylation in Stem Cell Self-Renewal and Differentiation

Cellular function is shaped by transcriptional and post-transcriptional mechanisms, including alternative polyadenylation (APA). By directly controlling 3'- untranslated region (UTR) length and the selection of the last exon, APA regulates up to 70% of all cellular transcripts influencing RNA stability, output, and protein isoform expression. Cell-state-dependent 3'-UTR shortening has been identified as a hallmark of cellular proliferation. Hence, quiescent/dormant stem cells are characterized by long 3'-UTRs, whereas proliferative stem/progenitor cells exhibit 3'-UTR shortening. Here, the latest studies analyzing the role of APA in regulating stem cell state, self-renewal, differentiation, and metabolism are reviewed. The new role of APA in controlling stem cell fate opens novel potential therapeutic avenues in the field of regenerative medicine.

The relationship between the gut microbiome and host gene expression: a review

Despite the growing knowledge surrounding host-microbiome interactions, we are just beginning to understand how the gut microbiome influences-and is influenced by-host gene expression. Here, we review recent literature that intersects these two fields, summarizing themes across studies. Work in model organisms, human biopsies, and cell culture demonstrate that the gut microbiome is an important regulator of several host pathways relevant for disease, including immune development and energy metabolism, and vice versa. The gut microbiome remodels host chromatin, causes differential splicing, alters the epigenetic landscape, and directly interrupts host signaling cascades. Emerging techniques like single-cell RNA sequencing and organoid generation have the potential to refine our understanding of the relationship between the gut microbiome and host gene expression in the future. By intersecting microbiome and host gene expression, we gain a window into the physiological processes important for fostering the extensive cross-kingdom interactions and ultimately our health.

Single-Cell Transcriptome Analysis of Human Adipose-Derived Stromal Cells Identifies a Contractile Cell Subpopulation

The potential for human adipose-derived stromal cells (hASCs) to be used as a therapeutic product is being assessed in multiple clinical trials. However, much is still to be learned about these cells before they can be used with confidence in the clinical setting. An inherent characteristic of hASCs that is not well understood is their heterogeneity. The aim of this exploratory study was to characterize the heterogeneity of freshly isolated hASCs after two population doublings (P2) using single-cell transcriptome analysis. A minimum of two subpopulations were identified at P2. A major subpopulation was identified as contractile cells which, based on gene expression patterns, are likely to be pericytes and/or vascular smooth muscle cells (vSMCs).

Uncovering transcriptional dark matter via gene annotation independent single-cell RNA sequencing analysis

Conventional scRNA-seq expression analyses rely on the availability of a high quality genome annotation. Yet, as we show here with scRNA-seq experiments and analyses spanning human, mouse, chicken, mole rat, lemur and sea urchin, genome annotations are often incomplete, in particular for organisms that are not routinely studied. To overcome this hurdle, we created a scRNA-seq analysis routine that recovers biologically relevant transcriptional activity beyond the scope of the best available genome annotation by performing scRNA-seq analysis on any region in the genome for which transcriptional products are detected. Our tool generates a single-cell expression matrix for all transcriptionally active regions (TARs), performs single-cell TAR expression analysis to identify biologically significant TARs, and then annotates TARs using gene homology analysis. This procedure uses single-cell expression analyses as a filter to direct annotation efforts to biologically significant transcripts and thereby uncovers biology to which scRNA-seq would otherwise be in the dark.

Conventional single-cell RNA sequencing analysis rely on genome annotations that may be incomplete or inaccurate especially for understudied organisms. Here the authors present a bioinformatic tool that leverages single-cell data to uncover biologically relevant transcripts beyond the best available genome annotation.

Imputation of spatially-resolved transcriptomes by graph-regularized tensor completion

High-throughput spatial-transcriptomics RNA sequencing (sptRNA-seq) based on in-situ capturing technologies has recently been developed to spatially resolve transcriptome-wide mRNA expressions mapped to the captured locations in a tissue sample. Due to the low RNA capture efficiency by in-situ capturing and the complication of tissue section preparation, sptRNA-seq data often only provides an incomplete profiling of the gene expressions over the spatial regions of the tissue. In this paper, we introduce a graph-regularized tensor completion model for imputing the missing mRNA expressions in sptRNA-seq data, namely FIST, Fast Imputation of Spatially-resolved transcriptomes by graph-regularized Tensor completion. We first model sptRNA-seq data as a 3-way sparse tensor in genes (p-mode) and the (x, y) spatial coordinates (x-mode and y-mode) of the observed gene expressions, and then consider the imputation of the unobserved entries or fibers as a tensor completion problem in Canonical Polyadic Decomposition (CPD) form. To improve the imputation of highly sparse sptRNA-seq data, we also introduce a protein-protein interaction network to add prior knowledge of gene functions, and a spatial graph to capture the the spatial relations among the capture spots. The tensor completion model is then regularized by a Cartesian product graph of protein-protein interaction network and the spatial graph to capture the high-order relations in the tensor. In the experiments, FIST was tested on ten 10x Genomics Visium spatial transcriptomic datasets of different tissue sections with cross-validation among the known entries in the imputation. FIST significantly outperformed the state-of-the-art methods for single-cell RNAseq data imputation. We also demonstrate that both the spatial graph and PPI network play an important role in improving the imputation. In a case study, we further analyzed the gene clusters obtained from the imputed gene expressions to show that the imputations by FIST indeed capture the spatial characteristics in the gene expressions and reveal functions that are highly relevant to three different kinds of tissues in mouse kidney.

Biological tissues are composed of different types of structurally organized cell units playing distinct functional roles. The exciting new spatial gene expression profiling methods have enabled the analysis of spatially resolved transcriptomes to understand the spatial and functional characteristics of these cells in the context of eco-environment of tissue. Due to the technical limitations, spatial transcriptomics data suffers from only sparsely measured mRNAs by in-situ capture and possibly missing spots in tissue regions that entirely failed fixing and permeabilizing RNAs. Our method, FIST (Fast Imputation of Spatially-resolved transcriptomes by graph-regularized Tensor completion), focuses on the spatial and high-sparsity nature of spatial transcriptomics data by modeling the data as a 3-way gene-by-(x, y)-location tensor and a product graph of a spatial graph and a protein-protein interaction network. Our comprehensive evaluation of FIST on ten 10x Genomics Visium spatial genomics datasets and comparison with the methods for single-cell RNA sequencing data imputation demonstrate that FIST is a better method more suitable for spatial gene expression imputation. Overall, we found FIST a useful new method for analyzing spatially resolved gene expressions based on novel modeling of spatial and functional information.

Microdroplet-based one-step RT-PCR for ultrahigh throughput single-cell multiplex gene expression analysis and rare cell detection

Gene expression analysis of individual cells enables characterization of heterogeneous and rare cell populations, yet widespread implementation of existing single-cell gene analysis techniques has been hindered due to limitations in scale, ease, and cost. Here, we present a novel microdroplet-based, one-step reverse-transcriptase polymerase chain reaction (RT-PCR) platform and demonstrate the detection of three targets simultaneously in over 100,000 single cells in a single experiment with a rapid read-out. Our customized reagent cocktail incorporates the bacteriophage T7 gene 2.5 protein to overcome cell lysate-mediated inhibition and allows for one-step RT-PCR of single cells encapsulated in nanoliter droplets. Fluorescent signals indicative of gene expressions are analyzed using a probabilistic deconvolution method to account for ambient RNA and cell doublets and produce single-cell gene signature profiles, as well as predict cell frequencies within heterogeneous samples. We also developed a simulation model to guide experimental design and optimize the accuracy and precision of the assay. Using mixtures of in vitro transcripts and murine cell lines, we demonstrated the detection of single RNA molecules and rare cell populations at a frequency of 0.1%. This low cost, sensitive, and adaptable technique will provide an accessible platform for high throughput single-cell analysis and enable a wide range of research and clinical applications.

Multi-omics approaches to improve malaria therapy

Malaria contributes to the most widespread infectious diseases worldwide. Even though current drugs are commercially available, the ever-increasing drug resistance problem by malaria parasites poses new challenges in malaria therapy. Hence, searching for efficient therapeutic strategies is of high priority in malaria control. In recent years, multi-omics technologies have been extensively applied to provide a more holistic view of functional principles and dynamics of biological mechanisms. We briefly review multi-omics technologies and focus on recent malaria progress conducted with the help of various omics methods. Then, we present up-to-date advances for multi-omics approaches in malaria. Next, we describe resistance phenomena to established antimalarial drugs and underlying mechanisms. Finally, we provide insight into novel multi-omics approaches, new drugs and vaccine developments and analyze current gaps in multi-omics research. Although multi-omics approaches have been successfully used in malaria studies, they are still limited. Many gaps need to be filled to bridge the gap between basic research and treatment of malaria patients. Multi-omics approaches will foster a better understanding of the molecular mechanisms of Plasmodium that are essential for the development of novel drugs and vaccines to fight this disastrous disease.

Tracing production instability in a clonally derived CHO cell line using single-cell transcriptomics

A variety of mechanisms including transcriptional silencing, gene copy loss, and increased susceptibility to cellular stress have been associated with a sudden or gradual loss of monoclonal antibody (mAb) production in Chinese hamster ovary (CHO) cell lines. In this study, we utilized single-cell RNA-seq (scRNA-seq) to study a clonally derived CHO cell line that underwent production instability leading to a dramatic reduction of the levels of mAb produced. From the scRNA-seq data, we identified subclusters associated with variations in the mAb transgenes and observed that heavy chain gene expression was significantly lower than that of the light chain across the population. Using trajectory inference, the evolution of the cell line was reconstructed and was found to correlate with a reduction in heavy and light chain gene expression. Genes encoding for proteins involved in the response to oxidative stress and apoptosis were found to increase in expression as cells progressed along the trajectory. Future studies of CHO cell lines using this technology have the potential to dramatically enhance our understanding of the characteristics underpinning efficient manufacturing performance as well as product quality.

Benchmarking Computational Doublet-Detection Methods for Single-Cell RNA Sequencing Data

In single-cell RNA sequencing (scRNA-seq), doublets form when two cells are encapsulated into one reaction volume. The existence of doublets, which appear to be-but are not-real cells, is a key confounder in scRNA-seq data analysis. Computational methods have been developed to detect doublets in scRNA-seq data; however, the scRNA-seq field lacks a comprehensive benchmarking of these methods, making it difficult for researchers to choose an appropriate method for specific analyses. We conducted a systematic benchmark study of nine cutting-edge computational doublet-detection methods. Our study included 16 real datasets, which contained experimentally annotated doublets, and 112 realistic synthetic datasets. We compared doublet-detection methods regarding detection accuracy under various experimental settings, impacts on downstream analyses, and computational efficiencies. Our results show that existing methods exhibited diverse performance and distinct advantages in different aspects. Overall, the DoubletFinder method has the best detection accuracy, and the cxds method has the highest computational efficiency. A record of this paper's transparent peer review process is included in the Supplemental Information.

Altered visual processing in the mdx52 mouse model of Duchenne muscular dystrophy

The mdx52 mouse model of Duchenne muscular dystrophy (DMD) is lacking exon 52 of the DMD gene that is located in a hotspot mutation region causing cognitive deficits and retinal anomalies in DMD patients. This deletion leads to the loss of the dystrophin proteins, Dp427, Dp260 and Dp140, while Dp71 is preserved. The flash electroretinogram (ERG) in mdx52 mice was previously characterized by delayed dark-adapted b-waves. A detailed description of functional ERG changes and visual performances in mdx52 mice is, however, lacking. Here an extensive full-field ERG repertoire was applied in mdx52 mice and WT littermates to analyze retinal physiology in scotopic, mesopic and photopic conditions in response to flash, sawtooth and/or sinusoidal stimuli. Behavioral contrast sensitivity was assessed using quantitative optomotor response (OMR) to sinusoidally modulated luminance gratings at 100% or 50% contrast. The mdx52 mice exhibited reduced amplitudes and delayed implicit times in dark-adapted ERG flash responses, particularly in their b-wave and oscillatory potentials, and diminished amplitudes of light-adapted flash ERGs. ERG responses to sawtooth stimuli were also diminished and delayed for both mesopic and photopic conditions in mdx52 mice and the first harmonic amplitudes to photopic sine-wave stimuli were smaller at all temporal frequencies. OMR indices were comparable between genotypes at 100% contrast but significantly reduced in mdx52 mice at 50% contrast. The complex ERG alterations and disturbed contrast vision in mdx52 mice include features observed in DMD patients and suggest altered photoreceptor-to-bipolar cell transmission possibly affecting contrast sensitivity. The mdx52 mouse is a relevant model to appraise the roles of retinal dystrophins and for preclinical studies related to DMD.

A heterogeneous cellular response to ionizing radiation revealed by single cell transcriptome sequencing

Our understanding on transcriptional regulation of tumour cells responding to ionizing radiation (IR) has mostly come from bulk sequencing. However, due to the heterogeneity of tumour, how each individual cell responds to IR differently is unclear. We report here a heterogeneous cellular response to IR by single cell transcriptome sequencing. We utilized the barcoded Smart-seq2 single cell transcriptome sequencing technology in breast cancer cell line MDA-MB-231 both without and with IR treatment. To further understand how ATM, a major hub protein required for an optimal DNA damage response, affected the heterogeneous IR response, we also knocked down ATM gene for single cell transcriptome sequencing. Single cell t-SNE analysis showed four clusters of cells responding to IR in distinctive ways: Cluster 1 changed the least; Cluster 2 responded to IR by upregulating ribosome associated genes, while Cluster 4 upregulated both ribosome and G1/S phase associated genes; Cluster 3 was a new cluster, which appeared only in irradiated cells. In the absence of ATM kinase, cells displayed much less transcriptional changes after IR. And Cluster 4 in wild-type cells, which had the greatest change in response to IR, was not present in the ATM knock-down cells. We also selected three IR-induced genes for functional validation in both MDA-MB-231 and an additional breast cancer cell line to demonstrate their importance in radiation sensitivity. Taken together, our single cell transcriptome analysis has revealed a heterogeneous cellular response to DNA damage induced by IR and identified potential biomarkers of radiation sensitivity.