Summary: We consider penalized linear regression, especially for “large \(p\), small \(n\)” problems, for which the relationships among predictors are described a priori by a network. A class of motivating examples includes modeling a phenotype through gene expression profiles while accounting for coordinated functioning of genes in the form of biological pathways or networks. To incorporate the prior knowledge of the similar effect sizes of neighboring predictors in a network, we propose a grouped penalty based on the \(L_{\gamma}\)-norm that smoothes the regression coefficients of the predictors over the network. The main feature of the proposed method is its ability to automatically realize grouped variable selection and exploit grouping effects. We also discuss effects of the choices of the \(\gamma \) and some weights inside the \(L_{\gamma}\) -norm. Simulation studies demonstrate the superior finite-sample performance of the proposed method as compared to Lasso, elastic net, and a recently proposed network-based method. The new method performs best in variable selection across all simulation set-ups considered. For illustration, the method is applied to a microarray data set to predict survival times for some glioblastoma patients using a gene expression data set and a gene network compiled from some Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways.