Linking protein to phenotype with Mendelian Randomization detects 38 proteins with causal roles in human diseases and traits.

To efficiently transform genetic associations into drug targets requires evidence that a particular gene, and its encoded protein, contribute causally to a disease. To achieve this, we employ a three-step proteome-by-phenome Mendelian Randomization (MR) approach. In step one, 154 protein quantitative trait loci (pQTLs) were identified and independently replicated. From these pQTLs, 64 replicated locally-acting variants were used as instrumental variables for proteome-by-phenome MR across 846 traits (step two). When its assumptions are met, proteome-by-phenome MR, is equivalent to simultaneously running many randomized controlled trials. Step 2 yielded 38 proteins that significantly predicted variation in traits and diseases in 509 instances. Step 3 revealed that amongst the 271 instances from GeneAtlas (UK Biobank), 77 showed little evidence of pleiotropy (HEIDI), and 92 evidence of colocalization (eCAVIAR). Results were wide ranging: including, for example, new evidence for a causal role of tyrosine-protein phosphatase non-receptor type substrate 1 (SHPS1; SIRPA) in schizophrenia, and a new finding that intestinal fatty acid binding protein (FABP2) abundance contributes to the pathogenesis of cardiovascular disease. We also demonstrated confirmatory evidence for the causal role of four further proteins (FGF5, IL6R, LPL, LTA) in cardiovascular disease risk.

Mortality and cancer incidence in carriers of constitutional t(11;22)(q23;q11) translocations: A prospective study.

The constitutional t(11;22)(q23;q11) translocation is the only recurrent non-Robertsonian translocation known in humans. Carriers are phenotypically normal and are usually referred for cytogenetic testing because of multiple miscarriages, infertility, or having aneuploidy in offspring. A breast cancer predisposition has been suggested, but previous studies have been small and had methodological shortcomings. We therefore conducted a long-term prospective study of cancer and mortality risk in carriers. We followed 65 male and 101 female carriers of t(11;22)(q23;q11) diagnosed in cytogenetic laboratories in Britain during 1976-2005 for cancer and deaths for an average of 21.4 years per subject. Standardised mortality (SMR) and incidence (SIR) ratios were calculated comparing the numbers of observed events with those expected from national age-, sex-, country- and calendar-period-specific population rates. Cancer incidence was borderline significantly raised for cancer overall (SIR = 1.56, 95% CI: 0.98-2.36, n = 22), and significantly raised for invasive breast cancer (SIR = 2.74, 95% CI: 1.18-5.40, n = 8) and in situ breast cancer (SIR = 13.0, 95% CI: 3.55-33.4, n = 4). Breast cancer risks were particularly increased at ages <50 (SIR = 4.37, 95% CI: 1.42-10.2 for invasive, SIR = 22.8, 95% CI: 2.76-82.5 for in situ). Mortality was borderline significantly raised for breast cancer (SMR = 4.82, 95% CI: 0.99-14.1) but not significantly raised for other cancers or causes. Individuals diagnosed with t(11;22)(q23;q11) appear to be at several-fold increased breast cancer risk, with the greatest risks at premenopausal ages. Further research is required to understand the genetic mechanism involving 11q23 and 22q11 and there may be a need for enhanced breast cancer surveillance among female carriers.

Mortality and Cancer Incidence in Carriers of Balanced Robertsonian Translocations: A National Cohort Study.

A balanced robertsonian translocation (rob) results from fusion of 2 acrocentric chromosomes. Carriers are phenotypically normal and are often diagnosed because of recurrent miscarriages, infertility, or aneuploid offspring. Mortality and site-specific cancer risks in carriers have not been prospectively investigated. We followed 1,987 carriers diagnosed in Great Britain for deaths and cancer risk, over an average of 24.1 years. Standardized mortality and incidence ratios were calculated comparing the number of observed events against population rates. Overall mortality was higher for carriers diagnosed before age 15 years (standardized mortality ratio (SMR) = 2.00, 95% confidence interval (CI): 1.09, 3.35), similar for those diagnosed aged 15-44 years (SMR = 1.06, 95% CI: 0.86-1.28), and lower for those diagnosed aged 45-84 years (SMR = 0.81, 95% CI: 0.68, 0.95). Cancer incidence was higher for non-Hodgkin lymphoma (standardized incidence ratio (SIR) = 1.90, 95% CI: 1.01, 3.24) and childhood leukemia (SIR = 14.5, 95% CI: 1.75, 52.2), the latter particularly in rob(15;21) carriers (SIR = 447.8, 95% CI: 11.3, 2,495). Rob(13;14) carriers had a higher breast cancer risk (SIR = 1.58, 95% CI: 1.12, 2.15). Mortality risks relative to the population in diagnosed carriers depend on age at cytogenetic diagnosis, possibly reflecting age-specific cytogenetic referral reasons. Carriers might be at greater risk of childhood leukemia and non-Hodgkin lymphoma and those diagnosed with rob(13;14) of breast cancer.

Variegated yet non-random rod and cone photoreceptor disease patterns in RPGR-ORF15-associated retinal degeneration.

Mutations in the ORF15 exon of the RPGR gene cause a common form of X-linked retinitis pigmentosa, which often results in severe loss of vision. In dogs and mice, gene augmentation therapy has been shown to arrest the progressive degeneration of rod and cone photoreceptors. However, the distribution of potentially treatable photoreceptors across the human retinas and the rate of degeneration are not known. Here, we have defined structural and functional features of the disease in 70 individuals with ORF15 mutations. We also correlated the features observed in patients with those of three Rpgr-mutant (Rpgr-ko, Rd9, and Rpgr-cko) mice. In patients, there was pronounced macular disease. Across the retina, rod and cone dysfunction showed a range of patterns and a spectrum of severity between individuals, but a high symmetry was observed between eyes of each individual. Genotype was not related to disease expression. In the Rpgr-ko mice, there were intra-retinal differences in rhodopsin and cone opsin trafficking. In Rd9 and Rpgr-cko mice, retinal degeneration showed inter-ocular symmetry. Longitudinal results in patients revealed localized rod and cone dysfunction with progression rates of 0.8 to 1.3 log per decade in sensitivity loss. Relatively retained rod and cone photoreceptors in mid- and far-peripheral temporal-inferior and nasal-inferior visual field regions should be good targets for future localized gene therapies in patients.

The genomic, biochemical, and cellular responses of the retina in inherited photoreceptor degenerations and prospects for the treatment of these disorders.

The association of more than 140 genes with human photoreceptor degenerations, together with studies of animal models of these monogenic diseases, has provided great insight into their pathogenesis. Here we review the responses of the retina to photoreceptor mutations, including mechanisms of photoreceptor death. We discuss the roles of oxidative metabolism, mitochondrial reactive oxygen species, metabolic stress, protein misfolding, and defects in ciliary proteins, as well as the responses of Müller glia, microglia, and the retinal vasculature. Finally, we report on potential pharmacologic and biologic therapies, the critical role of histopathology as a prerequisite to treatment, and the exciting promise of gene therapy in animal models and in phase 1 trials in humans.

MR-PheWAS: exploring the causal effect of SUA level on multiple disease outcomes by using genetic instruments in UK Biobank.

OBJECTIVES: We aimed to investigate the role of serum uric acid (SUA) level in a broad spectrum of disease outcomes using data for 120 091 individuals from UK Biobank. METHODS: We performed a phenome-wide association study (PheWAS) to identify disease outcomes associated with SUA genetic risk loci. We then implemented conventional Mendelianrandomisation (MR) analysis to investigate the causal relevance between SUA level and disease outcomes identified from PheWAS. We next applied MR Egger analysis to detect and account for potential pleiotropy, which conventional MR analysis might mistake for causality, and used the HEIDI (heterogeneity in dependent instruments) test to remove cross-phenotype associations that were likely due to genetic linkage. RESULTS: Our PheWAS identified 25 disease groups/outcomes associated with SUA genetic risk loci after multiple testing correction (P<8.57e-05). Our conventional MR analysis implicated a causal role of SUA level in three disease groups: inflammatory polyarthropathies (OR=1.22, 95% CI 1.11 to 1.34), hypertensive disease (OR=1.08, 95% CI 1.03 to 1.14) and disorders of metabolism (OR=1.07, 95% CI 1.01 to 1.14); and four disease outcomes: gout (OR=4.88, 95% CI 3.91 to 6.09), essential hypertension (OR=1.08, 95% CI 1.03 to 1.14), myocardial infarction (OR=1.16, 95% CI 1.03 to 1.30) and coeliac disease (OR=1.41, 95% CI 1.05 to 1.89). After balancing pleiotropic effects in MR Egger analysis, only gout and its encompassing disease group of inflammatory polyarthropathies were considered to be causally associated with SUA level. Our analysis highlighted a locus (ATXN2/S2HB3) that may influence SUA level and multiple cardiovascular and autoimmune diseases via pleiotropy. CONCLUSIONS: Elevated SUA level is convincing to cause gout and inflammatory polyarthropathies, and might act as a marker for the wider range of diseases with which it associates. Our findings support further investigation on the clinical relevance of SUA level with cardiovascular, metabolic, autoimmune and respiratory diseases.

Genomic prediction of complex human traits: relatedness, trait architecture and predictive meta-models.

We explore the prediction of individuals' phenotypes for complex traits using genomic data. We compare several widely used prediction models, including Ridge Regression, LASSO and Elastic Nets estimated from cohort data, and polygenic risk scores constructed using published summary statistics from genome-wide association meta-analyses (GWAMA). We evaluate the interplay between relatedness, trait architecture and optimal marker density, by predicting height, body mass index (BMI) and high-density lipoprotein level (HDL) in two data cohorts, originating from Croatia and Scotland. We empirically demonstrate that dense models are better when all genetic effects are small (height and BMI) and target individuals are related to the training samples, while sparse models predict better in unrelated individuals and when some effects have moderate size (HDL). For HDL sparse models achieved good across-cohort prediction, performing similarly to the GWAMA risk score and to models trained within the same cohort, which indicates that, for predicting traits with moderately sized effects, large sample sizes and familial structure become less important, though still potentially useful. Finally, we propose a novel ensemble of whole-genome predictors with GWAMA risk scores and demonstrate that the resulting meta-model achieves higher prediction accuracy than either model on its own. We conclude that although current genomic predictors are not accurate enough for diagnostic purposes, performance can be improved without requiring access to large-scale individual-level data. Our methodologically simple meta-model is a means of performing predictive meta-analysis for optimizing genomic predictions and can be easily extended to incorporate multiple population-level summary statistics or other domain knowledge.

Homozygous loss-of-function variants in European cosmopolitan and isolate populations.

Homozygous loss of function (HLOF) variants provide a valuable window on gene function in humans, as well as an inventory of the human genes that are not essential for survival and reproduction. All humans carry at least a few HLOF variants, but the exact number of inactivated genes that can be tolerated is currently unknown­as are the phenotypic effects of losing function for most human genes. Here, we make use of 1432 whole exome sequences from five European populations to expand the catalogue of known human HLOF mutations; after stringent filtering of variants in our dataset, we identify a total of 173 HLOF mutations, 76 (44%) of which have not been observed previously. We find that population isolates are particularly well suited to surveys of novel HLOF genes because individuals in such populations carry extensive runs of homozygosity, which we show are enriched for novel, rare HLOF variants. Further, we make use of extensive phenotypic data to show that most HLOFs, ascertained in population-based samples, appear to have little detectable effect on the phenotype. On the contrary, we document several genes directly implicated in disease that seem to tolerate HLOF variants. Overall HLOF genes are enriched for olfactory receptor function and are expressed in testes more often than expected, consistent with reduced purifying selection and incipient pseudogenisation.

Rare and low-frequency variants and their association with plasma levels of fibrinogen, FVII, FVIII, and vWF.

Fibrinogen, coagulation factor VII (FVII), and factor VIII (FVIII) and its carrier von Willebrand factor (vWF) play key roles in hemostasis. Previously identified common variants explain only a small fraction of the trait heritabilities, and additional variations may be explained by associations with rarer variants with larger effects. The aim of this study was to identify low-frequency (minor allele frequency [MAF] ≥0.01 and <0.05) and rare (MAF <0.01) variants that influence plasma concentrations of these 4 hemostatic factors by meta-analyzing exome chip data from up to 76,000 participants of 4 ancestries. We identified 12 novel associations of low-frequency (n = 2) and rare (n = 10) variants across the fibrinogen, FVII, FVIII, and vWF traits that were independent of previously identified associations. Novel loci were found within previously reported genes and had effect sizes much larger than and independent of previously identified common variants. In addition, associations at KCNT1, HID1, and KATNB1 identified new candidate genes related to hemostasis for follow-up replication and functional genomic analysis. Newly identified low-frequency and rare-variant associations accounted for modest amounts of trait variance and therefore are unlikely to increase predicted trait heritability but provide new information for understanding individual variation in hemostasis pathways.

