Functional Variants in Notch Pathway Genes NCOR2, NCSTN, and MAML2 Predict Survival of Patients with Cutaneous Melanoma.

BACKGROUND: The Notch signaling pathway is constitutively activated in human cutaneous melanoma to promote growth and aggressive metastatic potential of primary melanoma cells. Therefore, genetic variants in Notch pathway genes may affect the prognosis of cutaneous melanoma patients. METHODS: We identified 6,256 SNPs in 48 Notch genes in 858 cutaneous melanoma patients included in a previously published cutaneous melanoma genome-wide association study dataset. Multivariate and stepwise Cox proportional hazards regression and false-positive report probability corrections were performed to evaluate associations between putative functional SNPs and cutaneous melanoma disease-specific survival. Receiver operating characteristic curve was constructed, and area under the curve was used to assess the classification performance of the model. RESULTS: Four putative functional SNPs of Notch pathway genes had independent and joint predictive roles in survival of cutaneous melanoma patients. The most significant variant was NCOR2 rs2342924 T>C (adjusted HR, 2.71; 95% confidence interval, 1.73-4.23; Ptrend = 9.62 × 10(-7)), followed by NCSTN rs1124379 G>A, NCOR2 rs10846684 G>A, and MAML2 rs7953425 G>A (Ptrend = 0.005, 0.005, and 0.013, respectively). The receiver operating characteristic analysis revealed that area under the curve was significantly increased after adding the combined unfavorable genotype score to the model containing the known clinicopathologic factors. CONCLUSIONS: Our results suggest that SNPs in Notch pathway genes may be predictors of cutaneous melanoma disease-specific survival. IMPACT: Our discovery offers a translational potential for using genetic variants in Notch pathway genes as a genotype score of biomarkers for developing an improved prognostic assessment and personalized management of cutaneous melanoma patients.

Integrated pathway and epistasis analysis reveals interactive effect of genetic variants at TERF1 and AFAP1L2 loci on melanoma risk.

Genome-wide association studies (GWASs) have characterized 13 loci associated with melanoma, which only account for a small part of melanoma risk. To identify new genes with too small an effect to be detected individually but which collectively influence melanoma risk and/or show interactive effects, we used a two-step analysis strategy including pathway analysis of genome-wide SNP data, in a first step, and epistasis analysis within significant pathways, in a second step. Pathway analysis, using the gene-set enrichment analysis (GSEA) approach and the gene ontology (GO) database, was applied to the outcomes of MELARISK (3,976 subjects) and MDACC (2,827 subjects) GWASs. Cross-gene SNP-SNP interaction analysis within melanoma-associated GOs was performed using the INTERSNP software. Five GO categories were significantly enriched in genes associated with melanoma (false discovery rate ≤ 5% in both studies): response to light stimulus, regulation of mitotic cell cycle, induction of programmed cell death, cytokine activity and oxidative phosphorylation. Epistasis analysis, within each of the five significant GOs, showed significant evidence for interaction for one SNP pair at TERF1 and AFAP1L2 loci (pmeta-int  = 2.0 × 10(-7) , which met both the pathway and overall multiple-testing corrected thresholds that are equal to 9.8 × 10(-7) and 2.0 × 10(-7) , respectively) and suggestive evidence for another pair involving correlated SNPs at the same loci (pmeta-int  = 3.6 × 10(-6) ). This interaction has important biological relevance given the key role of TERF1 in telomere biology and the reported physical interaction between TERF1 and AFAP1L2 proteins. This finding brings a novel piece of evidence for the emerging role of telomere dysfunction into melanoma development.

Common variants at the MHC locus and at chromosome 16q24.1 predispose to Barrett's esophagus

Barrett's esophagus is an increasingly common disease that is strongly associated with reflux of stomach acid and usually a hiatus hernia, and it strongly predisposes to esophageal adenocarcinoma (EAC), a tumor with a very poor prognosis. We report the first genome-wide association study on Barrett's esophagus, comprising 1,852 UK cases and 5,172 UK controls in the discovery stage and 5,986 cases and 12,825 controls in the replication stage. Variants at two loci were associated with disease risk: chromosome 6p21, rs9257809 (P(combined) = 4.09 × 10(-9); odds ratio (OR) = 1.21, 95% confidence interval (CI) =1.13-1.28), within the major histocompatibility complex locus, and chromosome 16q24, rs9936833 (P(combined) = 2.74 × 10(-10); OR = 1.14, 95% CI = 1.10-1.19), for which the closest protein-coding gene is FOXF1, which is implicated in esophageal development and structure. We found evidence that many common variants of small effect contribute to genetic susceptibility to Barrett's esophagus and that SNP alleles predisposing to obesity also increase risk for Barrett's esophagus.

