Random-effects model aimed at discovering associations in meta-analysis of genome-wide association studies

A GWAS of angiotensin-converting enzyme inhibitor-induced angioedema in a South African population

Background

Angiotensin-converting enzyme inhibitor-induced angioedema (AE-ACEI) is a life-threatening adverse event; globally, it is the most common cause of emergency presentations with angioedema. Several genome-wide association studies (GWASs) have found genomic associations with AE-ACEI. However, despite African Americans having a 5-fold increased risk of AE-ACEI, there are no published GWASs from Africa.

Objective

The aim of this study was to conduct a GWAS of AE-ACEI in a South African population and perform a meta-analysis with an African American and European American population.

Methods

The GWAS included 202 South African adults with a history of AE-ACEI and 513 controls without angioedema following angiotensin-converting enzyme inhibitor (ACEI) treatment for at least 2 years. A meta-analysis was conducted with GWAS summary statistics from an African American and European American cohort (from the Vanderbilt-Marshfield cohort, which consisted of 174 case patients and 489 controls).

Results

No single-nucleotide polymorphisms (SNPs) attained genome-wide significance; however, 26 SNPs in the postimputation standard GWAS of the South African cohort and 73 SNPs in the meta-analysis attained suggestive thresholds (P < 5.0 × 10-06). Some of these SNPs were found to be located close to the genes PRKCQ (protein kinase C theta), RAD51B (RAD51 Paralog B), and RIMS1 (regulating synaptic membrane exocytosis 1), which were previously linked with drug-induced angioedema, and also close to the CSMD1 (CUB and sushi multiple domains 1) gene, which has been linked to ACEI cough, providing replication at the gene level but with novel lead SNPs. The study also replicated SNP rs500766 on chromosome 10, which was previously found to be associated with AE-ACEI.

Conclusions

Our results highlight the importance of African populations for detection of novel variants in replication studies. Further increased sampling across the continent and matched functional work are needed to confirm the importance of genetic variation in understanding the biology of AE-ACEI.

Blunted Niacin Skin Flushing Response in Schizophrenia: A Meta-analysis

BACKGROUND AND HYPOTHESIS

The multifactorial pathogenesis of schizophrenia (SZ) hinders the diagnosis and treatment of this disorder. Niacin skin flushing response (NSFR) has been identified as an endophenotype for SZ, but the proportion of blunted NSFR (BNR) varied between studies. This study aims to clarify the relationship between NSFR and SZ through a meta-analysis.

STUDY DESIGN

PubMed, Embase, Web of Science, Cochrane, and Scopus databases were searched for articles published until May 2024, and 32 studies were eligible. Using random-effects models, we examined the characteristics of NSFR in SZ, including the reaction degree, speed, sensitivity, and risk and prevalence of BNR. Subgroup analyses and regression analyses were performed to investigate the relevant effect factors of NSFR.

STUDY RESULTS

The reaction degree (SMD = -0.90; CI, -1.08 to -0.72), speed (SMD = 0.64; CI, 0.02-1.25), and sensitivity (SMD = 0.89; CI, 0.49-1.29) of NSFR was significantly reduced in SZ compared to healthy controls (HC). Moreover, we observed a positive association between BNR and SZ (OR = 8.50; CI, 5.93-12.19). The overall prevalence of BNR was 58.5% in SZ (CI, 49.3%-67.8%) compared to 11.8% in HC (CI, 7.7%-15.9%). In addition, NSFR detection method, geographical regions, and age were found to have effects on reaction degree and prevalence of BNR.

CONCLUSIONS

This study confirmed a significantly abnormal NSFR and higher prevalence of BNR in SZ, which highlights the potential facilitation of the diagnosis and personalized intervention of SZ subgroups. In addition, the study points to a need to establish a standardized method for NSFR assessment.

An integrative systems-biology approach defines mechanisms of Alzheimer's disease neurodegeneration

Despite years of intense investigation, the mechanisms underlying neuronal death in Alzheimer's disease, remain incompletely understood. To define relevant pathways, we conducted an unbiased, genome-scale forward genetic screen for age-associated neurodegeneration in Drosophila. We also measured proteomics, phosphoproteomics, and metabolomics in Drosophila models of Alzheimer's disease and identified Alzheimer's genetic variants that modify gene expression in disease-vulnerable neurons in humans. We then used a network model to integrate these data with previously published Alzheimer's disease proteomics, lipidomics and genomics. Here, we computationally predict and experimentally confirm how HNRNPA2B1 and MEPCE enhance toxicity of the tau protein, a pathological feature of Alzheimer's disease. Furthermore, we demonstrated that the screen hits CSNK2A1 and NOTCH1 regulate DNA damage in Drosophila and human stem cell-derived neural progenitor cells. Our study identifies candidate pathways that could be targeted to ameliorate neurodegeneration in Alzheimer's disease.

Methodological opportunities in genomic data analysis to advance health equity

The causes and consequences of inequities in genomic research and medicine are complex and widespread. However, it is widely acknowledged that underrepresentation of diverse populations in human genetics research risks exacerbating existing health disparities. Efforts to improve diversity are ongoing, but an often-overlooked source of inequity is the choice of analytical methods used to process, analyse and interpret genomic data. This choice can influence all areas of genomic research, from genome-wide association studies and polygenic score development to variant prioritization and functional genomics. New statistical and machine learning techniques to understand, quantify and correct for the impact of biases in genomic data are emerging within the wider genomic research and genomic medicine ecosystems. At this crucial time point, it is important to clarify where improvements in methods and practices can, or cannot, have a role in improving equity in genomics. Here, we review existing approaches to promote equity and fairness in statistical analysis for genomics, and propose future methodological developments that are likely to yield the most impact for equity.

Investigating the shared genetic architecture between obesity and depression: a large-scale genomewide cross-trait analysis

Introduction

Increasing evidence suggests that individuals with obesity are at a higher risk of developing depression, and conversely, depression can contribute to the onset of obesity, creating a detrimental cycle. This study aims to investigate the potential shared biological pathways between obesity and depression by examining genetic correlations, identifying common polymorphisms, and conducting cross-trait genetic analyses.

Methods

We assessed the genetic correlation between obesity and depression using linkage disequilibrium score regression and high-density lipoprotein levels. We combined two different sources of obesity data using METAL and employed bidirectional Mendelian randomization to determine the causal relationship between obesity and depression. Additionally, we conducted multivariate trait analysis using the MTAG method to improve statistical robustness and identify novel genetic associations. Furthermore, we performed a thorough investigation of independent risk loci using GCTA-COJO, PLACO, MAGMA, POPS, and SMR, integrating different QTL information and methods to further identify risk genes and proteins.

Results

Our analysis revealed genetic correlations and bidirectional positive causal relationships between obesity and depression, highlighting shared risk SNP (rs10789340). We identified RPL31P12, NEGR1, and DCC as common risk genes for obesity and depression. Using the BLISS method, we identified SCG3 and FLRT2 as potential drug targets.

Limitation

Most of our data sources are from Europe, which may limit the generalization of our findings to other ethnic populations.

Conclusion

This study demonstrates the genetic causal relationship and common risk SNPs, genes, proteins, and pathways between obesity and depression. These findings contribute to a deeper understanding of their pathogenesis and the identification of potential therapeutic targets.

Genome-wide association study for lung cancer in 6531 African Americans reveals new susceptibility loci

Despite lung cancer affecting all races and ethnicities, disparities are observed in incidence and mortality rates among different ethnic groups in the United States. Non-Hispanic African Americans had a high incidence rate of lung cancer at 55.8 per 100 000 people, as well as the highest death rate at 37.2 per 100 000 people from 2016 to 2020. While previous genome-wide association studies (GWAS) have identified over 45 susceptibility risk loci that influence lung cancer development, few GWAS have investigated the etiology of lung cancer in African Americans. To address this gap in knowledge, we conducted GWAS of lung cancer focused on studying African Americans, comprising 2267 lung cancer cases and 4264 controls. We identified three loci associated with lung cancer, one with lung adenocarcinoma, and four with lung squamous cell carcinoma in this population at the genomic-wide significance level. Among them, three novel loci were identified near VWF at 12p13.31 for overall lung cancer and GACAT3 at 2p24.3 and LMAN1L at 15q24.1 for lung squamous cell carcinoma. In addition, we confirmed previously reported risk loci with known or new lead variants near CHRNA5 at 15q25.1 and CYP2A6 at 19q13.2 associated with lung cancer and TRIP13 at 5p15.33 and ERC1 at 12p13.33 associated with lung squamous cell carcinoma. Further multi-step functional analyses shed light on biological mechanisms underlying these associations of lung cancer in this population. Our study highlights the importance of ancestry-specific studies for the potential alleviation of lung cancer burden in African Americans.

