Plasma proteomic associations with genetics and health in the UK Biobank.

The Pharma Proteomics Project is a precompetitive biopharmaceutical consortium characterizing the plasma proteomic profiles of 54,219 UK Biobank participants. Here we provide a detailed summary of this initiative, including technical and biological validations, insights into proteomic disease signatures, and prediction modelling for various demographic and health indicators. We present comprehensive protein quantitative trait locus (pQTL) mapping of 2,923 proteins that identifies 14,287 primary genetic associations, of which 81% are previously undescribed, alongside ancestry-specific pQTL mapping in non-European individuals. The study provides an updated characterization of the genetic architecture of the plasma proteome, contextualized with projected pQTL discovery rates as sample sizes and proteomic assay coverages increase over time. We offer extensive insights into trans pQTLs across multiple biological domains, highlight genetic influences on ligand-receptor interactions and pathway perturbations across a diverse collection of cytokines and complement networks, and illustrate long-range epistatic effects of ABO blood group and FUT2 secretor status on proteins with gastrointestinal tissue-enriched expression. We demonstrate the utility of these data for drug discovery by extending the genetic proxied effects of protein targets, such as PCSK9, on additional endpoints, and disentangle specific genes and proteins perturbed at loci associated with COVID-19 susceptibility. This public-private partnership provides the scientific community with an open-access proteomics resource of considerable breadth and depth to help to elucidate the biological mechanisms underlying proteo-genomic discoveries and accelerate the development of biomarkers, predictive models and therapeutics1.

Integrated analyses of growth differentiation factor-15 concentration and cardiometabolic diseases in humans.

Growth differentiation factor-15 (GDF15) is a stress response cytokine that is elevated in several cardiometabolic diseases and has attracted interest as a potential therapeutic target. To further explore the association of GDF15 with human disease, we conducted a broad study into the phenotypic and genetic correlates of GDF15 concentration in up to 14,099 individuals. Assessment of 772 traits across 6610 participants in FINRISK identified associations of GDF15 concentration with a range of phenotypes including all-cause mortality, cardiometabolic disease, respiratory diseases and psychiatric disorders, as well as inflammatory markers. A meta-analysis of genome-wide association studies (GWAS) of GDF15 concentration across three different assay platforms (n=14,099) confirmed significant heterogeneity due to a common missense variant (rs1058587; p.H202D) in GDF15, potentially due to epitope-binding artefacts. After conditioning on rs1058587, statistical fine mapping identified four independent putative causal signals at the locus. Mendelian randomisation (MR) analysis found evidence of a causal relationship between GDF15 concentration and high-density lipoprotein (HDL) but not body mass index (BMI). Using reverse MR, we identified a potential causal association of BMI on GDF15 (IVW pFDR = 0.0040). Taken together, our data derived from human population cohorts do not support a role for moderately elevated GDF15 concentrations as a causal factor in human cardiometabolic disease but support its role as a biomarker of metabolic stress.

Population-scale analysis of common and rare genetic variation associated with hearing loss in adults.

To better understand the genetics of hearing loss, we performed a genome-wide association meta-analysis with 125,749 cases and 469,497 controls across five cohorts. We identified 53/c loci affecting hearing loss risk, including common coding variants in COL9A3 and TMPRSS3. Through exome sequencing of 108,415 cases and 329,581 controls, we observed rare coding associations with 11 Mendelian hearing loss genes, including additive effects in known hearing loss genes GJB2 (Gly12fs; odds ratio [OR] = 1.21, P = 4.2 × 10-11) and SLC26A5 (gene burden; OR = 1.96, P = 2.8 × 10-17). We also identified hearing loss associations with rare coding variants in FSCN2 (OR = 1.14, P = 1.9 × 10-15) and KLHDC7B (OR = 2.14, P = 5.2 × 10-30). Our results suggest a shared etiology between Mendelian and common hearing loss in adults. This work illustrates the potential of large-scale exome sequencing to elucidate the genetic architecture of common disorders where both common and rare variation contribute to risk.

An effector index to predict target genes at GWAS loci.

