A common NFKB1 variant detected through antibody analysis in UK Biobank predicts risk of infection and allergy.

Infectious agents contribute significantly to the global burden of diseases through both acute infection and their chronic sequelae. We leveraged the UK Biobank to identify genetic loci that influence humoral immune response to multiple infections. From 45 genome-wide association studies in 9,611 participants from UK Biobank, we identified NFKB1 as a locus associated with quantitative antibody responses to multiple pathogens, including those from the herpes, retro-, and polyoma-virus families. An insertion-deletion variant thought to affect NFKB1 expression (rs28362491), was mapped as the likely causal variant and could play a key role in regulation of the immune response. Using 121 infection- and inflammation-related traits in 487,297 UK Biobank participants, we show that the deletion allele was associated with an increased risk of infection from diverse pathogens but had a protective effect against allergic disease. We propose that altered expression of NFKB1, as a result of the deletion, modulates hematopoietic pathways and likely impacts cell survival, antibody production, and inflammation. Taken together, we show that disruptions to the tightly regulated immune processes may tip the balance between exacerbated immune responses and allergy, or increased risk of infection and impaired resolution of inflammation.

Estimating disease heritability from complex pedigrees allowing for ascertainment and covariates.

We propose TetraHer, a method for estimating the liability heritability of binary phenotypes. TetraHer has five key features. First, it can be applied to data from complex pedigrees that contain multiple types of relationships. Second, it can correct for ascertainment of cases or controls. Third, it can accommodate covariates. Fourth, it can model the contribution of common environment. Fifth, it produces a likelihood that can be used for significance testing. We first demonstrate the validity of TetraHer on simulated data. We then use TetraHer to estimate liability heritability for 229 codes from the tenth International Classification of Diseases (ICD-10). We identify 107 codes with significant heritability (p < 0.05/229), which can be used in future analyses for investigating the genetic architecture of human diseases.

Evaluation of polygenic scoring methods in five biobanks shows larger variation between biobanks than methods and finds benefits of ensemble learning.

Methods of estimating polygenic scores (PGSs) from genome-wide association studies are increasingly utilized. However, independent method evaluation is lacking, and method comparisons are often limited. Here, we evaluate polygenic scores derived via seven methods in five biobank studies (totaling about 1.2 million participants) across 16 diseases and quantitative traits, building on a reference-standardized framework. We conducted meta-analyses to quantify the effects of method choice, hyperparameter tuning, method ensembling, and the target biobank on PGS performance. We found that no single method consistently outperformed all others. PGS effect sizes were more variable between biobanks than between methods within biobanks when methods were well tuned. Differences between methods were largest for the two investigated autoimmune diseases, seropositive rheumatoid arthritis and type 1 diabetes. For most methods, cross-validation was more reliable for tuning hyperparameters than automatic tuning (without the use of target data). For a given target phenotype, elastic net models combining PGS across methods (ensemble PGS) tuned in the UK Biobank provided consistent, high, and cross-biobank transferable performance, increasing PGS effect sizes (ß coefficients) by a median of 5.0% relative to LDpred2 and MegaPRS (the two best-performing single methods when tuned with cross-validation). Our interactively browsable online-results and open-source workflow prspipe provide a rich resource and reference for the analysis of polygenic scoring methods across biobanks.

A scalable and robust variance components method reveals insights into the architecture of gene-environment interactions underlying complex traits.

