The community persistence threshold was kept at 5 to remove the unstable clusters. This resulted in a hierarchical network with communities as nodes and their hierarchical relationships as edges. We identified significant communities to study network connectedness and the number of hierarchical modules of the overall network. For functional enrichment, the communities were based on the significance of overlap with gene sets from curated gene ontology databases. The Louvain method based on network modularity (Q) is defined as [1]}:
\({\displaystyle Q={\frac{1}{2m}}\sum \limits _{ij}{\bigg [}A_{ij}-{\frac{k_{i}k_{j}}{2m}}{\bigg ]}\delta (c_{i},c_{j})}\) 
