An Excess of Risk-Increasing Low-Frequency Variants Can Be a Signal of Polygenic Inheritance in Complex Diseases

In most complex diseases, much of the heritability remains unaccounted for by common variants. It has been postulated that lower-frequency variants contribute to the remaining heritability. Here, we describe a method to test for polygenic inheritance from lower-frequency variants by using GWAS summary association statistics. We explored scenarios with many causal low-frequency variants and showed that there is more power to detect risk variants than to detect protective variants, resulting in an increase in the ratio of detected risk to protective variants (R/P ratio). Such an excess can also occur if risk variants are present and kept at lower frequencies because of negative selection. The R/P ratio can be falsely elevated because of reasons unrelated to polygenic inheritance, such as uneven sample sizes or asymmetric population stratification, so precautions to correct for these confounders are essential. We tested our method on published GWAS results and observed a strong signal in some diseases (schizophrenia and type 2 diabetes) but not others. We also explored the shared genetic component in overlapping phenotypes related to inflammatory bowel disease (Crohn disease [CD] and ulcerative colitis [UC]) and diabetic nephropathy (macroalbuminuria and end-stage renal disease [ESRD]). Although the signal was still present when both CD and UC were jointly analyzed, the signal was lost when macroalbuminuria and ESRD were jointly analyzed, suggesting that these phenotypes should best be studied separately. Thus, our method may also help guide the design of future genetic studies of various traits and diseases.

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.

FLT1 Genetic Variation Predisposes to Neovascular AMD in Ethnically Diverse Populations and Alters Systemic FLT1 Expression

Purpose: Current understanding of the genetic risk factors for age-related macular degeneration (AMD) is not sufficiently predictive of the clinical course. The VEGF pathway is a key therapeutic target for treatment of neovascular AMD; however, risk attributable to genetic variation within pathway genes is unclear. We sought to identify single nucleotide polymorphisms (SNPs) associated with AMD within the VEGF pathway.
Methods: Using a tagSNP, direct sequencing and meta-analysis approach within four ethnically diverse cohorts, we identified genetic risk present in FLT1, though not within other VEGF pathway genes KDR, VEGFA, or VASH1. We used ChIP and ELISA in functional analysis.
Results: The FLT1 SNPs rs9943922, rs9508034, rs2281827, rs7324510, and rs9513115 were significantly associated with increased risk of neovascular AMD. Each association was more significant after meta-analysis than in any one of the four cohorts. All associations were novel, within noncoding regions of FLT1 that do not tag for coding variants in linkage disequilibrium. Analysis of soluble FLT1 demonstrated higher expression in unaffected individuals homozygous for the FLT1 risk alleles rs9943922 (P = 0.0086) and rs7324510 (P = 0.0057). In silico analysis suggests that these variants change predicted splice sites and RNA secondary structure, and have been identified in other neovascular pathologies. These data were supported further by murine chromatin immunoprecipitation demonstrating that FLT1 is a target of Nr2e3, a nuclear receptor gene implicated in regulating an AMD pathway.
Conclusions: Although exact variant functions are not known, these data demonstrate relevancy across ethnically diverse genetic backgrounds within our study and, therefore, hold potential for global efficacy.

Rare and common variants in extracellular matrix gene Fibrillin 2 (FBN2) are associated with macular degeneration

