Systems biology provides a comprehensive approach to biomarker discovery and biological hypothesis building. Indeed, it allows to jointly consider the statistical association between gene variation and a phenotype, and the biological context of each gene, represented as a network. In this work, we study six network methods which identify subnetworks with high association scores to a phenotype. Specifically, we examine their utility to discover new biomarkers for breast cancer susceptibility by interrogating a genome-wide association study (GWAS) focused on French women with a family history of breast cancer and tested negative for pathogenic variants in BRCA1 and BRCA2. We perform an in-depth benchmarking of the methods with regards to size of the solution subnetwork, their utility as biomarkers, and the stability and the runtime of the methods. By trading statistical stringency for biological meaningfulness, most network methods give more compelling results than standard SNP-and gene-level analyses, recovering causal subnetworks tightly related to cancer susceptibility. For instance, we show a general alteration of the neighborhood of COPS5, a gene related to multiple hallmarks of cancer. Importantly, we find a significantly large overlap between the genes in the solution networks and the genes significantly associated in the largest GWAS on susceptibility to breast cancer. Yet, network methods are notably unstable, producing di erent results when the input data changes slightly. To account for that, we produce a stable consensus subnetwork, formed by the most consistently selected genes. The stable consensus is composed of 68 genes, enriched in known breast cancer susceptibility genes (BLM, CASP8, CASP10, DNAJC1, FGFR2, MRPS30, and SLC4A7, Fisher's exact test P-value = 3 ◊ 10 ≠4 ) and occupying more central positions in the network than average. The network seems organized around CUL3, encoding an ubiquitin ligase related protein that regulates the protein levels of 1 several genes involved in cancer progression. In conclusion, this article shows the pertinence of network-based analyses to tackle known issues with GWAS, namely lack of statistical power and of interpretable solutions. Project-agnostic implementations of each of the network methods are available at https://github.com/hclimente/gwas-tools to facilitate their application to other GWAS datasets.
NetworksGene network The statistical frameworks of the di erent network methods are compatible with any type of network (protein interactions, gene coexpression, regulatory, etc.). Yet, we used proteinprotein interaction networks (PPIN) for all of them except SConES, as they are interpretable, well characterized, and they were designed to run e ciently on networks of their size. We built our PPIN from both binary and co-complex interactions stored in the HINT database (release April 2019) [17]. Unless otherwise specified, we used only interactions coming from high-throughput experiments, leaving out targeted studies that might bias the topology of the networ...