From multitude to singularity: An up-to-date overview of scRNA-seq data generation and analysis

Mitigating ambient RNA and doublets effects on single cell transcriptomics analysis in cancer research

In cancer biology, where understanding the tumor microenvironment at high resolution is vital, ambient RNA contamination becomes a considerable problem. This hinders accurate delineation of intratumoral heterogeneity, complicates the identification of potential biomarkers, and decelerates advancements in precision oncology. To solve this problem, several computational approaches are created to determine the ambient RNA contribution from scRNA-seq datasets. Techniques like SoupX and DecontX assist in assessing and eliminating ambient RNA contamination from primary gene expression profiles. Practical solutions like CellBender employ deep learning techniques to concurrently address ambient RNA contamination and background noise, offering a contemporary end-to-end strategy for data preparation. This high-quality, reliable data enables clinicians and researchers to make effective decisions that will help ensure interventions are rooted in reproducible evidence, giving hope for developing more effective targeted therapies.

Therapeutic and prognostic values of ferroptosis signature in glioblastoma

Ferroptosis is a regulated cell death process that results in decreased tumor growth and aggressiveness when targeted in various cancer cells. Studying the impact of ferroptosis in glioblastoma (GBM) will provide important knowledge about tumor biology and potential treatment strategies. The high metabolic activity resulting in ROS production, iron content and active lipid metabolism of glioblastoma cells make them particularly susceptible to ferroptosis. Single-cell RNA sequencing reveals the molecular signature of GBM and its tumor microenvironment, introducing ferroptosis-related biomarkers pathways and drug resistance mechanisms to enhance treatment outcomes for GBM patients. The relationship between ferroptosis and the immune landscape in GBM is complex and can have either positive or negative effects. These effects can be identified through single-cell RNA sequencing to develop targeted chemo-, radio- and immuno- therapies against glioma stem cells and tumor-supportive immune cells. Additionally, the implication of oncolytic virotherapy in combination with ferroptosis induction can lead to improved treatment of GBM in a clinical setting.

Single-cell RNA sequencing in ovarian cancer: Current progress and future prospects

Ovarian cancer is one of the most prevalent gynaecological malignancies. The rapid development of single-cell RNA sequencing (scRNA-seq) has allowed scientists to use this technique to study ovarian cancer development, heterogeneity, and tumour environment. Although multiple original research articles have reported the use of scRNA-seq in understanding ovarian cancer and how therapy resistance occurs, there is a lack of a comprehensive review that could summarize the findings from multiple studies. Therefore, this review aimed to fill this gap by comparing and summarizing the results from different studies that have used scRNA-seq in understanding ovarian cancer development, heterogeneity, tumour microenvironment, and treatment resistance. This review will begin with an overview of scRNA-seq workflow, followed by a discussion of various applications of scRNA-seq in studying ovarian cancer. Next, the limitations and future directions of scRNA-seq in ovarian cancer research will be presented.

High-throughput gene expression analysis with TempO-LINC sensitively resolves complex brain, lung and kidney heterogeneity at single-cell resolution

We report the development and performance of a novel genomics platform, TempO-LINC, for conducting high-throughput transcriptomic analysis on single cells and nuclei. TempO-LINC works by adding cell-identifying molecular barcodes onto highly selective and high-sensitivity gene expression probes within fixed cells, without having to first generate cDNA. Using an instrument-free combinatorial indexing approach, all probes within the same fixed cell receive an identical barcode, enabling the reconstruction of single-cell gene expression profiles across as few as several hundred cells and up to 100,000 + cells per sample. The TempO-LINC approach is easily scalable based on the number of barcodes and rounds of barcoding performed; however, for the experiments reported in this study, the assay utilized over 5.3 million unique barcodes. TempO-LINC offers a robust protocol for fixing and banking cells and displays high-sensitivity gene detection from multiple diverse sample types. We show that TempO-LINC has a multiplet rate of less than 1.1% and a cell capture rate of ~ 50%. Although the assay can accurately profile the whole transcriptome (19,683 human, 21,400 mouse and 21,119 rat genes), it can be targeted to measure only actionable/informative genes and molecular pathways of interest - thereby reducing sequencing requirements. In this study, we applied TempO-LINC to profile the transcriptomes of more than 90,000 cells across multiple species and sample types, including nuclei from mouse lung, kidney and brain tissues. The data demonstrated the ability to identify and annotate more than 50 unique cell populations and positively correlate expression of cell type-specific molecular markers within them. TempO-LINC is a robust new single-cell technology that is ideal for large-scale applications/studies with high data quality.

