Fast set-based association analysis using summary data from GWAS identifies novel gene loci for human complex traits

Integrating brain imaging features and genomic profiles for the subtyping of major depression

BACKGROUND

Precise stratification of patients into homogeneous disease subgroups could address the heterogeneity of phenotypes and enhance understanding of the pathophysiology underlying specific subtypes. Existing literature on subtyping patients with major depressive disorder (MDD) mainly utilized clinical features only. Genomic and imaging data may improve subtyping, but advanced methods are required due to the high dimensionality of features.

METHODS

We propose a novel disease subtyping framework for MDD by integrating brain structural features, genotype-predicted expression levels in brain tissues, and clinical features. Using a multi-view biclustering approach, we classify patients into clinically and biologically homogeneous subgroups. Additionally, we propose approaches to identify causally relevant genes for clustering.

RESULTS

We verified the reliability of the subtyping model by internal and external validation. High prediction strengths (PS) (average PS: 0.896, minimum: 0.854), a measure of generalizability of the derived clusters in independent datasets, support the validity of our approach. External validation using patient outcome variables (treatment response and hospitalization risks) confirmed the clinical relevance of the identified subgroups. Furthermore, subtype-defining genes overlapped with known susceptibility genes for MDD and were involved in relevant biological pathways. In addition, drug repositioning analysis based on these genes prioritized promising candidates for subtype-specific treatments.

CONCLUSIONS

Our approach successfully stratified MDD patients into subgroups with distinct clinical prognoses. The identification of biologically and clinically meaningful subtypes may enable more personalized treatment strategies. This study also provides a framework for disease subtyping that can be extended to other complex disorders.

Genetic correlation, pleiotropic loci and shared risk genes between major depressive disorder and gastrointestinal tract disorders

BACKGROUND

Major depressive disorder (MDD) is associated with gastrointestinal tract (GIT) disorders, while genetic correlation, pleiotropic loci and shared risk genes remain to be explored.

METHODS

Leveraging genome-wide association study statistics for MDD (n = 170,756), peptic ulcer disease (PUD; n = 16,666), gastroesophageal reflux disease (GORD; n = 54,854), PUD and/or GORD and/or medications (PGM; n = 90,175), irritable bowel syndrome (IBS; n = 28,518), and inflammatory bowel disease (IBD; n = 7045), we determined global and local genetic correlations, identified pleiotropic loci, performed gene-level evaluations, and inferred causal associations using bidirectional Mendelian randomization.

RESULTS

We found global correlation of MDD with PUD (rg = 0.444, P = 3.135 × 10-24), GORD (rg = 0.459, P = 2.568 × 10-65), PGM (rg = 0.498, P = 6.094 × 10-114), IBS (rg = 0.621, P = 2.483 × 10-63), and IBD (rg = 0.171, P = 1.824 × 10-5). We identified 12 locally correlated regions between MDD and GIT disorders except for IBD, and one shared region (chr11:111985737-113,103,996) for PGM, GORD, and IBS. We found one pleiotropic locus for PUD, 12 for GORD, 30 for PGM, eight for IBS, and seven for IBD, and five shared loci (rs138786869, rs2284189, rs3130063, rs35789010, rs7568369) for GORD and PGM. We respectively observed 14 and 20 overlapping genes for MDD-GORD and MDD-PGM. We showed genetic liabilities to GORD, PGM, and IBS causally increase MDD risk, while all reverse causalities are significant.

CONCLUSIONS

Our work identifies genetic architectures shared between MDD and GIT disorders, contributes genetic insights to understand depression in the context of gut-brain interactions, and provides potential targets to treat gastrointestinal symptoms in depressive patients.

Mapping variants in thyroid hormone transporter MCT8 to disease severity by genomic, phenotypic, functional, structural and deep learning integration

Predicting and quantifying phenotypic consequences of genetic variants in rare disorders is a major challenge, particularly pertinent for 'actionable' genes such as thyroid hormone transporter MCT8 (encoded by the X-linked SLC16A2 gene), where loss-of-function (LoF) variants cause a rare neurodevelopmental and (treatable) metabolic disorder in males. The combination of deep phenotyping data with functional and computational tests and with outcomes in population cohorts, enabled us to: (i) identify the genetic aetiology of divergent clinical phenotypes of MCT8 deficiency with genotype-phenotype relationships present across survival and 24 out of 32 disease features; (ii) demonstrate a mild phenocopy in ~400,000 individuals with common genetic variants in MCT8; (iii) assess therapeutic effectiveness, which did not differ among LoF-categories; (iv) advance structural insights in normal and mutated MCT8 by delineating seven critical functional domains; (v) create a pathogenicity-severity MCT8 variant classifier that accurately predicted pathogenicity (AUC:0.91) and severity (AUC:0.86) for 8151 variants. Our information-dense mapping provides a generalizable approach to advance multiple dimensions of rare genetic disorders.

Genomic Characterisation of the Relationship and Causal Links Between Vascular Calcification, Alzheimer's Disease, and Cognitive Traits

Background/Objectives: Observational studies suggest a link between vascular calcification and dementia or cognitive decline, but the evidence is conflicting, and the underlying mechanisms are unclear. Here, we investigate the shared genetic and causal relationships of vascular calcification-coronary artery calcification (CAC) and abdominal aortic calcification (AAC)-with Alzheimer's disease (AD), and five cognitive traits. Methods: We analyse large-scale genome-wide association studies (GWAS) summary statistics, using well-regarded methods, including linkage disequilibrium score regression (LDSC), Mendelian randomisation (MR), pairwise GWAS (GWAS-PW), and gene-based association analysis. Results: Our findings reveal a nominally significant positive genome-wide genetic correlation between CAC and AD, which becomes non-significant after excluding the APOE region. CAC and AAC demonstrate significant negative correlations with cognitive performance and educational attainment. MR found no causal association between CAC or AAC and AD or cognitive traits, except for a bidirectional borderline-significant association between AAC and fluid intelligence scores. Pairwise-GWAS analysis identifies no shared causal SNPs (posterior probability of association [PPA]3 < 0.5). However, we find pleiotropic loci (PPA4 > 0.9), particularly on chromosome 19, with gene association analyses revealing significant genes in shared regions, including APOE, TOMM40, NECTIN2, and APOC1. Moreover, we identify suggestively significant loci (PPA4 > 0.5) on chromosomes 1, 6, 7, 9 and 19, implicating pleiotropic genes, including NAV1, IPO9, PHACTR1, UFL1, FHL5, and FOCAD. Conclusions: Current findings reveal limited genetic correlation and no significant causal associations of CAC and AAC with AD or cognitive traits. However, significant pleiotropic loci, particularly at the APOE region, highlight the complex interplay between vascular calcification and neurodegenerative processes. Given APOE's roles in lipid metabolism, neuroinflammation, and vascular integrity, its involvement may link vascular and neurodegenerative disorders, pointing to potential targets for further investigation.

Multiomic single-cell profiling identifies critical regulators of postnatal brain

Human brain development spans from embryogenesis to adulthood, with dynamic gene expression controlled by cell-type-specific cis-regulatory element activity and three-dimensional genome organization. To advance our understanding of postnatal brain development, we simultaneously profiled gene expression and chromatin accessibility in 101,924 single nuclei from four brain regions across ten donors, covering five key postnatal stages from infancy to late adulthood. Using this dataset and chromosome conformation capture data, we constructed enhancer-based gene regulatory networks to identify cell-type-specific regulators of brain development and interpret genome-wide association study loci for ten main brain disorders. Our analysis connected 2,318 cell-specific loci to 1,149 unique genes, representing 41% of loci linked to the investigated traits, and highlighted 55 genes influencing several disease phenotypes. Pseudotime analysis revealed distinct stages of postnatal oligodendrogenesis and their regulatory programs. These findings provide a comprehensive dataset of cell-type-specific gene regulation at critical timepoints in postnatal brain development.

MARSweb: a fully automated web service for set-based association testing

BACKGROUND

Despite the successes in GWAS, there is still a large gap between the known heritability and the part explained by the SNPs identified by GWAS. Set-based analysis is one of the approaches that has tried to identify associations between multiple variants in a locus a trait, leveraging allelic heterogeneity to increase power in association testing. MARS is a set-based analysis method that integrates likelihood ratio test with a recently developed fine mapping technique to accurately account for causal status of variants in a risk locus. Unfortunately, due to its complex running process, time complexity, and the requirement of high-performance computing resources, it is not widely used.

