Genetic Variation and Stroke Recovery: The STRONG Study.

BACKGROUND: Genetic association studies can reveal biology and treatment targets but have received limited attention for stroke recovery. STRONG (Stroke, Stress, Rehabilitation, and Genetics) was a prospective, longitudinal (1-year), genetic study in adults with stroke at 28 US stroke centers. The primary aim was to examine the association that candidate genetic variants have with (1) motor/functional outcomes and (2) stress-related outcomes. METHODS: For motor/functional end points, 3 candidate gene variants (ApoE ε4, BDNF [brain-derived neurotrophic factor], and a dopamine polygenic score) were analyzed for associations with change in grip strength (3 months-baseline), function (3-month Stroke Impact Scale-Activities of Daily Living), mood (3-month Patient Health Questionnaire-8), and cognition (12-month telephone-Montreal Cognitive Assessment). For stress-related outcomes, 7 variants (serotonin transporter gene-linked promoter region, ACE [angiotensin-converting enzyme], oxytocin receptor, FKBP5 [FKBP prolyl isomerase 5], FAAH [fatty acid amide hydrolase], BDNF, and COMT [catechol-O-methyltransferase]) were assessed for associations with posttraumatic stress disorder ([PTSD]; PTSD Primary Care Scale) and depression (Patient Health Questionnaire-8) at 6 and 12 months; stress-related genes were examined as a function of poststroke stress level. Statistical models (linear, negative binomial, or Poisson regression) were based on response variable distribution; all included stroke severity, age, sex, and ancestry as covariates. Stroke subtype was explored secondarily. Data were Holm-Bonferroni corrected. A secondary replication analysis tested whether the rs1842681 polymorphism (identified in the GISCOME study [Genetics of Ischaemic Stroke Functional Outcome]) was related to 3-month modified Rankin Scale score in STRONG. RESULTS: The 763 enrollees were 63.1±14.9 (mean±SD) years of age, with a median initial National Institutes of Health Stroke Scale score of 4 (interquartile range, 2-9); outcome data were available in n=515 at 3 months, n=500 at 6 months, and n=489 at 12 months. At 1 year poststroke, the rs6265 (BDNF) variant was associated with poorer cognition (0.9-point lower telephone-Montreal Cognitive Assessment score, P=1×10-5). For stress-related outcomes, rs4291 (ACE) and rs324420 (FAAH) were risk factors linking increased poststroke stress with higher 1-year depression and PTSD symptoms (P<0.05), while rs4680 (COMT) linked poststroke stress with lower 1-year depression and PTSD. Findings were unchanged when considering stroke subtype. STRONG replicated GISCOME: rs1842681 was associated with lower 3-month modified Rankin Scale score (P=3.2×10-5). CONCLUSIONS: This study identified genetic associations with cognitive function, depression, and PTSD 1 year poststroke. Genetic susceptibility to PTSD and depressive symptoms varied according to the amount of poststroke stress, underscoring the critical role of lived experiences in recovery. Together, the results suggest that genetic association studies provide insights into the biology of stroke recovery in humans.

Electroencephalography Might Improve Diagnosis of Acute Stroke and Large Vessel Occlusion.

BACKGROUND AND PURPOSE: Clinical methods have incomplete diagnostic value for early diagnosis of acute stroke and large vessel occlusion (LVO). Electroencephalography is rapidly sensitive to brain ischemia. This study examined the diagnostic utility of electroencephalography for acute stroke/transient ischemic attack (TIA) and for LVO. METHODS: Patients (n=100) with suspected acute stroke in an emergency department underwent clinical exam then electroencephalography using a dry-electrode system. Four models classified patients, first as acute stroke/TIA or not, then as acute stroke with LVO or not: (1) clinical data, (2) electroencephalography data, (3) clinical+electroencephalography data using logistic regression, and (4) clinical+electroencephalography data using a deep learning neural network. Each model used a training set of 60 randomly selected patients, then was validated in an independent cohort of 40 new patients. RESULTS: Of 100 patients, 63 had a stroke (43 ischemic/7 hemorrhagic) or TIA (13). For classifying patients as stroke/TIA or not, the clinical data model had area under the curve=62.3, whereas clinical+electroencephalography using deep learning neural network model had area under the curve=87.8. Results were comparable for classifying patients as stroke with LVO or not. CONCLUSIONS: Adding electroencephalography data to clinical measures improves diagnosis of acute stroke/TIA and of acute stroke with LVO. Rapid acquisition of dry-lead electroencephalography is feasible in the emergency department and merits prehospital evaluation.

