Single-cell Samples Research Articles

Single-Cell Lipidomics: An Automated and Accessible Microfluidic Workflow Validated by Capillary Sampling.

We report the first demonstration of a microfluidics-based approach to measure lipids in single living cells using widely available liquid chromatography mass spectrometry (LC-MS) instrumentation. The method enables the rapid sorting of live cells into liquid chambers formed on standard Petri dishes and their subsequent dispensing into vials for analysis using LC-MS. This approach facilitates automated sampling, data acquisition, and analysis and carries the additional advantage of chromatographic separation, aimed at reducing matrix effects present in shotgun lipidomics approaches. We demonstrate that our method detects comparable numbers of features at around 200 lipids in populations of single cells versus established live single-cell capillary sampling methods and with greater throughput, albeit with the loss of spatial resolution. We also show the importance of optimization steps in addressing challenges from lipid contamination, especially in blanks, and demonstrate a 75% increase in the number of lipids identified. This work opens up a novel, accessible, and high-throughput way to obtain single-cell lipid profiles and also serves as an important validation of single-cell lipidomics through the use of different sampling methods.

Artificial intelligence-enabled multipurpose smart detection in active-matrix electrowetting-on-dielectric digital microfluidics

An active-matrix electrowetting-on-dielectric (AM-EWOD) system integrates hundreds of thousands of active electrodes for sample droplet manipulation, which can enable simultaneous, automatic, and parallel on-chip biochemical reactions. A smart detection system is essential for ensuring a fully automatic workflow and online programming for the subsequent experimental steps. In this work, we demonstrated an artificial intelligence (AI)-enabled multipurpose smart detection method in an AM-EWOD system for different tasks. We employed the U-Net model to quantitatively evaluate the uniformity of the applied droplet-splitting methods. We used the YOLOv8 model to monitor the droplet-splitting process online. A 97.76% splitting success rate was observed with 18 different AM-EWOD chips. A 99.982% model precision rate and a 99.980% model recall rate were manually verified. We employed an improved YOLOv8 model to detect single-cell samples in nanolitre droplets. Compared with manual verification, the model achieved 99.260% and 99.193% precision and recall rates, respectively. In addition, single-cell droplet sorting and routing experiments were demonstrated. With an AI-based smart detection system, AM-EWOD has shown great potential for use as a ubiquitous platform for implementing true lab-on-a-chip applications.

Mugen-UMAP: UMAP visualization and clustering of mutated genes in single-cell DNA sequencing data

BackgroundThe application of Uniform Manifold Approximation and Projection (UMAP) for dimensionality reduction and visualization has revolutionized the analysis of single-cell RNA expression and population genetics. However, its potential in single-cell DNA sequencing data analysis, particularly for visualizing gene mutation information, has not been fully explored.ResultsWe introduce Mugen-UMAP, a novel Python-based program that extends UMAP’s utility to single-cell DNA sequencing data. This innovative tool provides a comprehensive pipeline for processing gene annotation files of single-cell somatic single-nucleotide variants and metadata to the visualization of UMAP projections for identifying clusters, along with various statistical analyses. Employing Mugen-UMAP, we analyzed whole-exome sequencing data from 365 single-cell samples across 12 non-small cell lung cancer (NSCLC) patients, revealing distinct clusters associated with histological subtypes of NSCLC. Moreover, to demonstrate the general utility of Mugen-UMAP, we applied the program to 9 additional single-cell WES datasets from various cancer types, uncovering interesting patterns of cell clusters that warrant further investigation. In summary, Mugen-UMAP provides a quick and effective visualization method to uncover cell cluster patterns based on the gene mutation information from single-cell DNA sequencing data.ConclusionsThe application of Mugen-UMAP demonstrates its capacity to provide valuable insights into the visualization and interpretation of single-cell DNA sequencing data. Mugen-UMAP can be found at https://github.com/tengchn/Mugen-UMAP

Detecting expressed genes in cell populations at the single-cell level with scGeneXpress.

Determining whether genes are expressed or not remains a challenge in single-cell RNAseq experiments due to their different expression spectra, which are influenced by genetics, the microenvironment and gene length. Current approaches for addressing this issue fail to provide a comprehensive landscape of expressed genes, since they neglect the inherent differences in the expression ranges and distributions of genes. Here, we present scGeneXpress, a method for detecting expressed genes in cell populations of single-cell RNAseq samples based on gene-specific reference distributions. We demonstrate that scGeneXpress accurately detects expressed cell markers and identity genes in 34 human and mouse tissues and can be employed to improve differential expression analysis of single-cell RNAseq data.

