Large-scale circulating proteome association study (CPAS) meta-analysis identifies circulating proteins and pathways predicting incident hip fractures

Hip fractures are associated with significant disability, high cost, and mortality. However, the exact biological mechanisms underlying susceptibility to hip fractures remain incompletely understood. In an exploratory search of the underlying biology as reflected through the circulating proteome, we performed a comprehensive Circulating Proteome Association Study (CPAS) meta-analysis for incident hip fractures. Analyses included 6430 subjects from two prospective cohort studies (Cardiovascular Health Study and Trøndelag Health Study) with circulating proteomics data (aptamer-based 5 K SomaScan version 4.0 assay; 4979 aptamers). Associations between circulating protein levels and incident hip fractures were estimated for each cohort using age and sex-adjusted Cox regression models. Participants experienced 643 incident hip fractures. Compared with the individual studies, inverse-variance weighted meta-analyses yielded more statistically significant associations, identifying 23 aptamers associated with incident hip fractures (conservative Bonferroni correction 0.05/4979, P < 1.0 × 10-5). The aptamers most strongly associated with hip fracture risk corresponded to two proteins of the growth hormone/insulin growth factor system (GHR and IGFBP2), as well as GDF15 and EGFR. High levels of several inflammation-related proteins (CD14, CXCL12, MMP12, ITIH3) were also associated with increased hip fracture risk. Ingenuity pathway analysis identified reduced LXR/RXR activation and increased acute phase response signaling to be overrepresented among those proteins associated with increased hip fracture risk. These analyses identified several circulating proteins and pathways consistently associated with incident hip fractures. These findings underscore the usefulness of the meta-analytic approach for comprehensive CPAS in a similar manner as has previously been observed for large-scale human genetic studies. Future studies should investigate the underlying biology of these potential novel drug targets.

Proportion of venous thromboembolism attributed to recognized prothrombotic genotypes in men and women

Background

Data on the proportion of venous thromboembolism (VTE) risk attributed to prothrombotic genotypes in men and women are limited.

Objectives

We aimed to estimate the population attributable fraction (PAF) of VTE for recognized, common prothrombotic genotypes in men and women using a population-based case cohort.

Methods

Cases with incident VTE (n = 1493) and a randomly sampled subcohort (n = 13,069) were derived from the Tromsø study (1994-2012) and the Trøndelag Health Study (1995-2008) cohorts. DNA samples were genotyped for 17 single-nucleotide polymorphisms (SNPs) previously associated with VTE. PAFs with 95% bias-corrected CIs (based on 10,000 bootstrap samples) were estimated for SNPs significantly associated with VTE, and a 6-SNP cumulative model was constructed for both sexes.

Results

In women, the individual PAFs for SNPs included in the cumulative model were 16.9% for ABO (rs8176719), 17.6% for F11 (rs2036914), 15.1% for F11 (rs2289252), 8.7% for FVL (rs6025), 6.0% for FGG (rs2066865), and 0.2% for F2 (rs1799963). The cumulative PAF for this 6-SNP model was 37.8%. In men, the individual PAFs for SNPs included in the cumulative model were 21.3% for ABO, 12.2% for F11 (rs2036914), 10.4% for F11 (rs2289252), 7.5% for FVL, 7.8% for FGG, and 1.1% for F2. This resulted in a cumulative PAF in men of 51.9%.

Conclusion

Our findings in a Norwegian population suggest that 52% and 38% of the VTEs can be attributed to known prothrombotic genotypes in men and women, respectively.

An atlas of genetic determinants of forearm fracture

Osteoporotic fracture is among the most common and costly of diseases. While reasonably heritable, its genetic determinants have remained elusive. Forearm fractures are the most common clinically recognized osteoporotic fractures with a relatively high heritability. To establish an atlas of the genetic determinants of forearm fractures, we performed genome-wide association analyses including 100,026 forearm fracture cases. We identified 43 loci, including 26 new fracture loci. Although most fracture loci associated with bone mineral density, we also identified loci that primarily regulate bone quality parameters. Functional studies of one such locus, at TAC4, revealed that Tac4-/- mice have reduced mechanical bone strength. The strongest forearm fracture signal, at WNT16, displayed remarkable bone-site-specificity with no association with hip fractures. Tall stature and low body mass index were identified as new causal risk factors for fractures. The insights from this atlas may improve fracture prediction and enable therapeutic development to prevent fractures.

Cluster Headache Genomewide Association Study and Meta-Analysis Identifies Eight Loci and Implicates Smoking as Causal Risk Factor

OBJECTIVE

The objective of this study was to aggregate data for the first genomewide association study meta-analysis of cluster headache, to identify genetic risk variants, and gain biological insights.

METHODS

A total of 4,777 cases (3,348 men and 1,429 women) with clinically diagnosed cluster headache were recruited from 10 European and 1 East Asian cohorts. We first performed an inverse-variance genomewide association meta-analysis of 4,043 cases and 21,729 controls of European ancestry. In a secondary trans-ancestry meta-analysis, we included 734 cases and 9,846 controls of East Asian ancestry. Candidate causal genes were prioritized by 5 complementary methods: expression quantitative trait loci, transcriptome-wide association, fine-mapping of causal gene sets, genetically driven DNA methylation, and effects on protein structure. Gene set and tissue enrichment analyses, genetic correlation, genetic risk score analysis, and Mendelian randomization were part of the downstream analyses.

RESULTS

The estimated single nucleotide polymorphism (SNP)-based heritability of cluster headache was 14.5%. We identified 9 independent signals in 7 genomewide significant loci in the primary meta-analysis, and one additional locus in the trans-ethnic meta-analysis. Five of the loci were previously known. The 20 genes prioritized as potentially causal for cluster headache showed enrichment to artery and brain tissue. Cluster headache was genetically correlated with cigarette smoking, risk-taking behavior, attention deficit hyperactivity disorder (ADHD), depression, and musculoskeletal pain. Mendelian randomization analysis indicated a causal effect of cigarette smoking intensity on cluster headache. Three of the identified loci were shared with migraine.

