Reverse GWAS: Using genetics to identify and model phenotypic subtypes

Red Blood Cell-Related Phenotype-Genotype Correlations in Chronic and Acute Critical Illnesses (Traumatic Brain Injury Cohort and COVID-19 Cohort)

Changes in red blood cell (RBC)-related parameters and anemia are common in both severe chronic and acute diseases. RBC-related phenotypes have a heritable component. However, it is unclear whether the contribution of genetic variability is pronounced when hematological parameters are affected by physiological stress. In this study, we analyzed RBC-related phenotypes and phenotype-genotype correlations in two exome-sequenced patient cohorts with or at a high risk for a critical illness: chronic TBI patients admitted for rehabilitation and patients with acute COVID-19. In the analysis of exome data, we focused on the cumulative effects of rare high-impact variants (qualifying variants, QVs) in specific gene sets, represented by Notch signaling pathway genes, based on the results of enrichment analysis in anemic TBI patients and three predefined gene sets for phenotypes of interest derived from GO, GWAS, and HPO resources. In both patient cohorts, anemia was associated with the cumulative effects of QVs in the GO (TBI: p = 0.0003, OR = 2.47 (1.54-4.88); COVID-19: p = 0.0004, OR = 2.12 (1.39-3.25)) and Notch pathway-derived (TBI: p = 0.0017, OR = 2.33 (1.35-4.02); COVID-19: p = 0.0012, OR =8.00 (1.79-35.74)) gene sets. In the multiple linear regression analysis, genetic variables contributed to RBC indices in patients with TBI. In COVID-19 patients, QVs in Notch pathway genes influenced RBC, HGB, and HCT levels, whereas genes from other sets influenced MCHC levels. Thus, in this exploratory study, exome data analysis yielded similar and different results in the two patient cohorts, supporting the view that genetic factors may contribute to RBC-related phenotypic performance in both severe chronic and acute health conditions.

Phenotype integration improves power and preserves specificity in biobank-based genetic studies of major depressive disorder

Biobanks often contain several phenotypes relevant to diseases such as major depressive disorder (MDD), with partly distinct genetic architectures. Researchers face complex tradeoffs between shallow (large sample size, low specificity/sensitivity) and deep (small sample size, high specificity/sensitivity) phenotypes, and the optimal choices are often unclear. Here we propose to integrate these phenotypes to combine the benefits of each. We use phenotype imputation to integrate information across hundreds of MDD-relevant phenotypes, which significantly increases genome-wide association study (GWAS) power and polygenic risk score (PRS) prediction accuracy of the deepest available MDD phenotype in UK Biobank, LifetimeMDD. We demonstrate that imputation preserves specificity in its genetic architecture using a novel PRS-based pleiotropy metric. We further find that integration via summary statistics also enhances GWAS power and PRS predictions, but can introduce nonspecific genetic effects depending on input. Our work provides a simple and scalable approach to improve genetic studies in large biobanks by integrating shallow and deep phenotypes.

Deep learning on electronic medical records identifies distinct subphenotypes of diabetic kidney disease driven by genetic variations in the Rho pathway

Kidney disease affects 50% of all diabetic patients; however, prediction of disease progression has been challenging due to inherent disease heterogeneity. We use deep learning to identify novel genetic signatures prognostically associated with outcomes. Using autoencoders and unsupervised clustering of electronic health record data on 1,372 diabetic kidney disease patients, we establish two clusters with differential prevalence of end-stage kidney disease. Exome-wide associations identify a novel variant in ARHGEF18, a Rho guanine exchange factor specifically expressed in glomeruli. Overexpression of ARHGEF18 in human podocytes leads to impairments in focal adhesion architecture, cytoskeletal dynamics, cellular motility, and RhoA/Rac1 activation. Mutant GEF18 is resistant to ubiquitin mediated degradation leading to pathologically increased protein levels. Our findings uncover the first known disease-causing genetic variant that affects protein stability of a cytoskeletal regulator through impaired degradation, a potentially novel class of expression quantitative trait loci that can be therapeutically targeted.

Identification of Novel Metabolic Subtypes Using Multi-Trait Limited Mixed Regression in the Chinese Population

