We further identified two hub genes (S100A8 and S100A9) that were significantly associated with CAVD by combining WGCNA module genes and DEGs in the public database. Here, S100A9 is linked to congenital bilateral aplasia of vas deferens from CFTR mutation.