INTERPRETATION

This first genomewide association study meta-analysis gives clues to the biological basis of cluster headache and indicates that smoking is a causal risk factor. ANN NEUROL 2023;94:713-726.

Multi-ancestry meta-analysis identifies 5 novel loci for ischemic stroke and reveals heterogeneity of effects between sexes and ancestries

Stroke is the second leading cause of death and disability worldwide. Stroke prevalence varies by sex and ancestry, possibly due to genetic heterogeneity between subgroups. We performed a genome-wide meta-analysis of 16 biobanks across multiple ancestries to study the genetics of ischemic stroke (60,176 cases, 1,310,725 controls) as part of the Global Biobank Meta-analysis Initiative (GBMI) and further combined the results with previously published MegaStroke. Five novel loci for ischemic stroke (LAMC1, CALCRL, PLSCR1, CDKN1A, and SWAP70) were identified after replication in four additional datasets. One previously reported locus showed significant ancestry heterogeneity (ABO), and one showed significant sex heterogeneity (ALDH2). The ALDH2 association was male specific (males p = 1.67e-24, females p = 0.126) and was additionally observed only in the East Asian ancestry (male) samples. These findings emphasize the need for more diverse datasets with large sample sizes to further understand the genetic predisposition of stroke in different ancestry and sex groups.

Genetic effects on the timing of parturition and links to fetal birth weight

The timing of parturition is crucial for neonatal survival and infant health. Yet, its genetic basis remains largely unresolved. We present a maternal genome-wide meta-analysis of gestational duration (n = 195,555), identifying 22 associated loci (24 independent variants) and an enrichment in genes differentially expressed during labor. A meta-analysis of preterm delivery (18,797 cases, 260,246 controls) revealed six associated loci and large genetic similarities with gestational duration. Analysis of the parental transmitted and nontransmitted alleles (n = 136,833) shows that 15 of the gestational duration genetic variants act through the maternal genome, whereas 7 act both through the maternal and fetal genomes and 2 act only via the fetal genome. Finally, the maternal effects on gestational duration show signs of antagonistic pleiotropy with the fetal effects on birth weight: maternal alleles that increase gestational duration have negative fetal effects on birth weight. The present study provides insights into the genetic effects on the timing of parturition and the complex maternal-fetal relationship between gestational duration and birth weight.

Genetic diversity fuels gene discovery for tobacco and alcohol use

Tobacco and alcohol use are heritable behaviours associated with 15% and 5.3% of worldwide deaths, respectively, due largely to broad increased risk for disease and injury1-4. These substances are used across the globe, yet genome-wide association studies have focused largely on individuals of European ancestries5. Here we leveraged global genetic diversity across 3.4 million individuals from four major clines of global ancestry (approximately 21% non-European) to power the discovery and fine-mapping of genomic loci associated with tobacco and alcohol use, to inform function of these loci via ancestry-aware transcriptome-wide association studies, and to evaluate the genetic architecture and predictive power of polygenic risk within and across populations. We found that increases in sample size and genetic diversity improved locus identification and fine-mapping resolution, and that a large majority of the 3,823 associated variants (from 2,143 loci) showed consistent effect sizes across ancestry dimensions. However, polygenic risk scores developed in one ancestry performed poorly in others, highlighting the continued need to increase sample sizes of diverse ancestries to realize any potential benefit of polygenic prediction.

Assessment of the genetic and clinical determinants of hip fracture risk: Genome-wide association and Mendelian randomization study

Hip fracture is the clinically most important fracture, but the genetic architecture of hip fracture is unclear. Here, we perform a large-scale hip fracture genome-wide association study meta-analysis and Mendelian randomization study using five cohorts from European biobanks. The results show that five genetic signals associate with hip fractures. Among these, one signal associates with falls, but not with bone mineral density (BMD), while four signals are in loci known to be involved in bone biology. Mendelian randomization analyses demonstrate a strong causal effect of decreased femoral neck BMD and moderate causal effects of Alzheimer's disease and having ever smoked regularly on risk of hip fractures. The substantial causal effect of decreased femoral neck BMD on hip fractures in both young and old subjects and in both men and women supports the use of change in femoral neck BMD as a surrogate outcome for hip fractures in clinical trials.

Genome-wide risk prediction of common diseases across ancestries in one million people

Polygenic risk scores (PRS) measure genetic disease susceptibility by combining risk effects across the genome. For coronary artery disease (CAD), type 2 diabetes (T2D), and breast and prostate cancer, we performed cross-ancestry evaluation of genome-wide PRSs in six biobanks in Europe, the United States, and Asia. We studied transferability of these highly polygenic, genome-wide PRSs across global ancestries, within European populations with different health-care systems, and local population substructures in a population isolate. All four PRSs had similar accuracy across European and Asian populations, with poorer transferability in the smaller group of individuals of African ancestry. The PRSs had highly similar effect sizes in different populations of European ancestry, and in early- and late-settlement regions with different recent population bottlenecks in Finland. Comparing genome-wide PRSs to PRSs containing a smaller number of variants, the highly polygenic, genome-wide PRSs generally displayed higher effect sizes and better transferability across global ancestries. Our findings indicate that in the populations investigated, the current genome-wide polygenic scores for common diseases have potential for clinical utility within different health-care settings for individuals of European ancestry, but that the utility in individuals of African ancestry is currently much lower.

Genome-Wide Association Study of 2,093 Cases With Idiopathic Polyneuropathy and 445,256 Controls Identifies First Susceptibility Loci

Background: About one third of patients with chronic polyneuropathy have no obvious underlying etiology and are classified as having idiopathic polyneuropathy. The lack of knowledge about pathomechanisms and predisposing factors limits the development of effective prevention and treatment for these patients. We report the first genome-wide association study (GWAS) of idiopathic polyneuropathy.