The aggregation and interaction of metabolic risk factors leads to highly heterogeneous pathogeneses, manifestations, and outcomes, hindering risk stratification and targeted management. To deconstruct the heterogeneity, we used baseline data from phase II of the Fangshan Family-Based Ischemic Stroke Study (FISSIC), and a total of 4632 participants were included. A total of 732 individuals who did not have any component of metabolic syndrome (MetS) were set as a reference group, while 3900 individuals with metabolic abnormalities were clustered into subtypes using multi-trait limited mixed regression (MFMR). Four metabolic subtypes were identified with the dominant characteristics of abdominal obesity, hypertension, hyperglycemia, and dyslipidemia. Multivariate logistic regression showed that the hyperglycemia-dominant subtype had the highest coronary heart disease (CHD) risk (OR: 6.440, 95% CI: 3.177-13.977) and that the dyslipidemia-dominant subtype had the highest stroke risk (OR: 2.450, 95% CI: 1.250-5.265). Exome-wide association studies (EWASs) identified eight SNPs related to the dyslipidemia-dominant subtype with genome-wide significance, which were located in the genes APOA5, BUD13, ZNF259, and WNT4. Functional analysis revealed an enrichment of top genes in metabolism-related biological pathways and expression in the heart, brain, arteries, and kidneys. Our findings provide directions for future attempts at risk stratification and evidence-based management in populations with metabolic abnormalities from a systematic perspective.

Discerning asthma endotypes through comorbidity mapping

Asthma is a heterogeneous, complex syndrome, and identifying asthma endotypes has been challenging. We hypothesize that distinct endotypes of asthma arise in disparate genetic variation and life-time environmental exposure backgrounds, and that disease comorbidity patterns serve as a surrogate for such genetic and exposure variations. Here, we computationally discover 22 distinct comorbid disease patterns among individuals with asthma (asthma comorbidity subgroups) using diagnosis records for >151 M US residents, and re-identify 11 of the 22 subgroups in the much smaller UK Biobank. GWASs to discern asthma risk loci for individuals within each subgroup and in all subgroups combined reveal 109 independent risk loci, of which 52 are replicated in multi-ancestry meta-analysis across different ethnicity subsamples in UK Biobank, US BioVU, and BioBank Japan. Fourteen loci confer asthma risk in multiple subgroups and in all subgroups combined. Importantly, another six loci confer asthma risk in only one subgroup. The strength of association between asthma and each of 44 health-related phenotypes also varies dramatically across subgroups. This work reveals subpopulations of asthma patients distinguished by comorbidity patterns, asthma risk loci, gene expression, and health-related phenotypes, and so reveals different asthma endotypes.

Ten challenges for clinical translation in psychiatric genetics

Genome-wide association studies have identified hundreds of robust genetic associations underlying psychiatric disorders and provided important biological insights into disease onset and progression. There is optimism that genetic findings will pave the way to precision psychiatry by facilitating the development of more effective treatments and the identification of groups of patients that these treatments should be targeted toward. However, there are several challenges that must be addressed before genetic findings can be translated into the clinic. In this Perspective, we highlight ten challenges for the field of psychiatric genetics, focused on the robust and generalizable detection of genetic risk factors, improved definition and assessment of psychopathology and achieving better clinical indicators. We discuss recent advancements in the field that will improve the explanatory and predictive power of genetic data and ultimately contribute to improving the management and treatment of patients with a psychiatric disorder.

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.

Enhancing Discovery of Genetic Variants for Posttraumatic Stress Disorder Through Integration of Quantitative Phenotypes and Trauma Exposure Information

BACKGROUND

Posttraumatic stress disorder (PTSD) is heritable and a potential consequence of exposure to traumatic stress. Evidence suggests that a quantitative approach to PTSD phenotype measurement and incorporation of lifetime trauma exposure (LTE) information could enhance the discovery power of PTSD genome-wide association studies (GWASs).

METHODS

