Estimating SNP-Based Heritability and Genetic Correlation in Case-Control Studies Directly and with Summary Statistics

Protocol to estimate the heritability of drug response with GxEMM and identify gene-drug interactions with TxEWAS

Identifying factors that affect treatment response is a central objective of clinical research. Here, we present a protocol to study the genetic architecture of response to commonly prescribed drugs using gene-context interaction techniques. We describe steps for estimating the heritability of drug response with gene-environment interaction mixed model (GxEMM) and identifying gene-drug interactions with gene-environment interaction transcriptome-wide association study (TxEWAS). While the protocol describes application to drug treatments, this framework can be used to characterize the genetic basis of any covariate's effect. For complete details on the use and execution of this protocol, please refer to Sadowski et al.1.

Characterizing the genetic architecture of drug response using gene-context interaction methods

Identifying factors that affect treatment response is a central objective of clinical research, yet the role of common genetic variation remains largely unknown. Here, we develop a framework to study the genetic architecture of response to commonly prescribed drugs in large biobanks. We quantify treatment response heritability for statins, metformin, warfarin, and methotrexate in the UK Biobank. We find that genetic variation modifies the primary effect of statins on LDL cholesterol (9% heritable) as well as their side effects on hemoglobin A1c and blood glucose (10% and 11% heritable, respectively). We identify dozens of genes that modify drug response, which we replicate in a retrospective pharmacogenomic study. Finally, we find that polygenic score (PGS) accuracy varies up to 2-fold depending on treatment status, showing that standard PGSs are likely to underperform in clinical contexts.

GWAS-by-subtraction reveals an IOP-independent component of primary open angle glaucoma

The etiology of primary open angle glaucoma is constituted by both intraocular pressure-dependent and intraocular pressure-independent mechanisms. However, GWASs of traits affecting primary open angle glaucoma through mechanisms independent of intraocular pressure remains limited. Here, we address this gap by subtracting the genetic effects of a GWAS for intraocular pressure from a GWAS for primary open angle glaucoma to reveal the genetic contribution to primary open angle glaucoma via intraocular pressure-independent mechanisms. Seventeen independent genome-wide significant SNPs were associated with the intraocular pressure-independent component of primary open angle glaucoma. Of these, 7 are located outside known normal tension glaucoma loci, 11 are located outside known intraocular pressure loci, and 2 are novel primary open angle glaucoma loci. The intraocular pressure-independent genetic component of primary open angle glaucoma is associated with glaucoma endophenotypes, while the intraocular pressure-dependent component is associated with blood pressure and vascular permeability. A genetic risk score for the intraocular pressure-independent component of primary open angle glaucoma is associated with 26 different retinal micro-vascular features, which contrasts with the genetic risk score for the intraocular pressure-dependent component. Increased understanding of these intraocular pressure-dependent and intraocular pressure-independent components provides insights into the pathogenesis of glaucoma.

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.

Using Genetics to Investigate Relationships between Phenotypes: Application to Endometrial Cancer

Genome-wide association studies (GWAS) have accelerated the exploration of genotype-phenotype associations, facilitating the discovery of replicable genetic markers associated with specific traits or complex diseases. This narrative review explores the statistical methodologies developed using GWAS data to investigate relationships between various phenotypes, focusing on endometrial cancer, the most prevalent gynecological malignancy in developed nations. Advancements in analytical techniques such as genetic correlation, colocalization, cross-trait locus identification, and causal inference analyses have enabled deeper exploration of associations between different phenotypes, enhancing statistical power to uncover novel genetic risk regions. These analyses have unveiled shared genetic associations between endometrial cancer and many phenotypes, enabling identification of novel endometrial cancer risk loci and furthering our understanding of risk factors and biological processes underlying this disease. The current status of research in endometrial cancer is robust; however, this review demonstrates that further opportunities exist in statistical genetics that hold promise for advancing the understanding of endometrial cancer and other complex diseases.

A parametric bootstrap approach for computing confidence intervals for genetic correlations with application to genetically determined protein-protein networks

Genetic correlation refers to the correlation between genetic determinants of a pair of traits. When using individual-level data, it is typically estimated based on a bivariate model specification where the correlation between the two variables is identifiable and can be estimated from a covariance model that incorporates the genetic relationship between individuals, e.g., using a pre-specified kinship matrix. Inference relying on asymptotic normality of the genetic correlation parameter estimates may be inaccurate when the sample size is low, when the genetic correlation is close to the boundary of the parameter space, and when the heritability of at least one of the traits is low. We address this problem by developing a parametric bootstrap procedure to construct confidence intervals for genetic correlation estimates. The procedure simulates paired traits under a range of heritability and genetic correlation parameters, and it uses the population structure encapsulated by the kinship matrix. Heritabilities and genetic correlations are estimated using the close-form, method of moment, Haseman-Elston regression estimators. The proposed parametric bootstrap procedure is especially useful when genetic correlations are computed on pairs of thousands of traits measured on the same exact set of individuals. We demonstrate the parametric bootstrap approach on a proteomics dataset from the Jackson Heart Study.

A scalable and robust variance components method reveals insights into the architecture of gene-environment interactions underlying complex traits

Understanding the contribution of gene-environment interactions (GxE) to complex trait variation can provide insights into disease mechanisms, explain sources of heritability, and improve genetic risk prediction. While large biobanks with genetic and deep phenotypic data hold promise for obtaining novel insights into GxE, our understanding of GxE architecture in complex traits remains limited. We introduce a method to estimate the proportion of trait variance explained by GxE (GxE heritability) and additive genetic effects (additive heritability) across the genome and within specific genomic annotations. We show that our method is accurate in simulations and computationally efficient for biobank-scale datasets. We applied our method to common array SNPs (MAF ≥1%), fifty quantitative traits, and four environmental variables (smoking, sex, age, and statin usage) in unrelated white British individuals in the UK Biobank. We found 68 trait-E pairs with significant genome-wide GxE heritability (p<0.05/200) with a ratio of GxE to additive heritability of ≈6.8% on average. Analyzing ≈8 million imputed SNPs (MAF ≥0.1%), we documented an approximate 28% increase in genome-wide GxE heritability compared to array SNPs. We partitioned GxE heritability across minor allele frequency (MAF) and local linkage disequilibrium (LD) values, revealing that, like additive allelic effects, GxE allelic effects tend to increase with decreasing MAF and LD. Analyzing GxE heritability near genes highly expressed in specific tissues, we find significant brain-specific enrichment for body mass index (BMI) and basal metabolic rate in the context of smoking and adipose-specific enrichment for waist-hip ratio (WHR) in the context of sex.

Shared Genetic Architectures between Coronary Artery Disease and Type 2 Diabetes Mellitus in East Asian and European Populations

