Bayesian testing of many hypotheses ×\times many genes: A study of sleep apnea

Substantial statistical research has recently been devoted to the analysis of large-scale microarray experiments which provide a measure of the simultaneous expression of thousands of genes in a particular condition. A typical goal is the comparison of gene expression between two conditions (e.g., diseased vs. nondiseased) to detect genes which show differential expression. Classical hypothesis testing procedures have been applied to this problem and more recent work has employed sophisticated models that allow for the sharing of information across genes. However, many recent gene expression studies have an experimental design with several conditions that requires an even more involved hypothesis testing approach. In this paper, we use a hierarchical Bayesian model to address the situation where there are many hypotheses that must be simultaneously tested for each gene. In addition to having many hypotheses within each gene, our analysis also addresses the more typical multiple comparison issue of testing many genes simultaneously. We illustrate our approach with an application to a study of genes involved in obstructive sleep apnea in humans.Comment: Published in at http://dx.doi.org/10.1214/09-AOAS241 the Annals of Applied Statistics (http://www.imstat.org/aoas/) by the Institute of Mathematical Statistics (http://www.imstat.org

A bioinformatic analysis identifies circadian expression of splicing factors and time-dependent alternative splicing events in the HD-MY-Z cell line

The circadian clock regulates key cellular processes and its dysregulation is associated to several pathologies including cancer. Although the transcriptional regulation of gene expression by the clock machinery is well described, the role of the clock in the regulation of post-transcriptional processes, including splicing, remains poorly understood. In the present work, we investigated the putative interplay between the circadian clock and splicing in a cancer context. For this, we applied a computational pipeline to identify oscillating genes and alternatively spliced transcripts in time-course high-throughput data sets from normal cells and tissues, and cancer cell lines. We investigated the temporal phenotype of clock-controlled genes and splicing factors, and evaluated their impact in alternative splice patterns in the Hodgkin Lymphoma cell line HD-MY-Z. Our data points to a connection between clock-controlled genes and splicing factors, which correlates with temporal alternative splicing in several genes in the HD-MY-Z cell line. These include the genes DPYD, SS18, VIPR1 and IRF4, involved in metabolism, cell cycle, apoptosis and proliferation. Our results highlight a role for the clock as a temporal regulator of alternative splicing, which may impact malignancy in this cellular model

A temporal precedence based clustering method for gene expression microarray data

Background: Time-course microarray experiments can produce useful data which can help in understanding the underlying dynamics of the system. Clustering is an important stage in microarray data analysis where the data is grouped together according to certain characteristics. The majority of clustering techniques are based on distance or visual similarity measures which may not be suitable for clustering of temporal microarray data where the sequential nature of time is important. We present a Granger causality based technique to cluster temporal microarray gene expression data, which measures the interdependence between two time-series by statistically testing if one time-series can be used for forecasting the other time-series or not. Results: A gene-association matrix is constructed by testing temporal relationships between pairs of genes using the Granger causality test. The association matrix is further analyzed using a graph-theoretic technique to detect highly connected components representing interesting biological modules. We test our approach on synthesized datasets and real biological datasets obtained for Arabidopsis thaliana. We show the effectiveness of our approach by analyzing the results using the existing biological literature. We also report interesting structural properties of the association network commonly desired in any biological system. Conclusions: Our experiments on synthesized and real microarray datasets show that our approach produces encouraging results. The method is simple in implementation and is statistically traceable at each step. The method can produce sets of functionally related genes which can be further used for reverse-engineering of gene circuits

SensibleSleep: A Bayesian Model for Learning Sleep Patterns from Smartphone Events

We propose a Bayesian model for extracting sleep patterns from smartphone events. Our method is able to identify individuals' daily sleep periods and their evolution over time, and provides an estimation of the probability of sleep and wake transitions. The model is fitted to more than 400 participants from two different datasets, and we verify the results against ground truth from dedicated armband sleep trackers. We show that the model is able to produce reliable sleep estimates with an accuracy of 0.89, both at the individual and at the collective level. Moreover the Bayesian model is able to quantify uncertainty and encode prior knowledge about sleep patterns. Compared with existing smartphone-based systems, our method requires only screen on/off events, and is therefore much less intrusive in terms of privacy and more battery-efficient

