Candidate Gene Analysis of Tooth Agenesis Identifies Novel Mutations in Six Genes and Suggests Significant Role for WNT and EDA Signaling and Allele Combinations

Failure to develop complete dentition, tooth agenesis, is a common developmental anomaly manifested most often as isolated but also as associated with many developmental syndromes. It typically affects third molars or one or few other permanent teeth but severe agenesis is also relatively prevalent. Here we report mutational analyses of seven candidate genes in a cohort of 127 probands with non-syndromic tooth agenesis. 82 lacked more than five permanent teeth excluding third molars, called as oligodontia. We identified 28 mutations, 17 of which were novel. Together with our previous reports, we have identified two mutations in MSX1, AXIN2 and EDARADD, five in PAX9, four in EDA and EDAR, and nine in WNT10A. They were observed in 58 probands (44%), with a mean number of missing teeth of 11.7 (range 4 to 34). Almost all of these probands had severe agenesis. Only few of the probands but several relatives with heterozygous genotypes of WNT10A or EDAR conformed to the common type of non-syndromic tooth agenesis, incisor-premolar hypodontia. Mutations in MSX1 and PAX9 affected predominantly posterior teeth, whereas both deciduous and permanent incisors were especially sensitive to mutations in EDA and EDAR. Many mutations in EDAR, EDARADD and WNT10A were present in several families. Biallelic or heterozygous genotypes of WNT10A were observed in 32 and hemizygous or heterozygous genotypes of EDA, EDAR or EDARADD in 22 probands. An EDARADD variant were in seven probands present together with variants in EDAR or WNT10A, suggesting combined phenotypic effects of alleles in distinct genes.     

WNT10A variants in relation to nonsyndromic hypodontia in eastern Slovak population

Nonsyndromic hypodontia is a congenital absence of less than six permanent teeth, with a most common subtype maxillary lateral incisor agenesis (MLIA). Mutations in several genes have been described in severe tooth agenesis. The aim of this study was to search for the variants in wingless-type MMTV-integration site family member (WNT10A), paired box 9 (PAX9) and axis inhibitor 2 (AXIN2) genes, and investigate their potential role in the pathogenesis of non-syndromic hypodontia. Clinical examination and panoramic radiograph were performed in the cohort of 60 unrelated Slovak patients of Caucasian origin with nonsyndromic hypodontia including 37 MLIA cases and 48 healthy controls. Genomic DNA was isolated from buccal swabs and Sanger sequencing of WNT10A, PAX9 and AXIN2 was performed. Altogether, we identified 23 single-nucleotide variants, of which five were novel. We have found three rare nonsynonymous variants in WNT10A (p.Gly165Arg; p.Gly213Ser and p.Phe228Ile) in eight (13.33%) of 60 patients. Analysis showed potentially damaged WNT10A variant p.Phe228Ile predominantly occurred only in MLIA patients, and with a dominant form of tooth agenesis (odds ratio \(({\hbox {OR}}_{\mathrm{dom}}) = 9.841\); \(P=0.045\); 95% confidence interval (CI) 0.492–196.701; \({\hbox {OR}}_{\mathrm{rec}} = 0.773\); \(P =1.000\); 95% CI 0.015–39.877). In addition, the WNT10A variant p.Phe228Ile showed a trend associated with familial nonsyndromic hypodontia (\(P =0.024\); OR = 1.20; 95% CI 0.97–1.48). After Bonferroni correction, these effects remained with borderline tendencies. Using a 3D WNT10A protein model, we demonstrated that the variant Phe228Ile changes the protein secondary structure. In PAX9 and AXIN2, common variants were detected. Our findings suggest that the identified WNT10A variant p.Phe228Ile could represent risk for the inherited nonsyndromic hypodontia underlying MLIA. However, further study in different populations is required.     

Loss-of-Function Mutations in the WNT Co-receptor LRP6 Cause Autosomal-Dominant Oligodontia

