Multiple Related Phenotypes Research Articles

Identifying genetic loci and phenomic associations of substance use traits: A multi-trait analysis of GWAS (MTAG) study.

Genome-wide association studies (GWAS) of opioid use disorder (OUD) and cannabis use disorder (CUD) have lagged behind those of alcohol use disorder (AUD) and smoking, where many more loci have been identified. We sought to identify novel loci for substance use traits (SUTs) in both African- (AFR) and European- (EUR) ancestry individuals to enhance our understanding of the traits' genetic architecture. We used multi-trait analysis of GWAS (MTAG) to analyze four SUTs in EUR subjects (OUD, CUD, AUD and smoking initiation [SMKinitiation]), and three SUTs in AFR subjects (OUD, AUD and smoking trajectory [SMKtrajectory]). We conducted gene-set and protein-protein interaction analyses and calculated polygenic risk scores (PRS) in two independent samples. This study was conducted in the United States. A total of 5692 EUR and 4918 AFR individuals in the Yale-Penn sample and 29 054 EUR and 10 265 AFR individuals in the Penn Medicine BioBank sample. MTAG identified genome-wide significant (GWS) single nucleotide polymorphisms (SNPs) for all four traits in EUR: 41 SNPs in 36 loci for OUD; 74 SNPs in 60 loci for CUD; 63 SNPs in 52 loci for AUD; and 183 SNPs in 144 loci for SMKinitiation. MTAG also identified GWS SNPs in AFR: 2 SNPs in 2 loci for OUD; 3 SNPs in 3 loci for AUD; and 1 SNP in 1 locus for SMKtrajectory. In the Yale-Penn sample, the MTAG-derived PRS consistently yielded more significant associations with both the corresponding substance use disorder diagnosis and multiple related phenotypes than the GWAS-derived PRS. Multi-trait analysis of genome-wide association studies boosted the number of loci found for substance use traits, identifying genes not previously linked to any substance, and increased the power of polygenic risk scores. Multi-trait analysis of genome-wide association studies can be used to identify novel associations for substance use, especially those for which the samples are smaller than those for historically legal substances.

Key phenotypes related to wheat grain yield in a two‐site multicultivar test

AbstractYield performances of wheat (Triticum aestivum L.) cultivars are decided by the combined effects of key phenotypes. However, it is unclear which phenotypes are the key phenotypes that determine yield. We conducted a two‐site multicultivar test in Henan and Anhui in China during the 2016‒2017 wheat growing season. The objectives were to assess the relationships between grain yield and other related phenotypes and reveal the key phenotypes related to grain yield. Grain yield and nineteen phenotypes were assessed in a field trial consisting of 28 wheat cultivars. Multivariate statistical analysis methods (e.g., principal component analysis, partial least squares model, structural equation model) were used to reveal the relationships between grain yield and other related phenotypes. Wheat phenotypes changed as a result of the differences in environment (site) and genotype (cultivar). Eleven phenotypes were identified as the key phenotypes closely related to grain yield: biomass, harvest index, spike density, leaf nitrate reductase, net photosynthetic rate, canopy temperature, leaf chlorophyll, leaf area index, stomatal conductance, carboxylation efficiency, and transpiration rate. For analyzing the mechanisms of yield differences, multiple related phenotypes (rather than a single phenotype) should be considered.

Spatial rank-based multifactor dimensionality reduction to detect gene\u2013gene interactions for multivariate phenotypes

BackgroundIdentifying interaction effects between genes is one of the main tasks of genome-wide association studies aiming to shed light on the biological mechanisms underlying complex diseases. Multifactor dimensionality reduction (MDR) is a popular approach for detecting gene–gene interactions that has been extended in various forms to handle binary and continuous phenotypes. However, only few multivariate MDR methods are available for multiple related phenotypes. Current approaches use Hotelling’s T2 statistic to evaluate interaction models, but it is well known that Hotelling’s T2 statistic is highly sensitive to heavily skewed distributions and outliers.ResultsWe propose a robust approach based on nonparametric statistics such as spatial signs and ranks. The new multivariate rank-based MDR (MR-MDR) is mainly suitable for analyzing multiple continuous phenotypes and is less sensitive to skewed distributions and outliers. MR-MDR utilizes fuzzy k-means clustering and classifies multi-locus genotypes into two groups. Then, MR-MDR calculates a spatial rank-sum statistic as an evaluation measure and selects the best interaction model with the largest statistic. Our novel idea lies in adopting nonparametric statistics as an evaluation measure for robust inference. We adopt tenfold cross-validation to avoid overfitting. Intensive simulation studies were conducted to compare the performance of MR-MDR with current methods. Application of MR-MDR to a real dataset from a Korean genome-wide association study demonstrated that it successfully identified genetic interactions associated with four phenotypes related to kidney function. The R code for conducting MR-MDR is available at https://github.com/statpark/MR-MDR.ConclusionsIntensive simulation studies comparing MR-MDR with several current methods showed that the performance of MR-MDR was outstanding for skewed distributions. Additionally, for symmetric distributions, MR-MDR showed comparable power. Therefore, we conclude that MR-MDR is a useful multivariate non-parametric approach that can be used regardless of the phenotype distribution, the correlations between phenotypes, and sample size.

