Bayesian inference of causal gene networks