Genome-wide association study identifies three loci associated with melanoma risk

The delicate balance of melanoma immunotherapy

The strategy of immune modulation for the treatment of cancer is being refined with the introduction of multiple new therapeutic agents into the clinic. Melanoma is a disease where many of these agents have demonstrated efficacy. The mechanisms of action of these agents exploit the counter-regulatory mechanisms of the immune response. However, these agents are also associated with immune-related adverse events (IRAEs), which represent tissue-specific inflammatory responses. These IRAEs highlight the delicate balance of immunologic homeostasis and, with some interventions, may occur more frequently in patients who sustain a therapeutic response. This review will discuss melanoma immunogenicity and immunotherapy. Furthermore, the spectrum and distinction between a reversible immune adverse event and autoimmunity will be highlighted.

The genetics of melanoma: recent advances

Cutaneous malignant melanoma results from the interplay of genetic, host, and environmental factors. Genetic factors implicated in melanoma etiology include inherited high-, intermediate-, and low-risk susceptibility genes as well as numerous somatic mutations in melanoma tumors. CDKN2A is the major high-risk melanoma susceptibility gene identified to date. Recent identification of low-risk loci has been accomplished predominantly through genome-wide association studies. Whole-exome and whole-genome studies have identified numerous genes somatically altered in melanoma tumors and highlighted a higher mutation load in melanoma tumors compared with those in other cancers. This higher load is believed to be attributable to the preponderance of cytosine-to-thymine nucleotide substitutions as a result of UV radiation exposure. Technological advances, particularly next-generation sequencing, have increased the opportunities for germline and somatic gene discovery in melanoma and are opening up new avenues for understanding melanoma pathogenesis as well as leading to new opportunities for treatment.

Hypothesis driven single nucleotide polymorphism search (HyDn-SNP-S)

The advent of complete-genome genotyping across phenotype cohorts has provided a rich source of information for bioinformaticians. However the search for SNPs from this data is generally performed on a study-by-study case without any specific hypothesis of the location for SNPs that are predictive for the phenotype. We have designed a method whereby very large SNP lists (several gigabytes in size), combining several genotyping studies at once, can be sorted and traced back to their ultimate consequence in protein structure. Given a working hypothesis, researchers are able to easily search whole genome genotyping data for SNPs that link genetic locations to phenotypes. This allows a targeted search for correlations between phenotypes and potentially relevant systems, rather than utilizing statistical methods only. HyDn-SNP-S returns results that are less data dense, allowing more thorough analysis, including haplotype analysis. We have applied our method to correlate DNA polymerases to cancer phenotypes using four of the available cancer databases in dbGaP. Logistic regression and derived haplotype analysis indicates that ~80SNPs, previously overlooked, are statistically significant. Derived haplotypes from this work link POLL to breast cancer and POLG to prostate cancer with an increase in incidence of 3.01- and 9.6-fold, respectively. Molecular dynamics simulations on wild-type and one of the SNP mutants from the haplotype of POLL provide insights at the atomic level on the functional impact of this cancer related SNP. Furthermore, HyDn-SNP-S has been designed to allow application to any system. The program is available upon request from the authors.

Variants at the 9p21 locus and melanoma risk

Background

The influence of variants at the 9p21 locus on melanoma risk has been reported through investigation of CDKN2A variants through candidate gene approach as well as by genome wide association studies (GWAS).

Methods

In the present study we genotyped, 25 SNPs that tag 273 variants on chromosome 9p21 in 837 melanoma cases and 1154 controls from Spain. Ten SNPs were selected based on previous associations, reported in GWAS, with either melanocytic nevi or melanoma risk or both. The other 15 SNPs were selected to fine map the CDKN2A gene region.

Results

All the 10 variants selected from the GWAS showed statistically significant association with melanoma risk. Statistically significant association with melanoma risk was also observed for the carriers of the variant T-allele of rs3088440 (540 C>T) at the 3’ UTR of CDKN2A gene with an OR 1.52 (95% CI 1.14-2.04). Interaction analysis between risk associated polymorphisms and previously genotyped MC1R variants, in the present study, did not show any statistically significant association. Statistical significant association was observed for the interaction between phototypes and the rs10811629 (located in intron 5 of MTAP). The strongest association was observed between the homozygous carrier of the A–allele and phototype II with an OR of 15.93 (95% CI 5.34-47.54).

Conclusions

Our data confirmed the association of different variants at chromosome 9p21 with melanoma risk and we also found an association of a variant with skin phototypes.

Melanoma. Shall we move away from the sun and focus more on embryogenesis, body weight and longevity?

There are many observations regarding the behaviour of melanoma which points away from sunshine as the main cause of this tumour. Incidence data shows that the increase is mostly seen for thin melanomas which cannot be attributed to sun exposure but increasing screening over the last 20 years. Melanoma behaves in a similar fashion all over the world regarding age of onset, gender differences and histological subtypes. An excess of naevi is the strongest risk factor for melanoma and their appearance and involution throughout life, and the differences in naevus distribution according to gender is giving us a lot of clues about melanoma biology. Melanoma like all cancers is a complex disease with the involvement of many common and low penetrance genes many of them involved in pigmentation and naevogenesis but these only explain a very small portion of melanoma susceptibility. Genes involved in melanocyte differentiation early on in embryogenesis are also becoming relevant for melanoma initiation and progression. Reduced senescence and longevity as well as body weight and energy expenditure are also relevant for melanoma susceptibility. These observations with links between melanoma and non-sun related phenotypes as well as gene discoveries should help to assess the relative contribution of genetic and environmental factors in its causation.

Multiple neurofibromas as the presenting feature of familial atypical multiple malignant melanoma (FAMMM) syndrome

Mutations in the cyclin-dependent kinase inhibitor-2A (CDKN2A) gene have been associated with a number of malignancies, most notably cutaneous malignant melanoma (CMM). Mutations in this gene have also been associated with pancreatic cancer and breast cancer, as well as astrocytomas and other nervous system tumors (NST). Among NST, rare solitary internal neurofibromas have been reported, but multiple cutaneous neurofibromas have only been described in two families. In the first family, the affected individuals all carried a heterozygous G>C mutation at the splice acceptor site of intron 1 resulting in skipping of CDKN2A exon 2, while the affected individuals in the second family had a deletion that encompassed the whole CDKN2A/CDKN2B/ANRIL locus. We now report on a proposita presenting with multiple biopsy-proven cutaneous neurofibromas and a solitary spinal neurofibroma found to have a deletion of 14 nucleotides in exon 2 of CDKN2A, providing further evidence that p14, p16, and/or ANRIL are specifically involved in the pathogenesis of neurofibromas as a feature of the familial atypical multiple malignant melanoma spectrum.