RESULTS

To address these issues, we proposed a fully automated web-based analysis service, MARSweb. By providing a web service, we minimized the effort required for initial configuration. Additionally, users can perform analyses by simply uploading their data without needing to familiarize themselves with intricate analysis procedures. Furthermore, it facilitates easier interpretation of results by integrating advanced visualization tools. We confirmed the performance of MARSweb by detecting eGenes and performing pathway analysis of the genes using a Yeast Dataset.

CONCLUSIONS

MARSweb is a web-based analysis service that fully automates set-based analysis. It offers an intuitive user interface, making complex analyses more accessible while significantly reducing processing time for enhanced efficiency. MARSweb is available for use at http://cblab.dongguk.edu/MARSweb and its source code is available at https://github.com/DGU-CBLAB/MARSweb .

Gene-based Hardy-Weinberg equilibrium test using genotype count data: application to six types of cancers

BACKGROUND

An alternative approach to investigate associations between genetic variants and disease is to examine deviations from the Hardy-Weinberg equilibrium (HWE) in genotype frequencies within a case population, instead of case-control association analysis. The HWE analysis requires disease cases and demonstrates a notable ability in mapping recessive variants. Allelic heterogeneity is a common phenomenon in diseases. While gene-based case-control association analysis successfully incorporates this heterogeneity, there are no such approaches for HWE analysis. Therefore, we proposed a gene-based HWE test (gene-HWT) by aggregating single-nucleotide polymorphism (SNP)-level HWE test statistics in a gene to address allelic heterogeneity.

RESULTS

This method used only genotype count data and publicly available linkage disequilibrium information and has a very low computational cost. Extensive simulations demonstrated that gene-HWT effectively controls the type I error at a low significance level and outperforms SNP-level HWE test in power when there are multiple causal variants within a gene. Using gene-HWT, we analyzed genotype count data from a genome-wide association study of six cancer types in Japanese individuals and suggest DGKE and ANO3 as potential germline factors in colorectal cancer. Furthermore, FSTL4 was suggested through a combined analysis across the six cancer types, with particularly notable associations observed in colorectal and prostate cancers.

CONCLUSIONS

These findings indicate the potential of gene-HWT to elucidate the genetic basis of complex diseases, including cancer.

Major Depressive Disorder Working Group of the Psychiatric Genomics Consortium. Trans-ancestry genome-wide study of depression identifies 697 associations implicating cell types and pharmacotherapies

In a genome-wide association study (GWAS) meta-analysis of 688,808 individuals with major depression (MD) and 4,364,225 controls from 29 countries across diverse and admixed ancestries, we identify 697 associations at 635 loci, 293 of which are novel. Using fine-mapping and functional tools, we find 308 high-confidence gene associations and enrichment of postsynaptic density and receptor clustering. A neural cell-type enrichment analysis utilizing single-cell data implicates excitatory, inhibitory, and medium spiny neurons and the involvement of amygdala neurons in both mouse and human single-cell analyses. The associations are enriched for antidepressant targets and provide potential repurposing opportunities. Polygenic scores trained using European or multi-ancestry data predicted MD status across all ancestries, explaining up to 5.8% of MD liability variance in Europeans. These findings advance our global understanding of MD and reveal biological targets that may be used to target and develop pharmacotherapies addressing the unmet need for effective treatment.

Genetic Association of Juvenile Idiopathic Arthritis With Adult Rheumatic Disease

Importance

Patients with juvenile idiopathic arthritis (JIA) may develop adult rheumatic diseases later in life, and prolonged or recurrent disease activity is often associated with substantial disability; therefore, it is important to identify patients with JIA at high risk of developing adult rheumatic diseases and provide specialized attention and preventive care to them.

Objective

To elucidate the full extent of the genetic association of JIA with adult rheumatic diseases, to improve treatment strategies and patient outcomes for patients at high risk of developing long-term rheumatic diseases.

Design, Setting, and Participants

In this genetic association study of 4 disease genome-wide association study (GWAS) cohorts from 2013 to 2024 (JIA, rheumatoid arthritis [RA], systemic lupus erythematosus [SLE], and systemic sclerosis [SSc]), patients in the JIA cohort were recruited from the US, Australia, and Norway (with a UK cohort included in the meta-analyzed cohort), while patients in the other 3 cohorts were recruited from US and Western European countries. All analyses were conducted between September 2023 and April 2024.

Exposures

Genetic associations.

Main Outcomes and Measures

Genetic correlations and shared genomic loci between JIA and adult rheumatic diseases. Genetic correlation analyses and cross-trait meta-analysis were conducted on the JIA cohort and the summary statistics of the GWASs from adult rheumatic diseases (RA, SLE, and SSc). Mendelian randomization analyses were also conducted.

Results

This study included 33 207 patients across the 4 cohorts, with 4550 patients in the meta-analyzed JIA cohort (JIA cohort: 1485 patients with arthritis onset before 16 years; 1017 female [68.5%]; 10 352 controls; UK cohort: 3305 patients with JIA; 9196 controls), 143 61 patients in the RA cohort, 5201 patients in the SLE cohort; and 9095 patients in the SSc cohort. After the GWAS result of the JIA cohort was meta-analyzed with the UK JIA cohort, there was a total of 4550 JIA cases and 18 446 controls. The analysis revealed a significant global correlation between JIA and adult rheumatic diseases, with 84 regions harboring signals associated with multiple diseases. Cross-trait analyses uncovered novel disease loci and 20 loci associated with JIA and adult diseases. Mendelian randomization analysis revealed the significant association of 11 proteins with rheumatic disorders. Both shared, organ-specific, and disease-specific critical cell types were highlighted.

Conclusions and Relevance

In this genetic association study, there was significant genetic overlap between JIA and adult rheumatic diseases. These findings may help to refine JIA classification, risk stratification, and therapeutic strategy of repurposing adult disease drugs for pediatric patients with similar mechanisms.

Single-cell multi-modal integrative analyses highlight functional dynamic gene regulatory networks directing human cardiac development

Illuminating the precise stepwise genetic programs directing cardiac development provides insights into the mechanisms of congenital heart disease and strategies for cardiac regenerative therapies. Here, we integrate in vitro and in vivo human single-cell multi-omic studies with high-throughput functional genomic screening to reveal dynamic, cardiac-specific gene regulatory networks (GRNs) and transcriptional regulators during human cardiomyocyte development. Interrogating developmental trajectories reconstructed from single-cell data unexpectedly reveal divergent cardiomyocyte lineages with distinct gene programs based on developmental signaling pathways. High-throughput functional genomic screens identify key transcription factors from inferred GRNs that are functionally relevant for cardiomyocyte lineages derived from each pathway. Notably, we discover a critical heat shock transcription factor 1 (HSF1)-mediated cardiometabolic GRN controlling cardiac mitochondrial/metabolic function and cell survival, also observed in fetal human cardiomyocytes. Overall, these multi-modal genomic studies enable the systematic discovery and validation of coordinated GRNs and transcriptional regulators controlling the development of distinct human cardiomyocyte populations.

Genome-wide meta-analysis conducted in three large biobanks expands the genetic landscape of lumbar disc herniations

Given that lumbar disc herniation (LDH) is a prevalent spinal condition that causes significant individual suffering and societal costs, the genetic basis of LDH has received relatively little research. Our aim is to increase understanding of the genetic factors influencing LDH. We perform a genome-wide association analysis (GWAS) of LDH in the FinnGen project and in Estonian and UK biobanks, followed by a genome-wide meta-analysis to combine the results. In the meta-analysis, we identify 41 loci that have not been associated with LDH in prior studies on top of the 23 known risk loci. We detect LDH-associated loci in the vicinity of genes related to inflammation, disc-related structures, and synaptic transmission. Overall, our research contributes to a deeper understanding of the genetic factors behind LDH, potentially paving the way for the development of new therapeutics, prevention methods, and treatments for symptomatic LDH in the future.

Unraveling metabolic stress response in dairy cows: Genetic control of plasma biomarkers throughout lactation and the transition period

