Zhen Zhen Zhao

A Genome-Wide Association Study Identifies Novel Alleles Associated with Hair Color and Skin Pigmentation

We conducted a multi-stage genome-wide association study of natural hair color in more than 10,000 men and women of European ancestry from the United States and Australia. An initial analysis of 528,173 single nucleotide polymorphisms (SNPs) genotyped on 2,287 women identified IRF4 and SLC24A4 as loci highly associated with hair color, along with three other regions encompassing known pigmentation genes. We confirmed these associations in 7,028 individuals from three additional studies. Across these four studies, SLC24A4 rs12896399 and IRF4 rs12203592 showed strong associations with hair color, with p = 6.0×10−62 and p = 7.46×10−127, respectively. The IRF4 SNP was also associated with skin color (p = 6.2×10−14), eye color (p = 6.1×10−13), and skin tanning response to sunlight (p = 3.9×10−89). A multivariable analysis pooling data from the initial GWAS and an additional 1,440 individuals suggested that the association between rs12203592 and hair color was independent of rs1540771, a SNP between the IRF4 and EXOC2 genes previously found to be associated with hair color. After adjustment for rs12203592, the association between rs1540771 and hair color was not significant (p = 0.52). One variant in the MATP gene was associated with hair color. A variant in the HERC2 gene upstream of the OCA2 gene showed the strongest and independent association with hair color compared with other SNPs in this region, including three previously reported SNPs. The signals detected in a region around the MC1R gene were explained by MC1R red hair color alleles. Our results suggest that the IRF4 and SLC24A4 loci are associated with human hair color and skin pigmentation.

A Single SNP in an Evolutionary Conserved Region within Intron 86 of the HERC2 Gene Determines Human Blue-Brown Eye Color

We have previously demonstrated that haplotypes of three single nucleotide polymorphisms (SNPs) within the first intron of the OCA2 gene are extremely strongly associated with variation in human eye color. In the present work, we describe additional fine association mapping of eye color SNPs in the intergenic region upstream of OCA2 and within the neighboring HERC2 (hect domain and RLD2) gene. We screened an additional 92 SNPs in 300-3000 European individuals and found that a single SNP in intron 86 of HERC2, rs12913832, predicted eye color significantly better (ordinal logistic regression R(2) = 0.68, association LOD = 444) than our previous best OCA2 haplotype. Comparison of sequence alignments of multiple species showed that this SNP lies in the center of a short highly conserved sequence and that the blue-eye-associated allele (frequency 78%) breaks up this conserved sequence, part of which forms a consensus binding site for the helicase-like transcription factor (HLTF). We were also able to demonstrate the OCA2 R419Q, rs1800407, coding SNP acts as a penetrance modifier of this new HERC2 SNP for eye color, and somewhat independently, of melanoma risk. We conclude that the conserved region around rs12913832 represents a regulatory region controlling constitutive expression of OCA2 and that the C allele at rs12913832 leads to decreased expression of OCA2, particularly within iris melanocytes, which we postulate to be the ultimate cause of blue eye color.

Common sequence variants on 20q11.22 confer melanoma susceptibility

We conducted a genome-wide association pooling study for cutaneous melanoma and performed validation in samples totaling 2,019 cases and 2,105 controls. Using pooling, we identified a new melanoma risk locus on chromosome 20 (rs910873 and rs1885120), with replication in two further samples (combined P < 1 × 10−15). The per allele odds ratio was 1.75 (1.53, 2.01), with evidence for stronger association in early-onset cases.

Genome-wide association study identifies variants at 9p21 and 22q13 associated with development of cutaneous nevi

A high melanocytic nevi count is the strongest known risk factor for cutaneous melanoma. We conducted a genome-wide association study for nevus count using 297,108 SNPs in 1,524 twins, with validation in an independent cohort of 4,107 individuals. We identified strongly associated variants in MTAP, a gene adjacent to the familial melanoma susceptibility locus CDKN2A on 9p21 (rs4636294, combined P = 3.4 × 10−15), as well as in PLA2G6 on 22q13.1 (rs2284063, combined P = 3.4 × 10−8). In addition, variants in these two loci showed association with melanoma risk in 3,131 melanoma cases from two independent studies, including rs10757257 at 9p21, combined P = 3.4 × 10−8, OR = 1.23 (95% CI = 1.15–1.30) and rs132985 at 22q13.1, combined P = 2.6 × 10−7, OR = 1.23 (95% CI = 1.15–1.30). This provides the first report of common variants associated to nevus number and demonstrates association of these variants with melanoma susceptibility.

