Liability-scale heritability estimation for biobank studies of low-prevalence disease

Genetic association studies using disease liabilities from deep neural networks

The case-control study is a widely used method for investigating the genetic underpinnings of binary traits. However, long-term, prospective cohort studies often grapple with absent or evolving health-related outcomes. Here, we propose two methods, liability and meta, for conducting genome-wide association studies (GWASs) that leverage disease liabilities calculated from deep patient phenotyping. Analyzing 38 common traits in ∼300,000 UK Biobank participants, we identified an increased number of loci in comparison to the number identified by the conventional case-control approach, and there were high replication rates in larger external GWASs. Further analyses confirmed the disease specificity of the genetic architecture; the meta method demonstrated higher robustness when phenotypes were imputed with low accuracy. Additionally, polygenic risk scores based on disease liabilities more effectively predicted newly diagnosed cases in the 2022 dataset, which were controls in the earlier 2019 dataset. Our findings demonstrate that integrating high-dimensional phenotypic data into deep neural networks enhances genetic association studies while capturing disease-relevant genetic architecture.

MUTATE: a human genetic atlas of multiorgan artificial intelligence endophenotypes using genome-wide association summary statistics

Artificial intelligence (AI) has been increasingly integrated into imaging genetics to provide intermediate phenotypes (i.e. endophenotypes) that bridge the genetics and clinical manifestations of human disease. However, the genetic architecture of these AI endophenotypes remains largely unexplored in the context of human multiorgan system diseases. Using publicly available genome-wide association study summary statistics from the UK Biobank (UKBB), FinnGen, and the Psychiatric Genomics Consortium, we comprehensively depicted the genetic architecture of 2024 multiorgan AI endophenotypes (MAEs). We comparatively assessed the single-nucleotide polymorphism-based heritability, polygenicity, and natural selection signatures of 2024 MAEs using methods commonly used in the field. Genetic correlation and Mendelian randomization analyses reveal both within-organ relationships and cross-organ interconnections. Bi-directional causal relationships were established between chronic human diseases and MAEs across multiple organ systems, including Alzheimer's disease for the brain, diabetes for the metabolic system, asthma for the pulmonary system, and hypertension for the cardiovascular system. Finally, we derived polygenic risk scores for the 2024 MAEs for individuals not used to calculate MAEs and returned these to the UKBB. Our findings underscore the promise of the MAEs as new instruments to ameliorate overall human health. All results are encapsulated into the MUlTiorgan AI endophenoTypE genetic atlas and are publicly available at https://labs-laboratory.com/mutate.

Genome-wide association analysis provides insights into the molecular etiology of dilated cardiomyopathy

Dilated cardiomyopathy (DCM) is a leading cause of heart failure and cardiac transplantation. We report a genome-wide association study and multi-trait analysis of DCM (14,256 cases) and three left ventricular traits (36,203 UK Biobank participants). We identified 80 genomic risk loci and prioritized 62 putative effector genes, including several with rare variant DCM associations (MAP3K7, NEDD4L and SSPN). Using single-nucleus transcriptomics, we identify cellular states, biological pathways, and intracellular communications that drive pathogenesis. We demonstrate that polygenic scores predict DCM in the general population and modify penetrance in carriers of rare DCM variants. Our findings may inform the design of genetic testing strategies that incorporate polygenic background. They also provide insights into the molecular etiology of DCM that may facilitate the development of targeted therapeutics.

Converting estimated breeding values from the observed to probability scale for health traits

Dairy cattle health traits are paramount from a welfare and economic viewpoint, and modern breeding programs therefore prioritize the genetic improvement of these traits. Estimated breeding values for health traits are published as the probability of animals staying healthy. They are obtained using threshold models, which assume that the observed binary phenotype (i.e., healthy or sick) is dictated by an underlying normally distributed liability exceeding or not exceeding a threshold. This methodology requires significant computing time and faces convergence challenges, as it implies a nonlinear system of equations. Linear models have more straightforward computations and provide a robust approximation to threshold models; thus, they could be used to overcome these challenges. However, linear models yield estimated breeding values on the observed scale, requiring an approximation to the liability scale analogous to that from threshold models to later obtain the estimated breeding values on the probability scale. In addition, the robustness of the approximation of linear to threshold models depends on the amount of information and the incidence of the trait, with extreme incidence (i.e., ≤5%) deviating from optimal approximation. Our objective was to test a transformation from the observed to the liability, and then to the probability scale, in the genetic evaluation of health traits with moderate and very low (extreme) incidence. Data comprised displaced abomasum (5.1 million), ketosis (3.6 million), lameness (5 million), and mastitis (6.3 million) records from a Holstein population with a pedigree of 6 million animals, of which 1.7 million were genotyped. Univariate threshold and linear models were performed to predict breeding values. The agreement between estimated breeding values on the probability scale derived from threshold and linear models was assessed using Spearman rank correlations and comparison of estimated breeding values distributions. Correlations were at least 0.95, and estimated breeding value distributions almost entirely overlapped for all the traits but displaced abomasum, the trait with the lowest incidence (2%). Computing time was ∼3 times longer for threshold than for linear models. In this Holstein population, the approximation was suboptimal for a trait with extreme incidence (2%). However, when the incidence was ≥6%, the approximation was robust, and its use is recommended along with linear models for analyzing categorical traits in large populations to ease the computational burden.