Breeding animals able to effectively respond to stress could be a long-term, sustainable, and affordable strategy to improve resilience and welfare in livestock systems. In the present study, the concentrations of 29 plasma biomarkers were used as candidate endophenotypes for metabolic stress response in single-SNP, gene- and haplotype-based GWAS using 739 healthy lactating Italian Holstein cows and 88,271 variants. Significant genetic associations were found in all the 3 GWAS approaches for plasma γ-glutamyl transferase concentration on BTA17, for paraoxonase on BTA4, and for alkaline phosphatase and zinc on BTA2. On these chromosomes, single-SNP and gene-based chromosome-wide association studies were performed, confirming GWAS findings. The signals identified for paraoxonase, γ-glutamyl transferase, and alkaline phosphatase were in proximity to the genes coding for them. The heritability of these 4 biomarkers ranged from moderate to high (from 0.39 to 0.54). Plasma biomarkers are known to undergo large changes in concentration during metabolic stress in the transition period, with an interindividual variability in the rate of change and recovery time. Genetics may account in part for these differences. To assess this, we studied a subset of 139 periparturient cows homozygous at 3 SNPs known to be respectively associated with concentration of plasma ceruloplasmin, paraoxonase, and γ-glutamyl transferase. We compared the immune-metabolic profile measured in plasma at -7, +5, and +30 d relative to calving between groups of opposite homozygotes. A significant effect of the genotype was found on paraoxonase and γ-glutamyl transferase plasma concentration at all the 3 time points. No evidence for genotype effect was detected for ceruloplasmin. Understanding the genetic control underlying metabolic stress response may suggest new approaches to foster resilience in dairy cows.

Investigating genotype by environment interaction for beef cattle fertility traits in commercial herds in northern Australia with multi-trait analysis

BACKGROUND

Genotype by environment interactions (GxE) affect a range of production traits in beef cattle. Quantifying the effect of GxE in commercial and multi-breed herds is challenging due to unknown genetic linkage between animals across environment levels. The primary aim of this study was to use multi-trait models to investigate GxE for three heifer fertility traits, corpus luteum (CL) presence, first pregnancy and second pregnancy, in a large tropical beef multibreed dataset (n = 21,037). Environmental levels were defined by two different descriptors, burden of heat load (temperature humidity index, THI) and nutritional availability (based on mean average daily gain for the herd, ADWG). To separate the effects of genetic linkage and real GxE across the environments, 1000 replicates of a simulated phenotype were generated by simulating QTL effects with no GxE onto real marker genotypes from the population, to determine the genetic correlations that could be expected across environments due to the existing genetic linkage only. Correlations from the real phenotypes were then compared to the empirical distribution under the null hypothesis from the simulated data. By adopting this approach, this study attempted to establish if low genetic correlations between environmental levels were due to GxE or insufficient genetic linkage between animals in each environmental level.

RESULTS

The correlations (being less than <0.8) for the real phenotypes were indicative of GxE for CL presence between ADWG environmental levels and in pregnancy traits. However, none of the correlations for CL presence or first pregnancy between ADWG levels were below the 5th percentile value for the empirical distribution under the null hypothesis from the simulated data. Only one statistically significant (P < 0.05) indication of GxE for first pregnancy was found between THI environmental levels, where rg = 0.28 and 5th percentile value = 0.29, and this result was marginal.

CONCLUSIONS

Only one case of statistically significant GxE for fertility traits was detected for first pregnancy between THI environmental levels 2 and 3. Other initial indications of GxE that were observed from the real phenotypes did not prove significant when compared to an empirical null distribution from simulated phenotypes. The lack of compelling evidence of GxE indicates that direct selection for fertility traits can be made accurately, using a single evaluation, regardless of environment.

The goldmine of GWAS summary statistics: a systematic review of methods and tools

Genome-wide association studies (GWAS) have revolutionized our understanding of the genetic architecture of complex traits and diseases. GWAS summary statistics have become essential tools for various genetic analyses, including meta-analysis, fine-mapping, and risk prediction. However, the increasing number of GWAS summary statistics and the diversity of software tools available for their analysis can make it challenging for researchers to select the most appropriate tools for their specific needs. This systematic review aims to provide a comprehensive overview of the currently available software tools and databases for GWAS summary statistics analysis. We conducted a comprehensive literature search to identify relevant software tools and databases. We categorized the tools and databases by their functionality, including data management, quality control, single-trait analysis, and multiple-trait analysis. We also compared the tools and databases based on their features, limitations, and user-friendliness. Our review identified a total of 305 functioning software tools and databases dedicated to GWAS summary statistics, each with unique strengths and limitations. We provide descriptions of the key features of each tool and database, including their input/output formats, data types, and computational requirements. We also discuss the overall usability and applicability of each tool for different research scenarios. This comprehensive review will serve as a valuable resource for researchers who are interested in using GWAS summary statistics to investigate the genetic basis of complex traits and diseases. By providing a detailed overview of the available tools and databases, we aim to facilitate informed tool selection and maximize the effectiveness of GWAS summary statistics analysis.

Sex effects on DNA methylation affect discovery in epigenome-wide association study of schizophrenia

Sex differences in the epidemiology and clinical characteristics of schizophrenia are well-known; however, the molecular mechanisms underlying these differences remain unclear. Further, the potential advantages of sex-stratified meta-analyses of epigenome-wide association studies (EWAS) of schizophrenia have not been investigated. Here, we performed sex-stratified EWAS meta-analyses to investigate whether sex stratification improves discovery, and to identify differentially methylated regions (DMRs) in schizophrenia. Peripheral blood-derived DNA methylation data from 1519 cases of schizophrenia (male n = 989, female n = 530) and 1723 controls (male n = 997, female n = 726) from three publicly available datasets, and the TOP cohort were meta-analyzed to compare sex-specific, sex-stratified, and sex-adjusted EWAS. The predictive power of each model was assessed by polymethylation score (PMS). The number of schizophrenia-associated differentially methylated positions identified was higher for the sex-stratified model than for the sex-adjusted one. We identified 20 schizophrenia-associated DMRs in the sex-stratified analysis. PMS from sex-stratified analysis outperformed that from sex-adjusted analysis in predicting schizophrenia. Notably, PMSs from the sex-stratified and female-only analyses, but not those from sex-adjusted or the male-only analyses, significantly predicted schizophrenia in males. The findings suggest that sex-stratified EWAS meta-analyses improve the identification of schizophrenia-associated epigenetic changes and highlight an interaction between sex and schizophrenia status on DNA methylation. Sex-specific DNA methylation may have potential implications for precision psychiatry and the development of stratified treatments for schizophrenia.

From genetic associations to genes: methods, applications, and challenges

Genome-wide association studies (GWASs) have identified numerous genetic loci associated with human traits and diseases. However, pinpointing the causal genes remains a challenge, which impedes the translation of GWAS findings into biological insights and medical applications. In this review, we provide an in-depth overview of the methods and technologies used for prioritizing genes from GWAS loci, including gene-based association tests, integrative analysis of GWAS and molecular quantitative trait loci (xQTL) data, linking GWAS variants to target genes through enhancer-gene connection maps, and network-based prioritization. We also outline strategies for generating context-dependent xQTL data and their applications in gene prioritization. We further highlight the potential of gene prioritization in drug repurposing. Lastly, we discuss future challenges and opportunities in this field.

Abdominal aortic aneurysm and cardiometabolic traits share strong genetic susceptibility to lipid metabolism and inflammation

Abdominal aortic aneurysm has a high heritability and often co-occurs with other cardiometabolic disorders, suggesting shared genetic susceptibility. We investigate this commonality leveraging recent GWAS studies of abdominal aortic aneurysm and 32 cardiometabolic traits. We find significant genetic correlations between abdominal aortic aneurysm and 21 of the cardiometabolic traits investigated, including causal relationships with coronary artery disease, hypertension, lipid traits, and blood pressure. For each trait pair, we identify shared causal variants, genes, and pathways, revealing that cholesterol metabolism and inflammation are shared most prominently. Additionally, we show the tissue and cell type specificity in the shared signals, with strong enrichment across traits in the liver, arteries, adipose tissues, macrophages, adipocytes, and fibroblasts. Finally, we leverage drug-gene databases to identify several lipid-lowering drugs and antioxidants with high potential to treat abdominal aortic aneurysm with comorbidities. Our study provides insight into the shared genetic mechanism between abdominal aortic aneurysm and cardiometabolic traits, and identifies potential targets for pharmacological intervention.

Identifying BMI-associated genes via a genome-wide multi-omics integrative approach using summary data

OBJECTIVE

This study aims to identify BMI-associated genes by integrating aggregated summary information from different omics data.

METHODS

