RESUMEN
Multivariate analysis is becoming central in studies investigating high-throughput molecular data, yet, some important features of these data are seldom explored. Here, we present MANOCCA (Multivariate Analysis of Conditional CovAriance), a powerful method to test for the effect of a predictor on the covariance matrix of a multivariate outcome. The proposed test is by construction orthogonal to tests based on the mean and variance and is able to capture effects that are missed by both approaches. We first compare the performances of MANOCCA with existing correlation-based methods and show that MANOCCA is the only test correctly calibrated in simulation mimicking omics data. We then investigate the impact of reducing the dimensionality of the data using principal component analysis when the sample size is smaller than the number of pairwise covariance terms analysed. We show that, in many realistic scenarios, the maximum power can be achieved with a limited number of components. Finally, we apply MANOCCA to 1000 healthy individuals from the Milieu Interieur cohort, to assess the effect of health, lifestyle and genetic factors on the covariance of two sets of phenotypes, blood biomarkers and flow cytometry-based immune phenotypes. Our analyses identify significant associations between multiple factors and the covariance of both omics data.
Asunto(s)
Análisis de Componente Principal , Humanos , Análisis Multivariante , Biología Computacional/métodos , Fenotipo , Algoritmos , Genómica/métodos , Biomarcadores/sangre , Simulación por ComputadorRESUMEN
Primary aldosteronism affects up to 10% of hypertensive patients and is responsible for treatment resistance and increased cardiovascular risk. Here we perform a genome-wide association study in a discovery cohort of 562 cases and 950 controls and identify three main loci on chromosomes 1, 13 and X; associations on chromosome 1 and 13 are replicated in a second cohort and confirmed by a meta-analysis involving 1162 cases and 3296 controls. The association on chromosome 13 is specific to men and stronger in bilateral adrenal hyperplasia than aldosterone producing adenoma. Candidate genes located within the two loci, CASZ1 and RXFP2, are expressed in human and mouse adrenals in different cell clusters. Their overexpression in adrenocortical cells suppresses mineralocorticoid output under basal and stimulated conditions, without affecting cortisol biosynthesis. Our study identifies the first risk loci for primary aldosteronism and highlights new mechanisms for the development of aldosterone excess.
Asunto(s)
Neoplasias de la Corteza Suprarrenal , Adenoma Corticosuprarrenal , Hiperaldosteronismo , Neoplasias de la Corteza Suprarrenal/genética , Neoplasias de la Corteza Suprarrenal/cirugía , Adrenalectomía , Adenoma Corticosuprarrenal/genética , Adenoma Corticosuprarrenal/cirugía , Aldosterona , Animales , Proteínas de Unión al ADN/genética , Estudio de Asociación del Genoma Completo , Humanos , Hiperaldosteronismo/genética , Masculino , Ratones , Factores de Transcripción/genéticaRESUMEN
Genome-Wide Association Studies (GWAS) explain only a small fraction of heritability for most complex human phenotypes. Genomic heritability estimates the variance explained by the SNPs on the whole genome using mixed models and accounts for the many small contributions of SNPs in the explanation of a phenotype. This paper approaches heritability from a machine learning perspective, and examines the close link between mixed models and ridge regression. Our contribution is two-fold. First, we propose estimating genomic heritability using a predictive approach via ridge regression and Generalized Cross Validation (GCV). We show that this is consistent with classical mixed model based estimation. Second, we derive simple formulae that express prediction accuracy as a function of the ratio n p , where n is the population size and p the total number of SNPs. These formulae clearly show that a high heritability does not imply an accurate prediction when p > n. Both the estimation of heritability via GCV and the prediction accuracy formulae are validated using simulated data and real data from UK Biobank.