Neurodegenerative diseases affecting the macula constitute a major cause of incurable vision loss and exhibit considerable clinical and genetic heterogeneity, from early-onset monogenic disease to multifactorial late-onset age-related macular degeneration (AMD). As part of our continued efforts to define genetic causes of macular degeneration, we performed whole exome sequencing in four individuals of a two-generation family with autosomal dominant maculopathy and identified a rare variant p.Glu1144Lys in Fibrillin 2 (FBN2), a glycoprotein of the elastin-rich extracellular matrix (ECM). Sanger sequencing validated the segregation of this variant in the complete pedigree, including two additional affected and one unaffected individual. Sequencing of 192 maculopathy patients revealed additional rare variants, predicted to disrupt FBN2 function. We then undertook additional studies to explore the relationship of FBN2 to macular disease. We show that FBN2 localizes to Bruch's membrane and its expression appears to be reduced in aging and AMD eyes, prompting us to examine its relationship with AMD. We detect suggestive association of a common FBN2 non-synonymous variant, rs154001 (p.Val965Ile) with AMD in 10,337 cases and 11,174 controls (OR=1.10; p-value=3.79×10(-5)). Thus, it appears that rare and common variants in a single gene - FBN2 - can contribute to Mendelian and complex forms of macular degeneration. Our studies provide genetic evidence for a key role of elastin microfibers and Bruch's membrane in maintaining blood-retina homeostasis and establish the importance of studying orphan diseases for understanding more common clinical phenotypes.

Whole-mitochondrial genome sequencing in primary open-angle glaucoma using massively parallel sequencing identifies novel and known pathogenic variants

Purpose:The aim of this study was to determine whether mutations in mitochondrial DNA play a role in high-pressure primary open-angle glaucoma (OMIM 137760) by analyzing new data from massively parallel sequencing of mitochondrial DNA.
Methods:Glaucoma patients with high-tension primary open-angle glaucoma and ethnically matched and age-matched control subjects without glaucoma were recruited. The entire human mitochondrial genome was amplified in two overlapping fragments by long-range polymerase chain reaction and used as a template for massively parallel sequencing on an Ion Torrent Personal Genome Machine. All variants were confirmed by conventional Sanger sequencing.
Results:Whole-mitochondrial genome sequencing was performed in 32 patients with primary open-angle glaucoma from India (n = 16) and Ireland (n = 16). In 16 of the 32 patients with primary open-angle glaucoma (50% of cases), there were 22 mitochondrial DNA mutations consisting of 7 novel mutations and 8 previously reported disease-associated sequence variants. Eight of 22 (36.4%) of the mitochondrial DNA mutations were in complex I mitochondrial genes.
Conclusion:Massively parallel sequencing using the Ion Torrent Personal Genome Machine with confirmation by Sanger sequencing detected a pathogenic mitochondrial DNA mutation in 50% of the primary open-angle glaucoma cohort. Our findings support the emerging concept that mitochondrial dysfunction results in the development of glaucoma and, more specifically, that complex I defects play a significant role in primary open-angle glaucoma pathogenesis.

FAD binding overcomes defects in activity and stability displayed by cancer-associated variants of human NQO1

NAD(P)H quinone oxidoreductase 1 is involved in antioxidant defence and protection from cancer, stabilizing the apoptosis regulator p53 towards degradation. Here, we studied the enzymological, biochemical and biophysical properties of two cancer-associated variants (p.R139W and p.P187S). Both variants (especially p.187S) have lower thermal stability and greater susceptibility to proteolysis compared to the wild-type. p.P187S also has reduced activity due to a lower binding affinity for the FAD cofactor as assessed by activity measurements and direct titrations. Native gel electrophoresis and dynamic light scattering also suggest that p.P187S has a higher tendency to populate unfolded states under native conditions. Detailed thermal stability studies showed that all variants irreversibly denature causing dimer dissociation, while addition of FAD restores the stability of the polymorphic forms to wild-type levels. The kinetic destabilization induced by polymorphisms as well as the kinetic protection exerted by FAD was confirmed by measuring denaturation kinetics at temperatures close to physiological. Our data suggest that the main molecular mechanisms associated with these cancer-related variants are their low binding affinity for FAD and/or kinetic instability. Thus, pharmacological chaperones may be useful in the treatment of patients bearing these polymorphisms.

The metastability of human UDP-galactose 4'-epimerase (GALE) is increased by variants associated with type III galactosemia but decreased by substrate and cofactor binding