A two-step approach to testing overall effect of gene-environment interaction for multiple phenotypes.

While gene-environment (GxE) interactions contribute importantly to many different phenotypes, detecting such interactions requires well-powered studies and has proven difficult. To address this, we combine two approaches to improve GxE power: simultaneously evaluating multiple phenotypes and using a two-step analysis approach. Previous work shows that the power to identify a main genetic effect can be improved by simultaneously analyzing multiple related phenotypes. For a univariate phenotype, two-step methods produce higher power for detecting a GxE interaction compared to single step analysis. Therefore, we propose a two-step approach to test for an overall GxE effect for multiple phenotypes. Using simulations we demonstrate that, when more than one phenotype has GxE effect (i.e. GxE pleiotropy), our approach offers substantial gain in power (18-43%) to detect an aggregate-level GxE effect for a multivariate phenotype compared to an analogous two-step method to identify GxE effect for a univariate phenotype. We applied the proposed approach to simultaneously analyze three lipids, LDL, HDL and Triglyceride with the frequency of alcohol consumption as environmental factor in the UK Biobank. The method identified two loci with an overall GxE effect on the vector of lipids, one of which was missed by the competing approaches. We provide an R package MPGE implementing the proposed approach which is available from CRAN: https://cran.r-project.org/web/packages/MPGE/index.html. Supplementary data are available at Bioinformatics online.

A functional U-statistic method for association analysis of sequencing data.

Although sequencing studies hold great promise for uncovering novel variants predisposing to human diseases, the high dimensionality of the sequencing data brings tremendous challenges to data analysis. Moreover, for many complex diseases (e.g., psychiatric disorders) multiple related phenotypes are collected. These phenotypes can be different measurements of an underlying disease, or measurements characterizing multiple related diseases for studying common genetic mechanism. Although jointly analyzing these phenotypes could potentially increase the power of identifying disease-associated genes, the different types of phenotypes pose challenges for association analysis. To address these challenges, we propose a nonparametric method, functional U-statistic method (FU), for multivariate analysis of sequencing data. It first constructs smooth functions from individuals' sequencing data, and then tests the association of these functions with multiple phenotypes by using a U-statistic. The method provides a general framework for analyzing various types of phenotypes (e.g., binary and continuous phenotypes) with unknown distributions. Fitting the genetic variants within a gene using a smoothing function also allows us to capture complexities of gene structure (e.g., linkage disequilibrium, LD), which could potentially increase the power of association analysis. Through simulations, we compared our method to the multivariate outcome score test (MOST), and found that our test attained better performance than MOST. In a real data application, we apply our method to the sequencing data from Minnesota Twin Study (MTS) and found potential associations of several nicotine receptor subunit (CHRN) genes, including CHRNB3, associated with nicotine dependence and/or alcohol dependence.

Analysis of multiple related phenotypes in genome-wide association studies.

Most genome-wide association studies (GWAS) have been conducted by focusing on one phenotype of interest for identifying genetic variants associated with common complex phenotypes. However, despite many successful results from GWAS, only a small number of genetic variants tend to be identified and replicated given a very stringent genome-wide significance criterion, and explain only a small fraction of phenotype heritability. In order to improve power by using more information from data, we propose an alternative multivariate approach, which considers multiple related phenotypes simultaneously. We demonstrate through computer simulation that the multivariate approach can improve power for detecting disease-predisposing genetic variants and pleiotropic variants that have simultaneous effects on multiple related phenotypes. We apply the multivariate approach to a GWA dataset of 8,842 Korean individuals genotyped for 327,872 SNPs, and detect novel genetic variants associated with metabolic syndrome related phenotypes. Considering several related phenotype simultaneously, the multivariate approach provides not only more powerful results than the conventional univariate approach but also clue to identify pleiotropic genes that are important to the pathogenesis of many related complex phenotypes.

A method for analyzing multiple continuous phenotypes in rare variant association studies allowing for flexible correlations in variant effects.