Multiple pigmentation gene polymorphisms account for a substantial proportion of risk of cutaneous malignant melanoma.

We have previously described the role of red hair (melanocortin-1 receptor, MC1R) and blue eye (oculocutaneous albinism type II, OCA2) gene polymorphisms in modulating the risk of cutaneous malignant melanoma (CMM) in a highly sun-exposed population of European descent. A number of recent studies, including genome-wide association studies, have identified numerous polymorphisms controlling human hair, eye, and skin color. In this paper, we test a selected set of polymorphisms in pigmentation loci (ASIP (Agouti signalling protein, nonagouti homolog (mouse) gene), TYR (tyrosinase), TYRP1 (tyrosinase-related protein 1), MC1R, OCA2, IRF4 (interferon regulatory factor 4), SLC24A4 (solute carrier family 24, member 4), and SLC45A2 (solute carrier family 45, member 2)) for association with CMM risk in a large Australian population-based case-control study. Variants in IRF4 and SLC24A4, despite being strongly associated with pigmentation in our sample, did not modify CMM risk, but the other six did. Three single nucleotide polymorphisms (rs28777, rs35391, and rs16891982) in the MATP gene (SLC45A2) exhibited the strongest crude association with risk, but this was attenuated to approximately the same effect size as that of a MC1R red hair color allele by controlling for ancestry of cases and controls. We also detected significant epistatic interactions between SLC45A2 and OCA2 alleles, and MC1R and ASIP alleles. Overall, these measured variants account for 12% of the familial risk of CMM in our population.

Genome-wide association study identifies a new melanoma susceptibility locus at 1q21.3

We performed a genome-wide association study of melanoma in a discovery cohort of 2,168 Australian individuals with melanoma and 4,387 control individuals. In this discovery phase, we confirm several previously characterized melanoma-associated loci at MC1R, ASIP and MTAP–CDKN2A. We selected variants at nine loci for replication in three independent case-control studies (Europe: 2,804 subjects with melanoma, 7,618 control subjects; United States 1: 1,804 subjects with melanoma, 1,026 control subjects; United States 2: 585 subjects with melanoma, 6,500 control subjects). The combined meta-analysis of all case-control studies identified a new susceptibility locus at 1q21.3 (rs7412746, P = 9.0 × 10−11, OR in combined replication cohorts of 0.89 (95% CI 0.85–0.95)). We also show evidence suggesting that melanoma associates with 1q42.12 (rs3219090, P = 9.3 × 10−8). The associated variants at the 1q21.3 locus span a region with ten genes, and plausible candidate genes for melanoma susceptibility include ARNT and SETDB1. Variants at the 1q21.3 locus do not seem to be associated with human pigmentation or measures of nevus density.

Highly cost-efficient genome-wide association studies using DNA pools and dense SNP arrays

Genome-wide association (GWA) studies to map genes for complex traits are powerful yet costly. DNA-pooling strategies have the potential to dramatically reduce the cost of GWA studies. Pooling using Affymetrix arrays has been proposed and used but the efficiency of these arrays has not been quantified. We compared and contrasted Affymetrix Genechip HindIII and Illumina HumanHap300 arrays on the same DNA pools and showed that the HumanHap300 arrays are substantially more efficient. In terms of effective sample size, HumanHap300-based pooling extracts >80% of the information available with individual genotyping (IG). In contrast, Genechip HindIII-based pooling only extracts ∼30% of the available information. With HumanHap300 arrays concordance with IG data is excellent. Guidance is given on best study design and it is shown that even after taking into account pooling error, one stage scans can be performed for >100-fold reduced cost compared with IG. With appropriately designed two stage studies, IG can provide confirmation of pooling results whilst still providing ∼20-fold reduction in total cost compared with IG-based alternatives. The large cost savings with Illumina HumanHap300-based pooling imply that future studies need only be limited by the availability of samples and not cost.

Association of Helicobacter pylori Infection With Reduced Risk for Esophageal Cancer Is Independent of Environmental and Genetic Modifiers