Coronary artery disease (CAD) is a common comorbidity of type 2 diabetes mellitus (T2DM). However, the pathophysiology connecting these two phenotypes remains to be further understood. Combined analysis in multi-ethnic populations can help contribute to deepening our understanding of biological mechanisms caused by shared genetic loci. We applied genetic correlation analysis and then performed conditional and joint association analyses in Chinese, Japanese, and European populations to identify the genetic variants jointly associated with CAD and T2DM. Next, the associations between genes and the two traits were also explored. Finally, fine-mapping and functional enrichment analysis were employed to identify the potential causal variants and pathways. Genetic correlation results indicated significant genetic overlap between CAD and T2DM in the three populations. Over 10,000 shared signals were identified, and 587 were shared by East Asian and European populations. Fifty-six novel shared genes were found to have significant effects on both CAD and T2DM. Most loci were fine-mapped to plausible causal variant sets. Several similarities and differences of the involved genes in GO terms and KEGG pathways were revealed across East Asian and European populations. These findings highlight the importance of immunoregulation, neuroregulation, heart development, and the regulation of glucose metabolism in shared etiological mechanisms between CAD and T2DM.

LDER-GE estimates phenotypic variance component of gene-environment interactions in human complex traits accurately with GE interaction summary statistics and full LD information

Gene-environment (GE) interactions are essential in understanding human complex traits. Identifying these interactions is necessary for deciphering the biological basis of such traits. In this study, we review state-of-art methods for estimating the proportion of phenotypic variance explained by genome-wide GE interactions and introduce a novel statistical method Linkage-Disequilibrium Eigenvalue Regression for Gene-Environment interactions (LDER-GE). LDER-GE improves the accuracy of estimating the phenotypic variance component explained by genome-wide GE interactions using large-scale biobank association summary statistics. LDER-GE leverages the complete Linkage Disequilibrium (LD) matrix, as opposed to only the diagonal squared LD matrix utilized by LDSC (Linkage Disequilibrium Score)-based methods. Our extensive simulation studies demonstrate that LDER-GE performs better than LDSC-based approaches by enhancing statistical efficiency by ~23%. This improvement is equivalent to a sample size increase of around 51%. Additionally, LDER-GE effectively controls type-I error rate and produces unbiased results. We conducted an analysis using UK Biobank data, comprising 307 259 unrelated European-Ancestry subjects and 966 766 variants, across 217 environmental covariate-phenotype (E-Y) pairs. LDER-GE identified 34 significant E-Y pairs while LDSC-based method only identified 23 significant E-Y pairs with 22 overlapped with LDER-GE. Furthermore, we employed LDER-GE to estimate the aggregated variance component attributed to multiple GE interactions, leading to an increase in the explained phenotypic variance with GE interactions compared to considering main genetic effects only. Our results suggest the importance of impacts of GE interactions on human complex traits.

Heritability and genetic correlations for sleep apnea, insomnia, and hypersomnia in a large clinical biobank

RATIONALE

Comorbid insomnia and sleep apnea is reported to have worse outcomes than either condition alone. The local genetic correlations of these disorders are unknown.

OBJECTIVES

To identify local genomic regions with heritability for clinically diagnosed sleep apnea and insomnia, and to identify local genetic correlations between these disorders and/or hypersomnia.

METHODS

Fifty thousand two hundred seventeen patients of European ancestry were examined. Global and local heritability and genetic correlations for independent regions were calculated, adjusting for obesity and other covariates.

RESULTS

Sleep apnea and insomnia were significantly globally heritable and had 118 and 168 genetic regions with local heritability p-values <.05, respectively. One region had a significant genetic correlation for sleep apnea and hypersomnia (p-value = 9.85 × 10-4).

CONCLUSIONS

Clinically diagnosed sleep apnea and insomnia have minimal shared genetic architecture, supporting genetically distinct comorbid insomnia and sleep apnea components. However, additional correlated regions may be identified with additional sample size and methodological improvements.

A parametric bootstrap approach for computing confidence intervals for genetic correlations with application to genetically-determined protein-protein networks

Genetic correlation refers to the correlation between genetic determinants of a pair of traits. When using individual-level data, it is typically estimated based on a bivariate model specification where the correlation between the two variables is identifiable and can be estimated from a covariance model that incorporates the genetic relationship between individuals, e.g., using a pre-specified kinship matrix. Inference relying on asymptotic normality of the genetic correlation parameter estimates may be inaccurate when the sample size is low, when the genetic correlation is close to the boundary of the parameter space, and when the heritability of at least one of the traits is low. We address this problem by developing a parametric bootstrap procedure to construct confidence intervals for genetic correlation estimates. The procedure simulates paired traits under a range of heritability and genetic correlation parameters, and it uses the population structure encapsulated by the kinship matrix. Heritabilities and genetic correlations are estimated using the close-form, method of moment, Haseman-Elston regression estimators. The proposed parametric bootstrap procedure is especially useful when genetic correlations are computed on pairs of thousands of traits measured on the same exact set of individuals. We demonstrate the parametric bootstrap approach on a proteomics dataset from the Jackson Heart Study.

Genetics of hair graying with age

Hair graying is an early and obvious phenotypic and physiological trait with age in humans. Several recent advances in molecular biology and genetics have increased our understanding of the mechanisms of hair graying, which elucidate genes related to the synthesis, transport, and distribution of melanin in hair follicles, as well as genes regulating these processes above. Therefore, we review these advances and examine the trends in the genetic aspects of hair graying from enrichment theory, Genome-Wide association studies, whole exome sequencing, gene expression studies, and animal models for hair graying with age, aiming to overview the changes in hair graying at the genetic level and establish the foundation for future research. Meanwhile, by summarizing the genetics, it's of great value to explore the possible mechanism, treatment, or even prevention of hair graying with age.

Step by step: towards a better understanding of the genetic architecture of Alzheimer's disease

Alzheimer's disease (AD) is considered to have a large genetic component. Our knowledge of this component has progressed over the last 10 years, thanks notably to the advent of genome-wide association studies and the establishment of large consortia that make it possible to analyze hundreds of thousands of cases and controls. The characterization of dozens of chromosomal regions associated with the risk of developing AD and (in some loci) the causal genes responsible for the observed disease signal has confirmed the involvement of major pathophysiological pathways (such as amyloid precursor protein metabolism) and opened up new perspectives (such as the central role of microglia and inflammation). Furthermore, large-scale sequencing projects are starting to reveal the major impact of rare variants - even in genes like APOE - on the AD risk. This increasingly comprehensive knowledge is now being disseminated through translational research; in particular, the development of genetic risk/polygenic risk scores is helping to identify the subpopulations more at risk or less at risk of developing AD. Although it is difficult to assess the efforts still needed to comprehensively characterize the genetic component of AD, several lines of research can be improved or initiated. Ultimately, genetics (in combination with other biomarkers) might help to redefine the boundaries and relationships between various neurodegenerative diseases.

The Musical Abilities, Pleiotropy, Language, and Environment (MAPLE) Framework for Understanding Musicality-Language Links Across the Lifespan

Using individual differences approaches, a growing body of literature finds positive associations between musicality and language-related abilities, complementing prior findings of links between musical training and language skills. Despite these associations, musicality has been often overlooked in mainstream models of individual differences in language acquisition and development. To better understand the biological basis of these individual differences, we propose the Musical Abilities, Pleiotropy, Language, and Environment (MAPLE) framework. This novel integrative framework posits that musical and language-related abilities likely share some common genetic architecture (i.e., genetic pleiotropy) in addition to some degree of overlapping neural endophenotypes, and genetic influences on musically and linguistically enriched environments. Drawing upon recent advances in genomic methodologies for unraveling pleiotropy, we outline testable predictions for future research on language development and how its underlying neurobiological substrates may be supported by genetic pleiotropy with musicality. In support of the MAPLE framework, we review and discuss findings from over seventy behavioral and neural studies, highlighting that musicality is robustly associated with individual differences in a range of speech-language skills required for communication and development. These include speech perception-in-noise, prosodic perception, morphosyntactic skills, phonological skills, reading skills, and aspects of second/foreign language learning. Overall, the current work provides a clear agenda and framework for studying musicality-language links using individual differences approaches, with an emphasis on leveraging advances in the genomics of complex musicality and language traits.