For region-based sequencing data, power to detect genetic associations can be improved through analysis of multiple related phenotypes. With this motivation, we propose a novel test to detect association simultaneously between a set of rare variants, such as those obtained by sequencing in a small genomic region, and multiple continuous phenotypes. We allow arbitrary correlations among the phenotypes and build on a linear mixed model by assuming the effects of the variants follow a multivariate normal distribution with a zero mean and a specific covariance matrix structure. In order to account for the unknown correlation parameter in the covariance matrix of the variant effects, a data-adaptive variance component test based on score-type statistics is derived. As our approach can calculate the P-value analytically, the proposed test procedure is computationally efficient. Broad simulations and an application to the UK10K project show that our proposed multivariate test is generally more powerful than univariate tests, especially when there are pleiotropic effects or highly correlated phenotypes.

A multivariate genome-wide association analysis of 10 LDL subfractions, and their response to statin treatment, in 1868 Caucasians.

We conducted a genome-wide association analysis of 7 subfractions of low density lipoproteins (LDLs) and 3 subfractions of intermediate density lipoproteins (IDLs) measured by gradient gel electrophoresis, and their response to statin treatment, in 1868 individuals of European ancestry from the Pharmacogenomics and Risk of Cardiovascular Disease study. Our analyses identified four previously-implicated loci (SORT1, APOE, LPA, and CETP) as containing variants that are very strongly associated with lipoprotein subfractions (log10Bayes Factor > 15). Subsequent conditional analyses suggest that three of these (APOE, LPA and CETP) likely harbor multiple independently associated SNPs. Further, while different variants typically showed different characteristic patterns of association with combinations of subfractions, the two SNPs in CETP show strikingly similar patterns - both in our original data and in a replication cohort - consistent with a common underlying molecular mechanism. Notably, the CETP variants are very strongly associated with LDL subfractions, despite showing no association with total LDLs in our study, illustrating the potential value of the more detailed phenotypic measurements. In contrast with these strong subfraction associations, genetic association analysis of subfraction response to statins showed much weaker signals (none exceeding log10Bayes Factor of 6). However, two SNPs (in APOE and LPA) previously-reported to be associated with LDL statin response do show some modest evidence for association in our data, and the subfraction response proles at the LPA SNP are consistent with the LPA association, with response likely being due primarily to resistance of Lp(a) particles to statin therapy. An additional important feature of our analysis is that, unlike most previous analyses of multiple related phenotypes, we analyzed the subfractions jointly, rather than one at a time. Comparisons of our multivariate analyses with standard univariate analyses demonstrate that multivariate analyses can substantially increase power to detect associations. Software implementing our multivariate analysis methods is available at http://stephenslab.uchicago.edu/software.html.

CAPE: An R Package for Combined Analysis of Pleiotropy and Epistasis

Contemporary genetic studies are revealing the genetic complexity of many traits in humans and model organisms. Two hallmarks of this complexity are epistasis, meaning gene-gene interaction, and pleiotropy, in which one gene affects multiple phenotypes. Understanding the genetic architecture of complex traits requires addressing these phenomena, but interpreting the biological significance of epistasis and pleiotropy is often difficult. While epistasis reveals dependencies between genetic variants, it is often unclear how the activity of one variant is specifically modifying the other. Epistasis found in one phenotypic context may disappear in another context, rendering the genetic interaction ambiguous. Pleiotropy can suggest either redundant phenotype measures or gene variants that affect multiple biological processes. Here we present an R package, R/cape, which addresses these interpretation ambiguities by implementing a novel method to generate predictive and interpretable genetic networks that influence quantitative phenotypes. R/cape integrates information from multiple related phenotypes to constrain models of epistasis, thereby enhancing the detection of interactions that simultaneously describe all phenotypes. The networks inferred by R/cape are readily interpretable in terms of directed influences that indicate suppressive and enhancing effects of individual genetic variants on other variants, which in turn account for the variance in quantitative traits. We demonstrate the utility of R/cape by analyzing a mouse backcross, thereby discovering novel epistatic interactions influencing phenotypes related to obesity and diabetes. R/cape is an easy-to-use, platform-independent R package and can be applied to data from both genetic screens and a variety of segregating populations including backcrosses, intercrosses, and natural populations. The package is freely available under the GPL-3 license at http://cran.r-project.org/web/packages/cape.

Yunis-Varón Syndrome Is Caused by Mutations in FIG4, Encoding a Phosphoinositide Phosphatase