Methods: Cases with idiopathic polyneuropathy and healthy controls were identified by linkage to hospital records. We performed genome-wide association studies using genetic data from two large population-based health studies, the Trøndelag Health Study (HUNT, 1,147 cases and 62,204 controls) and UK Biobank (UKB, 946 cases and 383,052 controls). In a two-stage analysis design, we first treated HUNT as a discovery cohort and UK Biobank as a replication cohort. Secondly, we combined the two studies in a meta-analysis. Downstream analyses included genetic correlation to other traits and diseases. We specifically examined previously reported risk loci, and genes known to cause hereditary polyneuropathy.

Results: No replicable risk loci were identified in the discovery-replication stage, in line with the limited predicted power of this approach. When combined in a meta-analysis, two independent loci reached statistical significance (rs7294354 in B4GALNT3, P-value 4.51 × 10⁻⁸) and (rs147738081 near NR5A2, P-value 4.75 × 10⁻⁸). Idiopathic polyneuropathy genetically correlated with several anthropometric measures, most pronounced for height, and with several pain-related traits.

Conclusions: In this first GWAS of idiopathic polyneuropathy we identify two risk-loci that indicate possible pathogenetic mechanisms. Future collaborative efforts are needed to replicate and expand on these findings.

Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations

Osteoarthritis affects over 300 million people worldwide. Here, we conduct a genome-wide association study meta-analysis across 826,690 individuals (177,517 with osteoarthritis) and identify 100 independently associated risk variants across 11 osteoarthritis phenotypes, 52 of which have not been associated with the disease before. We report thumb and spine osteoarthritis risk variants and identify differences in genetic effects between weight-bearing and non-weight-bearing joints. We identify sex-specific and early age-at-onset osteoarthritis risk loci. We integrate functional genomics data from primary patient tissues (including articular cartilage, subchondral bone, and osteophytic cartilage) and identify high-confidence effector genes. We provide evidence for genetic correlation with phenotypes related to pain, the main disease symptom, and identify likely causal genes linked to neuronal processes. Our results provide insights into key molecular players in disease processes and highlight attractive drug targets to accelerate translation.

A genome-wide association study with 1,126,563 individuals identifies new risk loci for Alzheimer's disease

Late-onset Alzheimer's disease is a prevalent age-related polygenic disease that accounts for 50-70% of dementia cases. Currently, only a fraction of the genetic variants underlying Alzheimer's disease have been identified. Here we show that increased sample sizes allowed identification of seven previously unidentified genetic loci contributing to Alzheimer's disease. This study highlights microglia, immune cells and protein catabolism as relevant to late-onset Alzheimer's disease, while identifying and prioritizing previously unidentified genes of potential interest. We anticipate that these results can be included in larger meta-analyses of Alzheimer's disease to identify further genetic variants that contribute to Alzheimer's pathology.

Prothrombotic genotypes and risk of venous thromboembolism in occult cancer

BACKGROUND

Studies have reported that the combination of some prothrombotic genotypes and overt cancer yields a synergistic effect on VTE risk. Whether individual prothrombotic genotypes or number of risk alleles in a genetic risk score (GRS) affect VTE risk in occult cancer have not been addressed. The aim of this study was to investigate the joint effect of five prothrombotic genotypes and occult cancer on VTE risk.

METHODS

Cases with incident VTE (n = 1566) and a subcohort (n = 14,537) were sampled from the Scandinavian Thrombosis and Cancer Cohort (1993-2012). Five single nucleotide polymorphisms previously reported in a GRS were genotyped: ABO (rs8176719), F5 (rs6025), F2 (rs1799963), FGG (rs2066865) and F11 (rs2036914). Hazard ratios (HRs) for VTE by individual SNPs and GRS were estimated according to non-cancer and occult cancer (one year preceding a cancer diagnosis) exposure.

RESULTS

Occult cancer occurred in 1817 subjects, and of these, 93 experienced a VTE. The VTE risk was 4-fold higher (HR 4.05, 95% CI 3.28-5.00) in subjects with occult cancer compared with those without cancer. Among subjects with occult cancer, those with VTE had a higher proportion of prothrombotic and advanced cancers than those without VTE. The VTE risk increased according to individual prothrombotic genotypes and GRS in cancer-free subjects, while no such effect was observed in subjects with occult cancer (HR for ≥4 versus ≤1 risk alleles in GRS: 1.14, 95% CI 0.61-2.11).

CONCLUSIONS

Five well-established prothrombotic genotypes, individually or combined, were not associated with increased risk of VTE in individuals with occult cancer.

Joint Effect of Multiple Prothrombotic Genotypes and Obesity on the Risk of Incident Venous Thromboembolism

BACKGROUND

 The impact of the combination of obesity and multiple prothrombotic genotypes on venous thromboembolism (VTE) risk remains unclear.

OBJECTIVE

 To investigate the joint effect of obesity and a genetic risk score (GRS) composed of established prothrombotic single nucleotide polymorphisms (SNPs) on VTE risk using a population-based case-cohort.

METHODS

 Cases with incident VTE (n = 1,470) and a subcohort (n = 12,826) were derived from the Tromsø Study (1994-2012) and the Trøndelag Health Study (HUNT) (1995-2008). Participants were genotyped for ABO (rs8176719), F5 (rs6025), F2 (rs1799963), FGG (rs2066865), and F11 (rs2036914) SNPs. Age- and sex-adjusted hazard ratios (HRs) were estimated according to body mass index (BMI) categories and number of risk alleles for individual SNPs and the GRS (0-1, 2, 3, ≥4 alleles).