Detection of positive selection on depression-associated genes

Although depression significantly impacts fitness, some hypotheses suggest that it may offer a survival benefit. However, there has been limited systematic investigation into the selection pressures acting on genes associated with depression at the genomic level. Here, we conducted comparative genomic analyses and computational molecular evolutionary analyses on 320 depression-associated genes at two levels, i.e., across the primate phylogeny (long timescale selection) and in modern human populations (recent selection). We identified seven genes under positive selection in the human lineage, and 46 genes under positive selection in modern human populations. Most positively selected variants in modern human populations were at UTR regions and non-coding exons, indicating the importance of gene expression regulation in the evolution of depression-associated genes. Positively selected genes are not only related to immune responses, but also function in reproduction and dietary adaptation. Notably, the proportion of depression-associated genes under positive selection was significantly higher than the positively selected genes at the genome-wide average level in African, East Asian, and South Asian populations. We also identified two positively selected loci that happened to be associated with depression in the South Asian population. Our study revealed that depression-associated genes are subject to varying selection pressures across different populations. We suggest that, in precision medicine-particularly in gene therapy-it is crucial to consider the specific functions of genes within distinct populations.

Effects of exercise intervention on exercise capacity and cardiopulmonary function in patients with atrial fibrillation: A randomized controlled trial systematic review and meta-analysis

BACKGROUND

Atrial fibrillation (AF) is a common cardiac arrhythmia that significantly impacts the cardiopulmonary function and quality of life of patients. Despite various treatment strategies, non-pharmacological interventions, particularly exercise interventions, have gained attention in recent years.

OBJECTIVE

Through systematic review and meta-analysis, this study explores the impact of physical activity on the exercise capacity and quality of life of AF patients. It assesses the safety, clinical outcomes, and physiological mechanisms of exercise intervention in the treatment of AF.

METHODS

The systematic review and individual patient data (IPD) meta-analysis method were employed, following the PRISMA-IPD guidelines, for literature selection, data extraction, and quality assessment. The analysis focused on the impact of exercise on the cardiopulmonary function and quality of life of AF patients in randomized controlled trials.

RESULTS

A total of 12 randomized controlled trials involving 287 AF patients were included. Meta-analysis demonstrated a significant improvement in the 6-minute walk test capacity (SMD=87.87, 95% CI [42.23, 133.51]), static heart rate improvement (SMD=-7.63, 95% CI [-11.42, -3.85]), and cardiopulmonary function enhancement (SMD=2.37, 95% CI [0.96, 3.77]) due to exercise. There was also a significant improvement in the quality of life (SMD=0.720, 95% CI [0.038, 1.402]).

CONCLUSION

Exercise has a significant effect on improving exercise capacity and cardiopulmonary function in patients with atrial fibrillation. Particularly, high-intensity exercise training has a more significant impact on improving cardiopulmonary function and exercise capacity, emphasizing the importance of personalized exercise plans in enhancing the cardiopulmonary health of AF patients. Further research is needed to explore the effects of exercise on improving the quality of life in the future.

PROSPERO ID

CRD2023493917.

Translational genomics of osteoarthritis in 1,962,069 individuals

Osteoarthritis is the third most rapidly growing health condition associated with disability, after dementia and diabetes1. By 2050, the total number of patients with osteoarthritis is estimated to reach 1 billion worldwide2. As no disease-modifying treatments exist for osteoarthritis, a better understanding of disease aetiopathology is urgently needed. Here we perform a genome-wide association study meta-analyses across up to 489,975 cases and 1,472,094 controls, establishing 962 independent associations, 513 of which have not been previously reported. Using single-cell multiomics data, we identify signal enrichment in embryonic skeletal development pathways. We integrate orthogonal lines of evidence, including transcriptome, proteome and epigenome profiles of primary joint tissues, and implicate 700 effector genes. Within these, we find rare coding-variant burden associations with effect sizes that are consistently higher than common frequency variant associations. We highlight eight biological processes in which we find convergent involvement of multiple effector genes, including the circadian clock, glial-cell-related processes and pathways with an established role in osteoarthritis (TGFβ, FGF, WNT, BMP and retinoic acid signalling, and extracellular matrix organization). We find that 10% of the effector genes express a protein that is the target of approved drugs, offering repurposing opportunities, which can accelerate translation.

Cross-ancestry and sex-stratified genome-wide association analyses of amygdala and subnucleus volumes

The amygdala is a small but critical multi-nucleus structure for emotion, cognition and neuropsychiatric disorders. Although genetic associations with amygdala volumetric traits have been investigated in sex-combined European populations, cross-ancestry and sex-stratified analyses are lacking. Here we conducted cross-ancestry and sex-stratified genome-wide association analyses for 21 amygdala volumetric traits in 6,923 Chinese and 48,634 European individuals. We identified 191 variant-trait associations (P < 2.38 × 10-9), including 47 new associations (12 new loci) in sex-combined univariate analyses and seven additional new loci in sex-combined and sex-stratified multivariate analyses. We identified 12 ancestry-specific and two sex-specific associations. The identified genetic variants include 16 fine-mapped causal variants and regulate amygdala and fetal brain gene expression. The variants were enriched for brain development and colocalized with mood, cognition and neuropsychiatric disorders. These results indicate that cross-ancestry and sex-stratified genetic association analyses may provide insight into the genetic architectures of amygdala and subnucleus volumes.

Epigenome-Wide Association Study of Asthma Exacerbations in Europeans

BACKGROUND

Asthma exacerbations (AEs) represent the major contributor to the global asthma burden. Although genetic and environmental factors have been associated with AEs, the role of epigenetics remains uncovered.

OBJECTIVE

This study aimed to identify whole blood DNA methylation (DNAm) markers associated with AEs in Europeans.

METHODS

DNAm was assessed in 406 blood samples from Spanish individuals using the Infinium MethylationEPIC microarray (Illumina). An epigenome-wide association study was conducted to test the association of DNAm with AEs at differentially methylated positions, regions, and epigenetic modules. CpGs suggestively associated with AEs (false discovery rate [FDR] < 0.1) were followed up for replication in 222 European individuals, and the genome-wide significance (p < 9 × 10-8) was declared after meta-analyzing the discovery and replication samples. Additional assessment was performed using nasal tissue DNAm data from 155 Spanish individuals. The effects of genetic variation on DNAm were assessed through cis-methylation quantitative trait loci (meQTL) analysis. Enrichment analyses of previous EWAS signals were conducted.

RESULTS

Four CpGs were associated with AEs, and two were replicated and reached genomic significance in the meta-analysis (annotated to ZBTB16 and BAIAP2). Of those, CpG cg25345365 (ZBTB16) was cross-tissue validated in nasal epithelium (p= 0.003) and associated with five independent meQTLs (FDR < 0.05). Additionally, four differentially methylated regions and one module were significantly associated with AEs. Enrichment analyses revealed an overrepresentation of prior epigenetic associations with prenatal and environmental exposures, immune-mediated diseases, and mortality.

CONCLUSIONS

DNAm in whole blood and nasal samples may contribute to AEs in Europeans, capturing genetic and environmental risk factors.

Assessment and ascertainment in psychiatric molecular genetics: challenges and opportunities for cross-disorder research

Psychiatric disorders are highly comorbid, heritable, and genetically correlated [1-4]. The primary objective of cross-disorder psychiatric genetics research is to identify and characterize both the shared genetic factors that contribute to convergent disease etiologies and the unique genetic factors that distinguish between disorders [4, 5]. This information can illuminate the biological mechanisms underlying comorbid presentations of psychopathology, improve nosology and prediction of illness risk and trajectories, and aid the development of more effective and targeted interventions. In this review we discuss how estimates of comorbidity and identification of shared genetic loci between disorders can be influenced by how disorders are measured (phenotypic assessment) and the inclusion or exclusion criteria in individual genetic studies (sample ascertainment). Specifically, the depth of measurement, source of diagnosis, and time frame of disease trajectory have major implications for the clinical validity of the assessed phenotypes. Further, biases introduced in the ascertainment of both cases and controls can inflate or reduce estimates of genetic correlations. The impact of these design choices may have important implications for large meta-analyses of cohorts from diverse populations that use different forms of assessment and inclusion criteria, and subsequent cross-disorder analyses thereof. We review how assessment and ascertainment affect genetic findings in both univariate and multivariate analyses and conclude with recommendations for addressing them in future research.

Novel Parkinson's Disease Genetic Risk Factors Within and Across European Populations

Introduction

We conducted a meta-analysis of Parkinson's disease genome-wide association study summary statistics, stratified by source (clinically-recruited case-control cohorts versus population biobanks) and by general European versus European isolate ancestries. This study included 63,555 cases, 17,700 proxy cases with a family history of Parkinson's disease, and 1,746,386 controls, making it the largest investigation of Parkinson's disease genetic risk to date.