We conducted a meta-analysis to leverage information from a genome-wide association study (n = 339 224), a transcriptome-wide association study (n = 5619), and an epigenome-wide association study (n = 3743). We prioritized the significant genes with a machine learning-based method, netWAS, which borrows information from adipose tissue-specific interaction networks. We also used the brain-specific network in netWAS to investigate genes potentially involved in brain-adipose interaction.

RESULTS

We identified 195 genes that were significantly associated with BMI through meta-analysis. The netWAS analysis narrowed down the list to 21 genes in adipose tissue. Among these 21 genes, six genes, including FUS, STX4, CCNT2, FUBP1, NDUFS3, and RAPSN, were not reported to be BMI-associated in PubMed or GWAS Catalog. We also identified 11 genes that were significantly associated with BMI in both adipose and whole brain tissues.

CONCLUSION

This study integrated three types of omics data and identified a group of genes that have not previously been reported to be associated with BMI. This strategy could provide new insights for future studies to identify molecular mechanisms contributing to BMI regulation.

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.

Network propagation for GWAS analysis: a practical guide to leveraging molecular networks for disease gene discovery

MOTIVATION

Genome-wide association studies (GWAS) have enabled large-scale analysis of the role of genetic variants in human disease. Despite impressive methodological advances, subsequent clinical interpretation and application remains challenging when GWAS suffer from a lack of statistical power. In recent years, however, the use of information diffusion algorithms with molecular networks has led to fruitful insights on disease genes.

RESULTS

We present an overview of the design choices and pitfalls that prove crucial in the application of network propagation methods to GWAS summary statistics. We highlight general trends from the literature, and present benchmark experiments to expand on these insights selecting as case study three diseases and five molecular networks. We verify that the use of gene-level scores based on GWAS P-values offers advantages over the selection of a set of 'seed' disease genes not weighted by the associated P-values if the GWAS summary statistics are of sufficient quality. Beyond that, the size and the density of the networks prove to be important factors for consideration. Finally, we explore several ensemble methods and show that combining multiple networks may improve the network propagation approach.

HOPX-associated molecular programs control cardiomyocyte cell states underpinning cardiac structure and function

Genomic regulation of cardiomyocyte differentiation is central to heart development and function. This study uses genetic loss-of-function human-induced pluripotent stem cell-derived cardiomyocytes to evaluate the genomic regulatory basis of the non-DNA-binding homeodomain protein HOPX. We show that HOPX interacts with and controls cardiac genes and enhancer networks associated with diverse aspects of heart development. Using perturbation studies in vitro, we define how upstream cell growth and proliferation control HOPX transcription to regulate cardiac gene programs. We then use cell, organoid, and zebrafish regeneration models to demonstrate that HOPX-regulated gene programs control cardiomyocyte function in development and disease. Collectively, this study mechanistically links cell signaling pathways as upstream regulators of HOPX transcription to control gene programs underpinning cardiomyocyte identity and function.

Effects of genetics on odor perception: Can a quick smell test effectively screen everyone?

SCENTinel, a rapid smell test designed to screen for olfactory disorders, including anosmia (no ability to smell an odor) and parosmia (distorted sense of smell), measures 4 components of olfactory function: detection, intensity, identification, and pleasantness. Each test card contains one of 9 odorant mixtures. Some people born with genetic insensitivities to specific odorants (i.e. specific anosmia) may fail the test if they cannot smell an odorant but otherwise have a normal sense of smell. However, using odorant mixtures has largely been found to prevent this from happening. To better understand whether genetic differences affect SCENTinel test results, we asked genetically informative adult participants (twins or triplets, N = 630; singletons, N = 370) to complete the SCENTinel test. A subset of twins (n = 304) also provided a saliva sample for genotyping. We examined data for differences between the 9 possible SCENTinel odors; effects of age, sex, and race on SCENTinel performance, test-retest variability; and heritability using both structured equation modeling and SNP-based statistical methods. None of these strategies provided evidence for specific anosmia for any of the odors, but ratings of pleasantness were, in part, genetically determined (h2 = 0.40) and were nominally associated with alleles of odorant receptors (e.g. OR2T33 and OR1G1; P < 0.001). These results provide evidence that using odorant mixtures protected against effects of specific anosmia for ratings of intensity but that ratings of pleasantness showed effects of inheritance, possibly informed by olfactory receptor genotypes.

Sex bias in celiac disease: XWAS and monocyte eQTLs in women identify TMEM187 as a functional candidate gene

BACKGROUND

Celiac disease (CeD) is an immune-mediated disorder that develops in genetically predisposed individuals upon gluten consumption. HLA risk alleles explain 40% of the genetic component of CeD, so there have been continuing efforts to uncover non-HLA loci that can explain the remaining heritability. As in most autoimmune disorders, the prevalence of CeD is significantly higher in women. Here, we investigated the possible involvement of the X chromosome on the sex bias of CeD.

METHODS

We performed a X chromosome-wide association study (XWAS) and a gene-based association study in women from the CeD Immunochip (7062 cases, 5446 controls). We also constructed a database of X chromosome cis-expression quantitative trait loci (eQTLs) in monocytes from unstimulated (n = 226) and lipopolysaccharide (LPS)-stimulated (n = 130) female donors and performed a Summary-data-based MR (SMR) analysis to integrate XWAS and eQTL information. We interrogated the expression of the potentially causal gene (TMEM187) in peripheral blood mononuclear cells (PBMCs) from celiac patients at onset, on a gluten-free diet, potential celiac patients and non-celiac controls.

RESULTS

The XWAS and gene-based analyses identified 13 SNPs and 25 genes, respectively, 22 of which had not been previously associated with CeD. The X chromosome cis-eQTL analysis found 18 genes with at least one cis-eQTL in naïve female monocytes and 8 genes in LPS-stimulated female monocytes, 2 of which were common to both situations and 6 were unique to LPS stimulation. SMR identified a potentially causal association of TMEM187 expression in naïve monocytes with CeD in women, regulated by CeD-associated, eQTL-SNPs rs7350355 and rs5945386. The CeD-risk alleles were correlated with lower TMEM187 expression. These results were replicated using eQTLs from LPS-stimulated monocytes. We observed higher levels of TMEM187 expression in PBMCs from female CeD patients at onset compared to female non-celiac controls, but not in male CeD individuals.

CONCLUSION

Using X chromosome genotypes and gene expression data from female monocytes, SMR has identified TMEM187 as a potentially causal candidate in CeD. Further studies are needed to understand the implication of the X chromosome in the higher prevalence of CeD in women.

Biliary atresia is associated with polygenic susceptibility in ciliogenesis and planar polarity effector genes

BACKGROUND & AIMS

Biliary atresia (BA) is poorly understood and leads to liver transplantation (LT), with the requirement for and associated risks of lifelong immunosuppression, in most children. We performed a genome-wide association study (GWAS) to determine the genetic basis of BA.

METHODS

We performed a GWAS in 811 European BA cases treated with LT in US, Canadian and UK centers, and 4,654 genetically matched controls. Whole-genome sequencing of 100 cases evaluated synthetic association with rare variants. Functional studies included whole liver transcriptome analysis of 64 BA cases and perturbations in experimental models.

RESULTS

A GWAS of common single nucleotide polymorphisms (SNPs), i.e. allele frequencies >1%, identified intronic SNPs rs6446628 in AFAP1 with genome-wide significance (p = 3.93E-8) and rs34599046 in TUSC3 at sub-threshold genome-wide significance (p = 1.34E-7), both supported by credible peaks of neighboring SNPs. Like other previously reported BA-associated genes, AFAP1 and TUSC3 are ciliogenesis and planar polarity effectors (CPLANE). In gene-set-based GWAS, BA was associated with 6,005 SNPs in 102 CPLANE genes (p = 5.84E-15). Compared with non-CPLANE genes, more CPLANE genes harbored rare variants (allele frequency <1%) that were assigned Human Phenotype Ontology terms related to hepatobiliary anomalies by predictive algorithms, 87% vs. 40%, p <0.0001. Rare variants were present in multiple genes distinct from those with BA-associated common variants in most BA cases. AFAP1 and TUSC3 knockdown blocked ciliogenesis in mouse tracheal cells. Inhibition of ciliogenesis caused biliary dysgenesis in zebrafish. AFAP1 and TUSC3 were expressed in fetal liver organoids, as well as fetal and BA livers, but not in normal or disease-control livers. Integrative analysis of BA-associated variants and liver transcripts revealed abnormal vasculogenesis and epithelial tube formation, explaining portal vein anomalies that co-exist with BA.

CONCLUSIONS

