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.

Association of common gene variants in the WNT/β-catenin pathway with colon cancer recurrence

Wnt/β-catenin signaling has a central role in the development and progression of most colon cancers (CCs). Germline variants in Wnt/β-catenin pathway genes may result in altered gene function and/or activity, thereby causing inter-individual differences in relation to tumor recurrence capacity and chemoresistance. We investigated germline polymorphisms in a comprehensive panel of Wnt/β-catenin pathway genes to predict time to tumor recurrence (TTR) in patients with stage III and high-risk stage II CC. A total of 234 patients treated with 5-fluorouracil-based chemotherapy were included in this study. Whole-blood samples were analyzed for putative functional germline polymorphisms in SFRP3, SFRP4, DKK2, DKK3, Axin2, APC, TCF7L2, WNT5B, CXXC4, NOTCH2 and GLI1 genes by PCR-based restriction fragment-length polymorphism or direct DNA sequencing. Polymorphisms with statistical significance were validated in an independent study cohort. The minor allele of WNT5B rs2010851 T>G was significantly associated with a shorter TTR (10.7 vs 4.9 years; hazard ratio: 2.48; 95% CI, 0.96-6.38; P=0.04) in high-risk stage II CC patients. This result remained significant in multivariate Cox's regression analysis. This study shows that the WNT5B germline variant rs2010851 was significantly identified as a stage-dependent prognostic marker for CC patients after 5-fluorouracil-based adjuvant therapy.

Common cancer stem cell gene variants predict colon cancer recurrence

PURPOSE: Recent evidence suggests that cancer stem cells (CSC) are responsible for key elements of colon cancer progression and recurrence. Germline variants in CSC genes may result in altered gene function and/or activity, thereby causing interindividual differences in a patient's tumor recurrence capacity and chemoresistance. We investigated germline polymorphisms in a comprehensive panel of CSC genes to predict time to tumor recurrence (TTR) in patients with stage III and high-risk stage II colon cancer.

EXPERIMENTAL DESIGN: 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 polymorphisms in genes that have been previously associated with colon CSC (CD44, Prominin-1, DPP4, EpCAM, ALCAM, Msi-1, ITGB1, CD24, LGR5, and ALDH1A1) by PCR-RFLP or direct DNA-sequencing.

RESULTS: The minor alleles of CD44 rs8193 C>T, ALCAM rs1157 G>A, and LGR5 rs17109924 T>C were significantly associated with increased TTR (9.4 vs. 5.4 years; HR, 0.51; 95% CI: 0.35-0.93; P = 0.022; 11.3 vs. 5.7 years; HR, 0.56; 95% CI: 0.33-0.94; P = 0.024, and 10.7 vs. 5.7 years; HR, 0.33; 95% CI: 0.12-0.90; P = 0.023, respectively) and remained significant in the multivariate analysis stratified by ethnicity. In recursive partitioning, a specific gene variant profile including LGR5 rs17109924, CD44 rs8193, and ALDH1A1 rs1342024 represented a high-risk subgroup with a median TTR of 1.7 years (HR, 6.71, 95% CI: 2.71-16.63, P < 0.001).

CONCLUSION: This is the first study identifying common germline variants in colon CSC genes as independent prognostic markers for stage III and high-risk stage II colon cancer patients.

Prediction of breast cancer risk based on profiling with common genetic variants

BACKGROUND: Data for multiple common susceptibility alleles for breast cancer may be combined to identify women at different levels of breast cancer risk. Such stratification could guide preventive and screening strategies. However, empirical evidence for genetic risk stratification is lacking. METHODS: We investigated the value of using 77 breast cancer-associated single nucleotide polymorphisms (SNPs) for risk stratification, in a study of 33 673 breast cancer cases and 33 381 control women of European origin. We tested all possible pair-wise multiplicative interactions and constructed a 77-SNP polygenic risk score (PRS) for breast cancer overall and by estrogen receptor (ER) status. Absolute risks of breast cancer by PRS were derived from relative risk estimates and UK incidence and mortality rates. RESULTS: There was no strong evidence for departure from a multiplicative model for any SNP pair. Women in the highest 1% of the PRS had a three-fold increased risk of developing breast cancer compared with women in the middle quintile (odds ratio [OR] = 3.36, 95% confidence interval [CI] = 2.95 to 3.83). The ORs for ER-positive and ER-negative disease were 3.73 (95% CI = 3.24 to 4.30) and 2.80 (95% CI = 2.26 to 3.46), respectively. Lifetime risk of breast cancer for women in the lowest and highest quintiles of the PRS were 5.2% and 16.6% for a woman without family history, and 8.6% and 24.4% for a woman with a first-degree family history of breast cancer. CONCLUSIONS: The PRS stratifies breast cancer risk in women both with and without a family history of breast cancer. The observed level of risk discrimination could inform targeted screening and prevention strategies. Further discrimination may be achievable through combining the PRS with lifestyle/environmental factors, although these were not considered in this report.