No evidence that ACE2 or TMPRSS2 drive population disparity in COVID risks

Early in the SARS-CoV2 pandemic, in this journal, Hou et al. (BMC Med 18:216, 2020) interpreted public genotype data, run through functional prediction tools, as suggesting that members of particular human populations carry potentially COVID-risk-increasing variants in genes ACE2 and TMPRSS2 far more often than do members of other populations. Beyond resting on predictions rather than clinical outcomes, and focusing on variants too rare to typify population members even jointly, their claim mistook a well known artifact (that large samples reveal more of a population's variants than do small samples) as if showing real and congruent population differences for the two genes, rather than lopsided population sampling in their shared source data. We explain that artifact, and contrast it with empirical findings, now ample, that other loci shape personal COVID risks far more significantly than do ACE2 and TMPRSS2-and that variation in ACE2 and TMPRSS2 per se unlikely exacerbates any net population disparity in the effects of such more risk-informative loci.

Estimating disease heritability from complex pedigrees allowing for ascertainment and covariates

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

Genome-wide association study meta-analysis supports association between MUC1 and ectopic pregnancy

STUDY QUESTION

Can we identify genetic variants associated with ectopic pregnancy by undertaking the first genome-wide association study (GWAS) leveraging two large-scale biobank initiatives?

SUMMARY ANSWER

We identified two novel genome-wide significant associations with ectopic pregnancy, highlighting MUC1 (mucin 1) as the most plausible affected gene.

WHAT IS KNOWN ALREADY

Ectopic pregnancy is an important cause of maternal morbidity and mortality worldwide. Despite being a common early pregnancy complication, the genetic predisposition to this condition remains understudied and no large scale genetic studies have been performed so far.

STUDY DESIGN, SIZE, DURATION

A GWAS meta-analysis including 7070 women with ectopic pregnancy and 248 810 controls from Estonian Biobank and the FinnGen study.

PARTICIPANTS/MATERIALS, SETTING, METHODS

We identified ectopic pregnancy cases from national registers by ICD (International Classification of Disease) codes (ICD-10 O00), and all remaining women were considered controls. We carried out standard GWAS meta-analysis and additionally annotated GWAS signals, analysed co-localization with quantitative trait loci, estimated genetic correlations and identified associated phenotypes to characterize the genetic signals, as well as to analyse the genetic and phenotypic relationships with the condition.

MAIN RESULTS AND THE ROLE OF CHANCE

We identified two genome-wide significant loci on chromosomes 1 (rs4971091, P = 5.32×10-9) and 10 (rs11598956, P = 2.41×10-8) potentially associated with ectopic pregnancy. Follow-up analyses propose MUC1, which codes for an epithelial glycoprotein with an important role in barrier function, as the most likely candidate gene for the association on chromosome 1. We also characterize the phenotypic and genetic correlations with other phenotypes, identifying a genetic correlation with smoking and diseases of the (genito)urinary and gastrointestinal system, and phenotypic correlations with various reproductive health diagnoses, reflecting the previously known epidemiological associations.

LARGE SCALE DATA

The GWAS meta-analysis summary statistics are available from the GWAS Catalogue (GCST90272883).

LIMITATIONS, REASONS FOR CAUTION

The main limitation is that the findings are based on European-based ancestry populations, with limited data on other populations, and we only captured maternal genomes. Additionally, further larger meta-analysis or independent studies are needed to validate these findings.

WIDER IMPLICATIONS OF THE FINDINGS

This study encourages the use of large-scale genetic datasets to unravel genetic factors linked to ectopic pregnancy, which is difficult to study in experimental settings. Increased sample size might bring additional genetic factors associating with ectopic pregnancy and inform its heritability. Altogether, our results provide more insight into the biology of ectopic pregnancy and, accordingly, the biological processes governing embryo implantation.

STUDY FUNDING/COMPETING INTEREST(S)

N.P.G. was supported by MATER Marie Sklodowska-Curie which received funding from the European Union's Horizon 2020 research and innovation program under grant agreement No. 813707. This study was funded by European Union through the European Regional Development Fund Project No. 2014-2020.4.01.15-0012 GENTRANSMED. Computations were performed in the High-Performance Computing Center of University of Tartu. The authors declare no competing interests.

Accurate and efficient estimation of local heritability using summary statistics and the linkage disequilibrium matrix

Existing SNP-heritability estimators that leverage summary statistics from genome-wide association studies (GWAS) are much less efficient (i.e., have larger standard errors) than the restricted maximum likelihood (REML) estimators which require access to individual-level data. We introduce a new method for local heritability estimation-Heritability Estimation with high Efficiency using LD and association Summary Statistics (HEELS)-that significantly improves the statistical efficiency of summary-statistics-based heritability estimator and attains comparable statistical efficiency as REML (with a relative statistical efficiency >92%). Moreover, we propose representing the empirical LD matrix as the sum of a low-rank matrix and a banded matrix. We show that this way of modeling the LD can not only reduce the storage and memory cost, but also improve the computational efficiency of heritability estimation. We demonstrate the statistical efficiency of HEELS and the advantages of our proposed LD approximation strategies both in simulations and through empirical analyses of the UK Biobank data.

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.