In this study, a multi-omics approach was used to identify 10 CAVD-related biomarkers (COL5A1, COL5A2, THBS2, MFAP5, BTG2, COL1A1, COL1A2, MXRA5, LUM, CD34) and to develop an exploratory risk model. Here, COL1A1 is linked to congenital bilateral aplasia of vas deferens from CFTR mutation.