RESULTS

 The combination of obesity (BMI ≥ 30kg/m²) and risk alleles, either as individual SNPs or as a GRS, had an additive effect on VTE risk (i.e., no biological interaction). Obese subjects who were carriers of ≥4 risk alleles had a 2.85-fold (95% confidence interval [CI]: 2.05-3.96) increased risk of overall VTE compared with those with BMI <25 kg/m² and 0 to 1 risk allele. However, in subgroups, the combination of obesity and ≥4 risk alleles was more pronounced for deep vein thrombosis (DVT) (HR: 3.20; 95% CI: 2.09-4.90) and unprovoked VTE (HR: 3.82; 95% CI: 2.25-6.47), suggesting a supra-additive effect.

CONCLUSION

 Our findings indicate that the combination of obesity and GRS has an additive effect on the risk of overall VTE. However, it may have a supra-additive effect on the risk of DVT and unprovoked VTE.

Loss-of-function genomic variants highlight potential therapeutic targets for cardiovascular disease

Pharmaceutical drugs targeting dyslipidemia and cardiovascular disease (CVD) may increase the risk of fatty liver disease and other metabolic disorders. To identify potential novel CVD drug targets without these adverse effects, we perform genome-wide analyses of participants in the HUNT Study in Norway (n = 69,479) to search for protein-altering variants with beneficial impact on quantitative blood traits related to cardiovascular disease, but without detrimental impact on liver function. We identify 76 (11 previously unreported) presumed causal protein-altering variants associated with one or more CVD- or liver-related blood traits. Nine of the variants are predicted to result in loss-of-function of the protein. This includes ZNF529:p.K405X, which is associated with decreased low-density-lipoprotein (LDL) cholesterol (P = 1.3 × 10-8) without being associated with liver enzymes or non-fasting blood glucose. Silencing of ZNF529 in human hepatoma cells results in upregulation of LDL receptor and increased LDL uptake in the cells. This suggests that inhibition of ZNF529 or its gene product should be prioritized as a novel candidate drug target for treating dyslipidemia and associated CVD.

Combined effects of five prothrombotic genotypes and cancer on the risk of a first venous thromboembolic event

BACKGROUND

The role of combined prothrombotic genotypes in cancer-related venous thromboembolism (VTE) is scarcely studied. We aimed to investigate the impact of a 5-single nucleotide polymorphism (SNP) score on the risk of VTE in patients with and without cancer using a population-based case-cohort.

METHODS

Cases with a first VTE (n = 1493) and a subcohort (n = 13 072) were derived from the Tromsø Study (1994-2012) and the Nord-Trøndelag Health Study (1995-2008). Five SNPs previously reported as a risk score were genotyped: ABO (rs8176719), F5 (rs6025), F2 (rs1799963), FGG (rs2066865), and F11 (rs2036914). Hazard ratios (HRs) for VTE were estimated according to cancer status and the number of risk alleles in the 5-SNP score (0-1, 2-3, and ≥4 alleles).

RESULTS

During a median follow-up of 12.3 years, 1496 individuals were diagnosed with cancer, of whom 232 experienced VTE. The VTE risk increased with the number of risk alleles in the 5-SNP score among subjects without and with cancer. In cancer-free subjects, the HR was 2.17 (95% confidence interval [CI] 1.79-2.62) for ≥4 versus 0-1 risk alleles. In cancer patients, the corresponding HR was 1.93 (95% CI 1.28-2.91). The combination of cancer and ≥4 risk alleles yielded a 17-fold (HR 17.1, 95% CI 12.5-23.4) higher risk of VTE compared with cancer-free subjects with 0-1 risk alleles.

CONCLUSION

The risk of VTE increases with the number of prothrombotic risk alleles in subjects with and without cancer, and the combination of prothrombotic risk alleles and cancer leads to a highly elevated risk of VTE.

MEPE loss-of-function variant associates with decreased bone mineral density and increased fracture risk

A major challenge in genetic association studies is that most associated variants fall in the non-coding part of the human genome. We searched for variants associated with bone mineral density (BMD) after enriching the discovery cohort for loss-of-function (LoF) mutations by sequencing a subset of the Nord-Trøndelag Health Study, followed by imputation in the remaining sample (N = 19,705), and identified ten known BMD loci. However, one previously unreported variant, LoF mutation in MEPE, p.(Lys70IlefsTer26, minor allele frequency [MAF] = 0.8%), was associated with decreased ultradistal forearm BMD (P-value = 2.1 × 10⁻¹⁸), and increased osteoporosis (P-value = 4.2 × 10⁻⁵) and fracture risk (P-value = 1.6 × 10⁻⁵). The MEPE LoF association with BMD and fractures was further evaluated in 279,435 UK (MAF = 0.05%, heel bone estimated BMD P-value = 1.2 × 10⁻¹⁶, any fracture P-value = 0.05) and 375,984 Icelandic samples (MAF = 0.03%, arm BMD P-value = 0.12, forearm fracture P-value = 0.005). Screening for the MEPE LoF mutations before adulthood could potentially prevent osteoporosis and fractures due to the lifelong effect on BMD observed in the study. A key implication for precision medicine is that high-impact functional variants missing from the publicly available cosmopolitan panels could be clinically more relevant than polygenic risk scores.

Bone mineral density (BMD) is associated with fracture risk and many genetic loci with small effect sizes have been discovered by genome-wide association studies (GWAS). Here, the authors discover a large-effect rare loss-of-function genetic variant for BMD in the MEPE gene in the Norwegian HUNT study which replicates in the UK Biobank.

GWAS of thyroid stimulating hormone highlights pleiotropic effects and inverse association with thyroid cancer