Seven new loci associated with age-related macular degeneration

Age-related macular degeneration (AMD) is a common cause of blindness in older individuals. To accelerate the understanding of AMD biology and help design new therapies, we executed a collaborative genome-wide association study, including >17,100 advanced AMD cases and >60,000 controls of European and Asian ancestry. We identified 19 loci associated at P <5 × 10(-8). These loci show enrichment for genes involved in the regulation of complement activity, lipid metabolism, extracellular matrix remodeling and angiogenesis. Our results include seven loci with associations reaching P <5 × 10(-8) for the first time, near the genes COL8A1-FILIP1L, IER3-DDR1, SLC16A8, TGFBR1, RAD51B, ADAMTS9 and B3GALTL. A genetic risk score combining SNP genotypes from all loci showed similar ability to distinguish cases and controls in all samples examined. Our findings provide new directions for biological, genetic and therapeutic studies of AMD.

Five authors reply

We read with interest the comments offered by Drs. Hughes and Bradley (1) on our systematic review (2). Four single nucleotide polymorphisms (SNPs), rs9332739 and rs547154 in the complement component 2 gene (C2) and rs4151667 and rs641153 in the complement factor B gene (CFB), were pooled. Hughes and Bradley point out that we omitted the most common variant, rs12614. In fact, rs12614 is in high linkage disequilibrium (LD) with rs641153, which was included, and the major allele of both of these SNPs is in the range of 90% (population code, CEU, in the International HapMap Project (http://hapmap.ncbi.nlm.nih.gov/)). Moreover, our review was initiated in September 2010, at which point only 4 studies had published associations with rs12614, whereas 14 studies (n = 11,378) were available for rs641153. While it is true that both SNPs are better analyzed as a haplotype, these data were simply not available for pooling.
Hughes and Bradley also point out that we obtained and pooled new data that were not previously published. While it is recommended that contact with authors be completed as part of a comprehensive meta-analysis, we acknowledge that these additional data were not previously published and peer reviewed and, hence, do not have the same level of transparency. However, given that sample collections often increase over time and that the instrumentation for genotyping is continually improving, we thought that it would be advantageous to use the most recent information; this is a subjective decision.
We also agree that the allele frequencies given by Kaur et al. (3) were exactly opposite to those expected and were suggestive of strand flipping. However, we specifically queried this with the lead author on 2 separate occasions and were assured it was not.
Hughes and Bradley do make an interesting suggestion that SNPs in high LD should be used as a gauge of genotyping quality in HuGE reviews. This is an interesting idea but difficult to put into practice as the r2 parameter they propose as a measure of LD has some unusual properties. Although r2 is a measure of LD, it is also linked to the allele frequency; even small differences in allele frequencies between 2 linked SNPs can reduce the r2 dramatically. Wray (4) explored these effects and found that, at a baseline allele frequency of 10%, even a difference in allele frequency between 2 SNPs as small as 2% can drop the r2 value below 0.8. This degree of allele frequency difference is consistent with what could be expected for sampling error. Furthermore, when we look at 2 linked dialleleic SNPs, giving 4 possible haplotypes, the absence of 1 haplotype dramatically reduces r2, despite the 2 loci being in high LD as measured by D'. In fact, this is the situation for rs12614 and rs641153, where the low frequency of 1 haplotype means that the r2 is 0.01 but the D' is 1.
Hughes and Bradley also suggest consideration of genotype call rate restrictions as an inclusion criterion for metaanalysis. This would be more appropriate when focusing on genetic variants per se, as considered within the context of a genome-wide association study or other specific genetic analysis where large numbers of SNPs are evaluated (5).
The concerns raised by Hughes and Bradley reflect the limited ability of a meta-analysis based on summary data to tease out inconsistencies best identified at the individual level. We agree that SNPs in LD should be evaluated, but this will not necessarily be straightforward. A move to make genetic data sets publicly available, as in the Database of Genotypes and Phenotypes (http://www.ncbi.nlm.nih.gov/ gap), is a step in the right direction for greater transparency.

A Variant in LDLR is Associated with Abdominal Aortic Aneurysm

Background—Abdominal aortic aneurysm (AAA) is a common cardiovascular disease among older people and demonstrates significant heritability. In contrast to similar complex diseases, relatively few genetic associations with AAA have been confirmed. We reanalysed our genome-wide study and carried through to replication suggestive discovery associations at a lower level of significance.

Methods and Results—A genome-wide association study was conducted using 1,830 cases from the UK, New Zealand and Australia with infra-renal aorta diameter =30mm or ruptured AAA and 5,435 unscreened controls from the 1958 Birth Cohort and National Blood Service cohort from the Wellcome Trust Case Control Consortium. Eight suggestive associations with P<1x10-4 were carried through to in silico replication in 1,292 AAA cases and 30,503 controls. One SNP associated with P<0.05 after Bonferroni correction in the in silico study underwent further replication (706 AAA cases and 1,063 controls from the UK, 507 AAA cases and 199 controls from Denmark and 885 AAA cases and 1,000 controls from New Zealand). Low density lipoprotein receptor (LDLR) rs6511720 A, was significantly associated overall and in three of five individual replication studies. The full study showed an association that reached genome-wide significance (odds ratio 0.76; 95% confidence interval 0.70 to 0.83; P=2.08x10-10).

Conclusions—LDLR rs6511720 is associated with abdominal aortic aneurysm. This finding is consistent with established effects of this variant on coronary artery disease. Shared aetiological pathways with other cardiovascular diseases may present novel opportunities for preventative and therapeutic strategies for AAA.

Chromosome 2q31. 1 associates with ESRD in women with type 1 diabetes

Sex and genetic variation influence the risk of developing diabetic nephropathy and ESRD in patients with type 1 diabetes. We performed a genome-wide association study in a cohort of 3652 patients from the Finnish Diabetic Nephropathy (FinnDiane) Study with type 1 diabetes to determine whether sex-specific genetic risk factors for ESRD exist. A common variant, rs4972593 on chromosome 2q31. 1, was associated with ESRD in women (P

Genetics of New-Onset Diabetes after Transplantation

New-onset diabetes after transplantation is a common complication that reduces recipient survival. Research in renal transplant recipients has suggested that pancreatic ß-cell dysfunction, as opposed to insulin resistance, may be the key pathologic process. In this study, clinical and genetic factors associated with new-onset diabetes after transplantation were identified in a white population. A joint analysis approach, with an initial genome-wide association study in a subset of cases followed by de novo genotyping in the complete case cohort, was implemented to identify single-nucleotide polymorphisms (SNPs) associated with the development of new-onset diabetes after transplantation. Clinical variables associated with the development of diabetes after renal transplantation included older recipient age, female sex, and percentage weight gain within 12 months of transplantation. The genome-wide association study identified 26 SNPs associated with new-onset diabetes after transplantation; this association was validated for eight SNPs (rs10484821, rs7533125, rs2861484, rs11580170, rs2020902, rs1836882, rs198372, and rs4394754) by de novo genotyping. These associations remained significant after multivariate adjustment for clinical variables. Seven of these SNPs are associated with genes implicated in ß-cell apoptosis. These results corroborate recent clinical evidence implicating ß-cell dysfunction in the pathophysiology of new-onset diabetes after transplantation and support the pursuit of therapeutic strategies to protect ß cells in the post-transplant period.

Variants in the 1q21 risk region are associated with a visual endophenotype of autism and schizophrenia

Deficits in sensitivity to visual stimuli of low spatial frequency and high temporal frequency (so-called frequency-doubled gratings) have been demonstrated both in schizophrenia and in autism spectrum disorder (ASD). Such basic perceptual functions are ideal candidates for molecular genetic study, because the underlying neural mechanisms are well characterized; but they have sometimes been overlooked in favor of cognitive and neurophysiological endophenotypes, for which neural substrates are often unknown. Here, we report a genome-wide association study of a basic visual endophenotype associated with psychological disorder. Sensitivity to frequency-doubled gratings was measured in 1060 healthy young adults, and analyzed for association with genotype using linear regression at 642758 single nucleotide polymorphism (SNP) markers. A significant association (P=7.9×10) was found with the SNP marker rs1797052, situated in the 5′-untranslated region of PDZK1; each additional copy of the minor allele was associated with an increase in sensitivity equivalent to more than half a standard deviation. A permutation procedure, which accounts for multiple testing, showed that the association was significant at the α=0.005 level. The region on chromosome 1q21.1 surrounding PDZK1 is an established susceptibility locus both for schizophrenia and for ASD, mirroring the common association of the visual endophenotype with the two disorders. PDZK1 interacts with N-methyl-d-aspartate receptors and neuroligins, which have been implicated in the etiologies of schizophrenia and ASD. These findings suggest that perceptual abnormalities observed in two different disorders may be linked by common genetic elements. © 2013 John Wiley & Sons Ltd and International Behavioural and Neural Genetics Society.

Individual differences provide psychophysical evidence for separate on- and off-pathways deriving from short-wave cones

Distinct neural populations carry signals from short-wave (S) cones. We used individual differences to test whether two types of pathways, those that receive excitatory input (S+) and those that receive inhibitory input (S-), contribute independently to psychophysical performance. We also conducted a genome-wide association study (GWAS) to look for genetic correlates of the individual differences. Our psychophysical test was based on the Cambridge Color Test, but detection thresholds were measured separately for S-cone spatial increments and decrements. Our participants were 1060 healthy adults aged 16-40. Test-retest reliabilities for thresholds were good (ρ=0.64 for S-cone increments, 0.67 for decrements and 0.73 for the average of the two). "Regression scores," isolating variability unique to incremental or decremental sensitivity, were also reliable (ρ=0.53 for increments and ρ=0.51 for decrements). The correlation between incremental and decremental thresholds was ρ=0.65. No genetic markers reached genome-wide significance (p-7). We identified 18 "suggestive" loci (p-5). The significant test-retest reliabilities show stable individual differences in S-cone sensitivity in a normal adult population. Though a portion of the variance in sensitivity is shared between incremental and decremental sensitivity, over 26% of the variance is stable across individuals, but unique to increments or decrements, suggesting distinct neural substrates. Some of the variability in sensitivity is likely to be genetic. We note that four of the suggestive associations found in the GWAS are with genes that are involved in glucose metabolism or have been associated with diabetes.

Novel genetic susceptibility loci for diabetic end-stage renal disease identified through robust naive Bayes classification

Aims/hypothesis: Diabetic nephropathy is a major diabetic complication, and diabetes is the leading cause of end-stage renal disease (ESRD). Family studies suggest a hereditary component for diabetic nephropathy. However, only a few genes have been associated with diabetic nephropathy or ESRD in diabetic patients. Our aim was to detect novel genetic variants associated with diabetic nephropathy and ESRD. Methods: We exploited a novel algorithm, ‘Bag of Naive Bayes’, whose marker selection strategy is complementary to that of conventional genome-wide association models based on univariate association tests. The analysis was performed on a genome-wide association study of 3,464 patients with type 1 diabetes from the Finnish Diabetic Nephropathy (FinnDiane) Study and subsequently replicated with 4,263 type 1 diabetes patients from the Steno Diabetes Centre, the All Ireland-Warren 3-Genetics of Kidneys in Diabetes UK collection (UK–Republic of Ireland) and the Genetics of Kidneys in Diabetes US Study (GoKinD US). Results: Five genetic loci (WNT4/ZBTB40-rs12137135, RGMA/MCTP2-rs17709344, MAPRE1P2-rs1670754, SEMA6D/SLC24A5-rs12917114 and SIK1-rs2838302) were associated with ESRD in the FinnDiane study. An association between ESRD and rs17709344, tagging the previously identified rs12437854 and located between the RGMA and MCTP2 genes, was replicated in independent case–control cohorts. rs12917114 near SEMA6D was associated with ESRD in the replication cohorts under the genotypic model (p < 0.05), and rs12137135 upstream of WNT4 was associated with ESRD in Steno. Conclusions/interpretation: This study supports the previously identified findings on the RGMA/MCTP2 region and suggests novel susceptibility loci for ESRD. This highlights the importance of applying complementary statistical methods to detect novel genetic variants in diabetic nephropathy and, in general, in complex diseases.

Biological insights from 108 schizophrenia-associated genetic loci

Schizophrenia is a highly heritable disorder. Genetic risk is conferred by a large number of alleles, including common alleles of small effect that might be detected by genome-wide association studies. Here we report a multi-stage schizophrenia genome-wide association study of up to 36,989 cases and 113,075 controls. We identify 128 independent associations spanning 108 conservatively defined loci that meet genome-wide significance, 83 of which have not been previously reported. Associations were enriched among genes expressed in brain, providing biological plausibility for the findings. Many findings have the potential to provide entirely new insights into aetiology, but associations at DRD2 and several genes involved in glutamatergic neurotransmission highlight molecules of known and potential therapeutic relevance to schizophrenia, and are consistent with leading pathophysiological hypotheses. Independent of genes expressed in brain, associations were enriched among genes expressed in tissues that have important roles in immunity, providing support for the speculated link between the immune system and schizophrenia.

Genetic studies of body mass index yield new insights for obesity biology

Obesity is heritable and predisposes to many diseases. To understand the genetic basis of obesity better, here we conduct a genome-wide association study and Metabochip meta-analysis of body mass index (BMI), a measure commonly used to define obesity and assess adiposity, in up to 339,224 individuals. This analysis identifies 97 BMI-associated loci (P < 5 × 10(-8)), 56 of which are novel. Five loci demonstrate clear evidence of several independent association signals, and many loci have significant effects on other metabolic phenotypes. The 97 loci account for ∼2.7% of BMI variation, and genome-wide estimates suggest that common variation accounts for >20% of BMI variation. Pathway analyses provide strong support for a role of the central nervous system in obesity susceptibility and implicate new genes and pathways, including those related to synaptic function, glutamate signalling, insulin secretion/action, energy metabolism, lipid biology and adipogenesis.

Polymorphisms near TBX5 and GDF7 are associated with increased risk for Barrett's esophagus

BACKGROUND & AIMS: Barrett's esophagus (BE) increases the risk of esophageal adenocarcinoma (EAC). We found the risk to be BE has been associated with single nucleotide polymorphisms (SNPs) on chromosome 6p21 (within the HLA region) and on 16q23, where the closest protein-coding gene is FOXF1. Subsequently, the Barrett's and Esophageal Adenocarcinoma Consortium (BEACON) identified risk loci for BE and esophageal adenocarcinoma near CRTC1 and BARX1, and within 100 kb of FOXP1. We aimed to identify further SNPs that increased BE risk and to validate previously reported associations.

METHODS: We performed a genome-wide association study (GWAS) to identify variants associated with BE and further analyzed promising variants identified by BEACON by genotyping 10,158 patients with BE and 21,062 controls.

RESULTS: We identified 2 SNPs not previously associated with BE: rs3072 (2p24.1; odds ratio [OR] = 1.14; 95% CI: 1.09-1.18; P = 1.8 × 10(-11)) and rs2701108 (12q24.21; OR = 0.90; 95% CI: 0.86-0.93; P = 7.5 × 10(-9)). The closest protein-coding genes were respectively GDF7 (rs3072), which encodes a ligand in the bone morphogenetic protein pathway, and TBX5 (rs2701108), which encodes a transcription factor that regulates esophageal and cardiac development. Our data also supported in BE cases 3 risk SNPs identified by BEACON (rs2687201, rs11789015, and rs10423674). Meta-analysis of all data identified another SNP associated with BE and esophageal adenocarcinoma: rs3784262, within ALDH1A2 (OR = 0.90; 95% CI: 0.87-0.93; P = 3.72 × 10(-9)).

CONCLUSIONS: We identified 2 loci associated with risk of BE and provided data to support a further locus. The genes we found to be associated with risk for BE encode transcription factors involved in thoracic, diaphragmatic, and esophageal development or proteins involved in the inflammatory response.

Independent genomewide screens identify the tumor suppressor VTRNA2-1 as a human epiallele responsive to periconceptional environment

Background: Interindividual epigenetic variation that occurs systemically must be established prior to gastrulation in the very early embryo and, because it is systemic, can be assessed in easily biopsiable tissues. We employ two independent genome-wide approaches to search for such variants.

Results: First, we screen for metastable epialleles by performing genomewide bisulfite sequencing in peripheral blood lymphocyte (PBL) and hair follicle DNA from two Caucasian adults. Second, we conduct a genomewide screen for genomic regions at which PBL DNA methylation is affected by season of conception in rural Gambia. Remarkably, both approaches identify the genomically imprinted VTRNA2-1 as a top environmentally responsive epiallele. We demonstrate systemic and stochastic interindividual variation in DNA methylation at the VTRNA2-1 differentially methylated region in healthy Caucasian and Asian adults and show, in rural Gambians, that periconceptional environment affects offspring VTRNA2-1 epigenotype, which is stable over at least 10 years. This unbiased screen also identifies over 100 additional candidate metastable epialleles, and shows that these are associated with cis genomic features including transposable elements.

Conclusions: The non-coding VTRNA2-1 transcript (also called nc886) is a putative tumor suppressor and modulator of innate immunity. Thus, these data indicating environmentally induced loss of imprinting at VTRNA2-1 constitute a plausible causal pathway linking early embryonic environment, epigenetic alteration, and human disease. More broadly, the list of candidate metastable epialleles provides a resource for future studies of epigenetic variation and human disease.