Dissecting Complex Traits Using Omics Data: A Review on the Linear Mixed Models and Their Application in GWAS

Genome-wide association study (GWAS) is the most popular approach to dissecting complex traits in plants, humans, and animals. Numerous methods and tools have been proposed to discover the causal variants for GWAS data analysis. Among them, linear mixed models (LMMs) are widely used statistical methods for regulating confounding factors, including population structure, resulting in increased computational proficiency and statistical power in GWAS studies. Recently more attention has been paid to pleiotropy, multi-trait, gene-gene interaction, gene-environment interaction, and multi-locus methods with the growing availability of large-scale GWAS data and relevant phenotype samples. In this review, we have demonstrated all possible LMMs-based methods available in the literature for GWAS. We briefly discuss the different LMM methods, software packages, and available open-source applications in GWAS. Then, we include the advantages and weaknesses of the LMMs in GWAS. Finally, we discuss the future perspective and conclusion. The present review paper would be helpful to the researchers for selecting appropriate LMM models and methods quickly for GWAS data analysis and would benefit the scientific society.

Cross-trait assortative mating is widespread and inflates genetic correlation estimates

The observation of genetic correlations between disparate human traits has been interpreted as evidence of widespread pleiotropy. Here, we introduce cross-trait assortative mating (xAM) as an alternative explanation. We observe that xAM affects many phenotypes and that phenotypic cross-mate correlation estimates are strongly associated with genetic correlation estimates (R²=74%). We demonstrate that existing xAM plausibly accounts for substantial fractions of genetic correlation estimates and that previously reported genetic correlation estimates between some pairs of psychiatric disorders are congruent with xAM alone. Finally, we provide evidence for a history of xAM at the genetic level using cross-trait even/odd chromosome polygenic score correlations. Together, our results demonstrate that previous reports have likely overestimated the true genetic similarity between many phenotypes.

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

Theory for liability-scale models of the underlying genetic basis of complex disease provides an important way to interpret, compare, and understand results generated from biological studies. In particular, through estimation of the liability-scale heritability (LSH), liability models facilitate an understanding and comparison of the relative importance of genetic and environmental risk factors that shape different clinically important disease outcomes. Increasingly, large-scale biobank studies that link genetic information to electronic health records, containing hundreds of disease diagnosis indicators that mostly occur infrequently within the sample, are becoming available. Here, we propose an extension of the existing liability-scale model theory suitable for estimating LSH in biobank studies of low-prevalence disease. In a simulation study, we find that our derived expression yields lower mean square error (MSE) and is less sensitive to prevalence misspecification as compared to previous transformations for diseases with ≤2% population prevalence and LSH of ≤0.45, especially if the biobank sample prevalence is less than that of the wider population. Applying our expression to 13 diagnostic outcomes of ≤3% prevalence in the UK Biobank study revealed important differences in LSH obtained from the different theoretical expressions that impact the conclusions made when comparing LSH across disease outcomes. This demonstrates the importance of careful consideration for estimation and prediction of low-prevalence disease outcomes and facilitates improved inference of the underlying genetic basis of ≤2% population prevalence diseases, especially where biobank sample ascertainment results in a healthier sample population.

Open problems in human trait genetics

Genetic studies of human traits have revolutionized our understanding of the variation between individuals, and yet, the genetics of most traits is still poorly understood. In this review, we highlight the major open problems that need to be solved, and by discussing these challenges provide a primer to the field. We cover general issues such as population structure, epistasis and gene-environment interactions, data-related issues such as ancestry diversity and rare genetic variants, and specific challenges related to heritability estimates, genetic association studies, and polygenic risk scores. We emphasize the interconnectedness of these problems and suggest promising avenues to address them.

A review of SNP heritability estimation methods

Over the past decade, statistical methods have been developed to estimate single nucleotide polymorphism (SNP) heritability, which measures the proportion of phenotypic variance explained by all measured SNPs in the data. Estimates of SNP heritability measure the degree to which the available genetic variants influence phenotypes and improve our understanding of the genetic architecture of complex phenotypes. In this article, we review the recently developed and commonly used SNP heritability estimation methods for continuous and binary phenotypes from the perspective of model assumptions and parameter optimization. We primarily focus on their capacity to handle multiple phenotypes and longitudinal measurements, their ability for SNP heritability partition and their use of individual-level data versus summary statistics. State-of-the-art statistical methods that are scalable to the UK Biobank dataset are also elucidated in detail.

SNP-based heritability and selection analyses: Improved models and new results

Complex-trait genetics has advanced dramatically through methods to estimate the heritability tagged by SNPs, both genome-wide and in genomic regions of interest such as those defined by functional annotations. The models underlying many of these analyses are inadequate, and consequently many SNP-heritability results published to date are inaccurate. Here, we review the modelling issues, both for analyses based on individual genotype data and association test statistics, highlighting the role of a low-dimensional model for the heritability of each SNP. We use state-of-art models to present updated results about how heritability is distributed with respect to functional annotations in the human genome, and how it varies with allele frequency, which can reflect purifying selection. Our results give finer detail to the picture that has emerged in recent years of complex trait heritability widely dispersed across the genome. Confounding due to population structure remains a problem that summary statistic analyses cannot reliably overcome. Also see the video abstract here: https://youtu.be/WC2u03V65MQ.

Predicting Physical Appearance from DNA Data-Towards Genomic Solutions

The idea of forensic DNA intelligence is to extract from genomic data any information that can help guide the investigation. The clues to the externally visible phenotype are of particular practical importance. The high heritability of the physical phenotype suggests that genetic data can be easily predicted, but this has only become possible with less polygenic traits. The forensic community has developed DNA-based predictive tools by employing a limited number of the most important markers analysed with targeted massive parallel sequencing. The complexity of the genetics of many other appearance phenotypes requires big data coupled with sophisticated machine learning methods to develop accurate genomic predictors. A significant challenge in developing universal genomic predictive methods will be the collection of sufficiently large data sets. These should be created using whole-genome sequencing technology to enable the identification of rare DNA variants implicated in phenotype determination. It is worth noting that the correctness of the forensic sketch generated from the DNA data depends on the inclusion of an age factor. This, however, can be predicted by analysing epigenetic data. An important limitation preventing whole-genome approaches from being commonly used in forensics is the slow progress in the development and implementation of high-throughput, low DNA input sequencing technologies. The example of palaeoanthropology suggests that such methods may possibly be developed in forensics.

Fast estimation of genetic correlation for biobank-scale data