A common major histocompatibility complex class II allele HLA-DQB1* 0301 is present in clinical variants of pemphigoid.

Bullous pemphigoid (BP) is an autoimmune subepidermal blistering disease seen primarily in elderly persons. It is characterized clinically by the development of tense bullae and by the presence of an antibasement membrane antibody. In BP, the antigens involved in the autoimmunity are epidermal basement membrane peptides BPAg1 and BPAg2. We have compared high resolution typing of major histocompatibility complex class II loci (HLA-DRB1, DQB1) in 21 patients with BP, 17 with ocular cicatricial pemphigoid (OCP), and 22 with oral pemphigoid (OP) to a panel of 218 haplotypes of normal individuals. We found that the three diseases (BP, OCP, and OP) have significant association with DQB1*0301 (P = 0.005, P < 0.0001, and P = 0.001, respectively). The frequencies of alleles DQB1*0302, 0303, and 06, which share a specific amino acid sequence from position 71 to 77 (Thr-Arg-Ala-Glu-Leu-Val-Thr) were also increased (P = 0.01). We suggest that an identical major histocompatibility complex class II allele (DQB1*0301) is a common marker for enhanced susceptibility and that the same amino acid residues in positions 71-77 (DQB1*0301, -0302, -0305, -0602, -0603 alleles) are found in patients with BP, OCP and OP. Our findings propose that the autoimmune response in the three different clinical variants of pemphigoid, involves the recognition by T cells of a class II region of DQB1, bound to a peptide from the basement membrane of conjunctiva, oral mucosa, and skin.

Common E protein determinants for attenuation of glycosaminoglycan-binding variants of Japanese encephalitis and West Nile viruses

Natural isolates and laboratory strains of West Nile virus (WNV) and Japanese encephalitis virus (JEV) were attenuated for neuroinvasiveness in mouse models for flavivirus encephalitis by serial passage in human adenocarcinoma (SW13) cells. The passage variants displayed a small-plaque phenotype, augmented affinity for heparin-Sepharose, and a marked increase in specific infectivity for SW13 cells relative to the respective parental viruses, while the specific infectivity for Vero cells was not altered. Therefore, host cell adaptation of passage variants was most likely a consequence of altered receptor usage for virus attachment-entry with the involvement of cell surface glycosaminoglycans (GAG) in this process. In vivo blood clearance kinetics of the passage variants was markedly faster and viremia was reduced relative to the parental viruses, suggesting that affinity for GAG (ubiquitously present on cell surfaces and extracellular matrices) is a key determinant for the neuroinvasiveness of encephalitic flaviviruses. A difference in pathogenesis between WNV and JEV, which was reflected in more efficient growth in the spleen and liver of the WNV parent and passage variants, accounted for a less pronounced loss of neuroinvasiveness of GAG binding variants of WNV than JEV. Single gain-of-net-positive-charge amino acid changes at E protein residue 49, 138, 306, or 389/390, putatively positioned in two clusters on the virion surface, define molecular determinants for GAG binding and concomitant virulence attenuation that are shared by the JEV serotype flaviviruses.

No evidence for association of ataxia-telangiectasia mutated gene T2119C and C3161G amino acid substitution variants with risk of breast cancer

Background There is evidence that certain mutations in the double-strand break repair pathway ataxia-telangiectasia mutated gene act in a dominant-negative manner to increase the risk of breast cancer. There are also some reports to suggest that the amino acid substitution variants T2119C Ser707Pro and C3161G Pro1054Arg may be associated with breast cancer risk. We investigate the breast cancer risk associated with these two nonconservative amino acid substitution variants using a large Australian population-based case–control study. Methods The polymorphisms were genotyped in more than 1300 cases and 600 controls using 5' exonuclease assays. Case–control analyses and genotype distributions were compared by logistic regression. Results The 2119C variant was rare, occurring at frequencies of 1.4 and 1.3% in cases and controls, respectively (P = 0.8). There was no difference in genotype distribution between cases and controls (P = 0.8), and the TC genotype was not associated with increased risk of breast cancer (adjusted odds ratio = 1.08, 95% confidence interval = 0.59–1.97, P = 0.8). Similarly, the 3161G variant was no more common in cases than in controls (2.9% versus 2.2%, P = 0.2), there was no difference in genotype distribution between cases and controls (P = 0.1), and the CG genotype was not associated with an increased risk of breast cancer (adjusted odds ratio = 1.30, 95% confidence interval = 0.85–1.98, P = 0.2). This lack of evidence for an association persisted within groups defined by the family history of breast cancer or by age. Conclusion The 2119C and 3161G amino acid substitution variants are not associated with moderate or high risks of breast cancer in Australian women.