A GWAS on PTSD symptoms was performed in 51 cohorts followed by a fixed-effects meta-analysis (N = 182,199 European ancestry participants). A GWAS of LTE burden was performed in the UK Biobank cohort (N = 132,988). Genetic correlations were evaluated with linkage disequilibrium score regression. Multivariate analysis was performed using Multi-Trait Analysis of GWAS. Functional mapping and annotation of leading loci was performed with FUMA. Replication was evaluated using the Million Veteran Program GWAS of PTSD total symptoms.

RESULTS

GWASs of PTSD symptoms and LTE burden identified 5 and 6 independent genome-wide significant loci, respectively. There was a 72% genetic correlation between PTSD and LTE. PTSD and LTE showed largely similar patterns of genetic correlation with other traits, albeit with some distinctions. Adjusting PTSD for LTE reduced PTSD heritability by 31%. Multivariate analysis of PTSD and LTE increased the effective sample size of the PTSD GWAS by 20% and identified 4 additional loci. Four of these 9 PTSD loci were independently replicated in the Million Veteran Program.

CONCLUSIONS

Through using a quantitative trait measure of PTSD, we identified novel risk loci not previously identified using prior case-control analyses. PTSD and LTE have a high genetic overlap that can be leveraged to increase discovery power through multivariate methods.

Vascular Injury in the Zebrafish Tail Modulates Blood Flow and Peak Wall Shear Stress to Restore Embryonic Circular Network

Mechano-responsive signaling pathways enable blood vessels within a connected network to structurally adapt to partition of blood flow between organ systems. Wall shear stress (WSS) modulates endothelial cell proliferation and arteriovenous specification. Here, we study vascular regeneration in a zebrafish model by using tail amputation to disrupt the embryonic circulatory loop (ECL) at 3 days post fertilization (dpf). We observed a local increase in blood flow and peak WSS in the Segmental Artery (SeA) immediately adjacent to the amputation site. By manipulating blood flow and WSS via changes in blood viscosity and myocardial contractility, we show that the angiogenic Notch-ephrinb2 cascade is hemodynamically activated in the SeA to guide arteriogenesis and network reconnection. Taken together, ECL amputation induces changes in microvascular topology to partition blood flow and increase WSS-mediated Notch-ephrinb2 pathway, promoting new vascular arterial loop formation and restoring microcirculation.

Mendelian randomization highlights significant difference and genetic heterogeneity in clinically diagnosed Alzheimer's disease GWAS and self-report proxy phenotype GWAX

BACKGROUND

Until now, Mendelian randomization (MR) studies have investigated the causal association of risk factors with Alzheimer's disease (AD) using large-scale AD genome-wide association studies (GWAS), GWAS by proxy (GWAX), and meta-analyses of GWAS and GWAX (GWAS+GWAX) datasets. However, it currently remains unclear about the consistency of MR estimates across these GWAS, GWAX, and GWAS+GWAX datasets.

METHODS

Here, we first selected 162 independent educational attainment genetic variants as the potential instrumental variables (N = 405,072). We then selected one AD GWAS dataset (N = 63,926), two AD GWAX datasets (N = 314,278 and 408,942), and three GWAS+GWAX datasets (N = 388,324, 455,258, and 472,868). Finally, we conducted a MR analysis to evaluate the impact of educational attainment on AD risk across these datasets. Meanwhile, we tested the genetic heterogeneity of educational attainment genetic variants across these datasets.

RESULTS

In AD GWAS dataset, MR analysis showed that each SD increase in years of schooling (about 3.6 years) was significantly associated with 29% reduced AD risk (OR=0.71, 95% CI: 0.60-0.84, and P=1.02E-04). In AD GWAX dataset, MR analysis highlighted that each SD increase in years of schooling significantly increased 84% AD risk (OR=1.84, 95% CI: 1.59-2.13, and P=4.66E-16). Meanwhile, MR analysis suggested the ambiguous findings in AD GWAS+GWAX datasets. Heterogeneity test indicated evidence of genetic heterogeneity in AD GWAS and GWAX datasets.

CONCLUSIONS

