Multivariate Association Test Research Articles

A copula-based set-variant association test for bivariate continuous, binary or mixed phenotypes.

In genome wide association studies (GWAS), researchers are often dealing with dichotomous and non-normally distributed traits, or a mixture of discrete-continuous traits. However, most of the current region-based methods rely on multivariate linear mixed models (mvLMMs) and assume a multivariate normal distribution for the phenotypes of interest. Hence, these methods are not applicable to disease or non-normally distributed traits. Therefore, there is a need to develop unified and flexible methods to study association between a set of (possibly rare) genetic variants and non-normal multivariate phenotypes. Copulas are multivariate distribution functions with uniform margins on the [0, 1] interval and they provide suitable models to deal with non-normality of errors in multivariate association studies. We propose a novel unified and flexible copula-based multivariate association test (CBMAT) for discovering association between a genetic region and a bivariate continuous, binary or mixed phenotype. We also derive a data-driven analytic p-value procedure of the proposed region-based score-type test. Through simulation studies, we demonstrate that CBMAT has well controlled type I error rates and higher power to detect associations compared with other existing methods, for discrete and non-normally distributed traits. At last, we apply CBMAT to detect the association between two genes located on chromosome 11 and several lipid levels measured on 1477 subjects from the ASLPAC study.

Identifying pleiotropic genes for major psychiatric disorders with GWAS summary statistics using multivariate adaptive association tests

BackgroundGenome wide association studies (GWAS) have discovered a few of single nucleotide polymorphisms (SNPs) related to major psychiatric disorders. However, it is not completely clear which genes play a pleiotropic role in multiple disorders. The study aimed to identify the pleiotropic genes across five psychiatric disorders using multivariate adaptive association tests. MethodsSummary statistics of five psychiatric disorders were downloaded from Psychiatric Genomics Consortium. We applied linkage disequilibrium score regression (LDSC) to estimate genetic correlation and conducted tissue and cell type specificity analyses based on Multi-marker Analysis of GenoMic Annotation (MAGMA). Then, we identified the pleiotropic genes using MTaSPUsSet and aSPUs tests. We ultimately performed the functional analysis for pleiotropic genes. ResultsWe confirmed the significant genetic correlation and brain tissue and neuron specificity among five disorders. 100 pleiotropic genes were detected to be significantly associated with five psychiatric disorders, of which 55 were novel genes. These genes were functionally enriched in neuron differentiation and synaptic transmission. LimitationsThe effect direction of pleiotropic genes couldn't be distinguished due to without individual-level data. ConclusionWe identified pleiotropic genes using multivariate adaptive association tests and explored their biological function. The findings may provide novel insight into the development and implementation of prevention and treatment as well as targeted drug discovery in practice.

Abstract P3-13-02: A genome-wide association study of mammographic texture variation