BACKGROUND & AIMS Infection with Helicobacter pylori is associated with reduced risk of esophageal adenocarcinoma (EAC), but it is not clear whether this reduction is modified by genotype, other host characteristics, or environmental factors. Furthermore, little is known about the association between H pylori and adenocarcinomas of the esophagogastric junction (EGJAC) or squamous cell carcinomas (ESCC). We sought to measure the association between H pylori infection and esophageal cancer and identify potential modifiers. METHODS In an Australian, population-based, case-control study, we compared the prevalence of H pylori seropositivity and single nucleotide polymorphisms in interleukin (IL)-1B (-31, -511) and tumor necrosis factor (TNF)-alpha (-308, -238) among 260 EAC, 298 EGJAC, and 208 ESCC patients and 1346 controls. To estimate relative risks, we calculated odds ratios (OR) and 95% confidence intervals (CI) using multivariable logistic regression in the entire sample and within strata of phenotypic and genotypic risk factors. RESULTS H pylori infection was associated with significantly reduced risks of EAC (OR, 0.45; 95% CI: 0.30-0.67) and EGJAC (OR, 0.41; 95% CI: 0.27-0.60) but not ESCC (OR, 1.04; 95% CI: 0.71-1.50). For each cancer subtype, risks were of similar magnitude across strata of reflux frequency and smoking status. We found no evidence that polymorphisms in IL-1B or TNF-alpha modified the association between H pylori and EAC or EGJAC. CONCLUSIONS H pylori infection is inversely associated with risks of EAC and EGJAC (but not ESCC); the reduction in risk is similar across subgroups of potential modifiers.

IRF4 Variants Have Age-Specific Effects on Nevus Count and Predispose to Melanoma

High melanocytic nevus count is a strong predictor of melanoma risk. A GWAS of nevus count in Australian adolescent twins identified an association of nevus count with the interferon regulatory factor 4 gene (IRF4 [p = 6 x 10(-9)]). There was a strong genotype-by-age interaction, which was replicated in independent UK samples of adolescents and adults. The rs12203592(*)T allele was associated with high nevus counts and high freckling scores in adolescents, but with low nevus counts and high freckling scores in adults. The rs12203592(*)T increased counts of flat (compound and junctional) nevi in Australian adolescent twins, but decreased counts of raised (intradermal) nevi. In combined analysis of melanoma case-control data from Australia, the UK, and Sweden, the rs12203592(*)C allele was associated with melanoma (odds ratio [OR] 1.15, p = 4 x 10(-3)), most significantly on the trunk (OR = 1.33, p = 2.5 x 10(-5)). The melanoma association was corroborated in a GWAS performed by the GenoMEL consortium for an adjacent SNP, rs872071 (rs872071(*)T: OR 1.14, p = 0.0035; excluding Australian, the UK, and Swedish samples typed at rs12203592: OR 1.08, p = 0.08).

ATG16L1 T300A Shows Strong Associations With Disease Subgroups in a Large Australian IBD Population: Further Support for Significant Disease Heterogeneity

OBJECTIVES:Crohns disease (CD) and ulcerative colitis (UC) are the two most common forms of inflammatory bowel disease (IBD), representing a significant health-care burden. A variant in the autophagy gene ATG16L1 (T300A) has been newly identified as a CD susceptibility locus by genome-wide association. Our aim was to assess the contribution of T300A in determining disease susceptibility and phenotype in two independent Australian IBD cohorts and explore the relationship between T300A and known CD risk factors (NOD2[ nucleotide-binding oligomerization domain containing 2] status and smoking).METHODS:In total, 669 CD and 543 UC cases, and 1,244 controls (study 1), 154 CD cases and 420 controls (study 2), and 702 unaffected parents from both groups were genotyped. We conducted case–control and family association analyses, and investigated relationships between T300A and disease subgroups and between NOD2 status and cigarette smoking (CD only).RESULTS:The strong association between CD and T300A was confirmed (P < 0.001), with a two-fold increase in disease risk associated with the GG genotype (odds ratio [OR] 1.96, 95% confidence interval [CI] 1.49–2.58), while ileal CD risk was almost three-fold (OR 2.73, CI 1.87–4.0). ATG16L1 and NOD2 were found to contribute independently to CD risk. A greater than seven-fold increased CD risk was observed for current smokers with a GG genotype (vs nonsmoking AA genotype; P < 0.001, OR 7.65, CI 4.21–13.91). A significant inverse association was found between T300A and UC (P= 0.002). This was strongest for patients with extensive, severe disease.CONCLUSIONS:We confirm the strong association between T300A and CD, specifically ileal subphenotype, and also report the first strong association of this variant with UC.