Metabolomic pattern associated with physical sequelae in patients presenting with respiratory symptoms validates the aestivation concept in dehydrated patients.

BACKGROUND: Hypertonic dehydration is associated with muscle wasting and synthesis of organic osmolytes. We recently showed a metabolic shift to amino acid production and urea cycle activation in COVID-19, consistent with the aestivation response. The aim of the present investigation was to validate the metabolic shift and development of long-term physical outcome in the non-COVID cohort of the Biobanque Québécoise de la COVID-19 (BQC19). METHODS: We included 824 patients from BQC19, where of 571 patients had data of dehydration in the form of estimated osmolality (eOSM = 2Na+2K+glucose+urea), and 284 patients had metabolome data and long-term follow-up. We correlated the degree of dehydration to mortality, invasive mechanical ventilation, acute kidney injury, and long-term symptoms. RESULTS: As found in the COVID cohort, higher eOSM correlated with higher proportion of urea and glucose of total eOSM and an enrichment of amino acids compared to other metabolites. Sex stratified analysis indicated that women may show a weaker aestivation response. More severe dehydration was associated with mortality, invasive mechanical ventilation, and acute kidney injury during the acute illness. Importantly, more severe dehydration was associated with physical long-term symptoms but not mental long-term symptoms after adjustment for age, sex, and disease severity. CONCLUSIONS: Patients with water deficit in the form of increased eOSM tend to have more severe disease and experience more physical symptoms after an acute episode of care. This is associated with amino acid and urea production indicating dehydration induced muscle wasting.

Circulating metabolite abundances associated with risks of bipolar disorder, schizophrenia, and depression: a Mendelian randomization study.

BACKGROUND: Preventive measures and treatments for psychiatric disorders are limited. Circulating metabolites are potential candidates for biomarker and therapeutic target identification, given their measurability and essential roles in biological processes. METHODS: Leveraging large-scale genome-wide association studies, we conducted Mendelian randomization (MR) analyses to assess the associations between circulating metabolite abundances and the risks of bipolar disorder, schizophrenia, and depression. Genetic instruments were selected for 94 metabolites measured in the Canadian Longitudinal Study of Aging (N=8,299) cohort. We repeated MR analyses based on the UK Biobank, INTERVAL, and EPIC-Norfolk studies. RESULTS: After validating MR assumptions and colocalization evidence, we found that a one standard deviation (SD) increase in genetically predicted circulating abundances of eicosapentaenoate (EPA) and docosapentaenoate (n3 DPA) was associated with odds ratios (ORs) of 0.72 (95% CI: 0.65-0.79) and 0.63 (95% CI: 0.55-0.72) for bipolar disorder, respectively. Genetically increased Ω-3 unsaturated fatty acids abundance and Ω-3-to-total fatty acids ratio, as well as genetically decreased Ω-6-to-Ω-3 ratio were negatively associated with the risk of bipolar disorder in the UK Biobank. Genetically increased circulating abundances of three N-acetyl-amino acids were associated with an increased risk of schizophrenia with a maximum OR of 1.31 (95% CI: 1.18-1.44) per one SD increase. Furthermore, a one SD increase in genetically predicted circulating abundance of hypotaurine was associated with an OR of 0.85 (95% CI: 0.78-0.93) for depression. CONCLUSIONS: The biological mechanisms underlying Ω-3 unsaturated fatty acids, NAT8-catalyzed N-acetyl-amino acids, and hypotaurine warrant exploration to identify new biomarkers and potential therapeutic targets.

HLA allele-calling using multi-ancestry whole-exome sequencing from the UK Biobank identifies 129 novel associations in 11 autoimmune diseases.