Drug development and biological discovery require effective strategies to map existing genetic associations to causal genes. To approach this problem, we selected 12 common diseases and quantitative traits for which highly powered genome-wide association studies (GWAS) were available. For each disease or trait, we systematically curated positive control gene sets from Mendelian forms of the disease and from targets of medicines used for disease treatment. We found that these positive control genes were highly enriched in proximity of GWAS-associated single-nucleotide variants (SNVs). We then performed quantitative assessment of the contribution of commonly used genomic features, including open chromatin maps, expression quantitative trait loci (eQTL), and chromatin conformation data. Using these features, we trained and validated an Effector Index (Ei), to map target genes for these 12 common diseases and traits. Ei demonstrated high predictive performance, both with cross-validation on the training set, and an independently derived set for type 2 diabetes. Key predictive features included coding or transcript-altering SNVs, distance to gene, and open chromatin-based metrics. This work outlines a simple, understandable approach to prioritize genes at GWAS loci for functional follow-up and drug development, and provides a systematic strategy for prioritization of GWAS target genes.

Genome-wide analysis of 102,084 migraine cases identifies 123 risk loci and subtype-specific risk alleles.

Migraine affects over a billion individuals worldwide but its genetic underpinning remains largely unknown. Here, we performed a genome-wide association study of 102,084 migraine cases and 771,257 controls and identified 123 loci, of which 86 are previously unknown. These loci provide an opportunity to evaluate shared and distinct genetic components in the two main migraine subtypes: migraine with aura and migraine without aura. Stratification of the risk loci using 29,679 cases with subtype information indicated three risk variants that seem specific for migraine with aura (in HMOX2, CACNA1A and MPPED2), two that seem specific for migraine without aura (near SPINK2 and near FECH) and nine that increase susceptibility for migraine regardless of subtype. The new risk loci include genes encoding recent migraine-specific drug targets, namely calcitonin gene-related peptide (CALCA/CALCB) and serotonin 1F receptor (HTR1F). Overall, genomic annotations among migraine-associated variants were enriched in both vascular and central nervous system tissue/cell types, supporting unequivocally that neurovascular mechanisms underlie migraine pathophysiology.

Exome sequencing and analysis of 454,787 UK Biobank participants.

A major goal in human genetics is to use natural variation to understand the phenotypic consequences of altering each protein-coding gene in the genome. Here we used exome sequencing1 to explore protein-altering variants and their consequences in 454,787 participants in the UK Biobank study2. We identified 12 million coding variants, including around 1 million loss-of-function and around 1.8 million deleterious missense variants. When these were tested for association with 3,994 health-related traits, we found 564 genes with trait associations at P ≤ 2.18 × 10-11. Rare variant associations were enriched in loci from genome-wide association studies (GWAS), but most (91%) were independent of common variant signals. We discovered several risk-increasing associations with traits related to liver disease, eye disease and cancer, among others, as well as risk-lowering associations for hypertension (SLC9A3R2), diabetes (MAP3K15, FAM234A) and asthma (SLC27A3). Six genes were associated with brain imaging phenotypes, including two involved in neural development (GBE1, PLD1). Of the signals available and powered for replication in an independent cohort, 81% were confirmed; furthermore, association signals were generally consistent across individuals of European, Asian and African ancestry. We illustrate the ability of exome sequencing to identify gene-trait associations, elucidate gene function and pinpoint effector genes that underlie GWAS signals at scale.

Sequencing of 640,000 exomes identifies GPR75 variants associated with protection from obesity.

Large-scale human exome sequencing can identify rare protein-coding variants with a large impact on complex traits such as body adiposity. We sequenced the exomes of 645,626 individuals from the United Kingdom, the United States, and Mexico and estimated associations of rare coding variants with body mass index (BMI). We identified 16 genes with an exome-wide significant association with BMI, including those encoding five brain-expressed G protein-coupled receptors (CALCR, MC4R, GIPR, GPR151, and GPR75). Protein-truncating variants in GPR75 were observed in ~4/10,000 sequenced individuals and were associated with 1.8 kilograms per square meter lower BMI and 54% lower odds of obesity in the heterozygous state. Knock out of Gpr75 in mice resulted in resistance to weight gain and improved glycemic control in a high-fat diet model. Inhibition of GPR75 may provide a therapeutic strategy for obesity.

