Statistical Analysis of Multiple Phenotypes in Genetic Epidemiologic Studies: From Cross-Phenotype Associations to Pleiotropy

Identifying pleiotropic genes via the composite test amidst the complexity of polygenic traits

Identifying the causal relationship between genotype and phenotype is essential to expanding our understanding of the gene regulatory network spanning the molecular level to perceptible traits. A pleiotropic gene can act as a central hub in the network, influencing multiple outcomes. Identifying such a gene involves testing under a composite null hypothesis where the gene is associated with, at most, one trait. Traditional methods such as meta-analyses of top-hit $P$-values and sequential testing of multiple traits have been proposed, but these methods fail to consider the background of genome-wide signals. Since Huang's composite test produces uniformly distributed $P$-values for genome-wide variants under the composite null, we propose a gene-level pleiotropy test that entails combining the aforementioned method with the aggregated Cauchy association test. A polygenic trait involves multiple genes with different functions to co-regulate mechanisms. We show that polygenicity should be considered when identifying pleiotropic genes; otherwise, the associations polygenic traits initiate will give rise to false positives. In this study, we constructed gene-trait functional modules using the results of the proposed pleiotropy tests. Our analysis suite was implemented as an R package PGCtest. We demonstrated the proposed method with an application study of the Taiwan Biobank database and identified functional modules comprising specific genes and their co-regulated traits.

Multivariate genome-wide associations for immune traits in two maternal pig lines

BACKGROUND

Immune traits are considered to serve as potential biomarkers for pig's health. Medium to high heritabilities have been observed for some of the immune traits suggesting genetic variability of these phenotypes. Consideration of previously established genetic correlations between immune traits can be used to identify pleiotropic genetic markers. Therefore, genome-wide association study (GWAS) approaches are required to explore the joint genetic foundation for health biomarkers. Usually, GWAS explores phenotypes in a univariate (uv), trait-by-trait manner. Besides two uv GWAS methods, four multivariate (mv) GWAS approaches were applied on combinations out of 22 immune traits for Landrace (LR) and Large White (LW) pig lines.

RESULTS

In total 433 (LR: 351, LW: 82) associations were identified with the uv approach implemented in PLINK and a Bayesian linear regression uv approach (BIMBAM) software. Single Nucleotide Polymorphisms (SNPs) that were identified with both uv approaches (n = 32) were mostly associated with immune traits such as haptoglobin, red blood cell characteristics and cytokines, and were located in protein-coding genes. Mv GWAS approaches detected 647 associations for different mv immune trait combinations which were summarized to 133 Quantitative Trait Loci (QTL). SNPs for different trait combinations (n = 66) were detected with more than one mv method. Most of these SNPs are associated with red blood cell related immune trait combinations. Functional annotation of these QTL revealed 453 immune-relevant protein-coding genes. With uv methods shared markers were not observed between the breeds, whereas mv approaches were able to detect two conjoint SNPs for LR and LW. Due to unmapped positions for these markers, their functional annotation was not clarified.

CONCLUSIONS

This study evaluated the joint genetic background of immune traits in LR and LW piglets through the application of various uv and mv GWAS approaches. In comparison to uv methods, mv methodologies identified more significant associations, which might reflect the pleiotropic background of the immune system more accurately. In genetic research of complex traits, the SNP effects are generally small. Furthermore, one genetic variant can affect several correlated immune traits at the same time, termed pleiotropy. As mv GWAS methods consider strong dependencies among traits, the power to detect SNPs can be boosted. Both methods revealed immune-relevant potential candidate genes. Our results indicate that one single test is not able to detect all the different types of genetic effects in the most powerful manner and therefore, the methods should be applied complementary.

Simultaneous modeling of multivariate heterogeneous responses and heteroskedasticity via a two-stage composite likelihood

Multivariate heterogeneous responses and heteroskedasticity have attracted increasing attention in recent years. In genome-wide association studies, effective simultaneous modeling of multiple phenotypes would improve statistical power and interpretability. However, a flexible common modeling system for heterogeneous data types can pose computational difficulties. Here we build upon a previous method for multivariate probit estimation using a two-stage composite likelihood that exhibits favorable computational time while retaining attractive parameter estimation properties. We extend this approach to incorporate multivariate responses of heterogeneous data types (binary and continuous), and possible heteroskedasticity. Although the approach has wide applications, it would be particularly useful for genomics, precision medicine, or individual biomedical prediction. Using a genomics example, we explore statistical power and confirm that the approach performs well for hypothesis testing and coverage percentages under a wide variety of settings. The approach has the potential to better leverage genomics data and provide interpretable inference for pleiotropy, in which a locus is associated with multiple traits.

Benchmarking of univariate pleiotropy detection methods applied to epilepsy