Genetic correlation is an important parameter in efforts to understand the relationships among complex traits. Current methods that analyze individual genotype data for estimating genetic correlation are challenging to scale to large datasets. Methods that analyze summary data, while being computationally efficient, tend to yield estimates of genetic correlation with reduced precision. We propose SCORE (scalable genetic correlation estimator), a randomized method of moments estimator of genetic correlation that is both scalable and accurate. SCORE obtains more precise estimates of genetic correlations relative to summary-statistic methods that can be applied at scale; it achieves a 44% reduction in standard error relative to LD-score regression (LDSC) and a 20% reduction relative to high-definition likelihood (HDL) (averaged over all simulations). The efficiency of SCORE enables computation of genetic correlations on the UK Biobank dataset, consisting of ≈300 K individuals and ≈500 K SNPs, in a few h (orders of magnitude faster than methods that analyze individual data, such as GCTA). Across 780 pairs of traits in 291,273 unrelated white British individuals in the UK Biobank, SCORE identifies significant genetic correlation between 200 additional pairs of traits over LDSC (beyond the 245 pairs identified by both).

Early-Life Injuries and the Development of Attention-Deficit/Hyperactivity Disorder

Objective: To estimate phenotypic and familial association between early-life injuries and attention-deficit/hyperactivity disorder (ADHD) and the genetic contribution to the association using polygenic risk score for ADHD (PRS-ADHD) and genetic correlation analyses. Methods: Children born in Denmark between 1995-2010 (n = 786,543) were followed from age 5 years until a median age of 14 years (interquartile range: 10-18 years). Using ICD-10 diagnoses, we estimated hazard ratios (HRs) and absolute risks of ADHD by number of hospital/emergency ward-treated injuries by age 5. In a subset of ADHD cases and controls born 1995 to 2005 who had genetic data available (n = 16,580), we estimated incidence rate ratios (IRRs) for the association between PRS-ADHD and number of injuries before age 5 and the genetic correlation between ADHD and any injury before age 5. Results: Injuries were associated with ADHD (HR = 1.61; 95% CI, 1.55-1.66) in males (HR = 1.59; 1.53-1.65) and females (HR = 1.65; 1.54-1.77), with a dose-response relationship with number of injuries. The absolute ADHD risk by age 15 was 8.4% (3+ injuries) vs 3.1% (no injuries). ADHD was also associated with injuries in relatives, with a stronger association in first- than second-degree relatives. PRS-ADHD was marginally associated with the number of injuries in the general population (IRR = 1.06; 1.00-1.14), with a genetic correlation of 0.53 (0.21-0.85). Conclusions: Early-life injuries in individuals and their relatives were associated with a diagnosis of ADHD. However, even in children with the most injuries, more than 90% were not diagnosed with ADHD by age 15. Despite a low positive predictive value and that the impact of unmeasured factors such as parental behavior remains unclear, results indicate that the association is partly explained by genetics, suggesting that early-life injuries may represent or herald early behavioral manifestations of ADHD.

Exploring the genetic heterogeneity in major depression across diagnostic criteria

Major depressive disorder (MDD) is defined differently across genetic research studies and this may be a key source of heterogeneity. While previous literature highlights differences between minimal and strict phenotypes, the components contributing to this heterogeneity have not been identified. Using the cardinal symptoms (depressed mood/anhedonia) as a baseline, we build MDD phenotypes using five components-(1) five or more symptoms, (2) episode duration, (3) functional impairment, (4) episode persistence, and (5) episode recurrence-to determine the contributors to such heterogeneity. Thirty-two depression phenotypes which systematically incorporate different combinations of MDD components were created using the mental health questionnaire data within the UK Biobank. SNP-based heritabilities and genetic correlations with three previously defined major depression phenotypes were calculated (Psychiatric Genomics Consortium (PGC) defined depression, 23andMe self-reported depression and broad depression) and differences between estimates analysed. All phenotypes were heritable (h2SNP range: 0.102-0.162) and showed substantial genetic correlations with other major depression phenotypes (Rg range: 0.651-0.895 (PGC); 0.652-0.837 (23andMe); 0.699-0.900 (broad depression)). The strongest effect on SNP-based heritability was from the requirement for five or more symptoms (1.4% average increase) and for a long episode duration (2.7% average decrease). No significant differences were noted between genetic correlations. While there is some variation, the two cardinal symptoms largely reflect the genetic aetiology of phenotypes incorporating more MDD components. These components may index severity, however, their impact on heterogeneity in genetic results is likely to be limited.

Genetics of substance use disorders in the era of big data

Substance use disorders (SUDs) are conditions in which the use of legal or illegal substances, such as nicotine, alcohol or opioids, results in clinical and functional impairment. SUDs and, more generally, substance use are genetically complex traits that are enormously costly on an individual and societal basis. The past few years have seen remarkable progress in our understanding of the genetics, and therefore the biology, of substance use and abuse. Various studies - including of well-defined phenotypes in deeply phenotyped samples, as well as broadly defined phenotypes in meta-analysis and biobank samples - have revealed multiple risk loci for these common traits. A key emerging insight from this work establishes a biological and genetic distinction between quantity and/or frequency measures of substance use (which may involve low levels of use without dependence), versus symptoms related to physical dependence.

Utility of polygenic embryo screening for disease depends on the selection strategy

Polygenic risk scores (PRSs) have been offered since 2019 to screen in vitro fertilization embryos for genetic liability to adult diseases, despite a lack of comprehensive modeling of expected outcomes. Here we predict, based on the liability threshold model, the expected reduction in complex disease risk following polygenic embryo screening for a single disease. A strong determinant of the potential utility of such screening is the selection strategy, a factor that has not been previously studied. When only embryos with a very high PRS are excluded, the achieved risk reduction is minimal. In contrast, selecting the embryo with the lowest PRS can lead to substantial relative risk reductions, given a sufficient number of viable embryos. We systematically examine the impact of several factors on the utility of screening, including: variance explained by the PRS, number of embryos, disease prevalence, parental PRSs, and parental disease status. We consider both relative and absolute risk reductions, as well as population-averaged and per-couple risk reductions, and also examine the risk of pleiotropic effects. Finally, we confirm our theoretical predictions by simulating 'virtual' couples and offspring based on real genomes from schizophrenia and Crohn's disease case-control studies. We discuss the assumptions and limitations of our model, as well as the potential emerging ethical concerns.

Comparison of methods for estimating genetic correlation between complex traits using GWAS summary statistics

Genetic correlation is the correlation of phenotypic effects by genetic variants across the genome on two phenotypes. It is an informative metric to quantify the overall genetic similarity between complex traits, which provides insights into their polygenic genetic architecture. Several methods have been proposed to estimate genetic correlation based on data collected from genome-wide association studies (GWAS). Due to the easy access of GWAS summary statistics and computational efficiency, methods only requiring GWAS summary statistics as input have become more popular than methods utilizing individual-level genotype data. Here, we present a benchmark study for different summary-statistics-based genetic correlation estimation methods through simulation and real data applications. We focus on two major technical challenges in estimating genetic correlation: marker dependency caused by linkage disequilibrium (LD) and sample overlap between different studies. To assess the performance of different methods in the presence of these two challenges, we first conducted comprehensive simulations with diverse LD patterns and sample overlaps. Then we applied these methods to real GWAS summary statistics for a wide spectrum of complex traits. Based on these experiments, we conclude that methods relying on accurate LD estimation are less robust in real data applications due to the imprecision of LD obtained from reference panels. Our findings offer guidance on how to choose appropriate methods for genetic correlation estimation in post-GWAS analysis.