The human leukocyte antigen (HLA) region on chromosome 6 is strongly associated with many immune-mediated and infection-related diseases. Due to its highly polymorphic nature and complex linkage disequilibrium patterns, traditional genetic association studies of single nucleotide polymorphisms do not perform well in this region. Instead, the field has adopted the assessment of the association of HLA alleles (i.e., entire HLA gene haplotypes) with disease. Often based on genotyping arrays, these association studies impute HLA alleles, decreasing accuracy and thus statistical power for rare alleles and in non-European ancestries. Here, we use whole-exome sequencing (WES) from 454,824 UK Biobank (UKB) participants to directly call HLA alleles using the HLA-HD algorithm. We show this method is more accurate than imputing HLA alleles and harness the improved statistical power to identify 360 associations for 11 auto-immune phenotypes (at least 129 likely novel), leading to better insights into the specific coding polymorphisms that underlie these diseases. We show that HLA alleles with synonymous variants, often overlooked in HLA studies, can significantly influence these phenotypes. Lastly, we show that HLA sequencing may improve polygenic risk scores accuracy across ancestries. These findings allow better characterization of the role of the HLA region in human disease.

Identifying Rare Genetic Determinants for Improved Polygenic Risk Prediction of Bone Mineral Density and Fracture Risk.

Osteoporosis and fractures severely impact the elderly population. Polygenic risk scores for bone mineral density have demonstrated potential clinical utility. However, the value of rare genetic determinants in risk prediction has not been assessed. With whole-exome sequencing data from 436,824 UK Biobank participants, we assigned White British ancestry individuals into a training data set (n = 317,434) and a test data set (n = 74,825). In the training data set, we developed a common variant-based polygenic risk score for heel ultrasound speed of sound (SOS). Next, we performed burden testing to identify genes harboring rare determinants of bone mineral density, targeting influential rare variants with predicted high deleteriousness. We constructed a genetic risk score, called ggSOS, to incorporate influential rare variants in significant gene burden masks into the common variant-based polygenic risk score. We assessed the predictive performance of ggSOS in the White British test data set, as well as in populations of non-White British European (n = 18,885), African (n = 7165), East Asian (n = 2236), South Asian (n = 9829), and other admixed (n = 1481) ancestries. Twelve genes in pivotal regulatory pathways of bone homeostasis harbored influential rare variants associated with SOS (p < 5.5 × 10-7 ), including AHNAK, BMP5, CYP19A1, FAM20A, FBXW5, KDM5B, KREMEN1, LGR4, LRP5, SMAD6, SOST, and WNT1. Among 4013 (5.4%) individuals in the test data set carrying these variants, a one standard deviation decrease in ggSOS was associated with 1.35-fold (95% confidence interval [CI] 1.16-1.57) increased hazard of major osteoporotic fracture. However, compared with a common variant-based polygenic risk score (C-index = 0.641), ggSOS had only marginally improved prediction accuracy in identifying at-risk individuals (C-index = 0.644), with overlapping confidence intervals. Similarly, ggSOS did not demonstrate substantially improved predictive performance in non-European ancestry populations. In summary, modeling the effects of rare genetic determinants may assist polygenic prediction of fracture risk among carriers of influential rare variants. Nonetheless, improved clinical utility is not guaranteed for population-level risk screening. © 2023 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

Converging evidence from exome sequencing and common variants implicates target genes for osteoporosis.

In this study, we leveraged the combined evidence of rare coding variants and common alleles to identify therapeutic targets for osteoporosis. We undertook a large-scale multiancestry exome-wide association study for estimated bone mineral density, which showed that the burden of rare coding alleles in 19 genes was associated with estimated bone mineral density (P < 3.6 × 10-7). These genes were highly enriched for a set of known causal genes for osteoporosis (65-fold; P = 2.5 × 10-5). Exome-wide significant genes had 96-fold increased odds of being the top ranked effector gene at a given GWAS locus (P = 1.8 × 10-10). By integrating proteomics Mendelian randomization evidence, we prioritized CD109 (cluster of differentiation 109) as a gene for which heterozygous loss of function is associated with higher bone density. CRISPR-Cas9 editing of CD109 in SaOS-2 osteoblast-like cell lines showed that partial CD109 knockdown led to increased mineralization. This study demonstrates that the convergence of common and rare variants, proteomics and CRISPR can highlight new bone biology to guide therapeutic development.

Predicting ExWAS findings from GWAS data: a shorter path to causal genes.

