A subset-based approach improves power and interpretation for the combined analysis of genetic association studies of heterogeneous traits

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.

Genome-wide association analysis reveals the associations of NPHP4, TYW1-AUTS2 and SEMA6D for Behçet's disease and HLA-B*46:01 for its intestinal involvement

BACKGROUND

Intestinal involvement in Behçet's disease (BD) is associated with poor prognosis and is more prevalent in East Asian than in Mediterranean populations. Identifying the genetic causes of intestinal BD is important for understanding the pathogenesis and for appropriate treatment of BD patients.

METHODS

We performed genome-wide association studies (GWAS) and imputation/replication genotyping of human leukocyte antigen (HLA) alleles for 1,689 Korean and Turkish patients with BD (including 379 patients with intestinal BD) and 2,327 healthy controls, followed by replication using 593 Japanese patients with BD (101 patients with intestinal BD) and 737 healthy controls. Stratified cross-phenotype analyses were performed for 1) overall BD, 2) intestinal BD, and 3) intestinal BD without association of overall BD.

RESULTS

We identified three novel genome-wide significant susceptibility loci including NPHP4 (rs74566205; P=1.36 × 10-8), TYW1-AUTS2 (rs60021986; P=1.14 × 10-9), and SEMA6D (rs4143322; P=5.54 × 10-9) for overall BD, and a new association with HLA-B*46:01 for intestinal BD (P=1.67 × 10-8) but not for BD without intestinal involvement. HLA peptide binding analysis revealed that Mycobacterial peptides, have a stronger binding affinity to HLA-B*46:01 compared to the known risk allele HLA-B*51:01.

CONCLUSIONS

HLA-B*46:01 is associated with the development of intestinal BD; NPHP4, TYW1-AUTS2, and SEMA6D are susceptibility loci for overall BD.

Combining Transdiagnostic and Disorder-Level GWAS Enhances Precision of Psychiatric Genetic Risk Profiles in a Multi-Ancestry Sample

The etiology of substance use disorders (SUDs) and psychiatric disorders reflects a combination of both transdiagnostic (i.e., common) and disorder-level (i.e., independent) genetic risk factors. We applied genomic structural equation modeling to examine these genetic factors across SUDs, psychotic, mood, and anxiety disorders using genome-wide association studies (GWAS) of European- (EUR) and African-ancestry (AFR) individuals. In EUR individuals, transdiagnostic genetic factors represented SUDs (143 lead single nucleotide polymorphisms [SNPs]), psychotic (162 lead SNPs), and mood/anxiety disorders (112 lead SNPs). We identified two novel SNPs for mood/anxiety disorders that have probable regulatory roles on FOXP1, NECTIN3, and BTLA genes. In AFR individuals, genetic factors represented SUDs (1 lead SNP) and psychiatric disorders (no significant SNPs). The SUD factor lead SNP, although previously significant in EUR- and cross-ancestry GWAS, is a novel finding in AFR individuals. Shared genetic variance accounted for overlap between SUDs and their psychiatric comorbidities, with second-order GWAS identifying up to 12 SNPs not significantly associated with either first-order factor in EUR individuals. Finally, common and independent genetic effects showed different associations with psychiatric, sociodemographic, and medical phenotypes. For example, the independent components of schizophrenia and bipolar disorder had distinct associations with affective and risk-taking behaviors, and phenome-wide association studies identified medical conditions associated with tobacco use disorder independent of the broader SUDs factor. Thus, combining transdiagnostic and disorder-level genetic approaches can improve our understanding of co-occurring conditions and increase the specificity of genetic discovery, which is critical for psychiatric disorders that demonstrate considerable symptom and etiological overlap.

Dissecting Schizophrenia Biology Using Pleiotropy with Cognitive Genomics

Given the increasingly large number of loci discovered by psychiatric GWAS, specification of the key biological pathways underlying these loci has become a priority for the field. We have previously leveraged the pleiotropic genetic relationships between schizophrenia and two cognitive phenotypes (educational attainment and cognitive task performance) to differentiate two subsets of illness-relevant SNPs: (1) those with "concordant" alleles, which are associated with reduced cognitive ability/education and increased schizophrenia risk; and (2) those with "discordant" alleles linked to reduced educational and/or cognitive levels but lower schizophrenia susceptibility. In the present study, we extend our prior work, utilizing larger input GWAS datasets and a more powerful statistical approach to pleiotropic meta-analysis, the Pleiotropic Locus Exploration and Interpretation using Optimal test (PLEIO). Our pleiotropic meta-analysis of schizophrenia and the two cognitive phenotypes revealed 768 significant loci (159 novel). Among these, 347 loci harbored concordant SNPs, 270 encompassed discordant SNPs, and 151 "dual" loci contained concordant and discordant SNPs. Competitive gene-set analysis using MAGMA related concordant SNP loci with neurodevelopmental pathways (e.g., neurogenesis), whereas discordant loci were associated with mature neuronal synaptic functions. These distinctions were also observed in BrainSpan analysis of temporal enrichment patterns across developmental periods, with concordant loci containing more prenatally expressed genes than discordant loci. Dual loci were enriched for genes related to mRNA translation initiation, representing a novel finding in the schizophrenia literature.

Subset-based method for cross-tissue transcriptome-wide association studies improves power and interpretability

Integrating results from genome-wide association studies (GWASs) and studies of molecular phenotypes such as gene expressions can improve our understanding of the biological functions of trait-associated variants and can help prioritize candidate genes for downstream analysis. Using reference expression quantitative trait locus (eQTL) studies, several methods have been proposed to identify gene-trait associations, primarily based on gene expression imputation. To increase the statistical power by leveraging substantial eQTL sharing across tissues, meta-analysis methods aggregating such gene-based test results across multiple tissues or contexts have been developed as well. However, most existing meta-analysis methods have limited power to identify associations when the gene has weaker associations in only a few tissues and cannot identify the subset of tissues in which the gene is "activated." For this, we developed a cross-tissue subset-based transcriptome-wide association study (CSTWAS) meta-analysis method that improves power under such scenarios and can extract the set of potentially associated tissues. To improve applicability, CSTWAS uses only GWAS summary statistics and pre-computed correlation matrices to identify a subset of tissues that have the maximal evidence of gene-trait association. Through numerical simulations, we found that CSTWAS can maintain a well-calibrated type-I error rate, improves power especially when there is a small number of associated tissues for a gene-trait association, and identifies an accurate associated tissue set. By analyzing GWAS summary statistics of three complex traits and diseases, we demonstrate that CSTWAS could identify biological meaningful signals while providing an interpretation of disease etiology by extracting a set of potentially associated tissues.

Characterizing the polygenic overlap and shared loci between rheumatoid arthritis and cardiovascular diseases

BACKGROUND

Despite substantial research revealing that patients with rheumatoid arthritis (RA) have excessive morbidity and mortality of cardiovascular disease (CVD), the mechanism underlying this association has not been fully known. This study aims to systematically investigate the phenotypic and genetic correlation between RA and CVD.

METHODS

Based on UK Biobank, we conducted two cohort studies to evaluate the phenotypic relationships between RA and CVD, including atrial fibrillation (AF), coronary artery disease (CAD), heart failure (HF), and stroke. Next, we used linkage disequilibrium score regression, Local Analysis of [co]Variant Association, and bivariate causal mixture model (MiXeR) methods to examine the genetic correlation and polygenic overlap between RA and CVD, using genome-wide association summary statistics. Furthermore, we explored specific shared genetic loci by conjunctional false discovery rate analysis and association analysis based on subsets.

RESULTS

Compared with the general population, RA patients showed a higher incidence of CVD (hazard ratio [HR] = 1.21, 95% confidence interval [CI]: 1.15-1.28). We observed positive genetic correlations of RA with AF and stroke, and a mixture of negative and positive local genetic correlations underlying the global genetic correlation for CAD and HF, with 13 ~ 33% of shared genetic variants for these trait pairs. We further identified 23 pleiotropic loci associated with RA and at least one CVD, including one novel locus (rs7098414, TSPAN14, 10q23.1). Genes mapped to these shared loci were enriched in immune and inflammatory-related pathways, and modifiable risk factors, such as high diastolic blood pressure.

CONCLUSIONS

This study revealed the shared genetic architecture of RA and CVD, which may facilitate drug target identification and improved clinical management.

A cross-trait study of lung cancer and its related respiratory diseases based on large-scale exome sequencing population

Background

Genome-wide association studies (GWASs) explain the genetic susceptibility between diseases and common variants. Nevertheless, with the appearance of large-scale sequencing profiles, we could explore the rare coding variants in disease pathogenesis.

Methods

