Identification of interactions between molecular features (e.g. mutation, gene expression change) and gross phenotypes in diseases and other biological processes is one of the important challenges in genomic research. Popular approaches such as GSEA are limited to hypothesis tests of bivariate association. However, a specific phenotype is often dependent upon multiple molecular features. It is thus worth considering all possible interactions jointly for a more precise and realistic representation of the cellular network. In this article, a semiparametric copula model is developed to jointly model genotypes, pathways and phenotypes to accomplish this object. A two-step procedure for reconstruction of the network is described. Simulation studies indicate that the method is effective and accurate for the network reconstruction. Application using NCI60 cancer cell line data identifies several subsets of molecular features that jointly perform as the predictors of clinical phenotypes. The copula model is expected to have a broad impact on biomedical research, ranging from cancer treatment to disease prevention.