Miniaturized protein profiling permits targeted signaling pathway analysis in individual circulating tumor cells to improve personalized treatment

BackgroundTraditional genomic profiling and mutation analysis of single cells like Circulating Tumor Cells (CTCs) fails to capture post-translational and functional alterations of proteins, often leading to limited treatment efficacy. To overcome this gap, we developed a miniaturized ‘protein analysis on the single cell level’ workflow—baptized ZeptoCTC. It integrates established technologies for single-cell isolation with sensitive Reverse Phase Protein Array (RPPA) analysis, thus enabling the comprehensive assessment of multiple protein expression and activation in individual CTCs.MethodsThe ZeptoCTC workflow involves several critical steps. Firstly, individual cells are labeled and isolated. This is followed by cell lysis and the printing of true single cell lysate preparations onto a ZeptoChip using a modified micromanipulator, CellCelector™. The printed lysates then undergo fluorescence immunoassay RPPA protein detection using a ZeptoReader. Finally, signal quantification is carried out with Image J software, ensuring precise measurement of multiple protein levels.ResultsThe efficacy of ZeptoCTC was demonstrated through various applications. Initially, it was used for measuring EpCAM protein expression, a standard marker for CTC detection, revealing higher levels in single MCF-7 over MDA-MB-231 tumor cells. Furthermore, in Capivasertib (Akt-inhibitor)-treated MCF-7 single cells, ZeptoCTC detected a 2-fold increase in the pAkt/Akt ratio compared to control cells, and confirmed co-performed bulk-cell western blot analysis results. Notably, when applied to individual CTCs from metastasized breast cancer patients, ZeptoCTC revealed significant differences in protein activation levels, particularly in measured pAkt and pErk levels, compared to patient-matched WBCs. Moreover, it successfully differentiated between CTCs from patients with different Akt1 genotypes, highlighting its potential to determine the activation status of druggable cancer driving proteins for individual and targeted treatment decision making.ConclusionsThe ZeptoCTC workflow represents a valuable tool in single cell cancer research, crucial for personalized medicine. It permits detailed analysis of key proteins and their activation status of targeted, cancer-driven signaling pathways in single cell samples, aiding in understanding tumor response, progression, and treatment efficacy beyond bulk analysis. The method significantly advances clinical investigations in cancer, improving treatment precision and effectiveness. The workflow will be applicable to protein analysis on other types of single cells like relevant in stem cell, neuropathology and hemopoietic cell research.

Parallel sample processing for mass spectrometry-based single cell proteomics

BackgroundSingle cell mass spectrometry (scMS) has shown great promise for label free proteomics analysis recently. To present single cell samples for proteomics analysis by MS is not a trivial task. Existing methods rely on robotic liquid handlers to scale up sample preparation throughput. The cost associated with specialized equipment hinders the broad adoption of these workflows, and the sequential sample processing nature limits the ultimate throughput. ResultsIn this work, we report a parallel sample processing workflow that can simultaneously process 10 single cells without the need of robotic liquid handlers for scMS. This method utilized 3D printed microfluidic devices to form reagent arrays on a glass slide, and a magnetic beads-based streamlined sample processing workflow to present peptides for LC-MS detection. We optimized key operational parameters of the method and demonstrated the quantification consistency among 10 parallel processed samples. Finally, the utility of the method in differentiating cell lines and studying the proteome change induced by drug treatment were demonstrated. SignificanceThe present method allows parallel sample processing for single cells without the need of expensive liquid handlers, which has great potential to further improve throughput and decrease the barrier for single cell proteomics.

Deep Learning Powers Protein Identification from Precursor MS Information.

Proteome analysis currently heavily relies on tandem mass spectrometry (MS/MS), which does not fully utilize MS1 features, as many precursors remain unselected for MS/MS fragmentation, especially in the cases of low abundance samples and wide abundance dynamic range samples. Therefore, leveraging MS1 features as a complement to MS/MS has become an attractive option to improve the coverage of feature identification. Herein, we propose MonoMS1, an approach combining deep learning-based retention time, ion mobility, detectability prediction, and logistic regression-based scoring for MS1 feature identification. The approach achieved a significant increase in MS1 feature identification based on an E. coli data set. Application of MonoMS1 to data sets with wide dynamic range, such as human serum proteome samples, and with low sample abundance, such as single-cell proteome samples, enabled substantial complementation of MS/MS-based peptide and protein identification. This method opens a new avenue for proteomic analysis and can boost proteomic research on complex samples.