Tooth agenesis is one of the most common developmental anomalies in man. Oligodontia, a severe form of tooth agenesis, occurs both as an isolated anomaly and as a syndromal feature. We performed exome sequencing on 20 unrelated individuals with apparent non-syndromic oligodontia and failed to detect mutations in genes previously associated with oligodontia. In three of the probands, we detected heterozygous variants in LRP6, and sequencing of additional oligodontia-affected individuals yielded one additional mutation in LRP6. Three mutations (c.1144_1145dupAG [p.Ala383Glyfs^∗8], c.1779dupT [p.Glu594^∗], and c.2224_2225dupTT [p.Leu742Phefs^∗7]) are predicted to truncate the protein, whereas the fourth (c.56C>T [p.Ala19Val]) is a missense variant of a conserved residue located at the cleavage site of the protein’s signal peptide. All four affected individuals harboring a LRP6 mutation had a family history of tooth agenesis. LRP6 encodes a transmembrane cell-surface protein that functions as a co-receptor with members from the Frizzled protein family in the canonical Wnt/β-catenin signaling cascade. In this same pathway, WNT10A was recently identified as a major contributor in the etiology of non-syndromic oligodontia. We show that the LRP6 missense variant (c.56C>T) results in altered glycosylation and improper subcellular localization of the protein, resulting in abrogated activation of the Wnt pathway. Our results identify LRP6 variants as contributing to the etiology of non-syndromic autosomal-dominant oligodontia and suggest that this gene is a candidate for screening in DNA diagnostics.     

Involvement of and Interaction between WNT10A and EDA Mutations in Tooth Agenesis Cases in the Chinese Population

Background

Dental agenesis is the most common, often heritable, developmental anomaly in humans. Although WNT10A gene mutations are known to cause rare syndromes associated with tooth agenesis, including onycho-odontodermal dysplasia (OODD), Schöpf-Schulz-Passarge syndrome (SSPS), hypohidrotic ectodermal dysplasia (HED), and more than half of the cases of isolated oligodontia recently, the genotype-phenotype correlations and the mode of inheritance of WNT10A mutations remain unclear. The phenotypic expression with WNT10A mutations shows a high degree of variability, suggesting that other genes might function with WNT10A in regulating ectodermal organ development. Moreover, the involvement of mutations in other genes, such as EDA, which is also associated with HED and isolated tooth agenesis, is not clear. Therefore, we hypothesized that EDA mutations interact with WNT10A mutations to play a role in tooth agenesis. Additionally, EDA, EDAR, and EDARADD encode signaling molecules in the Eda/Edar/NF-κB signaling pathways, we also checked EDAR and EDARADD in this study.

Methods

WNT10A, EDA, EDAR and EDARADD were sequenced in 88 patients with isolated oligodontia and 26 patients with syndromic tooth agenesis. The structure of two mutated WNT10A and two mutated EDA proteins was analyzed.

Results

Digenic mutations of both WNT10A and EDA were identified in 2 of 88 (2.27%) isolated oligodontia cases and 4 of 26 (15.38%) syndromic tooth agenesis cases. No mutation in EDAR or EDARADD gene was found.

Conclusions

WNT10A and EDA digenic mutations could result in oligodontia and syndromic tooth agenesis in the Chinese population. Moreover, our results will greatly expand the genotypic spectrum of tooth agenesis.     

An Aberrant Splice Acceptor Site Due to a Novel Intronic Nucleotide Substitution in MSX1 Gene Is the Cause of Congenital Tooth Agenesis in a Japanese Family

Congenital tooth agenesis is caused by mutations in the MSX1, PAX9, WNT10A, or AXIN2 genes. Here, we report a Japanese family with nonsyndromic tooth agenesis caused by a novel nucleotide substitution in the intronic region between exons 1 and 2 of the MSX1 gene. Because the mutation is located 9 bp before exon 2 (c.452-9G>A), we speculated that the nucleotide substitution would generate an abnormal splice site. Using cDNA analysis of an immortalized patient blood cell, we confirmed that an additional 7-nucleotide sequence was inserted at the splice junction between exons 1 and 2 (c.451_452insCCCTCAG). The consequent frameshift generated a homeodomain-truncated MSX1 (p.R151fsX20). We then studied the subcellular localization of truncated MSX1 protein in COS cells, and observed that it had a whole cell distribution more than a nuclear localization, compared to that of wild-type protein. This result suggests a deletion of the nuclear localization signal, which is mapped to the MSX1 homeodomain. These results indicate that this novel intronic nucleotide substitution is the cause of tooth agenesis in this family. To date, most MSX1 variants isolated from patients with tooth agenesis involve single amino acid substitutions in the highly conserved homeodomain or deletion mutants caused by frameshift or nonsense mutations. We here report a rare case of an intronic mutation of the MSX1 gene responsible for human tooth agenesis. In addition, the missing tooth patterns were slightly but significantly different between an affected monozygotic twin pair of this family, showing that epigenetic or environmental factors also affect the phenotypic variations of missing teeth among patients with nonsyndromic tooth agenesis caused by an MSX1 haploinsufficiency.     