We estimated the genetic correlation of nine respiratory diseases and lung cancer in the UK Biobank (UKB) by linkage disequilibrium score regression (LDSC). Then, we performed exome-wide association studies at single-variant level and gene-level for lung cancer and lung cancer-related respiratory diseases using the whole-exome sequencing (WES) data of 427,934 European participants. Cross-trait meta-analysis was conducted by association analysis based on subsets (ASSET) to identify the pleiotropic variants, while in-silico functional analysis was performed to explore their function. Causal mediation analysis was used to explore whether these pleiotropic variants lead to lung cancer is mediated by affecting the chronic respiratory diseases.

Results

Five respiratory diseases [emphysema, pneumonia, asthma, chronic obstructive pulmonary disease (COPD), and fibrosis] were genetically correlated with lung cancer. We identified 102 significant independent variants at single-variant levels for lung cancer and five lung cancer-related diseases. 15:78590583:G>A (missense variant in CHRNA5) was shared in lung cancer, emphysema, and COPD. Meanwhile, 14 significant genes and 87 suggestive genes were identified in gene-based association tests, including HSD3B7 (lung cancer), SRSF2 (pneumonia), TNXB (asthma), TERT (fibrosis), MOSPD3 (emphysema). Based on the cross-trait meta-analysis, we detected 145 independent pleiotropic variants. We further identified abundant pathways with significant enrichment effects, demonstrating that these pleiotropic genes were functional. Meanwhile, the proportion of mediation effects of these variants ranged from 6 to 23 (emphysema: 23%; COPD: 20%; pneumonia: 20%; fibrosis: 7%; asthma: 6%) through these five respiratory diseases to the incidence of lung cancer.

Conclusions

The identified shared genetic variants, genes, biological pathways, and potential intermediate causal pathways provide a basis for further exploration of the relationship between lung cancer and respiratory diseases.

Identifying risk loci for obsessive-compulsive disorder and shared genetic component with schizophrenia: A large-scale multi-trait association analysis with summary statistics

Due to limited samples, no genetic loci have been identified for obsessive-compulsive disorder (OCD) in genome-wide association studies. Additionally, although co-morbidities between OCD and schizophrenia (SCZ) were observed, their common genetic etiology was not completely known. Here, we conducted a comprehensive investigation regarding the genetic architecture of OCD and the common genetic foundation shared by OCD and SCZ using summary statistics data (2688 cases and 7037 controls for OCD; 53,386 cases and 77,258 controls for SCZ). We discovered significant genetic correlation between OCD and SCZ (r̂g=0.296, P = 2.82 × 10-11). We then performed two multi-trait association analyses to detect OCD-associated loci and colocalization analysis to detect causal variants. Parallel gene-level analyses were also implemented. We identified 323 OCD-relevant variants located within 12 loci, with four loci shared the same causal variants between OCD and SCZ. Further, the gene-level analyses discovered 8 OCD-associated genes. Finally, multiple functional analyses at both SNP and gene levels showed that these genetic association signals had significant enrichments in the regions of left ventricle and anterior cingulate cortex, and suggested an important role of pathways involving regulation of telomere maintenance, histone phosphorylation, and GnRH secretion. Overall, this study identified new genetic loci for OCD and provided substantial evidence supporting common genetic foundation underlying OCD and SCZ. The findings advanced our understanding of genetic architecture and pathophysiology of OCD as well as shedding light on shared genetic etiology of the two disorders.

Identifying risk loci for FTD and shared genetic component with ALS: A large-scale multitrait association analysis

Current genome-wide association studies of frontotemporal dementia (FTD) are underpowered due to limited samples. Further, common genetic etiologies between FTD and amyotrophic lateral sclerosis (ALS) remain unknown. Using the largest summary statistics of FTD (3526 cases and 9402 controls) and ALS (27,205 cases and 110,881 controls), we found a significant genetic correlation between them (rˆg = 0.637, P = 0.032) and identified 190 FTD-related variants within 5 loci (3p22.1, 5q35.1, 9p21.2, 19p13.11, and 20q13.13). Among these, ALS and FTD had causal variants in 9p21.2 and 19p13.11. Moreover, MOBP (3p22.1), C9orf72 (9p21.2), MOB3B (9p21.2), UNC13A (19p13.11), SLC9A8 (20q13.13), SNAI1 (20q13.13), and SPATA2 (20q13.13) were discovered by both SNP- and gene-level analyses, which together discovered 15 FTD-associated genes, with 10 not detected before (IFNK, RNF114, SLC9A8, SPATA2, SNAI1, SCFD1, POLDIP2, TMEM97, G2E3, and PIGW). Functional analyses showed these genes were enriched in heart left ventricle, kidney cortex, and some brain regions. Overall, this study provides insights into genetic determinants of FTD and shared genetic etiology underlying FTD and ALS.

Cross-ancestry genetic investigation of schizophrenia, cannabis use disorder, and tobacco smoking

Individuals with schizophrenia frequently experience co-occurring substance use, including tobacco smoking and heavy cannabis use, and substance use disorders. There is interest in understanding the extent to which these relationships are causal, and to what extent shared genetic factors play a role. We explored the relationships between schizophrenia (Scz), cannabis use disorder (CanUD), and ever-regular tobacco smoking (Smk) using the largest available genome-wide studies of these phenotypes in individuals of African and European ancestries. All three phenotypes were positively genetically correlated (rgs = 0.17 - 0.62). Causal inference analyses suggested the presence of horizontal pleiotropy, but evidence for bidirectional causal relationships was also found between all three phenotypes even after correcting for horizontal pleiotropy. We identified 439 pleiotropic loci in the European ancestry data, 150 of which were novel (i.e., not genome-wide significant in the original studies). Of these pleiotropic loci, 202 had lead variants which showed convergent effects (i.e., same direction of effect) on Scz, CanUD, and Smk. Genetic variants convergent across all three phenotypes showed strong genetic correlations with risk-taking, executive function, and several mental health conditions. Our results suggest that both horizontal pleiotropy and causal mechanisms may play a role in the relationship between CanUD, Smk, and Scz, but longitudinal, prospective studies are needed to confirm a causal relationship.

PASTRY: achieving balanced power for detecting risk and protective minor alleles in meta-analysis of association studies with overlapping subjects

BACKGROUND

Meta-analysis is a statistical method that combines the results of multiple studies to increase statistical power. When multiple studies participating in a meta-analysis utilize the same public dataset as controls, the summary statistics from these studies become correlated. To solve this challenge, Lin and Sullivan proposed a method to provide an optimal test statistic adjusted for the correlation. This method quickly became the standard practice. However, we identified an unexpected power asymmetry phenomenon in this standard framework. This can lead to unbalanced power for detecting protective minor alleles and risk minor alleles.

RESULTS

We found that the power asymmetry of the current framework is mainly due to the errors in approximating the correlation term. We then developed a meta-analysis method based on an accurate correlation estimator, called PASTRY (A method to avoid Power ASymmeTRY). PASTRY outperformed the standard method on both simulated and real datasets in terms of the power symmetry.

CONCLUSIONS

Our findings suggest that PASTRY can help to alleviate the power asymmetry problem. PASTRY is available at https://github.com/hanlab-SNU/PASTRY .

mixWAS: An efficient distributed algorithm for mixed-outcomes genome-wide association studies

Genome-wide association studies (GWAS) have been instrumental in identifying genetic associations for various diseases and traits. However, uncovering genetic underpinnings among traits beyond univariate phenotype associations remains a challenge. Multi-phenotype associations (MPA), or genetic pleiotropy, offer important insights into shared genes and pathways among traits, enhancing our understanding of genetic architectures of complex diseases. GWAS of biobank-linked electronic health record (EHR) data are increasingly being utilized to identify MPA among various traits and diseases. However, methodologies that can efficiently take advantage of distributed EHR to detect MPA are still lacking. Here, we introduce mixWAS, a novel algorithm that efficiently and losslessly integrates multiple EHRs via summary statistics, allowing the detection of MPA among mixed phenotypes while accounting for heterogeneities across EHRs. Simulations demonstrate that mixWAS outperforms the widely used MPA detection method, Phenome-wide association study (PheWAS), across diverse scenarios. Applying mixWAS to data from seven EHRs in the US, we identified 4,534 MPA among blood lipids, BMI, and circulatory diseases. Validation in an independent EHR data from UK confirmed 97.7% of the associations. mixWAS fundamentally improves the detection of MPA and is available as a free, open-source software.

Subset scanning for multi-trait analysis using GWAS summary statistics

MOTIVATION

Multi-trait analysis has been shown to have greater statistical power than single-trait analysis. Most of the existing multi-trait analysis methods only work with a limited number of traits and usually prioritize high statistical power over identifying relevant traits, which heavily rely on domain knowledge.

RESULTS