Quantifying sleep architecture dynamics and individual differences using big data and Bayesian networks

The pattern of sleep stages across a night (sleep architecture) is influenced by biological, behavioral, and clinical variables. However, traditional measures of sleep architecture such as stage proportions, fail to capture sleep dynamics. Here we quantify the impact of individual differences on the dynamics of sleep architecture and determine which factors or set of factors best predict the next sleep stage from current stage information. We investigated the influence of age, sex, body mass index, time of day, and sleep time on static (e.g. minutes in stage, sleep efficiency) and dynamic measures of sleep architecture (e.g. transition probabilities and stage duration distributions) using a large dataset of 3202 nights from a non-clinical population. Multi-level regressions show that sex effects duration of all Non-Rapid Eye Movement (NREM) stages, and age has a curvilinear relationship for Wake After Sleep Onset (WASO) and slow wave sleep (SWS) minutes. Bayesian network modeling reveals sleep architecture depends on time of day, total sleep time, age and sex, but not BMI. Older adults, and particularly males, have shorter bouts (more fragmentation) of Stage 2, SWS, and they transition less frequently to these stages. Additionally, we showed that the next sleep stage and its duration can be optimally predicted by the prior 2 stages and age. Our results demonstrate the potential benefit of big data and Bayesian network approaches in quantifying static and dynamic architecture of normal sleep

Modeling circadian and sleep-homeostatic effects on short-term interval timing

Short-term interval timing i.e., perception and action relating to durations in the seconds range, has been suggested to display time-of-day as well as wake dependent fluctuations due to circadian and sleep-homeostatic changes to the rate at which an underlying pacemaker emits pulses; pertinent human data being relatively sparse and lacking in consistency however, the phenomenon remains elusive and its mechanism poorly understood. To better characterize the putative circadian and sleep-homeostatic effects on interval timing and to assess the ability of a pacemaker-based mechanism to account for the data, we measured timing performance in eighteen young healthy male subjects across two epochs of sustained wakefulness of 38.67 h each, conducted prior to (under entrained conditions) and following (under free-running conditions) a 28 h sleep-wake schedule, using the methods of duration estimation and duration production on target intervals of 10 and 40 s. Our findings of opposing oscillatory time courses across both epochs of sustained wakefulness that combine with increasing and, respectively, decreasing, saturating exponential change for the tasks of estimation and production are consistent with the hypothesis that a pacemaker emitting pulses at a rate controlled by the circadian oscillator and increasing with time awake determines human short-term interval timing; the duration-specificity of this pattern is interpreted as reflecting challenges to maintaining stable attention to the task that progressively increase with stimulus magnitude and thereby moderate the effects of pacemaker-rate changes on overt behavior

Precision medicine and artificial intelligence : a pilot study on deep learning for hypoglycemic events detection based on ECG

Tracking the fluctuations in blood glucose levels is important for healthy subjects and crucial diabetic patients. Tight glucose monitoring reduces the risk of hypoglycemia, which can result in a series of complications, especially in diabetic patients, such as confusion, irritability, seizure and can even be fatal in specific conditions. Hypoglycemia affects the electrophysiology of the heart. However, due to strong inter-subject heterogeneity, previous studies based on a cohort of subjects failed to deploy electrocardiogram (ECG)-based hypoglycemic detection systems reliably. The current study used personalised medicine approach and Artificial Intelligence (AI) to automatically detect nocturnal hypoglycemia using a few heartbeats of raw ECG signal recorded with non-invasive, wearable devices, in healthy individuals, monitored 24 hours for 14 consecutive days. Additionally, we present a visualisation method enabling clinicians to visualise which part of the ECG signal (e.g., T-wave, ST-interval) is significantly associated with the hypoglycemic event in each subject, overcoming the intelligibility problem of deep-learning methods. These results advance the feasibility of a real-time, non-invasive hypoglycemia alarming system using short excerpts of ECG signal