GWAS has identified thousands of loci associated with disease, yet the causal genes within these loci remain largely unknown. Identifying these causal genes would enable deeper understanding of the disease and assist in genetics-based drug development. Exome-wide association studies (ExWAS) are more expensive but can pinpoint causal genes offering high-yield drug targets, yet suffer from a high false-negative rate. Several algorithms have been developed to prioritize genes at GWAS loci, such as the Effector Index (Ei), Locus-2-Gene (L2G), Polygenic Prioritization score (PoPs), and Activity-by-Contact score (ABC) and it is not known if these algorithms can predict ExWAS findings from GWAS data. However, if this were the case, thousands of associated GWAS loci could potentially be resolved to causal genes. Here, we quantified the performance of these algorithms by evaluating their ability to identify ExWAS significant genes for nine traits. We found that Ei, L2G, and PoPs can identify ExWAS significant genes with high areas under the precision recall curve (Ei: 0.52, L2G: 0.37, PoPs: 0.18, ABC: 0.14). Furthermore, we found that for every unit increase in the normalized scores, there was an associated 1.3-4.6-fold increase in the odds of a gene reaching exome-wide significance (Ei: 4.6, L2G: 2.5, PoPs: 2.1, ABC: 1.3). Overall, we found that Ei, L2G, and PoPs can anticipate ExWAS findings from widely available GWAS results. These techniques are therefore promising when well-powered ExWAS data are not readily available and can be used to anticipate ExWAS findings, allowing for prioritization of genes at GWAS loci.

Increasing serum iron levels and their role in the risk of infectious diseases: a Mendelian randomization approach.

OBJECTIVES: Increased iron stores have been associated with elevated risks of different infectious diseases, suggesting that iron supplementation may increase the risk of infections. However, these associations may be biased by confounding or reverse causation. This is important, since up to 19% of the population takes iron supplementation. We used Mendelian randomization (MR) to bypass these biases and estimate the causal effect of iron on infections. METHODS: As instrumental variables, we used genetic variants associated with iron biomarkers in two genome-wide association studies (GWASs) of European ancestry participants. For outcomes, we used GWAS results from the UK Biobank, FinnGen, the COVID-19 Host Genetics Initiative or 23andMe, for seven infection phenotypes: 'any infections', combined, COVID-19 hospitalization, candidiasis, pneumonia, sepsis, skin and soft tissue infection (SSTI) and urinary tract infection (UTI). RESULTS: Most of our analyses showed increasing iron (measured by its biomarkers) was associated with only modest changes in the odds of infectious outcomes, with all 95% odds ratios confidence intervals within the 0.88 to 1.26 range. However, for the three predominantly bacterial infections (sepsis, SSTI, UTI), at least one analysis showed a nominally elevated risk with increased iron stores (P <0.05). CONCLUSION: Using MR, we did not observe an increase in risk of most infectious diseases with increases in iron stores. However for bacterial infections, higher iron stores may increase odds of infections. Hence, using genetic variation in iron pathways as a proxy for iron supplementation, iron supplements are likely safe on a population level, but we should continue the current practice of conservative iron supplementation during bacterial infections or in those at high risk of developing them.

Proteome-wide Mendelian randomization implicates nephronectin as an actionable mediator of the effect of obesity on COVID-19 severity.

Obesity is a major risk factor for Coronavirus disease (COVID-19) severity; however, the mechanisms underlying this relationship are not fully understood. As obesity influences the plasma proteome, we sought to identify circulating proteins mediating the effects of obesity on COVID-19 severity in humans. Here, we screened 4,907 plasma proteins to identify proteins influenced by body mass index using Mendelian randomization. This yielded 1,216 proteins, whose effect on COVID-19 severity was assessed, again using Mendelian randomization. We found that an s.d. increase in nephronectin (NPNT) was associated with increased odds of critically ill COVID-19 (OR = 1.71, P = 1.63 × 10-10). The effect was driven by an NPNT splice isoform. Mediation analyses supported NPNT as a mediator. In single-cell RNA-sequencing, NPNT was expressed in alveolar cells and fibroblasts of the lung in individuals who died of COVID-19. Finally, decreasing body fat mass and increasing fat-free mass were found to lower NPNT levels. These findings provide actionable insights into how obesity influences COVID-19 severity.

Genomic atlas of the plasma metabolome prioritizes metabolites implicated in human diseases.

