METHODS: To understand the contribution of the X chromosome in NPC susceptibility, we conducted an X chromosome-wide association analysis on 1615 NPC patients and 1025 healthy controls of Guangdong Chinese, followed by two validation analyses in Taiwan Chinese (n = 562) and Malaysian Chinese (n = 716).

RESULTS: Firstly, the proportion of variance of X-linked loci over phenotypic variance was estimated in the discovery samples, which revealed that the phenotypic variance explained by X chromosome polymorphisms was estimated to be 12.63% (non-dosage compensation model) in males, as compared with 0.0001% in females. This suggested that the contribution of X chromosome to the genetic variance of NPC should not be neglected. Secondly, association analysis revealed that rs5927056 in DMD gene achieved X chromosome-wide association significance in the discovery sample (OR = 0.81, 95% CI 0.73-0.89, P = 1.49 × 10-5). Combined analysis revealed rs5927056 for DMD gene with suggestive significance (P = 9.44 × 10-5). Moreover, the female-specific association of rs5933886 in ARHGAP6 gene (OR = 0.62, 95%CI: 0.47-0.81, P = 4.37 × 10-4) was successfully replicated in Taiwan Chinese (P = 1.64 × 10-2). rs5933886 also showed nominally significant gender × SNP interaction in both Guangdong (P = 6.25 × 10-4) and Taiwan datasets (P = 2.99 × 10-2).

METHODS: We conducted a meta-analysis of four NPC GWAS among Chinese individuals (2,152 cases; 3,740 controls). Forty-three noteworthy findings outside the MHC region were identified and targeted for replication in a pooled analysis of four independent case-control studies across three regions in Asia (4,716 cases; 5,379 controls). A meta-analysis that combined results from the initial GWA and replication studies was performed.

RESULTS: In the combined meta-analysis, rs31489, located within the CLPTM1L/TERT region on chromosome 5p15.33, was strongly associated with NPC (OR = 0.81; P value 6.3 × 10(-13)). Our results also provide support for associations reported from published NPC GWAS-rs6774494 (P = 1.5 × 10(-12); located in the MECOM gene region), rs9510787 (P = 5.0 × 10(-10); located in the TNFRSF19 gene region), and rs1412829/rs4977756/rs1063192 (P = 2.8 × 10(-8), P = 7.0 × 10(-7), and P = 8.4 × 10(-7), respectively; located in the CDKN2A/B gene region).

CONCLUSIONS: We have identified a novel association between genetic variation in the CLPTM1L/TERT region and NPC. Supporting our finding, rs31489 and other SNPs in this region have been reported to be associated with multiple cancer sites, candidate-based studies have reported associations between polymorphisms in this region and NPC, the TERT gene has been shown to be important for telomere maintenance and has been reported to be overexpressed in NPC, and an EBV protein expressed in NPC (LMP1) has been reported to modulate TERT expression/telomerase activity.

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: We did a genome-wide association study of NKTCL in multiple populations from east Asia. We recruited a discovery cohort of 700 cases with NKTCL and 7752 controls without NKTCL of Han Chinese ancestry from 19 centres in southern, central, and northern regions of China, and four independent replication samples including 717 cases and 12 650 controls. Three of these independent samples (451 cases and 5301 controls) were from eight centres in the same regions of southern, central, and northern China, and the fourth (266 cases and 7349 controls) was from 11 centres in Hong Kong, Taiwan, Singapore, and South Korea. All cases had primary NKTCL that was confirmed histopathologically, and matching with controls was based on geographical region and self-reported ancestry. Logistic regression analysis was done independently by geographical regions, followed by fixed-effect meta-analyses, to identify susceptibility loci. Bioinformatic approaches, including expression quantitative trait loci, binding motif and transcriptome analyses, and biological experiments were done to fine-map and explore the functional relevance of genome-wide association loci to the development of NKTCL.

FINDINGS: Genetic data were gathered between Jan 1, 2008, and Jan 23, 2019. Meta-analysis of all samples (a total of 1417 cases and 20 402 controls) identified two novel loci significantly associated with NKTCL: IL18RAP on 2q12.1 (rs13015714; p=2·83 × 10-16; odds ratio 1·39 [95% CI 1·28-1·50]) and HLA-DRB1 on 6p21.3 (rs9271588; 9·35 × 10-26 1·53 [1·41-1·65]). Fine-mapping and experimental analyses showed that rs1420106 at the promoter of IL18RAP was highly correlated with rs13015714, and the rs1420106-A risk variant had an upregulatory effect on IL18RAP expression. Cell growth assays in two NKTCL cell lines (YT and SNK-6 cells) showed that knockdown of IL18RAP inhibited cell proliferation by cell cycle arrest in NKTCL cells. Haplotype association analysis showed that haplotype 47F-67I was associated with reduced risk of NKTCL, whereas 47Y-67L was associated with increased risk of NKTCL. These two positions are component parts of the peptide-binding pocket 7 (P7) of the HLA-DR heterodimer, suggesting that these alterations might account for the association at HLA-DRB1, independent of the previously reported HLA-DPB1 variants.

INTERPRETATION: Our findings provide new insights into the development of NKTCL by showing the importance of inflammation and immune regulation through the IL18-IL18RAP axis and antigen presentation involving HLA-DRB1, which might help to identify potential therapeutic targets. Taken in combination with additional genetic and other risk factors, our results could potentially be used to stratify people at high risk of NKTCL for targeted prevention.