Evolutionary State-Space Model and Its Application to Time-Frequency Analysis of Local Field Potentials.

We propose an evolutionary state space model (E-SSM) for analyzing high dimensional brain signals whose statistical properties evolve over the course of a non-spatial memory experiment. Under E-SSM, brain signals are modeled as mixtures of components (e.g., AR(2) process) with oscillatory activity at pre-defined frequency bands. To account for the potential non-stationarity of these components (since the brain responses could vary throughout the entire experiment), the parameters are allowed to vary over epochs. Compared with classical approaches such as independent component analysis and filtering, the proposed method accounts for the entire temporal correlation of the components and accommodates non-stationarity. For inference purpose, we propose a novel computational algorithm based upon using Kalman smoother, maximum likelihood and blocked resampling. The E-SSM model is applied to simulation studies and an application to a multi-epoch local field potentials (LFP) signal data collected from a non-spatial (olfactory) sequence memory task study. The results confirm that our method captures the evolution of the power for different components across different phases in the experiment and identifies clusters of electrodes that behave similarly with respect to the decomposition of different sources. These findings suggest that the activity of different electrodes does change over the course of an experiment in practice; treating these epoch recordings as realizations of an identical process could lead to misleading results. In summary, the proposed method underscores the importance of capturing the evolution in brain responses over the study period.

Neural network gradient Hamiltonian Monte Carlo.

Hamiltonian Monte Carlo is a widely used algorithm for sampling from posterior distributions of complex Bayesian models. It can efficiently explore high-dimensional parameter spaces guided by simulated Hamiltonian flows. However, the algorithm requires repeated gradient calculations, and these computations become increasingly burdensome as data sets scale. We present a method to substantially reduce the computation burden by using a neural network to approximate the gradient. First, we prove that the proposed method still maintains convergence to the true distribution though the approximated gradient no longer comes from a Hamiltonian system. Second, we conduct experiments on synthetic examples and real data to validate the proposed method.

Prostatectomy-based validation of combined urine and plasma test for predicting high grade prostate cancer.

BACKGROUND: Distinguishing between low- and high-grade prostate cancers (PCa) is important, but biopsy may underestimate the actual grade of cancer. We have previously shown that urine/plasma-based prostate-specific biomarkers can predict high grade PCa. Our objective was to determine the accuracy of a test using cell-free RNA levels of biomarkers in predicting prostatectomy results. METHODS: This multicenter community-based prospective study was conducted using urine/blood samples collected from 306 patients. All recruited patients were treatment-naïve, without metastases, and had been biopsied, designated a Gleason Score (GS) based on biopsy, and assigned to prostatectomy prior to participation in the study. The primary outcome measure was the urine/plasma test accuracy in predicting high grade PCa on prostatectomy compared with biopsy findings. Sensitivity and specificity were calculated using standard formulas, while comparisons between groups were performed using the Wilcoxon Rank Sum, Kruskal-Wallis, Chi-Square, and Fisher's exact test. RESULTS: GS as assigned by standard 10-12 core biopsies was 3 + 3 in 90 (29.4%), 3 + 4 in 122 (39.8%), 4 + 3 in 50 (16.3%), and > 4 + 3 in 44 (14.4%) patients. The urine/plasma assay confirmed a previous validation and was highly accurate in predicting the presence of high-grade PCa (Gleason ≥3 + 4) with sensitivity between 88% and 95% as verified by prostatectomy findings. GS was upgraded after prostatectomy in 27% of patients and downgraded in 12% of patients. CONCLUSIONS: This plasma/urine biomarker test accurately predicts high grade cancer as determined by prostatectomy with a sensitivity at 92-97%, while the sensitivity of core biopsies was 78%.