Metabolic processes can influence disease risk and provide therapeutic targets. By conducting genome-wide association studies of 1,091 blood metabolites and 309 metabolite ratios, we identified associations with 690 metabolites at 248 loci and associations with 143 metabolite ratios at 69 loci. Integrating metabolite-gene and gene expression information identified 94 effector genes for 109 metabolites and 48 metabolite ratios. Using Mendelian randomization (MR), we identified 22 metabolites and 20 metabolite ratios having estimated causal effect on 12 traits and diseases, including orotate for estimated bone mineral density, α-hydroxyisovalerate for body mass index and ergothioneine for inflammatory bowel disease and asthma. We further measured the orotate level in a separate cohort and demonstrated that, consistent with MR, orotate levels were positively associated with incident hip fractures. This study provides a valuable resource describing the genetic architecture of metabolites and delivers insights into their roles in common diseases, thereby offering opportunities for therapeutic targets.

Circulating Proteins Influencing Psychiatric Disease: A Mendelian Randomization Study.

BACKGROUND: There is a pressing need for novel drug targets for psychiatric disorders. Circulating proteins are potential candidates because they are relatively easy to measure and modulate and play important roles in signaling. METHODS: We performed two-sample Mendelian randomization analyses to estimate the associations between circulating protein abundances and risk of 10 psychiatric disorders. Genetic variants associated with 1611 circulating protein abundances identified in 6 large-scale proteomic studies were used as genetic instruments. Effects of the circulating proteins on psychiatric disorders were estimated by Wald ratio or inverse variance-weighted ratio tests. Horizontal pleiotropy, colocalization, and protein-altering effects were examined to validate the assumptions of Mendelian randomization. RESULTS: Nine circulating protein-to-disease associations withstood multiple sensitivity analyses. Among them, 2 circulating proteins had associations replicated in 3 proteomic studies. A 1 standard deviation increase in the genetically predicted circulating TIMP4 level was associated with a reduced risk of anorexia nervosa (minimum odds ratio [OR] = 0.83; 95% CI, 0.76-0.91) and bipolar disorder (minimum OR = 0.88; 95% CI, 0.82-0.94). A 1 standard deviation increase in the genetically predicted circulating ESAM level was associated with an increased risk of schizophrenia (maximum OR = 1.32; 95% CI, 1.22-1.43). In addition, 58 suggestive protein-to-disease associations warrant validation with observational or experimental evidence. For instance, a 1 standard deviation increase in the ERLEC1-201-to-ERLEC1-202 splice variant ratio was associated with a reduced risk of schizophrenia (OR = 0.94; 95% CI, 0.90-0.97). CONCLUSIONS: Prioritized circulating proteins appear to influence the risk of psychiatric disease and may be explored as intervention targets.

Dehydration is associated with production of organic osmolytes and predicts physical long-term symptoms after COVID-19: a multicenter cohort study.

BACKGROUND: We have previously shown that iatrogenic dehydration is associated with a shift to organic osmolyte production in the general ICU population. The aim of the present investigation was to determine the validity of the physiological response to dehydration known as aestivation and its relevance for long-term disease outcome in COVID-19. METHODS: The study includes 374 COVID-19 patients from the Pronmed cohort admitted to the ICU at Uppsala University Hospital. Dehydration data was available for 165 of these patients and used for the primary analysis. Validation was performed in Biobanque Québécoise de la COVID-19 (BQC19) using 1052 patients with dehydration data. Dehydration was assessed through estimated osmolality (eOSM = 2Na + 2 K + glucose + urea), and correlated to important endpoints including death, invasive mechanical ventilation, acute kidney injury, and long COVID-19 symptom score grouped by physical or mental. RESULTS: Increasing eOSM was correlated with increasing role of organic osmolytes for eOSM, while the proportion of sodium and potassium of eOSM were inversely correlated to eOSM. Acute outcomes were associated with pronounced dehydration, and physical long-COVID was more strongly associated with dehydration than mental long-COVID after adjustment for age, sex, and disease severity. Metabolomic analysis showed enrichment of amino acids among metabolites that showed an aestivating pattern. CONCLUSIONS: Dehydration during acute COVID-19 infection causes an aestivation response that is associated with protein degradation and physical long-COVID. TRIAL REGISTRATION: The study was registered à priori (clinicaltrials.gov: NCT04316884 registered on 2020-03-13 and NCT04474249 registered on 2020-06-29).