A general approach for haplotype phasing across the full spectrum of relatedness.

Many existing cohorts contain a range of relatedness between genotyped individuals, either by design or by chance. Haplotype estimation in such cohorts is a central step in many downstream analyses. Using genotypes from six cohorts from isolated populations and two cohorts from non-isolated populations, we have investigated the performance of different phasing methods designed for nominally 'unrelated' individuals. We find that SHAPEIT2 produces much lower switch error rates in all cohorts compared to other methods, including those designed specifically for isolated populations. In particular, when large amounts of IBD sharing is present, SHAPEIT2 infers close to perfect haplotypes. Based on these results we have developed a general strategy for phasing cohorts with any level of implicit or explicit relatedness between individuals. First SHAPEIT2 is run ignoring all explicit family information. We then apply a novel HMM method (duoHMM) to combine the SHAPEIT2 haplotypes with any family information to infer the inheritance pattern of each meiosis at all sites across each chromosome. This allows the correction of switch errors, detection of recombination events and genotyping errors. We show that the method detects numbers of recombination events that align very well with expectations based on genetic maps, and that it infers far fewer spurious recombination events than Merlin. The method can also detect genotyping errors and infer recombination events in otherwise uninformative families, such as trios and duos. The detected recombination events can be used in association scans for recombination phenotypes. The method provides a simple and unified approach to haplotype estimation, that will be of interest to researchers in the fields of human, animal and plant genetics.

Enrichment of pathogenic alleles in the brittle cornea gene, ZNF469, in keratoconus.

Keratoconus, a common inherited ocular disorder resulting in progressive corneal thinning, is the leading indication for corneal transplantation in the developed world. Genome-wide association studies have identified common SNPs 100 kb upstream of ZNF469 strongly associated with corneal thickness. Homozygous mutations in ZNF469 and PR domain-containing protein 5 (PRDM5) genes result in brittle cornea syndrome (BCS) Types 1 and 2, respectively. BCS is an autosomal recessive generalized connective tissue disorder associated with extreme corneal thinning and a high risk of corneal rupture. Some individuals with heterozygous PRDM5 mutations demonstrate a carrier ocular phenotype, which includes a mildly reduced corneal thickness, keratoconus and blue sclera. We hypothesized that heterozygous variants in PRDM5 and ZNF469 predispose to the development of isolated keratoconus. We found a significant enrichment of potentially pathologic heterozygous alleles in ZNF469 associated with the development of keratoconus (P = 0.00102) resulting in a relative risk of 12.0. This enrichment of rare potentially pathogenic alleles in ZNF469 in 12.5% of keratoconus patients represents a significant mutational load and highlights ZNF469 as the most significant genetic factor responsible for keratoconus identified to date.

Nine loci for ocular axial length identified through genome-wide association studies, including shared loci with refractive error.

Refractive errors are common eye disorders of public health importance worldwide. Ocular axial length (AL) is the major determinant of refraction and thus of myopia and hyperopia. We conducted a meta-analysis of genome-wide association studies for AL, combining 12,531 Europeans and 8,216 Asians. We identified eight genome-wide significant loci for AL (RSPO1, C3orf26, LAMA2, GJD2, ZNRF3, CD55, MIP, and ALPPL2) and confirmed one previously reported AL locus (ZC3H11B). Of the nine loci, five (LAMA2, GJD2, CD55, ALPPL2, and ZC3H11B) were associated with refraction in 18 independent cohorts (n = 23,591). Differential gene expression was observed for these loci in minus-lens-induced myopia mouse experiments and human ocular tissues. Two of the AL genes, RSPO1 and ZNRF3, are involved in Wnt signaling, a pathway playing a major role in the regulation of eyeball size. This study provides evidence of shared genes between AL and refraction, but importantly also suggests that these traits may have unique pathways.

Photoreceptor degeneration: genetic and mechanistic dissection of a complex trait.