To handle diseases and traits with obscure etiology, we developed TraitScan, a powerful and fast algorithm that identifies potential pleiotropic traits from a moderate or large number of traits (e.g. dozens to thousands) and tests the association between one genetic variant and the selected traits. TraitScan can handle either individual-level or summary-level GWAS data. We evaluated TraitScan using extensive simulations and found that it outperformed existing methods in terms of both testing power and trait selection when sparsity was low or modest. We then applied it to search for traits associated with Ewing Sarcoma, a rare bone tumor with peak onset in adolescence, among 754 traits in UK Biobank. Our analysis revealed a few promising traits worthy of further investigation, highlighting the use of TraitScan for more effective multi-trait analysis as biobanks emerge. We also extended TraitScan to search and test association with a polygenic risk score and genetically imputed gene expression.

AVAILABILITY AND IMPLEMENTATION

Our algorithm is implemented in an R package "TraitScan" available at https://github.com/RuiCao34/TraitScan.

Isoform-level transcriptome-wide association uncovers genetic risk mechanisms for neuropsychiatric disorders in the human brain

Methods integrating genetics with transcriptomic reference panels prioritize risk genes and mechanisms at only a fraction of trait-associated genetic loci, due in part to an overreliance on total gene expression as a molecular outcome measure. This challenge is particularly relevant for the brain, in which extensive splicing generates multiple distinct transcript-isoforms per gene. Due to complex correlation structures, isoform-level modeling from cis-window variants requires methodological innovation. Here we introduce isoTWAS, a multivariate, stepwise framework integrating genetics, isoform-level expression and phenotypic associations. Compared to gene-level methods, isoTWAS improves both isoform and gene expression prediction, yielding more testable genes, and increased power for discovery of trait associations within genome-wide association study loci across 15 neuropsychiatric traits. We illustrate multiple isoTWAS associations undetectable at the gene-level, prioritizing isoforms of AKT3, CUL3 and HSPD1 in schizophrenia and PCLO with multiple disorders. Results highlight the importance of incorporating isoform-level resolution within integrative approaches to increase discovery of trait associations, especially for brain-relevant traits.

A cross-disorder study to identify causal relationships, shared genetic variants, and genes across 21 digestive disorders

Digestive disorders are a significant contributor to the global burden of disease and seriously affect human quality of life. Research has already confirmed the presence of pleiotropic genetic loci among digestive disorders, and studies have explored shared genetic factors among pan-cancers, including various malignant digestive disorders. However, most cross-phenotype studies within the digestive tract system have been limited to a few traits, with no systematic coverage of common benign and malignant digestive disorders. Here, we analyzed data from the UK Biobank to investigate 21 digestive disorders, exploring the genetic correlations and causal relationships between diseases, as well as the common genetic factors and potential biological pathways driving these relationships. Our findings confirmed the extensive genetic correlation and causal relationship between digestive disorders, providing important insights into the genetic etiology, causality, disease prevention, and clinical treatment of diseases.

Genome-Wide Assessment of Shared Genetic Architecture Between Rheumatoid Arthritis and Cardiovascular Diseases

