METAL: fast and efficient meta-analysis of genomewide association scans

Genome-Wide Association Study of Temporomandibular Disorder-Related Pain in Finnish Populations

BACKGROUND

Temporomandibular disorders (TMD) are multifactorial musculoskeletal pain and dysfunctions in temporomandibular joints (TMJs) and masticatory muscles. Genetic factors play a role in TMD-related pain, but only a few genome-wide association studies (GWAS) have been conducted.

OBJECTIVE

The aim of this GWAS was to explore genetic factors associated with painful TMD in Finnish populations.

METHODS

Data from two epidemiological surveys, the Northern Finland Birth Cohort 1966 (NFBC1966) and the Health 2000 Survey in Finland, including altogether 468 cases and 6833 controls, were used. Case definition was based on pain on palpation of masticatory muscles and/or TMJs. GWASs of the whole data and stratified by sex were conducted from both cohorts using additive models, followed by meta-analysis of the two cohorts. Replications of the previously reported TMD risk loci (rs73460075, DMD; rs4794106, SGCA; rs73271865, SP4; rs60249166, RXP2; rs1531554, BAHCCI; rs5862730, OTUD4/SMAD1; rs10092633, SFRP1; rs34612513, SOX14/CLDN18; rs878962, TSPAN9) were also investigated.

RESULTS

Four genome-wide significant loci were found in sex-stratified analysis of NFBC1966, including associations at three loci in males (rs1023114, PRIM2, p = 5 × 10-9; rs4244867, ALG10, p = 3 × 10-8; rs79841648, ADCYAP1, p = 4 × 10-9) and one locus in females (rs148476652, DNER, p = 4 × 10-9). However, the results could not be replicated in the Health 2000 Survey or in the meta-analysis of these two cohorts. The previous TMD GWAS associations did not replicate in our data either.

CONCLUSION

Several TMD pain-associated variants were found in sex-stratified analysis of NFBC1966, suggesting the role of neuroendocrine stress responses and central nervous system. These findings need to be confirmed in future studies.

Genetics of C-Peptide and Age at Diagnosis in Type 1 Diabetes

ARTICLE HIGHLIGHTS

Identified genetic loci for C-peptide and type 1 diabetes (T1D) age at diagnosis (AAD) explain only a small proportion of their variation. We aimed to identify additional genetic loci associated with C-peptide and AAD. Some HLA allele/haplotypes associated with T1D also contributed to variability of C-peptide and AAD, whereas outside the HLA region, T1D loci were mostly not associated with C-peptide or AAD. Genetic variation within CTSH can affect AAD. There is still residual heritability of C-peptide and AAD outside of HLA that could benefit from larger meta-genome-wide association studies.

Genetic risk for trait aggression and alcohol use predict unique facets of alcohol-related aggression

OBJECTIVE

A propensity for aggression or alcohol use may be associated with alcohol-related aggression. Previous research has shown genetic overlap between alcohol use and aggression but has not looked at how alcohol-related aggression may be uniquely influenced by genetic risk for aggression or alcohol use. The present study examined the associations of genetic risk for trait aggression, alcohol use, and alcohol use disorder (AUD) with alcohol-related aggression using a polygenic risk score (PRS) approach.

METHOD

Using genome-wide association study summary statistics, PRSs were created for trait aggression, alcohol consumption, and AUD. These PRSs were used to predict the phenotype of alcohol-related aggression among drinkers in two independent samples: the University of California at San Francisco (UCSF) Family Alcoholism Study (n = 1,162) and the National Longitudinal Study of Adolescent to Adult Health (Add Health; n = 4,291).

RESULTS

There were significant associations between the AUD PRS and lifetime alcohol-related aggression in the UCSF study sample. Additionally, the trait aggression PRS was associated with three or more experiences of hitting anyone else and getting into physical fights while under the influence of alcohol, along with a composite score of three or more experiences of alcohol-related aggression, in the UCSF study sample. No significant associations were observed in the Add Health sample. Limited sex-specific genetic effects were observed.

CONCLUSIONS

These results provide preliminary evidence that genetic influences underlying alcohol use and aggression are uniquely associated with alcohol-related aggression and suggest that these associations may differ by type and frequency of alcohol-related aggression incidents. (PsycInfo Database Record (c) 2025 APA, all rights reserved).

Genome-wide association study of prostate-specific antigen levels in 392,522 men identifies new loci and improves prediction across ancestry groups

We conducted a multiancestry genome-wide association study of prostate-specific antigen (PSA) levels in 296,754 men (211,342 European ancestry, 58,236 African ancestry, 23,546 Hispanic/Latino and 3,630 Asian ancestry; 96.5% of participants were from the Million Veteran Program). We identified 318 independent genome-wide significant (P ≤ 5 × 10-8) variants, 184 of which were novel. Most demonstrated evidence of replication in an independent cohort (n = 95,768). Meta-analyzing discovery and replication (n = 392,522) identified 447 variants, of which a further 111 were novel. Out-of-sample variance in PSA explained by our genome-wide polygenic risk scores ranged from 11.6% to 16.6% for European ancestry, 5.5% to 9.5% for African ancestry, 13.5% to 18.2% for Hispanic/Latino and 8.6% to 15.3% for Asian ancestry and decreased with increasing age. Midlife genetically adjusted PSA levels were more strongly associated with overall and aggressive prostate cancer than unadjusted PSA levels. Our study highlights how including proportionally more participants from underrepresented populations improves genetic prediction of PSA levels, offering potential to personalize prostate cancer screening.

Identification of 16 novel Alzheimer's disease loci using multi-ancestry meta-analyses

INTRODUCTION

Alzheimer's disease (AD) is the most prevalent form of dementia. While many AD-associated genetic determinants have been identified, few studies have analyzed individuals of non-European ancestry.

METHODS

We conducted a multi-ancestry genome-wide association study (GWAS) of clinically diagnosed AD and AD-by-proxy using whole genome sequencing data from the National Institute on Aging Genetics of Alzheimer's Disease Data Storage Site (NIAGADS), National Institute of Mental Health, UK Biobank (UKB), and All of Us (AoU) consisting of 49,149 cases (12,074 clinically diagnosed and 37,075 AD-by-proxy) and 383,225 controls. Nearly half of NIAGADS and AoU participants were of non-European ancestry.

RESULTS

For clinically diagnosed AD, we identified 14 new loci-five common (FBN2/SCL27A6, AC090115.1, DYM, KCNG1/AL121785.1, TIAM1) and nine rare (VWA5B1, RNU6-755P/LMX1A, MOB1A, MORC1-AS1, LINC00989, PDE4D, RNU2-49P/CDO1, NEO1, and SLC35G3/AC022916.1). Meta-analysis of UKB and AoU AD-by-proxy cases yielded two new rare loci (RPL23/LASP1 and CEBPA/AC008738.6), also nominally significant in NIAGADS.

DISCUSSION

In summary, we provide evidence for 16 novel AD loci and advocate for more studies using whole genome sequencing-based GWAS of diverse cohorts.

HIGHLIGHTS

We used whole-genome sequencing data from large and diverse cohorts. We found novel genome-wide association study findings based on whole-genome data. We performed a multiancestry meta-analysis and incorporated results from underrepresented groups.

Distinct Genetic Risk Profile in Aortic Stenosis Compared With Coronary Artery Disease

Importance

Aortic stenosis (AS) and coronary artery disease (CAD) frequently coexist. However, it is unknown which genetic and cardiovascular risk factors might be AS-specific and which could be shared between AS and CAD.

Objective

To identify genetic risk loci and cardiovascular risk factors with AS-specific associations.

Design, Setting, and Participants

This was a genomewide association study (GWAS) of AS adjusted for CAD with participants from the European Consortium for the Genetics of Aortic Stenosis (EGAS) (recruited 2000-2020), UK Biobank (recruited 2006-2010), Estonian Biobank (recruited 1997-2019), and FinnGen (recruited 1964-2019). EGAS participants were collected from 7 sites across Europe. All participants were of European ancestry, and information on comorbid CAD was available for all participants. Follow-up analyses with GWAS data on cardiovascular traits and tissue transcriptome data were also performed. Data were analyzed from October 2022 to July 2023.

Exposures

Genetic variants.

Main Outcomes and Measures

Cardiovascular traits associated with AS adjusted for CAD. Replication was performed in 2 independent AS GWAS cohorts.

Results

A total of 18 792 participants with AS and 434 249 control participants were included in this GWAS adjusted for CAD. The analysis found 17 AS risk loci, including 5 loci with novel and independently replicated associations (RNF114A, AFAP1, PDGFRA, ADAMTS7, HAO1). Of all 17 associated loci, 11 were associated with risk specifically for AS and were not associated with CAD (ALPL, PALMD, PRRX1, RNF144A, MECOM, AFAP1, PDGFRA, IL6, TPCN2, NLRP6, HAO1). Concordantly, this study revealed only a moderate genetic correlation of 0.15 (SE, 0.05) between AS and CAD (P = 1.60 × 10-3). Mendelian randomization revealed that serum phosphate was an AS-specific risk factor that was absent in CAD (AS: odds ratio [OR], 1.20; 95% CI, 1.11-1.31; P = 1.27 × 10-5; CAD: OR, 0.97; 95% CI 0.94-1.00; P = .04). Mendelian randomization also found that blood pressure, body mass index, and cholesterol metabolism had substantially lesser associations with AS compared with CAD. Pathway and transcriptome enrichment analyses revealed biological processes and tissues relevant for AS development.

Conclusions and Relevance

This GWAS adjusted for CAD found a distinct genetic risk profile for AS at the single-marker and polygenic level. These findings provide new targets for future AS research.

Novel risk loci encompassing genes influencing STAT3, GPCR, and oxidative stress signaling are associated with co-morbid GERD and COPD

