Learning large-scale gene regulatory networks from single cell transcriptomic data using multivariate information theory