Genetic determinants of polygenic prediction accuracy within a population.

Genomic risk prediction is on the emerging path toward personalized medicine. However, the accuracy of polygenic prediction varies strongly in different individuals. Based on up to 352,277 European ancestry participants in the UK Biobank, we constructed polygenic risk scores for 15 physiological and biochemical quantitative traits. We identified a total of 185 polygenic prediction variability quantitative trait loci for 11 traits by Levene's test among 254,376 unrelated individuals. We validated the effects of prediction variability quantitative trait loci using an independent test set of 58,927 individuals. For instance, a score aggregating 51 prediction variability quantitative trait locus variants for triglycerides had the strongest Spearman correlation of 0.185 (P-value <1.0 × 10-300) with the squared prediction errors. We found a strong enrichment of complex genetic effects conferred by prediction variability quantitative trait loci compared to risk loci identified in genome-wide association studies, including 89 prediction variability quantitative trait loci exhibiting dominance effects. Incorporation of dominance effects into polygenic risk scores significantly improved polygenic prediction for triglycerides, low-density lipoprotein cholesterol, vitamin D, and platelet. In conclusion, we have discovered and profiled genetic determinants of polygenic prediction variability for 11 quantitative biomarkers. These findings may assist interpretation of genomic risk prediction in various contexts and encourage novel approaches for constructing polygenic risk scores with complex genetic effects.

Capturing additional genetic risk from family history for improved polygenic risk prediction.

Family history of complex traits may reflect transmitted rare pathogenic variants, intra-familial shared exposures to environmental and lifestyle factors, as well as a common genetic predisposition. We developed a latent factor model to quantify trait heritability in excess of that captured by a common variant-based polygenic risk score, but inferable from family history. For 941 children in the Avon Longitudinal Study of Parents and Children cohort, a joint predictor combining a polygenic risk score for height and mid-parental height was able to explain ~55% of the total variance in sex-adjusted adult height z-scores, close to the estimated heritability. Marginal yet consistent risk prediction improvements were also achieved among ~400,000 European ancestry participants for 11 complex diseases in the UK Biobank. Our work showcases a paradigm for risk calculation, and supports incorporation of family history into polygenic risk score-based genetic risk prediction models.

Identifying Causes of Fracture Beyond Bone Mineral Density: Evidence From Human Genetics.

New therapies may help to prevent osteoporotic fractures other than through increasing bone mineral density (BMD). Because fracture risk has an important genetic component, we aim to identify loci increasing fracture risk that do not decrease BMD, using a recently-proposed structural equation model adapted to remove genetic influences of BMD on fracture risk. We used summary statistics of the largest genome-wide association studies (GWASs) for BMD and for fracture in these analyses. We next estimated the genetic correlation between the non-BMD or BMD-related genetic effects and other clinical risk factors for fracture. Last, based on white British participants in the UK Biobank, we conducted genetic risk score analyses to assess whether the aggregated genetic effects conferred increased major osteoporotic fracture risk. We found that only three loci affecting fracture risk exhibited genetic effects not mediated by BMD: SOST, CPED1-WNT16, and RSPO3, while these three loci simultaneously conferred BMD-related effects. No strong genetic associations between non-BMD or BMD-related effects and 16 clinical risk factors were observed. However, non-BMD effects might be genetic correlated with hip bone size. In the UK Biobank, a 1 standard deviation (1-SD) increase in the non-BMD genetic risk score conferred an odds ratio of 1.17 for incident major osteoporotic fracture, compared to 1.29 by a BMD-related genetic risk score. Our study suggests that the majority of common genetic predisposition toward fracture risk acts upon BMD. Although non-BMD genetic effects may exist, they are not strongly correlated with most traditional clinical risk factors. Risk loci harboring non-BMD genetic effects may influence other perspectives of bone quality, or confer effects that existing GWASs fail to capture, but they demonstrate weaker impact on fracture risk than BMD-related genetic effects. These findings suggest that most successful drug development programs for osteoporosis should focus on pathways identified through BMD-associated loci. © 2022 American Society for Bone and Mineral Research (ASBMR).

Polygenic risk score as a possible tool for identifying familial monogenic causes of complex diseases.