Multiple analysis of three common genetic alterations associated with thrombophilia

We have previously reported the use of a novel mini-sequencing protocol for detection of the factor V Leiden variant, the first nucleotide change (FNC) technology. This technology is based on a single nucleotide extension of a primer, which is hybridized immediately adjacent to the site of mutation. The extended nucleotide that carries a reporter molecule (fluorescein) has the power to discriminate the genotype at the site of mutation. More recently, the prothrombin 20210 and thermolabile methylene tetrahydrofolate reductase (MTHFR) 677 variants have been identified as possible risk factors associated with thrombophilia. This study describes the use of the FNC technology in a combined assay to detect factor V, prothrombin and MTHFR variants in a population of Australian blood donors, and describes the objective numerical methodology used to determine genotype cut-off values for each genetic variation. Using FNC to test 500 normal blood donors, the incidence of Factor V Leiden was 3.6% (all heterozygous), that of prothrombin 20210 was 2.8% (all heterozygous) and that of MTHFR was 10% (homozygous). The combined FNC technology offers a simple, rapid, automatable DNA-based test for the detection of these three important mutations that are associated with familial thrombophilia. (C) 2000 Lippincott Williams and Wilkins.

Investigation of lymphotoxin α genetic variants in migraine

Migraine is a common neurological disease with a genetic basis affecting approximately 12% of the population. Pain during a migraine attack is associated with activation of the trigeminal nerve system, which carries pain signals from the meninges and the blood vessels infusing the meninges to the trigeminal nucleus in the brain stem. The release of inflammatory mediators following cortical spreading depression (CSD) may further promote and sustain the activation and sensitization of meningeal nociceptors, inducing the persistent throbbing headache characterised in migraine. Lymphotoxin α (LTA) is a cytokine secreted by lymphocytes and is a member of the tumour necrosis factor (TNF) family. Genetic variation with the TNF and LTA genes may contribute to threshold brain excitability, propagation of neuronal hyperexcitability and thus initiation and maintenance of a migraine attack. Three LTA variants rs2009658, rs2844482 and rs2229094 were identified in a recent pGWAS study conducted in the Norfolk Island population as being potentially implicated in migraine with nominally significant p values of p = 0.0093, p = 0.0088 and p = 0.033 respectively. To determine whether these SNPs played a role in migraine in a general outbred population these SNPs were gentoyped in a large case control Australian Caucasian population and tested for association with migraine. All three SNPs showed no association in our cohort (p > 0.05). Validation of GWAS data in independent case-controls cohorts is essential to establish risk validity within specific population groups. The importance of cytokines in modulating neural inflammation and pain threshold in addition to other studies showing associations between TNF-α and SNPs in the LTA gene with migraine, suggests that LTA could be an important factor contributing to migraine. Although the present study did not support a role for the tested LTA variants in migraine, investigation of other variants within the LTA gene is still warranted.

A genome-wide analysis of 'Bounty' descendants implicates several novel variants in migraine susceptibility