Understanding the contribution of gene-environment interactions (GxE) to complex trait variation can provide insights into disease mechanisms, explain sources of heritability, and improve genetic risk prediction. While large biobanks with genetic and deep phenotypic data hold promise for obtaining novel insights into GxE, our understanding of GxE architecture in complex traits remains limited. We introduce a method to estimate the proportion of trait variance explained by GxE (GxE heritability) and additive genetic effects (additive heritability) across the genome and within specific genomic annotations. We show that our method is accurate in simulations and computationally efficient for biobank-scale datasets. We applied our method to common array SNPs (MAF ≥1%), fifty quantitative traits, and four environmental variables (smoking, sex, age, and statin usage) in unrelated white British individuals in the UK Biobank. We found 68 trait-E pairs with significant genome-wide GxE heritability (p<0.05/200) with a ratio of GxE to additive heritability of ≈6.8% on average. Analyzing ≈8 million imputed SNPs (MAF ≥0.1%), we documented an approximate 28% increase in genome-wide GxE heritability compared to array SNPs. We partitioned GxE heritability across minor allele frequency (MAF) and local linkage disequilibrium (LD) values, revealing that, like additive allelic effects, GxE allelic effects tend to increase with decreasing MAF and LD. Analyzing GxE heritability near genes highly expressed in specific tissues, we find significant brain-specific enrichment for body mass index (BMI) and basal metabolic rate in the context of smoking and adipose-specific enrichment for waist-hip ratio (WHR) in the context of sex.

FastBiCmrMLM: a fast and powerful compressed variance component mixed logistic model for big genomic case-control genome-wide association study.

Large sample datasets have been regarded as the primary basis for innovative discoveries and the solution to missing heritability in genome-wide association studies. However, their computational complexity cannot consider all comprehensive effects and all polygenic backgrounds, which reduces the effectiveness of large datasets. To address these challenges, we included all effects and polygenic backgrounds in a mixed logistic model for binary traits and compressed four variance components into two. The compressed model combined three computational algorithms to develop an innovative method, called FastBiCmrMLM, for large data analysis. These algorithms were tailored to sample size, computational speed, and reduced memory requirements. To mine additional genes, linkage disequilibrium markers were replaced by bin-based haplotypes, which are analyzed by FastBiCmrMLM, named FastBiCmrMLM-Hap. Simulation studies highlighted the superiority of FastBiCmrMLM over GMMAT, SAIGE and fastGWA-GLMM in identifying dominant, small α (allele substitution effect), and rare variants. In the UK Biobank-scale dataset, we demonstrated that FastBiCmrMLM could detect variants as small as 0.03% and with α ≈ 0. In re-analyses of seven diseases in the WTCCC datasets, 29 candidate genes, with both functional and TWAS evidence, around 36 variants identified only by the new methods, strongly validated the new methods. These methods offer a new way to decipher the genetic architecture of binary traits and address the challenges outlined above.

Promises and Challenges of populational Proteomics in Health and Disease.

Advances in proteomic assay technologies have significantly increased coverage and throughput, enabling recent increases in the number of large-scale population-based proteomic studies of human plasma and serum. Improvements in multiplexed protein assays have facilitated the quantification of thousands of proteins over a large dynamic range, a key requirement for detecting the lowest-ranging, and potentially the most disease-relevant, blood-circulating proteins. In this perspective, we examine how populational proteomic datasets in conjunction with other concurrent omic measures can be leveraged to better understand the genomic and non-genomic correlates of the soluble proteome, constructing biomarker panels for disease prediction, among others. Mass spectrometry workflows are discussed as they are becoming increasingly competitive with affinity-based array platforms in terms of speed, cost, and proteome coverage due to advances in both instrumentation and workflows. Despite much success, there remain considerable challenges such as orthogonal validation and absolute quantification. We also highlight emergent challenges associated with study design, analytical considerations, and data integration as population-scale studies are run in batches and may involve longitudinal samples collated over many years. Lastly, we take a look at the future of what the nascent next-generation proteomic technologies might provide to the analysis of large sets of blood samples, as well as the difficulties in designing large-scale studies that will likely require participation from multiple and complex funding sources and where data sharing, study designs, and financing must be solved.

Exploring pleiotropy in Mendelian randomisation analyses: What are genetic variants associated with 'cigarette smoking initiation' really capturing?