Genetic architecture of four smoking behaviors using partitioned SNP heritability

BACKGROUND AND AIMS

Although genome-wide association studies have identified many loci that influence smoking behaviors, much of the genetic variance remains unexplained. We characterized the genetic architecture of four smoking behaviors using single nucleotide polymorphism (SNP) heritability (h²_SNP ). This is an estimate of narrow-sense heritability specifically estimating the proportion of phenotypic variation due to causal variants (CVs) tagged by SNPs.

DESIGN

Partitioned h²_SNP analysis of smoking behavior traits.

SETTING

UK Biobank.

PARTICIPANTS

UK Biobank participants of European ancestry. The number of participants varied depending on the trait, from 54 792 to 323 068.

MEASUREMENTS

Smoking initiation, age of initiation, cigarettes per day (CPD; count, log-transformed, binned and dichotomized into heavy versus light) and smoking cessation with imputed genome-wide SNPs.

FINDINGS

We estimated that, in aggregate, approximately 18% of the phenotypic variance in smoking initiation was captured by imputed SNPs [h²_SNP = 0.18, standard error (SE) = 0.01] and 12% [SE = 0.02] for smoking cessation, both of which were more than twice the previously reported estimates. Estimated age of initiation (h²_SNP  = 0.05, SE = 0.01) and binned CPD (h²_SNP  = 0.1, SE = 0.01) were substantially below published twin-based h² of 50%. CPD encoding influenced estimates, with dichotomized CPD h²_SNP  = 0.28. There was no evidence of dominance genetic variance for any trait.

CONCLUSION

A biobank study of smoking behavior traits suggested that the phenotypic variance explained by SNPs of smoking initiation, age of initiation, cigarettes per day and smoking cessation is modest overall.

Analysis of genetic differences between psychiatric disorders: exploring pathways and cell types/tissues involved and ability to differentiate the disorders by polygenic scores

Although displaying genetic correlations, psychiatric disorders are clinically defined as categorical entities as they each have distinguishing clinical features and may involve different treatments. Identifying differential genetic variations between these disorders may reveal how the disorders differ biologically and help to guide more personalized treatment. Here we presented a statistical framework and comprehensive analysis to identify genetic markers differentially associated with various psychiatric disorders/traits based on GWAS summary statistics, covering 18 psychiatric traits/disorders and 26 comparisons. We also conducted comprehensive analysis to unravel the genes, pathways and SNP functional categories involved, and the cell types and tissues implicated. We also assessed how well one could distinguish between psychiatric disorders by polygenic risk scores (PRS). SNP-based heritabilities (h2snp) were significantly larger than zero for most comparisons. Based on current GWAS data, PRS have mostly modest power to distinguish between psychiatric disorders. For example, we estimated that AUC for distinguishing schizophrenia from major depressive disorder (MDD), bipolar disorder (BPD) from MDD and schizophrenia from BPD were 0.694, 0.602 and 0.618, respectively, while the maximum AUC (based on h2snp) were 0.763, 0.749 and 0.726, respectively. We also uncovered differences in each pair of studied traits in terms of their differences in genetic correlation with comorbid traits. For example, clinically defined MDD appeared to more strongly genetically correlated with other psychiatric disorders and heart disease, when compared to non-clinically defined depression in UK Biobank. Our findings highlight genetic differences between psychiatric disorders and the mechanisms involved. PRS may help differential diagnosis of selected psychiatric disorders in the future with larger GWAS samples.

Estimating heritability and its enrichment in tissue-specific gene sets in admixed populations

It is important to study the genetics of complex traits in diverse populations. Here, we introduce covariate-adjusted linkage disequilibrium (LD) score regression (cov-LDSC), a method to estimate SNP-heritability (${\boldsymbol{h}}_{\boldsymbol{g}}^{\mathbf{2}})$ and its enrichment in homogenous and admixed populations with summary statistics and in-sample LD estimates. In-sample LD can be estimated from a subset of the genome-wide association studies samples, allowing our method to be applied efficiently to very large cohorts. In simulations, we show that unadjusted LDSC underestimates ${\boldsymbol{h}}_{\boldsymbol{g}}^{\mathbf{2}}$ by 10-60% in admixed populations; in contrast, cov-LDSC is robustly accurate. We apply cov-LDSC to genotyping data from 8124 individuals, mostly of admixed ancestry, from the Slim Initiative in Genomic Medicine for the Americas study, and to approximately 161 000 Latino-ancestry individuals, 47 000 African American-ancestry individuals and 135 000 European-ancestry individuals, as classified by 23andMe. We estimate ${\boldsymbol{h}}_{\boldsymbol{g}}^{\mathbf{2}}$ and detect heritability enrichment in three quantitative and five dichotomous phenotypes, making this, to our knowledge, the most comprehensive heritability-based analysis of admixed individuals to date. Most traits have high concordance of ${\boldsymbol{h}}_{\boldsymbol{g}}^{\mathbf{2}}$ and consistent tissue-specific heritability enrichment among different populations. However, for age at menarche, we observe population-specific heritability estimates of ${\boldsymbol{h}}_{\boldsymbol{g}}^{\mathbf{2}}$. We observe consistent patterns of tissue-specific heritability enrichment across populations; for example, in the limbic system for BMI, the per-standardized-annotation effect size $ \tau $* is 0.16 ± 0.04, 0.28 ± 0.11 and 0.18 ± 0.03 in the Latino-, African American- and European-ancestry populations, respectively. Our approach is a powerful way to analyze genetic data for complex traits from admixed populations.

Quantifying the contribution of dominance deviation effects to complex trait variation in biobank-scale data

The proportion of variation in complex traits that can be attributed to non-additive genetic effects has been a topic of intense debate. The availability of biobank-scale datasets of genotype and trait data from unrelated individuals opens up the possibility of obtaining precise estimates of the contribution of non-additive genetic effects. We present an efficient method to estimate the variation in a complex trait that can be attributed to additive (additive heritability) and dominance deviation (dominance heritability) effects across all genotyped SNPs in a large collection of unrelated individuals. Over a wide range of genetic architectures, our method yields unbiased estimates of additive and dominance heritability. We applied our method, in turn, to array genotypes as well as imputed genotypes (at common SNPs with minor allele frequency [MAF] > 1%) and 50 quantitative traits measured in 291,273 unrelated white British individuals in the UK Biobank. Averaged across these 50 traits, we find that additive heritability on array SNPs is 21.86% while dominance heritability is 0.13% (about 0.48% of the additive heritability) with qualitatively similar results for imputed genotypes. We find no statistically significant evidence for dominance heritability (p<0.05/50 accounting for the number of traits tested) and estimate that dominance heritability is unlikely to exceed 1% for the traits analyzed. Our analyses indicate a limited contribution of dominance heritability to complex trait variation.

A unified framework for cross-population trait prediction by leveraging the genetic correlation of polygenic traits