Chronic obstructive pulmonary disease (COPD) is a leading cause of death globally. Gastroesophageal reflux disease (GERD) is a common comorbidity in COPD associated with worse pulmonary symptoms, reduced quality of life, and increased exacerbations and hospitalizations. GERD treatment in COPD is associated with a lower risk of exacerbations and mortality; however, it is not clear whether these findings can be attributed to aging populations where both diseases are likely to co-occur or reflect shared etiology. To test for the influence of common etiology in both diseases, we aimed to identify shared genetic etiology between GERD and COPD. We performed the first whole-genome sequence association analysis of comorbid GERD and COPD in 12,438 multi-ancestry participants. The co-heritability of GERD and COPD was 39.7% (h2 = 0.397, SE = 0.074) and we identified several ancestry-independent loci associated with co-morbid GERD and COPD (within LINC02493 and FRYL) known to be involved in oxidative stress and G protein-coupled receptor (GPCR) signaling mechanisms. We found several loci associated with co-morbid GERD and COPD previously associated with GERD or COPD individually, including HCG17, which plays a role in oxidative stress mechanisms. Gene set enrichment identified GPCR signaling pathways in co-morbid GERD and COPD loci. Rare variants in ZFP42, encoding key regulators of the IL6/STAT3 pathway, have been previously implicated with GI disorders and were associated with co-morbid GERD and COPD. We identified common genetic etiology for GERD in COPD which begins to provide a mechanistic foundation for the potential therapeutic utility of STAT3, oxidation, and GPCR signaling pathway modulators in both GERD and COPD.

Genetic Discovery and Risk Prediction for Type 1 Diabetes in Individuals Without High-Risk HLA-DR3/DR4 Haplotypes

OBJECTIVE

More than 10% of patients with type 1 diabetes (T1D) do not have high-risk HLA-DR3 or -DR4 haplotypes with distinct clinical features, such as later onset and reduced insulin dependence. We aimed to identify genetic drivers of T1D in the absence of DR3/DR4 and improve prediction of T1D risk in these individuals.

RESEARCH DESIGN AND METHODS

We performed T1D association and fine-mapping analyses in 12,316 non-DR3/DR4 samples. Next, we performed heterogeneity tests to examine differences in T1D risk variants in individuals without versus those with DR3/DR4 haplotypes. We further assessed genome-wide differences in gene regulatory element and biological pathway enrichments between the non-DR3/DR4 and DR3/DR4 cohorts. Finally, we developed a genetic risk score (GRS) to predict T1D in individuals without DR3/DR4 and compared with an existing T1D GRS.

RESULTS

A total of 18 T1D risk variants in non-DR3/DR4 samples were identified. Risk variants at the MHC and multiple other loci genome wide had heterogeneity in effects on T1D dependent on DR3/DR4 status, and non-DR3/DR4 T1D had evidence for a greater polygenic burden. T1D-associated variants in non-DR3/DR4 were more enriched for regulatory elements and pathways involved in antigen presentation, innate immunity, and β-cells and depleted in T cells compared with DR3/DR4. A non-DR3/DR4 GRS outperformed an existing risk score GRS2 in discriminating non-DR3/DR4 T1D from no diabetes (area under the curve 0.867; P = 7.48 × 10-32) and type 2 diabetes (0.907; P = 4.94 × 10-44).

CONCLUSIONS

In total, we identified heterogeneity in T1D genetic risk dependent on high-risk HLA-DR3/DR4 haplotype, which uncovers disease mechanisms and enables more accurate prediction of T1D across the HLA background.

The genetic architecture of hip shape and its role in the development of hip osteoarthritis and fracture

OBJECTIVES

Hip shape is thought to be an important causal risk factor for hip osteoarthritis and fracture. We aimed to identify genetic determinants of hip shape and use these to assess causal relationships with hip osteoarthritis.

METHODS

Statistical hip shape modelling was used to derive 10 hip shape modes (HSMs) from DXA images in UK Biobank and Shanghai Changfeng cohorts (ntotal = 43 485). Genome-wide association study meta-analyses were conducted for each HSM. Two-sample Mendelian randomisation (MR) was used to estimate causal effects between HSM and hip osteoarthritis using hip fracture as a positive control.

RESULTS

Analysis of the first 10 HSMs identified 203 independent association signals (P < 5 × 10-9). Hip shape SNPs were also associated (P < 2.5 × 10-4) with hip osteoarthritis (n = 26) and hip fracture (n = 4). Fine mapping implicated SMAD3 and PLEC as candidate genes that may be involved in the development of hip shape and hip osteoarthritis. MR analyses suggested there was no causal effect between any HSM and hip osteoarthritis, however there was evidence that HSM2 (more obtuse neck-shaft angle) and HSM4 (wider femoral neck) have a causal effect on hip fracture (ORIVW method 1.27 [95% CI 1.12-1.44], P = 1.79 × 10-4 and ORIVW 0.74 [0.65-0.84], P = 7.60 × 10-6 respectively).

CONCLUSIONS

We report the largest hip shape GWAS meta-analysis that identifies hundreds of novel loci, some of which are also associated with hip osteoarthritis and hip fracture. MR analyses suggest hip shape may not cause hip osteoarthritis but is implicated in hip fractures. Consequently, interventions targeting hip shape in older adults to prevent hip osteoarthritis may prove ineffective.

Network Mendelian randomisation analysis deciphers protein pathways linking type 2 diabetes and gastrointestinal disease

AIMS

The molecular mechanisms underlying the association between type 2 diabetes (T2D) and gastrointestinal (GI) disease are unclear. To identify protein pathways, we conducted a two-stage network Mendelian randomisation (MR) study.

MATERIALS AND METHODS

Genetic instruments for T2D were obtained from a large-scale summary-level genome-wide meta-analysis. Genetic associations with blood protein levels were obtained from three genome-wide association studies on plasma proteins (i.e. the deCODE study as the discovery and the UKB-PPP and Fenland studies as the replication). Summary-level data on 10 GI diseases were derived from genome-wide meta-analysis of the UK Biobank and FinnGen. MR and colocalisation analyses were performed. Pathways were constructed according to the directionality of total and indirect effects, and corresponding proportional mediation was estimated. Druggability assessments were conducted across four databases to prioritise protein mediators.

RESULTS

Genetic liability to T2D was associated with 69 proteins in the discovery protein dataset after multiple testing corrections. All associations were replicated at the nominal significance level. Among T2D-associated proteins, genetically predicted levels of nine proteins were associated with at least one of the GI diseases. Genetically predicted levels of SULT2A1 (odds ratio = 1.98, 95% CI 1.80-2.18), and ADH1B (odds ratio = 2.05, 95% CI 1.43-2.94) were associated with cholelithiasis and cirrhosis respectively. SULT2A1 and cholelithiasis (PH4 = 0.996) and ADH1B and cirrhosis (PH4 = 0.931) have strong colocalisation support, accounting for the mediation proportion of 72.8% (95% CI 45.7-99.9) and 42.9% (95% CI 15.5-70.4) respectively.

CONCLUSIONS

The study identified some proteins mediating T2D-GI disease associations, which provided biological insights into the underlying pathways.

GWAS of CRP response to statins further supports the role of APOE in statin response: A GIST consortium study

Statins are first-line treatments in the primary and secondary prevention of cardiovascular disease. Clinical studies show statins act independently of lipid-lowering mechanisms to decrease C-reactive protein (CRP), an inflammation marker. We aim to elucidate genetic loci associated with CRP statin response. CRP statin response is the change in log-CRP between off-treatment and on-treatment measurements. Cohort-level Genome-Wide Association Studies (GWAS) of CRP response were performed using 1000 Genomes imputed data, testing ∼10 million common genetic variants. GWAS meta-analysis combined results from seven cohorts and clinical trials totalling 14,070 statin-treated individuals of European ancestry within the GIST consortium. Secondary analyses included statin-by-placebo interaction analyses, and lookups in African ancestry cohorts. Our GWAS identified two genome-wide significant (P < 5e-8) loci: APOE and HNF1A for CRP statin response corrected for baseline CRP. The missense lead variant rs429358 at APOE, contributing to the APOE-E4 haplotype, is a risk locus for dyslipidaemia, Alzheimer's and coronary artery disease (CAD). The HNF1A locus is associated with diabetes, cholesterol levels, and CAD. Both loci are also associated with baseline CRP levels, and neither locus achieved a significant (P < 0.05) result from the statin v. placebo interaction meta-analysis using randomized clinical trial data. However, the interaction result (P-int=0.09) for APOE was suggestive and possibly underpowered. The APOE-E4 signal may therefore be associated with both CRP and LDL-cholesterol statin response. Combined with suggestions in the literature that APOE also leads to differential statin benefit in Alzheimer's, the APOE locus warrants further investigation for potential genetic effects on healthcare with statin treatment.

Multi-ancestry meta-analysis of keloids uncovers novel susceptibility loci in diverse populations

Keloids are raised scars that grow beyond original wound boundaries, resulting in pain and disfigurement. Reasons for keloid development are not well-understood, and current treatment options are limited. Keloids are more likely to occur in darker-skinned individuals of African and Asian descent than in Europeans. We performed a genome-wide association study (GWAS) examining keloid risk across and within continental ancestry groups, incorporating 7,837 cases and 1,593,009 controls. We detected 21 novel independent loci in the multi-ancestry analysis, including several previously associated with fibroproliferative disorders. Heritability estimates were 6%, 21%, and 34% for the European, East Asian, and African ancestry analyses, respectively. Genetically predicted gene expression and colocalization analyses identified 27 gene-tissue pairs, including nine in skin and fibroblasts. Pathway analysis implicated integrin signaling and upstream regulators involved in cancer, fibrosis, and sex hormone signaling. This investigation nearly quintuples the number of keloid-associated risk loci, illuminating biological processes in keloid pathology.

Multi-ancestry genome-wide meta-analysis with 472,819 individuals identifies 32 novel risk loci for psoriasis

BACKGROUND

Psoriasis is a common chronic, recurrent, immune-mediated disease involved in the skin or joints or both. However, deeper insight into the genetic susceptibility of psoriasis is still unclear.

METHODS

Here we performed the largest multi-ancestry meta-analysis of genome-wide association study including 28,869 psoriasis cases and 443,950 healthy controls.

RESULTS