Comprehensive review on single-cell RNA sequencing: A new frontier in Alzheimer's disease research

Alzheimer's disease (AD) is a multifaceted neurodegenerative condition marked by gradual cognitive deterioration and the loss of neurons. While conventional bulk RNA sequencing techniques have shed light on AD pathology, they frequently obscure the cellular diversity within brain tissues. The advent of single-cell RNA sequencing (scRNA-seq) has transformed our capability to analyze the cellular composition of AD, allowing for the detection of unique cell populations, rare cell types, and gene expression alterations at an individual cell level. This review examines the use of scRNA-seq in AD research, focusing on its contributions to understanding cellular diversity, disease progression, and potential therapeutic targets. We discuss key technological innovations, data analysis techniques, and challenges associated with scRNA-seq in studying AD. Furthermore, we highlight recent studies that have utilized scRNA-seq to identify novel biomarkers, uncover disease-associated pathways, and elucidate the role of non-neuronal cells, such as microglia and astrocytes, in AD pathogenesis. By providing a comprehensive overview of advancements in scRNA-seq for unraveling cellular heterogeneity in AD, this review highlights the transformative impact of scRNA-seq on our comprehension of disease mechanisms and the creation of targeted treatments.

High-throughput gene expression analysis with TempO-LINC sensitively resolves complex brain, lung and kidney heterogeneity at single-cell resolution

We report the development and performance of a novel genomics platform, TempO-LINC, for conducting high-throughput transcriptomic analysis on single cells and nuclei. TempO-LINC works by adding cell-identifying molecular barcodes onto highly selective and high-sensitivity gene expression probes within fixed cells, without having to first generate cDNA. Using an instrument-free combinatorial-indexing approach, all probes within the same fixed cell receive an identical barcode, enabling the reconstruction of single-cell gene expression profiles across as few as several hundred cells and up to 100,000+ cells per run. The TempO-LINC approach is easily scalable based on the number of barcodes and rounds of barcoding performed; however, for the experiments reported in this study, the assay utilized over 5.3 million unique barcodes. TempO-LINC has a robust protocol for fixing and banking cells and displays high-sensitivity gene detection from multiple diverse sample types. We show that TempO-LINC has an observed multiplet rate of less than 1.1% and a cell capture rate of ~50%. Although the assay can accurately profile the whole transcriptome (19,683 human or 21,400 mouse genes), it can be targeted to measure only actionable/informative genes and molecular pathways of interest - thereby reducing sequencing requirements. In this study, we applied TempO-LINC to profile the transcriptomes of 89,722 cells across multiple sample types, including nuclei from mouse lung, kidney and brain tissues. The data demonstrated the ability to identify and annotate at least 50 unique cell populations and positively correlate expression of cell type-specific molecular markers within them. TempO-LINC is a robust new single-cell technology that is ideal for large-scale applications/studies across thousands of samples with high data quality.

Single-cell proteomics by mass spectrometry: Advances and implications in cancer research

Cancer harbours extensive proteomic heterogeneity. Inspired by the prior success of single-cell RNA sequencing (scRNA-seq) in characterizing minute transcriptomics heterogeneity in cancer, researchers are now actively searching for information regarding the proteomics counterpart. Therefore recently, single-cell proteomics by mass spectrometry (SCP) has rapidly developed into state-of-the-art technology to cater the need. This review aims to summarize application of SCP in cancer research, while revealing current development progress of SCP technology. The review also aims to contribute ideas into research gaps and future directions, ultimately promoting the application of SCP in cancer research.

scEM: A New Ensemble Framework for Predicting Cell Type Composition Based on scRNA-Seq Data