We highlighted significant difference and genetic heterogeneity in clinically diagnosed AD GWAS and self-report proxy phenotype GWAX. Our MR findings are consistent with recent findings in AD genetic variants. Hence, the GWAX and GWAS+GWAX findings and MR findings from GWAX and GWAS+GWAX should be carefully interpreted and warrant further investigation using the AD GWAS dataset.

Recent innovations and in-depth aspects of post-genome wide association study (Post-GWAS) to understand the genetic basis of complex phenotypes

In the past decade, the high throughput and low cost of sequencing/genotyping approaches have led to the accumulation of a large amount of data from genome-wide association studies (GWASs). The first aim of this review is to highlight how post-GWAS analysis can be used make sense of the obtained associations. Novel directions for integrating GWAS results with other resources, such as somatic mutation, metabolite-transcript, and transcriptomic data, are also discussed; these approaches can help us move beyond each individual data point and provide valuable information about complex trait genetics. In addition, cross-phenotype association tests, when the loci detected by GWASs have significant associations with multiple traits, are reviewed to provide biologically informative results for use in real-time applications. This review also discusses the challenges of identifying interactions between genetic mutations (epistasis) and mutations of loci affecting more than one trait (pleiotropy) as underlying causes of cross-phenotype associations; these challenges can be overcome using post-GWAS analysis. Genetic similarities between phenotypes that can be revealed using post-GWAS analysis are also discussed. In summary, different methodologies of post-GWAS analysis are now available, enhancing the value of information obtained from GWAS results, and facilitating application in both humans and nonhuman species. However, precise methods still need to be developed to overcome challenges in the field and uncover the genetic underpinnings of complex traits.

Transcriptome-wide association study of treatment-resistant depression and depression subtypes for drug repurposing

Major depressive disorder (MDD) is the single largest contributor to global disability and up to 20-30% of patients do not respond to at least two antidepressants (treatment-resistant depression, TRD). This study leveraged imputed gene expression in TRD to perform a drug repurposing analysis. Among those with MDD, we defined TRD as having at least two antidepressant switches according to primary care records in UK Biobank (UKB). We performed a transcriptome-wide association study (TWAS) of TRD (n = 2165) vs healthy controls (n = 11,188) using FUSION and gene expression levels from 21 tissues. We identified compounds with opposite gene expression signatures (ConnectivityMap data) compared to our TWAS results using the Kolmogorov-Smirnov test, Spearman and Pearson correlation. As symptom patterns are routinely assessed in clinical practice and could be used to provide targeted treatments, we identified MDD subtypes associated with TRD in UKB and analysed them using the same pipeline described for TRD. Anxious MDD (n = 14,954) and MDD with weight gain (n = 4697) were associated with TRD. In the TWAS, two genes were significantly dysregulated (TMEM106B and ATP2A1 for anxious and weight gain MDD, respectively). A muscarinic receptor antagonist was identified as top candidate for repurposing in TRD; inhibition of heat shock protein 90 was the main mechanism of action identified for anxious MDD, while modulators of metabolism such as troglitazone showed promising results for MDD with weight gain. This was the first TWAS of TRD and associated MDD subtypes. Our results shed light on possible pharmacological approaches in individuals with difficult-to-treat depression.

Multitrait GWAS to connect disease variants and biological mechanisms

Genome-wide association studies (GWASs) have uncovered a wealth of associations between common variants and human phenotypes. Here, we present an integrative analysis of GWAS summary statistics from 36 phenotypes to decipher multitrait genetic architecture and its link with biological mechanisms. Our framework incorporates multitrait association mapping along with an investigation of the breakdown of genetic associations into clusters of variants harboring similar multitrait association profiles. Focusing on two subsets of immunity and metabolism phenotypes, we then demonstrate how genetic variants within clusters can be mapped to biological pathways and disease mechanisms. Finally, for the metabolism set, we investigate the link between gene cluster assignment and the success of drug targets in randomized controlled trials.

Genetic Influences on Disease Subtypes