We identified 74 genome-wide significant loci for psoriasis. Of 74 loci, 32 were novel psoriasis risk loci. Across 74 loci, 801 likely causal genes are indicated and 164 causal genes are prioritized. SNP-based heritability analyses demonstrated that common variants explain 15% of genetic risk for psoriasis. Gene-set analyses and the genetic correlation revealed that psoriasis-related genes have the positive correlations with autoimmune diseases such as ulcerative colitis, inflammatory bowel diseases, and Crohn's disease. Gene-drug interaction analysis suggested that psoriasis-associated genes overlapped with targets of current medications for psoriasis. Finally, we used the multi-ancestry meta-analysis to explore drug repurposing and the potential targets for psoriasis.

CONCLUSIONS

We identified 74 genome-wide significant loci for psoriasis. Based on 74 loci, we provided new biological insights to the etiology of psoriasis. Of clinical interest, we gave some hints for 76 potential targets and drug repurposing for psoriasis.

Exploring the role of circulating proteins in multiple myeloma risk: a Mendelian randomization study

Multiple myeloma (MM) is an incurable blood cancer with unclear aetiology. Proteomics is a valuable tool in exploring mechanisms of disease. We investigated the causal relationship between circulating proteins and MM risk, using two of the largest cohorts with proteomics data to-date. We performed bidirectional two-sample Mendelian randomization (MR; forward MR = causal effect estimation of proteins and MM risk; reverse MR = causal effect estimation of MM risk and proteins). Summary statistics for plasma proteins were obtained from genome-wide association studies performed using SomaLogic (N = 35,559; deCODE) and Olink (N = 34,557; UK Biobank; UKB) proteomic platforms and for MM risk from a meta-analysis of UKB and FinnGen (case = 1649; control = 727,247) or FinnGen only (case = 1085; control = 271,463). Cis-SNPs associated with protein levels were used to instrument circulating proteins. We evaluated proteins for the consistency of directions of effect across MR analyses (with 95% confidence intervals not overlapping the null) and corroborating evidence from genetic colocalization. In the forward MR, 994 (SomaLogic) and 1570 (Olink) proteins were instrumentable. 440 proteins were analysed in both deCODE and UKB; 302 (69%) of these showed consistent directions of effect in the forward MR. Seven proteins had 95% confidence intervals (CIs) that did not overlap the null in both forward MR analyses and did not have evidence for an effect in the reverse direction: higher levels of dermatopontin (DPT), beta-crystallin B1 (CRYBB1), interleukin-18-binding protein (IL18BP) and vascular endothelial growth factor receptor 2 (KDR) and lower levels of odorant-binding protein 2b (OBP2B), glutamate-cysteine ligase regulatory subunit (GCLM) and gamma-crystallin D (CRYGD) were implicated in increasing MM risk. Evidence from genetic colocalization did not meet our threshold for a shared causal signal between any of these proteins and MM risk (h4 < 0.8). Our results highlight seven circulating proteins which may be involved in MM risk. Although evidence from genetic colocalization suggests these associations may not be robust to the effects of horizontal pleiotropy, these proteins may be useful markers of MM risk. Future work should explore the utility of these proteins in disease prediction or prevention using proteomic data from patients with MM or precursor conditions.

Reassessing the link between adiposity and head and neck cancer: a Mendelian randomization study

Background

Adiposity has been associated with an increased risk of head and neck cancer (HNC). Although body mass index (BMI) has been inversely associated with HNC risk among smokers, this is likely due to confounding. Previous Mendelian randomization (MR) studies could not fully discount causality between adiposity and HNC due to limited statistical power. Hence, we aimed to revisit this using the largest genome-wide association study (GWAS) of HNC available, which has more granular data on HNC subsites.

Methods

We assessed the genetically predicted effects of BMI (N=806,834), waist-to-hip ratio (WHR; N=697,734) and waist circumference (N=462,166) on the risk of HNC (N=12,264 cases and 19,259 controls) and its subsites (oral, laryngeal, hypopharyngeal and oropharyngeal cancers) using a two-sample MR framework. We used the inverse variance weighted (IVW) MR approach and multiple sensitivity analyses including the weighted median, weighted mode, MR-Egger, MR-PRESSO, and CAUSE approaches. We also used multivariable MR (MVMR) to explore the direct effects of the adiposity measures on HNC, while accounting for smoking behaviour, a well-known HNC risk factor.

Results

In univariable MR, higher genetically predicted BMI increased the risk of overall HNC (IVW OR=1.17 per 1 standard deviation [1-SD] higher BMI, 95% CI 1.02-1.34, p=0.03), with no heterogeneity across subsites (Q p=0.78). However, the effect was not consistent in sensitivity analyses. The IVW effect was attenuated when smoking was included in the MVMR model (OR accounting for comprehensive smoking index=0.96 per 1-SD higher BMI, 95% CI 0.80-1.15, p=0.64) and CAUSE indicated the IVW results could be biased by correlated pleiotropy. Furthermore, we did not find a link between genetically predicted WHR (IVW OR=1.05 per 1-SD higher WHR, 95% CI 0.89-1.24, p=0.53) or waist circumference and HNC risk (IVW OR=1.01 per 1-SD higher waist circumference, 95% CI 0.85-1.21, p=0.87).

Conclusions

Our findings suggest that adiposity does not play a role in HNC risk.

Investigating the Contribution of Coding Variants in Alcohol Use Disorder Using Whole-Exome Sequencing Across Ancestries

BACKGROUND

Alcohol use disorder (AUD) is a leading cause of death and disability worldwide. There has been substantial progress in identifying genetic variants that underlie AUD. However, whole-exome sequencing studies of AUD have been hampered by the lack of available samples.

METHODS

We analyzed whole-exome sequencing data of 4530 samples from the Yale-Penn cohort and 469,835 samples from the UK Biobank, which represent an unprecedented resource for exploring the contribution of coding variants in AUD. After quality control, 1750 African-ancestry (1142 cases) and 2039 European-ancestry (1420 cases) samples from the Yale-Penn and 6142 African-ancestry (130 cases), 415,617 European-ancestry (12,861‬ cases), and 4607 South Asian (130 cases) samples from the UK Biobank cohorts were included in the analyses.

RESULTS

We confirmed the well-known functional variant rs1229984 in ADH1B (p = 4.88 × 10-31) and several other variants in ADH1C. Gene-based collapsing tests that considered the high allelic heterogeneity revealed the previously unreported genes CNST (p = 1.19 × 10-6), attributable to rare variants with allele frequency < 0.001, and IFIT5 (p = 3.74 × 10-6), driven by the burden of both common and rare loss-of-function and missense variants.

CONCLUSIONS

This study extends our understanding of the genetic architecture of AUD by providing insights into the contribution of rare coding variants, separately and convergently with common variants in AUD.

Genetic associations between orexin genes and phenotypes related to behavioral regulation in humans, including substance use

The hypothalamic neuropeptide system of orexin (hypocretin) neurons provides projections throughout the neuraxis and has been linked to sleep regulation, feeding and motivation for salient rewards including drugs of abuse. However, relatively little has been done to examine genes associated with orexin signaling and specific behavioral phenotypes in humans. Here, we tested for association of twenty-seven genes involved in orexin signaling with behavioral phenotypes in humans. We tested the full gene set, functional subsets, and individual genes involved in orexin signaling. Our primary phenotype of interest was Externalizing, a composite factor comprised of behaviors and disorders associated with reward-seeking, motivation, and behavioral regulation. We also tested for association with additional phenotypes that have been related to orexin regulation in model organism studies, including alcohol consumption, problematic alcohol use, daytime sleepiness, insomnia, cigarettes per day, smoking initiation, and body mass index. The composite set of 27 genes corresponding to orexin function was highly associated with Externalizing, as well as with alcohol consumption, insomnia, cigarettes per day, smoking initiation and BMI. In addition, all gene subsets (except the OXR2/HCRTR2 subset) were associated with Externalizing. BMI was significantly associated with all gene subsets. The "validated factors for PPOX/HCRT" and "PPOX/HCRT upregulation" gene subsets also were associated with alcohol consumption. Individually, 8 genes showed a strong association with Externalizing, 12 with BMI, 7 with smoking initiation, 3 with alcohol consumption, and 2 with problematic alcohol use, after correction for multiple testing. This study indicates that orexin genes are associated with multiple behaviors and disorders related to self-regulation in humans. This is consistent with prior work in animals that implicated orexin signaling in motivational activation induced by salient stimuli, and supports the hypothesis that orexin signaling is an important potential therapeutic target for numerous behavioral disorders.

Associations between common genetic variants and income provide insights about the socio-economic health gradient

We conducted a genome-wide association study on income among individuals of European descent (N = 668,288) to investigate the relationship between socio-economic status and health disparities. We identified 162 genomic loci associated with a common genetic factor underlying various income measures, all with small effect sizes (the Income Factor). Our polygenic index captures 1-5% of income variance, with only one fourth due to direct genetic effects. A phenome-wide association study using this index showed reduced risks for diseases including hypertension, obesity, type 2 diabetes, depression, asthma and back pain. The Income Factor had a substantial genetic correlation (0.92, s.e. = 0.006) with educational attainment. Accounting for the genetic overlap of educational attainment with income revealed that the remaining genetic signal was linked to better mental health but reduced physical health and increased risky behaviours such as drinking and smoking. These findings highlight the complex genetic influences on income and health.

Genetic insights into CRP levels in Indian adolescents: confirming adult genetic associations

CRP is a biomarker of acute inflammation linked to metabolic complications. Given the rising prevalence of these conditions in India, we investigated the genetic basis of CRP levels in Indian adolescents, an underrepresented group in genetic studies, to identify early markers of metabolic risk. We performed a two-phased genome-wide association study (GWAS; N = 5052) and an independent Exome-wide association study (ExWAS; N = 4547), to identify both common and rare genetic variants associated with CRP levels. The study identified intergenic variants near CRP and CRPP1 genes, and APOC1 gene as the key regulators of CRP levels establishing the universality of these associations. The GWAS identified the variant rs4247360 (PITPNC1) to be associated at a suggestive significance. The ExWAS single variant association identified novel associations in genes FGL1 (rs35431851), C19orf45 (rs608144, rs475923, rs484870), TRAPPC12 (rs11686212) and KIAA0087 (rs17153822). The SKATO analysis of the rare variants highlighted the role of loss of function and missense variants in genes EPS15, CCDC15, ZNF286A, ELF1, B3GNT8, ZNF850, MAP2, and PSG2. The GWAS and ExWAS in the present study validated the association of 56 variants previously reported for CRP levels. The meta-analysis with the CRP GWAS earlier reported in Indian adults revealed the shared genetic architecture of CRP levels across age groups. The gene-set enrichment analysis highlighted the role of CRP-associated genes in inflammatory and cardiometabolic pathways. The study enhances understanding of genetic predispositions to inflammation and metabolic disorders confirming known associations, identifying novel loci, and validating shared genetic architecture across age-groups, guiding targeted prevention for at-risk youth.

