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 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.