The development of polygenic risk scores (PRSs) has proved useful to stratify the general European population into different risk groups. However, PRSs are less accurate in non-European populations due to genetic differences across different populations. To improve the prediction accuracy in non-European populations, we propose a cross-population analysis framework for PRS construction with both individual-level (XPA) and summary-level (XPASS) GWAS data. By leveraging trans-ancestry genetic correlation, our methods can borrow information from the Biobank-scale European population data to improve risk prediction in the non-European populations. Our framework can also incorporate population-specific effects to further improve construction of PRS. With innovations in data structure and algorithm design, our methods provide a substantial saving in computational time and memory usage. Through comprehensive simulation studies, we show that our framework provides accurate, efficient, and robust PRS construction across a range of genetic architectures. In a Chinese cohort, our methods achieved 7.3%-198.0% accuracy gain for height and 19.5%-313.3% accuracy gain for body mass index (BMI) in terms of predictive R2 compared to existing PRS approaches. We also show that XPA and XPASS can achieve substantial improvement for construction of height PRSs in the African population, suggesting the generality of our framework across global populations.

Boosting heritability: estimating the genetic component of phenotypic variation with multiple sample splitting

Background

Heritability is a central measure in genetics quantifying how much of the variability observed in a trait is attributable to genetic differences. Existing methods for estimating heritability are most often based on random-effect models, typically for computational reasons. The alternative of using a fixed-effect model has received much more limited attention in the literature.

Results

In this paper, we propose a generic strategy for heritability inference, termed as “boosting heritability”, by combining the advantageous features of different recent methods to produce an estimate of the heritability with a high-dimensional linear model. Boosting heritability uses in particular a multiple sample splitting strategy which leads in general to a stable and accurate estimate. We use both simulated data and real antibiotic resistance data from a major human pathogen, Sptreptococcus pneumoniae, to demonstrate the attractive features of our inference strategy.

Conclusions

Boosting is shown to offer a reliable and practically useful tool for inference about heritability.

Multiple measures of depression to enhance validity of major depressive disorder in the UK Biobank

BACKGROUND

The UK Biobank contains data with varying degrees of reliability and completeness for assessing depression. A third of participants completed a Mental Health Questionnaire (MHQ) containing the gold-standard Composite International Diagnostic Interview (CIDI) criteria for assessing mental health disorders.

AIMS

To investigate whether multiple observations of depression from sources other than the MHQ can enhance the validity of major depressive disorder (MDD).

METHOD

In participants who did not complete the MHQ, we calculated the number of other depression measures endorsed, for example from hospital episode statistics and interview data. We compared cases defined this way with CIDI-defined cases for several estimates: the variance explained by polygenic risk scores (PRS), area under the curve attributable to PRS, single nucleotide polymorphisms (SNPs)-based heritability and genetic correlations with summary statistics from the Psychiatric Genomics Consortium MDD genome-wide association study.

RESULTS

The strength of the genetic contribution increased with the number of measures endorsed. For example, SNP-based heritability increased from 7% in participants who endorsed only one measure of depression, to 21% in those who endorsed four or five measures of depression. The strength of the genetic contribution to cases defined by at least two measures approximated that for CIDI-defined cases. Most genetic correlations between UK Biobank and the Psychiatric Genomics Consortium MDD study exceeded 0.7, but there was variability between pairwise comparisons.

CONCLUSIONS

Multiple measures of depression can serve as a reliable approximation for case status where the CIDI measure is not available, indicating sample size can be optimised using the entire suite of UK Biobank data.

Accurate genetic and environmental covariance estimation with composite likelihood in genome-wide association studies

Genetic and environmental covariances between pairs of complex traits are important quantitative measurements that characterize their shared genetic and environmental architectures. Accurate estimation of genetic and environmental covariances in genome-wide association studies (GWASs) can help us identify common genetic and environmental factors associated with both traits and facilitate the investigation of their causal relationship. Genetic and environmental covariances are often modeled through multivariate linear mixed models. Existing algorithms for covariance estimation include the traditional restricted maximum likelihood (REML) method and the recent method of moments (MoM). Compared to REML, MoM approaches are computationally efficient and require only GWAS summary statistics. However, MoM approaches can be statistically inefficient, often yielding inaccurate covariance estimates. In addition, existing MoM approaches have so far focused on estimating genetic covariance and have largely ignored environmental covariance estimation. Here we introduce a new computational method, GECKO, for estimating both genetic and environmental covariances, that improves the estimation accuracy of MoM while keeping computation in check. GECKO is based on composite likelihood, relies on only summary statistics for scalable computation, provides accurate genetic and environmental covariance estimates across a range of scenarios, and can accommodate SNP annotation stratified covariance estimation. We illustrate the benefits of GECKO through simulations and applications on analyzing 22 traits from five large-scale GWASs. In the real data applications, GECKO identified 50 significant genetic covariances among analyzed trait pairs, resulting in a twofold power gain compared to the previous MoM method LDSC. In addition, GECKO identified 20 significant environmental covariances. The ability of GECKO to estimate environmental covariance in addition to genetic covariance helps us reveal strong positive correlation between the genetic and environmental covariance estimates across trait pairs, suggesting that common pathways may underlie the shared genetic and environmental architectures between traits.

Reviewing the genetics of heterogeneity in depression: operationalizations, manifestations and etiologies

With progress in genome-wide association studies of depression, from identifying zero hits in ~16 000 individuals in 2013 to 223 hits in more than a million individuals in 2020, understanding the genetic architecture of this debilitating condition no longer appears to be an impossible task. The pressing question now is whether recently discovered variants describe the etiology of a single disease entity. There are a myriad of ways to measure and operationalize depression severity, and major depressive disorder as defined in the Diagnostic and Statistical Manual of Mental Disorders-5 can manifest in more than 10 000 ways based on symptom profiles alone. Variations in developmental timing, comorbidity and environmental contexts across individuals and samples further add to the heterogeneity. With big data increasingly enabling genomic discovery in psychiatry, it is more timely than ever to explicitly disentangle genetic contributions to what is likely 'depressions' rather than depression. Here, we introduce three sources of heterogeneity: operationalization, manifestation and etiology. We review recent efforts to identify depression subtypes using clinical and data-driven approaches, examine differences in genetic architecture of depression across contexts, and argue that heterogeneity in operationalizations of depression is likely a considerable source of inconsistency. Finally, we offer recommendations and considerations for the field going forward.

Heritability jointly explained by host genotype and microbiome: will improve traits prediction?

As we observe the $70$th anniversary of the publication by Robertson that formalized the notion of 'heritability', geneticists remain puzzled by the problem of missing/hidden heritability, where heritability estimates from genome-wide association studies (GWASs) fall short of that from twin-based studies. Many possible explanations have been offered for this discrepancy, including existence of genetic variants poorly captured by existing arrays, dominance, epistasis and unaccounted-for environmental factors; albeit these remain controversial. We believe a substantial part of this problem could be solved or better understood by incorporating the host's microbiota information in the GWAS model for heritability estimation and may also increase human traits prediction for clinical utility. This is because, despite empirical observations such as (i) the intimate role of the microbiome in many complex human phenotypes, (ii) the overlap between genetic variants associated with both microbiome attributes and complex diseases and (iii) the existence of heritable bacterial taxa, current GWAS models for heritability estimate do not take into account the contributory role of the microbiome. Furthermore, heritability estimate from twin-based studies does not discern microbiome component of the observed total phenotypic variance. Here, we summarize the concept of heritability in GWAS and microbiome-wide association studies, focusing on its estimation, from a statistical genetics perspective. We then discuss a possible statistical method to incorporate the microbiome in the estimation of heritability in host GWAS.