Disease classification, or nosology, was historically driven by careful examination of clinical features of patients. As technologies to measure and understand human phenotypes advanced, so too did classifications of disease, and the advent of genetic data has led to a surge in genetic subtyping in the past decades. Although the fundamental process of refining disease definitions and subtypes is shared across diverse fields, each field is driven by its own goals and technological expertise, leading to inconsistent and conflicting definitions of disease subtypes. Here, we review several classical and recent subtypes and subtyping approaches and provide concrete definitions to delineate subtypes. In particular, we focus on subtypes with distinct causal disease biology, which are of primary interest to scientists, and subtypes with pragmatic medical benefits, which are of primary interest to physicians. We propose genetic heterogeneity as a gold standard for establishing biologically distinct subtypes of complex polygenic disease. We focus especially on methods to find and validate genetic subtypes, emphasizing common pitfalls and how to avoid them.

Leveraging eQTLs to identify individual-level tissue of interest for a complex trait

Genetic predisposition for complex traits often acts through multiple tissues at different time points during development. As a simple example, the genetic predisposition for obesity could be manifested either through inherited variants that control metabolism through regulation of genes expressed in the brain, or that control fat storage through dysregulation of genes expressed in adipose tissue, or both. Here we describe a statistical approach that leverages tissue-specific expression quantitative trait loci (eQTLs) corresponding to tissue-specific genes to prioritize a relevant tissue underlying the genetic predisposition of a given individual for a complex trait. Unlike existing approaches that prioritize relevant tissues for the trait in the population, our approach probabilistically quantifies the tissue-wise genetic contribution to the trait for a given individual. We hypothesize that for a subgroup of individuals the genetic contribution to the trait can be mediated primarily through a specific tissue. Through simulations using the UK Biobank, we show that our approach can predict the relevant tissue accurately and can cluster individuals according to their tissue-specific genetic architecture. We analyze body mass index (BMI) and waist to hip ratio adjusted for BMI (WHRadjBMI) in the UK Biobank to identify subgroups of individuals whose genetic predisposition act primarily through brain versus adipose tissue, and adipose versus muscle tissue, respectively. Notably, we find that these individuals have specific phenotypic features beyond BMI and WHRadjBMI that distinguish them from random individuals in the data, suggesting biological effects of tissue-specific genetic contribution for these traits.

A model and test for coordinated polygenic epistasis in complex traits

Interactions between genetic variants-epistasis-is pervasive in model systems and can profoundly impact evolutionary adaption, population disease dynamics, genetic mapping, and precision medicine efforts. In this work, we develop a model for structured polygenic epistasis, called coordinated epistasis (CE), and prove that several recent theories of genetic architecture fall under the formal umbrella of CE. Unlike standard epistasis models that assume epistasis and main effects are independent, CE captures systematic correlations between epistasis and main effects that result from pathway-level epistasis, on balance skewing the penetrance of genetic effects. To test for the existence of CE, we propose the even-odd (EO) test and prove it is calibrated in a range of realistic biological models. Applying the EO test in the UK Biobank, we find evidence of CE in 18 of 26 traits spanning disease, anthropometric, and blood categories. Finally, we extend the EO test to tissue-specific enrichment and identify several plausible tissue-trait pairs. Overall, CE is a dimension of genetic architecture that can capture structured, systemic forms of epistasis in complex human traits.

Polygenic risk modeling with latent trait-related genetic components

Polygenic risk models have led to significant advances in understanding complex diseases and their clinical presentation. While polygenic risk scores (PRS) can effectively predict outcomes, they do not generally account for disease subtypes or pathways which underlie within-trait diversity. Here, we introduce a latent factor model of genetic risk based on components from Decomposition of Genetic Associations (DeGAs), which we call the DeGAs polygenic risk score (dPRS). We compute DeGAs using genetic associations for 977 traits and find that dPRS performs comparably to standard PRS while offering greater interpretability. We show how to decompose an individual's genetic risk for a trait across DeGAs components, with examples for body mass index (BMI) and myocardial infarction (heart attack) in 337,151 white British individuals in the UK Biobank, with replication in a further set of 25,486 non-British white individuals. We find that BMI polygenic risk factorizes into components related to fat-free mass, fat mass, and overall health indicators like physical activity. Most individuals with high dPRS for BMI have strong contributions from both a fat-mass component and a fat-free mass component, whereas a few "outlier" individuals have strong contributions from only one of the two components. Overall, our method enables fine-scale interpretation of the drivers of genetic risk for complex traits.