With the advent of single-cell RNA sequencing (scRNA-seq) technology, many scRNA-seq data have become available, providing an unprecedented opportunity to explore cellular composition and heterogeneity. Recently, many computational algorithms for predicting cell type composition have been developed, and these methods are typically evaluated on different datasets and performance metrics using diverse techniques. Consequently, the lack of comprehensive and standardized comparative analysis makes it difficult to gain a clear understanding of the strengths and weaknesses of these methods. To address this gap, we reviewed 20 cutting-edge unsupervised cell type identification methods and evaluated these methods comprehensively using 24 real scRNA-seq datasets of varying scales. In addition, we proposed a new ensemble cell-type identification method, named scEM, which learns the consensus similarity matrix by applying the entropy weight method to the four representative methods are selected. The Louvain algorithm is adopted to obtain the final classification of individual cells based on the consensus matrix. Extensive evaluation and comparison with 11 other similarity-based methods under real scRNA-seq datasets demonstrate that the newly developed ensemble algorithm scEM is effective in predicting cellular type composition.

Discovery of a Novel Intron in US10/US11/US12 of HSV-1 Strain 17

Herpes Simplex Virus type 1 (HSV-1) infects humans and causes a variety of clinical manifestations. Many HSV-1 genomes have been sequenced with high-throughput sequencing technologies and the annotation of these genome sequences heavily relies on the known genes in reference strains. Consequently, the accuracy of reference strain annotation is critical for future research and treatment of HSV-1 infection. In this study, we analyzed RNA-Seq data of HSV-1 from NCBI databases and discovered a novel intron in the overlapping coding sequence (CDS) of US10 and US11, and the 3' UTR of US12 in strain 17, a commonly used HSV-1 reference strain. To comprehensively understand the shared US10/US11/US12 intron structure, we used US11 as a representative and surveyed all US11 gene sequences from the NCBI nt/nr database. A total of 193 high-quality US11 sequences were obtained, of which 186 sequences have a domain of uninterrupted tandemly repeated RXP (Arg-X-Pro) in the C-terminus half of the protein. In total, 97 of the 186 sequences encode US11 protein with the same length of the mature US11 in strain 17:26 of them have the same structure of US11 and can be spliced as in strain 17; 71 of them have transcripts that are the same as mature US11 mRNA in strain 17. In total, 76 US11 gene sequences have either canonical or known noncanonical intron border sequences and may be spliced like strain 17 and obtain mature US11 CDS with the same length. If not spliced, they will have extra RXP repeats. A tandemly repeated RXP domain was proposed to be essential for US11 to bind with RNA and other host factors. US10 protein sequences from the same strains have also been studied. The results of this study show that even a frequently used reference organism may have errors in widely used databases. This study provides accurate annotation of the US10, US11, and US12 gene structure, which will build a more solid foundation to study expression regulation of the function of these genes.

Bioinformatics for wet-lab scientists: practical application in sequencing analysis

BACKGROUND

Genomics data is available to the scientific community after publication of research projects and can be investigated for a multitude of research questions. However, in many cases deposited data is only assessed and used for the initial publication, resulting in valuable resources not being exploited to their full depth. MAIN: A likely reason for this is that many wetlab-based researchers are not formally trained to apply bioinformatic tools and may therefore assume that they lack the necessary experience to do so themselves. In this article, we present a series of freely available, predominantly web-based platforms and bioinformatic tools that can be combined in analysis pipelines to interrogate different types of next-generation sequencing data. Additionally to the presented exemplary route, we also list a number of alternative tools that can be combined in a mix-and-match fashion. We place special emphasis on tools that can be followed and used correctly without extensive prior knowledge in programming. Such analysis pipelines can be applied to existing data downloaded from the public domain or be compared to the results of own experiments.

CONCLUSION

Integrating transcription factor binding to chromatin (ChIP-seq) with transcriptional output (RNA-seq) and chromatin accessibility (ATAC-seq) can not only assist to form a deeper understanding of the molecular interactions underlying transcriptional regulation but will also help establishing new hypotheses and pre-testing them in silico.

Cortical type: a conceptual tool for meaningful biological interpretation of high-throughput gene expression data in the human cerebral cortex