Thyroid stimulating hormone (TSH) is critical for normal development and metabolism. To better understand the genetic contribution to TSH levels, we conduct a GWAS meta-analysis at 22.4 million genetic markers in up to 119,715 individuals and identify 74 genome-wide significant loci for TSH, of which 28 are previously unreported. Functional experiments show that the thyroglobulin protein-altering variants P118L and G67S impact thyroglobulin secretion. Phenome-wide association analysis in the UK Biobank demonstrates the pleiotropic effects of TSH-associated variants and a polygenic score for higher TSH levels is associated with a reduced risk of thyroid cancer in the UK Biobank and three other independent studies. Two-sample Mendelian randomization using TSH index variants as instrumental variables suggests a protective effect of higher TSH levels (indicating lower thyroid function) on risk of thyroid cancer and goiter. Our findings highlight the pleiotropic effects of TSH-associated variants on thyroid function and growth of malignant and benign thyroid tumors.

Development and validation of a prediction model for incident hand osteoarthritis in the HUNT study

OBJECTIVE

To develop and externally validate prediction models for incident hand osteoarthritis (OA) in a large population-based cohort of middle aged and older men and women.

DESIGN

We included 17,153 men and 18,682 women from a population-based cohort, aged 35-70 years at baseline (1995-1997). Incident hand OA were obtained from diagnostic codes in the Norwegian National Patient Register (1995-2018). We studied whether a range of self-reported and clinically measured predictors could predict hand OA, using the Area Under the receiver-operating Curve (AUC) from logistic regression. External validation of an existing prediction model for male hand OA was tested on discrimination in a sample of men. Bootstrapping was used to avoid overfitting.

RESULTS

The model for men showed modest discriminatory ability (AUC = 0.67, 95% CI 0.62-0.71). Adding a genetic risk score did not improve prediction. Similar discrimination was observed in the model for women (AUC = 0.62, 95% CI 0.59-0.64). Prediction was not improved by adding a genetic risk score or hormonal and reproductive factors. Applying external validation, similar results were observed among men in HUNT (The Nord-Trøndelag Health Study) as in the developmental sample (AUC = 0.62, 95% CI 0.57-0.65).

CONCLUSION

We developed prediction models for incident hand OA in men and women. For women, the model included body mass index (BMI), heavy physical work, high physical activity and perceived poor health. The model showed moderate discrimination. For men, we have shown that a prediction model including BMI, education and information on sleep can predict incident hand OA in several populations with moderate discriminative ability.

Scalable generalized linear mixed model for region-based association tests in large biobanks and cohorts

With very large sample sizes, biobanks provide an exciting opportunity to identify genetic components of complex traits. To analyze rare variants, region-based multiple-variant aggregate tests are commonly used to increase power for association tests. However, because of the substantial computational cost, existing region-based tests cannot analyze hundreds of thousands of samples while accounting for confounders such as population stratification and sample relatedness. Here we propose a scalable generalized mixed-model region-based association test, SAIGE-GENE, that is applicable to exome-wide and genome-wide region-based analysis for hundreds of thousands of samples and can account for unbalanced case-control ratios for binary traits. Through extensive simulation studies and analysis of the HUNT study with 69,716 Norwegian samples and the UK Biobank data with 408,910 White British samples, we show that SAIGE-GENE can efficiently analyze large-sample data (N > 400,000) with type I error rates well controlled.

Fibrinogen gamma gene rs2066865 and risk of cancer-related venous thromboembolism

Venous thromboembolism (VTE) is a frequent complication in patients with cancer. Homozygous carriers of the fibrinogen gamma gene (FGG) rs2066865 have a moderately increased risk of VTE, but the effect of the FGG variant in cancer is unknown. We aimed to investigate the effect of the FGG variant and active cancer on the risk of VTE. Cases with incident VTE (n=640) and a randomly selected age-weighted sub-cohort (n=3,734) were derived from a population-based cohort (the Tromsø study). Cox-regression was used to estimate hazard ratios (HR) with 95% confidence intervals (CI) for VTE according to categories of cancer and FGG. In those without cancer, homozygosity at the FGG variant was associated with a 70% (HR 1.7, 95% CI: 1.2–2.3) increased risk of VTE compared to non-carriers. Cancer patients homozygous for the FGG variant had a two-fold (HR 2.0, 95% CI: 1.1–3.6) higher risk of VTE than cancer patients without the variant. Moreover, the six-months cumulative incidence of VTE among cancer patients was 6.4% (95% CI: 3.5–11.6) in homozygous carriers of FGG and 3.1% (95% CI: 2.3–4.7) in those without risk alleles. A synergistic effect was observed between rs2066865 and active cancer on the risk of VTE (synergy index: 1.81, 95% CI: 1.02–3.21, attributable proportion: 0.43, 95% CI: 0.11–0.74). In conclusion, homozygosity at the FGG variant and active cancer yielded a synergistic effect on the risk of VTE.

Impact of prothrombotic genotypes on the association between family history of myocardial infarction and venous thromboembolism

BACKGROUND

Family history of myocardial infarction (FHMI) is known to increase the risk of venous thromboembolism (VTE).

OBJECTIVES

To investigate the effect of prothrombotic genotypes on the association between FHMI and VTE in a case-cohort recruited from a general population.

METHODS

Cases with a first VTE (n = 1493) and a subcohort (n = 13 072) were sampled from the Tromsø study (1994-2012) and the Nord-Trøndelag health (HUNT) study (1995-2008). The DNA samples were genotyped for rs8176719 (ABO), rs6025 (F5), rs1799963 (F2), rs2066865 (FGG), and rs2036914 (F11). Participants with missing information on risk alleles (n = 175), FHMI (n = 2769), and BMI (n = 52) were excluded. Cox regression models were used to estimate hazard ratios (HRs) with 95% confidence intervals (CI) for VTE. To explore the role of prothrombotic genotypes for the association between FHMI and VTE, we (a) included the genotypes in the multivariable-adjusted models and (b) assessed the joint effects between FHMI and genotypes on VTE risk.

RESULTS

The FHMI was associated with a 1.3-fold increased risk of VTE (HR 1.32, 95% CI 1.16-1.50) and 1.5-fold increased risk of unprovoked VTE (HR 1.47, 95% CI 1.22-1.78). The risk of VTE by FHMI did not alter after adjustment for the five genotypes. The combination of FHMI and the different prothrombotic genotypes did not result in an excess VTE risk (i.e. no biological interaction).