Methods

Meta-analyses were performed using standard fixed and random effect models for the European sub-populations, the case-control studies, and the population biobanks separately. Finally, all of the European ancestries for all study types as well as proxy cases were combined in our final cross-European meta-analysis. We estimated heritable risk across ancestry groups, investigated tissue and cell-type enrichment, and prioritized risk genes using public data to facilitate functional follow-up efforts.

Results

The final combined cross-European meta-analysis identified 134 risk loci (59 novel), with a total of 157 independent signals, significantly expanding our understanding of Parkinson's disease risk. Multi-omic data integration revealed that expression of the nominated risk genes are highly enriched in brain tissues, particularly in neuronal and astrocyte cell types. Additionally, we prioritized 33 high-confidence genes across these 134 loci for future follow-up studies.

Conclusions

By integrating diverse European populations and leveraging harmonized data from the Global Parkinson's Genetics Program (GP2), we reveal new insight into the genetic architecture of Parkinson's disease. We identified a total of 134 risk loci, expanding the number of known loci associated with PD by approximately 24%. We also provided an initial layer of biological context to these results through follow-up analyses in an effort to facilitate follow-up studies and precision medicine efforts with the goal of advancing Parkinson's disease research.

Contribution of germline and somatic mutations to risk of neuromyelitis optica spectrum disorder

Neuromyelitis optica spectrum disorder (NMOSD) is a rare autoimmune disease characterized by optic neuritis and transverse myelitis, with an unclear genetic background. A genome-wide meta-analysis of NMOSD in Japanese individuals (240 patients and 50,578 controls) identified significant associations with the major histocompatibility complex region and a common variant close to CCR6 (rs12193698; p = 1.8 × 10-8, odds ratio [OR] = 1.73). In single-cell RNA sequencing (scRNA-seq) analysis (25 patients and 101 controls), the CCR6 risk variant showed disease-specific expression quantitative trait loci effects in CD4+ T (CD4T) cell subsets. Furthermore, we detected somatic mosaic chromosomal alterations (mCAs) in various autoimmune diseases and found that mCAs increase the risk of NMOSD (OR = 3.37 for copy number alteration). In scRNA-seq data, CD4T cells with 21q loss, a recurrently observed somatic event in NMOSD, showed dysregulation of type I interferon-related genes. Our integrated study identified novel germline and somatic mutations associated with NMOSD pathogenesis.

A genome-wide association study identifies a novel East Asian-specific locus for dementia with Lewy bodies in Japanese subjects

BACKGROUND

Dementia with Lewy bodies (DLB) is the second most common type of degenerative dementia in older patients. As with other multifactorial diseases, the pathogenesis results from interactions of environmental and genetic factors. The genetic basis of DLB is not yet fully understood. Recent genomic analyses of DLB in Caucasian cohorts identified genetic susceptibility loci for DLB, but the comprehensive genomic analysis in Asians was still not performed.

METHODS

We conducted a genome-wide association study (GWAS) in Japanese subjects (211 DLB cases and 6113 controls) to clarify the genetic architecture of DLB pathogenesis.

RESULTS

We identified the East Asian-specific DHTKD1 locus (rs138587229) on chromosome 10 with genome-wide significance (GWS; P = 3.27 × 10-8) and the ICOS/PARD3B locus on chromosome 2 with suggestive significance (P = 3.95 × 10-7) as novel DLB genetic risk loci. We also confirmed the APOE locus (rs429358, P < 5.0 × 10-8), a known risk locus for DLB and Alzheimer's disease in Caucasians. The DHTKD1 locus was associated with the gene expression of SEC61A2 and showed a causal relationship with cholinesterase levels. In a trans-ethnic meta-analysis that included Japanese, UK Biobank, and other Caucasian GWAS, we confirmed the risk for DLB at APOE and SNCA loci with GWS. Transcriptome-wide association analysis identified ZNF155 and ZNF284 in the brain cortex and GPRIN3 in the substantia nigra as putative causal genes for DLB.

CONCLUSIONS

This is the first GWAS for DLB in East Asians, and our findings provide new biological and clinical insights into the pathogenesis of DLB.

Structural variants linked to Alzheimer's disease and other common age-related clinical and neuropathologic traits

BACKGROUND

Alzheimer's disease (AD) is a complex neurodegenerative disorder with substantial genetic influence. While genome-wide association studies (GWAS) have identified numerous risk loci for late-onset AD (LOAD), the functional mechanisms underlying most of these associations remain unresolved. Large genomic rearrangements, known as structural variants (SVs), represent a promising avenue for elucidating such mechanisms within some of these loci.

METHODS

By leveraging data from two ongoing cohort studies of aging and dementia, the Religious Orders Study and Rush Memory and Aging Project (ROS/MAP), we performed genome-wide association analysis testing 20,205 common SVs from 1088 participants with whole genome sequencing (WGS) data. A range of Alzheimer's disease and other common age-related clinical and neuropathologic traits were examined.

RESULTS

First, we mapped SVs across 81 AD risk loci and discovered 22 SVs in linkage disequilibrium (LD) with GWAS lead variants and directly associated with the phenotypes tested. The strongest association was a deletion of an Alu element in the 3'UTR of the TMEM106B gene, in high LD with the respective AD GWAS locus and associated with multiple AD and AD-related disorders (ADRD) phenotypes, including tangles density, TDP-43, and cognitive resilience. The deletion of this element was also linked to lower TMEM106B protein abundance. We also found a 22-kb deletion associated with depression in ROS/MAP and bearing similar association patterns as GWAS SNPs at the IQCK locus. In addition, we leveraged our catalog of SV-GWAS to replicate and characterize independent findings in SV-based GWAS for AD and five other neurodegenerative diseases. Among these findings, we highlight the replication of genome-wide significant SVs for progressive supranuclear palsy (PSP), including markers for the 17q21.31 MAPT locus inversion and a 1483-bp deletion at the CYP2A13 locus, along with other suggestive associations, such as a 994-bp duplication in the LMNTD1 locus, suggestively linked to AD and a 3958-bp deletion at the DOCK5 locus linked to Lewy body disease (LBD) (P = 3.36 × 10-4).

CONCLUSIONS

While still limited in sample size, this study highlights the utility of including analysis of SVs for elucidating mechanisms underlying GWAS loci and provides a valuable resource for the characterization of the effects of SVs in neurodegenerative disease pathogenesis.

The prevalence of comorbid mental health difficulties in young people with chronic skin conditions: A systematic review and meta-analysis

Chronic skin conditions can have psychosocial and somatic implications, influencing well-being and quality of life. This systematic review and meta-analysis aimed to synthesise evidence on the prevalence of comorbid mental health difficulties in 0-25-year-olds with chronic skin conditions. A secondary aim included identifying factors associated with resilience. The narrative synthesis included 45 studies. Four meta-analyses were performed with moderate-high quality studies, one for each outcome: diagnosed mental disorders; mental health symptoms; suicidal behaviour; socio-emotional and behavioural difficulties. The pooled prevalence of diagnosed mental disorders was 1.2% (95% CI = 0.2-6.1); of mental health symptoms was 22.6% (95% CI = 18.9-26.7); of suicidal behaviour was 7.8% (95% CI = 1.4-3.1); of socio-emotional and behavioural difficulties was 20.9% (95% CI = 14.7-28.8). Findings demonstrate the pooled prevalence of comorbid mental health difficulties in youth with chronic skin conditions.

Genome-wide association analysis using multiple Atlantic salmon populations

BACKGROUND

In a previous study, we found low persistence of linkage disequilibrium (LD) phase across breeding populations of Atlantic salmon. Accordingly, we observed no increase in accuracy from combining these populations for genomic prediction. In this study, we aimed to examine if the same were true for detection power in genome-wide association studies (GWAS), in terms of reduction in p-values, and if the precision of mapping quantitative trait loci (QTL) would improve from such analysis. Since individual records may not always be available, e.g. due to proprietorship or confidentiality, we also compared mega-analysis and meta-analysis. Mega-analysis needs access to all individual records, whereas meta-analysis utilizes parameters, such as p-values or allele substitution effects, from multiple studies or populations. Furthermore, different methods for determining the presence or absence of independent or secondary signals, such as conditional association analysis, approximate conditional and joint analysis (COJO), and the clumping approach, were assessed.

RESULTS

Mega-analysis resulted in increased detection power, in terms of reduction in p-values, and increased precision, compared to the within-population GWAS. Only one QTL was detected using conditional association analysis, both within populations and in mega-analysis, while the number of QTL detected with COJO and the clumping approach ranged from 1 to 19. The allele substitution effect and -log10p-values obtained from mega-analysis were highly correlated with the corresponding values from various meta-analysis methods. Compared to mega-analysis, a higher detection power and reduced precision were obtained with the meta-analysis methods.