Computationally efficient whole-genome regression for quantitative and binary traits.

Genome-wide association analysis of cohorts with thousands of phenotypes is computationally expensive, particularly when accounting for sample relatedness or population structure. Here we present a novel machine-learning method called REGENIE for fitting a whole-genome regression model for quantitative and binary phenotypes that is substantially faster than alternatives in multi-trait analyses while maintaining statistical efficiency. The method naturally accommodates parallel analysis of multiple phenotypes and requires only local segments of the genotype matrix to be loaded in memory, in contrast to existing alternatives, which must load genome-wide matrices into memory. This results in substantial savings in compute time and memory usage. We introduce a fast, approximate Firth logistic regression test for unbalanced case-control phenotypes. The method is ideally suited to take advantage of distributed computing frameworks. We demonstrate the accuracy and computational benefits of this approach using the UK Biobank dataset with up to 407,746 individuals.

Genetics of 35 blood and urine biomarkers in the UK Biobank.

Clinical laboratory tests are a critical component of the continuum of care. We evaluate the genetic basis of 35 blood and urine laboratory measurements in the UK Biobank (n = 363,228 individuals). We identify 1,857 loci associated with at least one trait, containing 3,374 fine-mapped associations and additional sets of large-effect (>0.1 s.d.) protein-altering, human leukocyte antigen (HLA) and copy number variant (CNV) associations. Through Mendelian randomization (MR) analysis, we discover 51 causal relationships, including previously known agonistic effects of urate on gout and cystatin C on stroke. Finally, we develop polygenic risk scores (PRSs) for each biomarker and build 'multi-PRS' models for diseases using 35 PRSs simultaneously, which improved chronic kidney disease, type 2 diabetes, gout and alcoholic cirrhosis genetic risk stratification in an independent dataset (FinnGen; n = 135,500) relative to single-disease PRSs. Together, our results delineate the genetic basis of biomarkers and their causal influences on diseases and improve genetic risk stratification for common diseases.

Genome-wide association meta-analysis of nicotine metabolism and cigarette consumption measures in smokers of European descent.