The adaptive variant EDARV370A is associated with straight hair in East Asians

Hair straightness/curliness is a highly heritable trait amongst human populations. Previous studies have reported European specific genetic variants influencing hair straightness, but those in East Asians remain unknown. One promising candidate is a derived coding variant of the ectodysplasin A receptor (EDAR), EDARV370A (370A), associated with several phenotypic changes of epidermal appendages. One of the strongest signals of natural selection in human genomes, 370A, has risen to high prevalence in East Asian and Native American populations, whilst being almost absent in Europeans and Africans. This striking frequency distribution and the pleiotropic nature of 370A led us to pursue if hair straightness, another epidermal appendage-related phenotype, is affected by this variant. By studying 1,718 individuals from four distinctive East Asian populations (Han, Tibetan, Mongolian, and Li), we found a significant association between 370A and the straight hair type in the Han (p = 2.90 × 10^?6), Tibetan (p = 3.07 × 10^?2), and Mongolian (p = 1.03 × 10^?5) populations. Combining all the samples, the association is even stronger (p = 5.18 × 10^?10). The effect of 370A on hair straightness is additive, with an odds ratio of 2.05. The results indicate very different biological mechanisms of straight hair in Europe and Asia, and also present a more comprehensive picture of the phenotypic consequences of 370A, providing important clues into the potential adaptive forces shaping the evolution of this extraordinary genetic variant.     

Msx1 Gene Variant - its Association in Isolated Hypodontia: A Case Control Genetic study!!!

INTRODUCTION:

Non-syndromic tooth agenesis is a congenital anomaly with significant medical, psychological, and social ramifications. There is sufficient evidence to hypothesize that locus for this condition can be identified by candidate genes.

AIM OF THE STUDY:

The aim of this study was to test whether MSX1 671 T > C gene variant was involved in etiology of non-syndromic tooth agenesis in Raichur patients.

MATERIALS AND METHODS:

Blood samples were collected with informed consent from 50 subjects having non-syndromic tooth agenesis and 50 controls. Genomic deoxyribonucleic acid (DNA) was extracted from the blood samples, polymerase chain reaction (PCR) was performed, and restriction fragment length polymorphism (RFLP) was performed for digestion products that were evaluated.

RESULTS:

The results showed positive correlation between MSX1671 T > C gene variant and non-syndromic tooth agenesis in Raichur patients.

CONCLUSION:

MSX1 671 T > C gene variant may be a good screening marker for non-syndromic tooth agenesis in Raichur patients.     

WNT10A Mutations Are a Frequent Cause of a Broad Spectrum of Ectodermal Dysplasias with Sex-Biased Manifestation Pattern in Heterozygotes

Odonto-onycho-dermal dysplasia (OODD), a rare autosomal-recessive inherited form of ectodermal dysplasia including severe oligodontia, nail dystrophy, palmoplantar hyperkeratosis, and hyperhidrosis, was recently shown to be caused by a homozygous nonsense WNT10A mutation in three consanguineous Lebanese families. Here, we report on 12 patients, from 11 unrelated families, with ectodermal dysplasia caused by five previously undescribed WNT10A mutations. In this study, we show that (1) WNT10A mutations cause not only OODD but also other forms of ectodermal dysplasia, reaching from apparently monosymptomatic severe oligodontia to Schöpf-Schulz-Passarge syndrome, which is so far considered a unique entity by the findings of numerous cysts along eyelid margins and the increased risk of benign and malignant skin tumors; (2) WNT10A mutations are a frequent cause of ectodermal dysplasia and were found in about 9% of an unselected patient cohort; (3) about half of the heterozygotes (53.8%) show a phenotype manifestation, including mainly tooth and nail anomalies, which was not reported before in OODD; and (4) heterozygotes show a sex-biased manifestation pattern, with a significantly higher proportion of tooth anomalies in males than in females, which may implicate gender-specific differences of WNT10A expression.     

A genome-wide association study of severe teenage acne in European Americans

