Rxivist logo

Causal Network Inference from Gene Transcriptional Time Series Response to Glucocorticoids

By Jonathan Lu, Bianca Dumitrascu, Ian C McDowell, Brian Jo, Alejandro Barrera, Linda K Hong, Sarah M. Leichter, Timothy E. Reddy, Barbara E Engelhardt

Posted 24 Mar 2019
bioRxiv DOI: 10.1101/587170

Gene regulatory network inference is essential to uncover complex relationships among gene pathways and inform downstream experiments, ultimately paving the way for regulatory network re-engineering. Network inference from transcriptional time series data requires accurate, interpretable, and efficient determination of causal relationships among thousands of genes. Here, we develop Bootstrap Elastic net regression from Time Series (BETS), a statistical framework based on Granger causality for the recovery of a directed gene network from transcriptional time series data. BETS uses elastic net regression and stability selection from bootstrapped samples to infer causal relationships among genes. BETS is highly parallelized, enabling efficient analysis of large transcriptional data sets. We show competitive accuracy on a community benchmark, the DREAM4 100-gene network inference challenge, where BETS is one of the fastest among methods of similar performance but additionally infers whether the causal effects are activating or inhibitory. We apply BETS to transcriptional time series data of 2,768 differentially-expressed genes from A549 cells exposed to glucocorticoids over a period of 12 hours. We identify a network of 2,768 genes and 31,945 directed edges (FDR ≥ 0.2). We validate inferred causal network edges using two external data sources: overexpression experiments on the same glucocorticoid system, and genetic variants associated with inferred edges in primary lung tissue in the Genotype-Tissue Expression (GTEx) v6 project. BETS is freely available as an open source software package at https://github.com/lujonathanh/BETS.

Download data

  • Downloaded 1,106 times
  • Download rankings, all-time:
    • Site-wide: 9,762 out of 89,516
    • In bioinformatics: 1,613 out of 8,443
  • Year to date:
    • Site-wide: 5,549 out of 89,516
  • Since beginning of last month:
    • Site-wide: 4,200 out of 89,516

Altmetric data

Downloads over time

Distribution of downloads per paper, site-wide


Sign up for the Rxivist weekly newsletter! (Click here for more details.)