Genetic variants used as instruments for exposures in Mendelian randomisation (MR) analyses may have horizontal pleiotropic effects (i.e., influence outcomes via pathways other than through the exposure), which can undermine the validity of results. We examined the extent of this using smoking behaviours as an example. We first ran a phenome-wide association study in UK Biobank, using a smoking initiation genetic instrument. From the most strongly associated phenotypes, we selected those we considered could either plausibly or not plausibly be caused by smoking. We examined associations between genetic instruments for smoking initiation, smoking heaviness and lifetime smoking and these phenotypes in UK Biobank and the Avon Longitudinal Study of Parents and Children (ALSPAC). We conducted negative control analyses among never smokers, including children. We found evidence that smoking-related genetic instruments were associated with phenotypes not plausibly caused by smoking in UK Biobank and (to a lesser extent) ALSPAC. We observed associations with phenotypes among never smokers. Our results demonstrate that smoking-related genetic risk scores are associated with unexpected phenotypes that are less plausibly downstream of smoking. This may reflect horizontal pleiotropy in these genetic risk scores, and we would encourage researchers to exercise caution this when using these and genetic risk scores for other complex behavioural exposures. We outline approaches that could be taken to consider this and overcome issues caused by potential horizontal pleiotropy, for example, in genetically informed causal inference analyses (e.g., MR) it is important to consider negative control outcomes and triangulation approaches, to avoid arriving at incorrect conclusions.

Are trait-associated genes clustered together in a gene network?

Genome-wide association studies (GWAS) have provided an abundance of information about the genetic variants and their loci that are associated to complex traits and diseases. However, due to linkage disequilibrium (LD) and noncoding regions of loci, it remains a challenge to pinpoint the causal genes. Gene network-based approaches, paired with network diffusion methods, have been proposed to prioritize causal genes and to boost statistical power in GWAS based on the assumption that trait-associated genes are clustered in a gene network. Due to the difficulty in mapping trait-associated variants to genes in GWAS, this assumption has never been directly or rigorously tested empirically. On the other hand, whole exome sequencing (WES) data focuses on the protein-coding regions, directly identifying trait-associated genes. In this study, we tested the assumption by leveraging the recently available exome-based association statistics from the UK Biobank WES data along with two types of networks. We found that almost all trait-associated genes were significantly more proximal to each other than randomly selected genes within both networks. These results support the assumption that trait-associated genes are clustered in gene networks, which can be further leveraged to boost the power of GWAS such as by introducing less stringent p value thresholds.

The AORTA Gene score for detection and risk stratification of ascending aortic dilation.

BACKGROUND AND AIMS: This study assessed whether a model incorporating clinical features and a polygenic score for ascending aortic diameter would improve diameter estimation and prediction of adverse thoracic aortic events over clinical features alone. METHODS: Aortic diameter estimation models were built with a 1.1 million-variant polygenic score (AORTA Gene) and without it. Models were validated internally in 4394 UK Biobank participants and externally in 5469 individuals from Mass General Brigham (MGB) Biobank, 1298 from the Framingham Heart Study (FHS), and 610 from All of Us. Model fit for adverse thoracic aortic events was compared in 401 453 UK Biobank and 164 789 All of Us participants. RESULTS: AORTA Gene explained more of the variance in thoracic aortic diameter compared to clinical factors alone: 39.5% (95% confidence interval 37.3%-41.8%) vs. 29.3% (27.0%-31.5%) in UK Biobank, 36.5% (34.4%-38.5%) vs. 32.5% (30.4%-34.5%) in MGB, 41.8% (37.7%-45.9%) vs. 33.0% (28.9%-37.2%) in FHS, and 34.9% (28.8%-41.0%) vs. 28.9% (22.9%-35.0%) in All of Us. AORTA Gene had a greater area under the receiver operating characteristic curve for identifying diameter ≥ 4 cm: 0.836 vs. 0.776 (P < .0001) in UK Biobank, 0.808 vs. 0.767 in MGB (P < .0001), 0.856 vs. 0.818 in FHS (P < .0001), and 0.827 vs. 0.791 (P = .0078) in All of Us. AORTA Gene was more informative for adverse thoracic aortic events in UK Biobank (P = .0042) and All of Us (P = .049). CONCLUSIONS: A comprehensive model incorporating polygenic information and clinical risk factors explained 34.9%-41.8% of the variation in ascending aortic diameter, improving the identification of ascending aortic dilation and adverse thoracic aortic events compared to clinical risk factors.

Incorporating functional annotation with bilevel continuous shrinkage for polygenic risk prediction.