Smoking behaviors, including amount smoked, smoking cessation, and tobacco-related diseases, are altered by the rate of nicotine clearance. Nicotine clearance can be estimated using the nicotine metabolite ratio (NMR) (ratio of 3'hydroxycotinine/cotinine), but only in current smokers. Advancing the genomics of this highly heritable biomarker of CYP2A6, the main metabolic enzyme for nicotine, will also enable investigation of never and former smokers. We performed the largest genome-wide association study (GWAS) to date of the NMR in European ancestry current smokers (n = 5185), found 1255 genome-wide significant variants, and replicated the chromosome 19 locus. Fine-mapping of chromosome 19 revealed 13 putatively causal variants, with nine of these being highly putatively causal and mapping to CYP2A6, MAP3K10, ADCK4, and CYP2B6. We also identified a putatively causal variant on chromosome 4 mapping to TMPRSS11E and demonstrated an association between TMPRSS11E variation and a UGT2B17 activity phenotype. Together the 14 putatively causal SNPs explained ~38% of NMR variation, a substantial increase from the ~20 to 30% previously explained. Our additional GWASs of nicotine intake biomarkers showed that cotinine and smoking intensity (cotinine/cigarettes per day (CPD)) shared chromosome 19 and chromosome 4 loci with the NMR, and that cotinine and a more accurate biomarker, cotinine + 3'hydroxycotinine, shared a chromosome 15 locus near CHRNA5 with CPD and Pack-Years (i.e., cumulative exposure). Understanding the genetic factors influencing smoking-related traits facilitates epidemiological studies of smoking and disease, as well as assists in optimizing smoking cessation support, which in turn will reduce the enormous personal and societal costs associated with smoking.

An expanded analysis framework for multivariate GWAS connects inflammatory biomarkers to functional variants and disease.

Multivariate methods are known to increase the statistical power to detect associations in the case of shared genetic basis between phenotypes. They have, however, lacked essential analytic tools to follow-up and understand the biology underlying these associations. We developed a novel computational workflow for multivariate GWAS follow-up analyses, including fine-mapping and identification of the subset of traits driving associations (driver traits). Many follow-up tools require univariate regression coefficients which are lacking from multivariate results. Our method overcomes this problem by using Canonical Correlation Analysis to turn each multivariate association into its optimal univariate Linear Combination Phenotype (LCP). This enables an LCP-GWAS, which in turn generates the statistics required for follow-up analyses. We implemented our method on 12 highly correlated inflammatory biomarkers in a Finnish population-based study. Altogether, we identified 11 associations, four of which (F5, ABO, C1orf140 and PDGFRB) were not detected by biomarker-specific analyses. Fine-mapping identified 19 signals within the 11 loci and driver trait analysis determined the traits contributing to the associations. A phenome-wide association study on the 19 representative variants from the signals in 176,899 individuals from the FinnGen study revealed 53 disease associations (p < 1 × 10-4). Several reported pQTLs in the 11 loci provided orthogonal evidence for the biologically relevant functions of the representative variants. Our novel multivariate analysis workflow provides a powerful addition to standard univariate GWAS analyses by enabling multivariate GWAS follow-up and thus promoting the advancement of powerful multivariate methods in genomics.

Functionally informed fine-mapping and polygenic localization of complex trait heritability.

Fine-mapping aims to identify causal variants impacting complex traits. We propose PolyFun, a computationally scalable framework to improve fine-mapping accuracy by leveraging functional annotations across the entire genome-not just genome-wide-significant loci-to specify prior probabilities for fine-mapping methods such as SuSiE or FINEMAP. In simulations, PolyFun + SuSiE and PolyFun + FINEMAP were well calibrated and identified >20% more variants with a posterior causal probability >0.95 than identified in their nonfunctionally informed counterparts. In analyses of 49 UK Biobank traits (average n = 318,000), PolyFun + SuSiE identified 3,025 fine-mapped variant-trait pairs with posterior causal probability >0.95, a >32% improvement versus SuSiE. We used posterior mean per-SNP heritabilities from PolyFun + SuSiE to perform polygenic localization, constructing minimal sets of common SNPs causally explaining 50% of common SNP heritability; these sets ranged in size from 28 (hair color) to 3,400 (height) to 2 million (number of children). In conclusion, PolyFun prioritizes variants for functional follow-up and provides insights into complex trait architectures.

Polygenic Hyperlipidemias and Coronary Artery Disease Risk.

BACKGROUND: Hyperlipidemia is a highly heritable risk factor for coronary artery disease (CAD). While monogenic familial hypercholesterolemia associates with severely increased CAD risk, it remains less clear to what extent a high polygenic load of a large number of LDL (low-density lipoprotein) cholesterol (LDL-C) or triglyceride (TG)-increasing variants associates with increased CAD risk. METHODS: We derived polygenic risk scores (PRSs) with ≈6M variants separately for LDL-C and TG with weights from a UK Biobank-based genome-wide association study with ≈324K samples. We evaluated the impact of polygenic hypercholesterolemia and hypertriglyceridemia to lipid levels in 27 039 individuals from the National FINRISK Study (FINRISK) cohort and to CAD risk in 135 638 individuals (13 753 CAD cases) from the FinnGen project (FinnGen). RESULTS: In FINRISK, median LDL-C was 3.39 (95% CI, 3.38-3.40) mmol/L, and it ranged from 2.87 (95% CI, 2.82-2.94) to 3.78 (95% CI, 3.71-3.83) mmol/L between the lowest and highest 5% of the LDL-C PRS distribution. Median TG was 1.19 (95% CI, 1.18-1.20) mmol/L, ranging from 0.97 (95% CI, 0.94-1.00) to 1.55 (95% CI, 1.48-1.61) mmol/L with the TG PRS. In FinnGen, comparing the highest 5% of the PRS to the lowest 95%, CAD odds ratio was 1.36 (95% CI, 1.24-1.49) for the LDL-C PRS and 1.31 (95% CI, 1.19-1.43) for the TG PRS. These estimates were only slightly attenuated when adjusting for a CAD PRS (odds ratio, 1.26 [95% CI, 1.16-1.38] for LDL-C and 1.24 [95% CI, 1.13-1.36] for TG PRS). CONCLUSIONS: The CAD risk associated with a high polygenic load for lipid-increasing variants was proportional to their impact on lipid levels and partially overlapping with a CAD PRS. In contrast with a PRS for CAD, the lipid PRSs point to known and directly modifiable risk factors providing additional guidance for clinical translation.

Interrogation of human hematopoiesis at single-cell and single-variant resolution.

Widespread linkage disequilibrium and incomplete annotation of cell-to-cell state variation represent substantial challenges to elucidating mechanisms of trait-associated genetic variation. Here we perform genetic fine-mapping for blood cell traits in the UK Biobank to identify putative causal variants. These variants are enriched in genes encoding proteins in trait-relevant biological pathways and in accessible chromatin of hematopoietic progenitors. For regulatory variants, we explore patterns of developmental enhancer activity, predict molecular mechanisms, and identify likely target genes. In several instances, we localize multiple independent variants to the same regulatory element or gene. We further observe that variants with pleiotropic effects preferentially act in common progenitor populations to direct the production of distinct lineages. Finally, we leverage fine-mapped variants in conjunction with continuous epigenomic annotations to identify trait-cell type enrichments within closely related populations and in single cells. Our study provides a comprehensive framework for single-variant and single-cell analyses of genetic associations.

Non-operative meningiomas: long-term follow-up of 136 patients.

BACKGROUND: Improving access to neuroradiology investigations has led to an increased rate of diagnosis of incidental meningiomas. METHOD: A cohort of 136 incidental meningioma patients collected by a single neurosurgeon in a single neurosurgical centre is retrospectively analysed between 2002 and 2016. Demographic data, imaging and clinical features are presented. The radiological factors associated with meningiomas progression are also presented. RESULTS: The mean age at diagnosis was 65 (range, 33-94) years. Univariate analysis showed oedema was most strongly correlated with progression (p = 0.010) followed by hyperintensity in T2-weighted (T2W) MRI (p = 0.029) and in Flair-T2W MRI (p = 0.017). Isointensity in Flair-T2W MRI (0.004) was most strongly correlated with non-progression of the meningioma followed by calcification (p = 0.007), older age (p = 0.087), hypointensity in Flair-T2W MRI (p = 0.014) sequences and in T2W MRI (p = 0.096). In multivariate analysis, the strongest radiological factor predictive of progression was peritumoural oedema (p = 0.016) and that of non-progression was calcification (p = 0.002). At the end of the median follow-up (FU) of 43 (range, 4-150) months, 109 (80%) patients remained clinically stable, 13 (10%) became symptomatic and 14 (10%) showed clinical and radiological progression. CONCLUSIONS: One hundred and nine (80%) patients remained stable at the end of FU. Peritumoural oedema was predictive of meningiomas progression. Further prospective study is needed to identify the combination of factors which can predict the meningioma progression for an early surgery or early discharge.

Prospects of Fine-Mapping Trait-Associated Genomic Regions by Using Summary Statistics from Genome-wide Association Studies.

During the past few years, various novel statistical methods have been developed for fine-mapping with the use of summary statistics from genome-wide association studies (GWASs). Although these approaches require information about the linkage disequilibrium (LD) between variants, there has not been a comprehensive evaluation of how estimation of the LD structure from reference genotype panels performs in comparison with that from the original individual-level GWAS data. Using population genotype data from Finland and the UK Biobank, we show here that a reference panel of 1,000 individuals from the target population is adequate for a GWAS cohort of up to 10,000 individuals, whereas smaller panels, such as those from the 1000 Genomes Project, should be avoided. We also show, both theoretically and empirically, that the size of the reference panel needs to scale with the GWAS sample size; this has important consequences for the application of these methods in ongoing GWAS meta-analyses and large biobank studies. We conclude by providing software tools and by recommending practices for sharing LD information to more efficiently exploit summary statistics in genetics research.

biMM: efficient estimation of genetic variances and covariances for cohorts with high-dimensional phenotype measurements.

SUMMARY: Genetic research utilizes a decomposition of trait variances and covariances into genetic and environmental parts. Our software package biMM is a computationally efficient implementation of a bivariate linear mixed model for settings where hundreds of traits have been measured on partially overlapping sets of individuals. AVAILABILITY AND IMPLEMENTATION: Implementation in R freely available at www.iki.fi/mpirinen . CONTACT: matti.pirinen@helsinki.fi. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genome-Wide Meta-Analysis of Sciatica in Finnish Population.

Sciatica or the sciatic syndrome is a common and often disabling low back disorder in the working-age population. It has a relatively high heritability but poorly understood molecular mechanisms. The Finnish population is a genetic isolate where small founder population and bottleneck events have led to enrichment of certain rare and low frequency variants. We performed here the first genome-wide association (GWAS) and meta-analysis of sciatica. The meta-analysis was conducted across two GWAS covering 291 Finnish sciatica cases and 3671 controls genotyped and imputed at 7.7 million autosomal variants. The most promising loci (p<1x10-6) were replicated in 776 Finnish sciatica patients and 18,489 controls. We identified five intragenic variants, with relatively low frequencies, at two novel loci associated with sciatica at genome-wide significance. These included chr9:14344410:I (rs71321981) at 9p22.3 (NFIB gene; p = 1.30x10-8, MAF = 0.08) and four variants at 15q21.2: rs145901849, rs80035109, rs190200374 and rs117458827 (MYO5A; p = 1.34x10-8, MAF = 0.06; p = 2.32x10-8, MAF = 0.07; p = 3.85x10-8, MAF = 0.06; p = 4.78x10-8, MAF = 0.07, respectively). The most significant association in the meta-analysis, a single base insertion rs71321981 within the regulatory region of the transcription factor NFIB, replicated in an independent Finnish population sample (p = 0.04). Despite identifying 15q21.2 as a promising locus, we were not able to replicate it. It was differentiated; the lead variants within 15q21.2 were more frequent in Finland (6-7%) than in other European populations (1-2%). Imputation accuracies of the three significantly associated variants (chr9:14344410:I, rs190200374, and rs80035109) were validated by genotyping. In summary, our results suggest a novel locus, 9p22.3 (NFIB), which may be involved in susceptibility to sciatica. In addition, another locus, 15q21.2, emerged as a promising one, but failed to replicate.

Reduced genomic tumor heterogeneity after neoadjuvant chemotherapy is related to favorable outcome in patients with esophageal adenocarcinoma.

Neoadjuvant chemo(radio)therapy followed by surgery is the standard of care for patients with locally advanced resectable esophageal adenocarcinoma (EAC). There is increasing evidence that drug resistance might be related to genomic heterogeneity. We investigated whether genomic tumor heterogeneity is different after cytotoxic chemotherapy and is associated with EAC patient survival. We used arrayCGH and a quantitative assessment of the whole genome DNA copy number aberration patterns ('DNA copy number entropy') to establish the level of genomic tumor heterogeneity in 80 EAC treated with neoadjuvant chemotherapy followed by surgery (CS group) or surgery alone (S group). The association between DNA copy number entropy, clinicopathological variables and survival was investigated.DNA copy number entropy was reduced after chemotherapy, even if there was no morphological evidence of response to therapy (p<0.001). Low DNA copy number entropy was associated with improved survival in the CS group (p=0.011) but not in the S group (p=0.396).Our results suggest that cytotoxic chemotherapy reduces DNA copy number entropy, which might be a more sensitive tumor response marker than changes in the morphological tumor phenotype. The use of DNA copy number entropy in clinical practice will require validation of our results in a prospective study.