Multi-ancestry genome-wide association analyses: a comparison of meta- and mega-analyses in the Hyperglycemia and Adverse Pregnancy Outcome (HAPO) study

BACKGROUND

There is increasing need for effective incorporation of high-dimensional genetics data from individuals with varied ancestry in genome-wide association (GWAS) analyses. Classically, multi-ancestry GWAS analyses are performed using statistical meta-analysis to combine results conducted within homogeneous ancestry groups. The emergence of cosmopolitan reference panels makes collective preprocessing of GWAS data possible, but impact on downstream GWAS results in a mega-analysis framework merits investigation. We utilized GWAS data from the multi-national Hyperglycemia and Adverse Pregnancy Outcome Study to investigate differences in GWAS findings using a homogeneous ancestry meta-analysis versus a heterogeneous ancestry mega-analysis pipeline. Maternal fasting and 1-hr glucose and metabolomics measured during a 2-hr 75-gram oral glucose tolerance test during early third trimester pregnancy were evaluated as phenotypes.

RESULTS

For the homogeneous ancestry meta-analysis pipeline, variant data were prepared by identifying sets of individuals with similar ancestry and imputing to ancestry-specific reference panels. GWAS was conducted within each ancestry group and results were combined using random-effects meta-analysis. For the heterogeneous ancestry mega-analysis pipeline, data for all individuals were collectively imputed to the Trans-Omics for Precision Medicine (TOPMed) cosmopolitan reference panel, and GWAS was conducted using a unified mega-analysis. The meta-analysis pipeline identified genome-wide significant associations for 15 variants in a region close to GCK on chromosome 7 with maternal fasting glucose and no significant findings for 1-hr glucose. Associations in this same region were identified using the mega-analysis pipeline, along with a well-documented association at MTNR1B on chromosome 11 with both fasting and 1-hr maternal glucose. For metabolomics analyses, the number of significant findings in the heterogeneous ancestry mega-analysis far exceeded those from the homogeneous ancestry meta-analysis and confirmed many previously documented associations, but genomic inflation factors were much more variable.

CONCLUSIONS

For multi-ancestry GWAS, heterogeneous ancestry mega-analysis generates a rich set of variants for analysis using a cosmopolitan reference panel and results in vastly more significant, biologically credible and previously documented associations than a homogeneous ancestry meta-analysis approach. Genomic inflation factors do indicate that findings from the mega-analysis pipeline may merit cautious interpretation and further follow-up.

Characterization of gene-environment interactions for vitamin D through variance quantitative trait loci: a UK Biobank-based genetic epidemiology study

BACKGROUND

Understanding gene-environment interactions associated with vitamin D status may refine nutrition and public health strategies for vitamin D deficiency. Recent methodological advances have enabled the identification of variance quantitative trait loci (vQTLs) where gene-environment interactions are enriched.

OBJECTIVES

The study aims to identify vQTLs for serum 25-hydroxy vitamin D (25OHD) concentrations and characterize potential gene-environment interactions of vQTLs.

METHODS

We conducted vQTL discovery for 25OHD using a newly developed quantile integral linear model in the UK Biobank individuals of European (N = 313,514), African (N = 7800), East Asian (N = 2146), and South Asian (N = 8771) ancestries, respectively. We tested for interactions between the identified vQTL lead variants and 18 environmental, biological, or lifestyle factors, followed by multiple sensitivity analyses.

RESULTS

We identified 19 independent vQTL lead variants (P < 5 × 10-8) in the European ancestry population. No vQTLs were identified in the non-European ancestry populations, likely because of limited sample sizes. A total of 32 interactions were detected with a false discovery rate <0.05. Although known gene-season of measurement interactions were confirmed, additional interactions were identified involving modifiable risk factors, including time spent outdoors and body mass index. The magnitudes of these interactions were consistent within each locus upon adjusting for the season of measurement and other covariates. We also identified a gene-sex interaction at a vQTL that implicates DHCR7. Integrating transcript- and protein-level evidence, we found that the sex-differentiated genetic associations may act through sex-biased expression of DHCR7 isoforms in skin tissues because of alternative splicing.

CONCLUSIONS

Through the lens of vQTLs, we identified additional gene-environment interactions affecting vitamin D status in addition to the season of measurement. These findings may provide new insights into the etiology of vitamin D deficiency and encourage personalized prevention and management of associated diseases for at-risk individuals.

Associations of maternal night shift work during pregnancy with DNA methylation in offspring: a meta-analysis in the PACE consortium

BACKGROUND

Night shift work during pregnancy has been associated with differential DNA methylation in placental tissue, but no studies have explored this association in cord blood. We aimed to examine associations of maternal night shift work with cord blood DNA methylation.

METHODS

A total of 4487 mother-newborn pairs from 7 studies were included. Maternal night shift work during pregnancy was ascertained via questionnaires and harmonized into "any" versus "no". DNA methylation was measured in cord blood using the Illumina Infinium Methylation arrays. Robust linear regression models adjusted for relevant confounders were run in the individual cohorts, and results were meta-analyzed.

RESULTS

Maternal night shift work during pregnancy ranged from 3.4% to 26.3%. Three CpGs were differentially methylated in relation to maternal night shift work during pregnancy at a false discovery rate adjusted P < 0.05: cg10945885 (estimate (β) 0.38%, standard error (SE) 0.07), cg00773359 (β 0.25%, SE 0.05), and cg21836426 (β - 0.29%, SE 0.05). Associations of the identified CpGs were found in previous literature for gestational age and childhood and adolescent BMI. In a mouse model of prenatal jet lag exposure, information on offspring DNA methylation of ten homologous genes annotated to the 16 CpGs with P < 1 × 10-5 in our analysis was available, of which eight were associated (enrichment P: 1.62 × 10-11).

CONCLUSION

Maternal night shift work during pregnancy was associated with newborn DNA methylation at 3 CpGs. Top findings overlapped with those in a mouse model of gestational jet lag. This work strengthens evidence that DNA methylation could be a marker or mediator of impacts of circadian rhythm disturbances.

Epigenome-wide association study of objectively measured physical activity in peripheral blood leukocytes

BACKGROUND

Few studies have explored the association between DNA methylation and physical activity. The aim of this study was to evaluate the association of objectively measured hours of sedentary behavior (SB) and moderate physical activity (MPA) with DNA methylation. We further aimed to explore the association between SB or MPA related CpG sites and cardiometabolic traits, gene expression, and genetic variation.

RESULTS

For discovery, we performed cross sectional analyses in pregnant women from the Epigenetics in pregnancy (EPIPREG) sample with both DNA methylation (Illumina MethylationEPIC BeadChip) and objectively measured physical activity data (SenseWear™ Pro 3 armband) (European = 244, South Asian = 109). For EWAS of SB and MPA, two main models were designed: model (1) a linear mixed model adjusted for age, smoking, blood cell composition, including ancestry as random intercept, and model (2) which was additionally adjusted for the total number of steps per day. In model 1, we did not identify any CpG sites associated with neither SB nor MPA. In model 2, SB was positively associated (false discovery rate, FDR < 0.05) with two CpG sites within the VSX1 gene. Both CpG sites were positively associated with BMI and were associated with several genetic variants in cis. MPA was associated with 122 significant CpG sites at FDR < 0.05 (model 2). We further analyzed the ten most statistically significant MPA related CpG sites and found that they presented opposite associations with sedentary behavior and BMI. We were not able to replicate the SB and MPA-related CpG sites in the Avon Longitudinal Study of Parents and Children (ALSPAC). ALSPAC had available objectively measured physical activity data from Actigraph (without steps/day available) and leucocyte DNA methylation data collected during adolescence (n = 408, European).

CONCLUSION

This study suggests associations of objectively measured SB and MPA with maternal DNA methylation in peripheral blood leukocytes, that needs to be confirmed in larger samples of similar study design.

Matching Heterogeneous Cohorts by Projected Principal Components Reveals Two Novel Alzheimer's Disease-Associated Genes in the Hispanic Population

Alzheimer's disease (AD) is the most common form of dementia in elderly, affecting 6.9 million individuals in the United States. Some studies have suggested the prevalence of AD is greater in individuals who self-identify as Hispanic. Focused results are relevant for personalized and equitable clinical interventions. Ethnicity as a stratifying tool in genetic studies is often accompanied by genomic inflation due to heterogeneity. In this study, we report GWAS and meta-analyses conducted among NIAGADS subjects who self-identified as Hispanic and All of Us (AoU) sub-cohorts matched to that cohort, using projected genetically-derived principal components, with and without age and sex. In Hispanic NIAGADS subjects, we identified a common variant in PIEZO2 that was protective for AD with a p-value just beyond genome-wide significance (p = 5.4 * 10-8). Meta-analyses with genetically-matched AoU participants yielded three (two novel) genome-wide significant AD-associated loci based on rare lead variants: rs374043832 (RGS6/PSEN1), rs192423465 (ASPSCR1), and rs935208076 (GDAP2), which were also nominally significant in AoU sub-cohorts. We also show how genomic inflation can be mitigated in heterogeneous populations while increasing sample size and result generalizability.

Integrative genetics and multiomics analysis reveal mechanisms and therapeutic targets in vitiligo highlighting JAK STAT pathway regulation of CTSS