BA is associated with polygenic susceptibility in CPLANE genes. Rare variants contribute to polygenic risk in vulnerable pathways via unique genes.

IMPACT AND IMPLICATIONS

Liver transplantation is needed to cure most children born with biliary atresia, a poorly understood rare disease. Transplant immunosuppression increases the likelihood of life-threatening infections and cancers. To improve care by preventing this disease and its progression to transplantation, we examined its genetic basis. We find that this disease is associated with both common and rare mutations in highly specialized genes which maintain normal communication and movement of cells, and their organization into bile ducts and blood vessels during early development of the human embryo. Because defects in these genes also cause other birth defects, our findings could lead to preventive strategies to lower the incidence of biliary atresia and potentially other birth defects.

Investigating the shared genetic architecture of post-traumatic stress disorder and gastrointestinal tract disorders: a genome-wide cross-trait analysis

BACKGROUND

Observational studies suggest a correlation between post-traumatic stress disorder (PTSD) and gastrointestinal tract (GIT) disorders. However, the genetic overlap, causal relationships, and underlining mechanisms between PTSD and GIT disorders were absent.

METHODS

We obtained genome-wide association study statistics for PTSD (23 212 cases, 151 447 controls), peptic ulcer disease (PUD; 16 666 cases, 439 661 controls), gastroesophageal reflux disease (GORD; 54 854 cases, 401 473 controls), PUD and/or GORD and/or medications (PGM; 90 175 cases, 366 152 controls), irritable bowel syndrome (IBS; 28 518 cases, 426 803 controls), and inflammatory bowel disease (IBD; 7045 cases, 449 282 controls). We quantified genetic correlations, identified pleiotropic loci, and performed multi-marker analysis of genomic annotation, fast gene-based association analysis, transcriptome-wide association study analysis, and bidirectional Mendelian randomization analysis.

RESULTS

PTSD globally correlates with PUD (rg = 0.526, p = 9.355 × 10-7), GORD (rg = 0.398, p = 5.223 × 10-9), PGM (rg = 0.524, p = 1.251 × 10-15), and IBS (rg = 0.419, p = 8.825 × 10-6). Cross-trait meta-analyses identify seven genome-wide significant loci between PTSD and PGM (rs13107325, rs1632855, rs1800628, rs2188100, rs3129953, rs6973700, and rs73154693); three between PTSD and GORD (rs13107325, rs1632855, and rs3132450); one between PTSD and IBS/IBD (rs4937872 and rs114969413, respectively). Proximal pleiotropic genes are mainly enriched in immune response regulatory pathways, and in brain, digestive, and immune systems. Gene-level analyses identify five candidates: ABT1, BTN3A2, HIST1H3J, ZKSCAN4, and ZKSCAN8. We found significant causal effects of GORD, PGM, IBS, and IBD on PTSD. We observed no reverse causality of PTSD with GIT disorders, except for GORD.

CONCLUSIONS

PTSD and GIT disorders share common genetic architectures. Our work offers insights into the biological mechanisms, and provides genetic basis for translational research studies.

Multi-omics cannot replace sample size in genome-wide association studies

The integration of multi-omics information (e.g., epigenetics and transcriptomics) can be useful for interpreting findings from genome-wide association studies (GWAS). It has been suggested that multi-omics could circumvent or greatly reduce the need to increase GWAS sample sizes for novel variant discovery. We tested whether incorporating multi-omics information in earlier and smaller-sized GWAS boosts true-positive discovery of genes that were later revealed by larger GWAS of the same/similar traits. We applied 10 different analytic approaches to integrating multi-omics data from 12 sources (e.g., Genotype-Tissue Expression project) to test whether earlier and smaller GWAS of 4 brain-related traits (alcohol use disorder/problematic alcohol use, major depression/depression, schizophrenia, and intracranial volume/brain volume) could detect genes that were revealed by a later and larger GWAS. Multi-omics data did not reliably identify novel genes in earlier less-powered GWAS (PPV <0.2; 80% false-positive associations). Machine learning predictions marginally increased the number of identified novel genes, correctly identifying 1-8 additional genes, but only for well-powered early GWAS of highly heritable traits (i.e., intracranial volume and schizophrenia). Although multi-omics, particularly positional mapping (i.e., fastBAT, MAGMA, and H-MAGMA), can help to prioritize genes within genome-wide significant loci (PPVs = 0.5-1.0) and translate them into information about disease biology, it does not reliably increase novel gene discovery in brain-related GWAS. To increase power for discovery of novel genes and loci, increasing sample size is required.

Evaluating 17 methods incorporating biological function with GWAS summary statistics to accelerate discovery demonstrates a tradeoff between high sensitivity and high positive predictive value

Where sufficiently large genome-wide association study (GWAS) samples are not currently available or feasible, methods that leverage increasing knowledge of the biological function of variants may illuminate discoveries without increasing sample size. We comprehensively evaluated 17 functional weighting methods for identifying novel associations. We assessed the performance of these methods using published results from multiple GWAS waves across each of five complex traits. Although no method achieved both high sensitivity and positive predictive value (PPV) for any trait, a subset of methods utilizing pleiotropy and expression quantitative trait loci nominated variants with high PPV (>75%) for multiple traits. Application of functionally weighting methods to enhance GWAS power for locus discovery is unlikely to circumvent the need for larger sample sizes in truly underpowered GWAS, but these results suggest that applying functional weighting to GWAS can accurately nominate additional novel loci from available samples for follow-up studies.

Insights from Mendelian randomization and genetic correlation analyses into the relationship between endometriosis and its comorbidities

BACKGROUND

Endometriosis remains a poorly understood disease, despite its high prevalence and debilitating symptoms. The overlap in symptoms and the increased risk of multiple other traits in women with endometriosis is becoming increasingly apparent through epidemiological data. Genetic studies offer a method of investigating these comorbid relationships through the assessment of causal relationships with Mendelian randomization (MR), as well as identification of shared genetic variants and genes involved across traits. This has the capacity to identify risk factors for endometriosis as well as provide insight into the aetiology of disease.

OBJECTIVE AND RATIONALE

We aim to review the current literature assessing the relationship between endometriosis and other traits using genomic data, primarily through the methods of MR and genetic correlation. We critically examine the limitations of these studies in accordance with the assumptions of the utilized methods.

SEARCH METHODS

The PubMed database was used to search for peer-reviewed original research articles using the terms 'Mendelian randomization endometriosis' and '"genetic correlation" endometriosis'. Additionally, a Google Scholar search using the terms '"endometriosis" "mendelian randomization" "genetic correlation"' was performed. All relevant publications (n = 21) published up until 7 October 2022 were included in this review. Upon compilation of all traits with published MR and/or genetic correlation with endometriosis, additional epidemiological and genetic information on their comorbidity with endometriosis was sourced by searching for the trait in conjunction with 'endometriosis' on Google Scholar.

OUTCOMES

The association between endometriosis and multiple pain, gynaecological, cancer, inflammatory, gastrointestinal, psychological, and anthropometric traits has been assessed using MR analysis and genetic correlation analysis. Genetic correlation analyses provide evidence that genetic factors contributing to endometriosis are shared with multiple traits: migraine, uterine fibroids, subtypes of ovarian cancer, melanoma, asthma, gastro-oesophageal reflux disease, gastritis/duodenitis, and depression, suggesting the involvement of multiple biological mechanisms in endometriosis. The assessment of causality with MR has revealed several potential causes (e.g. depression) and outcomes (e.g. ovarian cancer and uterine fibroids) of a genetic predisposition to endometriosis; however, interpretation of these results requires consideration of potential violations of the MR assumptions.

WIDER IMPLICATIONS

Genomic studies have demonstrated that there is a molecular basis for the co-occurrence of endometriosis with other traits. Dissection of this overlap has identified shared genes and pathways, which provide insight into the biology of endometriosis. Thoughtful MR studies are necessary to ascertain causality of the comorbidities of endometriosis. Given the significant diagnostic delay of endometriosis of 7-11 years, determining risk factors is necessary to aid diagnosis and reduce the disease burden. Identification of traits for which endometriosis is a risk factor is important for holistic treatment and counselling of the patient. The use of genomic data to disentangle the overlap of endometriosis with other traits has provided insights into the aetiology of endometriosis.

Genomic characterisation of the overlap of endometriosis with 76 comorbidities identifies pleiotropic and causal mechanisms underlying disease risk