Yunis-Varón syndrome (YVS) is an autosomal-recessive disorder with cleidocranial dysplasia, digital anomalies, and severe neurological involvement. Enlarged vacuoles are found in neurons, muscle, and cartilage. By whole-exome sequencing, we identified frameshift and missense mutations of FIG4 in affected individuals from three unrelated families. FIG4 encodes a phosphoinositide phosphatase required for regulation of PI(3,5)P(2) levels, and thus endosomal trafficking and autophagy. In a functional assay, both missense substitutions failed to correct the vacuolar phenotype of Fig4-null mouse fibroblasts. Homozygous Fig4-null mice exhibit features of YVS, including neurodegeneration and enlarged vacuoles in neurons. We demonstrate that Fig4-null mice also have small skeletons with reduced trabecular bone volume and cortical thickness and that cultured osteoblasts accumulate large vacuoles. Our findings demonstrate that homozygosity or compound heterozygosity for null mutations of FIG4 is responsible for YVS, the most severe known human phenotype caused by defective phosphoinositide metabolism. In contrast, in Charcot-Marie-Tooth disease type 4J (also caused by FIG4 mutations), one of the FIG4 alleles is hypomorphic and disease is limited to the peripheral nervous system. This genotype-phenotype correlation demonstrates that absence of FIG4 activity leads to central nervous system dysfunction and extensive skeletal anomalies. Our results describe a role for PI(3,5)P(2) signaling in skeletal development and maintenance.

Increasing Power

Many of the phenotypic risk markers that are known, or suspected, to be involved in disease are interrelated. Given that many genes have pleiotropic effects, where the genetic variant is associated with multiple outcome variables, it is likely that multiple related phenotypes have a common genetic basis. At present, majority of genome-wide association studies analyze the effect of single nucleotide polymorphisms (SNPs) on individual phenotypes in univariate analyses. Where multiple phenotypes are examined, applying a traditional Bonferroni correction for multiple testing may be overly stringent, and does not take into account the common genetic origin of the correlated phenotypes. Inouye et al hypothesized that by analyzing multiple related phenotypes together in a multivariate genetic analysis, rather than individually in a typical univariate analysis, they would have increased power to detect genetic associations with disease-related phenotypes. The discovery effort focused on 2 separate population-based cohorts from Finland with genome-wide SNP genotypes, and extensive serum metabolite profiles. Using data from 6,600 individuals, the authors first identified 11 metabolic networks from the 130 metabolites profiled. This metabolic panel included metabolites …

Analysis of multiple phenotypes

The complex etiology of common diseases like cardiovascular disease, diabetes, hypertension, and rheumatoid arthritis has led investigators to focus on the genetics of correlated phenotypes and risk factors. Joint analysis of multiple disease-related phenotypes may reveal genes of pleiotropic effect and increase analytical power, but at the cost of increased analytical and computational complexity. All three data sets provided for analysis at the Genetic Analysis Workshop 16 offered multiple quantitative measures of phenotypes related to underlying disease processes as well as discrete measures of affection status. Participants in Group 6 addressed the challenges and possibilities of association analysis of these data sets on multiple levels, including phenotype definition and data reduction, multivariate approaches to gene discovery, analysis of causality and data structure, and development of predictive models. These approaches included combinations of continuous and discrete phenotypes, use of repeated measures in longitudinal data, and models that included multiple phenotypic measures and multiple single-nucleotide polymorphism variants. Most research teams regarded the use of multiple related phenotypes as a tool for increasing analytical power, as well as for clarifying the underlying biology of complex diseases.

Hierarchical modeling in association studies of multiple phenotypes

The genetic study of disease-associated phenotypes has become common because such phenotypes are often easier to measure and in many cases are under greater genetic control than the complex disease itself. Some disease-associated phenotypes are rare, however, making it difficult to evaluate their effects due to small informative sample sizes. In addition, analyzing numerous phenotypes introduces the issue of multiple comparisons. To address these issues, we have developed a hierarchical model (HM) for multiple phenotypes that provides more accurate effect estimates with a lower false-positive rate. We evaluated the validity and power of HM in association studies of multiple phenotypes using randomly selected cases and controls from the simulated data set in the Genetic Analysis Workshop 14. In particular, we first analyzed the association between each of the 12 subclinical phenotypes and single-nucleotide polymorphisms within the known causal loci using a conventional logistic regression model (LRM). Then we added a second-stage model by regressing all of the logistic coefficients of the phenotypes obtained from LRM on a Z matrix that incorporates the clinical correlation of the phenotypes. Specially, the 12 phenotypes were grouped into 3 clusters: 1) communally shared emotions; 2) behavioral related; and 3) anxiety related. A semi-Bayes HM effect estimate for each phenotype was calculated and compared with those from LRM. We observed that using HM to evaluate the association between SNPs and multiple related phenotypes slightly increased power for detecting the true associations and also led to fewer false-positive results.