Efficient variance components analysis across millions of genomes

While variance components analysis has emerged as a powerful tool in complex trait genetics, existing methods for fitting variance components do not scale well to large-scale datasets of genetic variation. Here, we present a method for variance components analysis that is accurate and efficient: capable of estimating one hundred variance components on a million individuals genotyped at a million SNPs in a few hours. We illustrate the utility of our method in estimating and partitioning variation in a trait explained by genotyped SNPs (SNP-heritability). Analyzing 22 traits with genotypes from 300,000 individuals across about 8 million common and low frequency SNPs, we observe that per-allele squared effect size increases with decreasing minor allele frequency (MAF) and linkage disequilibrium (LD) consistent with the action of negative selection. Partitioning heritability across 28 functional annotations, we observe enrichment of heritability in FANTOM5 enhancers in asthma, eczema, thyroid and autoimmune disorders.

Exploring the possibility of predicting human head hair greying from DNA using whole-exome and targeted NGS data

BACKGROUND

Greying of the hair is an obvious sign of human aging. In addition to age, sex- and ancestry-specific patterns of hair greying are also observed and the progression of greying may be affected by environmental factors. However, little is known about the genetic control of this process. This study aimed to assess the potential of genetic data to predict hair greying in a population of nearly 1000 individuals from Poland.

RESULTS

The study involved whole-exome sequencing followed by targeted analysis of 378 exome-wide and literature-based selected SNPs. For the selection of predictors, the minimum redundancy maximum relevance (mRMRe) method was used, and then two prediction models were developed. The models included age, sex and 13 unique SNPs. Two SNPs of the highest mRMRe score included whole-exome identified KIF1A rs59733750 and previously linked with hair loss FGF5 rs7680591. The model for greying vs. no greying prediction achieved accuracy of cross-validated AUC = 0.873. In the 3-grade classification cross-validated AUC equalled 0.864 for no greying, 0.791 for mild greying and 0.875 for severe greying. Although these values present fairly accurate prediction, most of the prediction information was brought by age alone. Genetic variants explained < 10% of hair greying variation and the impact of particular SNPs on prediction accuracy was found to be small.

CONCLUSIONS

The rate of changes in human progressive traits shows inter-individual variation, therefore they are perceived as biomarkers of the biological age of the organism. The knowledge on the mechanisms underlying phenotypic aging can be of special interest to the medicine, cosmetics industry and forensics. Our study improves the knowledge on the genetics underlying hair greying processes, presents prototype models for prediction and proves hair greying being genetically a very complex trait. Finally, we propose a four-step approach based on genetic and epigenetic data analysis allowing for i) sex determination; ii) genetic ancestry inference; iii) greying-associated SNPs assignment and iv) epigenetic age estimation, all needed for a final prediction of greying.

Transethnic Genome-Wide Association Study Provides Insights in the Genetic Architecture and Heritability of Long QT Syndrome

BACKGROUND

Long QT syndrome (LQTS) is a rare genetic disorder and a major preventable cause of sudden cardiac death in the young. A causal rare genetic variant with large effect size is identified in up to 80% of probands (genotype positive) and cascade family screening shows incomplete penetrance of genetic variants. Furthermore, a proportion of cases meeting diagnostic criteria for LQTS remain genetically elusive despite genetic testing of established genes (genotype negative). These observations raise the possibility that common genetic variants with small effect size contribute to the clinical picture of LQTS. This study aimed to characterize and quantify the contribution of common genetic variation to LQTS disease susceptibility.

METHODS

We conducted genome-wide association studies followed by transethnic meta-analysis in 1656 unrelated patients with LQTS of European or Japanese ancestry and 9890 controls to identify susceptibility single nucleotide polymorphisms. We estimated the common variant heritability of LQTS and tested the genetic correlation between LQTS susceptibility and other cardiac traits. Furthermore, we tested the aggregate effect of the 68 single nucleotide polymorphisms previously associated with the QT-interval in the general population using a polygenic risk score.

RESULTS

Genome-wide association analysis identified 3 loci associated with LQTS at genome-wide statistical significance (P<5×10-8) near NOS1AP, KCNQ1, and KLF12, and 1 missense variant in KCNE1(p.Asp85Asn) at the suggestive threshold (P<10-6). Heritability analyses showed that ≈15% of variance in overall LQTS susceptibility was attributable to common genetic variation (h2SNP 0.148; standard error 0.019). LQTS susceptibility showed a strong genome-wide genetic correlation with the QT-interval in the general population (rg=0.40; P=3.2×10-3). The polygenic risk score comprising common variants previously associated with the QT-interval in the general population was greater in LQTS cases compared with controls (P<10-13), and it is notable that, among patients with LQTS, this polygenic risk score was greater in patients who were genotype negative compared with those who were genotype positive (P<0.005).

CONCLUSIONS

This work establishes an important role for common genetic variation in susceptibility to LQTS. We demonstrate overlap between genetic control of the QT-interval in the general population and genetic factors contributing to LQTS susceptibility. Using polygenic risk score analyses aggregating common genetic variants that modulate the QT-interval in the general population, we provide evidence for a polygenic architecture in genotype negative LQTS.

Statistical Methods in Genome-Wide Association Studies

Since the initial success of genome-wide association studies (GWAS) in 2005, tens of thousands of genetic variants have been identified for hundreds of human diseases and traits. In a GWAS, genotype information at up to millions of genetic markers is collected from up to hundreds of thousands of individuals, together with their phenotype information. Several scientific goals can be accomplished through the analysis of GWAS data, including the identification of variants, genes, and pathways associated with diseases and traits of interest; the inference of the genetic architecture of these traits; and the development of genetic risk prediction models. In this review, we provide an overview of the statistical challenges in achieving these goals and recent progress in statistical methodology to address these challenges.

The Genetics of the Mood Disorder Spectrum: Genome-wide Association Analyses of More Than 185,000 Cases and 439,000 Controls

BACKGROUND

Mood disorders (including major depressive disorder and bipolar disorder) affect 10% to 20% of the population. They range from brief, mild episodes to severe, incapacitating conditions that markedly impact lives. Multiple approaches have shown considerable sharing of risk factors across mood disorders despite their diagnostic distinction.

METHODS

To clarify the shared molecular genetic basis of major depressive disorder and bipolar disorder and to highlight disorder-specific associations, we meta-analyzed data from the latest Psychiatric Genomics Consortium genome-wide association studies of major depression (including data from 23andMe) and bipolar disorder, and an additional major depressive disorder cohort from UK Biobank (total: 185,285 cases, 439,741 controls; nonoverlapping N = 609,424).

RESULTS

Seventy-three loci reached genome-wide significance in the meta-analysis, including 15 that are novel for mood disorders. More loci from the Psychiatric Genomics Consortium analysis of major depression than from that for bipolar disorder reached genome-wide significance. Genetic correlations revealed that type 2 bipolar disorder correlates strongly with recurrent and single-episode major depressive disorder. Systems biology analyses highlight both similarities and differences between the mood disorders, particularly in the mouse brain cell types implicated by the expression patterns of associated genes. The mood disorders also differ in their genetic correlation with educational attainment-the relationship is positive in bipolar disorder but negative in major depressive disorder.

