Graphical models are essential to describe networks among different factors.Various methods to construct graphs have been proposed.In practice,it is critical to assess the agreement between networks constructed under different treats or conditions.There is rather limited effort in this area.We propose a Bayesian method to build directed acyclic graphs(DAGs)based on a given order of nodes and simultaneously test the agreement between DAGs.The network construction and differential graphs testing are built upon the concept of variable selection.Simulations demonstrate the applicability of the method and a real data application to an epigenetic data is implemented to illustrate the method.