CONCLUSIONS

Our results show that combining multiple datasets or populations in a mega-analysis can increase detection power and mapping precision. With meta-analysis, a higher detection power was obtained compared to mega-analysis. However, care must be taken in the interpretation of the meta-analysis results from multiple populations because their test statistics might be inflated due to population structure or cryptic relatedness.

Trans-ancestry Genome-Wide Analyses in UK Biobank Yield Novel Risk Loci for Major Depression

Most genome-wide association studies (GWASs) of depression focus on broad, heterogeneous outcomes, limiting the discovery of genomic risk loci specific to major depressive disorder (MDD). Previous UK Biobank (UKB) studies had limited ability to pinpoint MDD-associated loci due to a smaller sample with strictly defined MDD outcomes and further exclusion of many participants based on ancestry or relatedness, significantly underutilizing this resource's potential for elucidating the genetic architecture of MDD. Here, we present novel genomic insights into MDD by fully utilizing existing UKB data through (1) a trans-ancestry GWAS pipeline using two complementary approaches controlling for population structure and relatedness and (2) an increased sample with MDD symptom-level data across two mental health assessments. We identified strict MDD outcomes among 211,535 participants, representing a 38% increase in eligible participants from prior studies with only one assessment. Ancestrally inclusive analyses yielded 61 genomic risk loci across depression phenotypes, compared to 47 in the analyses restricted to participants genetically similar to European ancestry. Fourteen of these loci, including five novel, were associated with strict MDD phenotypes, whereas only one locus has been previously reported in UKB. MDD-associated genomic loci and predicted gene expression levels showed little overlap with broad depression, indicating higher specificity. Notably, polygenic scores based on these results were significantly associated with depression diagnoses across ancestry groups in the All of Us Research Program, highlighting the shared genetic architecture across populations. While the trans-ancestry analyses, which included non-European participants, increased the number of associated loci, the discovery of non-European ancestry-specific loci was limited, underscoring the need for larger, globally representative studies of MDD. Importantly, beyond these results, our GWAS pipeline will facilitate inclusive analyses of other traits and disorders, helping improve statistical power, representation, and generalizability in genomic studies.

Sequence-based GWAS in 180,000 German Holstein cattle reveals new candidate variants for milk production traits

BACKGROUND

Milk production traits are complex and influenced by many genetic and environmental factors. Although extensive research has been performed for these traits, with many associations unveiled thus far, due to their crucial economic importance, complex genetic architecture, and the fact that causal variants in cattle are still scarce, there is a need for a better understanding of their genetic background. In this study, we aimed to identify new candidate loci associated with milk production traits in German Holstein cattle, the most important dairy breed in Germany and worldwide. For that purpose, 180,217 cattle were imputed to the sequence level and large-scale genome-wide association study (GWAS) followed by fine-mapping and evolutionary and functional annotation were carried out to identify and prioritize new association signals.

RESULTS

Using the imputed sequence data of a large cattle dataset, we identified 50,876 significant variants, confirming many known and identifying previously unreported candidate variants for milk (MY), fat (FY), and protein yield (PY). Genome-wide significant signals were fine-mapped with the Bayesian approach that determines the credible variant sets and generates the probability of causality for each signal. The variants with the highest probabilities of being causal were further classified using external information about the function and evolution, making the prioritization for subsequent validation experiments easier. The top potential causal variants determined with fine-mapping explained a large percentage of genetic variance compared to random ones; 178 variants explained 11.5%, 104 explained 7.7%, and 68 variants explained 3.9% of the variance for MY, FY, and PY, respectively, demonstrating the potential for causality.

CONCLUSIONS

Our findings proved the power of large samples and sequence-based GWAS in detecting new association signals. In order to fully exploit the power of GWAS, one should aim at very large samples combined with whole-genome sequence data. These can also come with both computational and time burdens, as presented in our study. Although milk production traits in cattle are comprehensively investigated, the genetic background of these traits is still not fully understood, with the potential for many new associations to be revealed, as shown. With constantly growing sample sizes, we expect more insights into the genetic architecture of milk production traits in the future.

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.

Using clustering of genetic variants in Mendelian randomization to interrogate the causal pathways underlying multimorbidity from a common risk factor

Mendelian randomization (MR) is an epidemiological approach that utilizes genetic variants as instrumental variables to estimate the causal effect of an exposure on a health outcome. This paper investigates an MR scenario in which genetic variants aggregate into clusters that identify heterogeneous causal effects. Such variant clusters are likely to emerge if they affect the exposure and outcome via distinct biological pathways. In the multi-outcome MR framework, where a shared exposure causally impacts several disease outcomes simultaneously, these variant clusters can provide insights into the common disease-causing mechanisms underpinning the co-occurrence of multiple long-term conditions, a phenomenon known as multimorbidity. To identify such variant clusters, we adapt the general method of agglomerative hierarchical clustering to multi-sample summary-data MR setup, enabling cluster detection based on variant-specific ratio estimates. Particularly, we tailor the method for multi-outcome MR to aid in elucidating the causal pathways through which a common risk factor contributes to multiple morbidities. We show in simulations that our "MR-AHC" method detects clusters with high accuracy, outperforming the existing methods. We apply the method to investigate the causal effects of high body fat percentage on type 2 diabetes and osteoarthritis, uncovering interconnected cellular processes underlying this multimorbid disease pair.

A Mixed-Effect Kernel Machine Regression Model for Integrative Analysis of Alpha Diversity in Microbiome Studies

Increasing evidence suggests that human microbiota plays a crucial role in many diseases. Alpha diversity, a commonly used summary statistic that captures the richness and/or evenness of the microbial community, has been associated with many clinical conditions. However, individual studies that assess the association between alpha diversity and clinical conditions often provide inconsistent results due to insufficient sample size, heterogeneous study populations and technical variability. In practice, meta-analysis tools have been applied to integrate data from multiple studies. However, these methods do not consider the heterogeneity caused by sequencing protocols, and the contribution of each study to the final model depends mainly on its sample size (or variance estimate). To combine studies with distinct sequencing protocols, a robust statistical framework for integrative analysis of microbiome datasets is needed. Here, we propose a mixed-effect kernel machine regression model to assess the association of alpha diversity with a phenotype of interest. Our approach readily incorporates the study-specific characteristics (including sequencing protocols) to allow for flexible modeling of microbiome effect via a kernel similarity matrix. Within the proposed framework, we provide three hypothesis testing approaches to answer different questions that are of interest to researchers. We evaluate the model performance through extensive simulations based on two distinct data generation mechanisms. We also apply our framework to data from HIV reanalysis consortium to investigate gut dysbiosis in HIV infection.

Effect of Epstein-Barr Virus infection on gene regulation in immune cells of patients with Immune-Mediated Diseases

It has been known that Epstein-Barr virus (EBV) can latently infect immune cells after the initial infection, and epidemiological studies have suggested its association with the onset of immune-mediated diseases (IMDs). However, the specific impact of EBV infection on IMDs pathology remains unclear. We quantified EBV load of B cell subsets (Naïve B cells, Unswitched memory B cells, Switched memory B cells, Double negative B cells, and Plasmablasts) in IMD patients as well as healthy control (HC) using bulk RNA sequencing data of 504 donors. The EBV load was clearly higher in IMD patients compared to HC. Furthermore, the wide range of EBV load in this dataset enabled us to assess the impact of EBV load on gene regulation. We found many examples of expression quantitative trait loci (eQTL) whose effects were associated with EBV load. Expression QTLs that exhibited larger effects with increasing EBV load were significantly overlapped with binding sites of transcription factors derived from the EBV genome. These EBV load-associated eQTLs exhibited high enrichment of systemic lupus erythematosus (SLE) GWAS signals, suggesting the mechanical link of EBV infection and the onset of the disease via gene regulation. These findings provide the first evidence of the influence of EBV infection on gene regulation in human primary cells and its association with the SLE onset and/or progression.

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).

The Impact of Ancestry on Genome-Wide Association Studies