CONCLUSIONS

The mood disorders share several genetic associations, and genetic studies of major depressive disorder and bipolar disorder can be combined effectively to enable the discovery of variants not identified by studying either disorder alone. However, we demonstrate several differences between these disorders. Analyzing subtypes of major depressive disorder and bipolar disorder provides evidence for a genetic mood disorders spectrum.

Quantitative genome-wide association study of six phenotypic subdomains identifies novel genome-wide significant variants in autism spectrum disorder

Autism spectrum disorders (ASD) are highly heritable and are characterized by deficits in social communication and restricted and repetitive behaviors. Twin studies on phenotypic subdomains suggest a differing underlying genetic etiology. Studying genetic variation explaining phenotypic variance will help to identify specific underlying pathomechanisms. We investigated the effect of common variation on ASD subdomains in two cohorts including >2500 individuals. Based on the Autism Diagnostic Interview-Revised (ADI-R), we identified and confirmed six subdomains with a SNP-based genetic heritability h²_SNP = 0.2-0.4. The subdomains nonverbal communication (NVC), social interaction (SI), and peer interaction (PI) shared genetic risk factors, while the subdomains of repetitive sensory-motor behavior (RB) and restricted interests (RI) were genetically independent of each other. The polygenic risk score (PRS) for ASD as categorical diagnosis explained 2.3-3.3% of the variance of SI, joint attention (JA), and PI, 4.5% for RI, 1.2% of RB, but only 0.7% of NVC. We report eight genome-wide significant hits-partially replicating previous findings-and 292 known and novel candidate genes. The underlying biological mechanisms were related to neuronal transmission and development. At the SNP and gene level, all subdomains showed overlap, with the exception of RB. However, no overlap was observed at the functional level. In summary, the ADI-R algorithm-derived subdomains related to social communication show a shared genetic etiology in contrast to restricted and repetitive behaviors. The ASD-specific PRS overlapped only partially, suggesting an additional role of specific common variation in shaping the phenotypic expression of ASD subdomains.

Genome-wide gene-environment analyses of major depressive disorder and reported lifetime traumatic experiences in UK Biobank

Depression is more frequent among individuals exposed to traumatic events. Both trauma exposure and depression are heritable. However, the relationship between these traits, including the role of genetic risk factors, is complex and poorly understood. When modelling trauma exposure as an environmental influence on depression, both gene-environment correlations and gene-environment interactions have been observed. The UK Biobank concurrently assessed Major Depressive Disorder (MDD) and self-reported lifetime exposure to traumatic events in 126,522 genotyped individuals of European ancestry. We contrasted genetic influences on MDD stratified by reported trauma exposure (final sample size range: 24,094-92,957). The SNP-based heritability of MDD with reported trauma exposure (24%) was greater than MDD without reported trauma exposure (12%). Simulations showed that this is not confounded by the strong, positive genetic correlation observed between MDD and reported trauma exposure. We also observed that the genetic correlation between MDD and waist circumference was only significant in individuals reporting trauma exposure (rg = 0.24, p = 1.8 × 10-7 versus rg = -0.05, p = 0.39 in individuals not reporting trauma exposure, difference p = 2.3 × 10-4). Our results suggest that the genetic contribution to MDD is greater when reported trauma is present, and that a complex relationship exists between reported trauma exposure, body composition, and MDD.

Dissecting the heritable risk of breast cancer: From statistical methods to susceptibility genes

Decades of research have shown that rare highly penetrant mutations can promote tumorigenesis, but it is still unclear whether variants observed at high-frequency in the broader population could modulate the risk of developing cancer. Genome-wide Association Studies (GWAS) have generated a wealth of data linking single nucleotide polymorphisms (SNPs) to increased cancer risk, but the effect of these mutations are usually subtle, leaving most of cancer heritability unexplained. Understanding the role of high-frequency mutations in cancer can provide new intervention points for early diagnostics, patient stratification and treatment in malignancies with high prevalence, such as breast cancer. Here we review state-of-the-art methods to study cancer heritability using GWAS data and provide an updated map of breast cancer susceptibility loci at the SNP and gene level.

Statistical methods for SNP heritability estimation and partition: A review

In GWAS studies, SNP heritability measures the proportion of phenotypic variance explained by all measured SNPs. Accurate estimation of SNP heritability can help us better understand the degree to which measured genetic variants influence phenotypes. Over the last decade, a variety of statistical methods and software tools have been developed for SNP heritability estimation with different data types including genotype array data, imputed genotype data, whole-genome sequencing data, RNA sequencing data, and bisulfite sequencing data. However, a thorough technical review of these methods, especially from a statistical and computational viewpoint, is currently missing. To fill this knowledge gap, we present a comprehensive review on a broad category of recently developed and commonly used SNP heritability estimation methods. We focus on their modeling assumptions; their interconnected relationships; their applicability to quantitative, binary and count phenotypes; their use of individual level data versus summary statistics, as well as their utility for SNP heritability partitioning. We hope that this review will serve as a useful reference for both methodologists who develop heritability estimation methods and practitioners who perform heritability analysis.

Estimating narrow-sense heritability using family data from admixed populations

Estimating total narrow-sense heritability in admixed populations remains an open question. In this work, we used extensive simulations to evaluate existing linear mixed-model frameworks for estimating total narrow-sense heritability in two population-based cohorts from Greenland, and compared the results with data from unadmixed individuals from Denmark. When our analysis focused on Greenlandic sib pairs, and under the assumption that shared environment among siblings has a negligible effect, the model with two relationship matrices, one capturing identity by descent and one capturing identity by state, returned heritability estimates close to the true simulated value, while using each of the two matrices alone led to downward biases. When phenotypes correlated with ancestry, heritability estimates were inflated. Based on these observations, we propose a PCA-based adjustment that recovers the true simulated heritability. We use this knowledge to estimate the heritability of ten quantitative traits from the two Greenlandic cohorts, and report differences such as lower heritability for height in Greenlanders compared with Europeans. In conclusion, narrow-sense heritability in admixed populations is best estimated when using a mixture of genetic relationship matrices on individuals with at least one first-degree relative included in the sample.

Minimal phenotyping yields genome-wide association signals of low specificity for major depression

Minimal phenotyping refers to the reliance on the use of a small number of self-reported items for disease case identification, increasingly used in genome-wide association studies (GWAS). Here we report differences in genetic architecture between depression defined by minimal phenotyping and strictly defined major depressive disorder (MDD): the former has a lower genotype-derived heritability that cannot be explained by inclusion of milder cases and a higher proportion of the genome contributing to this shared genetic liability with other conditions than for strictly defined MDD. GWAS based on minimal phenotyping definitions preferentially identifies loci that are not specific to MDD, and, although it generates highly predictive polygenic risk scores, the predictive power can be explained entirely by large sample sizes rather than by specificity for MDD. Our results show that reliance on results from minimal phenotyping may bias views of the genetic architecture of MDD and impede the ability to identify pathways specific to MDD.