Graphic Model-based Gene Regulatory Network Reconstruction using RNA Sequencing Count Data

Interactions between genes, such as regulations are best represented by gene regulatory networks (GRN). These are often constructed based on gene expression data. Few methods for the construction of GRN exist for RNA sequencing count data. One of the most used methods for microarray data is based on graphical Gaussian networks. Considering that count data have different distributions, a method assuming RNA sequencing counts distribute Poisson has been proposed recently. Nevertheless, it has been argued that the most likely distribution of RNA sequencing counts is not Poisson due to overdispersion. Therefore, the negative binomial distribution is much more likely. For this distribution, no model-based method for the construction of GRN has been proposed until now. Here, we present a graphical, model-based method for the construction of GRN assuming a negative binomial distribution of the RNA sequencing count data. The R code is available under request.

https://rajpub.com/index.php/jbt/article/view/8298