RESULTS: We identified a total of 644,225 SNPs in 131 neuropeptide genes in 6 worldwide population groups from a public database. Of these, 5163 SNPs that had ΔDAF |(African - non-African)| ≥ 0.20 were identified and fully annotated. A total of 20 outlier SNPs that included 19 missense SNPs with a moderate impact and one stop lost SNP with high impact, were identified in 16 neuropeptide genes. Our results indicate that an overall strong population differentiation was observed in the non-African populations that had a higher derived allele frequency for 15/20 of those SNPs. Highly differentiated SNPs in four genes were particularly striking: NPPA (rs5065) with high impact stop lost variant; CHGB (rs6085324, rs236150, rs236152, rs742710 and rs742711) with multiple moderate impact missense variants; IGF2 (rs10770125) and INS (rs3842753) with moderate impact missense variants that are in linkage disequilibrium. Phenotype and disease associations of these differentiated SNPs indicated their association with hypertension and diabetes and highlighted the pleiotropic effects of these neuropeptides and their role in maintaining physiological homeostasis in humans.

METHODS: We did a genome-wide association study of 189 patients with extranodal NKTCL, nasal type (WHO classification criteria; cases) and 957 controls from Guangdong province, southern China. We validated our findings in four independent case-control series, including 75 cases from Guangdong province and 296 controls from Hong Kong, 65 cases and 983 controls from Guangdong province, 125 cases and 1110 controls from Beijing (northern China), and 60 cases and 2476 controls from Singapore. We used imputation and conditional logistic regression analyses to fine-map the associations. We also did a meta-analysis of the replication series and of the entire dataset.

FINDINGS: Associations exceeding the genome-wide significance threshold (p<5 × 10(-8)) were seen at 51 single-nucleotide polymorphisms (SNPs) mapping to the class II MHC region on chromosome 6, with rs9277378 (located in HLA-DPB1) having the strongest association with NKTCL susceptibility (p=4·21 × 10(-19), odds ratio [OR] 1·84 [95% CI 1·61-2·11] in meta-analysis of entire dataset). Imputation-based fine-mapping across the class II MHC region suggests that four aminoacid residues (Gly84-Gly85-Pro86-Met87) in near-complete linkage disequilibrium at the edge of the peptide-binding groove of HLA-DPB1 could account for most of the association between the rs9277378*A risk allele and NKTCL susceptibility (OR 2·38, p value for haplotype 2·32 × 10(-14)). This association is distinct from MHC associations with Epstein-Barr virus infection.

INTERPRETATION: To our knowledge, this is the first time that a genetic variant conferring an NKTCL risk is noted at genome-wide significance. This finding underlines the importance of HLA-DP antigen presentation in the pathogenesis of NKTCL.

METHODS: Genomic DNAs were extracted from the blood samples followed by whole-genome sequencing. The reads were aligned to the reference human genome hg19 and variants in the CYP2D6 gene were analyzed. CYP2D6*5 and duplication of CYP2D6 were analyzed using previously established methods.

RESULTS: A total of 72 single nucleotide polymorphisms were identified. CYP2D6*1, *2, *4, *5, *10,*41, and duplication of the gene were found in the Orang Asli, whereby CYP2D6*2 and *41 alleles are reported for the first time in the Malaysian population.