The retina provides exquisitely sensitive vision that relies on the integrity of a uniquely vulnerable cell, the photoreceptor (PR). The genetic and mechanistic causes of retinal degeneration due to PR cell death--which occurs in conditions such as retinitis pigmentosa and age-related macular degeneration--are being successfully dissected. Over one hundred loci, some containing common variants but most containing rare variants, are implicated in the genetic architecture of this complex trait. This genetic heterogeneity results in equally diverse disease mechanisms that affect almost every aspect of PR function but converge on a common cell death pathway. Although genetic and mechanistic diversity creates challenges for therapy, some approaches--particularly gene-replacement therapy--are showing considerable promise.

Biological, clinical and population relevance of 95 loci for blood lipids.

Plasma concentrations of total cholesterol, low-density lipoprotein cholesterol, high-density lipoprotein cholesterol and triglycerides are among the most important risk factors for coronary artery disease (CAD) and are targets for therapeutic intervention. We screened the genome for common variants associated with plasma lipids in >100,000 individuals of European ancestry. Here we report 95 significantly associated loci (P < 5 x 10(-8)), with 59 showing genome-wide significant association with lipid traits for the first time. The newly reported associations include single nucleotide polymorphisms (SNPs) near known lipid regulators (for example, CYP7A1, NPC1L1 and SCARB1) as well as in scores of loci not previously implicated in lipoprotein metabolism. The 95 loci contribute not only to normal variation in lipid traits but also to extreme lipid phenotypes and have an impact on lipid traits in three non-European populations (East Asians, South Asians and African Americans). Our results identify several novel loci associated with plasma lipids that are also associated with CAD. Finally, we validated three of the novel genes-GALNT2, PPP1R3B and TTC39B-with experiments in mouse models. Taken together, our findings provide the foundation to develop a broader biological understanding of lipoprotein metabolism and to identify new therapeutic opportunities for the prevention of CAD.

Comparative performance of four methods for high-throughput glycosylation analysis of immunoglobulin G in genetic and epidemiological research.

The biological and clinical relevance of glycosylation is becoming increasingly recognized, leading to a growing interest in large-scale clinical and population-based studies. In the past few years, several methods for high-throughput analysis of glycans have been developed, but thorough validation and standardization of these methods is required before significant resources are invested in large-scale studies. In this study, we compared liquid chromatography, capillary gel electrophoresis, and two MS methods for quantitative profiling of N-glycosylation of IgG in the same data set of 1201 individuals. To evaluate the accuracy of the four methods we then performed analysis of association with genetic polymorphisms and age. Chromatographic methods with either fluorescent or MS-detection yielded slightly stronger associations than MS-only and multiplexed capillary gel electrophoresis, but at the expense of lower levels of throughput. Advantages and disadvantages of each method were identified, which should inform the selection of the most appropriate method in future studies.

Loci associated with N-glycosylation of human immunoglobulin G show pleiotropy with autoimmune diseases and haematological cancers.

Glycosylation of immunoglobulin G (IgG) influences IgG effector function by modulating binding to Fc receptors. To identify genetic loci associated with IgG glycosylation, we quantitated N-linked IgG glycans using two approaches. After isolating IgG from human plasma, we performed 77 quantitative measurements of N-glycosylation using ultra-performance liquid chromatography (UPLC) in 2,247 individuals from four European discovery populations. In parallel, we measured IgG N-glycans using MALDI-TOF mass spectrometry (MS) in a replication cohort of 1,848 Europeans. Meta-analysis of genome-wide association study (GWAS) results identified 9 genome-wide significant loci (P<2.27 × 10(-9)) in the discovery analysis and two of the same loci (B4GALT1 and MGAT3) in the replication cohort. Four loci contained genes encoding glycosyltransferases (ST6GAL1, B4GALT1, FUT8, and MGAT3), while the remaining 5 contained genes that have not been previously implicated in protein glycosylation (IKZF1, IL6ST-ANKRD55, ABCF2-SMARCD3, SUV420H1, and SMARCB1-DERL3). However, most of them have been strongly associated with autoimmune and inflammatory conditions (e.g., systemic lupus erythematosus, rheumatoid arthritis, ulcerative colitis, Crohn's disease, diabetes type 1, multiple sclerosis, Graves' disease, celiac disease, nodular sclerosis) and/or haematological cancers (acute lymphoblastic leukaemia, Hodgkin lymphoma, and multiple myeloma). Follow-up functional experiments in haplodeficient Ikzf1 knock-out mice showed the same general pattern of changes in IgG glycosylation as identified in the meta-analysis. As IKZF1 was associated with multiple IgG N-glycan traits, we explored biomarker potential of affected N-glycans in 101 cases with SLE and 183 matched controls and demonstrated substantial discriminative power in a ROC-curve analysis (area under the curve = 0.842). Our study shows that it is possible to identify new loci that control glycosylation of a single plasma protein using GWAS. The results may also provide an explanation for the reported pleiotropy and antagonistic effects of loci involved in autoimmune diseases and haematological cancer.

