- Open Access
X-chromosome association study reveals genetic susceptibility loci of nasopharyngeal carcinoma
- Xiao-Yu Zuo†1,
- Qi-Sheng Feng†1,
- Jian Sun1,
- Pan-Pan Wei1,
- Yoon-Ming Chin2,
- Yun-Miao Guo1,
- Yun-Fei Xia1,
- Bo Li3, 4,
- Xiao-Jun Xia1,
- Wei-Hua Jia1,
- Jian-Jun Liu5,
- Alan Soo-Beng Khoo6,
- Taisei Mushiroda7,
- Ching-Ching Ng†2Email author,
- Wen-Hui Su†8, 9Email author,
- Yi-Xin Zeng†1Email author and
- Jin-Xin Bei†1, 10Email author
© The Author(s). 2019
- Received: 9 January 2019
- Accepted: 27 February 2019
- Published: 25 March 2019
The male predominance in the incidence of nasopharyngeal carcinoma (NPC) suggests the contribution of the X chromosome to the susceptibility of NPC. However, no X-linked susceptibility loci have been examined by genome-wide association studies (GWASs) for NPC by far.
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).
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).
Our finding reveals new susceptibility loci at the X chromosome conferring risk of NPC and supports the value of including the X chromosome in large-scale association studies.
- Nasopharyngeal carcinoma
- Genetic susceptibility
- X chromosome
- Association study
- Male predominance
Nasopharyngeal carcinoma (NPC) is a human squamous-cell carcinoma that arises in the epithelium of the nasopharynx. The most remarkable feature of NPC is its restricted geographical distribution of incidence rate, which is about 30 per 100,000 in its prevalent regions including Southern China but rare in the western countries . The etiology of NPC is a multifactorial process, involving genetic, viral, and environmental factors . Large-scale case-control association studies and familial linkage studies have confirmed the genetic contribution to NPC predisposition, by revealing multiple susceptibility loci of NPC, such as HLA genes [3–6], TNFRSF19 , MECOM , GABBR1 , XRCC3 , ITGA9 , TERT-CLPTM1L [8–10], and CIITA . However, these explain only a small fraction of the NPC heritability [11, 12].
The predominant incidence of males is another feature of NPC, with a male/female ratio of two to threefolds in diverse populations [13, 14]. The age-standardized rates of NPC for males:females were 21.8:8.1 per 100,000 person-years in Hong Kong  and 28.93:11.26 in Sihui, Guangdong, where both cities are NPC prevalent regions in Southern China . Unequal exposures to environmental risk factors between males and females, such as smoking and diet habit, may partially contribute to the gender difference in the incidence of NPC . It has been suggested that some tumor suppressor genes may contribute to gender bias in cancer development by escaping from the inactivation of the X chromosome in females . Moreover, the involvement of the X chromosome has been implicated in the predisposition to NPC [18, 19]; however, the exact genetic components in the X chromosome have remained largely unexploited in NPC. By contrast, recent X chromosome association studies have identified novel risk loci at the sexual chromosome for follicular lymphoma , schizophrenia , systemic lupus erythematosus , etc., which partially explain the missing heritability of those complex diseases.
To explore the genetic contribution of the X chromosome to the gender difference of NPC incidence, we performed an X chromosome-wide association study, with the discovery data of 1615 cases and 1025 controls in Southern China, followed by two validation analyses involving a total of 543 cases and 735 controls from Taiwan Chinese and Malaysian Chinese, respectively.
The discovery sample included 1615 NPC cases and 1025 controls derived from a previous genome-wide association study (GWAS) in Southern Chinese . All cases were recruited through Sun Yat-sen University Cancer Center (SYSUCC) in Southern China during October 2005 and October 2007, whereas the controls were recruited from several physical examination centers in local communities in Guangdong during the same period. As for the validation stage, two independent case-control cohorts of Chinese descendants were included, with 266 NPC cases and 450 controls recruited from different states of Malaysia [4, 23], and 277 NPC cases and 285 controls recruited from the northern parts of Taiwan, respectively . The diagnosis of NPC was confirmed according to the World Health Organization (WHO) classification at each study site. The study was approved by the Institutional Review Board at SYSUCC. Informed consent was obtained from all participants. The information for all subjects is summarized in Additional file 1: Table S1.
Genotyping and quality control
Individual genotypes had been determined as described previously . In brief, the genomic DNA was extracted from peripheral blood sample using commercial (Qiagen) DNA extraction kit (Southern China and Taiwan studies) or conventional methods (Malaysia study), and the genotyping was conducted by using Illumina BeadChip arrays, according to the manufacturer’s protocols (Illumina, Inc., San Diego, CA, USA). For the discovery stage, quality checks were applied for each sample as described in the previous study , using autosomal SNPs and removing those with genotyping rate < 95%, excessive observed level of heterozygosity (departure from 3 standard deviation), cryptic relatedness, error or uncertainty in gender estimation, and population outliers through principal component analysis (PCA). A total of 1590 cases and 994 controls were retained for subsequent analyses. We extracted 18,133 genotyped SNPs at the X chromosome from the dataset and conducted further quality control filtering for the SNPs as suggested previously [24, 25]. As heterozygote of a SNP in a male should be deemed to be a genotyping/calling error, a fact of being haploidy for the X chromosome in a male, thus, we assigned a missing value for such a call in a male sample. The SNPs were removed if they met the following criteria : with genotyping rate < 92% for all samples or call rate difference > 3% or P < 1 × 10−5 between gender , with the difference of heterozygosity level > 5% or P < 1 × 10−3 between cases and controls in males , with minor allele frequency (MAF) < 1% for both males and females , deviation from Hardy-Weinberg equilibrium test P < 1 × 10−6 in females, or  with heterozygote genotypes found in 5% of males (due to genotyping error). For SNPs in pseudo-autosomal regions (PAR), we used the same filtering criteria as for the autosomal SNPs in our previous study . Finally, 6536 SNPs were removed, and 11,597 X chromosomal SNPs were used for subsequent analyses.
For the validation stage, genotypes were retrieved from the two GWAS studies [3, 4, 23]. Genotypes for unknown SNPs in the GWAS datasets were imputed by using non-PAR SNPs and IMPUTE2 program [26, 27], with 1000 Genomes Phase I integrated variant set (March 2012, Build 37) as the reference. Imputed SNPs with imputation score (INFO in the *.impute2_info output file) < 0.5 were considered low confidence and removed. Then, Gtool program (v0.7.5) was used to convert data from IMPUTE2 to PLINK. Genotypes with the threshold for calling genotypes (−threshold) < 0.9 were considered as low confidence and set as missing. Imputed SNPs were subjected to the same quality control procedure as that of the genotyped SNPs. As a result, 99,710 imputed loci were retained for downstream analysis. For the Malaysian Chinese dataset, 1000 Genomes Phase I interim set (June 2011, Build 37) was used as the reference. The genotypes of validating SNPs were retrieved from the Taiwan Chinese dataset as previously described .
Statistical and bioinformatics analysis
Association tests were performed by using PLINK implemented with logistic regression analysis under the generalized linear model [28, 29]. Genotypes of the X chromosome loci were coded as [0, 2] for males and [0, 1, 2] for females, accounting for the random inactivation on females and assuming activation of one of the two X chromosomes in females. Given that a proportion of SNPs can escape from random inactivation, genotype coding of [0, 1] for males and [0, 1, 2] for females was also tested in combined samples to account for the escaping from random inactivation. Gender was treated as a covariate, and the first 10 principal components (PCs) were also adjusted in the association test in the discovery sample. To investigate whether there are sex-specific effects underlying X chromosome loci, we performed association tests for males and females separately and also tests accounting for gender × SNP interaction. For combined analyses of discovery and replication datasets, meta-analyses were conducted in R (version 3.2.3) using the “metafor” package (version 1.9–8). The heterogeneity among datasets was first evaluated by I2 and Q test. I2 > 50% or P < 0.05 were considered as heterogeneous, and the random-effect model was applied; otherwise, the fixed-effect model would be used.
We used linear mixed model approach implemented in GCTA (version 1.24.7) [30, 31] to estimate the contribution of the X chromosome on the proportion of the NPC phenotypic variance in the discovery dataset. The equal variance in both sexes (EV), full-dosage compensation (FDC), and non-dosage compensation (NDC) models were applied . Only genotyped SNPs were used in this study to reduce the potential impact of subtle imputation uncertainty. For a comparison, we also estimated the variance explained by each autosome by using common SNPs passing QC and with MAF > 5%.
We used the simpleM method to address the effective number of independent association tests  since the Bonferroni correction for multiple testing was too conservative for the association test among loci with considerable linkage disequilibrium (LD). The simpleM method adopts a principal component analysis on the SNP correlation matrix to pick up the least “effective number” of tests accounting for ≥ 99.5% variance. In this study, the effective number was estimated as 2308, irrespective of genotyped or imputed loci. Therefore, we adopted the X chromosome-wide significance level as 0.05/2308 = 2.17 × 10−5. In the discovery stage, we chose a relatively relaxed significance threshold (P < 1 × 10−3) to select candidate SNPs for the follow-up replications, allowing inclusion of more SNPs and higher sensitivity of reproducible association signals. We used the R package “twoStageGwasPower” to calculate statistical power in this study, following the methods described previously . The key raw data of this study have been uploaded onto the Research Data Deposit (RDD; http://www.researchdata.org.cn/; Number: RDDB2019000532).
Furthermore, noncoding susceptible SNPs identified in this study were subjected to HaploReg (V4.1; https://pubs.broadinstitute.org/mammals/haploreg/haploreg.php)  to annotate their functional and regulatory potentials on chromatin state and protein-binding annotations from the Roadmap Epigenomics and ENCODE projects, sequence conservation across mammals, and the effect of SNPs on regulatory motifs. Expression quantitative trait locus (eQTL) effects of the SNPs were estimated by using GTEx portal, where samples were collected from 53 non-diseased tissue sites across nearly 1000 individuals and tissue-specific gene expression and regulation with particularly the correlation between SNP and the expression of nearby genes was analyzed and archived (V7; https://gtexportal.org/home/).
Genetic variance in X chromosome contributes to the risk of NPC
X chromosome-wide association analysis
Representative associations of SNPs on X chromosome for combined samples
Replication 1 (Malaysia)
Replication 2 (Taiwan)
1.49 × 10−5
9.47 × 10−1
3.33 × 10−1
9.44 × 10−5
299 kb 3′ of LOC101928201
2.86 × 10−5
6.58 × 10−1
3.44 × 10−1
3.01 × 10−4
7.20 × 10−5
5.58 × 10−1
8.86 × 10−1
3.98 × 10−1
261 kb 5′ of CXorf51A
4.02 × 10−4
1.43 × 10−1
7.04 × 10−1
7.30 × 10−4
7.69 × 10−4
6.44 × 10−1
6.86 × 10−1
1.95 × 10−3
108 kb 3′ of C1GALT1C1
3.86 × 10−4
4.89 × 10−1
6.97 × 10−1
6.17 × 10−1
43 kb 5′ of SMS
2.41 × 10−4
2.07 × 10−1
2.76 × 10−1
9.86 × 10−1
430 kb 5′ of BRDTP1
3.40 × 10−4
1.39 × 10−1
3.74 × 10−1
9.90 × 10−1
Representative gender-specific association on X chromosome in discovery cohort
299 kb 3′ of LOC101928201
2.50 × 10−5
6.05 × 10−1
2.34 × 10−1
4.51 × 10−5
7.02 × 10−1
2.58 × 10−1
427 kb 3′ of ZFX
1.21 × 10−4
2.17 × 10−1
1.20 × 10−2
1.17 × 10−4
4.66 × 10−2
7.62 × 10−1
400 kb 5′ of BRDTP1
1.20 × 10−4
1.77 × 10−1
1.20 × 10−2
47Kb 3′ of MAGEB18
8.65 × 10−1
5.02 × 10−5
2.40 × 10−3
1.40 × 10−1
1.42 × 10−4
1.19 × 10−4
8.01 × 10−1
1.67 × 10−4
7.57 × 10−2
4.10 × 10−1
4.37 × 10−4
6.25 × 10−4
177Kb 5′ of VENTXP1
2.58 × 10−1
2.32 × 10−4
3.21 × 10−3
1.29 × 10−1
1.48 × 10−4
8.35 × 10−5
20 kb 5′ of MAOA
7.88 × 10−3
2.12 × 10−3
1.22 × 10−4
143 kb 5′ of NAP1L3
2.00 × 10−1
3.49 × 10−4
1.86 × 10−4
374 kb 3′ of CT83
1.90 × 10−1
7.25 × 10−4
2.24 × 10−4
20 kb 3′ of MBTPS2
6.23 × 10−1
6.01 × 10−4
4.71 × 10−4
3.99 × 10−1
3.81 × 10−4
4.72 × 10−4
Validation and combined analyses
Functional annotations of the SNPs
In silico analyses were conducted to explore the functional potentials of the susceptibility SNPs on the X chromosomes. HaploReg revealed that the top significant rs5927056, located in the intron of DMD, might alter six regulatory motifs including AP-3, Evi-1_4, Hoxa10, Hoxb13, Hoxd10, and Pou2f2_known11, suggesting its regulatory potentials (Additional file 1: Table S7). HaploReg also revealed regulatory motifs (HMG-IY_1, Ik-3, NFKB_known5, Pou3f1, and STAT_known5) and DNase peaks at rs5933886, which is an intronic SNP of ARHGAP6 (Additional file 1: Table S8). Moreover, the eQTL analysis revealed a significant cis-eQTL effect of rs5933886 in the aorta artery sample from the GTEx database, which were collected from the ascending aorta or other thoracic regions (nonatherosclerotic; P = 5.88 × 10−3, Additional file 1: Table S8). These results imply that rs5933886 may act as a regulatory SNP and predispose NPC by regulating the expression of ARHGAP6.
The human X chromosome consists ~ 155 million base pairs and contains more than 1000 genes that are important for many biological processes [17, 30]. However, none of the previous GWAS of NPC included association analysis of variations at the X chromosome, overlooking potential susceptibility loci that could partially account for the missing heritability of NPC. Here, we provided, to the best of our knowledge, the first large-scale X chromosome-wide association study of NPC in the Chinese population, revealing the contribution of the X chromosome and suggestive susceptibility loci for NPC risk.
At the chromosome-wide level, we observed different contribution of the X chromosome to NPC heritability between genders. The genetic variance of NPC could be partially explained by X chromosome for males, with considerable portion next to that by the chromosome 6 harboring the well-known susceptible HLA loci. However, a subtle genetic contribution of the X chromosome to NPC was observed in females. The different genetic effects on males and females might explain the sexual dimorphism phenomena in NPC incidence; however, the mechanisms underlying this phenomenon remain unclear. Random inactivation in the X chromosome is a major characteristic for females, such that homozygotes in females exhibit a similar genetic effect as hemizygotes in males. It has been implicated that X-linked tumor suppressor genes may escape from X-inactivation, whereby the females carrying such two copies of functional genes have a reduced risk to develop tumor . Moreover, sex-specific regulatory variants are likely enriched in the X chromosome due to the sex-specific chromatin accessibility, implying different cumulative effects on gene reregulation between genders and thereby sex-biases in disease prevalence .
We identified a novel susceptibility locus at rs5927056 to be associated with NPC risk, surpassing X chromosome-wide significance in the discovery sample collection. Moreover, the combined analysis revealed a consistent association in males. Heterogeneous effect of the SNP was found in females (OR = 0.77, P = 0.04 in discovery data; OR = 0.99 and 1.00 in two replications), which could be partially explained by the inadequate statistical power because of limited female individuals included. Furthermore, the genetic variance of NPC risk to be explained by the X chromosome SNPs decreased with 8.3% after exclusion of rs5927056, suggesting a remarkable genetic effect of this locus on NPC susceptibility. rs5927056 is located within the intron of the DMD gene, a large gene that encodes the component of dystrophin-glycoprotein complex and bridges the inner cytoskeleton and the extracellular matrix. It has been demonstrated that DMD is sensitive to the replication stress on DNA damage and genome instability in tumor cells  and may act as a tumor suppressor involved in the development and progression of mesenchymal tumors . Alteration in DMD was shown to be associated with many non-myogenic tumors . Considering that rs5927056 alters some motifs in the region as revealed by HaploReg, we suspect that the haplotype of rs5927056 might regulate DMD expression and thus predispose individuals to NPC.
We also observed that rs5933886 in ARHGAP6 was associated with NPC risk, specifically in females. ARHGAP6 encodes Ras homology GTPase activation protein 6 and is involved in the regulation of actin polymerization at the plasma membrane during several cellular processes. ARHGAP6 may act as a tumor suppressor by inhibiting cell proliferation, migration, invasion, and adhesion of cervical carcinoma . Our bioinformatics analyses revealed that rs5933886 was located within a regulatory region and was associated with cis-eQTL effects on ARHGAP6, suggesting that it may confer a protective effect on NPC risk for females through transcriptional regulation of ARHGAP6.
Although our study adopted a case-control approach with more than 3900 samples across multiple centers, the sample size is still our limitation. Only one third of the heritability on X chromosome could be explained by the identified association signals in current study. Given the observed minor allele frequency and effect size, the power of association tests in discovery dataset varied from 0.097 (rs371000 in males) to 0.841 (rs6418572 in females; Additional file 1: Table S9) under the present significance level (α = 2.17 × 10−5), suggesting that we had certain power to discover some susceptibility loci underlying the X chromosome and meanwhile other potential susceptibility loci would be overlooked under the current sample size setting. Moreover, as compared to that of the discovery cohort, the relatively smaller sample sizes of the replication cohorts might be insufficient to establish associations with effect size smaller than that in the discovery due to Winner’s curse phenomenon. Therefore, further studies with larger cohorts are necessary to boost the power for replicating our findings as well as detecting missing susceptibility loci on X chromosome, in particular for SNPs with low minor allele frequency and weak effect size.
Taken together, we reported for the first time the large-scale X chromosome-wide association study for NPC. Our study illustrates the value of association test by including the X chromosome in finding the missing heritability of complex diseases other than the conventional GWAS, particularly for those diseases showing the remarkable difference in incidence rate between genders. We acknowledge that future association studies are warranted to validate our findings in independent cohorts with a large sample size and to identify more novel genetic variants on sexual chromosomes associated with NPC risk. Moreover, fine-mapping studies incorporating next-generation sequencing technology would be important to pinpoint the causal variants underlying the susceptibility loci.
Our study illustrates that the X chromosome confers different genetic effects on males and females and might explain the sexual dimorphism phenomena in NPC incidence. Our study reveals new susceptibility loci at the X chromosome conferring risk of NPC and supports the value of including the X chromosome in large-scale association studies. Further association studies using independent large cohorts are warranted to validate our findings.
We thank all the participants in the study, and staff members at the biobank of Sun Yat-Sen University Cancer Center (China) for their generous contribution in preparing patient samples.
This work was supported by the National Key Research and Development Program of China [2016YFC0902001]; the National High Technology Research and Development Program of China [2012AA02A206]; the National Natural Science Foundation of China [81572781, 81602478, 81372882, 81222035]; the China Postdoctoral Science Foundation [2015 M580759]; the Natural Science Foundation of Guangdong Province [2016A030310194]; Guangdong Innovative and Entrepreneurial Research Team Program [2016ZT06S638]; the National Program for Support of Top-Notch Young Professionals; Chang Jiang Scholars Program; the Special Support Program of Guangdong, and the Program for New Century Excellent Talents in University [NCET-11-0529].
Availability of data and materials
The datasets generated and/or analyzed during the current study are not publicly available due to restrictions apply to the availability of these data, which were used under license for the current study.
JXB and YXZ designed and directed the study. QSF, YMG, PPW, WHJ, YFX, and XYZ prepared the samples and obtained the genotypes for the Southern Chinese datasets. WHS prepared samples and obtained genotypes for the Taiwan Chinese dataset. ASK, TM, YMC, and CCN prepared samples and obtained genotypes for the Malaysian Chinese dataset. YMC analyzed data for Malaysian subjects. XYZ analyzed the data. XYZ, BL, JS, and JJL interpreted the data. XYZ wrote the report. All authors approved the final version.
Ethics approval and consent to participate
The study was approved by the Institutional Review Board at SYSUCC. Informed consent was obtained from all participants.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Curado MP, Edwards B, Shin HR, Storm H, Ferlay J, Heanue M, et al. Cancer incidence in five continents, Vol. IX. Lyon: IARC; 2007.Google Scholar
- Bei JX, Zuo XY, Liu WS, Guo YM, Zeng YX. Genetic susceptibility to the endemic form of NPC. Chin Clin Oncol. 2016;5(2):15.View ArticleGoogle Scholar
- Tse KP, Su WH, Chang KP, Tsang NM, Yu CJ, Tang P, et al. Genome-wide association study reveals multiple nasopharyngeal carcinoma-associated loci within the HLA region at chromosome 6p21.3. Am J Hum Genet. 2009;85(2):194–203.View ArticleGoogle Scholar
- Ng CC, Yew PY, Puah SM, Krishnan G, Yap LF, Teo SH, et al. A genome-wide association study identifies ITGA9 conferring risk of nasopharyngeal carcinoma. J Hum Genet. 2009;54(7):392–7.View ArticleGoogle Scholar
- Bei JX, Li Y, Jia WH, Feng BJ, Zhou G, Chen LZ, et al. A genome-wide association study of nasopharyngeal carcinoma identifies three new susceptibility loci. Nat Genet. 2010;42(7):599–603.View ArticleGoogle Scholar
- Tang M, Lautenberger JA, Gao X, Sezgin E, Hendrickson SL, Troyer JL, et al. The principal genetic determinants for nasopharyngeal carcinoma in China involve the HLA class I antigen recognition groove. PLoS Genet. 2012;8(11):e1003103.View ArticleGoogle Scholar
- Cui Q, Zuo XY, Lian YF, Feng QS, Xia YF, He CY, et al. Association between XRCC3 Thr241Met polymorphism and nasopharyngeal carcinoma risk: evidence from a large-scale case-control study and a meta-analysis. Tumour Biol. 2016;37(11):14825–30.View ArticleGoogle Scholar
- Cui Q, Feng QS, Mo HY, Sun J, Xia YF, Zhang H, et al. An extended genome-wide association study identifies novel susceptibility loci for nasopharyngeal carcinoma. Hum Mol Genet. 2016;25(16):3626–34.View ArticleGoogle Scholar
- Yee Ko JM, Dai W, Wun Wong EH, Kwong D, Tong Ng W, Lee A, et al. Multigene pathway-based analyses identify nasopharyngeal carcinoma risk associations for cumulative adverse effects of TERT-CLPTM1L and DNA double-strand breaks repair. Int J Cancer. 2014;135(7):1634–45.View ArticleGoogle Scholar
- Bei JX, Su WH, Ng CC, Yu K, Chin YM, Lou PJ, et al. A GWAS meta-analysis and replication study identifies a novel locus within CLPTM1L/TERT associated with nasopharyngeal carcinoma in individuals of Chinese ancestry. Cancer Epidemiol Biomark Prev. 2016;25(1):188–92.View ArticleGoogle Scholar
- Eichler EE, Flint J, Gibson G, Kong A, Leal SM, Moore JH, et al. Missing heritability and strategies for finding the underlying causes of complex disease. Nat Rev Genet. 2010;11(6):446–50.View ArticleGoogle Scholar
- Manolio TA, Collins FS, Cox NJ, Goldstein DB, Hindorff LA, Hunter DJ, et al. Finding the missing heritability of complex diseases. Nature. 2009;461(7265):747–53.View ArticleGoogle Scholar
- Chang ET, Adami HO. The enigmatic epidemiology of nasopharyngeal carcinoma. Cancer Epidemiol Biomark Prev. 2006;15(10):1765–77.View ArticleGoogle Scholar
- Yu MC, Yuan JM. Epidemiology of nasopharyngeal carcinoma. Semin Cancer Biol. 2002;12(6):421–9.View ArticleGoogle Scholar
- Xie SH, Yu IT, Tse LA, Mang OW, Yue L. Sex difference in the incidence of nasopharyngeal carcinoma in Hong Kong 1983–2008: suggestion of a potential protective role of oestrogen. Eur J Cancer. 2013;49(1):150–5.View ArticleGoogle Scholar
- Zhang LF, Li YH, Xie SH, Ling W, Chen SH, Liu Q, et al. Incidence trend of nasopharyngeal carcinoma from 1987 to 2011 in Sihui County, Guangdong Province, South China: an age-period-cohort analysis. Chin J Cancer. 2015;34(6):15.View ArticleGoogle Scholar
- Dunford A, Weinstock DM, Savova V, Schumacher SE, Cleary JP, Yoda A, et al. Tumor-suppressor genes that escape from X-inactivation contribute to cancer sex bias. Nat Genet. 2017;49(1):10–6.View ArticleGoogle Scholar
- Wee JT, Ha TC, Loong SL, Qian CN. Is nasopharyngeal cancer really a “Cantonese cancer”? Chin J Cancer. 2010;29(5):517–26.View ArticleGoogle Scholar
- Simons MJ. Nasopharyngeal carcinoma as a paradigm of cancer genetics. Chin J Cancer. 2011;30(2):79–84.View ArticleGoogle Scholar
- Conde L, Foo JN, Riby J, Liu J, Darabi H, Hjalgrim H, et al. X chromosome-wide association study of follicular lymphoma. Br J Haematol. 2013;162(6):858–62.View ArticleGoogle Scholar
- Wong EH, So HC, Li M, Wang Q, Butler AW, Paul B, et al. Common variants on Xq28 conferring risk of schizophrenia in Han Chinese. Schizophr Bull. 2014;40(4):777–86.View ArticleGoogle Scholar
- Zhu Z, Liang Z, Liany H, Yang C, Wen L, Lin Z, et al. Discovery of a novel genetic susceptibility locus on X chromosome for systemic lupus erythematosus. Arthritis Res Ther. 2015;17:349.View ArticleGoogle Scholar
- Chin YM, Mushiroda T, Takahashi A, Kubo M, Krishnan G, Yap LF, et al. HLA-A SNPs and amino acid variants are associated with nasopharyngeal carcinoma in Malaysian Chinese. Int J Cancer. 2015;136(3):678–87.PubMedGoogle Scholar
- Konig IR, Loley C, Erdmann J, Ziegler A. How to include chromosome X in your genome-wide association study. Genet Epidemiol. 2014;38(2):97–103.View ArticleGoogle Scholar
- Ziegler A. Genome-wide association studies: quality control and population-based measures. Genet Epidemiol. 2009;33(Suppl 1):S45–50.View ArticleGoogle Scholar
- Howie BN, Donnelly P, Marchini J. A flexible and accurate genotype imputation method for the next generation of genome-wide association studies. PLoS Genet. 2009;5(6):e1000529.View ArticleGoogle Scholar
- van Leeuwen EM, Kanterakis A, Deelen P, Kattenberg MV, Genome of the Netherlands C, Slagboom PE, et al. population-specific genotype imputations using minimac or IMPUTE2. Nat Protoc 2015;10(9):1285–1296.Google Scholar
- Clayton D. Testing for association on the X chromosome. Biostatistics. 2008;9(4):593–600.View ArticleGoogle Scholar
- Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.View ArticleGoogle Scholar
- Yang J, Manolio TA, Pasquale LR, Boerwinkle E, Caporaso N, Cunningham JM, et al. Genome partitioning of genetic variation for complex traits using common SNPs. Nat Genet. 2011;43(6):519–25.View ArticleGoogle Scholar
- Yang J, Lee SH, Goddard ME, Visscher PM. GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 2011;88(1):76–82.View ArticleGoogle Scholar
- Tukiainen T, Pirinen M, Sarin AP, Ladenvall C, Kettunen J, Lehtimaki T, et al. Chromosome X-wide association study identifies loci for fasting insulin and height and evidence for incomplete dosage compensation. PLoS Genet. 2014;10(2):e1004127.View ArticleGoogle Scholar
- Gao X, Starmer J, Martin ER. A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms. Genet Epidemiol. 2008;32(4):361–9.View ArticleGoogle Scholar
- Skol AD, Scott LJ, Abecasis GR, Boehnke M. Joint analysis is more efficient than replication-based analysis for two-stage genome-wide association studies. Nat Genet. 2006;38(2):209–13.View ArticleGoogle Scholar
- Ward LD, Kellis M. HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 2012;40(Database issue):D930–4.View ArticleGoogle Scholar
- Kukurba KR, Parsana P, Balliu B, Smith KS, Zappala Z, Knowles DA, et al. Impact of the X chromosome and sex on regulatory variation. Genome Res. 2016;26(6):768–77.View ArticleGoogle Scholar
- Durkin SG, Glover TW. Chromosome fragile sites. Annu Rev Genet. 2007;41:169–92.View ArticleGoogle Scholar
- Wang Y, Marino-Enriquez A, Bennett RR, Zhu M, Shen Y, Eilers G, et al. Dystrophin is a tumor suppressor in human cancers with myogenic programs. Nat Genet. 2014;46(6):601–6.View ArticleGoogle Scholar
- Luce LN, Abbate M, Cotignola J, Giliberto F. Non-myogenic tumors display altered expression of dystrophin (DMD) and a high frequency of genetic alterations. Oncotarget. 2017;8(1):145–55.Google Scholar
- Li J, Liu Y, Yin Y. Inhibitory effects of Arhgap6 on cervical carcinoma cells. Tumour Biol. 2016;37(2):1411–25.View ArticleGoogle Scholar