Geodesic Lagrangian Monte Carlo over the space of positive definite matrices: with application to Bayesian spectral density estimation.

We present geodesic Lagrangian Monte Carlo, an extension of Hamiltonian Monte Carlo for sampling from posterior distributions defined on general Riemannian manifolds. We apply this new algorithm to Bayesian inference on symmetric or Hermitian positive definite matrices. To do so, we exploit the Riemannian structure induced by Cartan's canonical metric. The geodesics that correspond to this metric are available in closed-form and-within the context of Lagrangian Monte Carlo-provide a principled way to travel around the space of positive definite matrices. Our method improves Bayesian inference on such matrices by allowing for a broad range of priors, so we are not limited to conjugate priors only. In the context of spectral density estimation, we use the (non-conjugate) complex reference prior as an example modeling option made available by the algorithm. Results based on simulated and real-world multivariate time series are presented in this context, and future directions are outlined.

What time is it? Deep learning approaches for circadian rhythms.

MOTIVATION: Circadian rhythms date back to the origins of life, are found in virtually every species and every cell, and play fundamental roles in functions ranging from metabolism to cognition. Modern high-throughput technologies allow the measurement of concentrations of transcripts, metabolites and other species along the circadian cycle creating novel computational challenges and opportunities, including the problems of inferring whether a given species oscillate in circadian fashion or not, and inferring the time at which a set of measurements was taken. RESULTS: We first curate several large synthetic and biological time series datasets containing labels for both periodic and aperiodic signals. We then use deep learning methods to develop and train BIO_CYCLE, a system to robustly estimate which signals are periodic in high-throughput circadian experiments, producing estimates of amplitudes, periods, phases, as well as several statistical significance measures. Using the curated data, BIO_CYCLE is compared to other approaches and shown to achieve state-of-the-art performance across multiple metrics. We then use deep learning methods to develop and train BIO_CLOCK to robustly estimate the time at which a particular single-time-point transcriptomic experiment was carried. In most cases, BIO_CLOCK can reliably predict time, within approximately 1 h, using the expression levels of only a small number of core clock genes. BIO_CLOCK is shown to work reasonably well across tissue types, and often with only small degradation across conditions. BIO_CLOCK is used to annotate most mouse experiments found in the GEO database with an inferred time stamp. AVAILABILITY AND IMPLEMENTATION: All data and software are publicly available on the CircadiOmics web portal: circadiomics.igb.uci.edu/ CONTACTS: fagostin@uci.edu or pfbaldi@uci.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An efficient Bayesian inference framework for coalescent-based nonparametric phylodynamics.

MOTIVATION: The field of phylodynamics focuses on the problem of reconstructing population size dynamics over time using current genetic samples taken from the population of interest. This technique has been extensively used in many areas of biology but is particularly useful for studying the spread of quickly evolving infectious diseases agents, e.g. influenza virus. Phylodynamic inference uses a coalescent model that defines a probability density for the genealogy of randomly sampled individuals from the population. When we assume that such a genealogy is known, the coalescent model, equipped with a Gaussian process prior on population size trajectory, allows for nonparametric Bayesian estimation of population size dynamics. Although this approach is quite powerful, large datasets collected during infectious disease surveillance challenge the state-of-the-art of Bayesian phylodynamics and demand inferential methods with relatively low computational cost. RESULTS: To satisfy this demand, we provide a computationally efficient Bayesian inference framework based on Hamiltonian Monte Carlo for coalescent process models. Moreover, we show that by splitting the Hamiltonian function, we can further improve the efficiency of this approach. Using several simulated and real datasets, we show that our method provides accurate estimates of population size dynamics and is substantially faster than alternative methods based on elliptical slice sampler and Metropolis-adjusted Langevin algorithm. AVAILABILITY AND IMPLEMENTATION: The R code for all simulation studies and real data analysis conducted in this article are publicly available at http://www.ics.uci.edu/∼slan/lanzi/CODES.html and in the R package phylodyn available at https://github.com/mdkarcher/phylodyn. CONTACT: S.Lan@warwick.ac.uk or babaks@uci.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Neural function, injury, and stroke subtype predict treatment gains after stroke.