Genome-wide association studies (GWAS) are an important tool for the study of complex disease genetics. Decisions regarding the quality control (QC) procedures employed as part of a GWAS can have important implications on the results and their biological interpretation. Many GWAS have been conducted predominantly in cohorts of European ancestry, but many initiatives aim to increase the representation of diverse ancestries in genetic studies. The question of how these data should be combined and the consequences that genetic variation across ancestry groups might have on GWAS results warrants further investigation. In this study, we focus on several commonly used methods for combining genetic data across diverse ancestry groups and the impact these decisions have on the outcome of GWAS summary statistics. We ran GWAS on two binary phenotypes using ancestry-specific, multi-ancestry mega-analysis, and meta-analysis approaches. We found that while multi-ancestry mega-analysis and meta-analysis approaches can aid in identifying signals shared across ancestries, they can diminish the signal of ancestry-specific associations and modify their effect sizes. These results demonstrate the potential impact on downstream post-GWAS analyses and follow-up studies. Decisions regarding how the genetic data are combined has the potential to mask important findings that might serve individuals of ancestries that have been historically underrepresented in genetic studies. New methods that consider ancestry-specific variants in conjunction with the shared variants need to be developed.

Impact of disease-associated chromatin accessibility QTLs across immune cell types and contexts

Only a third of immune-associated loci from genome-wide association studies (GWAS) colocalize with expression quantitative trait loci (eQTLs). To learn about causal genes and mechanisms at the remaining loci, we created a unified single-cell chromatin accessibility (scATAC-seq) map in peripheral blood comprising a total of 282,424 cells from 48 individuals. Clustering and topic modeling of scATAC data identified discrete cell-types and continuous cell states, which helped reveal disease-relevant cellular contexts, and allowed mapping of genetic effects on chromatin accessibility across these contexts. We identified 37,390 chromatin accessibility QTLs (caQTL) at 10% FDR across eight cell groups and observed extensive sharing of caQTLs across immune cell contexts, finding that fewer than 20% of caQTLs are specific to a single cell type. Notably, caQTLs colocalized with ∼50% more GWAS loci compared to eQTLs, helping to nominate putative causal genes for many unexplained loci. However, most GWAS-caQTL colocalizations had no detectable downstream regulatory effects on gene expression levels in the same cell type. We find evidence that the higher rates of colocalization between caQTLs and GWAS signals reflect missing disease-relevant cellular contexts among existing eQTL studies. Thus, there remains a pressing need for identifying disease-causing cellular contexts and for mapping gene regulatory variation in these cells.

Genome-wide association analyses identify distinct genetic architectures for age-related macular degeneration across ancestries

To effectively reduce vision loss due to age-related macular generation (AMD) on a global scale, knowledge of its genetic architecture in diverse populations is necessary. A critical element, AMD risk profiles in African and Hispanic/Latino ancestries, remains largely unknown. We combined data in the Million Veteran Program with five other cohorts to conduct the first multi-ancestry genome-wide association study of AMD and discovered 63 loci (30 novel). We observe marked cross-ancestry heterogeneity at major risk loci, especially in African-ancestry populations which demonstrate a primary signal in a major histocompatibility complex class II haplotype and reduced risk at the established CFH and ARMS2/HTRA1 loci. Dissecting local ancestry in admixed individuals, we find significantly smaller marginal effect sizes for CFH risk alleles in African ancestry haplotypes. Broadening efforts to include ancestrally distinct populations helped uncover genes and pathways that boost risk in an ancestry-dependent manner and are potential targets for corrective therapies.

Genome-wide association study and polygenic risk score analysis for schizophrenia in a Korean population

Although large-scale genome-wide association studies (GWASs) have revealed the genetic architecture of schizophrenia, these studies have mainly focused on populations of European ancestry. This study aimed to identify common genetic variants associated with schizophrenia in the Korean population and evaluate the performance of polygenic risk scores (PRSs) derived from large-scale GWASs across ancestries. In the Korean psychiatric GWAS project (KPGP), seven academic institutes and their affiliated hospitals across South Korea recruited a cohort of 1670 patients with DSM-IV-defined schizophrenia and 2271 healthy controls. A total of 6690,822 SNPs were tested for association with schizophrenia. We identified one previously unreported genome-wide significant locus rs2423464 (P = 2.83 × 10-11; odds ratio = 1.65; 95 % confidence interval = 1.43-1.91, minor allele frequency = 0.126). This variant was also associated with increased lysosomal-associated membrane protein family member 5 (LAMP5) gene expression. The PRS derived from the meta-analysis results of East Asian and European GWASs explained a larger proportion of the phenotypic variance in the Korean schizophrenia sample than the PRS of an East Asian or European GWAS. (R2 = 0.073 for meta-analysis; 0.028 for East Asian GWAS; 0.037 for European GWAS). GWASs involving diverse ethnic groups will expand our understanding of the genetic architecture of schizophrenia.

Shared genetic architecture between COVID-19 and irritable bowel syndrome: a large-scale genome-wide cross-trait analysis

BACKGROUND

It has been reported that COVID-19 patients have an increased risk of developing IBS; however, the underlying genetic mechanisms of these associations remain largely unknown. The aim of this study was to investigate potential shared SNPs, genes, proteins, and biological pathways between COVID-19 and IBS by assessing pairwise genetic correlations and cross-trait genetic analysis.

MATERIALS AND METHODS

We assessed the genetic correlation between three COVID-19 phenotypes and IBS using linkage disequilibrium score regression (LDSC) and high-definition likelihood (HDL) methods. Two different sources of IBS data were combined using METAL, and the Multi-trait analysis of GWAS (MTAG) method was applied for multi-trait analysis to enhance statistical robustness and discover new genetic associations. Independent risk loci were examined using genome-wide complex trait analysis (GCTA)-conditional and joint analysis (COJO), multi-marker analysis of genomic annotation (MAGMA), and functional mapping and annotation (FUMA), integrating various QTL information and methods to further identify risk genes and proteins. Gene set variation analysis (GSVA) was employed to compute pleiotropic gene scores, and combined with immune infiltration algorithms, IBS patients were categorized into high and low immune infiltration groups.

RESULTS

We found a positive genetic correlation between COVID-19 infection, COVID-19 hospitalization, and IBS. Subsequent multi-trait analysis identified nine significantly associated genomic loci. Among these, eight genetic variants were closely related to the comorbidity of IBS and COVID-19. The study also highlighted four genes and 231 proteins associated with the susceptibility to IBS identified through various analytical strategies and a stratification approach for IBS risk populations.

CONCLUSIONS

Our study reveals a shared genetic architecture between these two diseases, providing new insights into potential biological mechanisms and laying the groundwork for more effective interventions.

Accounting for heterogeneity due to environmental sources in meta-analysis of genome-wide association studies

Meta-analysis of genome-wide association studies (GWAS) across diverse populations offers power gains to identify loci associated with complex traits and diseases. Often heterogeneity in effect sizes across populations will be correlated with genetic ancestry and environmental exposures (e.g. lifestyle factors). We present an environment-adjusted meta-regression model (env-MR-MEGA) to detect genetic associations by adjusting for and quantifying environmental and ancestral heterogeneity between populations. In simulations, env-MR-MEGA has similar or greater association power than MR-MEGA, with notable gains when the environmental factor has a greater correlation with the trait than ancestry. In our analysis of low-density lipoprotein cholesterol in ~19,000 individuals across twelve sex-stratified GWAS from Africa, adjusting for sex, BMI, and urban status, we identify additional heterogeneity beyond ancestral effects for seven variants. Env-MR-MEGA provides an approach to account for environmental effects using summary-level data, making it a useful tool for meta-analyses without the need to share individual-level data.

Identifying Common Genetic Etiologies Between Inflammatory Bowel Disease and Related Immune-Mediated Diseases

BACKGROUND

Patients with inflammatory bowel disease (IBD) have an increased risk of developing immune-mediated diseases. However, the genetic basis of IBD is complex, and an integrated approach should be used to elucidate the complex genetic relationship between IBD and immune-mediated diseases.

METHODS

The genetic relationship between IBD and 16 immune-mediated diseases was examined using linkage disequilibrium score regression. GWAS data were synthesized from two IBD databases using the METAL, and multi-trait analysis of genome-wide association studies was performed to enhance statistical robustness and identify novel genetic associations. Independent risk loci were meticulously examined using conditional and joint genome-wide multi-trait analysis, multi-marker analysis of genomic annotation, and functional mapping and annotation of significant genetic loci, integrating the information of quantitative trait loci and different methodologies to identify risk-related genes and proteins.

RESULTS

The results revealed four immune-mediated diseases (AS, psoriasis, iridocyclitis, and PsA) with a significant relationship with IBD. The multi-trait analysis revealed 909 gene loci of statistical significance. Of these loci, 28 genetic variants were closely related to IBD, and 7 single-nucleotide polymorphisms represented novel independent risk loci. In addition, 14 genes and 514 proteins were found to be associated with susceptibility to immune-mediated diseases. Notably, IL1RL1 emerged as a key player, present within pleiotropic genes across multiple protein databases, highlighting its potential as a therapeutic target.

CONCLUSIONS

This study suggests that the common polygenic determinants between IBD and immune-mediated diseases are widely distributed across the genome. The findings not only support a shared genetic relationship between IBD and immune-mediated diseases but also provide novel therapeutic targets for these diseases.