Comorbid conditions can be driven by underlying pleiotropic and causal mechanisms that can provide insights into shared molecular and biological processes contributing to disease risk. Endometriosis is a chronic condition affecting one in nine women of reproductive age and poses many challenges including lengthy diagnostic delays and limited treatment efficacy owing to poor understanding of disease aetiology. To shed light on the underlying biological mechanisms and to identify potential risk factors, we examine the epidemiological and genomic relationship between endometriosis and its comorbidities. In the UK Biobank 292 ICD10 codes were epidemiologically correlated with endometriosis diagnosis, including gynaecological, immune, infection, pain, psychiatric, cancer, gastrointestinal, urinary, bone and cardiovascular traits. A subset of the identified comorbidities (n = 76) underwent follow-up genetic analysis. Whilst Mendelian randomisation suggested causality was not responsible for most comorbid relationships, 22 traits were genetically correlated with endometriosis, including pain, gynaecological and gastrointestinal traits, suggestive of a shared genetic background. Pleiotropic genetic variants and genes were identified using gene-based and colocalisation analysis. Shared genetic risk factors and potential target genes suggest a diverse collection of biological systems are involved in these comorbid relationships including coagulation factors, development of the female reproductive tract and cell proliferation. These findings highlight the diversity of traits with epidemiological and genomic overlap with endometriosis and implicate a key role for pleiotropy in the comorbid relationships.

Haplotype-Based Genome-Wide Association Analysis Using Exome Capture Assay and Digital Phenotyping Identifies Genetic Loci Underlying Salt Tolerance Mechanisms in Wheat

Soil salinity can impose substantial stress on plant growth and cause significant yield losses. Crop varieties tolerant to salinity stress are needed to sustain yields in saline soils. This requires effective genotyping and phenotyping of germplasm pools to identify novel genes and QTL conferring salt tolerance that can be utilised in crop breeding schemes. We investigated a globally diverse collection of 580 wheat accessions for their growth response to salinity using automated digital phenotyping performed under controlled environmental conditions. The results show that digitally collected plant traits, including digital shoot growth rate and digital senescence rate, can be used as proxy traits for selecting salinity-tolerant accessions. A haplotype-based genome-wide association study was conducted using 58,502 linkage disequilibrium-based haplotype blocks derived from 883,300 genome-wide SNPs and identified 95 QTL for salinity tolerance component traits, of which 54 were novel and 41 overlapped with previously reported QTL. Gene ontology analysis identified a suite of candidate genes for salinity tolerance, some of which are already known to play a role in stress tolerance in other plant species. This study identified wheat accessions that utilise different tolerance mechanisms and which can be used in future studies to investigate the genetic and genic basis of salinity tolerance. Our results suggest salinity tolerance has not arisen from or been bred into accessions from specific regions or groups. Rather, they suggest salinity tolerance is widespread, with small-effect genetic variants contributing to different levels of tolerance in diverse, locally adapted germplasm.

Complex genetic architecture underlying the plasticity of maize agronomic traits

Phenotypic plasticity is the ability of a given genotype to produce multiple phenotypes in response to changing environmental conditions. Understanding the genetic basis of phenotypic plasticity and establishing a predictive model is highly relevant to future agriculture under a changing climate. Here we report findings on the genetic basis of phenotypic plasticity for 23 complex traits using a diverse maize population planted at five sites with distinct environmental conditions. We found that latitude-related environmental factors were the main drivers of across-site variation in flowering time traits but not in plant architecture or yield traits. For the 23 traits, we detected 109 quantitative trait loci (QTLs), 29 for mean values, 66 for plasticity, and 14 for both parameters, and 80% of the QTLs interacted with latitude. The effects of several QTLs changed in magnitude or sign, driving variation in phenotypic plasticity. We experimentally validated one plastic gene, ZmTPS14.1, whose effect was likely mediated by the compensation effect of ZmSPL6 from a downstream pathway. By integrating genetic diversity, environmental variation, and their interaction into a joint model, we could provide site-specific predictions with increased accuracy by as much as 9.9%, 2.2%, and 2.6% for days to tassel, plant height, and ear weight, respectively. This study revealed a complex genetic architecture involving multiple alleles, pleiotropy, and genotype-by-environment interaction that underlies variation in the mean and plasticity of maize complex traits. It provides novel insights into the dynamic genetic architecture of agronomic traits in response to changing environments, paving a practical way toward precision agriculture.

StocSum: stochastic summary statistics for whole genome sequencing studies

Genomic summary statistics, usually defined as single-variant test results from genome-wide association studies, have been widely used to advance the genetics field in a wide range of applications. Applications that involve multiple genetic variants also require their correlations or linkage disequilibrium (LD) information, often obtained from an external reference panel. In practice, it is usually difficult to find suitable external reference panels that represent the LD structure for underrepresented and admixed populations, or rare genetic variants from whole genome sequencing (WGS) studies, limiting the scope of applications for genomic summary statistics. Here we introduce StocSum, a novel reference-panel-free statistical framework for generating, managing, and analyzing stochastic summary statistics using random vectors. We develop various downstream applications using StocSum including single-variant tests, conditional association tests, gene-environment interaction tests, variant set tests, as well as meta-analysis and LD score regression tools. We demonstrate the accuracy and computational efficiency of StocSum using two cohorts from the Trans-Omics for Precision Medicine Program. StocSum will facilitate sharing and utilization of genomic summary statistics from WGS studies, especially for underrepresented and admixed populations.

Vascular cells improve functionality of human cardiac organoids

Crosstalk between cardiac cells is critical for heart performance. Here we show that vascular cells within human cardiac organoids (hCOs) enhance their maturation, force of contraction, and utility in disease modeling. Herein we optimize our protocol to generate vascular populations in addition to epicardial, fibroblast, and cardiomyocyte cells that self-organize into in-vivo-like structures in hCOs. We identify mechanisms of communication between endothelial cells, pericytes, fibroblasts, and cardiomyocytes that ultimately contribute to cardiac organoid maturation. In particular, (1) endothelial-derived LAMA5 regulates expression of mature sarcomeric proteins and contractility, and (2) paracrine platelet-derived growth factor receptor β (PDGFRβ) signaling from vascular cells upregulates matrix deposition to augment hCO contractile force. Finally, we demonstrate that vascular cells determine the magnitude of diastolic dysfunction caused by inflammatory factors and identify a paracrine role of endothelin driving dysfunction. Together this study highlights the importance and role of vascular cells in organoid models.

Discovery of genomic loci associated with sleep apnea risk through multi-trait GWAS analysis with snoring

STUDY OBJECTIVES

Despite its association with severe health conditions, the etiology of sleep apnea (SA) remains understudied. This study sought to identify genetic variants robustly associated with SA risk.

METHODS

We performed a genome-wide association study (GWAS) meta-analysis of SA across five cohorts (NTotal = 523 366), followed by a multi-trait analysis of GWAS (multi-trait analysis of genome-wide association summary statistics [MTAG]) to boost power, leveraging the high genetic correlation between SA and snoring. We then adjusted our results for the genetic effects of body mass index (BMI) using multi-trait-based conditional and joint analysis (mtCOJO) and sought replication of lead hits in a large cohort of participants from 23andMe, Inc (NTotal = 1 477 352; Ncases = 175 522). We also explored genetic correlations with other complex traits and performed a phenome-wide screen for causally associated phenotypes using the latent causal variable method.

RESULTS

Our SA meta-analysis identified five independent variants with evidence of association beyond genome-wide significance. After adjustment for BMI, only one genome-wide significant variant was identified. MTAG analyses uncovered 49 significant independent loci associated with SA risk. Twenty-nine variants were replicated in the 23andMe GWAS adjusting for BMI. We observed genetic correlations with several complex traits, including multisite chronic pain, diabetes, eye disorders, high blood pressure, osteoarthritis, chronic obstructive pulmonary disease, and BMI-associated conditions.

CONCLUSION

Our study uncovered multiple genetic loci associated with SA risk, thus increasing our understanding of the etiology of this condition and its relationship with other complex traits.

LDAK-GBAT: Fast and powerful gene-based association testing using summary statistics

We present LDAK-GBAT, a tool for gene-based association testing using summary statistics from genome-wide association studies that is computationally efficient, produces well-calibrated p values, and is significantly more powerful than existing tools. LDAK-GBAT takes approximately 30 min to analyze imputed data (2.9M common, genic SNPs), requiring less than 10 Gb memory. It shows good control of type 1 error given an appropriate reference panel. Across 109 phenotypes (82 from the UK Biobank, 18 from the Million Veteran Program, and nine from the Psychiatric Genetics Consortium), LDAK-GBAT finds on average 19% (SE: 1%) more significant genes than the existing tool sumFREGAT-ACAT, with even greater gains in comparison with MAGMA, GCTA-fastBAT, sumFREGAT-SKAT-O, and sumFREGAT-PCA.