CONCLUSIONS

Our findings suggest that the risk of VTE by FHMI is not explained by rs8176719 (ABO), rs6025 (F5), rs1799963 (F2), rs2066865 (FGG), and rs2036914 (F11). The combination of FHMI with prothrombotic genotypes had an additive effect on VTE risk.

Is smoking heaviness causally associated with alcohol use? A Mendelian randomization study in four European cohorts

Background

Observational studies have shown that tobacco and alcohol use co-occur, but it is not clear whether this relationship is causal.

Methods

Using data from the Avon Longitudinal Study of Parents and Children (ALSPAC) and UK Biobank, we used observational methods to test the hypothesis that smoking heaviness increases alcohol consumption. Mendelian randomization (MR) analyses were then used to test the causal relationship between smoking heaviness and alcohol consumption using 55 967 smokers from four European studies [ALSPAC, The Nord-Trøndelag Health Study (HUNT), the Copenhagen General Population Study (CGPS) and UK Biobank]. MR analyses used rs1051730/rs16969968 as a genetic proxy for smoking heaviness.

Results

Observational results provided evidence of an association between cigarettes per day and weekly alcohol consumption (increase in units of alcohol per additional cigarette smoked per day = 0.10, 95% confidence interval (CI) 0.05 to 0.15, P ≤ 0.001 in ALSPAC; and 0.48, 95% CI 0.45 to 0.52, P ≤ 0.001 in UK Biobank). However, there was little evidence for an association between rs1051730/rs16969968 and units of alcohol consumed per week across ALSPAC, HUNT, CGPS and UK Biobank (standard deviation increase in units of alcohol per additional copy of the risk allele = –0.004, 95% CI –0.023 to 0.016, P=0.708, I² = 51.9%). We had 99% and 88% power to detect a change of 0.03 and 0.02 standard deviation units of alcohol per additional copy of the risk allele, respectively.

Conclusions

Previously reported associations between smoking and alcohol are unlikely to be causal, and may be the result of confounding and/or reverse causation. This has implications for public health research and intervention research.

Sex-specific and pleiotropic effects underlying kidney function identified from GWAS meta-analysis

Chronic kidney disease (CKD) is a growing health burden currently affecting 10–15% of adults worldwide. Estimated glomerular filtration rate (eGFR) as a marker of kidney function is commonly used to diagnose CKD. We analyze eGFR data from the Nord-Trøndelag Health Study and Michigan Genomics Initiative and perform a GWAS meta-analysis with public summary statistics, more than doubling the sample size of previous meta-analyses. We identify 147 loci (53 novel) associated with eGFR, including genes involved in transcriptional regulation, kidney development, cellular signaling, metabolism, and solute transport. Additionally, sex-stratified analysis identifies one locus with more significant effects in women than men. Using genetic risk scores constructed from these eGFR meta-analysis results, we show that associated variants are generally predictive of CKD with only modest improvements in detection compared with other known clinical risk factors. Collectively, these results yield additional insight into the genetic factors underlying kidney function and progression to CKD.

Estimated glomerular filtration rate (eGFR) is a measure of kidney function and used to characterize chronic kidney disease. Here, Graham et al. identify 53 novel loci for eGFR in a GWAS meta-analysis, a subset of which are associated with other common diseases, such as diabetes and hypertension, based on PheWAS.

Genome-wide Study of Atrial Fibrillation Identifies Seven Risk Loci and Highlights Biological Pathways and Regulatory Elements Involved in Cardiac Development

Atrial fibrillation (AF) is a common cardiac arrhythmia and a major risk factor for stroke, heart failure, and premature death. The pathogenesis of AF remains poorly understood, which contributes to the current lack of highly effective treatments. To understand the genetic variation and biology underlying AF, we undertook a genome-wide association study (GWAS) of 6,337 AF individuals and 61,607 AF-free individuals from Norway, including replication in an additional 30,679 AF individuals and 278,895 AF-free individuals. Through genotyping and dense imputation mapping from whole-genome sequencing, we tested almost nine million genetic variants across the genome and identified seven risk loci, including two novel loci. One novel locus (lead single-nucleotide variant [SNV] rs12614435; p = 6.76 × 10^-18) comprised intronic and several highly correlated missense variants situated in the I-, A-, and M-bands of titin, which is the largest protein in humans and responsible for the passive elasticity of heart and skeletal muscle. The other novel locus (lead SNV rs56202902; p = 1.54 × 10^-11) covered a large, gene-dense chromosome 1 region that has previously been linked to cardiac conduction. Pathway and functional enrichment analyses suggested that many AF-associated genetic variants act through a mechanism of impaired muscle cell differentiation and tissue formation during fetal heart development.

Heavier smoking increases coffee consumption: findings from a Mendelian randomization analysis

Background

There is evidence for a positive relationship between cigarette and coffee consumption in smokers. Cigarette smoke increases metabolism of caffeine, so this may represent a causal effect of smoking on caffeine intake.

Methods

We performed Mendelian randomization analyses in the UK Biobank (N = 114 029), the Norwegian HUNT study (N = 56 664) and the Copenhagen General Population Study (CGPS) (N = 78 650). We used the rs16969968 genetic variant as a proxy for smoking heaviness in all studies and rs4410790 and rs2472297 as proxies for coffee consumption in UK Biobank and CGPS. Analyses were conducted using linear regression and meta-analysed across studies.

Results