Genetic risk factors underlying white matter hyperintensities and cortical atrophy

White matter hyperintensities index structural abnormalities in the cerebral white matter, including axonal damage. The latter may promote atrophy of the cerebral cortex, a key feature of dementia. Here, we report a study of 51,065 individuals from 10 cohorts demonstrating that higher white matter hyperintensity volume associates with lower cortical thickness. The meta-GWAS of white matter hyperintensities-associated cortical 'atrophy' identifies 20 genome-wide significant loci, and enrichment in genes specific to vascular cell types, astrocytes, and oligodendrocytes. White matter hyperintensities-associated cortical 'atrophy' showed positive genetic correlations with vascular-risk traits and plasma biomarkers of neurodegeneration, and negative genetic correlations with cognitive functioning. 15 of the 20 loci regulated the expression of 54 genes in the cerebral cortex that, together with their co-expressed genes, were enriched in biological processes of axonal cytoskeleton and intracellular transport. The white matter hyperintensities-cortical thickness associations were most pronounced in cortical regions with higher expression of genes specific to excitatory neurons with long-range axons traversing through the white matter. The meta-GWAS-based polygenic risk score predicts vascular and all-cause dementia in an independent sample of 500,348 individuals. Thus, the genetics of white matter hyperintensities-related cortical atrophy involves vascular and neuronal processes and increases dementia risk.

Cross-ancestry meta-genome-wide association studies provide insights to the understanding of semen traits in pigs

Semen traits play a crucial role in pig reproduction and fertility. However, limited data availability hinder a comprehensive understanding of the genetic mechanisms underlying these traits. In this study, we integrated 597 299 ejaculates and 3 596 sequence data to identify genetic variants and candidate genes related to four semen traits, including sperm progressive motility (MOT), semen volume, sperm concentration (CON), and effective sperm count (SUM). A cross-ancestry meta-genome-wide association study was conducted to detect 163 lead single nucleotide polymorphisms (SNPs) associated with MOT, CON, and SUM. Subsequently, transcriptome-wide association studies and colocalisation analyses were integrated to identify 176 candidate genes, many of which have documented roles in spermatogenesis or male mammal semen traits. Our analysis highlighted the potential involvement of CSM5, PDZD9, and LDAF1 in regulating semen traits through multiple methods. Finally, to validate the function of significant SNPs, we performed genomic feature best linear unbiased prediction in 348 independent pigs using identified trait-related SNP subsets as genomic features. We found that integrating the top 0.1, 1, and 5% significant SNPs as genomic features could enhance genomic prediction accuracy for CON and MOT compared to traditional genomic best linear unbiased prediction. This study contributes to a comprehensive understanding of the genetic mechanisms of boar semen traits and provides insight for developing genomic selection models.

Single-Ancestry versus Multi-Ancestry Polygenic Risk Scores for CKD in Black American Populations

Key Points

The predictive performance of an African ancestry–specific polygenic risk score (PRS) was comparable to a European ancestry–derived PRS for kidney traits. However, multi-ancestry PRSs outperform single-ancestry PRSs in Black American populations. Predictive accuracy of PRSs for CKD was improved with the use of race-free eGFR.

Background

CKD is a risk factor of cardiovascular disease and early death. Recently, polygenic risk scores (PRSs) have been developed to quantify risk for CKD. However, African ancestry populations are underrepresented in both CKD genetic studies and PRS development overall. Moreover, European ancestry–derived PRSs demonstrate diminished predictive performance in African ancestry populations.

Methods

This study aimed to develop a PRS for CKD in Black American populations. We obtained score weights from a meta-analysis of genome-wide association studies for eGFR in the Million Veteran Program and Reasons for Geographic and Racial Differences in Stroke Study to develop an eGFR PRS. We optimized the PRS risk model in a cohort of participants from the Hypertension Genetic Epidemiology Network. Validation was performed in subsets of Black participants of the Trans-Omics in Precision Medicine Consortium and Genetics of Hypertension Associated Treatment Study.

Results

The prevalence of CKD—defined as stage 3 or higher—was associated with the PRS as a continuous predictor (odds ratio [95% confidence interval]: 1.35 [1.08 to 1.68]) and in a threshold-dependent manner. Furthermore, including APOL1 risk status—a putative variant for CKD with higher prevalence among those of sub-Saharan African descent—improved the score's accuracy. PRS associations were robust to sensitivity analyses accounting for traditional CKD risk factors, as well as CKD classification based on prior eGFR equations. Compared with previously published PRS, the predictive performance of our PRS was comparable with a European ancestry–derived PRS for kidney traits. However, single-ancestry PRSs were less predictive than multi-ancestry–derived PRSs.

Conclusions

In this study, we developed a PRS that was significantly associated with CKD with improved predictive accuracy when including APOL1 risk status. However, PRS generated from multi-ancestry populations outperformed single-ancestry PRS in our study.

An integrative systems-biology approach defines mechanisms of Alzheimer's disease neurodegeneration

Despite years of intense investigation, the mechanisms underlying neuronal death in Alzheimer's disease, the most common neurodegenerative disorder, remain incompletely understood. To define relevant pathways, we integrated the results of an unbiased, genome-scale forward genetic screen for age-associated neurodegeneration in Drosophila with human and Drosophila Alzheimer's disease-associated multi-omics. We measured proteomics, phosphoproteomics, and metabolomics in Drosophila models of Alzheimer's disease and identified Alzheimer's disease human genetic variants that modify expression in disease-vulnerable neurons. We used a network optimization approach to integrate these data with previously published Alzheimer's disease multi-omic data. We computationally predicted and experimentally demonstrated how HNRNPA2B1 and MEPCE enhance tau-mediated neurotoxicity. Furthermore, we demonstrated that the screen hits CSNK2A1 and NOTCH1 regulate DNA damage in Drosophila and human iPSC-derived neural progenitor cells. Our work identifies candidate pathways that could be targeted to ameliorate neurodegeneration in Alzheimer's disease.

Multi-ancestry GWAS meta-analyses of lung cancer reveal susceptibility loci and elucidate smoking-independent genetic risk

Lung cancer remains the leading cause of cancer mortality, despite declining smoking rates. Previous lung cancer GWAS have identified numerous loci, but separating the genetic risks of lung cancer and smoking behavioral susceptibility remains challenging. Here, we perform multi-ancestry GWAS meta-analyses of lung cancer using the Million Veteran Program cohort (approximately 95% male cases) and a previous study of European-ancestry individuals, jointly comprising 42,102 cases and 181,270 controls, followed by replication in an independent cohort of 19,404 cases and 17,378 controls. We then carry out conditional meta-analyses on cigarettes per day and identify two novel, replicated loci, including the 19p13.11 pleiotropic cancer locus in squamous cell lung carcinoma. Overall, we report twelve novel risk loci for overall lung cancer, lung adenocarcinoma, and squamous cell lung carcinoma, nine of which are externally replicated. Finally, we perform PheWAS on polygenic risk scores for lung cancer, with and without conditioning on smoking. The unconditioned lung cancer polygenic risk score is associated with smoking status in controls, illustrating a reduced predictive utility in non-smokers. Additionally, our polygenic risk score demonstrates smoking-independent pleiotropy of lung cancer risk across neoplasms and metabolic traits.

A comprehensive framework for trans-ancestry pathway analysis using GWAS summary data from diverse populations

As more multi-ancestry GWAS summary data become available, we have developed a comprehensive trans-ancestry pathway analysis framework that effectively utilizes this diverse genetic information. Within this framework, we evaluated various strategies for integrating genetic data at different levels-SNP, gene, and pathway-from multiple ancestry groups. Through extensive simulation studies, we have identified robust strategies that demonstrate superior performance across diverse scenarios. Applying these methods, we analyzed 6,970 pathways for their association with schizophrenia, incorporating data from African, East Asian, and European populations. Our analysis identified over 200 pathways significantly associated with schizophrenia, even after excluding genes near genome-wide significant loci. This approach substantially enhances detection efficiency compared to traditional single-ancestry pathway analysis and the conventional approach that amalgamates single-ancestry pathway analysis results across different ancestry groups. Our framework provides a flexible and effective tool for leveraging the expanding pool of multi-ancestry GWAS summary data, thereby improving our ability to identify biologically relevant pathways that contribute to disease susceptibility.

Micklesfield, Ramsay M, Sarah P, Peter J. A GWAS of ACE Inhibitor-Induced Angioedema in a South African Population

Background