mBAT-combo: A more powerful test to detect gene-trait associations from GWAS data

Gene-based association tests aggregate multiple SNP-trait associations into sets defined by gene boundaries and are widely used in post-GWAS analysis. A common approach for gene-based tests is to combine SNPs associations by computing the sum of χ² statistics. However, this strategy ignores the directions of SNP effects, which could result in a loss of power for SNPs with masking effects, e.g., when the product of two SNP effects and the linkage disequilibrium (LD) correlation is negative. Here, we introduce "mBAT-combo," a set-based test that is better powered than other methods to detect multi-SNP associations in the context of masking effects. We validate the method through simulations and applications to real data. We find that of 35 blood and urine biomarker traits in the UK Biobank, 34 traits show evidence for masking effects in a total of 4,273 gene-trait pairs, indicating that masking effects is common in complex traits. We further validate the improved power of our method in height, body mass index, and schizophrenia with different GWAS sample sizes and show that on average 95.7% of the genes detected only by mBAT-combo with smaller sample sizes can be identified by the single-SNP approach with a 1.7-fold increase in sample sizes. Eleven genes significant only in mBAT-combo for schizophrenia are confirmed by functionally informed fine-mapping or Mendelian randomization integrating gene expression data. The framework of mBAT-combo can be applied to any set of SNPs to refine trait-association signals hidden in genomic regions with complex LD structures.

Genome-Wide Integration of Genetic and Genomic Studies of Atopic Dermatitis: Insights into Genetic Architecture and Pathogenesis

Atopic dermatitis (AD) is a common heterogeneous, chronic, itching, and inflammatory skin disease. Genetic studies have identified multiple AD susceptibility genes. However, the genetic architecture of AD has not been elucidated. In this study, we conducted a large-scale meta-analysis of AD (35,647 cases and 1,013,885 controls) to characterize the genetic basis of AD. The heritability of AD in different datasets varied from 0.6 to 7.1%. We identified 31 previously unreported genes by integrating multiomics data. Among the 31 genes, MCL1 was identified as a potential treatment target for AD by mediating gene‒drug interactions. Tissue enrichment analyses and phenome-wide association study provided strong support for the role of the hemic and immune systems in AD. Across 1,207 complex traits and diseases, genetic correlations indicated that AD shared links with multiple respiratory phenotypes. The phenome-wide Mendelian randomization analysis (Mendelian randomization‒phenome-wide association study) revealed that the age of onset of diabetes exhibited a positive causal effect on AD (inverse-variance weighted β = 0.39, SEM = 0.09, P = 2.77 × 10^-5). Overall, these results provide important insights into the genetic architecture of AD and will lead to a more thorough and complete understanding of the molecular mechanisms underlying AD.

Genome Wide Association Study with Imputed Whole Genome Sequence Data Identifies a 431 kb Risk Haplotype on CFA18 for Congenital Laryngeal Paralysis in Alaskan Sled Dogs

Congenital laryngeal paralysis (CLP) is an inherited disorder that affects the ability of the dog to exercise and precludes it from functioning as a working sled dog. Though CLP is known to occur in Alaskan sled dogs (ASDs) since 1986, the genetic mutation underlying the disease has not been reported. Using a genome-wide association study (GWAS), we identified a 708 kb region on CFA 18 harboring 226 SNPs to be significantly associated with CLP. The significant SNPs explained 47.06% of the heritability of CLP. We narrowed the region to 431 kb through autozygosity mapping and found 18 of the 20 cases to be homozygous for the risk haplotype. Whole genome sequencing of two cases and a control ASD, and comparison with the genome of 657 dogs from various breeds, confirmed the homozygous status of the risk haplotype to be unique to the CLP cases. Most of the dogs that were homozygous for the risk allele had blue eyes. Gene annotation and a gene-based association study showed that the risk haplotype encompasses genes implicated in developmental and neurodegenerative disorders. Pathway analysis showed enrichment of glycoproteins and glycosaminoglycans biosynthesis, which play a key role in repairing damaged nerves. In conclusion, our results suggest an important role for the identified candidate region in CLP.

Cross-ancestry meta-analysis of opioid use disorder uncovers novel loci with predominant effects in brain regions associated with addiction

Despite an estimated heritability of ~50%, genome-wide association studies of opioid use disorder (OUD) have revealed few genome-wide significant loci. We conducted a cross-ancestry meta-analysis of OUD in the Million Veteran Program (N = 425,944). In addition to known exonic variants in OPRM1 and FURIN, we identified intronic variants in RABEPK, FBXW4, NCAM1 and KCNN1. A meta-analysis including other datasets identified a locus in TSNARE1. In total, we identified 14 loci for OUD, 12 of which are novel. Significant genetic correlations were identified for 127 traits, including psychiatric disorders and other substance use-related traits. The only significantly enriched cell-type group was CNS, with gene expression enrichment in brain regions previously associated with substance use disorders. These findings increase our understanding of the biological basis of OUD and provide further evidence that it is a brain disease, which may help to reduce stigma and inform efforts to address the opioid epidemic.

Single nucleotide polymorphism marker combinations for classifying Yeonsan Ogye chicken using a machine learning approach

Genetic analysis has great potential as a tool to differentiate between different species and breeds of livestock. In this study, the optimal combinations of single nucleotide polymorphism (SNP) markers for discriminating the Yeonsan Ogye chicken (Gallus gallus domesticus) breed were identified using high-density 600K SNP array data. In 3,904 individuals from 198 chicken breeds, SNP markers specific to the target population were discovered through a case-control genome-wide association study (GWAS) and filtered out based on the linkage disequilibrium blocks. Significant SNP markers were selected by feature selection applying two machine learning algorithms: Random Forest (RF) and AdaBoost (AB). Using a machine learning approach, the 38 (RF) and 43 (AB) optimal SNP marker combinations for the Yeonsan Ogye chicken population demonstrated 100% accuracy. Hence, the GWAS and machine learning models used in this study can be efficiently utilized to identify the optimal combination of markers for discriminating target populations using multiple SNP markers.

A comprehensive comparison of multilocus association methods with summary statistics in genome-wide association studies

BACKGROUND

Multilocus analysis on a set of single nucleotide polymorphisms (SNPs) pre-assigned within a gene constitutes a valuable complement to single-marker analysis by aggregating data on complex traits in a biologically meaningful way. However, despite the existence of a wide variety of SNP-set methods, few comprehensive comparison studies have been previously performed to evaluate the effectiveness of these methods.

RESULTS

We herein sought to fill this knowledge gap by conducting a comprehensive empirical comparison for 22 commonly-used summary-statistics based SNP-set methods. We showed that only seven methods could effectively control the type I error, and that these well-calibrated approaches had varying power performance under the simulation scenarios. Overall, we confirmed that the burden test was generally underpowered and score-based variance component tests (e.g., sequence kernel association test) were much powerful under the polygenic genetic architecture in both common and rare variant association analyses. We further revealed that two linkage-disequilibrium-free P value combination methods (e.g., harmonic mean P value method and aggregated Cauchy association test) behaved very well under the sparse genetic architecture in simulations and real-data applications to common and rare variant association analyses as well as in expression quantitative trait loci weighted integrative analysis. We also assessed the scalability of these approaches by recording computational time and found that all these methods can be scalable to biobank-scale data although some might be relatively slow.

CONCLUSION

In conclusion, we hope that our findings can offer an important guidance on how to choose appropriate multilocus association analysis methods in post-GWAS era. All the SNP-set methods are implemented in the R package called MCA, which is freely available at https://github.com/biostatpzeng/ .

Sodium-calcium exchanger-3 regulates pain "wind-up": From human psychophysics to spinal mechanisms