OBJECTIVE: This study was undertaken to better understand the high variability in response seen when treating human subjects with restorative therapies poststroke. Preclinical studies suggest that neural function, neural injury, and clinical status each influence treatment gains; therefore, the current study hypothesized that a multivariate approach incorporating these 3 measures would have the greatest predictive value. METHODS: Patients 3 to 6 months poststroke underwent a battery of assessments before receiving 3 weeks of standardized upper extremity robotic therapy. Candidate predictors included measures of brain injury (including to gray and white matter), neural function (cortical function and cortical connectivity), and clinical status (demographics/medical history, cognitive/mood, and impairment). RESULTS: Among all 29 patients, predictors of treatment gains identified measures of brain injury (smaller corticospinal tract [CST] injury), cortical function (greater ipsilesional motor cortex [M1] activation), and cortical connectivity (greater interhemispheric M1-M1 connectivity). Multivariate modeling found that best prediction was achieved using both CST injury and M1-M1 connectivity (r(2) = 0.44, p = 0.002), a result confirmed using Lasso regression. A threshold was defined whereby no subject with >63% CST injury achieved clinically significant gains. Results differed according to stroke subtype; gains in patients with lacunar stroke were best predicted by a measure of intrahemispheric connectivity. INTERPRETATION: Response to a restorative therapy after stroke is best predicted by a model that includes measures of both neural injury and function. Neuroimaging measures were the best predictors and may have an ascendant role in clinical decision making for poststroke rehabilitation, which remains largely reliant on behavioral assessments. Results differed across stroke subtypes, suggesting the utility of lesion-specific strategies.

Maternal Folate Concentration in Early Pregnancy and Newborn Telomere Length.

BACKGROUND/AIMS: Telomere biology plays a fundamental role in genomic integrity and cell physiology. The newborn setting of telomere length (TL) likely has important implications for telomere dynamics over the lifespan; however, its determinants are poorly understood. Folate is essential for DNA integrity. The maternal compartment is the only source of folate for the developing fetus. We, therefore, tested the hypothesis that variation in maternal folate during pregnancy is associated with newborn TL. METHODS: A prospective, longitudinal study was conducted in 119 mother-newborn dyads. Eligible mothers were enrolled at 9.5 (SD ±2.1) weeks gestation and followed through birth. Concentrations of maternal serum folate were measured in the first trimester of pregnancy. Newborn telomere length was measured in cord blood mononuclear cells (CBMC). RESULTS: After accounting for the effects of other established determinants of newborn TL, each 10 ng/ml increase in maternal total folate was associated with a 5.8% increase in median TL (p = 0.03). The median TL in newborns of mother in the lowest quartile of total folate levels was approximately 10% shorter than that of newborns of mothers in the highest folate quartile. CONCLUSIONS: Our findings suggest that fetal TL exhibits developmental plasticity, and provide evidence that maternal nutrition may exert a 'programming' effect on this system.

Maternal cortisol over the course of pregnancy and subsequent child amygdala and hippocampus volumes and affective problems.

Stress-related variation in the intrauterine milieu may impact brain development and emergent function, with long-term implications in terms of susceptibility for affective disorders. Studies in animals suggest limbic regions in the developing brain are particularly sensitive to exposure to the stress hormone cortisol. However, the nature, magnitude, and time course of these effects have not yet been adequately characterized in humans. A prospective, longitudinal study was conducted in 65 normal, healthy mother-child dyads to examine the association of maternal cortisol in early, mid-, and late gestation with subsequent measures at approximately 7 y age of child amygdala and hippocampus volume and affective problems. After accounting for the effects of potential confounding pre- and postnatal factors, higher maternal cortisol levels in earlier but not later gestation was associated with a larger right amygdala volume in girls (a 1 SD increase in cortisol was associated with a 6.4% increase in right amygdala volume), but not in boys. Moreover, higher maternal cortisol levels in early gestation was associated with more affective problems in girls, and this association was mediated, in part, by amygdala volume. No association between maternal cortisol in pregnancy and child hippocampus volume was observed in either sex. The current findings represent, to the best of our knowledge, the first report linking maternal stress hormone levels in human pregnancy with subsequent child amygdala volume and affect. The results underscore the importance of the intrauterine environment and suggest the origins of neuropsychiatric disorders may have their foundations early in life.