Despite the family aggregation of severe teenage acne, the genetic basis of this common skin condition remains unclear. We conducted a genome-wide association study (GWAS) on severe teenage acne in 928 European Americans. The SNP rs4133274 on chromosome 8q24 (72 kb upstream of MYC) revealed the most significant association with severe teenage acne (p value = 1.7 × 10^?6). The variant allele of this SNP (G allele) was associated with an increased risk of severe teenage acne with odds ratio of 4.01 (95 % confidence interval = 2.37–6.82). Upon further replication, our findings suggest new genetic basis of acne and may explain the association between acne and cancer risk observed in the epidemiological studies.     

Novel PAX9 and COL1A2 Missense Mutations Causing Tooth Agenesis and OI/DGI without Skeletal Abnormalities

Inherited dentin defects are classified into three types of dentinogenesis imperfecta (DGI) and two types of dentin dysplasia (DD). The genetic etiology of DD-I is unknown. Defects in dentin sialophosphoprotein (DSPP) cause DD type II and DGI types II and III. DGI type I is the oral manifestation of osteogenesis imperfecta (OI), a systemic disease typically caused by defects in COL1A1 or COL1A2. Mutations in MSX1, PAX9, AXIN2, EDA and WNT10A can cause non-syndromic familial tooth agenesis. In this study a simplex pattern of clinical dentinogenesis imperfecta juxtaposed with a dominant pattern of hypodontia (mild tooth agenesis) was evaluated, and available family members were recruited. Mutational analyses of the candidate genes for DGI and hypodontia were performed and the results validated. A spontaneous novel mutation in COL1A2 (c.1171G>A; p.Gly391Ser) causing only dentin defects and a novel mutation in PAX9 (c.43T>A; p.Phe15Ile) causing hypodontia were identified and correlated with the phenotypic presentations in the family. Bone radiographs of the proband’s dominant leg and foot were within normal limits. We conclude that when no DSPP mutation is identified in clinically determined isolated DGI cases, COL1A1 and COL1A2 should be considered as candidate genes. PAX9 mutation p.Phe15Ile within the N-terminal β-hairpin structure of the PAX9 paired domain causes tooth agenesis.     

Polymorphism of dental formula and segregation of its variants in a pedigree of Kerry Blue Terrier dogs

Polymorphism of the dental formula was analyzed in a complex pedigree of Kerry Blue Terrier. A lack of one or more lower premolars was observed in some dogs. Two different patterns of missing teeth were identified. One pattern consisted in agenesis of a second premolar, often in combination with agenesis of neighbor teeth, including the fourth premolar. In the second pattern, agenesis of a fourth premolar was expressed as an isolated abnormality. It was shown previously that the first pattern is inherited as a recessive trait with near complete penetrance. In this work, the major-gene control was demonstrated for the second pattern. This abnormality develops in 70–80% of mutant homozygotes and in no more than 20% of heterozygotes and wild-type homozygotes. It was shown that the two dentition abnormalities are controlled by different genes, which were designated LPA2 and LPA4 (Lower Premolar Agenesis).     

Association between common alcohol dehydrogenase gene (ADH) variants and schizophrenia and autism

Humans express at least seven alcohol dehydrogenase (ADH) isoforms that are encoded by ADH gene cluster (ADH7–ADH1C–ADH1B–ADH1A–ADH6–ADH4–ADH5) at chromosome 4. ADHs are key catabolic enzymes for retinol and ethanol. The functional ADH variants (mostly rare) have been implicated in alcoholism risk. In addition to catalyzing the oxidation of retinol and ethanol, ADHs may be involved in the metabolic pathways of several neurotransmitters that are implicated in the neurobiology of neuropsychiatric disorders. In the present study, we comprehensively examined the associations between common ADH variants [minor allele frequency (MAF) >0.05] and 11 neuropsychiatric and neurological disorders. A total of 50,063 subjects in 25 independent cohorts were analyzed. The entire ADH gene cluster was imputed across these 25 cohorts using the same reference panels. Association analyses were conducted, adjusting for multiple comparisons. We found 28 and 15 single nucleotide polymorphisms (SNPs), respectively, that were significantly associated with schizophrenia in African-Americans and autism in European-Americans after correction by false discovery rate (FDR) (q < 0.05); and 19 and 6 SNPs, respectively, that were significantly associated with these two disorders after region-wide correction by SNPSpD (8.9 × 10^?5 ≤ p  ≤ 0.0003 and 2.4 × 10^?5 ≤ p ≤ 0.0003, respectively). No variants were significantly associated with the other nine neuropsychiatric disorders, including alcohol dependence. We concluded that common ADH variants conferred risk for both schizophrenia in African-Americans and autism in European-Americans.     