Each additional cigarette per day consumed by current smokers was associated with higher coffee consumption (0.10 cups per day, 95% CI: 0.03, 0.17). There was weak evidence for an increase in tea consumption per additional cigarette smoked per day (0.04 cups per day, 95% CI: -0.002, 0.07). There was strong evidence that each additional copy of the minor allele of rs16969968 (which increases daily cigarette consumption) in current smokers was associated with higher coffee consumption (0.16 cups per day, 95% CI: 0.11, 0.20), but only weak evidence for an association with tea consumption (0.04 cups per day, 95% CI: -0.01, 0.09). There was no clear evidence that rs16969968 was associated with coffee or tea consumption in never or former smokers or that the coffee-related variants were associated with cigarette consumption.

Conclusions

Higher cigarette consumption causally increases coffee intake. This is consistent with faster metabolism of caffeine by smokers, but could also reflect a behavioural effect of smoking on coffee drinking.

The causal role of smoking on the risk of hip or knee replacement due to primary osteoarthritis: a Mendelian randomisation analysis of the HUNT study

OBJECTIVE

Smoking has been associated with a reduced risk of hip and knee osteoarthritis (OA) and subsequent joint replacement. The aim of the present study was to assess whether the observed association is likely to be causal.

METHOD

55,745 participants of a population-based cohort were genotyped for the rs1051730 C > T single-nucleotide polymorphism (SNP), a proxy for smoking quantity among smokers. A Mendelian randomization analysis was performed using rs1051730 as an instrument to evaluate the causal role of smoking on the risk of hip or knee replacement (combined as total joint replacement (TJR)). Association between rs1051730 T alleles and TJR was estimated by hazard ratios (HRs) and 95% confidence intervals (CIs). All analyses were adjusted for age and sex.

RESULTS

Smoking quantity (no. of cigarettes) was inversely associated with TJR (HR 0.97, 95% CI 0.97-0.98). In the Mendelian randomization analysis, rs1051730 T alleles were associated with reduced risk of TJR among current smokers (HR 0.84, 95% CI 0.76-0.98, per T allele), however we found no evidence of association among former (HR 0.97, 95% CI 0.88-1.07) and never smokers (HR 0.97, 95% CI 0.89-1.06). Neither adjusting for body mass index (BMI), cardiovascular disease (CVD) nor accounting for the competing risk of mortality substantially changed the results.

CONCLUSION

This study suggests that smoking may be causally associated with the reduced risk of TJR. Our findings add support to the inverse association found in previous observational studies. More research is needed to further elucidate the underlying mechanisms of this causal association.

Investigating the causal effect of smoking on hay fever and asthma: a Mendelian randomization meta-analysis in the CARTA consortium

Observational studies on smoking and risk of hay fever and asthma have shown inconsistent results. However, observational studies may be biased by confounding and reverse causation. Mendelian randomization uses genetic variants as markers of exposures to examine causal effects. We examined the causal effect of smoking on hay fever and asthma by using the smoking-associated single nucleotide polymorphism (SNP) rs16969968/rs1051730. We included 231,020 participants from 22 population-based studies. Observational analyses showed that current vs never smokers had lower risk of hay fever (odds ratio (OR) = 0·68, 95% confidence interval (CI): 0·61, 0·76; P < 0·001) and allergic sensitization (OR = 0·74, 95% CI: 0·64, 0·86; P < 0·001), but similar asthma risk (OR = 1·00, 95% CI: 0·91, 1·09; P = 0·967). Mendelian randomization analyses in current smokers showed a slightly lower risk of hay fever (OR = 0·958, 95% CI: 0·920, 0·998; P = 0·041), a lower risk of allergic sensitization (OR = 0·92, 95% CI: 0·84, 1·02; P = 0·117), but higher risk of asthma (OR = 1·06, 95% CI: 1·01, 1·11; P = 0·020) per smoking-increasing allele. Our results suggest that smoking may be causally related to a higher risk of asthma and a slightly lower risk of hay fever. However, the adverse events associated with smoking limit its clinical significance.

Effect of Smoking on Blood Pressure and Resting Heart Rate: A Mendelian Randomization Meta-Analysis in the CARTA Consortium

BACKGROUND

Smoking is an important cardiovascular disease risk factor, but the mechanisms linking smoking to blood pressure are poorly understood.

METHODS AND RESULTS

Data on 141 317 participants (62 666 never, 40 669 former, 37 982 current smokers) from 23 population-based studies were included in observational and Mendelian randomization meta-analyses of the associations of smoking status and smoking heaviness with systolic and diastolic blood pressure, hypertension, and resting heart rate. For the Mendelian randomization analyses, a genetic variant rs16969968/rs1051730 was used as a proxy for smoking heaviness in current smokers. In observational analyses, current as compared with never smoking was associated with lower systolic blood pressure and diastolic blood pressure and lower hypertension risk, but with higher resting heart rate. In observational analyses among current smokers, 1 cigarette/day higher level of smoking heaviness was associated with higher (0.21 bpm; 95% confidence interval 0.19; 0.24) resting heart rate and slightly higher diastolic blood pressure (0.05 mm Hg; 95% confidence interval 0.02; 0.08) and systolic blood pressure (0.08 mm Hg; 95% confidence interval 0.03; 0.13). However, in Mendelian randomization analyses among current smokers, although each smoking increasing allele of rs16969968/rs1051730 was associated with higher resting heart rate (0.36 bpm/allele; 95% confidence interval 0.18; 0.54), there was no strong association with diastolic blood pressure, systolic blood pressure, or hypertension. This would suggest a 7 bpm higher heart rate in those who smoke 20 cigarettes/day.

CONCLUSIONS

This Mendelian randomization meta-analysis supports a causal association of smoking heaviness with higher level of resting heart rate, but not with blood pressure. These findings suggest that part of the cardiovascular risk of smoking may operate through increasing resting heart rate.

Identification of lung cancer histology-specific variants applying Bayesian framework variant prioritization approaches within the TRICL and ILCCO consortia