Distinct Immune Homeostasis Remodeling Patterns after HLA-Matched and Haploidentical Transplantation.

Allogeneic hematopoietic stem cell transplantation (allo-HSCT) is a widely used treatment for a variety of hematopoietic disorders, and also provides a valuable platform for investigating the development of donor-derived immune cells in recipients post-HSCT. The immune system remodels from the donor to the recipient during allo-HSCT. However, little is known about the cell profile alterations as donor homeostasis rebalances to recipient homeostasis following HSCT. Here, multi-omics technology is applied at both the singlecell and bulk sample levels, as well as spectrum flow cytometry and fluorescent transgenic mouse models, to dissect the dynamics of the rebalanced homeostatic immune system in recipients after allo-HSCT. The data reveal that all immune subpopulations observed in donors are successfully restored in recipients, though with varying levels of abundance. The remodeling of immune homeostasis exhibits different patterns in HLA-matched and haploidentical HSCT, highlighting distinct biases in T cell reconstitution from the central and peripheral pathways. Furthermore, ZNF683 is critical for maintaining the persistence and quiescence of CD8 T-cell in haploidentical HSCT. The research can serve as a foundation for developing novel strategies to induce immune tolerance.

Integrative genomic analyses identify neuroblastoma risk genes involved in neuronal differentiation

Genome-Wide Association Studies (GWAS) have been decisive in elucidating the genetic predisposition of neuroblastoma (NB). The majority of genetic variants identified in GWAS are found in non-coding regions, suggesting that they can be causative of pathogenic dysregulations of gene expression. Nonetheless, pinpointing the potential causal genes within implicated genetic loci remains a major challenge. In this study, we integrated NB GWAS and expression Quantitative Trait Loci (eQTL) data from adrenal gland to identify candidate genes impacting NB susceptibility. We found that ZMYM1, CBL, GSKIP and WDR81 expression was dysregulated by NB predisposing variants. We further investigated the functional role of the identified genes through computational analysis of RNA sequencing (RNA-seq) data from single-cell and whole-tissue samples of NB, neural crest, and adrenal gland tissues, as well as through in vitro differentiation assays in NB cell cultures. Our results indicate that dysregulation of ZMYM1, CBL, GSKIP, WDR81 may lead to malignant transformation by affecting early and late stages of normal program of neuronal differentiation. Our findings enhance the understanding of how specific genes contribute to NB pathogenesis by highlighting their influence on neuronal differentiation and emphasizing the impact of genetic risk variants on the regulation of genes involved in critical biological processes.

Enhancement of redox cycling stability of Ni/GDC cermets for intermediate-temperature SOFC anodes through Ge incorporation in the ceramic phase

In this study, Ge4+ was incorporated into the ceramic matrix of nickel-based anode materials, synthesizing anode materials NiO-Gd0.1Ce0.9-xGexO1.95 (x = 0.01, 0.04, 0.07) and NiO-Ce0.9Gd0.1O2–δ. Corresponding anode support samples were also prepared, designated as A-CGDC, A-C1Ge, A-C4Ge, and A-C7Ge. Additionally, single cell samples using the newly co-doped anode materials were prepared, identified as C-CGDC, C-C1Ge, C-C4Ge, and C-C7Ge. The research indicates that after doping the ceramic phase framework with Ge4+, these composite anode materials exhibit the potential to withstand more redox cycles than C-CGDC without catastrophic mechanical structural damage or degradation in electrochemical performance. Particularly, after incorporating 7mol% Ge4+ into the ceramic backbone, the issue of 86 % strain accumulation in the anode support was mitigated. Additionally, the ion doping ratio is closely linked to the stability of the single cell samples, with the A-C7Ge sample able to withstand up to 11 extreme recycling sessions at 700 °C without complete structural failure, improving structural stability by approximately 57 % compared to A-C1Ge. In terms of electrochemical output performance, the impedance of C-C1Ge prior to complete failure was nearly double that before oxidation-reduction, with a maximum power density reduction of about 26.47 %. As the co-doping concentration in the ceramic phase reached 7mol%, the impedance increase during single cell cycle testing was about four times smaller, with C-C7Ge 's impedance before mechanical failure only 37 % higher than before cycling, and a reduction in maximum power density of just 11.11 %. Moreover, during consecutive cycling tests, the average damage to the open-circuit voltage was only 0.02V, showing almost no negative impact from the cycling operations, and it maintained efficient output performance after multiple cycles. Therefore, the newly synthesized nickel-based anode materials, modified via co-doping methods within the ceramic phase, can be considered as potential composite anode materials for intermediate-temperature solid oxide fuel cells.