Pleiotropy is a widespread phenomenon that may increase insight into the etiology of biological and disease traits. Since genome-wide association studies frequently provide information on a single trait only, only univariate pleiotropy detection methods are applicable, with yet unknown comparative performance. Here, we compared five such methods with respect to their ability to detect pleiotropy, including meta-analysis, ASSET, cFDR, CPBayes, and PLACO, by performing extended computer simulations that varied the underlying etiological model for pleiotropy for a pair of traits, including the number of causal variants, degree of traits' overlap, effect sizes as well as trait prevalence, and varying sample sizes. Our results indicate that ASSET provides the best trade-off between power and protection against false positives. We then applied ASSET to a previously published ILAE consortium dataset on complex epilepsies, comprising genetic generalized epilepsy and focal epilepsy cases and corresponding controls. We identified a novel candidate locus at 17q21.32 and confirmed locus 2q24.3, previously identified to act pleiotropically on both epilepsy subtypes by a mega-analysis. Functional annotation, tissue-specific expression and regulatory function analysis as well as Bayesian co-localization analysis corroborated this result, rendering 17q21.32 a worthwhile candidate for follow-up studies on pleiotropy in epilepsies. This article is protected by copyright. All rights reserved.

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.

Associating complex traits with genetic variants: polygenic risk scores, pleiotropy and endophenotypes

Genotype-phenotype causal modeling has evolved significantly since Johannsen's and Wright's original designs were published. The development of genomewide assays to interrogate and detect possible causal variants associated with complex traits has expanded the scope of genotype-phenotype research considerably. Clusters of causal variants discovered by genomewide assays and associated with complex traits have been used to develop polygenic risk scores to predict clinical diagnoses of multidimensional human disorders. However, genomewide investigations have met with many challenges to their research designs and statistical complexities which have hindered the reliability and validity of their predictions. Findings linked to differences in heritability estimates between causal clusters and complex traits among unrelated individuals remain a research area of some controversy. Causal models developed from case-control studies as opposed to experiments, as well as other issues concerning the genotype-phenotype causal model and the extent to which various forms of pleiotropy and the concept of the endophenotype add to its complexity, will be reviewed.

How Well Can Multivariate and Univariate GWAS Distinguish Between True and Spurious Pleiotropy?

Quantification of the simultaneous contributions of loci to multiple traits, a phenomenon called pleiotropy, is facilitated by the increased availability of high-throughput genotypic and phenotypic data. To understand the prevalence and nature of pleiotropy, the ability of multivariate and univariate genome-wide association study (GWAS) models to distinguish between pleiotropic and non-pleiotropic loci in linkage disequilibrium (LD) first needs to be evaluated. Therefore, we used publicly available maize and soybean genotypic data to simulate multiple pairs of traits that were either (i) controlled by quantitative trait nucleotides (QTNs) on separate chromosomes, (ii) controlled by QTNs in various degrees of LD with each other, or (iii) controlled by a single pleiotropic QTN. We showed that multivariate GWAS could not distinguish between QTNs in LD and a single pleiotropic QTN. In contrast, a unique QTN detection rate pattern was observed for univariate GWAS whenever the simulated QTNs were in high LD or pleiotropic. Collectively, these results suggest that multivariate and univariate GWAS should both be used to infer whether or not causal mutations underlying peak GWAS associations are pleiotropic. Therefore, we recommend that future studies use a combination of multivariate and univariate GWAS models, as both models could be useful for identifying and narrowing down candidate loci with potential pleiotropic effects for downstream biological experiments.

Discovery and Mediation Analysis of Cross-Phenotype Associations Between Asthma and Body Mass Index in 12q13.2

Twin studies suggest that shared genetics contributes to the comorbidity of asthma and obesity, but candidate-gene studies provide limited evidence of pleiotropy. We conducted genome-wide association analyses of asthma and body mass index (BMI; weight (kg)/height (m)2)) among 305,945 White British subjects recruited into the UK Biobank in 2006-2010. We searched for overlapping signals and conducted mediation analyses on genome-wide-significant cross-phenotype associations, assessing moderation by sex and age at asthma diagnosis, and adjusting for confounders of the asthma-BMI relationship. We identified a genome-wide-significant cross-phenotype association at rs705708 (asthma odds ratio = 1.05, 95% confidence interval: 1.03, 1.07; P = 7.20 × 10-9; and BMI β = -0.065, 95% confidence interval: -0.087, -0.042; P = 1.30 × 10-8). rs705708 resides on 12q13.2, which harbors 9 other asthma- and BMI-associated variants (all P < 5 × 10-5 for asthma; all but one P < 5 × 10-5 for BMI). Follow-up analyses of rs705708 show that most of the BMI association occurred independently of asthma, with consistent magnitude between men and women and persons with and without asthma, irrespective of age at diagnosis; the asthma association was stronger for childhood versus adult asthma; and both associations remained after confounder adjustment. This suggests that 12q13.2 displays pleiotropy for asthma and BMI. Upon further characterization, 12q13.2 might provide a target for interventions that simultaneously prevent or treat asthma and obesity.