Migraine is a common neurological disease with a complex genetic aetiology. The disease affects ~12% of the Caucasian population and females are three times more likely than males to be diagnosed. In an effort to identify loci involved in migraine susceptibility, we performed a pedigree-based genome-wide association study of the isolated population of Norfolk Island, which has a high prevalence of migraine. This unique population originates from a small number of British and Polynesian founders who are descendents of the Bounty mutiny and forms a very large multigenerational pedigree (Bellis et al.; Human Genetics, 124(5):543-5542, 2008). These population genetic features may facilitate disease gene mapping strategies (Peltonen et al.; Nat Rev Genet, 1(3):182-90, 2000. In this study, we identified a high heritability of migraine in the Norfolk Island population (h (2) = 0.53, P = 0.016). We performed a pedigree-based GWAS and utilised a statistical and pathological prioritisation approach to implicate a number of variants in migraine. An SNP located in the zinc finger protein 555 (ZNF555) gene (rs4807347) showed evidence of statistical association in our Norfolk Island pedigree (P = 9.6 × 10(-6)) as well as replication in a large independent and unrelated cohort with >500 migraineurs. In addition, we utilised a biological prioritisation to implicate four SNPs, in within the ADARB2 gene, two SNPs within the GRM7 gene and a single SNP in close proximity to a HTR7 gene. Association of SNPs within these neurotransmitter-related genes suggests a disrupted serotoninergic system that is perhaps specific to the Norfolk Island pedigree, but that might provide clues to understanding migraine more generally.

Functional analysis of missense variants in the TRESK (KCNK18) K+ channel

A loss of function mutation in the TRESK K2P potassium channel (KCNK18), has recently been linked with typical familial migraine with aura. We now report the functional characterisation of additional TRESK channel missense variants identified in unrelated patients. Several variants either had no apparent functional effect, or they caused a reduction in channel activity. However, the C110R variant was found to cause a complete loss of TRESK function, yet is present in both sporadic migraine and control cohorts, and no variation in KCNK18 copy number was found. Thus despite the previously identified association between loss of TRESK channel activity and migraine in a large multigenerational pedigree, this finding indicates that a single non-functional TRESK variant is not alone sufficient to cause typical migraine and highlights the genetic complexity of this disorder. Migraine is a common, disabling neurological disorder with a genetic, environmental and in some cases hormonal component. It is characterized by attacks of severe, usually unilateral and throbbing headache, can be accompanied by nausea, vomiting and photophobia and is clinically divided into two main subtypes, migraine with aura (MA) when a migraine is accompanied by transient and reversible focal neurological symptoms and migraine without aura (MO)1. The multifactorial and clinical heterogeneity of the disorder have considerably hindered the identification of common migraine susceptibility genes and most of our current understanding comes from the studies of familial hemiplegic migraine (FHM), a rare monogenic autosomal dominant form of MA2. So far, the three susceptibility genes that have been convincingly identified in FHM families all encode ion channels or transporters: CACNA1A encoding the α1 subunit of the Cav2.1 calcium channel3, SCN1A encoding the Nav1.1 sodium channel4 and ATP1A2 encoding the α2 subunit of the Na+/K+ pump5. It is believed that mutations in these genes may lead to increased efflux of glutamate and potassium in the synapse and thereby cause migraine by rendering the brain more susceptible to cortical spreading depression (CSD)6 which is thought to play a role in initiating a migraine attack7,8. However, these genes have not to date been implicated in common forms of migraine9. Nevertheless, current opinion suggests that typical migraine, like FHM, is also disorder of neuronal excitability, ion homeostasis and neurotransmitter release10,11,12. Mutations in the SLC4A4 gene encoding the sodium-bicarbonate cotransporter NBCe1, have recently been implicated in several different forms of migraine13, and a variety of genes involved in glutamate homeostasis (PGCP, MTDH14 and LRP115) and a cation channel (TRPM8)15 have also recently been implicated in migraine via genome-wide association studies. Ion channels are therefore highly likely to play an important role in the pathogenesis of typical migraine. TRESK (KCNK18), is a member of the two-pore domain (K2P) family of potassium channels involved in the control of cellular electrical excitability16. Regulation of TRESK activity by the calcium-dependent phosphatase calcineurin17, as well as its expression in dorsal root ganglia (DRG)18 and trigeminal ganglia (TG)19,20 has led to a proposed role for this channel in a variety of pain pathways. In a recent study, a frameshift mutation (F139Wfsx24) in TRESK was identified in a large multigenerational pedigree where it co-segregated perfectly with typical MA and a significant genome-wide linkage LOD score of 3.0. Furthermore, functional analysis revealed that this mutation caused a complete loss of TRESK function and that the truncated subunit was also capable of down regulating wild-type channel function. This therefore highlighted KCNK18 as potentially important candidate gene and suggested that TRESK dysfunction might play a possible role in the pathogenesis of familial migraine with visual aura20. Additional screening for KCNK18 mutations in unrelated sporadic migraine and control cohorts also identified a number of other missense variants; R10G, A34V, C110R, S231P and A233V20. The A233V variant was found only in the control cohort, whilst A34V was identified in a single Australian migraine proband for which family samples were not available, but it was not detected in controls. By contrast, the R10G, C110R, and S231P variants were found in both migraineurs and controls in both cohorts. In this study, we have investigated the functional effect of these variants to further probe the potential association of TRESK dysfunction with typical migraine.