Background Patients with rheumatoid arthritis (RA) have a 2- to 10-fold increased risk of cardiovascular disease (CVD), but the biological mechanisms and existence of causality underlying such associations remain to be investigated. We aimed to investigate the genetic associations and underlying mechanisms between RA and CVD by leveraging large-scale genomic data and genetic cross-trait analytic approaches. Methods and Results Within UK Biobank data, we examined the genetic correlation, shared genetics, and potential causality between RA (Ncases=6754, Ncontrols=452 384) and cardiovascular diseases (CVD, Ncases=44 238, Ncontrols=414 900) using linkage disequilibrium score regression, cross-trait meta-analysis, and Mendelian randomization. We observed significant genetic correlations of RA with myocardial infarction (rg:0.40 [95% CI, 0.24-0.56), angina (rg:0.42 [95% CI, 0.28-0.56]), coronary heart diseases (rg:0.41 [95% CI, 0.27-0.55]), and CVD (rg:0.43 [95% CI, 0.29-0.57]) after correcting for multiple testing (P<0.05/5). When stratified by frequent use of analgesics, we found increased genetic correlation between RA and CVD among participants without aspirin usage (rg:0.54 [95% CI, 0.30-0.78] for angina; Pvalue=6.69×10-6) and among participants with paracetamol usage (rg:0.75 [95% CI, 0.20-1.29] for myocardial infarction; Pvalue=8.90×10-3), whereas others remained similar. Cross-trait meta-analysis identified 9 independent shared loci between RA and CVD, including PTPN22 at chr1p13.2, BCL2L11 at chr2q13, and CCR3 at chr3p21.31 (Psingle trait<1×10-3 and Pmeta<5×10-8), highlighting potential shared pathogenesis including accelerating atherosclerosis, upregulating oxidative stress, and vascular permeability. Finally, Mendelian randomization estimates showed limited evidence of causality between RA and CVD. Conclusions Our results supported shared genetic pathogenesis rather than causality in explaining the observed association between RA and CVD. The identified shared genetic factors provided insights into potential novel therapeutic target for RA-CVD comorbidities.

Rare variants with large effects provide functional insights into the pathology of migraine subtypes, with and without aura

Migraine is a complex neurovascular disease with a range of severity and symptoms, yet mostly studied as one phenotype in genome-wide association studies (GWAS). Here we combine large GWAS datasets from six European populations to study the main migraine subtypes, migraine with aura (MA) and migraine without aura (MO). We identified four new MA-associated variants (in PRRT2, PALMD, ABO and LRRK2) and classified 13 MO-associated variants. Rare variants with large effects highlight three genes. A rare frameshift variant in brain-expressed PRRT2 confers large risk of MA and epilepsy, but not MO. A burden test of rare loss-of-function variants in SCN11A, encoding a neuron-expressed sodium channel with a key role in pain sensation, shows strong protection against migraine. Finally, a rare variant with cis-regulatory effects on KCNK5 confers large protection against migraine and brain aneurysms. Our findings offer new insights with therapeutic potential into the complex biology of migraine and its subtypes.

Trait selection strategy in multi-trait GWAS: Boosting SNPs discoverability

Since the first Genome-Wide Association Studies (GWAS), thousands of variant-trait associations have been discovered. However, the sample size required to detect additional variants using standard univariate association screening is increasingly prohibitive. Multi-trait GWAS offers a relevant alternative: it can improve statistical power and lead to new insights about gene function and the joint genetic architecture of human phenotypes. Although many methodological hurdles of multi-trait testing have been discussed, the strategy to select trait, among overwhelming possibilities, has been overlooked. In this study, we conducted extensive multi-trait tests using JASS (Joint Analysis of Summary Statistics) and assessed which genetic features of the analysed sets were associated with an increased detection of variants as compared to univariate screening. Our analyses identified multiple factors associated with the gain in the association detection in multi-trait tests. Together, these factors of the analysed sets are predictive of the gain of the multi-trait test (Pearson's ρ equal to 0.43 between the observed and predicted gain, P < 1.6 × 10-60). Applying an alternative multi-trait approach (MTAG, multi-trait analysis of GWAS), we found that in most scenarios but particularly those with larger numbers of traits, JASS outperformed MTAG. Finally, we benchmark several strategies to select set of traits including the prevalent strategy of selecting clinically similar traits, which systematically underperformed selecting clinically heterogenous traits or selecting sets that issued from our data-driven models. This work provides a unique picture of the determinant of multi-trait GWAS statistical power and outline practical strategies for multi-trait testing.

GWAS meta-analysis of over 29,000 people with epilepsy identifies 26 risk loci and subtype-specific genetic architecture

Epilepsy is a highly heritable disorder affecting over 50 million people worldwide, of which about one-third are resistant to current treatments. Here we report a multi-ancestry genome-wide association study including 29,944 cases, stratified into three broad categories and seven subtypes of epilepsy, and 52,538 controls. We identify 26 genome-wide significant loci, 19 of which are specific to genetic generalized epilepsy (GGE). We implicate 29 likely causal genes underlying these 26 loci. SNP-based heritability analyses show that common variants explain between 39.6% and 90% of genetic risk for GGE and its subtypes. Subtype analysis revealed markedly different genetic architectures between focal and generalized epilepsies. Gene-set analyses of GGE signals implicate synaptic processes in both excitatory and inhibitory neurons in the brain. Prioritized candidate genes overlap with monogenic epilepsy genes and with targets of current antiseizure medications. Finally, we leverage our results to identify alternate drugs with predicted efficacy if repurposed for epilepsy treatment.

GCPBayes pipeline: a tool for exploring pleiotropy at the gene level

Cross-phenotype association using gene-set analysis can help to detect pleiotropic genes and inform about common mechanisms between diseases. Although there are an increasing number of statistical methods for exploring pleiotropy, there is a lack of proper pipelines to apply gene-set analysis in this context and using genome-scale data in a reasonable running time. We designed a user-friendly pipeline to perform cross-phenotype gene-set analysis between two traits using GCPBayes, a method developed by our team. All analyses could be performed automatically by calling for different scripts in a simple way (using a Shiny app, Bash or R script). A Shiny application was also developed to create different plots to visualize outputs from GCPBayes. Finally, a comprehensive and step-by-step tutorial on how to use the pipeline is provided in our group's GitHub page. We illustrated the application on publicly available GWAS (genome-wide association studies) summary statistics data to identify breast cancer and ovarian cancer susceptibility genes. We have shown that the GCPBayes pipeline could extract pleiotropic genes previously mentioned in the literature, while it also provided new pleiotropic genes and regions that are worthwhile for further investigation. We have also provided some recommendations about parameter selection for decreasing computational time of GCPBayes on genome-scale data.

Interplay of chronic obstructive pulmonary disease and colorectal cancer development: unravelling the mediating role of fatty acids through a comprehensive multi-omics analysis

BACKGROUND

Chronic obstructive pulmonary disease (COPD) patients often exhibit gastrointestinal symptoms, A potential association between COPD and Colorectal Cancer (CRC) has been indicated, warranting further examination.

METHODS

In this study, we collected COPD and CRC data from the National Health and Nutrition Examination Survey, genome-wide association studies, and RNA sequence for a comprehensive analysis. We used weighted logistic regression to explore the association between COPD and CRC incidence risk. Mendelian randomization analysis was performed to assess the causal relationship between COPD and CRC, and cross-phenotype meta-analysis was conducted to pinpoint crucial loci. Multivariable mendelian randomization was used to uncover mediating factors connecting the two diseases. Our results were validated using both NHANES and GEO databases.

RESULTS

In our analysis of the NHANES dataset, we identified COPD as a significant contributing factor to CRC development. MR analysis revealed that COPD increased the risk of CRC onset and progression (OR: 1.16, 95% CI 1.01-1.36). Cross-phenotype meta-analysis identified four critical genes associated with both CRC and COPD. Multivariable Mendelian randomization suggested body fat percentage, omega-3, omega-6, and the omega-3 to omega-6 ratio as potential mediating factors for both diseases, a finding consistent with the NHANES dataset. Further, the interrelation between fatty acid-related modules in COPD and CRC was demonstrated via weighted gene co-expression network analysis and Kyoto Encyclopedia of Genes and Genomes enrichment results using RNA expression data.

CONCLUSIONS

This study provides novel insights into the interplay between COPD and CRC, highlighting the potential impact of COPD on the development of CRC. The identification of shared genes and mediating factors related to fatty acid metabolism deepens our understanding of the underlying mechanisms connecting these two diseases.

Topic modeling identifies novel genetic loci associated with multimorbidities in UK Biobank

Many diseases show patterns of co-occurrence, possibly driven by systemic dysregulation of underlying processes affecting multiple traits. We have developed a method (treeLFA) for identifying such multimorbidities from routine health-care data, which combines topic modeling with an informative prior derived from medical ontology. We apply treeLFA to UK Biobank data and identify a variety of topics representing multimorbidity clusters, including a healthy topic. We find that loci identified using topic weights as traits in a genome-wide association study (GWAS) analysis, which we validated with a range of approaches, only partially overlap with loci from GWASs on constituent single diseases. We also show that treeLFA improves upon existing methods like latent Dirichlet allocation in various ways. Overall, our findings indicate that topic models can characterize multimorbidity patterns and that genetic analysis of these patterns can provide insight into the etiology of complex traits that cannot be determined from the analysis of constituent traits alone.

Integrative analyses for the identification of idiopathic pulmonary fibrosis-associated genes and shared loci with other diseases

BACKGROUND

Although genome-wide association studies (GWAS) have identified many genomic regions associated with idiopathic pulmonary fibrosis (IPF), the causal genes and functions remain largely unknown. Many single-cell expression data have become available for IPF, and there is increasing evidence suggesting a shared genetic basis between IPF and other diseases.

METHODS

We conducted integrative analyses to improve the power of GWAS. First, we calculated global and local genetic correlations to identify IPF genetically associated traits and local regions. Then, we prioritised candidate genes contributing to local genetic correlation. Second, we performed transcriptome-wide association analysis (TWAS) of 44 tissues to identify candidate genes whose genetically predicted expression level is associated with IPF. To replicate our findings and investigate the regulatory role of the transcription factors (TF) in identified candidate genes, we first conducted the heritability enrichment analysis in TF binding sites. Then, we examined the enrichment of the TF target genes in cell-type-specific differentially expressed genes (DEGs) identified from single-cell expression data of IPF and healthy lung samples.

FINDINGS

We identified 12 candidate genes across 13 genomic regions using local genetic correlation, including the POT1 locus (p value=0.00041), which contained variants with protective effects on lung cancer but increasing IPF risk. We identified another 13 novel genes using TWAS. Two TFs, MAFK and SMAD2, showed significant enrichment in both partitioned heritability and cell-type-specific DEGs.

INTERPRETATION

Our integrative analysis identified new genes for IPF susceptibility and expanded the understanding of the complex genetic architecture and disease mechanism of IPF.

Rare variants confer shared susceptibility to gastrointestinal tract cancer risk

Background

Cancers arising within the gastrointestinal tract are complex disorders involving genetic events that cause the conversion of normal tissue to premalignant lesions and malignancy. Shared genetic features are reported in epithelial-based gastrointestinal cancers which indicate common susceptibility among this group of malignancies. In addition, the contribution of rare variants may constitute parts of genetic susceptibility.

Methods

A cross-cancer analysis of 38,171 shared rare genetic variants from genome-wide association assays was conducted, which included data from 3,194 cases and 1,455 controls across three cancer sites (esophageal, gastric and colorectal). The SNP-level association was performed by multivariate logistic regression analyses for single cancer, followed by association analysis for SubSETs (ASSET) to adjust the bias of overlapping controls. Gene-level analyses were conducted by SKAT-O, with multiple comparison adjustments by false discovery rate (FDR). Based on the significant genes indicated by SKATO analysis, pathways analysis was conducted using Gene Ontology (GO), the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome databases.

Results

Meta-analysis in three gastrointestinal (GI) cancers identified 13 novel susceptibility loci that reached genome-wide significance (P _ASSET< 5×10^-8). SKAT-O analysis revealed EXOC6, LRP5L and MIR1263/LINC01324 to be significant genes shared by GI cancers (P _adj<0.05, P _FDR<0.05). Furthermore, GO pathway analysis identified significant enrichment of synaptic transmission and neuron development pathways shared by all three cancer types.

Conclusion

Rare variants and the corresponding genes potentially contribute to shared susceptibility in different GI cancer types. The discovery of these novel variants and genes offers new insights for the carcinogenic mechanisms and missing heritability of GI cancers.

Shared genetics and causal relationships between major depressive disorder and COVID-19 related traits: a large-scale genome-wide cross-trait meta-analysis

Introduction

The comorbidity between major depressive disorder (MDD) and coronavirus disease of 2019 (COVID-19) related traits have long been identified in clinical settings, but their shared genetic foundation and causal relationships are unknown. Here, we investigated the genetic mechanisms behind COVID-19 related traits and MDD using the cross-trait meta-analysis, and evaluated the underlying causal relationships between MDD and 3 different COVID-19 outcomes (severe COVID-19, hospitalized COVID-19, and COVID-19 infection).

Methods

In this study, we conducted a comprehensive analysis using the most up-to-date and publicly available GWAS summary statistics to explore shared genetic etiology and the causality between MDD and COVID-19 outcomes. We first used genome-wide cross-trait meta-analysis to identify the pleiotropic genomic SNPs and the genes shared by MDD and COVID-19 outcomes, and then explore the potential bidirectional causal relationships between MDD and COVID-19 outcomes by implementing a bidirectional MR study design. We further conducted functional annotations analyses to obtain biological insight for shared genes from the results of cross-trait meta-analysis.

Results

We have identified 71 SNPs located on 25 different genes are shared between MDD and COVID-19 outcomes. We have also found that genetic liability to MDD is a causal factor for COVID-19 outcomes. In particular, we found that MDD has causal effect on severe COVID-19 (OR = 1.832, 95% CI = 1.037-3.236) and hospitalized COVID-19 (OR = 1.412, 95% CI = 1.021-1.953). Functional analysis suggested that the shared genes are enriched in Cushing syndrome, neuroactive ligand-receptor interaction.

Discussion

Our findings provide convincing evidence on shared genetic etiology and causal relationships between MDD and COVID-19 outcomes, which is crucial to prevention, and therapeutic treatment of MDD and COVID-19.

Pan-cancer and cross-population genome-wide association studies dissect shared genetic backgrounds underlying carcinogenesis

Integrating genomic data of multiple cancers allows de novo cancer grouping and elucidating the shared genetic basis across cancers. Here, we conduct the pan-cancer and cross-population genome-wide association study (GWAS) meta-analysis and replication studies on 13 cancers including 250,015 East Asians (Biobank Japan) and 377,441 Europeans (UK Biobank). We identify ten cancer risk variants including five pleiotropic associations (e.g., rs2076295 at DSP on 6p24 associated with lung cancer and rs2525548 at TRIM4 on 7q22 nominally associated with six cancers). Quantifying shared heritability among the cancers detects positive genetic correlations between breast and prostate cancer across populations. Common genetic components increase the statistical power, and the large-scale meta-analysis of 277,896 breast/prostate cancer cases and 901,858 controls identifies 91 newly genome-wide significant loci. Enrichment analysis of pathways and cell types reveals shared genetic backgrounds across said cancers. Focusing on genetically correlated cancers can contribute to enhancing our insights into carcinogenesis.

Genome-wide analyses characterize shared heritability among cancers and identify novel cancer susceptibility regions

BACKGROUND

The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci.

METHODS

We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci.

RESULTS

We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci.

CONCLUSIONS

Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types.

Identification of new risk loci shared across systemic vasculitides points towards potential target genes for drug repurposing

OBJECTIVES

The number of susceptibility loci currently associated with vasculitis is lower than in other immune-mediated diseases due in part to small cohort sizes, a consequence of the low prevalence of vasculitides. This study aimed to identify new genetic risk loci for the main systemic vasculitides through a comprehensive analysis of their genetic overlap.

METHODS

Genome-wide data from 8467 patients with any of the main forms of vasculitis and 29 795 healthy controls were meta-analysed using ASSET. Pleiotropic variants were functionally annotated and linked to their target genes. Prioritised genes were queried in DrugBank to identify potentially repositionable drugs for the treatment of vasculitis.

RESULTS

Sixteen variants were independently associated with two or more vasculitides, 15 of them representing new shared risk loci. Two of these pleiotropic signals, located close to CTLA4 and CPLX1, emerged as novel genetic risk loci in vasculitis. Most of these polymorphisms appeared to affect vasculitis by regulating gene expression. In this regard, for some of these common signals, potential causal genes were prioritised based on functional annotation, including CTLA4, RNF145, IL12B, IL5, IRF1, IFNGR1, PTK2B, TRIM35, EGR2 and ETS2, each of which has key roles in inflammation. In addition, drug repositioning analysis showed that several drugs, including abatacept and ustekinumab, could be potentially repurposed in the management of the analysed vasculitides.

CONCLUSIONS

We identified new shared risk loci with functional impact in vasculitis and pinpointed potential causal genes, some of which could represent promising targets for the treatment of vasculitis.

Genome-wide analysis of genetic pleiotropy and causal genes across three age-related ocular disorders

Age-related macular degeneration (AMD), cataract, and glaucoma are leading causes of blindness worldwide. Previous genome-wide association studies (GWASs) have revealed a variety of susceptible loci associated with age-related ocular disorders, yet the genetic pleiotropy and causal genes across these diseases remain poorly understood. By leveraging large-scale genetic and observational data from ocular disease GWASs and UK Biobank (UKBB), we found significant pairwise genetic correlations and consistent epidemiological associations among these ocular disorders. Cross-disease meta-analysis uncovered seven pleiotropic loci, three of which were replicated in an additional cohort. Integration of variants in pleiotropic loci and multiple single-cell omics data identified that Müller cells and astrocytes were likely trait-related cell types underlying ocular comorbidity. In addition, we comprehensively integrated eye-specific gene expression quantitative loci (eQTLs), epigenomic profiling, and 3D genome data to prioritize causal pleiotropic genes. We found that pleiotropic genes were essential in nerve development and eye pigmentation, and targetable by aflibercept and pilocarpine for the treatment of AMD and glaucoma. These findings will not only facilitate the mechanistic research of ocular comorbidities but also benefit the therapeutic optimization of age-related ocular diseases.

The genetic overlap between Alzheimer’s disease, amyotrophic lateral sclerosis, Lewy body dementia, and Parkinson’s disease

Neurodegenerative diseases are a group of disorders characterized by neuronal cell death causing a variety of physical and mental problems. While these disorders can be characterized by their phenotypic presentation within the nervous system, their aetiologies differ to varying degrees. The majority of previous genetic evidence for overlap between neurodegenerative diseases has been pairwise. In this study, we aimed to identify overlap between the 4 investigated neurodegenerative disorders (Alzheimer's disease, amyotrophic lateral sclerosis, Lewy body dementia, and Parkinson's disease) at the variant, gene, genomic locus, gene-set, cell, or tissue level, with specific interest in overlap between 3 or more diseases. Using local genetic correlation, we found 2 loci (TMEM175 and HLA) that were shared across 3 disorders. We also highlighted genes, genomic loci, gene sets, cell types, and tissue types which may be important to 2 or more disorders by analyzing the association of variants with a common factor estimated from the 4 disorders. Our study successfully highlighted genetic loci and tissues associated with 2 or more neurodegenerative diseases.

Shared genetics of psychiatric disorders and type 2 diabetes:a large-scale genome-wide cross-trait analysis

BACKGROUND

Individuals with psychiatric disorders have elevated rates of type 2 diabetes comorbidity. Although little is known about the shared genetics and causality of this association. Thus, we aimed to investigate shared genetics and causal link between different type 2 diabetes and psychiatric disorders.

METHODS

We conducted a large-scale genome-wide cross-trait association study(GWAS) to investigate genetic overlap between type 2 diabetes and anorexia nervosa, attention deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, major depressive disorder, obsessive-compulsive disorder, schizophrenia, anxiety disorders and Tourette syndrome. By post-GWAS functional analysis, we identify variants genes expression in various tissues. Enrichment pathways, potential protein interaction and mendelian randomization also provided to research the relationship between type 2 diabetes and psychiatric disorders.

RESULTS

We discovered that type 2 diabetes and psychiatric disorders had a significant correlation. We identified 138 related loci, 32 were novel loci. Post-GWAS analysis revealed that 86 differentially expressed genes were located in different brain regions and peripheral blood in type 2 diabetes and related psychiatric disorders. MAPK signaling pathway plays an important role in neural development and insulin signaling. In addition, there is a causal relationship between T2D and mental disorders. In PPI analysis, the central genes of the DEG PPI network were FTO and TCF7L2.

CONCLUSION

This large-scale genome-wide cross-trait analysis identified shared genetics andpotential causal links between type 2 diabetes and related psychiatric disorders, suggesting potential new biological functions in common among them.

Multivariate genome-wide association meta-analysis of over 1 million subjects identifies loci underlying multiple substance use disorders

Genetic liability to substance use disorders can be parsed into loci that confer general or substance-specific addiction risk. We report a multivariate genome-wide association meta-analysis that disaggregates general and substance-specific loci for published summary statistics of problematic alcohol use, problematic tobacco use, cannabis use disorder, and opioid use disorder in a sample of 1,025,550 individuals of European descent and 92,630 individuals of African descent. Nineteen independent SNPs were genome-wide significant (P < 5e-8) for the general addiction risk factor (addiction-rf), which showed high polygenicity. Across ancestries, PDE4B was significant (among other genes), suggesting dopamine regulation as a cross-substance vulnerability. An addiction-rf polygenic risk score was associated with substance use disorders, psychopathologies, somatic conditions, and environments associated with the onset of addictions. Substance-specific loci (9 for alcohol, 32 for tobacco, 5 for cannabis, 1 for opioids) included metabolic and receptor genes. These findings provide insight into genetic risk loci for substance use disorders that could be leveraged as treatment targets.

Investigation of convergent and divergent genetic influences underlying schizophrenia and alcohol use disorder

BACKGROUND

Alcohol use disorder (AUD) and schizophrenia (SCZ) frequently co-occur, and large-scale genome-wide association studies (GWAS) have identified significant genetic correlations between these disorders.

METHODS

We used the largest published GWAS for AUD (total cases = 77 822) and SCZ (total cases = 46 827) to identify genetic variants that influence both disorders (with either the same or opposite direction of effect) and those that are disorder specific.

RESULTS

We identified 55 independent genome-wide significant single nucleotide polymorphisms with the same direction of effect on AUD and SCZ, 8 with robust effects in opposite directions, and 98 with disorder-specific effects. We also found evidence for 12 genes whose pleiotropic associations with AUD and SCZ are consistent with mediation via gene expression in the prefrontal cortex. The genetic covariance between AUD and SCZ was concentrated in genomic regions functional in brain tissues (p = 0.001).

CONCLUSIONS

Our findings provide further evidence that SCZ shares meaningful genetic overlap with AUD.

Comparison of Low-Value Care Among Commercial and Medicaid Enrollees

BACKGROUND

Low-value healthcare is costly and inefficient and may adversely affect patient outcomes. Despite increases in low-value service use, little is known about how the receipt of low-value care differs across payers.

OBJECTIVE

To evaluate differences in the use of low-value care between patients with commercial versus Medicaid coverage.

DESIGN

Retrospective observational analysis of the 2017 Rhode Island All-payer Claims Database, estimating the probability of receiving each of 14 low-value services between commercial and Medicaid enrollees, adjusting for patient sociodemographic and clinical characteristics. Ensemble machine learning minimized the possibility of model misspecification.

PARTICIPANTS

Medicaid and commercial enrollees aged 18-64 with continuous coverage and an encounter at which they were at risk of receiving a low-value service.

INTERVENTION

Enrollment in Medicaid or Commercial insurance.

MAIN MEASURES

Use of one of 14 validated measures of low-value care.

KEY RESULTS

Among 110,609 patients, Medicaid enrollees were younger, had more comorbidities, and were more likely to be female than commercial enrollees. Medicaid enrollees had higher rates of use for 7 low-value care measures, and those with commercial coverage had higher rates for 5 measures. Across all measures of low-value care, commercial enrollees received more (risk difference [RD] 6.8 percentage points; CI: 6.6 to 7.0) low-value services than their counterparts with Medicaid. Commercial enrollees were also more likely to receive low-value services typically performed in the emergency room (RD 11.4 percentage points; CI: 10.7 to 12.2) and services that were less expensive (RD 15.3 percentage points; CI 14.6 to 16.0).

CONCLUSION

Differences in the provision of low-value care varied across measures, though average use was slightly higher among commercial than Medicaid enrollees. This difference was more pronounced for less expensive services indicating that financial incentives may not be the sole driver of low-value care.

A genome-wide cross-cancer meta-analysis highlights the shared genetic links of five solid cancers

Breast, ovarian, prostate, lung, and head/neck cancers are five solid cancers with complex interrelationships. However, the shared genetic factors of the five cancers were often revealed either by the combination of individual genome-wide association study (GWAS) approach or by the fixed-effect model-based meta-analysis approach with practically impossible assumptions. Here, we presented a random-effect model-based cross-cancer meta-analysis framework for identifying the genetic variants jointly influencing the five solid cancers. A comprehensive genetic correlation analysis (genome-wide, partitioned, and local) approach was performed by using GWAS summary statistics of the five cancers, and we observed three cancer pairs with significant genetic correlation: breast-ovarian cancer (r _g = 0.221, p = 0.0003), breast-lung cancer (r _g = 0.234, p = 7.6 × 10^-6), and lung-head/neck cancer (r _g = 0.652, p = 0.010). Furthermore, a random-effect model-based cross-trait meta-analysis was conducted for each significant cancer pair, and we found 27 shared genetic loci between breast and ovarian cancers, 18 loci between breast and lung cancers, and three loci between lung and head/neck cancers. Functional analysis indicates that the shared genes are enriched in human T-cell leukemia virus 1 infection (HTLV-1) and antigen processing and presentation (APP) pathways. Our study investigates the shared genetic links across five solid cancers and will help to reveal their potential molecular mechanisms.

Inflammatory and infectious upper respiratory diseases associate with 41 genomic loci and type 2 inflammation

Inflammatory and infectious upper respiratory diseases (ICD-10: J30-J39), such as diseases of the sinonasal tract, pharynx and larynx, are growing health problems yet their genomic similarity is not known. We analyze genome-wide association to eight upper respiratory diseases (61,195 cases) among 260,405 FinnGen participants, meta-analyzing diseases in four groups based on an underlying genetic correlation structure. Aiming to understand which genetic loci contribute to susceptibility to upper respiratory diseases in general and its subtypes, we detect 41 independent genome-wide significant loci, distinguishing impact on sinonasal or pharyngeal diseases, or both. Fine-mapping implicated non-synonymous variants in nine genes, including three linked to immune-related diseases. Phenome-wide analysis implicated asthma and atopic dermatitis at sinonasal disease loci, and inflammatory bowel diseases and other immune-mediated disorders at pharyngeal disease loci. Upper respiratory diseases also genetically correlated with autoimmune diseases such as rheumatoid arthritis, autoimmune hypothyroidism, and psoriasis. Finally, we associated separate gene pathways in sinonasal and pharyngeal diseases that both contribute to type 2 immunological reaction. We show shared heritability among upper respiratory diseases that extends to several immune-mediated diseases with diverse mechanisms, such as type 2 high inflammation.

DNA methylation signatures as biomarkers of socioeconomic position

This review article provides a framework for the use of deoxyribonucleic acid (DNA) methylation (DNAm) biomarkers to study the biological embedding of socioeconomic position (SEP) and summarizes the latest developments in the area. It presents the emerging literature showing associations between individual- and neighborhood-level SEP exposures and DNAm across the life course. In contrast to questionnaire-based methods of assessing SEP, we suggest that DNAm biomarkers may offer an accessible metric to study questions about SEP and health outcomes, acting as a personal dosimeter of exposure. However, further work remains in standardizing SEP measures across studies and evaluating consistency across domains, tissue types, and time periods. Meta-analyses of epigenetic associations with SEP are offered as one approach to confirm the replication of DNAm loci across studies. The development of DNAm biomarkers of SEP would provide a method for examining its impact on health outcomes in a more robust way, increasing the rigor of epidemiological studies.

Distinct germline genetic susceptibility profiles identified for common non-Hodgkin lymphoma subtypes

Lymphoma risk is elevated for relatives with common non-Hodgkin lymphoma (NHL) subtypes, suggesting shared genetic susceptibility across subtypes. To evaluate the extent of mutual heritability among NHL subtypes and discover novel loci shared among subtypes, we analyzed data from eight genome-wide association studies within the InterLymph Consortium, including 10,629 cases and 9505 controls. We utilized Association analysis based on SubSETs (ASSET) to discover loci for subsets of NHL subtypes and evaluated shared heritability across the genome using Genome-wide Complex Trait Analysis (GCTA) and polygenic risk scores. We discovered 17 genome-wide significant loci (P < 5 × 10-8) for subsets of NHL subtypes, including a novel locus at 10q23.33 (HHEX) (P = 3.27 × 10-9). Most subset associations were driven primarily by only one subtype. Genome-wide genetic correlations between pairs of subtypes varied broadly from 0.20 to 0.86, suggesting substantial heterogeneity in the extent of shared heritability among subtypes. Polygenic risk score analyses of established loci for different lymphoid malignancies identified strong associations with some NHL subtypes (P < 5 × 10-8), but weak or null associations with others. Although our analyses suggest partially shared heritability and biological pathways, they reveal substantial heterogeneity among NHL subtypes with each having its own distinct germline genetic architecture.

Genetic heterogeneity: Challenges, impacts, and methods through an associative lens

Genetic heterogeneity describes the occurrence of the same or similar phenotypes through different genetic mechanisms in different individuals. Robustly characterizing and accounting for genetic heterogeneity is crucial to pursuing the goals of precision medicine, for discovering novel disease biomarkers, and for identifying targets for treatments. Failure to account for genetic heterogeneity may lead to missed associations and incorrect inferences. Thus, it is critical to review the impact of genetic heterogeneity on the design and analysis of population level genetic studies, aspects that are often overlooked in the literature. In this review, we first contextualize our approach to genetic heterogeneity by proposing a high-level categorization of heterogeneity into "feature," "outcome," and "associative" heterogeneity, drawing on perspectives from epidemiology and machine learning to illustrate distinctions between them. We highlight the unique nature of genetic heterogeneity as a heterogeneous pattern of association that warrants specific methodological considerations. We then focus on the challenges that preclude effective detection and characterization of genetic heterogeneity across a variety of epidemiological contexts. Finally, we discuss systems heterogeneity as an integrated approach to using genetic and other high-dimensional multi-omic data in complex disease research.

Evidence that the pituitary gland connects type 2 diabetes mellitus and schizophrenia based on large-scale trans-ethnic genetic analyses

BACKGROUND

Previous studies on European (EUR) samples have obtained inconsistent results regarding the genetic correlation between type 2 diabetes mellitus (T2DM) and Schizophrenia (SCZ). A large-scale trans-ethnic genetic analysis may provide additional evidence with enhanced power.

OBJECTIVE

We aimed to explore the genetic basis for both T2DM and SCZ based on large-scale genetic analyses of genome-wide association study (GWAS) data from both East Asian (EAS) and EUR subjects.

METHODS

A range of complementary approaches were employed to cross-validate the genetic correlation between T2DM and SCZ at the whole genome, autosomes (linkage disequilibrium score regression, LDSC), loci (Heritability Estimation from Summary Statistics, HESS), and causal variants (MiXeR and Mendelian randomization, MR) levels. Then, genome-wide and transcriptome-wide cross-trait/ethnic meta-analyses were performed separately to explore the effective shared organs, cells and molecular pathways.

RESULTS

A weak genome-wide negative genetic correlation between SCZ and T2DM was found for the EUR (r_g = - 0.098, P = 0.009) and EAS (r_g =- 0.053 and P = 0.032) populations, which showed no significant difference between the EUR and EAS populations (P = 0.22). After Bonferroni correction, the r_g remained significant only in the EUR population. Similar results were obtained from analyses at the levels of autosomes, loci and causal variants. 25 independent variants were firstly identified as being responsible for both SCZ and T2DM. The variants associated with the two disorders were significantly correlated to the gene expression profiles in the brain (P = 1.1E-9) and pituitary gland (P = 1.9E-6). Then, 61 protein-coding and non-coding genes were identified as effective genes in the pituitary gland (P < 9.23E-6) and were enriched in metabolic pathways related to glutathione mediated arsenate detoxification and to D-myo-inositol-trisphosphate.

CONCLUSION

Here, we show that a negative genetic correlation exists between SCZ and T2DM at the whole genome, autosome, locus and causal variant levels. We identify pituitary gland as a common effective organ for both diseases, in which non-protein-coding effective genes, such as lncRNAs, may be responsible for the negative genetic correlation. This highlights the importance of molecular metabolism and neuroendocrine modulation in the pituitary gland, which may be responsible for the initiation of T2DM in SCZ patients.

Leveraging pleiotropy for joint analysis of genome-wide association studies with per trait interpretations

We introduce pleiotropic association test (PAT) for joint analysis of multiple traits using genome-wide association study (GWAS) summary statistics. The method utilizes the decomposition of phenotypic covariation into genetic and environmental components to create a likelihood ratio test statistic for each genetic variant. Though PAT does not directly interpret which trait(s) drive the association, a per trait interpretation of the omnibus p-value is provided through an extension to the meta-analysis framework, m-values. In simulations, we show PAT controls the false positive rate, increases statistical power, and is robust to model misspecifications of genetic effect. Additionally, simulations comparing PAT to three multi-trait methods, HIPO, MTAG, and ASSET, show PAT identified 15.3% more omnibus associations over the next best method. When these associations were interpreted on a per trait level using m-values, PAT had 37.5% more true per trait interpretations with a 0.92% false positive assignment rate. When analyzing four traits from the UK Biobank, PAT discovered 22,095 novel variants. Through the m-values interpretation framework, the number of per trait associations for two traits were almost tripled and were nearly doubled for another trait relative to the original single trait GWAS.

Shared genetics between classes of obesity and psychiatric disorders: A large-scale genome-wide cross-trait analysis

AIMS

Epidemiological studies demonstrate an association between classes of obesity and psychiatric disorders, although little is known about shared genetics and causality of association. Thus, we aimed to investigate shared genetics and causal link between different classes of obesity and psychiatric disorders.

METHODS

We used genome-wide association study (GWAS) summary data range from 9725 to 500,199 sample sizes of European descent, conducted a large-scale genome-wide cross-trait association study to investigate genetic overlap between the classes of obesity and anorexia nervosa, attention-deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, major depressive disorder, obsessive-compulsive disorder, schizophrenia, anxiety disorders and Tourette syndrome. We conducted transcriptome-wide association study analysis (TWAS) to identified variants regulated gene expression in those related disorders. Finally, pathway enrichment analysis to identified major pathways.

RESULTS

In the combined analysis, we replicated 211 previously reported loci and discovered 58 novel independent loci that were associated with all three classes of obesity and related psychiatric disorders. Functional analysis revealed that the identified variants regulated gene expression in major tissues belonging to exocrine/endocrine, digestive, circulatory, adipose, digestive, respiratory, and nervous systems, such as DCC, NEGR1, INO80E. Mendelian randomization analyses suggested that there may be a two-way or one-way causal relationship between obesity and psychiatric disorders.

CONCLUSION

This large-scale genome-wide cross-trait analysis identified shared genetics and potential causal links between classes of obesity and psychiatric disorders (attention deficit hyperactivity disorder, autism spectrum disorder, anorexia nervosa, major depressive disorder, schizophrenia, and obsessive-compulsive disorder). Such shared genetics suggests potential new biological functions in common among them.

Understanding the comorbidity between posttraumatic stress severity and coronary artery disease using genome-wide information and electronic health records

The association between coronary artery disease (CAD) and posttraumatic stress disorder (PTSD) contributes to the high morbidity and mortality observed for these conditions. To understand the dynamics underlying PTSD-CAD comorbidity, we investigated large-scale genome-wide association (GWA) statistics from the Million Veteran Program (MVP), the UK Biobank (UKB), the Psychiatric Genomics Consortium, and the CARDIoGRAMplusC4D Consortium. We observed a genetic correlation of CAD with PTSD case-control and quantitative outcomes, ranging from 0.18 to 0.32. To investigate possible cause-effect relationships underlying these genetic correlations, we performed a two-sample Mendelian randomization (MR) analysis, observing a significant bidirectional relationship between CAD and PTSD symptom severity. Genetically-determined PCL-17 (PTSD 17-item Checklist) total score was associated with increased CAD risk (odds ratio = 1.04; 95% confidence interval, 95% CI = 1.01-1.06). Conversely, CAD genetic liability was associated with reduced PCL-17 total score (beta = -0.42; 95% CI = -0.04 to -0.81). Because of these opposite-direction associations, we conducted a pleiotropic meta-analysis to investigate loci with concordant vs. discordant effects on PCL-17 and CAD, observing that concordant-effect loci were enriched for molecular pathways related to platelet amyloid precursor protein (beta = 1.53, p = 2.97 × 10-7) and astrocyte activation regulation (beta = 1.51, p = 2.48 × 10-6) while discordant-effect loci were enriched for biological processes related to lipid metabolism (e.g., triglyceride-rich lipoprotein particle clearance, beta = 2.32, p = 1.61 × 10-10). To follow up these results, we leveraged MVP and UKB electronic health records (EHR) to assess longitudinal changes in the association between CAD and posttraumatic stress severity. This EHR-based analysis highlighted that earlier CAD diagnosis is associated with increased PCL-total score later in life, while lower PCL total score was associated with increased risk of a later CAD diagnosis (Mann-Kendall trend test: MVP tau = 0.932, p < 2 × 10-16; UKB tau = 0.376, p = 0.005). In conclusion, both our genetically-informed analyses and our EHR-based follow-up investigation highlighted a bidirectional relationship between PTSD and CAD where multiple pleiotropic mechanisms are likely to be involved.

Substance use and common contributors to morbidity: A genetics perspective

Excessive substance use and substance use disorders (SUDs) are common, serious and relapsing medical conditions. They frequently co-occur with other diseases that are leading contributors to disability worldwide. While heavy substance use may potentiate the course of some of these illnesses, there is accumulating evidence suggesting common genetic architectures. In this narrative review, we focus on four heritable medical conditions - cardiometabolic disease, chronic pain, depression and COVID-19, which are commonly overlapping with, but not necessarily a direct consequence of, SUDs. We find persuasive evidence of underlying genetic liability that predisposes to both SUDs and chronic pain, depression, and COVID-19. For cardiometabolic disease, there is greater support for a potential causal influence of problematic substance use. Our review encourages de-stigmatization of SUDs and the assessment of substance use in clinical settings. We assert that identifying shared pathways of risk has high translational potential, allowing tailoring of treatments for multiple medical conditions.

Funding

SSR acknowledges T29KT0526, T32IR5226 and DP1DA054394; RLK acknowledges AA028292; AA acknowledges DA054869 & K02DA032573. The funders had no role in the conceptualization or writing of the paper.

Genome-wide association meta-analysis of 88,250 individuals highlights pleiotropic mechanisms of five ocular diseases in UK Biobank

BACKGROUND

Ocular diseases may exhibit common clinical symptoms and epidemiological comorbidity. However, the extent of pleiotropic mechanisms across ocular diseases remains unclear. We aim to examine shared genetic etiology in age-related macular degeneration (AMD), diabetic retinopathy (DR), glaucoma, retinal detachment (RD), and myopia.

METHODS

We analyzed genome-wide association analyses for the five ocular diseases in 43,877 cases and 44,373 controls of European ancestry from UK Biobank, estimated their genetic relationships (LDSC, GNOVA, and Genomic SEM), and identified pleiotropic loci (ASSET and METASOFT).

FINDINGS

The genetic correlation of common SNPs revealed a meaningful genetic structure within these diseases, identifying genetic correlations between AMD, DR, and glaucoma. Cross-trait meta-analysis identified 23 pleiotropic loci associated with at least two ocular diseases and 14 loci unique to individual disorders (non-pleiotropic). We found that the genes associated with these shared genetic loci are involved in neuron differentiation (P = 8.80 × 10^-6) and eye development systems (P = 3.86 × 10^-5), and single cell RNA sequencing data reveals their heightened gene expression from multipotent progenitors to other differentiated retinal cells during retina developmental process.

INTERPRETATION

These results highlighted the potential common genetic architectures among these ocular diseases and can deepen the understanding of the molecular mechanisms underlying the related diseases.

FUNDING

The National Natural Science Foundation of China (61871294), Zhejiang Provincial Natural Science Foundation of China (LR19C060001), and the Scientific Research Foundation for Talents of Wenzhou Medical University (QTJ18023).

Cross-ancestry genome-wide meta-analysis of 61,047 cases and 947,237 controls identifies new susceptibility loci contributing to lung cancer

To identify new susceptibility loci to lung cancer among diverse populations, we performed cross-ancestry genome-wide association studies in European, East Asian and African populations and discovered five loci that have not been previously reported. We replicated 26 signals and identified 10 new lead associations from previously reported loci. Rare-variant associations tended to be specific to populations, but even common-variant associations influencing smoking behavior, such as those with CHRNA5 and CYP2A6, showed population specificity. Fine-mapping and expression quantitative trait locus colocalization nominated several candidate variants and susceptibility genes such as IRF4 and FUBP1. DNA damage assays of prioritized genes in lung fibroblasts indicated that a subset of these genes, including the pleiotropic gene IRF4, potentially exert effects by promoting endogenous DNA damage.

An Atlas of Genetic Correlations and Genetically Informed Associations Linking Psychiatric and Immune-Related Phenotypes

Importance

Certain psychiatric and immune-related disorders are reciprocal risk factors. However, the nature of these associations is unclear.

Objective

To characterize the pleiotropy between psychiatric and immune-related traits, as well as risk factors of hypothesized relevance.

Design, Setting, and Participants

This genetic association study was conducted from July 10, 2020, to January 15, 2022. Analyses used genome-wide association (GWA) statistics related to 14 psychiatric traits; 13 immune-related phenotypes, ie, allergic, autoimmune, and inflammatory disorders; and 15 risk factors related to health-related behaviors, social determinants of health, and stress response. Genetically correlated psychiatric-immune pairs were assessed using 2-sample mendelian randomization (MR) with sensitivity analyses and multivariable adjustment for genetic associations of third variables. False discovery rate correction (Q value < .05) was applied for each analysis.

Exposures

Genetic associations.

Main Outcomes and Measures

Genetic correlations and MR association estimates with SEs and P values. A data-driven approach was used that did not test a priori planned hypotheses.

Results

A total of 44 genetically correlated psychiatric-immune pairs were identified, including 31 positive correlations (most consistently involving asthma, Crohn disease, hypothyroidism, and ulcerative colitis) and 13 negative correlations (most consistently involving allergic rhinitis and type 1 diabetes). Correlations with third variables were especially strong for psychiatric phenotypes. MR identified 7 associations of psychiatric phenotypes on immune-related phenotypes that were robust to multivariable adjustment, including the positive association of (1) the psychiatric cross-disorder phenotype with asthma (odds ratio [OR], 1.04; 95% CI, 1.02-1.06), Crohn disease (OR, 1.09; 95% CI, 1.05-1.14), and ulcerative colitis (OR, 1.09; 95% CI, 1.05-1.14); (2) major depression with asthma (OR, 1.25; 95% CI, 1.13-1.37); (3) schizophrenia with Crohn disease (OR, 1.12; 95% CI, 1.05-1.18) and ulcerative colitis (OR, 1.14; 95% CI, 1.07-1.21); and a negative association of risk tolerance with allergic rhinitis (OR, 0.77; 95% CI, 0.67-0.92).

Conclusions and Relevance

Results of this genetic association study suggest that genetic liability for psychiatric disorders was associated with liability for several immune disorders, suggesting that vertical pleiotropy related to behavioral traits (or correlated third variables) contributes to clinical associations observed in population-scale data.

Pinpointing novel risk loci for Lewy body dementia and the shared genetic etiology with Alzheimer's disease and Parkinson's disease: a large-scale multi-trait association analysis

BACKGROUND

The current genome-wide association study (GWAS) of Lewy body dementia (LBD) suffers from low power due to a limited sample size. In addition, the genetic determinants underlying LBD and the shared genetic etiology with Alzheimer's disease (AD) and Parkinson's disease (PD) remain poorly understood.

METHODS

Using the largest GWAS summary statistics of LBD to date (2591 cases and 4027 controls), late-onset AD (86,531 cases and 676,386 controls), and PD (33,674 cases and 449,056 controls), we comprehensively investigated the genetic basis of LBD and shared genetic etiology among LBD, AD, and PD. We first conducted genetic correlation analysis using linkage disequilibrium score regression (LDSC), followed by multi-trait analysis of GWAS (MTAG) and association analysis based on SubSETs (ASSET) to identify the trait-specific SNPs. We then performed SNP-level functional annotation to identify significant genomic risk loci paired with Bayesian fine-mapping and colocalization analysis to identify potential causal variants. Parallel gene-level analysis including GCTA-fastBAT and transcriptome-wide association analysis (TWAS) was implemented to explore novel LBD-associated genes, followed by pathway enrichment analysis to understand underlying biological mechanisms.

RESULTS

Pairwise LDSC analysis found positive genome-wide genetic correlations between LBD and AD (rg = 0.6603, se = 0.2001; P = 0.0010), between LBD and PD (rg = 0.6352, se = 0.1880; P = 0.0007), and between AD and PD (rg = 0.2136, se = 0.0860; P = 0.0130). We identified 13 significant loci for LBD, including 5 previously reported loci (1q22, 2q14.3, 4p16.3, 4q22.1, and 19q13.32) and 8 novel biologically plausible genetic associations (5q12.1, 5q33.3, 6p21.1, 8p23.1, 8p21.1, 16p11.2, 17p12, and 17q21.31), among which APOC1 (19q13.32), SNCA (4q22.1), TMEM175 (4p16.3), CLU (8p21.1), MAPT (17q21.31), and FBXL19 (16p11.2) were also validated by gene-level analysis. Pathway enrichment analysis of 40 common genes identified by GCTA-fastBAT and TWAS implicated significant role of neurofibrillary tangle assembly (GO:1902988, adjusted P = 1.55 × 10^-2).

CONCLUSIONS

Our findings provide novel insights into the genetic determinants of LBD and the shared genetic etiology and biological mechanisms of LBD, AD, and PD, which could benefit the understanding of the co-pathology as well as the potential treatment of these diseases simultaneously.

Large-scale genomic analyses reveal insights into pleiotropy across circulatory system diseases and nervous system disorders

Clinical and epidemiological studies have shown that circulatory system diseases and nervous system disorders often co-occur in patients. However, genetic susceptibility factors shared between these disease categories remain largely unknown. Here, we characterized pleiotropy across 107 circulatory system and 40 nervous system traits using an ensemble of methods in the eMERGE Network and UK Biobank. Using a formal test of pleiotropy, five genomic loci demonstrated statistically significant evidence of pleiotropy. We observed region-specific patterns of direction of genetic effects for the two disease categories, suggesting potential antagonistic and synergistic pleiotropy. Our findings provide insights into the relationship between circulatory system diseases and nervous system disorders which can provide context for future prevention and treatment strategies.

Shared components of heritability across genetically correlated traits

Most disease-associated genetic variants are pleiotropic, affecting multiple genetically correlated traits. Their pleiotropic associations can be mechanistically informative: if many variants have similar patterns of association, they may act via similar pleiotropic mechanisms, forming a shared component of heritability. We developed pleiotropic decomposition regression (PDR) to identify shared components and their underlying genetic variants. We validated PDR on simulated data and identified limitations of existing methods in recovering the true components. We applied PDR to three clusters of five to six traits genetically correlated with coronary artery disease (CAD), asthma, and type II diabetes (T2D), producing biologically interpretable components. For CAD, PDR identified components related to BMI, hypertension, and cholesterol, and it clarified the relationship among these highly correlated risk factors. We assigned variants to components, calculated their posterior-mean effect sizes, and performed out-of-sample validation. Our posterior-mean effect sizes pool statistical power across traits and substantially boost the correlation (r2) between true and estimated effect sizes (compared with the original summary statistics) by 94% and 70% for asthma and T2D out of sample, respectively, and by a predicted 300% for CAD.