A semiparametric Bayesian model for detecting synchrony among multiple neurons.

We propose a scalable semiparametric Bayesian model to capture dependencies among multiple neurons by detecting their cofiring (possibly with some lag time) patterns over time. After discretizing time so there is at most one spike at each interval, the resulting sequence of 1s (spike) and 0s (silence) for each neuron is modeled using the logistic function of a continuous latent variable with a gaussian process prior. For multiple neurons, the corresponding marginal distributions are coupled to their joint probability distribution using a parametric copula model. The advantages of our approach are as follows. The nonparametric component (i.e., the gaussian process model) provides a flexible framework for modeling the underlying firing rates, and the parametric component (i.e., the copula model) allows us to make inferences regarding both contemporaneous and lagged relationships among neurons. Using the copula model, we construct multivariate probabilistic models by separating the modeling of univariate marginal distributions from the modeling of a dependence structure among variables. Our method is easy to implement using a computationally efficient sampling algorithm that can be easily extended to high-dimensional problems. Using simulated data, we show that our approach could correctly capture temporal dependencies in firing rates and identify synchronous neurons. We also apply our model to spike train data obtained from prefrontal cortical areas.

Neonatal Outcomes after Maternal Biomarker-Guided Preterm Birth Intervention: The AVERT PRETERM Trial.

The AVERT PRETERM trial (NCT03151330) evaluated whether screening clinically low-risk pregnancies with a validated maternal blood biomarker test for spontaneous preterm birth (sPTB) risk, followed by preventive treatments for those screening positive, would improve neonatal outcomes compared to a clinically low-risk historical population that had received the usual care. Prospective arm participants with singleton non-anomalous pregnancies and no PTB history were tested for sPTB risk at 191/7-206/7 weeks' gestation and followed up with after neonatal discharge. Screen-positive individuals (≥16% sPTB risk) were offered vaginal progesterone (200 mg) and aspirin (81 mg) daily, with twice-weekly nurse phone calls. Co-primary outcomes were neonatal morbidity and mortality, measured using a validated composite index (NMI), and neonatal hospital length of stay (NNLOS). Endpoints were assessed using survival analysis and logistic regression in a modified intent-to-treat population comprising screen-negative individuals and screen-positive individuals accepting treatment. Of 1460 eligible participants, 34.7% screened positive; of these, 56.4% accepted interventions and 43.6% declined. Compared to historical controls, prospective arm neonates comprising mothers accepting treatment had lower NMI scores (odds ratio 0.81, 95% CI, 0.67-0.98, p = 0.03) and an 18% reduction in severe morbidity. NNLOS was shorter (hazard ratio 0.73, 95% CI, 0.58-0.92, p = 0.01), with a 21% mean stay decrease among neonates having the longest stays. Sensitivity analyses in the entire intent-to-treat population supported these findings. These results suggest that biomarker sPTB risk stratification and preventive interventions can ameliorate PTB complications in singleton, often nulliparous, pregnancies historically deemed low risk.

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.

Maternal psychosocial stress during pregnancy is associated with newborn leukocyte telomere length.

OBJECTIVE: In adults, one of the major determinants of leukocyte telomere length (LTL), a predictor of age-related diseases and mortality, is cumulative psychosocial stress exposure. More recently we reported that exposure to maternal psychosocial stress during intrauterine life is associated with LTL in young adulthood. The objective of the present study was to determine how early in life this effect of stress on LTL is apparent by quantifying the association of maternal psychosocial stress during pregnancy with newborn telomere length. STUDY DESIGN: In a prospective study of N = 27 mother-newborn dyads maternal pregnancy-specific stress was assessed in early gestation and cord blood peripheral blood mononuclear cells were subsequently collected and analyzed for LTL measurement. RESULTS: After accounting for the effects of potential determinants of newborn LTL (gestational age at birth, weight, sex, and exposure to antepartum obstetric complications), there was a significant, independent, linear effect of pregnancy-specific stress on newborn LTL that accounted for 25% of the variance in adjusted LTL (ß = -0.099; P = .04). CONCLUSION: Our finding provides the first preliminary evidence in human beings that maternal psychological stress during pregnancy may exert a "programming" effect on the developing telomere biology system that is already apparent at birth, as reflected by the setting of newborn LTL.