Meta-analysis of genome-wide association studies identifies six new Loci for serum calcium concentrations.

Calcium is vital to the normal functioning of multiple organ systems and its serum concentration is tightly regulated. Apart from CASR, the genes associated with serum calcium are largely unknown. We conducted a genome-wide association meta-analysis of 39,400 individuals from 17 population-based cohorts and investigated the 14 most strongly associated loci in ≤ 21,679 additional individuals. Seven loci (six new regions) in association with serum calcium were identified and replicated. Rs1570669 near CYP24A1 (P = 9.1E-12), rs10491003 upstream of GATA3 (P = 4.8E-09) and rs7481584 in CARS (P = 1.2E-10) implicate regions involved in Mendelian calcemic disorders: Rs1550532 in DGKD (P = 8.2E-11), also associated with bone density, and rs7336933 near DGKH/KIAA0564 (P = 9.1E-10) are near genes that encode distinct isoforms of diacylglycerol kinase. Rs780094 is in GCKR. We characterized the expression of these genes in gut, kidney, and bone, and demonstrate modulation of gene expression in bone in response to dietary calcium in mice. Our results shed new light on the genetics of calcium homeostasis.

ATR localizes to the photoreceptor connecting cilium and deficiency leads to severe photoreceptor degeneration in mice.

Ataxia-telangiectasia and Rad3 (ATR), a sensor of DNA damage, is associated with the regulation and control of cell division. ATR deficit is known to cause Seckel syndrome, characterized by severe proportionate short stature and microcephaly. We used a mouse model for Seckel disease to study the effect of ATR deficit on retinal development and function and we have found a new role for ATR, which is critical for the postnatal development of the photoreceptor (PR) layer in mouse retina. The structural and functional characterization of the ATR(+/s) mouse retinas displayed a specific, severe and early degeneration of rod and cone cells resembling some characteristics of human retinal degenerations. A new localization of ATR in the cilia of PRs and the fact that mutant mice have shorter cilia suggests that the PR degeneration here described results from a ciliary defect.

Genome-wide association study identifies genetic risk underlying primary rhegmatogenous retinal detachment.

Rhegmatogenous retinal detachment (RRD) is an important cause of vision loss and can potentially lead to blindness. The underlying pathogenesis is complex and incompletely understood. We applied a two-stage genetic association discovery phase followed by a replication phase in a combined total of 2833 RRD cases and 7871 controls. The discovery phase involved a genome-wide association scan of 867 affected individuals and 1953 controls from Scotland, followed by genotyping and testing 4347 highest ranking or candidate single nucleotide polymorphisms (SNPs) in independent sets of cases (1000) and controls (2912) of Dutch and British origin. None of the SNPs selected reached a Bonferroni-corrected threshold for significance (P < 1.27 × 10(-7)). The strongest association, for rs12960119 (P = 1.58 × 10(-7)) located within an intron of the SS18 gene. Further testing was carried out in independent case-control series from London (846 cases) and Croatia (120 cases). The combined meta-analysis identified one association reaching genome-wide significance for rs267738 (OR = 1.29, P = 2.11 × 10(-8)), a missense coding SNP and eQTL for CERS2 encoding the protein ceramide synthase 2. Several of the top signals showing suggestive significance in the combined meta-analysis encompassed genes with a documented role in cell adhesion or migration, including SS18, TIAM1, TSTA3 and LDB2, which warrant further investigation. This first genetic association study of RRD supports a polygenic component underlying RRD risk since 27.4% of the underlying RRD liability could be explained by the collective additive effects of the genotyped SNP from the discovery genome-wide scan.