Previous seed- and atlas-based structural covariance/connectivity analyses have demonstrated that patients with schizophrenia is accompanied by aberrant structural connection and abnormal topological organization. However, it remains unclear whether this disruption is present in unbiased whole-brain voxel-wise structural covariance networks (SCNs) and whether brain genetic expression variations are linked with network alterations. In this study, ninety-five patients with schizophrenia and 95 matched healthy controls were recruited and gray matter volumes were extracted from high-resolution structural magnetic resonance imaging scans. Whole-brain voxel-wise gray matter SCNs were constructed at the group level and were further analyzed by using graph theory method. Nonparametric permutation tests were employed for group comparisons. In addition, regression modes along with random effect analysis were utilized to explore the associations between structural network changes and gene expression from the Allen Human Brain Atlas. Compared with healthy controls, the patients with schizophrenia showed significantly increased structural covariance strength (SCS) in the right orbital part of superior frontal gyrus and bilateral middle frontal gyrus, while decreased SCS in the bilateral superior temporal gyrus and precuneus. The altered SCS showed reproducible correlations with the expression profiles of the gene classes involved in therapeutic targets and neurodevelopment. Overall, our findings not only demonstrate that the topological architecture of whole-brain voxel-wise SCNs is impaired in schizophrenia, but also provide evidence for the possible role of therapeutic targets and neurodevelopment-related genes in gray matter structural brain networks in schizophrenia.