Large-scale genome-wide association studies (GWAS) have likely uncovered all common variants at the GWAS significance level. Additional variants within the suggestive range (0.0001> P > 5×10(-8)) are, however, still of interest for identifying causal associations. This analysis aimed to apply novel variant prioritization approaches to identify additional lung cancer variants that may not reach the GWAS level. Effects were combined across studies with a total of 33456 controls and 6756 adenocarcinoma (AC; 13 studies), 5061 squamous cell carcinoma (SCC; 12 studies) and 2216 small cell lung cancer cases (9 studies). Based on prior information such as variant physical properties and functional significance, we applied stratified false discovery rates, hierarchical modeling and Bayesian false discovery probabilities for variant prioritization. We conducted a fine mapping analysis as validation of our methods by examining top-ranking novel variants in six independent populations with a total of 3128 cases and 2966 controls. Three novel loci in the suggestive range were identified based on our Bayesian framework analyses: KCNIP4 at 4p15.2 (rs6448050, P = 4.6×10(-7)) and MTMR2 at 11q21 (rs10501831, P = 3.1×10(-6)) with SCC, as well as GAREM at 18q12.1 (rs11662168, P = 3.4×10(-7)) with AC. Use of our prioritization methods validated two of the top three loci associated with SCC (P = 1.05×10(-4) for KCNIP4, represented by rs9799795) and AC (P = 2.16×10(-4) for GAREM, represented by rs3786309) in the independent fine mapping populations. This study highlights the utility of using prior functional data for sequence variants in prioritization analyses to search for robust signals in the suggestive range.

Stratification by smoking status reveals an association of CHRNA5-A3-B4 genotype with body mass index in never smokers

We previously used a single nucleotide polymorphism (SNP) in the CHRNA5-A3-B4 gene cluster associated with heaviness of smoking within smokers to confirm the causal effect of smoking in reducing body mass index (BMI) in a Mendelian randomisation analysis. While seeking to extend these findings in a larger sample we found that this SNP is associated with 0.74% lower body mass index (BMI) per minor allele in current smokers (95% CI -0.97 to -0.51, P = 2.00×10⁻¹⁰), but also unexpectedly found that it was associated with 0.35% higher BMI in never smokers (95% CI +0.18 to +0.52, P = 6.38×10⁻⁵). An interaction test confirmed that these estimates differed from each other (P = 4.95×10⁻¹³). This difference in effects suggests the variant influences BMI both via pathways unrelated to smoking, and via the weight-reducing effects of smoking. It would therefore be essentially undetectable in an unstratified genome-wide association study of BMI, given the opposite association with BMI in never and current smokers. This demonstrates that novel associations may be obscured by hidden population sub-structure. Stratification on well-characterized environmental factors known to impact on health outcomes may therefore reveal novel genetic associations.

We found that a single nucleotide polymorphism in the CHRNA5-A3-B4 gene cluster, which is known to influence smoking heaviness, is associated with lower body mass index (BMI) in current smokers, but higher BMI in never smokers. This difference in effects suggests that the variant influences BMI both via pathways other than smoking, and via the weight-reducing effects of smoking, in opposite directions. The overall effect on BMI would therefore be undetectable in an unstratified genome-wide association study, indicating that novel associations may be obscured by hidden population sub-structure.

Causal associations of tobacco smoking with cardiovascular risk factors: a Mendelian randomization analysis of the HUNT Study in Norway

BACKGROUND

Tobacco smoking has been associated with cardiovascular risk factors including adverse serum lipid levels, central obesity and higher resting heart rate, but lower blood pressure and body mass index (BMI). We used a Mendelian randomization approach to study whether these associations may be causal. If smoking affects cardiovascular risk factors then rs1051730 T alleles, predictors of increased smoking quantity, should be associated with cardiovascular risk factors among smokers, but not among never smokers.

METHODS

Among 56,625 participants of a population-based study, we estimated associations of rs1051730 T alleles with cardiovascular risk factors and examined whether the associations differed by smoking status.

RESULTS

Rs1051730 T alleles were associated with lower BMI and waist and hip circumferences and higher resting heart rate and estimated glomerular filtration rate (eGFR), and the associations were strongest among current smokers (P interaction 5×10(-9) to 0.01). Rs1051730 T alleles were associated with lower systolic blood pressure and pulse pressure and higher HDL cholesterol concentrations, but these associations did not robustly differ by smoking status. There were no convincing associations of rs1051730 T alleles with waist-hip ratio, diastolic blood pressure and non-fasting serum concentrations of non-HDL cholesterol, triglycerides, glucose and C-reactive protein.

CONCLUSIONS

This Mendelian randomization analysis provides evidence that smoking may cause lower BMI and waist and hip circumferences and higher resting heart rate and eGFR. The findings further suggest that smoking is not a major determinant of waist-hip ratio or adverse blood pressure, serum lipid or glucose levels.

Association between a 15q25 gene variant, nicotine-related habits, lung cancer and COPD among 56,307 individuals from the HUNT study in Norway

Genetic studies have shown an association between single-nucleotide polymorphisms on chromosome 15q25 and smoking-related traits and diseases, such as quantity of smoking, lung cancer and chronic obstructive pulmonary disease (COPD). A discussion has centred on the variants and their effects being directly disease related or indirect via nicotine addiction. To address these discrepancies, we genotyped the single-nucleotide polymorphism rs16969968 in the CHRNA5/A3/B4 gene cluster at chromosome 15q25, in 56 307 individuals from a large homogenous population-based cohort, the North Trøndelag Health Study (HUNT) in Norway. The variant was examined in relation to four different outcomes: lung cancer, loss of lung function equivalent to that of COPD, smoking behaviour and the use of smokeless tobacco (snus). Novel associations were found between rs16969968 and the motivational factor for starting to use snus, and the quantity of snus used. Our results also confirm and extend previous findings for associations between rs16969968 and lung cancer, loss of lung function equivalent to that of COPD, and smoking quantity. Our data suggest a role for rs16969968 in nicotine addiction, and the novel association with snus strengthens this observation.