Bayesian nonparametric variable selection as an exploratory tool for discovering differentially expressed genes.

High-throughput scientific studies involving no clear a priori hypothesis are common. For example, a large-scale genomic study of a disease may examine thousands of genes without hypothesizing that any specific gene is responsible for the disease. In these studies, the objective is to explore a large number of possible factors (e.g., genes) in order to identify a small number that will be considered in follow-up studies that tend to be more thorough and on smaller scales. A simple, hierarchical, linear regression model with random coefficients is assumed for case-control data that correspond to each gene. The specific model used will be seen to be related to a standard Bayesian variable selection model. Relatively large regression coefficients correspond to potential differences in responses for cases versus controls and thus to genes that might 'matter'. For large-scale studies, and using a Dirichlet process mixture model for the regression coefficients, we are able to find clusters of regression effects of genes with increasing potential effect or 'relevance', in relation to the outcome of interest. One cluster will always correspond to genes whose coefficients are in a neighborhood that is relatively close to zero and will be deemed least relevant. Other clusters will correspond to increasing magnitudes of the random/latent regression coefficients. Using simulated data, we demonstrate that our approach could be quite effective in finding relevant genes compared with several alternative methods. We apply our model to two large-scale studies. The first study involves transcriptome analysis of infection by human cytomegalovirus. The second study's objective is to identify differentially expressed genes between two types of leukemia.

Improving clinical efficiency in screening for cognitive impairment due to Alzheimer's.

INTRODUCTION: To reduce demands on expert time and improve clinical efficiency, we developed a framework to evaluate whether inexpensive, accessible data could accurately classify Alzheimer's disease (AD) clinical diagnosis and predict the likelihood of progression. METHODS: We stratified relevant data into three tiers: obtainable at primary care (low-cost), mostly available at specialty visits (medium-cost), and research-only (high-cost). We trained several machine learning models, including a hierarchical model, an ensemble model, and a clustering model, to distinguish between diagnoses of cognitively unimpaired, mild cognitive impairment, and dementia due to AD. RESULTS: All models showed viable classification, but the hierarchical and ensemble models outperformed the conventional model. Classifier "error" was predictive of progression rates, and cluster membership identified subgroups with high and low risk of progression within 1.5 to 3 years. DISCUSSION: Accessible, inexpensive clinical data can be used to guide AD diagnosis and are predictive of current and future disease states. HIGHLIGHTS: Classification performance using cost-effective features was accurate and robustHierarchical classification outperformed conventional multinomial classificationClassification labels indicated significant changes in conversion risk at follow-upA clustering-classification method identified subgroups at high risk of decline.

Fully Bayesian Autoencoders with Latent Sparse Gaussian Processes.

We present a fully Bayesian autoencoder model that treats both local latent variables and global decoder parameters in a Bayesian fashion. This approach allows for flexible priors and posterior approximations while keeping the inference costs low. To achieve this, we introduce an amortized MCMC approach by utilizing an implicit stochastic network to learn sampling from the posterior over local latent variables. Furthermore, we extend the model by incorporating a Sparse Gaussian Process prior over the latent space, allowing for a fully Bayesian treatment of inducing points and kernel hyperparameters and leading to improved scalability. Additionally, we enable Deep Gaussian Process priors on the latent space and the handling of missing data. We evaluate our model on a range of experiments focusing on dynamic representation learning and generative modeling, demonstrating the strong performance of our approach in comparison to existing methods that combine Gaussian Processes and autoencoders.

tauFisher accurately predicts circadian time from a single sample of bulk and single-cell 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 tau-Fisher, 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 out-standing performance in both bulk and single-cell transcriptomic data collected from multiple tissue types and experimental settings. Application of tauFisher at a cell-type level in a single-cell RNA-seq 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 RNA-seq data analysis, tauFisher is a promising tool that facilitates circadian medicine and research.