Vitiligo is a complex autoimmune disease characterized by the loss of melanocytes, leading to skin depigmentation. Despite advances in understanding its genetic and molecular basis, the precise mechanisms driving vitiligo remain elusive. Integrating multiple layers of omics data can provide a comprehensive view of disease pathogenesis and identify potential therapeutic targets. The study aims to delineate the genetic and molecular mechanisms of vitiligo pathogenesis using an integrative multiomics strategy. We focus on exploring the regulatory influence of the JAK/STAT pathway on Cathepsin S, a potential therapeutic target in vitiligo. Our GWAS-meta analysis pinpointed five druggable genes: ERBB3, RHOH, CDK10, MC1R, and NDUFAF3, and underwent drug target exploration and molecular docking. SMR analysis linked CTSS, CTSH, STX8, KIR2DL3, and GRHPR to vitiligo through pQTL and eQTL associations. Microarray and single-cell RNA-seq data showed differential expression of CTSS and STAT1/3 in vitiligo patients' blood and skin lesions. Our study offers novel perspectives on vitiligo's genetic and molecular basis, highlighting the JAK/STAT pathway's role in regulating CTSS for antigen processing in melanocytes. Further research is needed to confirm these results and assess the therapeutic potential of CTSS and related genes.

Polymorphisms within autophagy-related genes as susceptibility biomarkers for pancreatic cancer: A meta-analysis of three large European cohorts and functional characterization

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers with patients having unresectable or metastatic disease at diagnosis, with poor prognosis and very short survival. Given that genetic variation within autophagy-related genes influences autophagic flux and susceptibility to solid cancers, we decided to investigate whether 55,583 single nucleotide polymorphisms (SNPs) within 234 autophagy-related genes could influence the risk of developing PDAC in three large independent cohorts of European ancestry including 12,754 PDAC cases and 324,926 controls. The meta-analysis of these populations identified, for the first time, the association of the BIDrs9604789 variant with an increased risk of developing the disease (ORMeta = 1.31, p = 9.67 × 10-6). We also confirmed the association of TP63rs1515496 and TP63rs35389543 variants with PDAC risk (OR = 0.89, p = 6.27 × 10-8 and OR = 1.16, p = 2.74 × 10-5). Although it is known that BID induces autophagy and TP63 promotes cell growth, cell motility and invasion, we also found that carriers of the TP63rs1515496G allele had increased numbers of FOXP3+ Helios+ T regulatory cells and CD45RA+ T regulatory cells (p = 7.67 × 10-4 and p = 1.56 × 10-3), but also decreased levels of CD4+ T regulatory cells (p = 7.86 × 10-4). These results were in agreement with research suggesting that the TP63rs1515496 variant alters binding sites for FOXA1 and CTCF, which are transcription factors involved in modulating specific subsets of regulatory T cells. In conclusion, this study identifies BID as new susceptibility locus for PDAC and confirms previous studies suggesting that the TP63 gene is involved in the development of PDAC. This study also suggests new pathogenic mechanisms of the TP63 locus in PDAC.

The impact of Indigenous American-like ancestry on risk of acute lymphoblastic leukemia in Hispanic/Latino children

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer, with Hispanic/Latino children having a higher incidence of ALL than other racial/ethnic groups. Genetic variants, particularly ones found enriched in Indigenous American (IA)-like ancestry and inherited by Hispanics/Latinos, may contribute to this disparity. In this study, we characterized the impact of IA-like ancestry on overall ALL risk and the frequency and effect size of known risk alleles in a large cohort of self-reported Hispanic/Latino individuals. We also performed genome-wide admixture mapping analysis to identify potentially novel ALL risk loci. We found that global IA ancestry was positively associated with ALL risk, but the association was not significant after adjusting for socio-economic indicators. In a series of local ancestry analyses, we uncovered that at known ALL risk loci, increasing copies of the IA-like haplotype were positively and significantly associated with ALL case-control status. Further, the IA-like haplotype had ~1.33 times the odds of harboring the risk allele compared to non-IA-like haplotypes. We found no evidence of interaction between genotype and ancestry (local or global) in relation to ALL risk. Admixture mapping identified association signals on chromosomes 2 (2q21.2), 7 (7p12.2), 10 (10q21.2), and 15 (15q22.31); however, only the variants at 7p12.2 and 10q21.2 replicated in additional cohorts. Taken together, our results suggest that increased risk of ALL in Hispanic/Latino children may be conferred by higher frequency of risk alleles within IA-like ancestry, which can be leveraged as targets of new precision health strategies and therapeutics.

Effects of Genetic Risk and Lifestyle Habits on Gout: A Korean Cohort Study

BACKGROUND

Gout is a type of inflammatory arthritis caused by monosodium urate crystal deposits, and the prevalence of this condition has been increasing. This study aimed to determine the combined effects of genetic risk factors and lifestyle habits on gout, using data from a Korean cohort study. Identifying high-risk individuals in advance can help prevent gout and its associated disorders.

METHODS

We analyzed data from the Korean Genome and Epidemiology Study-Urban Health Examinees cohort (KoGES-HEXA). Genetic information of the participants was collected at baseline, and gout cases were identified based on patient statements. The polygenic risk score (PRS) was calculated using nine independent genome-wide association study datasets, and lifestyle factors and metabolic syndrome status were measured for each participant using the KoGES. Logistic regression models were used to estimate the odds ratios (ORs) for gout in relation to genetic risk, lifestyle habits, and metabolic health status, after adjusting for age and sex.

RESULTS

Among 44,605 participants, 617 were diagnosed with gout. Gout was associated with older age, higher body mass index, and higher prevalence of hypertension, diabetes, and hypertriglyceridemia. High PRS, unfavorable lifestyle habits, and poor metabolic profiles were significantly associated with an increased risk of gout. Compared with that in the low-genetic-risk and healthy lifestyle group or ideal metabolic profile group, the risk of gout was increased in the high-genetic-risk plus unfavorable lifestyle (OR, 3.64; 95% confidence interval [CI], 2.32-6.03) or poor metabolic profile (OR, 7.78; 95% CI, 4.61-13.40) group. Conversely, adherence to favorable lifestyle habits significantly reduced gout risk, especially in high-genetic-risk groups.

CONCLUSION

Genetic predisposition and unhealthy lifestyle habits significantly increase the risk of gout. Promoting healthy lifestyle habits is crucial to prevent the development of gout, particularly in individuals with high genetic susceptibility.

The genomic architecture of circulating cytokine levels points to drug targets for immune-related diseases

Circulating cytokines orchestrate immune reactions and are promising drug targets for immune-mediated and inflammatory diseases. Exploring the genetic architecture of circulating cytokine levels could yield key insights into causal mediators of human disease. Here, we performed genome-wide association studies (GWAS) for 40 circulating cytokines in meta-analyses of 74,783 individuals. We detected 359 significant associations between cytokine levels and variants in 169 independent loci, including 150 trans- and 19 cis-acting loci. Integration with transcriptomic data point to key regulatory mechanisms, such as the buffering function of the Atypical Chemokine Receptor 1 (ACKR1) acting as scavenger for multiple chemokines and the role of tumor necrosis factor receptor-associated factor 1 (TRAFD1) in modulating the cytokine storm triggered by TNF signaling. Applying Mendelian randomization (MR), we detected a network of complex cytokine interconnections with TNF-b, VEGF, and IL-1ra exhibiting pleiotropic downstream effects on multiple cytokines. Drug target cis-MR using 2 independent proteomics datasets paired with colocalization revealed G-CSF/CSF-3 and CXCL9/MIG as potential causal mediators of asthma and Crohn's disease, respectively, but also a potentially protective role of TNF-b in multiple sclerosis. Our results provide an overview of the genetic architecture of circulating cytokines and could guide the development of targeted immunotherapies.

Patient stratification by genetic risk in Alzheimer's disease is only effective in the presence of phenotypic heterogeneity

Case-only designs in longitudinal cohorts are a valuable resource for identifying disease-relevant genes, pathways, and novel targets influencing disease progression. This is particularly relevant in Alzheimer's disease (AD), where longitudinal cohorts measure disease "progression," defined by rate of cognitive decline. Few of the identified drug targets for AD have been clinically tractable, and phenotypic heterogeneity is an obstacle to both clinical research and basic science. In four cohorts (n = 7241), we performed genome-wide association studies (GWAS) and Mendelian randomization (MR) to discover novel targets associated with progression and assess causal relationships. We tested opportunities for patient stratification by deriving polygenic risk scores (PRS) for AD risk and severity and tested the value of these scores in predicting progression. Genome-wide association studies identified no loci associated with progression at genome-wide significance (α = 5×10-8); MR analyses provided no significant evidence of an association between cognitive decline in AD patients and protein levels in brain, cerebrospinal fluid (CSF), and plasma. Polygenic risk scores for AD risk did not reliably stratify fast from slow progressors; however, a deeper investigation found that APOE ε4 status predicts amyloid-β and tau positive versus negative patients (odds ratio for an additional APOE ε4 allele = 5.78 [95% confidence interval: 3.76-8.89], P<0.001) when restricting to a subset of patients with available CSF biomarker data. These results provided no evidence for large-effect, common-variant loci involved in the rate of memory decline, suggesting that patient stratification based on common genetic risk factors for progression may have limited utility. Where clinically relevant biomarkers suggest diagnostic heterogeneity, there is evidence that a priori identified genetic risk factors may have value in patient stratification. Mendelian randomization was less tractable due to the lack of large-effect loci, and future analyses with increased samples sizes are needed to replicate and validate our results.

A multi-ancestry genome-wide association study identifies novel candidate loci in the RARB gene associated with hypertensive disorders of pregnancy