BACKGROUND: Genetic variants can contribute differently to trait heritability by their functional categories, and recent studies have shown that incorporating functional annotation can improve the predictive performance of polygenic risk scores (PRSs). In addition, when only a small proportion of variants are causal variants, PRS methods that employ a Bayesian framework with shrinkage can account for such sparsity. It is possible that the annotation group level effect is also sparse. However, the number of PRS methods that incorporate both annotation information and shrinkage on effect sizes is limited. We propose a PRS method, PRSbils, which utilizes the functional annotation information with a bilevel continuous shrinkage prior to accommodate the varying genetic architectures both on the variant-specific level and on the functional annotation level. RESULTS: We conducted simulation studies and investigated the predictive performance in settings with different genetic architectures. Results indicated that when there was a relatively large variability of group-wise heritability contribution, the gain in prediction performance from the proposed method was on average 8.0% higher AUC compared to the benchmark method PRS-CS. The proposed method also yielded higher predictive performance compared to PRS-CS in settings with different overlapping patterns of annotation groups and obtained on average 6.4% higher AUC. We applied PRSbils to binary and quantitative traits in three real world data sources (the UK Biobank, the Michigan Genomics Initiative (MGI), and the Korean Genome and Epidemiology Study (KoGES)), and two sources of annotations: ANNOVAR, and pathway information from the Kyoto Encyclopedia of Genes and Genomes (KEGG), and demonstrated that the proposed method holds the potential for improving predictive performance by incorporating functional annotations. CONCLUSIONS: By utilizing a bilevel shrinkage framework, PRSbils enables the incorporation of both overlapping and non-overlapping annotations into PRS construction to improve the performance of genetic risk prediction. The software is available at https://github.com/styvon/PRSbils .

Causal association of blood lipids with all-cause and cause-specific mortality risk: a Mendelian randomization study.

Dyslipidemia has long been implicated in elevating mortality risk; yet, the precise associations between lipid traits and mortality remained undisclosed. Our study aimed to explore the causal effects of lipid traits on both all-cause and cause-specific mortality. One-sample Mendelian randomization (MR) with linear and nonlinear assumptions was conducted in a cohort of 407,951 European participants from the UK Biobank. Six lipid traits, consisting of low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), triglycerides, apolipoprotein A1 (ApoA1), apolipoprotein B (ApoB), and lipoprotein(a), were included to investigate the causal associations with mortality. Two-sample MR was performed to replicate the association between each lipid trait and all-cause mortality. Univariable MR results showed that genetically predicted higher ApoA1 was significantly associated with a decreased all-cause mortality risk (HR[95% CI]:0.93 [0.89-0.97], P value = 0.001), which was validated by the two-sample MR analysis. Higher lipoprotein(a) was associated with an increased risk of all-cause mortality (1.03 [1.01-1.04], P value = 0.002). Multivariable MR confirmed the direct causal effects of ApoA1 and lipoprotein(a) on all-cause mortality. Meanwhile, nonlinear MR found no evidence for nonlinearity between lipids and all-cause mortality. Our examination into cause-specific mortality revealed a suggestive inverse association between ApoA1 and cancer mortality, a significant positive association between lipoprotein(a) and cardiovascular disease mortality, and a suggestive positive association between lipoprotein(a) and digestive disease mortality. High LDL-C was associated with an increased risk of cardiovascular disease mortality but a decreased risk of neurodegenerative disease mortality. The findings suggest that implementing interventions to raise ApoA1 and decrease lipoprotein(a) levels may improve overall health outcomes and mitigate cancer and digestive disease mortality.

Metabolites, Healthy Lifestyle, and Polygenic Risk Score Associated with Upper Gastrointestinal Cancer: Findings from the UK Biobank Study.