PURPOSE: The study aimed to evaluate whether polygenic risk scores could be helpful in addition to family history for triaging individuals to undergo deep-depth diagnostic sequencing for identifying monogenic causes of complex diseases. METHODS: Among 44,550 exome-sequenced European ancestry UK Biobank participants, we identified individuals with a clinically reported or computationally predicted monogenic pathogenic variant for breast cancer, bowel cancer, heart disease, diabetes, or Alzheimer disease. We derived polygenic risk scores for these diseases. We tested whether a polygenic risk score could identify rare pathogenic variant heterozygotes among individuals with a parental disease history. RESULTS: Monogenic causes of complex diseases were more prevalent among individuals with a parental disease history than in the rest of the population. Polygenic risk scores showed moderate discriminative power to identify familial monogenic causes. For instance, we showed that prescreening the patients with a polygenic risk score for type 2 diabetes can prioritize individuals to undergo diagnostic sequencing for monogenic diabetes variants and reduce needs for such sequencing by up to 37%. CONCLUSION: Among individuals with a family history of complex diseases, those with a low polygenic risk score are more likely to have monogenic causes of the disease and could be prioritized to undergo genetic testing.

Cohort profile: genomic data for 26 622 individuals from the Canadian Longitudinal Study on Aging (CLSA).

PURPOSE: The Canadian Longitudinal Study on Aging (CLSA) Comprehensive cohort was established to provide unique opportunities to study the genetic and environmental contributions to human disease as well as ageing process. The aim of this report was to describe the genomic data included in CLSA. PARTICIPANTS: A total of 26 622 individuals from the CLSA Comprehensive cohort of men and women aged 45-85 recruited between 2010 and 2015 underwent genome-wide genotyping of DNA samples collected from blood. Comprehensive quality control metrics were measured for genetic markers and samples, respectively. The genotypes were imputed to the TOPMed reference panel. Sex chromosome abnormalities were identified by copy number profiling. Classical human leukocyte antigen gene haplotypes were imputed at two-field (four-digit). FINDINGS TO DATE: Of the 26 622 genotyped participants, 24 655 (92.6%) were identified as having European ancestry. These genomic data were linked to physical, lifestyle, medical, economic, environmental and psychosocial factors collected longitudinally in CLSA. The combined analysis, including CLSA genomic data, uncovered over 100 novel loci associated with key parameters to define glaucoma. The CLSA genomic dataset validated the contribution of a polygenic risk score to screen individuals with high fracture risk. It is also a valuable resource to directly identify common genetic variations associated with conditions related to complex traits. Taking advantage of the comprehensive interview and physical information collected in CLSA, this genomic dataset has been linked to psychosocial factors to investigate both the independent and interactive effects on cardiovascular disease. FUTURE PLANS: The CLSA overall is ongoing. Follow-up data will continue to be collected from participants in the current genomic subcohort, including the DNA methylation and metabolomic data. Ongoing studies focus on elucidating the role of genetic factors in cognitive decline and cardiovascular diseases. This genomic data resource is available on request through the CLSA data access application process.

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.

Genetically increased circulating FUT3 level leads to reduced risk of idiopathic pulmonary fibrosis: a Mendelian randomisation study.