Genetic factors related to pregnancy-related traits are understudied, especially in ancestrally diverse cohorts. To assess maternal contributions to hypertensive disorders of pregnancy (HDP), we performed a multi-ancestry genome-wide association study (GWAS) of HDP in data from the North Carolina-based Personalized Environment and Genes Study (PEGS) cohort with validation in the UK Biobank (UKBB). The GWAS revealed two maternal loci associated with HDP at the genome-wide significance level. The lead independent variants were rs114954125 on chromosome 2 (near LRP1B; odds ratio [OR] [95% confidence interval {CI}]): 2.96 [2.02-4.34]; p = 2.82 × 10-8) and rs61176331 on chromosome 3 (on RARB; OR (95% CI): 3.08 (2.12-4.48); p = 3.52 × 10-9). We validated the associations near RARB with a meta-analysis of PEGS and the UKBB. We also identified cis-expression quantitative trait loci in the candidate region associated with decreased RARB expression in macrophage cells exposed to Salmonella. Chromatin mapping in FUMA identified a significant interaction within chromosome 3's enhancer and open chromatin regions, with strong effects observed for RARB and H3P10 gene regulation in mesendoderm cells, mesenchymal stem cells, and trophoblast-like stem cells. We applied existing polygenic scores (PGS) for preeclampsia and gestational hypertension and found that the scores were significantly associated with HDP in PEGS. The findings demonstrate the power of multi-ancestry studies for genetic discovery and highlight the relationship between immune response, regulation, and HDP and the utility of PGS for risk prediction. (PEGS is registered at ClinicalTrials.gov: NCT00341237.).

Soluble Immune Checkpoint Protein and Lipid Network Associations with All-Cause Mortality Risk: Trans-Omics for Precision Medicine (TOPMed) Program

Adverse cardiovascular events are emerging with the use of immune checkpoint therapies in oncology. Using datasets in the Trans-Omics for Precision Medicine program (Multi-Ethnic Study of Atherosclerosis, Jackson Heart Study [JHS], and Framingham Heart Study), we examined the association of immune checkpoint plasma proteins with each other, their associated protein network with high-density lipoprotein cholesterol (HDL-C) and low-density lipoprotein cholesterol (LDL-C), and the association of HDL-C- and LDL-C-associated protein networks with all-cause mortality risk. Plasma levels of LAG3 and HAVCR2 showed statistically significant associations with mortality risk. Colocalization analysis using genome wide-association studies of HDL-C or LDL-C and protein quantitative trait loci from JHS and the Atherosclerosis Risk in Communities identified TFF3 rs60467699 and CD36 rs3211938 variants as significantly colocalized with HDL-C; in contrast, none colocalized with LDL-C. The measurement of plasma LAG3, HAVCR2, and associated proteins plus targeted genotyping may identify patients at increased mortality risk.

Genome-wide association study reveals shared and distinct genetic architecture of fatty acids and oxylipins in the Hispanic Community Health Study/Study of Latinos

Bioactive fatty acid-derived oxylipin molecules play key roles mediating inflammation and oxidative stress. Circulating levels of fatty acids and oxylipins are influenced by environmental and genetic factors; characterizing the genetic architecture of bioactive lipids could yield new insights into underlying biology. We performed a genome-wide association study (GWAS) of 81 fatty acids and oxylipins in 11,584 Hispanic Community Health Study/Study of Latinos (HCHS/SOL) participants with genetic and lipidomic data measured at study baseline (58.6% female, mean age = 46.1 years (standard deviation 13.8)). Additionally, given the effects of central obesity on inflammation, we examined interactions with waist circumference using two-degree-of-freedom joint tests. Thirty-three of the 81 oxylipins and fatty acids were significantly heritable (heritability range: 0-32.7%). Forty (49.4%) oxylipins and fatty acids had at least one genome-wide significant (p < 6.94E-11) variant resulting in 19 independent genetic loci. Six loci (lead variant minor allele frequency [MAF] range: 0.08-0.50), including desaturase-encoding FADS and OATP1B1 transporter protein-encoding SLCO1B1, exhibited associations with two or more fatty acids and oxylipins. At several of these loci, there was evidence of colocalization of the top variant across fatty acids and oxylipins. The remaining loci were only associated with one oxylipin or fatty acid and included several CYP loci. We also identified an additional rare variant (MAF = 0.002) near CARS2 in two-degree-of-freedom tests. Our analyses revealed shared and distinct genetic architecture underlying fatty acids and oxylipins, providing insights into genetic factors and motivating work to characterize these compounds and elucidate their roles in disease.

Genome-wide association studies in a large Korean cohort identify novel quantitative trait loci for 36 traits and illuminate their genetic architectures

Genome-wide association studies (GWAS) have been predominantly conducted in populations of European ancestry, limiting opportunities for biological discovery in diverse populations. We report GWAS findings from 153,950 individuals across 36 quantitative traits in the Korean Cancer Prevention Study-II (KCPS2) Biobank. We discovered 301 novel genetic loci in KCPS2, including an association between thyroid-stimulating hormone and CD36. Meta-analysis with the Korean Genome and Epidemiology Study, Biobank Japan, Taiwan Biobank, and UK Biobank identified 4,588 loci that were not significant in any contributing GWAS. We describe differences in genetic architectures across these East Asian and European samples. We also highlight East Asian specific associations, including a known pleiotropic missense variant in ALDH2, which fine-mapping identified as a likely causal variant for a diverse set of traits. Our findings provide insights into the genetic architecture of complex traits in East Asian populations and highlight how broadening the population diversity of GWAS samples can aid discovery.

Exploring the Role of Glycine Metabolism in Coronary Artery Disease: Insights from Human Genetics and Mouse Models

Background: Circulating glycine levels have been associated with reduced risk of coronary artery disease (CAD) in humans but these associations have not been observed in all studies. We evaluated whether the relationship between glycine levels and atherosclerosis was causal using genetic analyses in humans and feeding studies in mice. Methods: Serum glycine levels were evaluated for association with risk of CAD in the UK Biobank. Genetic determinants of glycine levels were identified through a genome-wide association study (GWAS) and used to evaluate the causal relationship between glycine and risk of CAD by Mendelian randomization (MR). A dietary supplementation study was carried out with atherosclerosis-prone apolipoprotein E deficient (ApoE-/-) mice to determine the effects of increased circulating glycine levels on cardiometabolic traits and aortic lesion formation. Results: Among 105,718 UK Biobank subjects, elevated serum glycine levels were associated with significantly reduced risk of prevalent CAD (Quintile 5 vs. Quintile 1 OR = 0.76, 95% CI 0.67-0.87; p < 0.0001) and incident CAD (Quintile 5 vs. Quintile 1 HR = 0.70, 95% CI 0.65-0.77; p < 0.0001) after adjustment for age, sex, ethnicity, anti-hypertensive and lipid-lowering medications, blood pressure, kidney function, and diabetes. A GWAS meta-analysis with 230,947 subjects identified 61 loci for glycine levels, of which 26 were novel. MR analyses provided modest evidence that genetically elevated glycine levels were causally associated with reduced systolic blood pressure and risk of type 2 diabetes, but did not provide significant evidence for an association with decreased risk of CAD. Glycine supplementation in mice had no effects on cardiometabolic traits or atherosclerotic lesion development. Conclusions: While expanding the genetic architecture of glycine metabolism, MR analyses and in vivo feeding studies did not provide evidence that the clinical association of this amino acid with atherosclerosis represents a causal relationship.

Identification of Endometriosis Pathophysiologic-Related Genes Based on Meta-Analysis and Bayesian Approach

Endometriosis is a complex disease with diverse etiologies, including hormonal, immunological, and environmental factors; however, its exact pathogenesis remains unknown. While surgical approaches are the diagnostic and therapeutic gold standard, identifying endometriosis-associated genes is a crucial first step. Five endometriosis-related gene expression studies were selected from the available datasets. Approximately, 14,167 genes common to these 5 datasets were analyzed for differential expression. Meta-analyses utilized fold-change values and standard errors obtained from each analysis, with the binomial and continuous datasets contributing to endometriosis presence and endometriosis severity meta-analysis, respectively. Approximately 160 genes showed significant results in both meta-analyses. For Bayesian analysis, endometriosis-related single nucleotide polymorphisms (SNPs), the human transcription factor catalog, uterine SNP-related gene expression, disease-gene databases, and interactome databases were utilized. Twenty-four genes, present in at least three or more databases, were identified. Network analysis based on Pearson's correlation coefficients revealed the HLA-DQB1 gene with both a high score in the Bayesian analysis and a central position in the network. Although ZNF24 had a lower score, it occupied a central position in the network, followed by other ZNF family members. Bayesian analysis identified genes with high confidence that could support discovering key diagnostic biomarkers and therapeutic targets for endometriosis.

GWAS highlights the neuronal contribution to multiple sclerosis susceptibility

Multiple Sclerosis (MS) is a chronic inflammatory and neurodegenerative disease affecting the brain and spinal cord. Genetic studies have identified many risk loci, that were thought to primarily impact immune cells and microglia. Here, we performed a multi-ancestry genome-wide association study with 20,831 MS and 729,220 control participants, identifying 236 susceptibility variants outside the Major Histocompatibility Complex, including four novel loci. We derived a polygenic score for MS and, optimized for European ancestry, it is informative for African-American and Latino participants. Integrating single-cell data from blood and brain tissue, we identified 76 genes affected by MS risk variants. Notably, while T cells showed the strongest enrichment, inhibitory neurons emerged as a key cell type. The expression of IL7 and STAT3 are affected only in inhibitory neurons, highlighting the importance of neuronal and glial dysfunction in MS susceptibility.

A meta-analysis of genome-wide association studies to identify candidate genes associated with feed efficiency traits in pigs