The interpretation of massive high-throughput gene expression data requires computational and biological analyses to identify statistically and biologically significant differences, respectively. There are abundant sources that describe computational tools for statistical analysis of massive gene expression data but few address data analysis for biological significance. In the present article we exemplify the importance of selecting the proper biological context in the human brain for gene expression data analysis and interpretation. For this purpose, we use cortical type as conceptual tool to make predictions about gene expression in areas of the human temporal cortex. We predict that the expression of genes related to glutamatergic transmission would be higher in areas of simpler cortical type, the expression of genes related to GABAergic transmission would be higher in areas of more complex cortical type, and the expression of genes related to epigenetic regulation would be higher in areas of simpler cortical type. Then, we test these predictions with gene expression data from several regions of the human temporal cortex obtained from the Allen Human Brain Atlas. We find that the expression of several genes shows statistically significant differences in agreement with the predicted gradual expression along the laminar complexity gradient of the human cortex, suggesting that simpler cortical types may have greater glutamatergic excitability and epigenetic turnover compared to more complex types; on the other hand, complex cortical types seem to have greater GABAergic inhibitory control compared to simpler types. Our results show that cortical type is a good predictor of synaptic plasticity, epigenetic turnover, and selective vulnerability in human cortical areas. Thus, cortical type can provide a meaningful context for interpreting high-throughput gene expression data in the human cerebral cortex.

RevGel-seq: instrument-free single-cell RNA sequencing using a reversible hydrogel for cell-specific barcoding

Progress in sample preparation for scRNA-seq is reported based on RevGel-seq, a reversible-hydrogel technology optimized for samples of fresh cells. Complexes of one cell paired with one barcoded bead are stabilized by a chemical linker and dispersed in a hydrogel in the liquid state. Upon gelation on ice the complexes are immobilized and physically separated without requiring nanowells or droplets. Cell lysis is triggered by detergent diffusion, and RNA molecules are captured on the adjacent barcoded beads for further processing with reverse transcription and preparation for cDNA sequencing. As a proof of concept, analysis of PBMC using RevGel-seq achieves results similar to microfluidic-based technologies when using the same original sample and the same data analysis software. In addition, a clinically relevant application of RevGel-seq is presented for pancreatic islet cells. Furthermore, characterizations carried out on cardiomyocytes demonstrate that the hydrogel technology readily accommodates very large cells. Standard analyses are in the 10,000-input cell range with the current gelation device, in order to satisfy common requirements for single-cell research. A convenient stopping point after two hours has been established by freezing at the cell lysis step, with full preservation of gene expression profiles. Overall, our results show that RevGel-seq represents an accessible and efficient instrument-free alternative, enabling flexibility in terms of experimental design and timing of sample processing, while providing broad coverage of cell types.

Analysis of candidate genes for cleft lip ± cleft palate using murine single-cell expression data

Introduction: Cleft lip ± cleft palate (CL/P) is one of the most common birth defects. Although research has identified multiple genetic risk loci for different types of CL/P (i.e., syndromic or non-syndromic forms), determining the respective causal genes and understanding the relevant functional networks remain challenging. The recent introduction of single-cell RNA sequencing (scRNA-seq) has provided novel opportunities to study gene expression patterns at cellular resolution. The aims of our study were to: (i) aggregate available scRNA-seq data from embryonic mice and provide this as a resource for the craniofacial community; and (ii) demonstrate the value of these data in terms of the investigation of the gene expression patterns of CL/P candidate genes. Methods and Results: First, two published scRNA-seq data sets from embryonic mice were re-processed, i.e., data representing the murine time period of craniofacial development: (i) facial data from embryonic day (E) E11.5; and (ii) whole embryo data from E9.5-E13.5 from the Mouse Organogenesis Cell Atlas (MOCA). Marker gene expression analyses demonstrated that at E11.5, the facial data were a high-resolution representation of the MOCA data. Using CL/P candidate gene lists, distinct groups of genes with specific expression patterns were identified. Among others we identified that a co-expression network including Irf6, Grhl3 and Tfap2a in the periderm, while it was limited to Irf6 and Tfap2a in palatal epithelia, cells of the ectodermal surface, and basal cells at the fusion zone. The analyses also demonstrated that additional CL/P candidate genes (e.g., Tpm1, Arid3b, Ctnnd1, and Wnt3) were exclusively expressed in Irf6+ facial epithelial cells (i.e., as opposed to Irf6- epithelial cells). The MOCA data set was finally used to investigate differences in expression profiles for candidate genes underlying different types of CL/P. These analyses showed that syndromic CL/P genes (syCL/P) were expressed in significantly more cell types than non-syndromic CL/P candidate genes (nsCL/P). Discussion: The present study illustrates how scRNA-seq data can empower research on craniofacial development and disease.