Previous metabolomics studies have highlighted the predictive value of metabolites on upper gastrointestinal (UGI) cancer, while most of them ignored the potential effects of lifestyle and genetic risk on plasma metabolites. This study aimed to evaluate the role of lifestyle and genetic risk in the metabolic mechanism of UGI cancer. Differential metabolites of UGI cancer were identified using partial least-squares discriminant analysis and the Wilcoxon test. Then, we calculated the healthy lifestyle index (HLI) score and polygenic risk score (PRS) and divided them into three groups, respectively. A total of 15 metabolites were identified as UGI-cancer-related differential metabolites. The metabolite model (AUC = 0.699) exhibited superior discrimination ability compared to those of the HLI model (AUC = 0.615) and the PRS model (AUC = 0.593). Moreover, subgroup analysis revealed that the metabolite model showed higher discrimination ability for individuals with unhealthy lifestyles compared to that with healthy individuals (AUC = 0.783 vs 0.684). Furthermore, in the genetic risk subgroup analysis, individuals with a genetic predisposition to UGI cancer exhibited the best discriminative performance in the metabolite model (AUC = 0.770). These findings demonstrated the clinical significance of metabolic biomarkers in UGI cancer discrimination, especially in individuals with unhealthy lifestyles and a high genetic risk.

Determinants of plasma levels of proglucagon and the metabolic impact of glucagon receptor signalling: a UK Biobank study.

AIMS/HYPOTHESES: Glucagon and glucagon-like peptide-1 (GLP-1) are derived from the same precursor; proglucagon, and dual agonists of their receptors are currently being explored for the treatment of obesity and metabolic dysfunction-associated steatotic liver disease (MASLD). Elevated levels of endogenous glucagon (hyperglucagonaemia) have been linked with hyperglycaemia in individuals with type 2 diabetes but are also observed in individuals with obesity and MASLD. GLP-1 levels have been reported to be largely unaffected or even reduced in similar conditions. We investigated potential determinants of plasma proglucagon and associations of glucagon receptor signalling with metabolic diseases based on data from the UK Biobank. METHODS: We used exome sequencing data from the UK Biobank for ~410,000 white participants to identify glucagon receptor variants and grouped them based on their known or predicted signalling. Data on plasma levels of proglucagon estimated using Olink technology were available for a subset of the cohort (~40,000). We determined associations of glucagon receptor variants and proglucagon with BMI, type 2 diabetes and liver fat (quantified by liver MRI) and performed survival analyses to investigate if elevated proglucagon predicts type 2 diabetes development. RESULTS: Obesity, MASLD and type 2 diabetes were associated with elevated plasma levels of proglucagon independently of each other. Baseline proglucagon levels were associated with the risk of type 2 diabetes development over a 14 year follow-up period (HR 1.13; 95% CI 1.09, 1.17; n=1562; p=1.3×10-12). This association was of the same magnitude across strata of BMI. Carriers of glucagon receptor variants with reduced cAMP signalling had elevated levels of proglucagon (ß 0.847; 95% CI 0.04, 1.66; n=17; p=0.04), and carriers of variants with a predicted frameshift mutation had higher levels of liver fat compared with the wild-type reference group (ß 0.504; 95% CI 0.03, 0.98; n=11; p=0.04). CONCLUSIONS/INTERPRETATION: Our findings support the suggestion that glucagon receptor signalling is involved in MASLD, that plasma levels of proglucagon are linked to the risk of type 2 diabetes development, and that proglucagon levels are influenced by genetic variation in the glucagon receptor, obesity, type 2 diabetes and MASLD. Determining the molecular signalling pathways downstream of glucagon receptor activation may guide the development of biased GLP-1/glucagon co-agonist with improved metabolic benefits. DATA AVAILABILITY: All coding is available through https://github.com/nicwin98/UK-Biobank-GCG.

Physical Activity, Mental Activity, and Risk of Incident Stroke: A Prospective Cohort Study.