Association between putative functional variants in the PSMB9 gene and risk of melanoma--re-analysis of published melanoma genome-wide association studies

To mine possibly hidden causal single-nucleotide polymorphisms (SNPs) of melanoma, we investigated the association of SNPs in 76 M/G1 transition genes with melanoma risk using our published genome-wide association study (GWAS) data set with 1804 melanoma cases and 1026 cancer-free controls. We found multiple SNPs with P < 0.01 and performed validation studies for 18 putative functional SNPs in PSMB9 in two other GWAS data sets. Two SNPs (rs1351383 and rs2127675) were associated with melanoma risk in the GenoMEL data set (P = 0.013 and 0.004, respectively), but failed in validation using the Australian data set. Genotype-phenotype analysis revealed these two SNPs were significantly correlated with mRNA expression level of PSMB9. Further experiments revealed that SNP rs2071480, which is in high LD with rs1351383 and rs2127675, may have a weak effect on the promoter activity of PSMB9. Taken together, our data suggested that functional variants in PSMB9 may contribute to melanoma susceptibility.

Robust methods for population stratification in genome wide association studies

Background

Genome-wide association studies can provide novel insights into diseases of interest, as well as to the responsiveness of an individual to specific treatments. In such studies, it is very important to correct for population stratification, which refers to allele frequency differences between cases and controls due to systematic ancestry differences. Population stratification can cause spurious associations if not adjusted properly. The principal component analysis (PCA) method has been relied upon as a highly useful methodology to adjust for population stratification in these types of large-scale studies. Recently, the linear mixed model (LMM) has also been proposed to account for family structure or cryptic relatedness. However, neither of these approaches may be optimal in properly correcting for sample structures in the presence of subject outliers.

Results

We propose to use robust PCA combined with k-medoids clustering to deal with population stratification. This approach can adjust for population stratification for both continuous and discrete populations with subject outliers, and it can be considered as an extension of the PCA method and the multidimensional scaling (MDS) method. Through simulation studies, we compare the performance of our proposed methods with several widely used stratification methods, including PCA and MDS. We show that subject outliers can greatly influence the analysis results from several existing methods, while our proposed robust population stratification methods perform very well for both discrete and admixed populations with subject outliers. We illustrate the new method using data from a rheumatoid arthritis study.

Conclusions

We demonstrate that subject outliers can greatly influence the analysis result in GWA studies, and propose robust methods for dealing with population stratification that outperform existing population stratification methods in the presence of subject outliers.

Colorful DNA polymorphisms in humans

In this review article we summarize current knowledge on how variation on the DNA level influences human pigmentation including color variation of iris, hair, and skin. We review recent progress in the field of human pigmentation genetics by focusing on the genes and DNA polymorphisms discovered to be involved in determining human pigmentation traits, their association with diseases particularly skin cancers, and their power to predict human eye, hair, and skin colors with potential utilization in forensic investigations.

NRAS and BRAF mutations in cutaneous melanoma and the association with MC1R genotype: findings from Spanish and Austrian populations

There is increasing epidemiologic and molecular evidence that cutaneous melanomas arise through multiple causal pathways. To further define the pathways to melanoma, we explored the relationship between germline and somatic mutations in a series of melanomas collected from 134 Spanish and 241 Austrian patients. Tumor samples were analyzed for melanocortin-1 receptor (MC1R) variants and mutations in the BRAF and NRAS genes. Detailed clinical data were systematically collected from patients. We found that NRAS-mutant melanomas were significantly more likely from older patients and BRAF-mutant melanomas were more frequent in melanomas from the trunk. We observed a nonsignificant association between germline MC1R status and somatic BRAF mutations in melanomas from trunk sites (odds ratio (OR) 1.8 (0.8-4.1), P=0.1), whereas we observed a significant inverse association between MC1R and BRAF for melanomas of the head and neck (OR 0.3 (0.1-0.8), P=0.02). This trend was observed in both the Spanish and Austrian populations.

A variant in FTO shows association with melanoma risk not due to BMI

We report the results of an association study of melanoma that is based on the genome-wide imputation of the genotypes of 1,353 cases and 3,566 controls of European origin conducted by the GenoMEL consortium. This revealed an association between several SNPs in intron 8 of the FTO gene, including rs16953002, which replicated using 12,313 cases and 55,667 controls of European ancestry from Europe, the USA and Australia (combined P = 3.6 × 10(-12), per-allele odds ratio for allele A = 1.16). In addition to identifying a new melanoma-susceptibility locus, this is to our knowledge the first study to identify and replicate an association with SNPs in FTO not related to body mass index (BMI). These SNPs are not in intron 1 (the BMI-related region) and exhibit no association with BMI. This suggests FTO's function may be broader than the existing paradigm that FTO variants influence multiple traits only through their associations with BMI and obesity.

Association between functional polymorphisms in genes involved in the MAPK signaling pathways and cutaneous melanoma risk

Genome-wide association studies (GWASs) have mainly focused on top significant single nucleotide polymorphisms (SNPs), most of which did not have clear biological functions but were just surrogates for unknown causal variants. Studying SNPs with modest association and putative functions in biologically plausible pathways has become one complementary approach to GWASs. To unravel the key roles of mitogen-activated protein kinase (MAPK) pathways in cutaneous melanoma (CM) risk, we re-evaluated the associations between 47 818 SNPs in 280 MAPK genes and CM risk using our published GWAS dataset with 1804 CM cases and 1026 controls. We initially found 105 SNPs with P ≤ 0.001, more than expected by chance, 26 of which were predicted to be putatively functional SNPs. The risk associations with 16 SNPs around DUSP14 (rs1051849) and a previous reported melanoma locus MAFF/PLA2G6 (proxy SNP rs4608623) were replicated in the GenoMEL dataset (P < 0.01) but failed in the Australian dataset. Meta-analysis showed that rs1051849 in the 3' untranslated regions of DUSP14 was associated with a reduced risk of melanoma (odds ratio = 0.89, 95% confidence interval: 0.82-0.96, P = 0.003, false discovery rate = 0.056). Further genotype-phenotype correlation analysis using the 90 HapMap lymphoblastoid cell lines from Caucasians showed significant correlations between two SNPs (rs1051849 and rs4608623) and messenger RNA expression levels of DUSP14 and MAFF (P = 0.025 and P = 0.010, respectively). Gene-based tests also revealed significant SNPs were over-represented in MAFF, PLA2G6, DUSP14 and other 16 genes. Our results suggest that functional SNPs in MAPK pathways may contribute to CM risk. Further studies are warranted to validate our findings.

MicroRNA-31 controls phenotypic modulation of human vascular smooth muscle cells by regulating its target gene cellular repressor of E1A-stimulated genes