Angiotensin-converting enzyme inhibitor-induced angioedema (AE-ACEI) is a life-threatening adverse event and, globally, the commonest cause of emergency presentations with angioedema. Several large genome-wide association studies (GWAS) have found genomic associations with AE-ACEI. However, despite African Americans having a 5-fold increased risk of AE-ACEI, there are no published GWAS from Africa. The aim of this study was to conduct a case-control GWAS of AE-ACEI in a South African population and perform a meta-analysis with an African American and European American population.

Methods

The GWAS included 202 South African adults with a history of AE-ACEI and 513 controls without angioedema following angiotensin-converting enzyme inhibitor (ACEI) treatment for at least 2 years. A meta-analysis was conducted with GWAS summary statistics from an African American and European American cohort (from Vanderbilt/Marshfield with 174 cases and 489 controls).

Results

No SNPs attained genome-wide significance. However, 26 SNPs in the post-imputation standard GWAS of the South African cohort and 37 SNPs in the meta-analysis were associated to AE-ACEI with suggestive threshold(p-value<5.0×10-06). Some of these SNPs were found to be located close to the genes PRKCQ and RIMS1, previously linked with drug-induced angioedema, and also close to the CSMD1 gene linked to ACEI cough, providing replication at the gene level, but with novel lead SNPs.

Conclusions

Our results highlight the importance of African populations to detect novel variants in replication studies. Further increased sampling across the continent and matched functional work are needed to confirm the importance of genetic variation in understanding the biology of AE-ACEI.

Genome-wide association analysis and admixture mapping in a Puerto Rican cohort supports an Alzheimer disease risk locus on chromosome 12

Introduction

Hispanic/Latino populations are underrepresented in Alzheimer Disease (AD) genetic studies. Puerto Ricans (PR), a three-way admixed (European, African, and Amerindian) population is the second-largest Hispanic group in the continental US. We aimed to conduct a genome-wide association study (GWAS) and comprehensive analyses to identify novel AD susceptibility loci and characterize known AD genetic risk loci in the PR population.

Materials and methods

Our study included Whole Genome Sequencing (WGS) and phenotype data from 648 PR individuals (345 AD, 303 cognitively unimpaired). We used a generalized linear-mixed model adjusting for sex, age, population substructure, and genetic relationship matrix. To infer local ancestry, we merged the dataset with the HGDP/1000G reference panel. Subsequently, we conducted univariate admixture mapping (AM) analysis.

Results

We identified suggestive signals within the SLC38A1 and SCN8A genes on chromosome 12q13. This region overlaps with an area of linkage of AD in previous studies (12q13) in independent data sets further supporting. Univariate African AM analysis identified one suggestive ancestral block (p = 7.2×10-6) located in the same region. The ancestry-aware approach showed that this region has both European and African ancestral backgrounds and both contributing to the risk in this region. We also replicated 11 different known AD loci -including APOE- identified in mostly European studies, which is likely due to the high European background of the PR population.

Conclusion

PR GWAS and AM analysis identified a suggestive AD risk locus on chromosome 12, which includes the SLC38A1 and SCN8A genes. Our findings demonstrate the importance of designing GWAS and ancestry-aware approaches and including underrepresented populations in genetic studies of AD.

Fine-mapping across diverse ancestries drives the discovery of putative causal variants underlying human complex traits and diseases

Genome-wide association studies (GWAS) of human complex traits or diseases often implicate genetic loci that span hundreds or thousands of genetic variants, many of which have similar statistical significance. While statistical fine-mapping in individuals of European ancestry has made important discoveries, cross-population fine-mapping has the potential to improve power and resolution by capitalizing on the genomic diversity across ancestries. Here we present SuSiEx, an accurate and computationally efficient method for cross-population fine-mapping. SuSiEx integrates data from an arbitrary number of ancestries, explicitly models population-specific allele frequencies and linkage disequilibrium patterns, accounts for multiple causal variants in a genomic region and can be applied to GWAS summary statistics. We comprehensively assessed the performance of SuSiEx using simulations. We further showed that SuSiEx improves the fine-mapping of a range of quantitative traits available in both the UK Biobank and Taiwan Biobank, and improves the fine-mapping of schizophrenia-associated loci by integrating GWAS across East Asian and European ancestries.

Hierarchical joint analysis of marginal summary statistics-Part I: Multipopulation fine mapping and credible set construction

Recent advancement in genome-wide association studies (GWAS) comes from not only increasingly larger sample sizes but also the shift in focus towards underrepresented populations. Multipopulation GWAS increase power to detect novel risk variants and improve fine-mapping resolution by leveraging evidence and differences in linkage disequilibrium (LD) from diverse populations. Here, we expand upon our previous approach for single-population fine-mapping through Joint Analysis of Marginal SNP Effects (JAM) to a multipopulation analysis (mJAM). Under the assumption that true causal variants are common across studies, we implement a hierarchical model framework that conditions on multiple SNPs while explicitly incorporating the different LD structures across populations. The mJAM framework can be used to first select index variants using the mJAM likelihood with different feature selection approaches. In addition, we present a novel approach leveraging the ideas of mediation to construct credible sets for these index variants. Construction of such credible sets can be performed given any existing index variants. We illustrate the implementation of the mJAM likelihood through two implementations: mJAM-SuSiE (a Bayesian approach) and mJAM-Forward selection. Through simulation studies based on realistic effect sizes and levels of LD, we demonstrated that mJAM performs well for constructing concise credible sets that include the underlying causal variants. In real data examples taken from the most recent multipopulation prostate cancer GWAS, we showed several practical advantages of mJAM over other existing multipopulation methods.

Biological Insights from Schizophrenia-associated Loci in Ancestral Populations

Large-scale genome-wide association studies of schizophrenia have uncovered hundreds of associated loci but with extremely limited representation of African diaspora populations. We surveyed electronic health records of 200,000 individuals of African ancestry in the Million Veteran and All of Us Research Programs, and, coupled with genotype-level data from four case-control studies, realized a combined sample size of 13,012 affected and 54,266 unaffected persons. Three genome-wide significant signals - near PLXNA4, PMAIP1, and TRPA1 - are the first to be independently identified in populations of predominantly African ancestry. Joint analyses of African, European, and East Asian ancestries across 86,981 cases and 303,771 controls, yielded 376 distinct autosomal loci, which were refined to 708 putatively causal variants via multi-ancestry fine-mapping. Utilizing single-cell functional genomic data from human brain tissue and two complementary approaches, transcriptome-wide association studies and enhancer-promoter contact mapping, we identified a consensus set of 94 genes across ancestries and pinpointed the specific cell types in which they act. We identified reproducible associations of schizophrenia polygenic risk scores with schizophrenia diagnoses and a range of other mental and physical health problems. Our study addresses a longstanding gap in the generalizability of research findings for schizophrenia across ancestral populations, underlining shared biological underpinnings of schizophrenia across global populations in the presence of broadly divergent risk allele frequencies.

Ancestry-aligned polygenic scores combined with conventional risk factors improve prediction of cardiometabolic outcomes in African populations

BACKGROUND

Cardiovascular diseases (CVD) are a major health concern in Africa. Improved identification and treatment of high-risk individuals can reduce adverse health outcomes. Current CVD risk calculators are largely unvalidated in African populations and overlook genetic factors. Polygenic scores (PGS) can enhance risk prediction by measuring genetic susceptibility to CVD, but their effectiveness in genetically diverse populations is limited by a European-ancestry bias. To address this, we developed models integrating genetic data and conventional risk factors to assess the risk of developing cardiometabolic outcomes in African populations.

METHODS

We used summary statistics from a genome-wide association meta-analysis (n = 14,126) in African populations to derive novel genome-wide PGS for 14 cardiometabolic traits in an independent African target sample (Africa Wits-INDEPTH Partnership for Genomic Research (AWI-Gen), n = 10,603). Regression analyses assessed relationships between each PGS and corresponding cardiometabolic trait, and seven CVD outcomes (CVD, heart attack, stroke, diabetes mellitus, dyslipidaemia, hypertension, and obesity). The predictive utility of the genetic data was evaluated using elastic net models containing multiple PGS (MultiPGS) and reference-projected principal components of ancestry (PPCs). An integrated risk prediction model incorporating genetic and conventional risk factors was developed. Nested cross-validation was used when deriving elastic net models to enhance generalisability.

RESULTS