Data acquisition approaches for single cell proteomics.

Single-cell proteomics (SCP) aims to characterize the proteome of individual cells, providing insights into complex biological systems. It reveals subtle differences in distinct cellular populations that bulk proteome analysis may overlook, which is essential for understanding disease mechanisms and developing targeted therapies. Mass spectrometry (MS) methods in SCP allow the identification and quantification of thousands of proteins from individual cells. Two major challenges in SCP are the limited material in single-cell samples necessitating highly sensitive analytical techniques and the efficient processing of samples, as each biological sample requires thousands of single cell measurements. This review discusses MS advancements to mitigate these challenges using data-dependent acquisition (DDA) and data-independent acquisition (DIA). Additionally, we examine the use of short liquid chromatography gradients and sample multiplexing methods that increase the sample throughput and scalability of SCP experiments. We believe these methods will pave the way for improving our understanding of cellular heterogeneity and its implications for systems biology.

Label-Free Single-Cell Cancer Classification from the Spatial Distribution of Adhesion Contact Kinetics.

There is an increasing need for simple-to-use, noninvasive, and rapid tools to identify and separate various cell types or subtypes at the single-cell level with sufficient throughput. Often, the selection of cells based on their direct biological activity would be advantageous. These steps are critical in immune therapy, regenerative medicine, cancer diagnostics, and effective treatment. Today, live cell selection procedures incorporate some kind of biomolecular labeling or other invasive measures, which may impact cellular functionality or cause damage to the cells. In this study, we first introduce a highly accurate single-cell segmentation methodology by combining the high spatial resolution of a phase-contrast microscope with the adhesion kinetic recording capability of a resonant waveguide grating (RWG) biosensor. We present a classification workflow that incorporates the semiautomatic separation and classification of single cells from the measurement data captured by an RWG-based biosensor for adhesion kinetics data and a phase-contrast microscope for highly accurate spatial resolution. The methodology was tested with one healthy and six cancer cell types recorded with two functionalized coatings. The data set contains over 5000 single-cell samples for each surface and over 12,000 samples in total. We compare and evaluate the classification using these two types of surfaces (fibronectin and noncoated) with different segmentation strategies and measurement timespans applied to our classifiers. The overall classification performance reached nearly 95% with the best models showing that our proof-of-concept methodology could be adapted for real-life automatic diagnostics use cases. The label-free measurement technique has no impact on cellular functionality, directly measures cellular activity, and can be easily tuned to a specific application by varying the sensor coating. These features make it suitable for applications requiring further processing of selected cells.

Inference of chromosome selection parameters and missegregation rate in cancer from DNA-sequencing data

Aneuploidy is frequently observed in cancers and has been linked to poor patient outcome. Analysis of aneuploidy in DNA-sequencing (DNA-seq) data necessitates untangling the effects of the Copy Number Aberration (CNA) occurrence rates and the selection coefficients that act upon the resulting karyotypes. We introduce a parameter inference algorithm that takes advantage of both bulk and single-cell DNA-seq cohorts. The method is based on Approximate Bayesian Computation (ABC) and utilizes CINner, our recently introduced simulation algorithm of chromosomal instability in cancer. We examine three groups of statistics to summarize the data in the ABC routine: (A) Copy Number-based measures, (B) phylogeny tip statistics, and (C) phylogeny balance indices. Using these statistics, our method can recover both the CNA probabilities and selection parameters from ground truth data, and performs well even for data cohorts of relatively small sizes. We find that only statistics in groups A and C are well-suited for identifying CNA probabilities, and only group A carries the signals for estimating selection parameters. Moreover, the low number of CNA events at large scale compared to cell counts in single-cell samples means that statistics in group B cannot be estimated accurately using phylogeny reconstruction algorithms at the chromosome level. As data from both bulk and single-cell DNA-sequencing techniques becomes increasingly available, our inference framework promises to facilitate the analysis of distinct cancer types, differentiation between selection and neutral drift, and prediction of cancer clonal dynamics.

Investigating the characteristics of mild intervertebral disc degeneration at various age stages using single-cell genomics.