Pleiotropy and Cross-Disorder Genetics Among Psychiatric Disorders

Genome-wide analyses of common and rare genetic variations have documented the heritability of major psychiatric disorders, established their highly polygenic genetic architecture, and identified hundreds of contributing variants. In recent years, these studies have illuminated another key feature of the genetic basis of psychiatric disorders: the important role and pervasive nature of pleiotropy. It is now clear that a substantial fraction of genetic influences on psychopathology transcend clinical diagnostic boundaries. In this review, we summarize evidence in psychiatry for pleiotropy at multiple levels of analysis: from overall genome-wide correlation to biological pathways and down to the level of individual loci. We examine underlying mechanisms of observed pleiotropy, including genetic effects on neurodevelopment, diverse actions of regulatory elements, mediated effects, and spurious associations of genomic variation with multiple phenotypes. We conclude with an exploration of the implications of pleiotropy for understanding the genetic basis of psychiatric disorders, informing nosology, and advancing the aims of precision psychiatry and genomic medicine.

Genetic Pleiotropy of Bone-Related Phenotypes: Insights from Osteoporosis

PURPOSE OF REVIEW

We summarize recent evidence on the shared genetics within and outside the musculoskeletal system (mostly related to bone density and osteoporosis).

RECENT FINDINGS

Osteoporosis is determined by an interplay between multiple genetic and environmental factors. Significant progress has been made regarding its genetic background revealing a number of robustly validated loci and respective pathways. However, pleiotropic factors affecting bone and other tissues are not well understood. The analytical methods proposed to test for potential associations between genetic variants and multiple phenotypes can be applied to bone-related data. A number of recent genetic studies have shown evidence of pleiotropy between bone density and other different phenotypes (traits, conditions, or diseases), within and outside the musculoskeletal system. Power benefits of combining correlated phenotypes, as well as unbiased discovery, make these studies promising. Studies in humans are supported by evidence from animal models. Drug development and repurposing should benefit from the pleiotropic approach. We believe that future studies should take into account shared genetics between the bone and related traits.

The E Is in the G: Gene-Environment-Trait Correlations and Findings From Genome-Wide Association Studies

Genome-wide association studies (GWASs) have shown that pleiotropy is widespread (i.e., the same genetic variants affect multiple traits) and that complex traits are polygenic (i.e., affected by many genetic variants with very small effect sizes). However, despite the growing number of GWASs, the possible contribution of gene-environment correlations (rGEs) to pleiotropy and polygenicity has been mostly ignored. rGEs can lead to environmentally mediated pleiotropy or gene-environment-trait correlations (rGETs), given that an environment that is affected by one genetically influenced phenotype, can in turn affect a different phenotype. By adding correlations with environmentally mediated genetic variants, rGETs can contribute to polygenicity. Socioeconomic status (SES) and the experience of stressful life events may, for example, be involved in rGETs. Both are genetically influenced and have been associated with a myriad of physical and mental disorders. As a result, GWASs of these disorders may find the genetic correlates of SES and stressful life events. Consequently, some of the genetic correlates of physical and mental disorders may be modified by public policy that affects environments such as SES and stressful life events. Thus, identifying rGETs can shed light on findings from GWASs and have important implications for public health.

Assessing pleiotropy and mediation in genetic loci associated with chronic obstructive pulmonary disease

Genetic association studies have increasingly recognized variant effects on multiple phenotypes. Chronic obstructive pulmonary disease (COPD) is a heterogeneous disease with environmental and genetic causes. Multiple genetic variants have been associated with COPD, many of which show significant associations to additional phenotypes. However, it is unknown if these associations represent biological pleiotropy or if they exist through correlation of related phenotypes ("mediated pleiotropy"). Using 6,670 subjects from the COPDGene study, we describe the association of known COPD susceptibility loci with other COPD-related phenotypes and distinguish if these act directly on the phenotypes (i.e., biological pleiotropy) or if the association is due to correlation (i.e., mediated pleiotropy). We identified additional associated phenotypes for 13 of 25 known COPD loci. Tests for pleiotropy between genotype and associated outcomes were significant for all loci. In cases of significant pleiotropy, we performed mediation analysis to test if SNPs had a direct association to phenotype. Most loci showed a mediated effect through the hypothesized causal pathway. However, many loci also had direct associations, suggesting causal explanations (i.e., emphysema leading to reduced lung function) are incomplete. Our results highlight the high degree of pleiotropy in complex disease-associated loci and provide novel insights into the mechanisms underlying COPD.