Characterizing the effect of background selection on the polygenicity of brain-related traits

BACKGROUND

Genome-wide association studies (GWAS) have demonstrated that psychopathology phenotypes are affected by many risk alleles with small effect (polygenicity). It is unclear how ubiquitously evolutionary pressures influence the genetic architecture of these traits.

METHODS

We partitioned SNP heritability to assess the contribution of background (BGS) and positive selection, Neanderthal local ancestry, functional significance, and genotype networks in 75 brain-related traits (8411 ≤ N ≤ 1,131,181, mean N = 205,289). We applied binary annotations by dichotomizing each measure based on top 2%, 1%, and 0.5% of all scores genome-wide. Effect size distribution features were calculated using GENESIS. We tested the relationship between effect size distribution descriptive statistics and natural selection. In a subset of traits, we explore the inclusion of diagnostic heterogeneity (e.g., number of diagnostic combinations and total symptoms) in the tested relationship.

RESULTS

SNP-heritability was enriched (false discovery rate q < 0.05) for loci with elevated BGS (7 phenotypes) and in genic (34 phenotypes) and loss-of-function (LoF)-intolerant regions (67 phenotypes). These effects were strongest in GWAS of schizophrenia (1.90-fold BGS, 1.16-fold genic, and 1.92-fold LoF), educational attainment (1.86-fold BGS, 1.12-fold genic, and 1.79-fold LoF), and cognitive performance (2.29-fold BGS, 1.12-fold genic, and 1.79-fold LoF). BGS (top 2%) significantly predicted effect size variance for trait-associated loci (σ2 parameter) in 75 brain-related traits (β = 4.39 × 10-5, p = 1.43 × 10-5, model r2 = 0.548). Considering the number of DSM-5 diagnostic combinations per psychiatric disorder improved model fit (σ2 ~ BTop2% × Genic × diagnostic combinations; model r2 = 0.661).

CONCLUSIONS

Brain-related phenotypes with larger variance in risk locus effect sizes are associated with loci under BGS. We show exploratory results suggesting that diagnostic complexity may also contribute to the increased polygenicity of psychiatric disorders.

Module Analysis Using Single-Patient Differential Expression Signatures Improves the Power of Association Studies for Alzheimer's Disease

The causal mechanism of Alzheimer's disease is extremely complex. Achieving great statistical power in association studies usually requires a large number of samples. In this work, we illustrated a different strategy to identify AD risk genes by clustering AD patients into modules based on their single-patient differential expression signatures. The evaluation suggested that our method could enrich AD patients with similar clinical manifestations. Applying this to a cohort of only 310 AD patients, we identified 174 AD risk loci at a strict threshold of empirical p < 0.05, while only two loci were identified using all the AD patients. As an evaluation, we collected 23 AD risk genes reported in a recent large-scale meta-analysis and found that 18 of them were rediscovered by association studies using clustered AD patients, while only three of them were rediscovered using all AD patients. Functional annotation suggested that AD-associated genetic variants mainly disturbed neuronal/synaptic function. Our results suggested module analysis helped to enrich AD patients affected by the common risk variants.

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.

Genetic risk factors of ME/CFS: a critical review

Myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) is a complex multisystem illness that lacks effective therapy and a biomedical understanding of its causes. Despite a prevalence of ∼0.2-0.4% and its high public health burden, and evidence that it has a heritable component, ME/CFS has not yet benefited from the advances in technology and analytical tools that have improved our understanding of many other complex diseases. Here we critically review existing evidence that genetic factors alter ME/CFS risk before concluding that most ME/CFS candidate gene associations are not replicated by the larger CFS cohort within the UK Biobank. Multiple genome-wide association studies of this cohort also have not yielded consistently significant associations. Ahead of upcoming larger genome-wide association studies, we discuss how these could generate new lines of enquiry into the DNA variants, genes and cell types that are causally involved in ME/CFS disease.

Leveraging correlations between variants in polygenic risk scores to detect heterogeneity in GWAS cohorts