Pig production is an agricultural sector of great economic and social relevance to Brazil and global markets. Feed efficiency traits directly influence the sustainability of pig production due to the economic impact of feed costs on the production system and the environmental footprint of the industry. Therefore, breeding for improved feed efficiency has been a target of worldwide pig breeding programs. Genome-wide association studies (GWAS) enable the assessment of the genetic background of complex traits, which contributes to a better understanding of the biological mechanisms regulating their phenotypic expression. In this context, the primary objective of this study was to identify and validate genomic regions and candidate genes associated with feed conversion ratio (FCR) and residual feed intake (RFI) in pigs based on a comprehensive systematic review and meta-analysis of GWAS. The METAL software was used to implement the meta-analysis and the Bonferroni multiple testing correction considering a significance threshold 0.05. The significant single nucleotide polymorphisms (SNPs) in the meta-analysis were used to identify candidate genes, followed by a functional genomic enrichment analysis. The systematic review identified 13 studies, of which 7 evaluated FCR, 3 evaluated RFI, and 3 studies investigated both traits, with 160 and 96 SNPs identified for FCR and RFI, respectively. After the meta-analysis, 145 markers were significantly associated with FCR and 90 with RFI. The gene annotation process resulted in 105 and 114 genes for FCR and RFI, respectively. The enrichment analysis for FCR resulted in 16 significant gene ontology (GO) terms, while 6 terms were identified for RFI. The main GO terms were actin cytoskeleton (GO_BP:0030036), membrane (GO_CC:0016020), integral components of the peroxisomal membrane (GO_CC:0005779), and carbohydrate-binding (GO_MF:0030246). The main candidate genes identified were MED18, PHACTR4, ABCC2, TRHDE, FRS2, FAR2 and FIS1 for FCR, and ADGRL2, ASGR1, ASGR2, and MAN2B1 for RFI. These findings contribute to a better understanding of the genetic mechanisms associated with feed efficiency traits in pigs, providing a foundation for future improvements in pig breeding programs.

Neuropathology-based approach reveals novel Alzheimer's Disease genes and highlights female-specific pathways and causal links to disrupted lipid metabolism: insights into a vicious cycle

Dementia refers to an umbrella phenotype of many different underlying pathologies with Alzheimer's disease (AD) being the most common type. Neuropathological examination remains the gold standard for accurate AD diagnosis, however, most that we know about AD genetics is based on Genome-Wide Association Studies (GWAS) of clinically defined AD. Such studies have identified multiple AD susceptibility variants with a significant portion of the heritability unexplained and highlighting the phenotypic and genetic heterogeneity of the clinically defined entity. Furthermore, despite women's increased susceptibility to dementia, there is a lack of sex-specific genetic studies and understanding of sex-specific background for the disorder. Here, we aim to tackle the heterogeneity of AD by specifically concentrating on neuropathological features and pursuing sex-specific analysis. We bring together 14 different genomic and neuropathology datasets (6960 individuals) and we integrate our GWAS findings with transcriptomic and phenotypic data aiming to also identify biomarkers for AD progression. We uncover novel genetic associations to AD neuropathology, including BIN1 and OPCML. Our sex-specific analysis points to a role for BIN1 specifically in women as well as novel AD loci including QRFPR and SGCZ. Post-GWAS analyses illuminate the functional and biological mechanisms underlying AD and reveal sex-specific differences. Finally, through PheWAS and Mendelian Randomization analysis, we identify causal links with AD neuropathology pointing to disrupted lipid metabolism, as well as impaired peripheral immune response and liver dysfunction as part of a vicious cycle that fuels neurodegeneration.

A genome-wide association study of high-sensitivity C-reactive protein in a large Korean population highlights its genetic relationship with cholesterol metabolism

High-sensitivity C-reactive protein (hsCRP) is a representative biomarker of systemic inflammation and is associated with numerous chronic diseases. To explore the biological pathways and functions underlying chronic inflammation, we conducted a genome-wide association study (GWAS) and several post-GWAS analyses of the hsCRP levels. This study was performed on data from 71,019 Koreans and is one of the largest East Asian studies. Overall, 69 independent single nucleotide polymorphisms (SNPs) were identified, including 13 novel variants. The implicated genes and pathways are primarily involved in cholesterol metabolism and the immune response. A phenome-wide association study was performed based on a polygenic risk score (PRS) constructed using 69 hsCRP-associated SNPs. Notably, the alleles associated with higher hsCRP levels appeared to be associated with lower low-density lipoprotein cholesterol levels (P = 1.69 × 10-33, β = -1.47) and higher γ -glutamyl transpeptidase (P = 8.30 × 10-8, β = 0.84). It suggests that increase in genetically determined hsCRP may contribute to a decrease in cholesterol level and a strong oxidative environment in the blood vessel. Thus, individuals with higher hsCRP-PRS may have a greater risk of cardiovascular diseases. Our findings suggest the genetic association between cholesterol and hsCRP, as well as the clinical importance of hsCRP-PRS for predicting the potential risk of cardiovascular diseases.

Unraveling the causal impact of smoking and its DNA methylation signatures on cardiovascular disease: Mendelian randomization and colocalization analysis

BACKGROUND

To explore the mechanisms linking smoking to cardiovascular diseases (CVDs) from an epigenetic perspective.

METHODS

Mendelian Randomization (MR) analysis was performed to assess the causal effects of smoking behavior and DNA methylation levels at smoking-related CpG sites on nine CVDs, including aortic aneurysm, atrial fibrillation, coronary atherosclerosis, coronary heart disease, heart failure, intracerebral hemorrhage, ischemic stroke, myocardial infarction, subarachnoid hemorrhage. Colocalization analysis was used to further identify key smoking-related CpG sites from the MR causal estimates. Reactome enrichment analysis was used to elucidate the potential mechanisms.

RESULTS

MR analysis indicates that smoking behaviors are significantly associated with an increased risk of nine CVDs (OR > 1, P < 0.05). Through MR and colocalization analysis, five key smoking-related CpG sites were ultimately determined. DNA methylation alteration at cg25313468 (located in the TSS1500 region of REST) is simultaneously associated with the risk of atrial fibrillation, coronary atherosclerosis, coronary heart disease, and myocardial infarction. Additionally, cg21647257 (located in the TSS200 region of CLIP3) is associated with the risk of atrial fibrillation; cg06197751 (located in SGEF gene body) and cg07520810 (located in ARID5B gene body) are associated with the risk of coronary atherosclerosis; cg16822035 (located in MCF2L gene body) is associated with the risk of myocardial infarction. Enrichment analysis suggests that phosphatase and tensin homologue (PTEN) may be involved in the downstream mechanisms of cg25313468 (REST).

CONCLUSION

This study uncovers the relationship between smoking, DNA methylation, and CVDs, providing new insights into the pathogenic effect of smoking on CVDs from an epigenetic perspective.

Genome-wide characterization of 54 urinary metabolites reveals molecular impact of kidney function

Dissecting the genetic mechanisms underlying urinary metabolite concentrations can provide molecular insights into kidney function and open possibilities for causal assessment of urinary metabolites with risk factors and disease outcomes. Proton nuclear magnetic resonance metabolomics provides a high-throughput means for urinary metabolite profiling, as widely applied for blood biomarker studies. Here we report a genome-wide association study meta-analysed for 3 European cohorts comprising 8,011 individuals, covering both people with type 1 diabetes and general population settings. We identify 54 associations (p < 9.3 × 10-10) for 19 of 54 studied metabolite concentrations. Out of these, 33 were not reported previously for relevant urinary or blood metabolite traits. Subsequent two-sample Mendelian randomization analysis suggests that estimated glomerular filtration rate causally affects 13 urinary metabolite concentrations whereas urinary ethanolamine, an initial precursor for phosphatidylcholine and phosphatidylethanolamine, was associated with higher eGFR lending support for a potential protective role. Our study provides a catalogue of genetic associations for 53 metabolites, enabling further investigation on how urinary metabolites are linked to human health.

Novel Genetic Loci in Early-Onset Gout Derived From Whole-Genome Sequencing of an Adolescent Gout Cohort

OBJECTIVE

Mechanisms underlying the adolescent-onset and early-onset gout are unclear. This study aimed to discover variants associated with early-onset gout.

METHODS

We conducted whole-genome sequencing in a discovery adolescent-onset gout cohort of 905 individuals (gout onset 12 to 19 years) to discover common and low-frequency single-nucleotide variants (SNVs) associated with gout. Candidate common SNVs were genotyped in an early-onset gout cohort of 2,834 individuals (gout onset ≤30 years old), and meta-analysis was performed with the discovery and replication cohorts to identify loci associated with early-onset gout. Transcriptome and epigenomic analyses, quantitative real-time polymerase chain reaction and RNA sequencing in human peripheral blood leukocytes, and knock-down experiments in human THP-1 macrophage cells investigated the regulation and function of candidate gene RCOR1.

RESULTS

In addition to ABCG2, a urate transporter previously linked to pediatric-onset and early-onset gout, we identified two novel loci (Pmeta < 5.0 × 10-8): rs12887440 (RCOR1) and rs35213808 (FSTL5-MIR4454). Additionally, we found associations at ABCG2 and SLC22A12 that were driven by low-frequency SNVs. SNVs in RCOR1 were linked to elevated blood leukocyte messenger RNA levels. THP-1 macrophage culture studies revealed the potential of decreased RCOR1 to suppress gouty inflammation.

CONCLUSION

This is the first comprehensive genetic characterization of adolescent-onset gout. The identified risk loci of early-onset gout mediate inflammatory responsiveness to crystals that could mediate gouty arthritis. This study will contribute to risk prediction and therapeutic interventions to prevent adolescent-onset gout.

metaGE: Investigating genotype x environment interactions through GWAS meta-analysis

Elucidating the genetic components of plant genotype-by-environment interactions is of key importance in the context of increasing climatic instability, diversification of agricultural practices and pest pressure due to phytosanitary treatment limitations. The genotypic response to environmental stresses can be investigated through multi-environment trials (METs). However, genome-wide association studies (GWAS) of MET data are significantly more complex than that of single environments. In this context, we introduce metaGE, a flexible and computationally efficient meta-analysis approach for jointly analyzing single-environment GWAS of any MET experiment. The metaGE procedure accounts for the heterogeneity of quantitative trait loci (QTL) effects across the environmental conditions and allows the detection of QTL whose allelic effect variations are strongly correlated to environmental cofactors. We evaluated the performance of the proposed methodology and compared it to two competing procedures through simulations. We also applied metaGE to two emblematic examples: the detection of flowering QTLs whose effects are modulated by competition in Arabidopsis and the detection of yield QTLs impacted by drought stresses in maize. The procedure identified known and new QTLs, providing valuable insights into the genetic architecture of complex traits and QTL effects dependent on environmental stress conditions. The whole statistical approach is available as an R package.

Interactions of genes with alcohol consumption affect insulin sensitivity and beta cell function

AIMS/HYPOTHESIS