Our African-specific PGS displayed significant but variable within- and cross- trait prediction (max.R2 = 6.8%, p = 1.86 × 10-173). Significantly associated PGS with dyslipidaemia included the PGS for total cholesterol (logOR = 0.210, SE = 0.022, p = 2.18 × 10-21) and low-density lipoprotein (logOR =  - 0.141, SE = 0.022, p = 1.30 × 10-20); with hypertension, the systolic blood pressure PGS (logOR = 0.150, SE = 0.045, p = 8.34 × 10-4); and multiple PGS associated with obesity: body mass index (max. logOR = 0.131, SE = 0.031, p = 2.22 × 10-5), hip circumference (logOR = 0.122, SE = 0.029, p = 2.28 × 10-5), waist circumference (logOR = 0.013, SE = 0.098, p = 8.13 × 10-4) and weight (logOR = 0.103, SE = 0.029, p = 4.89 × 10-5). Elastic net models incorporating MultiPGS and PPCs significantly improved prediction over MultiPGS alone. Models including genetic data and conventional risk factors were more predictive than conventional risk models alone (dyslipidaemia: R2 increase = 2.6%, p = 4.45 × 10-12; hypertension: R2 increase = 2.6%, p = 2.37 × 10-13; obesity: R2 increase = 5.5%, 1.33 × 10-34).

CONCLUSIONS

In African populations, CVD and associated cardiometabolic trait prediction models can be improved by incorporating ancestry-aligned PGS and accounting for ancestry. Combining PGS with conventional risk factors further enhances prediction over traditional models based on conventional factors. Incorporating data from target populations can improve the generalisability of international predictive models for CVD and associated traits in African populations.

Structural variants linked to Alzheimer's Disease and other common age-related clinical and neuropathologic traits

Advances have led to a greater understanding of the genetics of Alzheimer's Disease (AD). However, the gap between the predicted and observed genetic heritability estimates when using single nucleotide polymorphisms (SNPs) and small indel data remains. Large genomic rearrangements, known as structural variants (SVs), have the potential to account for this missing genetic heritability. By leveraging data from two ongoing cohort studies of aging and dementia, the Religious Orders Study and Rush Memory and Aging Project (ROS/MAP), we performed genome-wide association analysis testing around 20,000 common SVs from 1,088 participants with whole genome sequencing (WGS) data. A range of Alzheimer's Disease and Related Disorders (AD/ADRD) clinical and pathologic traits were examined. Given the limited sample size, no genome-wide significant association was found, but we mapped SVs across 81 AD risk loci and discovered 22 SVs in linkage disequilibrium (LD) with GWAS lead variants and directly associated with AD/ADRD phenotypes (nominal P < 0.05). The strongest association was a deletion of an Alu element in the 3'UTR of the TMEM106B gene. This SV was in high LD with the respective AD GWAS locus and was associated with multiple AD/ADRD phenotypes, including tangle density, TDP-43, and cognitive resilience. The deletion of this element was also linked to lower TMEM106B protein abundance. We also found a 22 kb deletion associated with depression in ROSMAP and bearing similar association patterns as AD GWAS SNPs at the IQCK locus. In addition, genome-wide scans allowed the identification of 7 SVs, with no LD with SNPs and nominally associated with AD/ADRD traits. This result suggests potentially new ADRD risk loci not discoverable using SNP data. Among these findings, we highlight a 5.6 kb duplication of coding regions of the gene C1orf186 at chromosome 1 associated with indices of cognitive impairment, decline, and resilience. While further replication in independent datasets is needed to validate these findings, our results support the potential roles of common structural variations in the pathogenesis of AD/ADRD.

Genome-Wide Association Study Meta-Analysis Uncovers Novel Genetic Variants Associated with Olfactory Dysfunction

IMPORTANCE

Olfactory dysfunction is among the earliest signs of many age-related neurodegenerative diseases and has been associated with increased mortality in older adults; however, its genetic basis remains largely unknown.

OBJECTIVE

To identify the genetic loci associated with olfactory dysfunction in the general population.

DESIGN SETTING AND PARTIICIPANTS

This genome-wide association study meta-analysis (GWMA) included participants of European ancestry (N = 22,730) enrolled in four different large population-based studies, followed by a multi-ancestry GWMA including participants of African ancestry (N = 1,030). The data analysis was performed from March 2023 through June 2024.

EXPOSURES

Genome-wide single nucleotide polymorphisms.

MAIN OUTCOMES AND MEASURES

Olfactory dysfunction was the outcome and assessed using a 12-item smell identification test.

RESULTS

GWMA revealed a novel genome-wide significant locus (tagged by rs11228623 at 11q12) associated with olfactory dysfunction. Gene-based analysis revealed a high enrichment for olfactory receptor genes in this region. Phenome-wide association studies demonstrated associations between genetic variants related to olfactory dysfunction and blood cell counts, kidney function, skeletal muscle mass, cholesterol levels and cardiovascular disease. Using individual-level data, we also confirmed and quantified the strength of these associations on a phenotypic level. Moreover, employing two-sample Mendelian Randomization analyses, we found evidence for causal associations between olfactory dysfunction and these phenotypes.

CONCLUSIONS

These findings provide novel insights into the genetic architecture of the sense of smell and highlight its importance for many aspects of human health.

Addressing overlapping sample challenges in genome-wide association studies: Meta-reductive approach

Polygenic risk scores (PRS) are instrumental in genetics, offering insights into an individual level genetic risk to a range of diseases based on accumulated genetic variations. These scores rely on Genome-Wide Association Studies (GWAS). However, precision in PRS is often challenged by the requirement of extensive sample sizes and the potential for overlapping datasets that can inflate PRS calculations. In this study, we present a novel methodology, Meta-Reductive Approach (MRA), that was derived algebraically to adjust GWAS results, aiming to neutralize the influence of select cohorts. Our approach recalibrates summary statistics using algebraic derivations. Validating our technique with datasets from Alzheimer disease studies, we showed that the summary statistics of the MRA and those derived from individual-level data yielded the exact same values. This innovative method offers a promising avenue for enhancing the accuracy of PRS, especially when derived from meta-analyzed GWAS data.

Methods for multiancestry genome-wide association study meta-analysis

Genome-wide association studies (GWAS) have significantly enhanced our understanding of the genetic basis of complex diseases. Despite the technological advancements, gaps in our understanding remain, partly due to small effect sizes and inadequate coverage of genetic variation. Multiancestry GWAS meta-analysis (MAGMA) addresses these challenges by integrating genetic data from diverse populations, thereby increasing power to detect loci and improving fine-mapping resolution to identify causal variants across different ancestry groups. This review provides an overview of the protocols, statistical methods, and software of MAGMA, as well as highlighting some challenges associated with this approach.

Genome-Wide Association Study Points to Novel Locus for Gilles de la Tourette Syndrome

BACKGROUND

Tourette syndrome (TS) is a childhood-onset neurodevelopmental disorder of complex genetic architecture and is characterized by multiple motor tics and at least one vocal tic persisting for more than 1 year.

METHODS

We performed a genome-wide meta-analysis integrating a novel TS cohort with previously published data, resulting in a sample size of 6133 individuals with TS and 13,565 ancestry-matched control participants.

RESULTS

We identified a genome-wide significant locus on chromosome 5q15. Integration of expression quantitative trait locus, Hi-C (high-throughput chromosome conformation capture), and genome-wide association study data implicated the NR2F1 gene and associated long noncoding RNAs within the 5q15 locus. Heritability partitioning identified statistically significant enrichment in brain tissue histone marks, while polygenic risk scoring of brain volume data identified statistically significant associations with right and left thalamus volumes and right putamen volume.

CONCLUSIONS

Our work presents novel insights into the neurobiology of TS, thereby opening up new directions for future studies.

Shared molecular mechanisms and transdiagnostic potential of neurodevelopmental disorders and immune disorders

The co-occurrence and familial clustering of neurodevelopmental disorders and immune disorders suggest shared genetic risk factors. Based on genome-wide association summary statistics from five neurodevelopmental disorders and four immune disorders, we conducted genome-wide, local genetic correlation and polygenic overlap analysis. We further performed a cross-trait GWAS meta-analysis. Pleotropic loci shared between the two categories of diseases were mapped to candidate genes using multiple algorithms and approaches. Significant genetic correlations were observed between neurodevelopmental disorders and immune disorders, including both positive and negative correlations. Neurodevelopmental disorders exhibited higher polygenicity compared to immune disorders. Around 50%-90% of genetic variants of the immune disorders were shared with neurodevelopmental disorders. The cross-trait meta-analysis revealed 154 genome-wide significant loci, including 8 novel pleiotropic loci. Significant associations were observed for 30 loci with both types of diseases. Pathway analysis on the candidate genes at these loci revealed common pathways shared by the two types of diseases, including neural signaling, inflammatory response, and PI3K-Akt signaling pathway. In addition, 26 of the 30 lead SNPs were associated with blood cell traits. Neurodevelopmental disorders exhibit complex polygenic architecture, with a subset of individuals being at a heightened genetic risk for both neurodevelopmental and immune disorders. The identification of pleiotropic loci has important implications for exploring opportunities for drug repurposing, enabling more accurate patient stratification, and advancing genomics-informed precision in the medical field of neurodevelopmental disorders.