BACKGROUND: Idiopathic pulmonary fibrosis (IPF) is a progressive, fatal fibrotic interstitial lung disease. Few circulating biomarkers have been identified to have causal effects on IPF. METHODS: To identify candidate IPF-influencing circulating proteins, we undertook an efficient screen of circulating proteins by applying a two-sample Mendelian randomisation (MR) approach with existing publicly available data. For instruments, we used genetic determinants of circulating proteins which reside cis to the encoded gene (cis-single nucleotide polymorphisms (SNPs)), identified by two genome-wide association studies (GWASs) in European individuals (3301 and 3200 subjects). We then applied MR methods to test if the levels of these circulating proteins influenced IPF susceptibility in the largest IPF GWAS (2668 cases and 8591 controls). We validated the MR results using colocalisation analyses to ensure that both the circulating proteins and IPF shared a common genetic signal. RESULTS: MR analyses of 834 proteins found that a 1 sd increase in circulating galactoside 3(4)-l-fucosyltransferase (FUT3) and α-(1,3)-fucosyltransferase 5 (FUT5) was associated with a reduced risk of IPF (OR 0.81, 95% CI 0.74-0.88; p=6.3×10-7 and OR 0.76, 95% CI 0.68-0.86; p=1.1×10-5, respectively). Sensitivity analyses including multiple cis-SNPs provided similar estimates both for FUT3 (inverse variance weighted (IVW) OR 0.84, 95% CI 0.78-0.91; p=9.8×10-6 and MR-Egger OR 0.69, 95% CI 0.50-0.97; p=0.03) and FUT5 (IVW OR 0.84, 95% CI 0.77-0.92; p=1.4×10-4 and MR-Egger OR 0.59, 95% CI 0.38-0.90; p=0.01). FUT3 and FUT5 signals colocalised with IPF signals, with posterior probabilities of a shared genetic signal of 99.9% and 97.7%, respectively. Further transcriptomic investigations supported the protective effects of FUT3 for IPF. CONCLUSIONS: An efficient MR scan of 834 circulating proteins provided evidence that genetically increased circulating FUT3 level is associated with reduced risk of IPF.

Clinically Relevant Circulating Protein Biomarkers for Type 1 Diabetes: Evidence From a Two-Sample Mendelian Randomization Study.

OBJECTIVE: To identify circulating proteins influencing type 1 diabetes susceptibility using Mendelian randomization (MR). RESEARCH DESIGN AND METHODS: We used a large-scale two-sample MR study, using cis genetic determinants (protein quantitative trait loci [pQTL]) of up to 1,611 circulating proteins from five large genome-wide association studies, to screen for causal associations of these proteins with type 1 diabetes risk in 9,684 case subjects with type 1 diabetes and 15,743 control subjects. Further, pleiotropy-robust MR methods were used in sensitivity analyses using both cis and trans-pQTL. RESULTS: We found that a genetically predicted SD increase in signal regulatory protein gamma (SIRPG) level was associated with increased risk of type 1 diabetes risk (MR odds ratio [OR] 1.66 [95% 1.36-2.03]; P = 7.1 × 10-7). The risk of type 1 diabetes increased almost twofold per genetically predicted standard deviation (SD) increase in interleukin-27 Epstein-Barr virus-induced 3 (IL27-EBI3) protein levels (MR OR 1.97 [95% CI 1.48-2.62]; P = 3.7 × 10-6). However, an SD increase in chymotrypsinogen B1 (CTRB1) was associated with decreased risk of type 1 diabetes (MR OR 0.84 [95% CI 0.77-0.90]; P = 6.1 × 10-6). Sensitivity analyses using MR methods testing for pleiotropy while including trans-pQTL showed similar results. While the MR-Egger suggested no pleotropic effect (P value MR-Egger intercept = 0.31), there was evidence of pleiotropy in MR-PRESSO (P value global test = 0.006). CONCLUSIONS: We identified three novel circulating protein biomarkers associated with type 1 diabetes risk using an MR approach. These biomarkers are promising targets for development of drugs and/or of screening tools for early prediction of type 1 diabetes.

Health Effects of Calcium: Evidence From Mendelian Randomization Studies.

Calcium is widely used in conjunction with vitamin D to prevent osteoporosis. The use of calcium supplementation is also promoted for its potential benefits in lowering the risk for metabolic syndromes and cancers. However, the causal link between calcium and various health outcomes remains unclear. This review focuses on the evidence from 24 Mendelian randomization (MR) studies that were designed to minimize bias from confounding and reverse causation. These MR studies evaluated the effect of lifelong genetically higher serum calcium levels on various health outcomes. Overall, available MR studies found no conclusive effects of serum calcium levels on bone mineral density and fracture, ischemic stroke and heart failure, cancers, type 2 diabetes, Parkinson disease, or offspring birth weight. However, a higher serum calcium concentration was reported to have estimated causal effects on increased risks for coronary artery disease (especially myocardial infarction), migraine, renal colic, allergy/adverse effect of penicillin, and reduced risks for osteoarthrosis and osteoarthritis. In conclusion, supplementation of calcium in individuals from the general population is not predicted to influence the risk of most investigated diseases to date. Moreover, long-term high serum calcium concentrations may result in adverse health outcomes. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.