Alcohol consumption has complex effects on diabetes and metabolic disease, but there is widespread heterogeneity within populations and the specific reasons are unclear. Genetic factors may play a role and warrant exploration. The aim of this study was to elucidate genetic variants modulating the impact of alcohol consumption on insulin sensitivity and pancreatic beta cell function within populations presenting normal glucose tolerance (NGT).

METHODS

We recruited 4194 volunteers in Nanjing, 854 in Jurong and an additional 5833 in Nanjing for Discovery cohorts 1 and 2 and a Validation cohort, respectively. We performed an OGTT on all participants, establishing a stringent NGT group, and then assessed insulin sensitivity and beta cell function. Alcohol consumption was categorised as abstinent, light-to-moderate (<210 g per week) or heavy (≥210 g per week). After excluding ineligible individuals, an exploratory genome-wide association study identified potential variants interacting with alcohol consumption in 1862 NGT individuals. These findings were validated in an additional cohort of 2169 NGT individuals. Cox proportional hazard regression was further employed to evaluate the effect of the interaction between the potential variants and alcohol consumption on the risk of type 2 diabetes within the UK Biobank cohort.

RESULTS

A significant correlation was observed between drinking levels and insulin sensitivity, accompanied by a consequent inverse relationship with insulin resistance and beta cell insulin secretion after adjusting for confounding factors in NGT individuals. However, no significant associations were noted in the disposition indexes. The interaction of variant rs56221195 with alcohol intake exhibited a pronounced effect on the liver insulin resistance index (LIRI) in the discovery set, corroborated in the validation set (combined p=1.32 × 10-11). Alcohol consumption did not significantly affect LIRI in rs56221195 wild-type (TT) carriers, but a strong negative association emerged in heterozygous (TA) and homozygous (AA) individuals. The rs56221195 variant also significantly interacts with alcohol consumption, influencing the total insulin secretion index INSR120 (the ratio of the AUC of insulin to glucose from 0 to 120 min) (p=2.06 × 10-9) but not disposition index. In the UK Biobank, we found a significant interaction between rs56221195 and alcohol consumption, which was linked to the risk of type 2 diabetes (HR 0.897, p=0.008).

CONCLUSIONS/INTERPRETATION

Our findings reveal the effects of the interaction of alcohol and rs56221195 on hepatic insulin sensitivity in NGT individuals. It is imperative to weigh potential benefits and detriments thoughtfully when considering alcohol consumption across diverse genetic backgrounds.

Transethnic analysis identifies SORL1 variants and haplotypes protective against Alzheimer's disease

INTRODUCTION

The SORL1 locus exhibits protective effects against Alzheimer's disease (AD) across ancestries, yet systematic studies in diverse populations are sparse.

METHODS

Logistic regression identified AD-associated SORL1 haplotypes in East Asian (N = 5249) and European (N = 8588) populations. Association analysis between SORL1 haplotypes and AD-associated traits or plasma biomarkers was conducted. The effects of non-synonymous mutations were assessed in cell-based systems.

RESULTS

Protective SORL1 variants/haplotypes were identified in the East Asian and European populations. Haplotype Hap_A showed a strong protective effect against AD in East Asians, linked to less severe AD phenotypes, higher SORL1 transcript levels, and plasma proteomic changes. A missense variant within Hap_A, rs2282647-C allele, was linked to a lower risk of AD and decreased expression of a truncated SORL1 protein isoform.

DISCUSSION

Our transethnic analysis revealed key SORL1 haplotypes that exert protective effects against AD, suggesting mechanisms of the protective role of SORL1 in AD.

HIGHLIGHTS

We examined the AD-protective mechanisms of SORL1 in the general population across diverse ancestral backgrounds by jointly analyzing data from three East Asian cohorts (ie, mainland China, Hong Kong, and Japan) and a European cohort. Comparative analysis unveiled key ethnic-specific SORL1 genetic variants and haplotypes. Among these, the SORL1 minor haplotype, Hap_A, emerged as the primary AD-protective factor in East Asians. Hap_A exerts significant AD-protective effects in both APOE ε4 carriers and non-carriers. SORL1 haplotype Hap_A is associated with cognitive function, brain volume, and the activity of specific neuronal and immune-related pathways closely connected to AD risk. Protective variants within Hap_A are linked to increased SORL1 expression in human tissues. We identified an isoform-specific missense variant in Hap_A that modifies the function and levels of a truncated SORL1 protein isoform that is poorly investigated.

Adipose tissue eQTL meta-analysis highlights the contribution of allelic heterogeneity to gene expression regulation and cardiometabolic traits

Complete characterization of the genetic effects on gene expression is needed to elucidate tissue biology and the etiology of complex traits. In the present study, we analyzed 2,344 subcutaneous adipose tissue samples and identified 34,774 conditionally distinct expression quantitative trait locus (eQTL) signals at 18,476 genes. Over half of eQTL genes exhibited at least two eQTL signals. Compared with primary eQTL signals, nonprimary eQTL signals had lower effect sizes, lower minor allele frequencies and less promoter enrichment; they corresponded to genes with higher heritability and higher tolerance for loss of function. Colocalization of eQTLs with genome-wide association study (GWAS) signals for 28 cardiometabolic traits identified 1,835 genes. Inclusion of nonprimary eQTL signals increased discovery of colocalized GWAS-eQTL signals by 46%. Furthermore, 21 genes with ≥2 colocalized GWAS-eQTL signals showed a mediating gene dosage effect on the GWAS trait. Thus, expanded eQTL identification reveals more mechanisms underlying complex traits and improves understanding of the complexity of gene expression regulation.

Association of exposure to second-hand smoke during childhood with blood DNA methylation

INTRODUCTION

By recent estimates, 40% of children worldwide are exposed to second-hand smoke (SHS), which has been associated with adverse health outcomes. While numerous studies have linked maternal smoking during pregnancy (MSDP) to widespread differences in child blood DNA methylation (DNAm), research specifically examining postnatal SHS exposure remains sparse. To address this gap, we conducted epigenome-wide meta-analyses to identify associations of postnatal SHS and child blood DNAm.

METHODS

Six cohorts from the Pregnancy And Childhood Epigenetics (PACE) Consortium (total N = 2,695), with SHS data and child blood DNAm (aged 7-9 years) measured with the Illumina 450K array were included in the meta-analysis. Linear regression models adjusted for covariates were fitted to examine the association between the number of household smokers in postnatal life (0, 1, 2+) and child blood DNAm. Sensitivity models without adjusting for MSDP and restricted to mothers who did not smoke during pregnancy were evaluated.

RESULTS

Our analysis revealed significant associations (False Discovery Rate < 0.05) between household postnatal SHS exposure and DNAm at 11 CpGs in exposed children. Nine CpGs were mapped to genes (MYO1G, FAM184B, CTDSPL2, LTBP3, PDE10A, and FIBCD1), while 2 CpGs were located in open sea regions. Notably, all except 2 CpGs (mapped to FIBCD1 and CTDSPL2) have previously been linked to either personal smoking habits or in utero exposure to smoking. The models restricted to non-smoking mothers provided similar results. Importantly, several of these CpGs and their associated genes are implicated in conditions exacerbated by or directly linked to SHS.

CONCLUSIONS

Our findings highlight the potential biological effects of SHS on blood DNAm. These findings support further research on epigenetic factors mediating deleterious effects of SHS on child health and call for public health policies aimed at reducing exposure, particularly in environments where children are present.

Refining breast cancer genetic risk and biology through multi-ancestry fine-mapping analyses of 192 risk regions

Genome-wide association studies have identified approximately 200 genetic risk loci for breast cancer, but the causal variants and target genes are mostly unknown. We sought to fine-map all known breast cancer risk loci using genome-wide association study data from 172,737 female breast cancer cases and 242,009 controls of African, Asian and European ancestry. We identified 332 independent association signals for breast cancer risk, including 131 signals not reported previously, and for 50 of them, we narrowed the credible causal variants down to a single variant. Analyses integrating functional genomics data identified 195 putative susceptibility genes, enriched in PI3K/AKT, TNF/NF-κB, p53 and Wnt/β-catenin pathways. Single-cell RNA sequencing or in vitro experiment data provided additional functional evidence for 105 genes. Our study uncovered large numbers of association signals and candidate susceptibility genes for breast cancer, uncovered breast cancer genetics and biology, and supported the value of including multi-ancestry data in fine-mapping analyses.

Shared genetic architecture and bidirectional clinical risks within the psycho-metabolic nexus

BACKGROUND

Increasing evidence suggests a complex interplay between psychiatric disorders and metabolic dysregulations. However, most research has been limited to specific disorder pairs, leaving a significant gap in our understanding of the broader psycho-metabolic nexus.

METHODS

This study leveraged large-scale cohort data and genome-wide association study (GWAS) summary statistics, covering 8 common psychiatric disorders and 43 metabolic traits. We introduced a comprehensive analytical strategy to identify shared genetic bases sequentially, from key genetic correlation regions to local pleiotropy and pleiotropic genes. Finally, we developed polygenic risk score (PRS) models to translate these findings into clinical applications.

FINDINGS

We identified significant bidirectional clinical risks between psychiatric disorders and metabolic dysregulations among 310,848 participants from the UK Biobank. Genetic correlation analysis confirmed 104 robust trait pairs, revealing 1088 key genomic regions, including critical hotspots such as chr3: 47588462-50387742. Cross-trait meta-analysis uncovered 388 pleiotropic single nucleotide variants (SNVs) and 126 shared causal variants. Among variants, 45 novel SNVs were associated with psychiatric disorders and 75 novel SNVs were associated with metabolic traits, shedding light on new targets to unravel the mechanism of comorbidity. Notably, RBM6, a gene involved in alternative splicing and cellular stress response regulation, emerged as a key pleiotropic gene. When psychiatric and metabolic genetic information were integrated, PRS models demonstrated enhanced predictive power.

INTERPRETATION

The study highlights the intertwined genetic and clinical relationships between psychiatric disorders and metabolic dysregulations, emphasising the need for integrated approaches in diagnosis and treatment.

FUNDING

The National Key Research and Development Program of China (2023YFC2506200, SHH). The National Natural Science Foundation of China (82273741, SY).