Phenotypic modulation of vascular smooth muscle cells (VSMCs) plays a critical role in the pathogenesis of a variety of proliferative vascular diseases. The cellular repressor of E1A-stimulated genes (CREG) has been shown to play an important role in phenotypic modulation of VSMCs. However, the mechanism regulating CREG upstream signaling remains unclear. MicroRNAs (miRNAs) have recently been found to play a critical role in cell differentiation via target-gene regulation. This study aimed to identify a miRNA that binds directly to CREG, and may thus be involved in CREG-mediated VSMC phenotypic modulation. Computational analysis indicated that miR-31 bound to the CREG mRNA 3' untranslated region (3'-UTR). miR-31 was upregulated in quiescent differentiated VSMCs and downregulated in proliferative cells stimulated by platelet-derived growth factor and serum starvation, demonstrating a negative relationship with the VSMC differentiation marker genes, smooth muscle α-actin, calponin and CREG. Using gain-of-function and loss-of-function approaches, CREG and VSMC differentiation marker gene expression levels were shown to be suppressed by a miR-31 mimic, but increased by a miR-31 inhibitor at both protein and mRNA levels. Notably, miR-31 overexpression or inhibition affected luciferase expression driven by the CREG 3'-UTR containing the miR-31 binding site. Furthermore, miR-31-mediated VSMC phenotypic modulation was inhibited in CREG-knockdown human VSMCs. We also determined miR-31 levels in the serum of patients with coronary artery disease (CAD), with or without in stent restenosis and in healthy controls. miR-31 levels were higher in the serum of CAD patients with restenosis compared to CAD patients without restenosis and in healthy controls. In summary, these data demonstrate that miR-31 not only directly binds to its target gene CREG and modulates the VSMC phenotype through this interaction, but also can be an important biomarker in diseases involving VSMC phenotypic modulation. These novel findings may have extensive implications for the diagnosis and therapy of a variety of proliferative vascular diseases.

Melanoma as a subsequent neoplasm in adult survivors of childhood cancer: a report from the childhood cancer survivor study

BACKGROUND

Childhood cancer survivors have a sixfold increased risk of developing subsequent neoplasms when compared to the general population. We sought to describe the occurrence of melanoma as a subsequent neoplasm among adult survivors of childhood cancer.

PATIENTS AND METHODS

Among 14,358 5-year survivors of childhood cancer diagnosed between 1970 and 1986, we calculated the cumulative incidence, standardized incidence ratio (SIR), and absolute excess risk (AER) of subsequent melanoma. Potential risk factors were assessed using a cause-specific hazards model.

RESULTS

Fifty-seven melanomas (46 invasive, 2 ocular, and 9 in situ) occurred in 51 survivors. The median time to the development of melanoma was 21.0 years (range: 5.6-35.4 years) and the median age at melanoma was 32.3 years (range: 10.9-49.0 years). Initial cancer diagnoses included soft tissue and bone sarcoma (n = 15), leukemia (13), lymphoma (14), central nervous system malignancy (5), Wilms tumor (3), and neuroblastoma (1). The cumulative incidence of first subsequent melanoma at 35 years from initial cancer diagnosis was 0.55% [95% confidence interval (CI): 0.37-0.73]. The SIR of subsequent invasive malignant melanoma of the skin was 2.42 (95% CI: 1.77-3.23), and the AER was 0.10 (95% CI: 0.05-0.15) per 1,000 person-years. No statistically significant associations were found between melanoma risk and family history of cancer, demographic, or treatment-related factors.

CONCLUSION

Survivors of childhood cancer have an approximate 2.5-fold increased risk of melanoma. Early screening and prevention strategies are warranted.

Resequencing and clinical associations of the 9p21.3 region: a comprehensive investigation in the Framingham heart study

BACKGROUND

9p21.3 is among the most strongly replicated regions for cardiovascular disease. There are few reports of sequencing the associated 9p21.3 interval. We set out to sequence the 9p21.3 region followed by a comprehensive study of genetic associations with clinical and subclinical cardiovascular disease and its risk factors, as well as with copy number variation and gene expression, in the Framingham Heart Study (FHS).

METHODS AND RESULTS

We sequenced 281 individuals (94 with myocardial infarction, 94 with high coronary artery calcium levels, and 93 control subjects free of elevated coronary artery calcium or myocardial infarction), followed by genotyping and association in >7000 additional FHS individuals. We assessed genetic associations with clinical and subclinical cardiovascular disease, risk factor phenotypes, and gene expression levels of the protein-coding genes CDKN2A and CDKN2B and the noncoding gene ANRIL in freshly harvested leukocytes and platelets. Within this large sample, we found strong associations of 9p21.3 variants with increased risk for myocardial infarction, higher coronary artery calcium levels, and larger abdominal aorta diameters and no evidence for association with traditional cardiovascular disease risk factors. No common protein-coding variation, variants in splice donor or acceptor sites, or copy number variation events were observed. By contrast, strong associations were observed between genetic variants and gene expression, particularly for a short isoform of ANRIL and for CDKN2B.

CONCLUSIONS

Our thorough genomic characterization of 9p21.3 suggests common variants likely account for observed disease associations and provides further support for the hypothesis that complex regulatory variation affecting ANRIL and CDKN2B gene expression may contribute to increased risk for clinically apparent and subclinical coronary artery disease and aortic disease.

Serum 25-hydroxyvitamin D3 levels and vitamin D receptor variants in melanoma patients from the Mediterranean area of Barcelona

BACKGROUND

Serum 25-hydroxyvitamin D3 (Vitamin D) insufficiency and single-nucleotide polymorphisms (SNPs) on its receptor, Vitamin D receptor (VDR), have been reported to be involved in melanoma susceptibility in populations mostly from northern countries.

OBJECTIVE

To investigate 25-hydroxyvitamin D3 levels and VDR SNPs in melanoma patients from sunny area of Barcelona, two studies were carried out. The first study evaluated the levels of Vitamin D at time of melanoma diagnosis and the second one analyzed the association between VDR genetic variants and risk of having a high nevus number, the strongest phenotypic risk factor for melanoma.

METHODS

The levels of 25-hydroxyvitamin D3 in 81 melanoma patients at diagnosis were measured. In a second group of melanoma patients, including 150 with low and 113 with high nevus number, 11 VDR SNPs were analyzed for their association with nevus number.

RESULTS

In the first study, 68% of patients had insufficient levels of 25-hydroxyvitamin D3 (<25 ng/ml). Autumn-winter months and fair phototype were associated with 25-hydroxyvitamin D3 insufficiency; after multivariate analysis, season of sampling remained the only independent predictor of 25-hydroxyvitamin D3 levels. In the second study, VDR variant rs2189480 (P = 0.006) was associated with risk of high nevus number whereas rs2239179 (P = 0.044) and rs7975128 (P = 0.0005) were protective against high nevus number. After Bonferroni adjustment only rs7975128 remained significant. In stratified analysis, SNP rs7975128 was found protective against multiple melanomas (P = 0.021).

