Genetic interactions or epistasis could make a substantial contribution to variation in human complex traits including longevity. However, detecting epistatic interactions in high dimensional datasets is difficult due to various reasons including multiple testing of correlated tests. We introduce a novel permutation strategy to the case-only analysis of gene by gene interaction using multiple SNPs. The method is applied to genes coding for Forkhead box O transcription factors which recently have been associated with human longevity across different populations hypothesizing that epistatic interaction in the regulation and expression of the FOXO gene family could contribute to the human longevity phenotype. Genotype data was collected from 1088 individuals from the Danish 1905 birth cohort aged over 92/93 years with 12 SNPs in the FOXO1a and 15 SNPs in the FOXO3a genes. Our analysis detected a joint effect between rs9486902 in FOXO3a and rs2701858 in FOXO1a that highly significantly contributes to human longevity (OR=3.23, 95% CI: 2.93-3.53) which is consistent in both males and females. Our results were compared with published studies and importance of our novel method and findings discussed. This article is protected by copyright. All rights reserved.