BACKGROUND: Cumulative evidence suggests a correlation between physical or mental activity and the risk of stroke. However, the combined impact of these activities on stroke onset remains unexplored. This study identified physical and mental activity patterns using principal component analysis and investigated their associations with risk of incident stroke in the general population. METHODS: Our study was sourced from the UK Biobank cohort between 2006 and 2010. Information on physical and mental-related activities were obtained through a touch-screen questionnaire. The incident stroke was diagnosed by physicians and subsequently verified through linkage to Hospital Episode Statistics. Principal component analysis was used to identify potential physical and mental activity patterns. Cox proportional hazard regression models were performed to calculate hazard ratios (HRs) and 95% CIs of incident stroke, adjusting for potential confounders. RESULTS: The initial UK Biobank cohort originally consisted of 502 411 individuals, of whom a total of 386 902 participants (aged 38-79 years) without any history of stroke at baseline were included in our study. During a median follow-up of 7.7 years, 6983 (1.8%) cases of stroke were documented. The mean age of the included participants was 55.9 years, and the proportion of women was 55.1%. We found that multiple individual items related to physical and mental activity showed significant associations with risk of stroke. We identified 4 patterns of physical activity and 3 patterns of mental activity using principal component analysis. The adherence to activity patterns of vigorous exercise, housework, and walking predominant patterns were associated with a lower risk of stroke by 17% (HR, 0.83 [95% CI, 0.78-0.89]; 20% (HR, 0.80 [95% CI, 0.75-0.85]; and 20% (HR, 0.80 [95% CI, 0.75-0.86), respectively. Additionally, the transportation predominant pattern (HR, 1.36 [95% CI, 1.28-1.45) and watching TV pattern (HR, 1.43 [95% CI, 1.33-1.53) were found to be significantly associated with a higher risk of stroke. These associations remained consistent across all subtypes of stroke. CONCLUSIONS: Activity patterns mainly related to frequent vigorous exercise, housework, and walking were associated with lower risks of stroke and all its subtypes. Our findings provide new insights for promoting suitable patterns of physical and mental activity for primary prevention of stroke.

Genetic variants affecting mitochondrial function provide further insights for kidney disease.

BACKGROUND: Chronic kidney disease (CKD) is a complex disorder that has become a high prevalence global health problem, with diabetes being its predominant pathophysiologic driver. Autosomal genetic variation only explains some of the predisposition to kidney disease. Variations in the mitochondrial genome (mtDNA) and nuclear-encoded mitochondrial genes (NEMG) are implicated in susceptibility to kidney disease and CKD progression, but they have not been thoroughly explored. Our aim was to investigate the association of variation in both mtDNA and NEMG with CKD (and related traits), with a particular focus on diabetes. METHODS: We used the UK Biobank (UKB) and UK-ROI, an independent collection of individuals with type 1 diabetes mellitus (T1DM) patients. RESULTS: Fourteen mitochondrial variants were associated with estimated glomerular filtration rate (eGFR) in UKB. Mitochondrial variants and haplogroups U, H and J were associated with eGFR and serum variables. Mitochondrial haplogroup H was associated with all the serum variables regardless of the presence of diabetes. Mitochondrial haplogroup X was associated with end-stage kidney disease (ESKD) in UKB. We confirmed the influence of several known NEMG on kidney disease and function and found novel associations for SLC39A13, CFL1, ACP2 or ATP5G1 with serum variables and kidney damage, and for SLC4A1, NUP210 and MYH14 with ESKD. The G allele of TBC1D32-rs113987180 was associated with higher risk of ESKD in patients with diabetes (OR:9.879; CI95%:4.440-21.980; P = 2.0E-08). In UK-ROI, AGXT2-rs71615838 and SURF1-rs183853102 were associated with diabetic nephropathies, and TFB1M-rs869120 with eGFR. CONCLUSIONS: We identified novel variants both in mtDNA and NEMG which may explain some of the missing heritability for CKD and kidney phenotypes. We confirmed the role of MT-ND5 and mitochondrial haplogroup H on renal disease (serum variables), and identified the MT-ND5-rs41535848G variant, along with mitochondrial haplogroup X, associated with higher risk of ESKD. Despite most of the associations were independent of diabetes, we also showed potential roles for NEMG in T1DM.

Genome-wide association study identifies novel susceptible loci and evaluation of polygenic risk score for chronic obstructive pulmonary disease in a Taiwanese population.

BACKGROUND: Chronic Obstructive Pulmonary Disease (COPD) describes a group of progressive lung diseases causing breathing difficulties. While COPD development typically involves a complex interplay between genetic and environmental factors, genetics play a role in disease susceptibility. This study used genome-wide association studies (GWAS) and polygenic risk score (PRS) to elucidate the genetic basis for COPD in Taiwanese patients. RESULTS: GWAS was performed on a Taiwanese COPD case-control cohort with a sample size of 5,442 cases and 17,681 controls. Additionally, the PRS was calculated and assessed in our target groups. GWAS results indicate that although there were no single nucleotide polymorphisms (SNPs) of genome-wide significance, prominent COPD susceptibility loci on or nearby genes such as WWTR1, EXT1, INTU, MAP3K7CL, MAMDC2, BZW1/CLK1, LINC01197, LINC01894, and CFAP95 (C9orf135) were identified, which had not been reported in previous studies. Thirteen susceptibility loci, such as CHRNA4, AFAP1, and DTWD1, previously reported in other populations were replicated and confirmed to be associated with COPD in Taiwanese populations. The PRS was determined in the target groups using the summary statistics from our base group, yielding an effective association with COPD (odds ratio [OR] 1.09, 95% confidence interval [CI] 1.02-1.17, p = 0.011). Furthermore, replication a previous lung function trait PRS model in our target group, showed a significant association of COPD susceptibility with PRS of Forced Expiratory Volume in one second (FEV1)/Forced Vital Capacity (FCV) (OR 0.89, 95% CI 0.83-0.95, p = 0.001). CONCLUSIONS: Novel COPD-related genes were identified in the studied Taiwanese population. The PRS model, based on COPD or lung function traits, enables disease risk estimation and enhances prediction before suffering. These results offer new perspectives on the genetics of COPD and serve as a basis for future research.

Association of Life's Essential 8 cardiovascular health with breast cancer incidence and mortality according to genetic susceptibility of breast cancer: a prospective cohort study.

BACKGROUND: Accumulating evidence suggests that cardiovascular diseases and breast cancer share a number of common risk factors, however, evidence on the association between cardiovascular health (CVH) and breast cancer is limited. The present study aimed to assess the association of CVH, defined by Life's Essential 8 (LE8) and genetic risk with breast cancer incidence and mortality among premenopausal and postmenopausal women. METHODS: We used data from the UK Biobank and conducted the multivariate Cox proportional-hazards models to examine associations of LE8 score and genetic risk with breast cancer incidence and mortality. Date on LE8 score was collected between 2006 and 2010 and composed of eight components, including behavioral metrics (diet, tobacco or nicotine exposure, physical activity, and sleep health), and biological metrics (body mass index, blood lipids, blood glucose, and blood pressure). The polygenic risk score (PRS) was calculated as the sum of effect sizes of individual genetic variants multiplied by the allele dosage. RESULTS: A total of 150,566 premenopausal and postmenopausal women were included. Compared to postmenopausal women with low LE8 score, those with high LE8 score were associated with 22% lower risk of breast cancer incidence (HR: 0.78, 95% CI: 0.70-0.87) and 43% lower risk of breast cancer mortality (HR: 0.57, 95% CI: 0.36-0.90). By contrast, we did not observe the significant association among premenopausal women. Further analyses stratified by PRS categories showed that high LE8 score was associated with 28% and 71% decreased risk of breast cancer incidence (HR: 0.72, 95% CI: 0.60-0.87) and mortality (HR: 0.29, 95% CI: 0.10-0.83) compared to low LE8 score among high genetic risk groups, but no significant associations were found among low genetic risk groups. Furthermore, compared with postmenopausal women with high LE8 score and low genetic risk, those with low LE8 score and high genetic risk were associated with increased risk of breast cancer incidence (HR: 6.26, 95% CI: 4.43-8.84). CONCLUSIONS: The present study suggests that better CVH is a protective factor for both breast cancer incidence and mortality among postmenopausal women. Moreover, the risk of developing breast cancer caused by high genetic susceptibility could be largely offset by better CVH.

Circulating white blood cell traits and colorectal cancer risk: A Mendelian randomisation study.

Observational studies have suggested a protective role for eosinophils in colorectal cancer (CRC) development and implicated neutrophils, but the causal relationships remain unclear. Here, we aimed to estimate the causal effect of circulating white blood cell (WBC) counts (N = ~550 000) for basophils, eosinophils, monocytes, lymphocytes and neutrophils on CRC risk (N = 52 775 cases and 45 940 controls) using Mendelian randomisation (MR). For comparison, we also examined this relationship using individual-level data from UK Biobank (4043 incident CRC cases and 332 773 controls) in a longitudinal cohort analysis. The inverse-variance weighted (IVW) MR analysis suggested a protective effect of increased basophil count and eosinophil count on CRC risk [OR per 1-SD increase: 0.88, 95% CI: 0.78-0.99, P = .04; OR: 0.93, 95% CI: 0.88-0.98, P = .01]. The protective effect of eosinophils remained [OR per 1-SD increase: 0.88, 95% CI: 0.80-0.97, P = .01] following adjustments for all other WBC subtypes, to account for genetic correlation between the traits, using multivariable MR. A protective effect of increased lymphocyte count on CRC risk was also found [OR: 0.84, 95% CI: 0.76-0.93, P = 6.70e-4] following adjustment. Consistent with MR results, a protective effect for eosinophils in the cohort analysis in the fully adjusted model [RR per 1-SD increase: 0.96, 95% CI: 0.93-0.99, P = .02] and following adjustment for the other WBC subtypes [RR: 0.96, 95% CI: 0.93-0.99, P = .001] was observed. Our study implicates peripheral blood immune cells, in particular eosinophils and lymphocytes, in CRC development, highlighting a need for mechanistic studies to interrogate these relationships.

A large study of metabolomics reveals common and distinct metabolic biomarkers for type 2 diabetes, coronary heart disease, and stroke.

We aimed at examining the shared and unique associations of metabolites with multiple cardiometabolic diseases (CMD), i.e. type 2 diabetes (T2D), coronary heart disease (CHD) and stroke. In this study, a total of 168 plasma metabolites were measured by targeted high-throughput nuclear magnetic resonance spectroscopy among 98,162 participants free of T2D, CHD, and stroke at baseline. Cox proportional hazard models estimated hazard ratios for one SD increase in metabolite concentration levels, and false discovery rate (at 10%) was used to correct for multiple comparisons. Over 12.1 years of follow-up on average, 3,463 T2D, 6,186 CHD, and 1,892 stroke events were recorded. Most lipoprotein metabolites were associated with risks of T2D and CHD but not with the risk of stroke, with stronger associations for T2D than for CHD. Phospholipids within intermediate-density lipoprotein or large low-density lipoprotein particles showed positive associations with CHD and inverse associations with T2D. Metabolites indicating very small very low-density lipoprotein, histidine, creatinine, albumin, and glycoprotein acetyls were associated with risks of all three conditions. This large-scale metabolomics study revealed common and distinct metabolic biomarkers for T2D, CHD and stroke, providing instrumental information to possibly implement precision medicine for preventing and treating these conditions.

Joint Modeling of Gene-Environment Correlations and Interactions Using Polygenic Risk Scores in Case-Control Studies.

Polygenic risk scores (PRS) are rapidly emerging as a way to measure disease risk by aggregating multiple genetic variants. Understanding the interplay of PRS with environmental factors is critical for interpreting and applying PRS in a wide variety of settings. We develop an efficient method for simultaneously modeling gene-environment correlations and interactions using PRS in case control studies. We use a logistic-normal regression modeling framework to specify the disease risk and PRS distribution in the underlying population and propose joint inference across the two models using the retrospective likelihood of the case-control data. Extensive simulation studies demonstrate the flexibility of the method in trading-off bias and efficiency for the estimation of various model parameters compared to the standard logistic regression or a case-only analysis for gene-environment interactions, or a control-only analysis for gene-environment correlations. Finally using simulated case-control data sets within the UK Biobank study, we demonstrate the power of our method for its ability to recover results from the full prospective cohort for the detection of an interaction between long-term oral contraceptive use and PRS on the risk of breast cancer. This method is computationally efficient and implemented in a user-friendly R package.