CONCLUSION

This study showed that even in Barcelona, a sunny Mediterranean area, 25-hydroxyvitamin D3 levels were sub-optimal in the majority of melanoma patients at diagnosis. The involvement of VDR in nevi and, in turn, in melanoma susceptibility has also been suggested. Larger studies are needed to confirm our findings.

Replication and predictive value of SNPs associated with melanoma and pigmentation traits in a Southern European case-control study

Background

Genetic association studies have revealed numerous polymorphisms conferring susceptibility to melanoma. We aimed to replicate previously discovered melanoma-associated single-nucleotide polymorphisms (SNPs) in a Greek case-control population, and examine their predictive value.

Methods

Based on a field synopsis of genetic variants of melanoma (MelGene), we genotyped 284 patients and 284 controls at 34 melanoma-associated SNPs of which 19 derived from GWAS. We tested each one of the 33 SNPs passing quality control for association with melanoma both with and without accounting for the presence of well-established phenotypic risk factors. We compared the risk allele frequencies between the Greek population and the HapMap CEU sample. Finally, we evaluated the predictive ability of the replicated SNPs.

Results

Risk allele frequencies were significantly lower compared to the HapMap CEU for eight SNPs (rs16891982 – SLC45A2, rs12203592 – IRF4, rs258322 – CDK10, rs1805007 – MC1R, rs1805008 - MC1R, rs910873 - PIGU, rs17305573- PIGU, and rs1885120 - MTAP) and higher for one SNP (rs6001027 – PLA2G6) indicating a different profile of genetic susceptibility in the studied population. Previously identified effect estimates modestly correlated with those found in our population (r = 0.72, P<0.0001). The strongest associations were observed for rs401681-T in CLPTM1L (odds ratio [OR] 1.60, 95% CI 1.22–2.10; P = 0.001), rs16891982-C in SCL45A2 (OR 0.51, 95% CI 0.34–0.76; P = 0.001), and rs1805007-T in MC1R (OR 4.38, 95% CI 2.03–9.43; P = 2×10⁻⁵). Nominally statistically significant associations were seen also for another 5 variants (rs258322-T in CDK10, rs1805005-T in MC1R, rs1885120-C in MYH7B, rs2218220-T in MTAP and rs4911442-G in the ASIP region). The addition of all SNPs with nominal significance to a clinical non-genetic model did not substantially improve melanoma risk prediction (AUC for clinical model 83.3% versus 83.9%, p = 0.66).

Conclusion

Overall, our study has validated genetic variants that are likely to contribute to melanoma susceptibility in the Greek population.

Common genetic variants in the 9p21 region and their associations with multiple tumours

BACKGROUND

The chromosome 9p21.3 region has been implicated in the pathogenesis of multiple cancers.

METHODS

We systematically examined up to 203 tagging SNPs of 22 genes on 9p21.3 (19.9-32.8 Mb) in eight case-control studies: thyroid cancer, endometrial cancer (EC), renal cell carcinoma, colorectal cancer (CRC), colorectal adenoma (CA), oesophageal squamous cell carcinoma (ESCC), gastric cardia adenocarcinoma and osteosarcoma (OS). We used logistic regression to perform single SNP analyses for each study separately, adjusting for study-specific covariates. We combined SNP results across studies by fixed-effect meta-analyses and a newly developed subset-based statistical approach (ASSET). Gene-based P-values were obtained by the minP method using the Adaptive Rank Truncated Product program. We adjusted for multiple comparisons by Bonferroni correction.

RESULTS

Rs3731239 in cyclin-dependent kinase inhibitors 2A (CDKN2A) was significantly associated with ESCC (P=7 × 10(-6)). The CDKN2A-ESCC association was further supported by gene-based analyses (Pgene=0.0001). In the meta-analyses by ASSET, four SNPs (rs3731239 in CDKN2A, rs615552 and rs573687 in CDKN2B and rs564398 in CDKN2BAS) showed significant associations with ESCC and EC (P<2.46 × 10(-4)). One SNP in MTAP (methylthioadenosine phosphorylase) (rs7023329) that was previously associated with melanoma and nevi in multiple genome-wide association studies was associated with CRC, CA and OS by ASSET (P=0.007).

CONCLUSION

Our data indicate that genetic variants in CDKN2A, and possibly nearby genes, may be associated with ESCC and several other tumours, further highlighting the importance of 9p21.3 genetic variants in carcinogenesis.

Assessment of the 9p21.3 locus in severity of coronary artery disease in the presence and absence of type 2 diabetes

BACKGROUND

The 9p21.3 locus is strongly associated with the risk of coronary artery disease (CAD) and with type 2 diabetes (T2D). We investigated the association of 9p21.3 variants with severity of CAD (defined by the number of vessel diseased [VD]) in the presence and absence of T2D.

METHODS

We tested 11 9p21.3-variants for association in a white Italian study (N = 2,908), and carried out replication in 2 independent white populations, a German study (N = 2,028) and a Canadian Study (N=950). SNP association and permutation analyses were conducted.

RESULTS

We identified two 9p21.3-variants, rs4977574 (P < 4×10(-4)) and rs2383207 (P < 1.5×10(-3)) that were associated with severity of CAD in subjects without T2D. Association of rs4977574 with severity of CAD was confirmed in the Canadian Study. Results from subgroup analysis among patients with T2D showed an interaction between rs10738610 and T2D with P = 4.82×10(-2). Further investigation showed that rs10738610 (P < 1.99×10(-2)) was found to be significantly associated with severity of CAD in subjects with T2D.

CONCLUSIONS

The 9p21.3 locus is significantly associated with severity of CAD. The number of associations of 9p21.3 variants with severity of CAD is variable to the presence and absence of T2D. In a CAD-susceptible region of 115 kb, there is only one variant associated with the severity of coronary vessel disease in the presence of type 2 diabetes.

Breed-predispositions to cancer in pedigree dogs

Cancer is a common problem in dogs and although all breeds of dog and crossbred dogs may be affected, it is notable that some breeds of pedigree dogs appear to be at increased risk of certain types of cancer suggesting underlying genetic predisposition to cancer susceptibility. Although the aetiology of most cancers is likely to be multifactorial, the limited genetic diversity seen in purebred dogs facilitates genetic linkage or association studies on relatively small populations as compared to humans, and by using newly developed resources, genome-wide association studies in dog breeds are proving to be a powerful tool for unravelling complex disorders. This paper will review the literature on canine breed susceptibility to histiocytic sarcoma, osteosarcoma, haemangiosarcoma, mast cell tumours, lymphoma, melanoma, and mammary tumours including the recent advances in knowledge through molecular genetic, cytogenetic, and genome wide association studies.