Evidence from both GWAS and clinical observation has suggested that certain psychiatric, metabolic, and autoimmune diseases are heterogeneous, comprising multiple subtypes with distinct genomic etiologies and Polygenic Risk Scores (PRS). However, the presence of subtypes within many phenotypes is frequently unknown. We present CLiP (Correlated Liability Predictors), a method to detect heterogeneity in single GWAS cohorts. CLiP calculates a weighted sum of correlations between SNPs contributing to a PRS on the case/control liability scale. We demonstrate mathematically and through simulation that among i.i.d. homogeneous cases generated by a liability threshold model, significant anti-correlations are expected between otherwise independent predictors due to ascertainment on the hidden liability score. In the presence of heterogeneity from distinct etiologies, confounding by covariates, or mislabeling, these correlation patterns are altered predictably. We further extend our method to two additional association study designs: CLiP-X for quantitative predictors in applications such as transcriptome-wide association, and CLiP-Y for quantitative phenotypes, where there is no clear distinction between cases and controls. Through simulations, we demonstrate that CLiP and its extensions reliably distinguish between homogeneous and heterogeneous cohorts when the PRS explains as low as 3% of variance on the liability scale and cohorts comprise 50, 000 − 100, 000 samples, an increasingly practical size for modern GWAS. We apply CLiP to heterogeneity detection in schizophrenia cohorts totaling > 50, 000 cases and controls collected by the Psychiatric Genomics Consortium. We observe significant heterogeneity in mega-analysis of the combined PGC data (p-value 8.54 × 0⁻⁴), as well as in individual cohorts meta-analyzed using Fisher’s method (p-value 0.03), based on significantly associated variants. We also apply CLiP-Y to detect heterogeneity in neuroticism in over 10, 000 individuals from the UK Biobank and detect heterogeneity with a p-value of 1.68 × 10⁻⁹. Scores were not significantly reduced when partitioning by known subclusters (“Depression” and “Worry”), suggesting that these factors are not the primary source of observed heterogeneity.

Several traits, such as bipolar disease, are known to be heterogeneous and comprise distinct subtypes with unique genomic associations. For other traits such as schizophrenia, heterogeneity may be suspected, but specific subtypes are less well characterized. Furthermore, conventional mixture model-based methods to detect subtypes in genome-wide association data struggle with the high polygenicity of complex traits. We propose CLiP (Correlated Liability Predictors), a method that does not identify subtype-specific effects, but is very well-powered to detect heterogeneity of any kind within the very weak signals of GWAS. CLiP serves as a method to flag particular phenotypes for potential further study into the genomic factors driving heterogeneity, as well as a means to evaluate the transferability of polygenic risk scores. We also develop extensions of CLiP applicable to scoring heterogeneity in quantitative phenotypes and quantitative predictors such as gene expression. We apply CLiP to scoring heterogeneity in schizophrenia cohorts from the Psychiatric Genomics Consortium and neuroticism in individuals in the UK Biobank and find significant heterogeneity in both phenotypes, warranting further study.

Emerging phenotyping strategies will advance our understanding of psychiatric genetics

Over the last decade, genome-wide association studies of psychiatric disorders have identified numerous significant loci. Whereas these studies initially depended on cohorts ascertained for specific disorders, there has been a gradual shift in the ascertainment strategy toward population-based cohorts for which both genotype and heterogeneous phenotypic information are available. One of the advantages of population-based cohorts is that, in addition to clinical diagnoses and various proxies for diagnoses ('minimal phenotyping'), many of them also provide non-clinical phenotypes, including putative endophenotypes, that can be used to study domains of normal function in addition to, or instead of, clinical diagnoses. By studying endophenotypes it is possible to both dissect psychiatric disorders ('splitting') and to combine multiple phenotypes ('clumping'), which can either reinforce or challenge traditional diagnostic categories. Such endophenotypes may also permit a deeper exploration of the neurobiology of psychiatric disorders. A coordinated effort to fully exploit the potential of endophenotypes is overdue.

Genetics and epigenetics of allergy

PURPOSE OF REVIEW