Type III galactosemia is an inherited disease caused by mutations which affect the activity of UDP-galactose 4'-epimerase (GALE). We evaluated the impact of four disease-associated variants (p.N34S, p.G90E, p.V94M and p.K161N) on the conformational stability and dynamics of GALE. Thermal denaturation studies showed that wild-type GALE denatures at temperatures close to physiological, and disease-associated mutations often reduce GALE's thermal stability. This denaturation is under kinetic control and results partly from dimer dissociation. The natural ligands, NAD(+) and UDP-glucose, stabilize GALE. Proteolysis studies showed that the natural ligands and disease-associated variations affect local dynamics in the N-terminal region of GALE. Proteolysis kinetics followed a two-step irreversible model in which the intact protein is cleaved at Ala38 forming a long-lived intermediate in the first step. NAD(+) reduces the rate of the first step, increasing the amount of undigested protein whereas UDP-glucose reduces the rate of the second step, increasing accumulation of the intermediate. Disease-associated variants affect these rates and the amounts of protein in each state. Our results also suggest communication between domains in GALE. We hypothesize that, in vivo, concentrations of natural ligands modulate GALE stability and that it should be possible to discover compounds which mimic the stabilising effects of the natural ligands overcoming mutation-induced destabilization.

Partitioning Heritability of Regulatory and Cell-Type-Specific Variants across 11 Common Diseases

Regulatory and coding variants are known to be enriched with associations identified by genome-wide association studies (GWASs) of complex disease, but their contributions to trait heritability are currently unknown. We applied variance-component methods to imputed genotype data for 11 common diseases to partition the heritability explained by genotyped SNPs () across functional categories (while accounting for shared variance due to linkage disequilibrium). Extensive simulations showed that in contrast to current estimates from GWAS summary statistics, the variance-component approach partitions heritability accurately under a wide range of complex-disease architectures. Across the 11 diseases DNaseI hypersensitivity sites (DHSs) from 217 cell types spanned 16% of imputed SNPs (and 24% of genotyped SNPs) but explained an average of 79% (SE = 8%) of  from imputed SNPs (5.1× enrichment; p = 3.7 × 10⁻¹⁷) and 38% (SE = 4%) of  from genotyped SNPs (1.6× enrichment, p = 1.0 × 10⁻⁴). Further enrichment was observed at enhancer DHSs and cell-type-specific DHSs. In contrast, coding variants, which span 1% of the genome, explained <10% of  despite having the highest enrichment. We replicated these findings but found no significant contribution from rare coding variants in independent schizophrenia cohorts genotyped on GWAS and exome chips. Our results highlight the value of analyzing components of heritability to unravel the functional architecture of common disease.

Shared common variants in prostate cancer and blood lipids

BACKGROUND: Epidemiological and clinical studies suggest comorbidity between prostate cancer (PCA) and cardiovascular disease (CVD) risk factors. However, the relationship between these two phenotypes is still not well understood. Here we sought to identify shared genetic loci between PCA and CVD risk factors.

METHODS: We applied a genetic epidemiology method based on conjunction false discovery rate (FDR) that combines summary statistics from different genome-wide association studies (GWAS), and allows identification of genetic overlap between two phenotypes. We evaluated summary statistics from large, multi-centre GWA studies of PCA (n=50 000) and CVD risk factors (n=200 000) [triglycerides (TG), low-density lipoprotein (LDL) cholesterol and high-density lipoprotein (HDL) cholesterol, systolic blood pressure, body mass index, waist-hip ratio and type 2 diabetes (T2D)]. Enrichment of single nucleotide polymorphisms (SNPs) associated with PCA and CVD risk factors was assessed with conditional quantile-quantile plots and the Anderson-Darling test. Moreover, we pinpointed shared loci using conjunction FDR.