Cis-acting noncoding RNAs: friends and foes

In recent years, the number and types of known functional noncoding RNAs have increased considerably. A subset of both short- and long-sized species are known to be involved in the cis regulation of target genes located at or near the same genomic locus. Their expression is often coordinated with that of neighboring protein-coding genes, and in many cases, related transcripts can influence each other at one step or another during their biogenesis. Here, we review the current literature, summarizing the existing knowledge about mammalian cis-acting RNAs and their impact on physiological and disease states.

Update on the Epidemiology of Melanoma

Cutaneous malignant melanoma (CMM) has been increasing steadily in incidence over the past 30 years. Recent studies have explored associations between CMM and varying physiologic risk factors, such as nevi or hair and eye color, in addition to historical features such as a personal history of nonmelanoma skin cancer (NMSC), childhood cancers, Parkinson's Disease, hormone exposure and family history of CMM. Genome-wide association studies have also uncovered many genetic determinants of CMM risk. Ultimately, ultraviolet (UV) radiation exposure remains the most important modifiable risk factor for CMM. Organ transplant recipients, and nonsteroidal anti-inflammatory usage may also play a role. While risk factors are important to identify, effective campaigns to reduce the burden of disease through early detection and prevention are essential. We present detailed data regarding these facets of care for the CMM patient, and provide an update on the epidemiology of CMM.

Insights into the genetic architecture of early stage age-related macular degeneration: a genome-wide association study meta-analysis

Genetic factors explain a majority of risk variance for age-related macular degeneration (AMD). While genome-wide association studies (GWAS) for late AMD implicate genes in complement, inflammatory and lipid pathways, the genetic architecture of early AMD has been relatively under studied. We conducted a GWAS meta-analysis of early AMD, including 4,089 individuals with prevalent signs of early AMD (soft drusen and/or retinal pigment epithelial changes) and 20,453 individuals without these signs. For various published late AMD risk loci, we also compared effect sizes between early and late AMD using an additional 484 individuals with prevalent late AMD. GWAS meta-analysis confirmed previously reported association of variants at the complement factor H (CFH) (peak P = 1.5×10⁻³¹) and age-related maculopathy susceptibility 2 (ARMS2) (P = 4.3×10⁻²⁴) loci, and suggested Apolipoprotein E (ApoE) polymorphisms (rs2075650; P = 1.1×10⁻⁶) associated with early AMD. Other possible loci that did not reach GWAS significance included variants in the zinc finger protein gene GLI3 (rs2049622; P = 8.9×10⁻⁶) and upstream of GLI2 (rs6721654; P = 6.5×10⁻⁶), encoding retinal Sonic hedgehog signalling regulators, and in the tyrosinase (TYR) gene (rs621313; P = 3.5×10⁻⁶), involved in melanin biosynthesis. For a range of published, late AMD risk loci, estimated effect sizes were significantly lower for early than late AMD. This study confirms the involvement of multiple established AMD risk variants in early AMD, but suggests weaker genetic effects on the risk of early AMD relative to late AMD. Several biological processes were suggested to be potentially specific for early AMD, including pathways regulating RPE cell melanin content and signalling pathways potentially involved in retinal regeneration, generating hypotheses for further investigation.

Molecular signatures of basal cell carcinoma susceptibility and pathogenesis: a genomic approach

Gene expression profiling can be useful for phenotypic classification, investigation of functional pathways, and to facilitate the search for disease risk genes through the integration of transcriptional data with available genomic information. To enhance our understanding of the genetic and molecular basis of basal cell carcinoma (BCC) we performed global gene expression analysis to generate a disease-associated transcriptional profile. A gene signature composed of 331 differentially expressed genes (DEGs) was generated from comparing 4 lesional and 4 site-matched control samples using Affymetrix Human Genome U95A microarrays. Hierarchical clustering based on the obtained gene signature separated the samples into their corresponding phenotype. Pathway analysis identified several significantly overrepresented pathways including PPAR-γ signaling, TGF-β signaling and lipid metabolism, as well as confirmed the importance of SHH and p53 in the pathogenesis of BCC. Comparison of our microarray data with previous microarray studies revealed 13 DEGs overlapping in 3 studies. Several of these overlapping genes function in lipid metabolism or are components of the extracellular matrix, suggesting the importance of these and related pathways in BCC pathogenesis. BCC-associated DEGs were mapped to previously reported BCC susceptibility loci including 1p36, 1q42, 5p13.3, 5p15 and 12q11-13. Our analysis also revealed transcriptional 'hot spots' on chromosome 5 which help to confirm (5p13 and 5p15) and suggest novel (5q11.2-14.3, 5q22.1-23.3 and 5q31-35.3) disease susceptibility loci/regions. Integrating microarray analyses with reported genetic information helps to confirm and suggest novel disease susceptibility loci/regions. Identification of these specific genomic and/or transcriptional targets may lead to novel diagnostic and therapeutic modalities.

Influence of genetic variants in type I interferon genes on melanoma survival and therapy

Melanoma is an immunogenic tumor; however, the efficacy of immune-therapy shows large inter-individual variation with possible influence of background genetic variation. In this study we report the influence of genetic polymorphisms in the type I interferon gene cluster on chromosome 9p22 on melanoma survival. We genotyped 625 melanoma patients recruited in an oncology center in Germany for 44 polymorphisms located on chromosome 9p22 that were informative for 299 polymorphisms and spanned 15 type I interferon genes. Our results showed associations between time to metastasis/survival and two linked (r² = 0.76) polymorphisms, rs10964859 (C>G) and rs10964862 (C>A). The rs10964859 polymorphism was located at 3′UTR and rs10964862 was 9.40 Kb towards 5′UTR of IFNW1 gene. The carriers of the variant alleles of the rs10964859 and rs10964862 polymorphisms were associated with a reduced disease-free survival. The validation of data in an independent group of 710 patients from Spain showed that the direction of the effect was similar. Stratification based on therapy showed that the adverse effect on metastasis development was statistically significant in the patients from Spain who did not receive any treatment and were homozygous for variant allele of rs10964862 (HR = 2.52, 95% CI 1.07–5.90; P = 0.03). Patients homozygous for rs10964859 (HR = 2.01, 95% CI 1.17–3.44; P = 0.01) and rs10964862 (HR 1.84, 95%CI 1.03–3.27, P = 0.04) were associated to increased risk of death following metastasis. GTCGACAA haplotype, found in 8.8% of the patients, was associated with an increased risk of death (HR 1.94, 95%CI 1.16–3.26, P = 0.01). In conclusion, our results identified genetic variants in interferon genes that influence melanoma progression and survival with modulation of effect due to treatment status.

