AbstractInflammatory bowel diseases (IBDs) are complex medical conditions in which the gut microbiota is attacked by the immune system of genetically predisposed subjects when exposed to yet unclear environmental factors. The complexity of this class of diseases makes them suitable to be represented and studied with network science. In this paper, the metagenomic data of control, Crohn’s disease, and ulcerative colitis subjects’ gut microbiota were investigated by representing this data as correlation networks and co-expression networks. We obtained correlation networks by calculating Pearson’s correlation between gene expression across subjects. A percolation-based procedure was used to threshold and binarize the adjacency matrices. In contrast, co-expression networks involved the construction of the bipartite subjects-genes networks and the monopartite genes-genes projection after binarization of the biadjacency matrix. Centrality measures and community detection were used on the so-built networks to mine data complexity and highlight possible biomarkers of the diseases. The main results were about the modules of Bacteroides, which were connected in the control subjects’ correlation network, Faecalibacterium prausnitzii, where co-enzyme A became central in IBD correlation networks and Escherichia coli, whose module has different patterns of integration within the whole network in the different diagnoses.