Prostate cancer represents a significant worldwide public health burden. Epidemiological and genetic epidemiological studies have consistently provided data supporting the existence of inherited prostate cancer susceptibility genes. Segregation analyses of prostate cancer suggest that a multigene model may best explain familial clustering of this disease. Therefore, modeling gene-gene interactions in linkage analysis may improve the power to detect chromosomal regions harboring these disease susceptibility genes. In this study, we systematically screened for prostate cancer linkage by modeling two-locus gene-gene interactions for all possible pairs of loci across the genome in 426 prostate cancer families from Johns Hopkins Hospital, University of Michigan, University of Umeå, and University of Tampere. We found suggestive evidence for an epistatic interaction for six sets of loci (target chromosome-wide/reference marker-specific P< or =0.0001). Evidence for these interactions was found in two independent subsets from within the 426 families. While the validity of these results requires confirmation from independent studies and the identification of the specific genes underlying this linkage evidence, our approach of systematically assessing gene-gene interactions across the entire genome represents a promising alternative approach for gene identification for prostate cancer.