Genetic variants at 6p21.1 and 7p15.3 are associated with risk of multiple cancers in Han Chinese

Cancer susceptibility loci identified in reported genome-wide association studies (GWAS) are often tumor-specific; however, evidence of pleiotropy of some genes/loci has also been observed and biologically plausible. We hypothesized that there are important regions in the genome harboring genetic variants associated with risk of multiple types of cancer. In the current study, we attempted to map genetic variants that have consistent effects on risk of multiple cancers using our existing genome-wide scan data of lung cancer, noncardia gastric cancer, and esophageal squamous-cell carcinoma with overall 5,368 cases and 4,006 controls (GWAS stage), followed by a further evaluation in additional 9,001 cases with one of these cancer types and 11,436 controls (replication stage). Five variants satisfying the criteria of pleiotropy with p values from 1.10 × 10(-8) to 8.96 × 10(-6) for genome-wide scans of three cancer types were further evaluated in the replication stage. We found consistent associations of rs2494938 at 6p21.1 and rs2285947 at 7p15.3 with these three cancers in both GWAS and replication stages. In combined samples of GWAS and replication stages, the minor alleles of rs2494938 and rs2285947 were significantly associated with an increased risk of the cancers (odds ratio [OR] = 1.15, 95% confidence interval [CI], 1.10-1.19 and OR = 1.17, 95% CI, 1.12-1.21), with the p values being 1.20 × 10(-12) and 1.26 × 10(-16), respectively, which are at a genome-wide significance level. Our findings highlight the potential importance of variants at 6p21.1 and 7p15.3 in the susceptibility to multiple cancers.

Genetic investigation into the endophenotypic status of central corneal thickness and optic disc parameters in relation to open-angle glaucoma

PURPOSE

To ascertain if single nucleotide polymorphisms (SNPs) involved in the determination of central corneal thickness, optic disc area, and vertical cup-to-disc ratio (VCDR) also are associated with open-angle glaucoma (OAG).

DESIGN

Retrospective case-control genetic association study.

METHODS

A total of 16 SNPs associated with central corneal thickness, optic disc area, and VCDR were genotyped in 876 OAG cases and 883 normal controls. To determine if the SNPs were also correlated with OAG severity, the cohort was stratified into advanced OAG (n = 326) and nonadvanced OAG (n = 550). Both the cases and controls were of European descent and were recruited from within Australia.

RESULTS

Two VCDR SNPs were found to be significantly associated with OAG after correction for multiple testing. The 2 SNPs were rs10483727, found adjacent to the SIX1 gene (P = 6.2 × 10(-06); odds ratio, 1.38; 95% confidence interval, 1.20 to 1.59), and rs1063192, found within the CDKN2B gene (P = 2.2 × 10(-05); odds ratio, 0.74; 95% confidence interval, 0.64 to 0.85). The CDKN2B variant rs1063192 also was found to be associated more strongly with advanced OAG.

CONCLUSIONS

The findings from this study indicate that variants influencing VCDR are also risk alleles for OAG in our Australian cohort of European descent. The identification of SIX1 and CDKN2B as susceptibility loci will assist in understanding the pathologic mechanisms involved in the development of OAG.

Variants affecting exon skipping contribute to complex traits

DNA variants that affect alternative splicing and the relative quantities of different gene transcripts have been shown to be risk alleles for some Mendelian diseases. However, for complex traits characterized by a low odds ratio for any single contributing variant, very few studies have investigated the contribution of splicing variants. The overarching goal of this study is to discover and characterize the role that variants affecting alternative splicing may play in the genetic etiology of complex traits, which include a significant number of the common human diseases. Specifically, we hypothesize that single nucleotide polymorphisms (SNPs) in splicing regulatory elements can be characterized in silico to identify variants affecting splicing, and that these variants may contribute to the etiology of complex diseases as well as the inter-individual variability in the ratios of alternative transcripts. We leverage high-throughput expression profiling to 1) experimentally validate our in silico predictions of skipped exons and 2) characterize the molecular role of intronic genetic variations in alternative splicing events in the context of complex human traits and diseases. We propose that intronic SNPs play a role as genetic regulators within splicing regulatory elements and show that their associated exon skipping events can affect protein domains and structure. We find that SNPs we would predict to affect exon skipping are enriched among the set of SNPs reported to be associated with complex human traits.

Alternative splicing is a common eukaryotic cellular mechanism that allows for the production of multiple proteins from one gene and occurs in 40%–90% of all human genes. Alternative splicing has been shown to be important for many critical biological processes, including development, evolution, and even psychological behavior. Additionally, alternative splicing has been associated with 15%–50% of human genetic diseases, including breast cancer; however, the precise mechanism by which genetic variations regulate this process remains to be fully elucidated. In this study, we develop an integrative approach that utilizes sequence-based analysis and genome-wide expression profiling to identify genetic variations that may affect alternative splicing. We also evaluate their enrichment among established disease-associated variations. Our study provides insights into the functionality of these variations and emphasizes their importance for complex human traits and diseases.

Long non-coding RNAs in cancer progression

Recent large-scale transcriptome analyses have revealed that transcription is spread throughout the mammalian genomes, yielding large numbers of transcripts, including long non-coding RNAs (lncRNAs) with little or no protein-coding capacity. Dozens of lncRNAs have been identified as biologically significant. In many cases, lncRNAs act as key molecules in the regulation of processes such as chromatin remodeling, transcription, and post-transcriptional processing. Several lncRNAs (e.g., MALAT1, HOTAIR, and ANRIL) are associated with human diseases, including cancer. Those lncRNAs associated with cancer are often aberrantly expressed. Although the underlying molecular mechanisms by which lncRNAs regulate cancer development are unclear, recent studies have revealed that such aberrant expression of lncRNAs affects the progression of cancers. In this review, we highlight recent findings regarding the roles of lncRNAs in cancer biology.

Network-based SNP meta-analysis identifies joint and disjoint genetic features across common human diseases

BACKGROUND

Genome-wide association studies (GWAS) have provided a large set of genetic loci influencing the risk for many common diseases. Association studies typically analyze one specific trait in single populations in an isolated fashion without taking into account the potential phenotypic and genetic correlation between traits. However, GWA data can be efficiently used to identify overlapping loci with analogous or contrasting effects on different diseases.

RESULTS

Here, we describe a new approach to systematically prioritize and interpret available GWA data. We focus on the analysis of joint and disjoint genetic determinants across diseases. Using network analysis, we show that variant-based approaches are superior to locus-based analyses. In addition, we provide a prioritization of disease loci based on network properties and discuss the roles of hub loci across several diseases. We demonstrate that, in general, agonistic associations appear to reflect current disease classifications, and present the potential use of effect sizes in refining and revising these agonistic signals. We further identify potential branching points in disease etiologies based on antagonistic variants and describe plausible small-scale models of the underlying molecular switches.

