Multifactor dimensionality reduction for interactions among genes revealed a 4-locus model of BMPR-2 rs6717924G/A; ALK-1 rs11169953C/T-rs706816T/C and 5-HTT rs6354A/C as the best disease predicting (P≤0.001); whereas HapEvolution analysis confirmed the alleles rs6717924A, rs4303700A and rs6354C as the best interacting (P≤0.01).