Abstract Background: Breast parenchymal texture features capture the patterns of texture variation on a mammogram. One such texture feature is called V, which reflects the gray scale variation on the image. Recent studies have shown that texture features, including V, are associated with breast cancer risk, independent of mammographic density (MD). However, our knowledge on the genetic basis of texture features and the shared genetic components between these features, MD, and breast cancer is limited. Methods: We conducted a genome-wide association study (GWAS) of breast texture variation, including four different assessments of V, in 7,040 European-ancestry women within the Nurses’ Health Studies and Mayo Mammography Health Study cohorts. Genotype data were imputed to the 1000 Genomes Phase 3 version 5 reference panel. The four V assessments, differ by amount of breast edge erosion and image resolution, were generated from digitized film mammograms using fully automated methods. We used linear regression to test the single-nucleotide polymorphism (SNP)-phenotype associations adjusting for age, body mass index (BMI), MD phenotypes (percent density, dense area, and nondense area), and the top four principal components. Multivariate phenotype association tests were performed to estimate summary test statistics of the four V assessments. We compared our results to GWAS of breast cancer and MD phenotypes, and identified shared susceptibility loci. We calculated genetic correlations of V measures with MD phenotypes, breast cancer, overall and by estrogen receptor status, BMI, childhood body fatness, age at menarche, and age at natural menopause. Aggregate tests were performed to assess the relationship between V measures and the above traits using genome-wide significant SNPs for these traits. Results: We identified three novel genome-wide significant loci that are associated with V measures: rs79670367 (8q24.22) in LINC01591, rs113174754 (12q22) near PGAM1P5, and rs138141444 (6q24.1) in ECT2L. 12q22 is a known locus for both MD (nondense area) and breast cancer. 8q24.22 and 6q24.1 have not previously been associated with MD or breast cancer. We identified four additional loci that are associated with MD or breast cancer, or both and are significantly associated with the V measures at the Bonferroni-corrected thresholds accounting for the number of MD or breast cancer SNPs tested: rs335189 (5q23.2) in PRDM6, rs13256025 (8p21.2) in EBF2, rs11836164 (12p12.1) in LOC105369705, and rs17817449 (16q12.2) in FTO, the directions of association are consistent with MD or breast cancer for all four loci. We observed significant genetic correlations between V measures and dense area (rg = 0.79, p = 6.0 × 10-5), percent density (rg = 0.73, p = 1.0 × 10-4), and BMI (rg = -0.36, p = 3.9 × 10-7). Additional significant relationships were observed for nondense area (z = -4.14, p = 3.5 × 10-5) and childhood body fatness (z = -4.91, p = 9.1 × 10-7) from the aggregate tests. Conclusions: This is the first GWAS of breast parenchymal texture features. These findings provide new insights into the genetic basis of texture variation and their associations with MD, breast cancer risk, and other breast cancer risk factors. Citation Format: Yuxi Liu, Xia Jiang, Constance Turman, Celine M Vachon, Rulla Tamimi, Peter Kraft. A genome-wide association study of mammographic texture variation [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P3-13-02.

Genetic Variants and Their Associations to Type 2 Diabetes Mellitus Complications in the United Arab Emirates.

AimType 2 Diabetes Mellitus (T2DM) is associated with microvascular complications, including diabetic retinopathy (DR), diabetic nephropathy (DNp), and diabetic peripheral neuropathy (DPN). In this study, we investigated genetic variations and Single Nucleotide Polymorphisms (SNPs) associated with DR, DNp, DPN and their combinations among T2DM patients of Arab origin from the United Arab Emirates, to establish the role of genes in the progression of microvascular diabetes complications.MethodsA total of 158 Emirati patients with T2DM were recruited in this study. The study population was divided into 8 groups based on the presence of single, dual, or all three complications. SNPs were selected for association analyses through a search of publicly available databases, specifically genome-wide association study (GWAS) catalog, infinome genome interpretation platform, and GWAS Central database. A multivariate logistic regression analysis and association test were performed to evaluate the association between 83 SNPs and DR, DNp, DPN, and their combinations.ResultsEighty-three SNPs were identified as being associated with T2DM and 18 SNPs had significant associations to one or more diabetes complications. The most strongly significant association for DR was rs3024997 SNP in the VEGFA gene. The top-ranked SNP for DPN was rs4496877 in the NOS3 gene. A trend towards association was detected at rs833068 and rs3024998 in the VEGFA gene with DR and rs743507 and rs1808593 in the NOS3 gene with DNp. For dual complications, the rs833061, rs833068 and rs3024997 in the VEGFA gene and the rs4149263 SNP in the ABCA1 gene were also with borderline association with DR/DNp and DPN/DNp, respectively. Diabetic with all of the complications was significantly associated with rs2230806 in the ABCA1 gene. In addition, the highly associated SNPs rs3024997 of the VEGFA gene and rs4496877 of the NOS3 gene were linked to DR and DPN after adjusting for the effects of other associated markers, respectively.ConclusionsThe present study reports associations of different genetic polymorphisms with microvascular complications and their combinations in Emirati T2DM patients, reporting new associations, and corroborating previous findings. Of interest is that some SNPs/genes were only present if multiple comorbidities were present and not associated with any single complication.

Multivariate genome wide association and network analysis of subcortical imaging phenotypes in Alzheimer\u2019s disease

BackgroundGenome-wide association studies (GWAS) have identified many individual genes associated with brain imaging quantitative traits (QTs) in Alzheimer’s disease (AD). However single marker level association discovery may not be able to address the underlying biological interactions with disease mechanism.ResultsIn this paper, we used the MGAS (Multivariate Gene-based Association test by extended Simes procedure) tool to perform multivariate GWAS on eight AD-relevant subcortical imaging measures. We conducted multiple iPINBPA (integrative Protein-Interaction-Network-Based Pathway Analysis) network analyses on MGAS findings using protein-protein interaction (PPI) data, and identified five Consensus Modules (CMs) from the PPI network. Functional annotation and network analysis were performed on the identified CMs. The MGAS yielded significant hits within APOE, TOMM40 and APOC1 genes, which were known AD risk factors, as well as a few new genes such as LAMA1, XYLB, HSD17B7P2, and NPEPL1. The identified five CMs were enriched by biological processes related to disorders such as Alzheimer’s disease, Legionellosis, Pertussis, and Serotonergic synapse.ConclusionsThe statistical power of coupling MGAS with iPINBPA was higher than traditional GWAS method, and yielded new findings that were missed by GWAS. This study provides novel insights into the molecular mechanism of Alzheimer’s Disease and will be of value to novel gene discovery and functional genomic studies.

The human gut microbiota in people with amyotrophic lateral sclerosis

Objective To characterize the gut microbiota in people with amyotrophic lateral sclerosis (ALS) relative to controls and to test the hypothesis that butyrate-producing bacteria are less abundant in the gastrointestinal tracts of people with ALS (PALS). Methods: We conducted a case–control study at Massachusetts General Hospital to compare the gut microbiota in people with ALS to that in controls. Metagenomic shotgun sequencing was performed on DNA extracted from stool samples of 66 people with ALS (PALS), 61 healthy controls (HC), and 12 neurodegenerative controls (NDC). Taxonomic metagenomic profiles were analyzed for shifts in the microbial community structure between the comparator groups using per-feature univariate and multivariate association tests. Results: The relative abundance of the dominant butyrate-producing bacteria Eubacterium rectale and Roseburia intestinalis was significantly lower in ALS patients compared to HC. Adjustment for age, sex, and constipation did not materially change the results. The total abundance of 8 dominant species capable of producing butyrate was also significantly lower in ALS compared to HC (p < 0.001). Conclusions: The levels of several butyrate-producing bacteria, which are important for gut integrity and regulation of inflammation, were lower in people with ALS compared to controls. These findings lend support to the inference that the gut microbiota could be a risk factor for ALS. Further investigations are warranted, preferably earlier in the disease with corresponding dietary collection and a longitudinal design.

MF-TOWmuT: Testing an optimally weighted combination of common and rare variants with multiple traits using family data.

With rapid advancements of sequencing technologies and accumulations of electronic health records, a large number of genetic variants and multiple correlated human complex traits have become available in many genetic association studies. Thus, it becomes necessary and important to develop new methods that can jointly analyze the association between multiple genetic variants and multiple traits. Compared with methods that only use a single marker or trait, the joint analysis of multiple genetic variants and multiple traits is more powerful since such an analysis can fully incorporate the correlation structure of genetic variants and/or traits and their mutual dependence patterns. However, most of existing methods that simultaneously analyze multiple genetic variants and multiple traits are only applicable to unrelated samples. We develop a new method called MF-TOWmuT to detect association of multiple phenotypes and multiple genetic variants in a genomic region with family samples. MF-TOWmuT is based on an optimally weighted combination of variants. Our method can be applied to both rare and common variants and both qualitative and quantitative traits. Our simulation results show that (1) the type I error of MF-TOWmuT is preserved; (2) MF-TOWmuT outperforms two existing methods such as Multiple Family-based Quasi-Likelihood Score Test and Multivariate Family-based Rare Variant Association Testin terms of power. We also illustrate the usefulness of MF-TOWmuT by analyzing genotypic and phenotipic data from the Genetics of Kidneys in Diabetes study. R program is available at https://github.com/gaochengPRC/MF-TOWmuT.

Comparison of methods for multivariate gene-based association tests for complex diseases using common variants

Complex diseases are usually associated with multiple correlated phenotypes, and the analysis of composite scores or disease status may not fully capture the complexity (or multidimensionality). Joint analysis of multiple disease-related phenotypes in genetic tests could potentially increase power to detect association of a disease with common SNPs (or genes). Gene-based tests are designed to identify genes containing multiple risk variants that individually are weakly associated with a univariate trait. We combined three multivariate association tests (O’Brien method, TATES, and MultiPhen) with two gene-based association tests (GATES and VEGAS) and compared performance (type I error and power) of six multivariate gene-based methods using simulated data. Data (n = 2000) for genetic sequence and correlated phenotypes were simulated by varying causal variant proportions and phenotype correlations for various scenarios. These simulations showed that two multivariate association tests (TATES and MultiPhen, but not O’Brien) paired with VEGAS have inflated type I error in all scenarios, while the three multivariate association tests paired with GATES have correct type I error. MultiPhen paired with GATES has higher power than competing methods if the correlations among phenotypes are low (r < 0.57). We applied these gene-based association methods to a GWAS dataset from the Alzheimer’s Disease Genetics Consortium containing three neuropathological traits related to Alzheimer disease (neuritic plaque, neurofibrillary tangles, and cerebral amyloid angiopathy) measured in 3500 autopsied brains. Gene-level significant evidence (P < 2.7 × 10−6) was identified in a region containing three contiguous genes (TRAPPC12, TRAPPC12-AS1, ADI1) using O’Brien and VEGAS. Gene-wide significant associations were not observed in univariate gene-based tests.

Multivariate association test for rare variants controlling for cryptic and family relatedness

AbstractIn genetic studies of complex diseases, multiple measures of related phenotypes are often collected. Jointly analyzing these phenotypes may improve statistical power to detect sets of rare variants affecting multiple traits. In this work, we consider association testing between a set of rare variants and multiple phenotypes in family‐based designs. We use a mixed linear model to express the correlations among the phenotypes and between related individuals. Given the many sources of correlations in this situation, deriving an appropriate test statistic is not straightforward. We derive a vector of score statistics, whose joint distribution is approximated using a copula. This allows us to have closed‐form expressions for the p‐values of several test statistics. A comprehensive simulation study and an application to Genetic Analysis Workshop 18 data highlight the gains associated with joint testing over univariate approaches, especially in the presence of pleiotropy or highly correlated phenotypes. The Canadian Journal of Statistics 47: 90–107; 2019 © 2018 Statistical Society of Canada

Loss-of-function mutation in inositol monophosphatase 1 (IMPA1) results in abnormal synchrony in resting-state EEG

BackgroundDysregulation of the inositol cycle is implicated in a wide variety of human diseases, including developmental defects and neurological diseases. A homozygous frameshift mutation in IMPA1, coding for the enzyme inositol monophosphatase 1 (IMPase), has recently been associated with severe intellectual disability (ID) in a geographically isolated consanguineous family in Northeastern Brazil (Figueredo et al., 2016). However, the neurophysiologic mechanisms that mediate the IMPA1 mutation and associated ID phenotype have not been characterized. To this end, resting EEG (eyes-open and eyes-closed) was collected from the Figueredo et al. pedigree. Quantitative EEG measures, including mean power, dominant frequency and dominant frequency variability, were investigated for allelic associations using multivariate family-based association test using generalized estimating equations.ResultsWe found that the IMPA1 mutation was associated with relative decreases in frontal theta band power as well as altered alpha-band variability with no regional specificity during the eyes-open condition. For the eyes-closed condition, there was altered dominant theta frequency variability in the central and parietal regions.ConclusionsThese findings represent the first human in vivo phenotypic assessment of brain function disturbances associated with a loss-of-function IMPA1 mutation, and thus an important first step towards an understanding the pathophysiologic mechanisms of intellectual disability associated with the mutation that affects this critical metabolic pathway.

Identifying pleiotropic genes in genome-wide association studies from related subjects using the linear mixed model and Fisher combination function

BackgroundA multivariate genome-wide association test is proposed for analyzing data on multivariate quantitative phenotypes collected from related subjects. The proposed method is a two-step approach. The first step models the association between the genotype and marginal phenotype using a linear mixed model. The second step uses the correlation between residuals of the linear mixed model to estimate the null distribution of the Fisher combination test statistic.ResultsThe simulation results show that the proposed method controls the type I error rate and is more powerful than the marginal tests across different population structures (admixed or non-admixed) and relatedness (related or independent). The statistical analysis on the database of the Study of Addiction: Genetics and Environment (SAGE) demonstrates that applying the multivariate association test may facilitate identification of the pleiotropic genes contributing to the risk for alcohol dependence commonly expressed by four correlated phenotypes.ConclusionsThis study proposes a multivariate method for identifying pleiotropic genes while adjusting for cryptic relatedness and population structure between subjects. The two-step approach is not only powerful but also computationally efficient even when the number of subjects and the number of phenotypes are both very large.

Genome-Wide Association Study of Psychosis Proneness in the Finnish Population.

The current study examined quantitative measures of psychosis proneness in a nonpsychotic population, in order to elucidate their underlying genetic architecture and to observe if there is any commonality to that already detected in the studies of individuals with overt psychotic conditions, such as schizophrenia and bipolar disorder. Heritability, univariate and multivariate genome-wide association (GWAs) tests, including a series of comprehensive gene-based association analyses, were developed in 4269 nonpsychotic persons participating in the Northern Finland Birth Cohort 1966 study with information on the following psychometric measures: Hypomanic Personality, Perceptual Aberration, Physical and Social Anhedonia (also known as Chapman's Schizotypia scales), and Schizoidia scale. Genome-wide genetic data was available for ~9.84 million SNPs. Heritability estimates ranged from 16% to 27%. Phenotypic, genetic and environmental correlations ranged from 0.04-0.43, 0.25-0.73, and 0.12-0.43, respectively. Univariate GWAs tests revealed an intronic SNP (rs12449097) at the TMC7 gene (16p12.3) that significantly associated (P = 3.485 × 10-8) with the hypomanic scale. Bivariate GWAs tests including the hypomanic and physical anhedonia scales suggested a further borderline significant SNP (rs188320715; P-value = 5.261 × 10-8, ~572 kb downstream the ARID1B gene at 6q25.3). Gene-based tests highlighted 20 additional genes of which 5 had previously been associated to schizophrenia and/or bipolar disorder: CSMD1, CCDC141, SLC1A2, CACNA1C, and SNAP25. Altogether the findings explained from 3.7% to 14.1% of the corresponding trait heritability. In conclusion, this study provides preliminary genomic evidence suggesting that qualitatively similar biological factors may underlie different psychosis proneness measures, some of which could further predispose to schizophrenia and bipolar disorder.

Identifying Pleiotropic Genes in Genome-Wide Association Studies for Multivariate Phenotypes with Mixed Measurement Scales.

We propose a multivariate genome-wide association test for mixed continuous, binary, and ordinal phenotypes. A latent response model is used to estimate the correlation between phenotypes with different measurement scales so that the empirical distribution of the Fisher’s combination statistic under the null hypothesis is estimated efficiently. The simulation study shows that our proposed correlation estimation methods have high levels of accuracy. More importantly, our approach conservatively estimates the variance of the test statistic so that the type I error rate is controlled. The simulation also shows that the proposed test maintains the power at the level very close to that of the ideal analysis based on known latent phenotypes while controlling the type I error. In contrast, conventional approaches–dichotomizing all observed phenotypes or treating them as continuous variables–could either reduce the power or employ a linear regression model unfit for the data. Furthermore, the statistical analysis on the database of the Study of Addiction: Genetics and Environment (SAGE) demonstrates that conducting a multivariate test on multiple phenotypes can increase the power of identifying markers that may not be, otherwise, chosen using marginal tests. The proposed method also offers a new approach to analyzing the Fagerström Test for Nicotine Dependence as multivariate phenotypes in genome-wide association studies.

A small-sample multivariate kernel machine test for microbiome association studies

High-throughput sequencing technologies have enabled large-scale studies of the role of the human microbiome in health conditions and diseases. Microbial community level association test, as a critical step to establish the connection between overall microbiome composition and an outcome of interest, has now been routinely performed in many studies. However, current microbiome association tests all focus on a single outcome. It has become increasingly common for a microbiome study to collect multiple, possibly related, outcomes to maximize the power of discovery. As these outcomes may share common mechanisms, jointly analyzing these outcomes can amplify the association signal and improve statistical power to detect potential associations. We propose the multivariate microbiome regression-based kernel association test (MMiRKAT) for testing association between multiple continuous outcomes and overall microbiome composition, where the kernel used in MMiRKAT is based on Bray-Curtis or UniFrac distance. MMiRKAT directly regresses all outcomes on the microbiome profiles via a semiparametric kernel machine regression framework, which allows for covariate adjustment and evaluates the association via a variance-component score test. Because most of the current microbiome studies have small sample sizes, a novel small-sample correction procedure is implemented in MMiRKAT to correct for the conservativeness of the association test when the sample size is small or moderate. The proposed method is assessed via simulation studies and an application to a real data set examining the association between host gene expression and mucosal microbiome composition. We demonstrate that MMiRKAT is more powerful than large sample based multivariate kernel association test, while controlling the type I error. A free implementation of MMiRKAT in R language is available at http://research.fhcrc.org/wu/en.html.

Joint analysis of multiple blood pressure phenotypes in GAW19 data by using a multivariate rare-variant association test.

IntroductionLarge-scale sequencing studies often measure many related phenotypes in addition to the genetic variants. Joint analysis of multiple phenotypes in genetic association studies may increase power to detect disease-associated loci.MethodsWe apply a recently developed multivariate rare-variant association test to the Genetic Analysis Workshop 19 data in order to test associations between genetic variants and multiple blood pressure phenotypes simultaneously. We also compare this multivariate test with a widely used univariate test that analyzes phenotypes separately.ResultsThe multivariate test identified 2 genetic variants that have been previously reported as associated with hypertension or coronary artery disease. In addition, our region-based analyses also show that the multivariate test tends to give smaller p values than the univariate test.ConclusionsHence, the multivariate test has potential to improve test power, especially when multiple phenotypes are correlated.

Multivariate Gene-Based Association Test on Family Data in MGAS

In analyses of unrelated individuals, the program multivariate gene-based association test by extended Simes (MGAS), which facilitates multivariate gene-based association testing, was shown to have correct Type I error rate and superior statistical power compared to other multivariate gene-based approaches. Here we show, through simulation, that MGAS can also be applied to data including genetically related subjects (e.g., family data), by using p value information obtained in Plink or in generalized estimating equations (with the 'exchangeable' working correlation matrix), both of which account for the family structure on a univariate single nucleotide polymorphism-based level by applying a sandwich correction of standard errors. We show that when applied to family-data, MGAS has correct Type I error rate, and given the details of the simulation setup, adequate power. Application of MGAS to seven eye measurement phenotypes showed statistically significant association with two genes that were not discovered in previous univariate analyses of a composite score. We conclude that MGAS is a useful and convenient tool for multivariate gene-based genome-wide association analysis in both unrelated and related individuals.

Function of cancer associated genes revealed by modern univariate and multivariate association tests.

Copy number variation (CNV) plays a role in pathogenesis of many human diseases, especially cancer. Several whole genome CNV association studies have been performed for the purpose of identifying cancer associated CNVs. Here we undertook a novel approach to whole genome CNV analysis, with the goal being identification of associations between CNV of different genes (CNV-CNV) across 60 human cancer cell lines. We hypothesize that these associations point to the roles of the associated genes in cancer, and can be indicators of their position in gene networks of cancer-driving processes. Recent studies show that gene associations are often non-linear and non-monotone. In order to obtain a more complete picture of all CNV associations, we performed omnibus univariate analysis by utilizing dCov, MIC, and HHG association tests, which are capable of detecting any type of association, including non-monotone relationships. For comparison we used Spearman and Pearson association tests, which detect only linear or monotone relationships. Application of dCov, MIC and HHG tests resulted in identification of twice as many associations compared to those found by Spearman and Pearson alone. Interestingly, most of the new associations were detected by the HHG test. Next, we utilized dCov's and HHG's ability to perform multivariate analysis. We tested for association between genes of unknown function and known cancer-related pathways. Our results indicate that multivariate analysis is much more effective than univariate analysis for the purpose of ascribing biological roles to genes of unknown function. We conclude that a combination of multivariate and univariate omnibus association tests can reveal significant information about gene networks of disease-driving processes. These methods can be applied to any large gene or pathway dataset, allowing more comprehensive analysis of biological processes.

MGAS: a powerful tool for multivariate gene-based genome-wide association analysis.

Motivation: Standard genome-wide association studies, testing the association between one phenotype and a large number of single nucleotide polymorphisms (SNPs), are limited in two ways: (i) traits are often multivariate, and analysis of composite scores entails loss in statistical power and (ii) gene-based analyses may be preferred, e.g. to decrease the multiple testing problem.Results: Here we present a new method, multivariate gene-based association test by extended Simes procedure (MGAS), that allows gene-based testing of multivariate phenotypes in unrelated individuals. Through extensive simulation, we show that under most trait-generating genotype–phenotype models MGAS has superior statistical power to detect associated genes compared with gene-based analyses of univariate phenotypic composite scores (i.e. GATES, multiple regression), and multivariate analysis of variance (MANOVA). Re-analysis of metabolic data revealed 32 False Discovery Rate controlled genome-wide significant genes, and 12 regions harboring multiple genes; of these 44 regions, 30 were not reported in the original analysis.Conclusion: MGAS allows researchers to conduct their multivariate gene-based analyses efficiently, and without the loss of power that is often associated with an incorrectly specified genotype–phenotype models.Availability and implementation: MGAS is freely available in KGG v3.0 (http://statgenpro.psychiatry.hku.hk/limx/kgg/download.php). Access to the metabolic dataset can be requested at dbGaP (https://dbgap.ncbi.nlm.nih.gov/). The R-simulation code is available from http://ctglab.nl/people/sophie_van_der_sluis.Contact: mxli@hku.hkSupplementary information: Supplementary data are available at Bioinformatics online.

A Family‐Based Robust Multivariate Association Test Using Maximum Statistic

For characterizing the genetic mechanisms of complex diseases familial data with multiple correlated quantitative traits are usually collected in genetic studies. To analyze such data, various multivariate tests have been proposed to investigate the association between the underlying disease genes and the multiple traits. Although these multivariate association tests may have better power performance than the univariate association tests, they suffer from loss of testing power when the genetic models of the putative genes are misspecified. To address the problem, in this paper we aim to develop a family-based robust multivariate association test. We will first establish the optimal multivariate score tests for the recessive, additive, and dominant genetic models. Based on these optimal tests, a maximum-type robust multivariate association test is then obtained. Simulations are conducted to compare the power of our method with that of other existing multivariate methods. The results show that the robust multivariate test does manifest the robustness in power over all plausible genetic models. A practical data set is applied to demonstrate the applicability of our approach. The results suggest that the robust multivariate test is more powerful than the robust univariate test when dealing with multiple quantitative traits.

Applying genome-wide gene-based expression quantitative trait locus mapping to study population ancestry and pharmacogenetics.

BackgroundGene-based analysis has become popular in genomic research because of its appealing biological and statistical properties compared with those of a single-locus analysis. However, only a few, if any, studies have discussed a mapping of expression quantitative trait loci (eQTL) in a gene-based framework. Neither study has discussed ancestry-informative eQTL nor investigated their roles in pharmacogenetics by integrating single nucleotide polymorphism (SNP)-based eQTL (s-eQTL) and gene-based eQTL (g-eQTL).ResultsIn this g-eQTL mapping study, the transcript expression levels of genes (transcript-level genes; T-genes) were correlated with the SNPs of genes (sequence-level genes; S-genes) by using a method of gene-based partial least squares (PLS). Ancestry-informative transcripts were identified using a rank-score-based multivariate association test, and ancestry-informative eQTL were identified using Fisher’s exact test. Furthermore, key ancestry-predictive eQTL were selected in a flexible discriminant analysis. We analyzed SNPs and gene expression of 210 independent people of African-, Asian- and European-descent. We identified numerous cis- and trans-acting g-eQTL and s-eQTL for each population by using PLS. We observed ancestry information enriched in eQTL. Furthermore, we identified 2 ancestry-informative eQTL associated with adverse drug reactions and/or drug response. Rs1045642, located on MDR1, is an ancestry-informative eQTL (P = 2.13E-13, using Fisher’s exact test) associated with adverse drug reactions to amitriptyline and nortriptyline and drug responses to morphine. Rs20455, located in KIF6, is an ancestry-informative eQTL (P = 2.76E-23, using Fisher’s exact test) associated with the response to statin drugs (e.g., pravastatin and atorvastatin). The ancestry-informative eQTL of drug biotransformation genes were also observed; cross-population cis-acting expression regulators included SPG7, TAP2, SLC7A7, and CYP4F2. Finally, we also identified key ancestry-predictive eQTL and established classification models with promising training and testing accuracies in separating samples from close populations.ConclusionsIn summary, we developed a gene-based PLS procedure and a SAS macro for identifying g-eQTL and s-eQTL. We established data archives of eQTL for global populations. The program and data archives are accessible at http://www.stat.sinica.edu.tw/hsinchou/genetics/eQTL/HapMapII.htm. Finally, the results from our investigations regarding the interrelationship between eQTL, ancestry information, and pharmacodynamics provide rich resources for future eQTL studies and practical applications in population genetics and medical genetics.