CONCLUSIONS

The observation that a surprisingly high fraction (>15%) of the SNPs considered in our study are associated both agonistically and antagonistically with related as well as unrelated disorders indicates that the molecular mechanisms influencing causes and progress of human diseases are in part interrelated. Genetic overlaps between two diseases also suggest the importance of the affected entities in the specific pathogenic pathways and should be investigated further.

SNP in the genome-wide association study hotspot on chromosome 9p21 confers susceptibility to diabetic nephropathy in type 1 diabetes

AIMS/HYPOTHESIS

Parental type 2 diabetes mellitus increases the risk of diabetic nephropathy in offspring with type 1 diabetes mellitus. Several single nucleotide polymorphisms (SNPs) that predispose to type 2 diabetes mellitus have recently been identified. It is, however, not known whether such SNPs also confer susceptibility to diabetic nephropathy in patients with type 1 diabetes mellitus.

METHODS

We genotyped nine SNPs associated with type 2 diabetes mellitus in genome-wide association studies in the Finnish population, and tested for their association with diabetic nephropathy as well as with severe retinopathy and cardiovascular disease in 2,963 patients with type 1 diabetes mellitus. Replication of significant SNPs was sought in 2,980 patients from three other cohorts.

RESULTS

In the discovery cohort, rs10811661 near gene CDKN2A/B was associated with diabetic nephropathy. The association remained after robust Bonferroni correction for the total number of tests performed in this study (OR 1.33 [95% CI 1.14, 1.56], p = 0.00045, p (36tests) = 0.016). In the meta-analysis, the combined result for diabetic nephropathy was significant, with a fixed effects p value of 0.011 (OR 1.15 [95% CI 1.02, 1.29]). The association was particularly strong when patients with end-stage renal disease were compared with controls (OR 1.35 [95% CI 1.13, 1.60], p = 0.00038). The same SNP was also associated with severe retinopathy (OR 1.37 [95% CI 1.10, 1.69] p = 0.0040), but the association did not remain after Bonferroni correction (p (36tests) = 0.14). None of the other selected SNPs was associated with nephropathy, severe retinopathy or cardiovascular disease.

CONCLUSIONS/INTERPRETATION

A SNP predisposing to type 2 diabetes mellitus, rs10811661 near CDKN2A/B, is associated with diabetic nephropathy in patients with type 1 diabetes mellitus.

Franklin H. Epstein Lecture. Both ends of the leash--the human links to good dogs with bad genes

For nearly 350 years, veterinary medicine and human medicine have been separate entities, with one geared toward the diagnosis and treatment in animals and the other toward parallel goals in the owners. However, that model no longer fits, since research on diseases of humans and companion animals has coalesced.^– The catalyst for this union has been the completion of the human genome sequence, coupled with draft sequence assemblies of genomes for companion animals.^, Here, we summarize the critical events in canine genetics and genomics that have led to this development, review major applications in canine health that will be of interest to human caregivers, and discuss expectations for the future.

Benefits of oral Polypodium Leucotomos extract in MM high-risk patients

BACKGROUND

UV radiation and the presence of melanocytic nevi are the main risk factors of sporadic melanoma (MM). Protection of skin by an oral photoprotective agent would have substantial benefits.

OBJECTIVE

We investigated the possible role of an oral Polypodium leucotomos (PL) extract to improve systemic photoprotection in patients at risk of skin cancer analyzing the ability to decrease UV-induced erythema. We also studied the interaction among MC1R polymorphisms and CDKN2A status with the minimal erythematous dose (MED) and their influence in the response after oral PL.

METHODS

A total of 61 patients (25 with familial and/or multiple MM, 20 with sporadic MM and 16 with atypical mole syndrome without history of MM) were exposed to varying doses of artificial UVB radiation without and after oral administration of a total dose of 1080 mg of PL.

RESULTS

Oral PL treatment significantly increased the MED mean in all group patients (0.123 to 0.161 J/cm(2) , p<0.05). Although not significant, we noticed a stronger effect of PL on the MED of patients with familial MM compared to those with MM (U=273, p=0.06). Among the patients with familial MM, those exhibiting a mutated CDKN2A and/or polymorphisms in MC1R had the bigger differences in response to treatment with PL.

LIMITATIONS

Reduced number of patients. No control population.

CONCLUSIONS

Administration of PL leads to a significant reduction of sensitivity to UVR (p<0.05) in all patients. Dark-eye patients and patients with higher UVR sensibility (lower basal MED) would be the most benefited from oral PL treatment.

The three M's: melanoma, microphthalmia-associated transcription factor and microRNA

Studies examining intratumor heterogeneity have indicated that several cancer types, including melanoma, can display phenotypic plasticity, corresponding to their capacity to undergo transient reversible cellular changes. Conceptual models constructed to explain the process of cancer propagation differ in their treatment of intratumor heterogeneity. Recent observations of reversible phenotypic heterogeneity in melanoma have led to the proposal of a novel 'phenotypic plasticity' model of cancer propagation. Microphthalmia-associated transcription factor (MITF), the melanocyte 'lineage-specific' transcription factor, has emerged as one of the central players in melanoma phenotypic plasticity. Here we discuss the conceptual models suggested to explain the relations between MITF and melanoma plasticity, in addition to the complex regulatory roles that MITF plays in melanocytes and melanoma development. Finally, we provide an in-depth literature survey of microRNAs (miRNAs) involved in MITF activity, melanoma propagation and metastasis, in addition to their potential use as agents of personalized therapy.

The chromosome 9p21.3 coronary heart disease risk allele is associated with altered gene expression in normal heart and vascular tissues

Genome-wide association studies have identified a coronary artery disease (CAD) risk locus in a non-coding region at 9p21.3, the nearest genes being CDKN2A and CDKN2B. To understand the pathways by which this locus might influence CAD susceptibility, we investigated associations between the 9p21.3 risk genotype and global gene expression in heart tissue from donors with no diagnosed heart disease (n = 108, predominant cause of death, cerebral vascular accident) and in carotid plaque (n = 106), aorta (n = 104) and mammary artery (n = 88) tissues from heart valve and carotid endarterectomy patients. Genotyping was performed with Taqman assays and Illumina arrays, and gene expression profiles generated with Affymetrix microarrays. Associations were analyzed with an additive genetic model. In heart tissue, 46 genes were putatively altered in association with the 9p21.3 risk allele (70% down-regulated, fold-change >1.1 per allele, p<0.05 adjusted for age, gender, ethnicity, cause of death). These genes were enriched for biomarkers of myocardial infarction (p = 1.53×10⁻⁹), response to wounding (p = 2.65×10⁻¹⁰) and inflammatory processes (p<1.97×10⁻⁷). Among the top 10 most down-regulated genes, 7 genes shared a set of transcription factor binding sites within conserved promoter regions (p<1.14×10⁻⁵), suggesting they may be co-regulated. Canonical pathway modelling of the most differentially expressed transcripts across all tissues (154 genes, 60% down-regulated, fold-change >1.1 per allele, p<0.01) showed that 75% of the genes could be transcriptionally regulated through the cell cycle G1 phase progression pathway (p<1.08×10⁻²⁵⁸), in which CDKN2A and CDKN2B play a regulatory role. These data suggest that the cell cycle G1 phase progression pathway is activated in individuals with the 9p21.3 risk allele. This may contribute to a proliferative phenotype that promotes adverse cardiac hypertrophy and vascular remodeling, leading to an increased CAD risk.