Introduction: Intervertebral disc degeneration often occurs in the elderly population, but in recent years, there has been an increasing incidence of disc degeneration in younger individuals, primarily with mild degeneration. Methods: In order to explore the underlying mechanisms of disc degeneration in both young and aging individuals, we collected four types of nucleus pulposus (NP) single-cell sequencing samples for analysis based on Pfirrmann grading: normal-young (NY) (Grade I), normal-old (NO) (Grade I), mild degenerative-young (MY) (Grade II-III), and mild degenerative-old (MO) (Grade II-III). Results: We found that most NP cells in NO and MY samples exhibited oxidative stress, which may be important pathogenic factors in NO and MY groups. On the other hand, NP cells in MO group exhibited endoplasmic reticulum stress. In terms of inflammation, myeloid cells were mainly present in the degenerative group, with the MY group showing a stronger immune response compared to the MO group. Interestingly, dendritic cells in the myeloid lineage played a critical role in the process of mild degeneration. Discussion: Our study investigated the molecular mechanisms of intervertebral disc degeneration from an age perspective, providing insights for improving treatment strategies for patients with disc degeneration at different age groups.

Accurate assembly of multiple RNA-seq samples with Aletsch.

High-throughput RNA sequencing has become indispensable for decoding gene activities, yet the challenge of reconstructing full-length transcripts persists. Traditional single-sample assemblers frequently produce fragmented transcripts, especially in single-cell RNA-seq data. While algorithms designed for assembling multiple samples exist, they encounter various limitations. We present Aletsch, a new assembler for multiple bulk or single-cell RNA-seq samples. Aletsch incorporates several algorithmic innovations, including a "bridging" system that can effectively integrate multiple samples to restore missed junctions in individual samples, and a new graph-decomposition algorithm that leverages "supporting" information across multiple samples to guide the decomposition of complex vertices. A standout feature of Aletsch is its application of a random forest model with 50 well-designed features for scoring transcripts. We demonstrate its robust adaptability across different chromosomes, datasets, and species. Our experiments, conducted on RNA-seq data from several protocols, firmly demonstrate Aletsch's significant outperformance over existing meta-assemblers. As an example, when measured with the partial area under the precision-recall curve (pAUC, constrained by precision), Aletsch surpasses the leading assemblers TransMeta by 22.9%-62.1% and PsiCLASS by 23.0%-175.5% on human datasets. Aletsch is freely available at https://github.com/Shao-Group/aletsch. Scripts that reproduce the experimental results of this manuscript is available at https://github.com/Shao-Group/aletsch-test.

Refraction color based contact detection for single cell micromanipulation: a cornerstone for precise automation

Abstract Single cell studies represent a cutting-edge field in life science. For cells growing on a soft substrate, a precise micromanipulation is indispensable in order for high-quality single-cell sampling. The foundation of micromanipulation is the accurate triaxial movement control, of which the major challenge lies in the z-axis positioning. Conventional methods rely on manual operation but are hampered by low efficiency. Resorting to additional sensors and cameras might help, whereas the system configuration becomes inevitably complicated. In this work, a vision-based approach utilizing interface refraction color is reported. By monitoring the color change induced by refraction during the contact process, we realized an accurate movement control of the end-effector. The reported method is capable of detecting the z-axis contact between the capillary tip and the substrate in a fully automated manner. Tested by 70 batches of a total of 3840 images, the proposed technique achieves an overall accuracy of 94.28% with precision of 2.73 μm, sufficing in the applications of single cell sampling. Detecting z-axis contact by interface refraction color marks an innovative computer vision method and it is proposed and validated in this work. Our method holds a promise of being a potent part of a close-loop control mechanism for automated single cell micromanipulation.

Increasing Proteome Coverage Through a Reduction in Analyte Complexity in Single-Cell Equivalent Samples.

The advancement of sophisticated instrumentation in mass spectrometry has catalyzed an in-depth exploration of complex proteomes. This exploration necessitates a nuanced balance in experimental design, particularly between quantitative precision and the enumeration of analytes detected. In bottom-up proteomics, a key challenge is that oversampling of abundant proteins can adversely affect the identification of a diverse array of unique proteins. This issue is especially pronounced in samples with limited analytes, such as small tissue biopsies or single-cell samples. Methods such as depletion and fractionation are suboptimal to reduce oversampling in single cell samples, and other improvements on LC and mass spectrometry technologies and methods have been developed to address the trade-off between precision and enumeration. We demonstrate that by using a monosubstrate protease for proteomic analysis of single-cell equivalent digest samples, an improvement in quantitative accuracy can be achieved, while maintaining high proteome coverage established by trypsin. This improvement is particularly vital for the field of single-cell proteomics, where single-cell samples with limited number of protein copies, especially in the context of low-abundance proteins, can benefit from considering analyte complexity. Considerations about analyte complexity, alongside chromatographic complexity, integration with data acquisition methods, and other factors such as those involving enzyme kinetics, will be crucial in the design of future single-cell workflows.