Allergic diseases are prototypic examples for gene × environment-wide interactions. This review considers the current evidence for genetic and epigenetic mechanisms in allergic diseases and highlights barriers and facilitators for the implementation of these novel tools both for research and clinical practice.

RECENT FINDINGS

The value of whole-genome sequencing studies and the use of polygenic risk score analysis in homogeneous well characterized populations are currently being tested. Epigenetic mechanisms are known to play a crucial role in the pathogenesis of allergic disorders, especially through mediating the effects of the environmental factors, well recognized risk modifiers. There is emerging evidence for the immune-modulatory role of probiotics through epigenetic changes. Direct or indirect targeting of epigenetic mechanisms affect expression of the genes favouring the development of allergic diseases and can improve tissue biology. The ability to specifically edit the epigenome, especially using the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)/Cas9 technology, holds the promise of enhancing understanding of how epigenetic modifications function and enabling manipulation of cell phenotype for research or therapeutic purposes.

SUMMARY

Additional research in the role of genetic and epigenetic mechanisms in relation to allergic diseases' endotypes is needed. An international project characterizing the human epigenome in relation to allergic diseases is warranted.

A Robust Method Uncovers Significant Context-Specific Heritability in Diverse Complex Traits

Gene-environment interactions (GxE) can be fundamental in applications ranging from functional genomics to precision medicine and is a conjectured source of substantial heritability. However, unbiased methods to profile GxE genome-wide are nascent and, as we show, cannot accommodate general environment variables, modest sample sizes, heterogeneous noise, and binary traits. To address this gap, we propose a simple, unifying mixed model for gene-environment interaction (GxEMM). In simulations and theory, we show that GxEMM can dramatically improve estimates and eliminate false positives when the assumptions of existing methods fail. We apply GxEMM to a range of human and model organism datasets and find broad evidence of context-specific genetic effects, including GxSex, GxAdversity, and GxDisease interactions across thousands of clinical and molecular phenotypes. Overall, GxEMM is broadly applicable for testing and quantifying polygenic interactions, which can be useful for explaining heritability and invaluable for determining biologically relevant environments.

Advances in asthma and allergic disease genetics: Is bigger always better?

This review focuses on genome-wide association studies (GWASs) of asthma and allergic diseases published between January 1, 2018, and June 30, 2019. During this time period, there were 38 GWASs reported in 19 articles, including the largest performed to date for many of these conditions. Overall, we learned that childhood-onset asthma is associated with the most independent loci compared with other defined groups of asthma and allergic disease cases; adult-onset asthma and moderate-to-severe asthma are associated with fewer genes, which are largely a subset of those associated with childhood-onset asthma. There is significant genetic overlap between asthma and allergic diseases, particularly with respect to childhood-onset asthma, which involves genes that reflect the importance of barrier function biology, and to HLA region genes, which are the most frequently associated genes overall in both groups of diseases. Although the largest GWASs in African American and Latino/Hispanic populations were reported during this period, they are still significantly underpowered compared with studies reported in populations of European ancestry, highlighting the need for larger studies, particularly in patients with childhood-onset asthma and allergic diseases, in these important populations that carry the greatest burden of disease.

A comprehensive study of metabolite genetics reveals strong pleiotropy and heterogeneity across time and context

Genetic studies of metabolites have identified thousands of variants, many of which are associated with downstream metabolic and obesogenic disorders. However, these studies have relied on univariate analyses, reducing power and limiting context-specific understanding. Here we aim to provide an integrated perspective of the genetic basis of metabolites by leveraging the Finnish Metabolic Syndrome In Men (METSIM) cohort, a unique genetic resource which contains metabolic measurements, mostly lipids, across distinct time points as well as information on statin usage. We increase effective sample size by an average of two-fold by applying the Covariates for Multi-phenotype Studies (CMS) approach, identifying 588 significant SNP-metabolite associations, including 228 new associations. Our analysis pinpoints a small number of master metabolic regulator genes, balancing the relative proportion of dozens of metabolite levels. We further identify associations to changes in metabolic levels across time as well as genetic interactions with statin at both the master metabolic regulator and genome-wide level.