RESULTS: We found the strongest enrichment of P-values in PCA was conditional on LDL and conditional on TG. In contrast, we found only weak enrichment conditional on HDL or conditional on the other traits investigated. Conjunction FDR identified altogether 17 loci; 10 loci were associated with PCA and LDL, 3 loci were associated with PCA and TG and additionally 4 loci were associated with PCA, LDL and TG jointly (conjunction FDR <0.01). For T2D, we detected one locus adjacent to HNF1B.

CONCLUSIONS: We found polygenic overlap between PCA predisposition and blood lipids, in particular LDL and TG, and identified 17 pleiotropic gene loci between PCA and LDL, and PCA and TG, respectively. These findings provide novel pathobiological insights and may have implications for trials using targeting lipid-lowering agents in a prevention or cancer setting.

Integrin genetic variants and stage-specific tumor recurrence in patients with stage II and III colon cancer

Integrins (ITGs) are key elements in cancer biology, regulating tumor growth, angiogenesis and lymphangiogenesis through interactions of the tumor cells with the microenvironment. Moving from the hypothesis that ITGs could have different effects in stage II and III colon cancer, we tested whether a comprehensive panel of germline single-nucleotide polymorphisms (SNPs) in ITG genes could predict stage-specific time to tumor recurrence (TTR). A total of 234 patients treated with 5-fluorouracil-based chemotherapy at the University of Southern California were included in this study. Whole-blood samples were analyzed for germline SNPs in ITG genes using PCR-restriction fragment length polymorphism or direct DNA sequencing. In the multivariable analysis, stage II colon cancer patients with at least one G allele for ITGB3 rs4642 had higher risk of recurrence (hazard ratio (HR)=4.027, 95% confidence interval (95% CI) 1.556-10.421, P=0.004). This association was also significant in the combined stage II-III cohort (HR=1.975, 95% CI 1.194-3.269, P=0.008). The predominant role of ITGB3 rs4642 in stage II diseases was confirmed using recursive partitioning, showing that ITGB3 rs4642 was the most important factor in stage II diseases. In contrast, in stage III diseases the combined analysis of ITGB1 rs2298141 and ITGA4 rs7562325 allowed to identify three distinct prognostic subgroups (P=0.009). The interaction between stage and the combined ITGB1 rs2298141 and ITGA4 rs7562325 on TTR was significant (P=0.025). This study identifies germline polymorphisms in ITG genes as independent stage-specific prognostic markers for stage II and III colon cancer. These data may help to select subgroups of patients who may benefit from ITG-targeted treatments.

Genetic variants of kinase suppressors of Ras (KSR1) to predict survival in patients with ERα-positive advanced breast cancer

In patients with breast cancer (BC), deregulation of estrogen receptor (ERα) activity may account for most resistance to endocrine therapies. Our previous study used a whole-human kinome siRNA screen to identify functional actors in ERα modulation and showed the implication of proteins kinase suppressors of ras (KSR1). From those findings we evaluated the clinical impact of KSR1 variants in patients with ERα+ BC treated with TAM. DNA was obtained from 222 patients with advanced ERα+ BC treated with TAM who had undergone surgery from 1981 to 2003. We selected three potentially functional relevant KSR1 polymorphisms; two within the 3'UTR (rs224190, rs1075952) and one in the coding exon 7 (rs2293180). The primary end points were overall survival (OS) and disease-free survival (DFS). After a 6.4-year median follow-up, patients carrying the rs2241906 TT genotype showed shorter DFS (2.1 vs 7.1 years, P=0.005) and OS (2.6 vs 8.4 years P=0.002) than those with the TC or TT genotypes. Those associations remained significant in the multivariable analysis adjusting age, lymph node status, LMTK3 and IGFR variants and HER2 status. The polymorphisms rs2241906 and rs1075952 were in linkage disequilibrium. No association was shown between rs2293180 and survival. Among the actors of ERα signaling, KSR1 rs2241906 variants may predict survival in patients with advanced ERα+ BC treated with adjuvant TAM.