Human Thromboxane A2 Receptor Genetic Variants: In Silico,In Vitro and “In Platelet” Analysis

Thromboxane and its receptor have emerged as key players in modulating vascular thrombotic events. Thus, a dysfunctional hTP genetic variant may protect against (hypoactivity) or promote (hyperactivity) vascular events, based upon its activity on platelets. After extensive in silico analysis, six hTP-α variants were selected (C⁶⁸S, V⁸⁰E, E⁹⁴V, A¹⁶⁰T, V¹⁷⁶E, and V²¹⁷I) for detailed biochemical studies based on structural proximity to key regions involved in receptor function and in silico predictions. Variant biochemical profiles ranged from severe instability (C⁶⁸S) to normal (V²¹⁷I), with most variants demonstrating functional alteration in binding, expression or activation (V⁸⁰E, E⁹⁴V, A¹⁶⁰T, and V¹⁷⁶E). In the absence of patient platelet samples, we developed and validated a novel megakaryocyte based system to evaluate human platelet function in the presence of detected dysfunctional genetic variants. Interestingly, variant V80E exhibited reduced platelet activation whereas A160T demonstrated platelet hyperactivity. This report provides the most comprehensive in silico, in vitro and “in platelet” evaluation of hTP variants to date and highlightscurrent inherent problems in evaluating genetic variants, with possible solutions. The study additionally provides clinical relevance to characterized dysfunctional hTP variants.     

Molecular genetics of <Emphasis Type="Italic">MARVELD2</Emphasis> and clinical phenotype in Pakistani and Slovak families segregating DFNB49 hearing loss

Pathogenic mutations of MARVELD2, encoding tricellulin, a tricelluar tight junction protein, cause autosomal recessive non-syndromic hearing loss (DFNB49) in families of Pakistan and Czech Roma origin. In fact, they are a significant cause of prelingual hearing loss in the Czech Roma, second only to GJB2 variants. Previously, we reported that mice homozygous for p.Arg497* variant of Marveld2 had a broad phenotypic spectrum, where defects were observed in the inner ear, heart, mandibular salivary gland, thyroid gland and olfactory epithelium. The current study describes the types and frequencies of MARVELD2 alleles and clinically reexamines members of DFNB49 families. We found that MARVELD2 variants are responsible for about 1.5 % (95 % CI 0.8–2.6) of non-syndromic hearing loss in our cohort of 800 Pakistani families. The c.1331+2T>C allele is recurrent. In addition, we identified a novel large deletion in a single family, which appears to have resulted from non-allelic homologous recombination between two similar Alu short interspersed elements. Finally, we observed no other clinical manifestations co-segregating with hearing loss in DFNB49 human families, and hypothesize that the additional abnormalities in the Marveld2 mutant mouse indicates a critical non-redundant function for tricellulin in other organ systems.     

Large multiethnic Candidate Gene Study for C-reactive protein levels: identification of a novel association at CD36 in African Americans

C-reactive protein (CRP) is a heritable biomarker of systemic inflammation and a predictor of cardiovascular disease (CVD). Large-scale genetic association studies for CRP have largely focused on individuals of European descent. We sought to uncover novel genetic variants for CRP in a multiethnic sample using the ITMAT Broad-CARe (IBC) array, a custom 50,000 SNP gene-centric array having dense coverage of over 2,000 candidate CVD genes. We performed analyses on 7,570 African Americans (AA) from the Candidate gene Association Resource (CARe) study and race-combined meta-analyses that included 29,939 additional individuals of European descent from CARe, the Women’s Health Initiative (WHI) and KORA studies. We observed array-wide significance (p < 2.2 × 10^?6) for four loci in AA, three of which have been reported previously in individuals of European descent (IL6R, p = 2.0 × 10^?6; CRP, p = 4.2 × 10^?71; APOE, p = 1.6 × 10^?6). The fourth significant locus, CD36 (p = 1.6 × 10^?6), was observed at a functional variant (rs3211938) that is extremely rare in individuals of European descent. We replicated the CD36 finding (p = 1.8 × 10^?5) in an independent sample of 8,041 AA women from WHI; a meta-analysis combining the CARe and WHI AA results at rs3211938 reached genome-wide significance (p = 1.5 × 10^?10). In the race-combined meta-analyses, 13 loci reached significance, including ten (CRP, TOMM40/APOE/APOC1, HNF1A, LEPR, GCKR, IL6R, IL1RN, NLRP3, HNF4A and BAZ1B/BCL7B) previously associated with CRP, and one (ARNTL) previously reported to be nominally associated with CRP. Two novel loci were also detected (RPS6KB1, p = 2.0 × 10^?6; CD36, p = 1.4 × 10^?6). These results highlight both shared and unique genetic risk factors for CRP in AA compared to populations of European descent.     