Variants at chromosome 20 (ASIP locus) and melanoma risk

Agouti signaling protein (ASIP) locus on chromosome 20q11 is implicated, as shown by genome-wide association studies, in phenotype variation and melanoma risk. We genotyped 837 melanoma cases and 1,154 controls for 21 single nucleotide polymorphisms (SNPs) informative for 495 polymorphisms at the locus. Our data showed an increased risk of melanoma (odds ratio [OR] 1.27, 95% confidence interval [95% CI] 1.03-1.57) in carriers of the rs4911414 variant, located 120 kb upstream of ASIP. The main effect of rs4911414, as reported previously, was in tandem with a 10 kb adjacent polymorphism rs1015362; two constituted risk-associated haplotype/diplotype. Except for rs1015363, none of the 12 tagging SNPs, genotyped to cover 239.9 kb region with polymorphisms linked to rs4911414 and rs1015362, were associated with melanoma. Our data confirmed a previous association of melanoma risk (OR 1.82, 95% CI 1.37-2.41) with rs4911442, located in intron 5 of the nuclear receptor coactivator 6 (NCOA6) gene. The rs910871, one of the six variants, genotyped to cover NCOA6, showed an association with melanoma risk (OR 1.33, 95% CI 1.04-1.70). Both, rs4911442 and rs910871 were in moderate linkage with a, previously reported, risk-associated rs910873 polymorphism. A haplotype from the variants within NCOA6 showed an association with risk of melanoma (OR 1.49, 95% CI 1.17-1.88). Interaction between risk-associated polymorphisms and previously genotyped melanocortin receptor 1 (MC1R) variants, in our study, was not statistically significant. Nevertheless, the carriers of the variant alleles over the background of MC1R variants were at a higher risk than the carriers not enriched for MC1R variants. Our data confirmed the association of different variants at chromosome 20q11 with melanoma risk.

Dysplastic nevi

Dysplastic nevi have been a subject of much debate since their original description in 1978. Although some question the biological potential of dysplastic nevi themselves, several studies have shown that their presence confers substantial risk for melanoma. In addition to predisposing patients to melanoma, dysplastic nevi have been shown to harbor genetic mutations, indicating their position on a continuum between banal nevi and melanomas. Dysplastic nevi are also clinically relevant as mimickers of melanoma, and can be challenging diagnostically. This article reviews the history, epidemiology, biology and genetics, clinical features, histopathologic features, and management guidelines for patients with these lesions.

Inherited variants in the MC1R gene and survival from cutaneous melanoma: a BioGenoMEL study

Inherited MC1R variants modulate MITF transcription factor signaling, which in turn affects tumor cell proliferation, apoptosis, and DNA repair. The aim of this BioGenoMEL collaborative study in 10 melanoma cohorts was to test the hypothesis that inherited variants thereby moderate survival expectation. A survival analysis in the largest cohort (Leeds) was carried out adjusting for factors known to impact on survival. The results were then compared with data from nine smaller cohorts. The absence of any consensus MC1R alleles was associated with a significantly lower risk of death in the Leeds set (HR, 0.64; 95% CI, 0.46-0.89) and overall in the 10 data sets (HR, 0.78; 95% CI, 0.65-0.94) with some support from the nine smaller data sets considered together (HR, 0.83; 95% CI, 0.67-1.04). The data are suggestive of a survival benefit for inherited MC1R variants in melanoma patients.

Glaucoma risk alleles at CDKN2B-AS1 are associated with lower intraocular pressure, normal-tension glaucoma, and advanced glaucoma

PURPOSE

Genetic variation at the 9p21 locus encompassing the CDKN2B-AS1, CDKN2A, and CDKN2B genes has been associated with primary open-angle glaucoma (POAG) in several independent studies. This study aimed to dissect the association further and to determine genotype-phenotype correlations between genetic variation at this locus and a range of glaucoma-related traits in a large cohort of POAG patients.

DESIGN

Comparative case series and case-control study.

PARTICIPANTS

One thousand four hundred thirty-two POAG patients and 595 unaffected controls recruited from 2 population-based and 2 cross-sectional studies.

METHODS

Each patient was genotyped at 9 single nucleotide polymorphisms (SNPs) previously associated with POAG at the 9p21 locus. Each SNP was assessed for association with each outcome measure using linear regression under an additive genetic model. Associated traits were explored further including adjustment for relevant covariates. Highest recorded intraocular pressure (IOP) also was analyzed both with and without correction for central corneal thickness (CCT) and was dichotomized into high-tension glaucoma and normal-tension glaucoma (NTG).

MAIN OUTCOME MEASURES

Intraocular pressure and vertical cup-to-disc ratio (VCDR).

RESULTS

Glaucoma risk alleles at 9p21, particularly, rs7049105 and rs10120688, were associated with the presence of both NTG and advanced POAG. The SNP rs10120688 was associated with greater VCDR after adjustment for covariates (P = 0.003; β = 0.016; standard error, 0.006). In addition, multiple SNPs in the region were associated with reduced IOP, before and after adjustment for CCT. The SNP most significantly associated with IOP was also rs10120688 (P = 0.001; β = -2.135; standard error, 0.634) after adjustment for covariates under an additive model. In a comparison of high-tension versus low-tension glaucoma, this SNP was also the most significantly associated, particularly when IOP was corrected for CCT before classification of the type of glaucoma (P = 0.0009; odds ratio, 0.63; 95% confidence interval, 0.48-0.83).

CONCLUSIONS

Patients with POAG carrying the glaucoma risk alleles at the 9p21 locus have larger VCDR and lower IOP than POAG patients without these alleles. Carriers of these alleles seem to be predisposed to POAG developing at lower IOP levels and exhibit stronger associations with NTG and advanced glaucoma phenotypes. This may be of relevance when setting target pressures in patients carrying these risk alleles.