Repeated application of noxious stimuli leads to a progressively increased pain perception; this temporal summation is enhanced in and predictive of clinical pain disorders. Its electrophysiological correlate is "wind-up," in which dorsal horn spinal neurons increase their response to repeated nociceptor stimulation. To understand the genetic basis of temporal summation, we undertook a GWAS of wind-up in healthy human volunteers and found significant association with SLC8A3 encoding sodium-calcium exchanger type 3 (NCX3). NCX3 was expressed in mouse dorsal horn neurons, and mice lacking NCX3 showed normal, acute pain but hypersensitivity to the second phase of the formalin test and chronic constriction injury. Dorsal horn neurons lacking NCX3 showed increased intracellular calcium following repetitive stimulation, slowed calcium clearance, and increased wind-up. Moreover, virally mediated enhanced spinal expression of NCX3 reduced central sensitization. Our study highlights Ca2+ efflux as a pathway underlying temporal summation and persistent pain, which may be amenable to therapeutic targeting.

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.

Multi-Trait Genetic Analysis Identifies Autoimmune Loci Associated with Cutaneous Melanoma

Genome-wide association studies (GWAS) have identified a number of risk loci for cutaneous melanoma. Cutaneous melanoma shares overlapping genetic risk (genetic correlation) with a number of other traits, including its risk factors such as sunburn propensity. This genetic correlation can be exploited to identify additional cutaneous melanoma risk loci by multitrait analysis of GWAS (MTAG). We used bivariate linkage disequilibrium-score regression score regression to identify traits that are genetically correlated with clinically confirmed cutaneous melanoma and then used publicly available GWAS for these traits in a multitrait analysis of GWAS. Multitrait analysis of GWAS allows GWAS to be combined while accounting for sample overlap and incomplete genetic correlation. We identified a total of 74 genome-wide independent loci, 19 of them were not previously reported in the input cutaneous melanoma GWAS meta-analysis. Of these loci, 55 were replicated (P < 0.05/74, Bonferroni-corrected P-value in two independent cutaneous melanoma replication cohorts from Melanoma Institute Australia and 23andMe, Inc. Among the, to our knowledge, previously unreported cutaneous melanoma loci are ones that have also been associated with autoimmune traits including rs715199 near LPP and rs10858023 near AP4B1. Our analysis indicates genetic correlation between traits can be leveraged to identify new risk genes for cutaneous melanoma.

Genome-wide analysis of genetic predisposition to common polygenic cancers

Lung, breast, prostate, and colorectal cancers are among the most common and fatal malignancies worldwide. They are mainly caused by multifactorial mechanisms and are genetically heterogeneous. We investigated the genetic architecture of these cancers through genome-wide association, pathway-based, and summary-based transcriptome-/methylome-wide association analyses using three independent cohorts. Our genome-wide association analyses identified the associations of 33 single-nucleotide polymorphisms (SNPs) at P < 5E - 06, of which 32 SNPs were not previously reported and did not have proxy variants within their ± 1 Mb flanking regions. Moreover, other polymorphisms mapped to their closest genes were not previously associated with the same cancers at P < 5E - 06. Our pathway enrichment analyses revealed associations of 32 pathways; mainly related to the immune system, DNA replication/transcription, and chromosomal organization; with the studied cancers. Also, 60 probes were associated with these cancers in our transcriptome-wide and methylome-wide analyses. The ± 1 Mb flanking regions of most probes had not attained P < 5E - 06 in genome-wide association studies. The genes corresponding to the significant probes can be considered as potential targets for further functional studies. Two genes (i.e., CDC14A and PMEL) demonstrated stronger evidence of associations with lung cancer as they had significant probes in both transcriptome-wide and methylome-wide association analyses. The novel cancer-associated SNPs and genes identified here would advance our understanding of the genetic heterogeneity of the common cancers.

A multi-level investigation of the genetic relationship between endometriosis and ovarian cancer histotypes

Endometriosis is associated with increased risk of epithelial ovarian cancers (EOCs). Using data from large endometriosis and EOC genome-wide association meta-analyses, we estimate the genetic correlation and evaluate the causal relationship between genetic liability to endometriosis and EOC histotypes, and identify shared susceptibility loci. We estimate a significant genetic correlation (rg) between endometriosis and clear cell (rg = 0.71), endometrioid (rg = 0.48), and high-grade serous (rg = 0.19) ovarian cancer, associations supported by Mendelian randomization analyses. Bivariate meta-analysis identified 28 loci associated with both endometriosis and EOC, including 19 with evidence for a shared underlying association signal. Differences in the shared risk suggest different underlying pathways may contribute to the relationship between endometriosis and the different histotypes. Functional annotation using transcriptomic and epigenomic profiles of relevant tissues/cells highlights several target genes. This comprehensive analysis reveals profound genetic overlap between endometriosis and EOC histotypes with valuable genomic targets for understanding the biological mechanisms linking the diseases.

Distinct biological ages of organs and systems identified from a multi-omics study

Biological age (BA) has been proposed to evaluate the aging status instead of chronological age (CA). Our study shows evidence that there might be multiple "clocks" within the whole-body system: systemic aging drivers/clocks overlaid with organ/tissue-specific counterparts. We utilize multi-omics data, including clinical tests, immune repertoire, targeted metabolomic molecules, gut microbiomes, physical fitness examinations, and facial skin examinations, to estimate the BA of different organs (e.g., liver, kidney) and systems (immune and metabolic system). The aging rates of organs/systems are diverse. People's aging patterns are different. We also demonstrate several applications of organs/systems BA in two independent datasets. Mortality predictions are compared among organs' BA in the dataset of the United States National Health and Nutrition Examination Survey. Polygenic risk score of BAs constructed in the Chinese Longitudinal Healthy Longevity Survey cohort can predict the possibility of becoming centenarian.

A flexible summary statistics-based colocalization method with application to the mucin cystic fibrosis lung disease modifier locus

Mucus obstruction is a central feature in the cystic fibrosis (CF) airways. A genome-wide association study (GWAS) of lung disease by the CF Gene Modifier Consortium (CFGMC) identified a significant locus containing two mucin genes, MUC20 and MUC4. Expression quantitative trait locus (eQTL) analysis using human nasal epithelia (HNE) from 94 CF-affected Canadians in the CFGMC demonstrated MUC4 eQTLs that mirrored the lung association pattern in the region, suggesting that MUC4 expression may mediate CF lung disease. Complications arose, however, with colocalization testing using existing methods: the locus is complex and the associated SNPs span a 0.2 Mb region with high linkage disequilibrium (LD) and evidence of allelic heterogeneity. We previously developed the Simple Sum (SS), a powerful colocalization test in regions with allelic heterogeneity, but SS assumed eQTLs to be present to achieve type I error control. Here we propose a two-stage SS (SS2) colocalization test that avoids a priori eQTL assumptions, accounts for multiple hypothesis testing and the composite null hypothesis, and enables meta-analysis. We compare SS2 to published approaches through simulation and demonstrate type I error control for all settings with the greatest power in the presence of high LD and allelic heterogeneity. Applying SS2 to the MUC20/MUC4 CF lung disease locus with eQTLs from CF HNE revealed significant colocalization with MUC4 (p = 1.31 × 10-5) rather than with MUC20. The SS2 is a powerful method to inform the responsible gene(s) at a locus and guide future functional studies. SS2 has been implemented in the application LocusFocus.

STXBP6 and B3GNT6 Genes are Associated With Selective IgA Deficiency

Immunoglobulin A Deficiency (IgAD) is a polygenic primary immune deficiency, with a strong genetic association to the human leukocyte antigen (HLA) region. Previous genome-wide association studies (GWAS) have identified five non-HLA risk loci (IFIH1, PVT1, ATG13-AMBRA1, AHI1 and CLEC16A). In this study, we investigated the genetic interactions between different HLA susceptibility haplotypes and non-MHC genes in IgAD. To do this, we stratified IgAD subjects and healthy controls based on HLA haplotypes (N = 10,993), and then performed GWAS to identify novel genetic regions contributing to IgAD susceptibility. After replicating previously published HLA risk haplotypes, we compared individuals carrying at least one HLA risk allele (HLA-B*08:01-DRB1*03:01-DQB1*02:01 or HLA-DRB1*07:01-DQB1*02:02 or HLA-DRB1*01-DQB1*05:01) with individuals lacking an HLA risk allele. Subsequently, we stratified subjects based on the susceptibility alleles/haplotypes and performed gene-based association analysis using 572,856 SNPs and 24,125 genes. A significant genome-wide association in STXBP6 (rs4097492; p = 7.63 × 10^-9) was observed in the cohort carrying at least one MHC risk allele. We also identified a significant gene-based association for B3GNT6 (P _Gene = 2.1 × 10^-6) in patients not carrying known HLA susceptibility alleles. Our findings indicate that the etiology of IgAD differs depending on the genetic background of HLA susceptibility haplotypes.