Picoliter Single-Cell Reactor for Proteome Profiling by In Situ Cell Lysis, Protein Immobilization, Digestion, and Droplet Transfer.

Global profiling of single-cell proteomes can reveal cellular heterogeneity, thus benefiting precision medicine. However, current mass spectrometry (MS)-based single-cell proteomic sample processing still faces technical challenges associated with processing efficiency and protein recovery. Herein, we present an innovative sample processing platform based on a picoliter single-cell reactor (picoSCR) for single-cell proteome profiling, which involves in situ protein immobilization and sample transfer. PicoSCR helped minimize surface adsorptive losses by downscaling the processing volume to 400 pL with a contact area of less than 0.4 mm2. Besides, picoSCR reached highly efficient cell lysis and digestion within 30 min, benefiting from optimal reagent and high reactant concentrations. Using the picoSCR-nanoLC-MS system, over 1400 proteins were identified from an individual HeLa cell using data-dependent acquisition mode. Proteins with copy number below 1000 were identified, demonstrating this system with a detection limit of 1.7 zmol. Furthermore, we profiled the proteome of circulating tumor cells (CTCs). Data are available via ProteomeXchange with the identifier PXD051468. Proteins associated with epithelial-mesenchymal transition and neutrophil extracellular traps formation (which are both related to tumor metastasis) were observed in all CTCs. The cellular heterogeneity was revealed by differences in signaling pathways within individual cells. These results highlighted the potential of the picoSCR platform to help discover new biomarkers and explore differences in biological processes between cells.

ScCross: efficient search for rare subpopulations across multiple single-cell samples.

Identifying rare cell types is an important task to capture the heterogeneity of single-cell data, such as scRNA-seq. The widespread availability of such data enables to aggregate multiple samples, corresponding for example to different donors, into the same study. Yet, such aggregated data is often subject to batch effects between samples. Clustering it therefore generally requires the use of data integration methods, which can lead to overcorrection, making the identification of rare cells difficult. We present scCross, a biclustering method identifying rare subpopulations of cells present across multiple single-cell samples. It jointly identifies a group of cells with specific marker genes by relying on a global sum criterion, computed over entire subpopulation of cells, rather than pairwise comparisons between individual cells. This proves robust with respect to the high variability of scRNA-seq data, in particular batch effects. We show through several case studies that scCross is able to identify rare subpopulations across multiple samples without performing prior data integration. Namely, it identifies a cilium subpopulation with potential new ciliary genes from lung cancer cells, which is not detected by typical alternatives. It also highlights rare subpopulations in human pancreas samples sequenced with different protocols, despite visible shifts in expression levels between batches. We further show that scCross outperforms typical alternatives at identifying a target rare cell type in a controlled experiment with artificially created batch effects. This shows the ability of scCross to efficiently identify rare cell subpopulations characterized by specific genes despite the presence of batch effects. The R and Scala implementation of scCross is freely available on GitHub, at https://github.com/agerniers/scCross/. A snapshot of the code and the data underlying this article are available on Zenodo, at https://zenodo.org/doi/10.5281/zenodo.10471063.

TauFisher predicts circadian time from a single sample of bulk and single-cell pseudobulk transcriptomic data

As the circadian clock regulates fundamental biological processes, disrupted clocks are often observed in patients and diseased tissues. Determining the circadian time of the patient or the tissue of focus is essential in circadian medicine and research. Here we present tauFisher, a computational pipeline that accurately predicts circadian time from a single transcriptomic sample by finding correlations between rhythmic genes within the sample. We demonstrate tauFisher’s performance in adding timestamps to both bulk and single-cell transcriptomic samples collected from multiple tissue types and experimental settings. Application of tauFisher at a cell-type level in a single-cell RNAseq dataset collected from mouse dermal skin implies that greater circadian phase heterogeneity may explain the dampened rhythm of collective core clock gene expression in dermal immune cells compared to dermal fibroblasts. Given its robustness and generalizability across assay platforms, experimental setups, and tissue types, as well as its potential application in single-cell RNAseq data analysis, tauFisher is a promising tool that facilitates circadian medicine and research.