Case–control association study of polymorphisms in the voltage-gated sodium channel genes SCN1A,SCN2A,SCN3A,SCN1B,and SCN2B and epilepsy

High-frequency action potentials are mediated by voltage-gated sodium channels, composed of one large α subunit and two small β subunits, encoded mainly by SCN1A, SCN2A, SCN3A, SCN1B, and SCN2B genes in the brain. These play a key role in epilepsy, with the most commonly mutated gene in epilepsy being SCN1A. We examined whether polymorphisms in the above genes affect epilepsy risk in 1,529 epilepsy patients and 1,935 controls from four ethnicities or locations: Malay, Indian, and Chinese, all from Malaysia, and Chinese from Hong Kong. Of patients, 19 % were idiopathic, 42 % symptomatic, and 40 % cryptogenic. We genotyped 43 polymorphisms: 27 in Hong Kong, 28 in Malaysia, and 12 in both locations. The strongest association with epilepsy was rs3812718, or SCN1A IVS5N+5G>A: odds ratio (OR) = 0.85 for allele G (p = 0.0009) and 0.73 for genotype GG versus AA (p = 0.003). The OR was between 0.76 and 0.87 for all ethnicities. Meta-analysis confirmed the association (OR = 0.81 and p = 0.002 for G, and OR = 0.67 and p = 0.007 for GG versus AA), which appeared particularly strong for Indians and for febrile seizures. Allele G affects splicing and speeds recovery from inactivation. Since SCN1A is preferentially expressed in inhibitory neurons, G may decrease epilepsy risk. SCN1A rs10188577 displayed OR = 1.20 for allele C (p = 0.003); SCN2A rs12467383 had OR = 1.16 for allele A (p = 0.01), and displayed linkage disequilibrium with rs2082366 (r ² = 0.67), whose genotypes tended toward association with SCN2A brain expression (p = 0.10). SCN1A rs2298771 was associated in Indians (OR = 0.56, p = 0.005) and SCN2B rs602594 with idiopathic epilepsy (OR = 0.62, p = 0.002). Therefore, sodium channel polymorphisms are associated with epilepsy.     

Functional polymorphism rs4072037 in MUC1 gene contributes to the susceptibility to gastric cancer: evidence from pooled 6,580 cases and 10,324 controls

Genome-wide association studies have reported a promising association of rs4072037 with gastric cancer (GC). This variant was associated with altered physiological function of MUC1 possibly by modulating promoter activity and alternative splicing of MUC1. However, the association results were inconclusive and estimate of the effect of this variant was not well evaluated. A meta-analysis by systematically reviewing relevant reports may facilitate to address these concerns. Association studies involving MUC1 rs4072037 polymorphism and GC risk were identified up to June 30, 2012. Odds ratio (OR) and 95 % confidence interval (CI) in additive model were estimated or extracted from each study. The pooled effect size was quantitatively synthesized using meta-analysis. Heterogeneity between studies was measured by the Q test and I ² statistic, and publication bias was evaluated by a funnel plot and the Egger’s test. A total of 10 independent case–control studies including 6,580 GC cases and 10,324 controls were included in this meta-analysis. Eight of the ten studies were Asian ethnicity and two European. The G allele of MUC1 rs4072037 was significantly associated with a decreased risk of GC (OR = 0.72, 95 % CI 0.68–0.77; P = 7.82 × 10^?25), as compared with A allele. Stratification for different ethnicity, tumor localization or type showed similar results. These findings represent important evidence for association of MUC1 rs4072037 variant with GC risk, and also provide a relatively reliable estimate of effect size. MUC1 is a strong candidate as a susceptibility gene of GC.