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

FaceAge, a deep learning system to estimate biological age from face photographs to improve prognostication: a model development and validation study

BACKGROUND

As humans age at different rates, physical appearance can yield insights into biological age and physiological health more reliably than chronological age. In medicine, however, appearance is incorporated into medical judgements in a subjective and non-standardised way. In this study, we aimed to develop and validate FaceAge, a deep learning system to estimate biological age from easily obtainable and low-cost face photographs.

METHODS

FaceAge was trained on data from 58 851 presumed healthy individuals aged 60 years or older: 56 304 individuals from the IMDb-Wiki dataset (training) and 2547 from the UTKFace dataset (initial validation). Clinical utility was evaluated on data from 6196 patients with cancer diagnoses from two institutions in the Netherlands and the USA: the MAASTRO, Harvard Thoracic, and Harvard Palliative cohorts FaceAge estimates in these cancer cohorts were compared with a non-cancerous reference cohort of 535 individuals. To assess the prognostic relevance of FaceAge, we performed Kaplan-Meier survival analysis and Cox modelling, adjusting for several clinical covariates. We also assessed the performance of FaceAge in patients with metastatic cancer receiving palliative treatment at the end of life by incorporating FaceAge into clinical prediction models. To evaluate whether FaceAge has the potential to be a biomarker for molecular ageing, we performed a gene-based analysis to assess its association with senescence genes.

FINDINGS

FaceAge showed significant independent prognostic performance in various cancer types and stages. Looking older was correlated with worse overall survival (after adjusting for covariates per-decade hazard ratio [HR] 1·151, p=0·013 in a pan-cancer cohort of n=4906; 1·148, p=0·011 in a thoracic cohort of n=573; and 1·117, p=0·021 in a palliative cohort of n=717). We found that, on average, patients with cancer looked older than their chronological age (mean increase of 4·79 years with respect to non-cancerous reference cohort, p<0·0001). We found that FaceAge can improve physicians' survival predictions in patients with incurable cancer receiving palliative treatments (from area under the curve 0·74 [95% CI 0·70-0·78] to 0·8 [0·76-0·83]; p<0·0001), highlighting the clinical use of the algorithm to support end-of-life decision making. FaceAge was also significantly associated with molecular mechanisms of senescence through gene analysis, whereas age was not.

INTERPRETATION

Our results suggest that a deep learning model can estimate biological age from face photographs and thereby enhance survival prediction in patients with cancer. Further research, including validation in larger cohorts, is needed to verify these findings in patients with cancer and to establish whether the findings extend to patients with other diseases. Subject to further testing and validation, approaches such as FaceAge could be used to translate a patient's visual appearance into objective, quantitative, and clinically valuable measures.

FUNDING

US National Institutes of Health and EU European Research Council.

Modulatory Neurotransmitter Genotypes Shape Dynamic Functional Connectome Reconfigurations

Dynamic reconfigurations of the functional connectome across different connectivity states are highly heritable, predictive of cognitive abilities, and linked to mental health. Despite their established heritability, the specific polymorphisms that shape connectome dynamics are largely unknown. Given the widespread regulatory impact of modulatory neurotransmitters on functional connectivity, we comprehensively investigated a large set of single nucleotide polymorphisms (SNPs) of their receptors, metabolic enzymes, and transporters in 674 healthy adult subjects (347 females) from the Human Connectome Project. Preregistered modulatory neurotransmitter SNPs and dynamic connectome features entered a Stability Selection procedure with resampling. We found that specific subsets of these SNPs explain individual differences in temporal phenotypes of fMRI-derived connectome dynamics for which we previously established heritability. Specifically, noradrenergic polymorphisms explained Fractional Occupancy, i.e., the proportion of time spent in each connectome state, and cholinergic polymorphisms explained Transition Probability, i.e., the probability to transition between state pairs, respectively. This work identifies specific genetic effects on connectome dynamics via the regulatory impact of modulatory neurotransmitter systems. Our observations highlight the potential of dynamic connectome features as endophenotypes for neurotransmitter-focused precision psychiatry.

A Novel Statistical Method for Unmasking Sex-Specific Genomics Signatures in Complex Traits

Genotype-phenotype association studies have advanced our understanding of complex traits but often overlook sex-specific genetic signals. The growing awareness of sex-specific influences on human traits and diseases necessitates tailored statistical methodologies to dissect these genetic intricacies. Rare genetic variants play a significant role in disease development, often exhibiting stronger per-allele effects than common variants. In sex-dimorphic analysis, traits are viewed as having two sex-specific subsets rather than being uniformly defined. Existing methods for gene-based analysis of rare variants across multiple traits can identify shared genetic signals but cannot reveal the specific subsets from which significant signals originate. This means that when a significant signal is detected, it remains unclear whether it arises from the male samples, female samples, or both. To address this limitation, we propose SubsetRV, a new methodology capable of identifying genes associated with specific traits or diseases in males, females, or both. SubsetRV can also be applied to broader applications in multiple traits analysis. Simulation studies have demonstrated SubsetRV's reliability, and real data analysis on bipolar disorder and schizophrenia has revealed potential sex-specific genetic signals. SubsetRV offers a valuable tool for identifying sex-specific genetic candidates, aiding in understanding disease mechanisms. An R package for SubsetRV is available on GitHub. It can be accessed directly through this link: https://github.com/Mansouri-S/SubsetRV.

GbyE: an integrated tool for genome widely association study and genome selection based on genetic by environmental interaction

BACKGROUND

The growth and development of organism were dependent on the effect of genetic, environment, and their interaction. In recent decades, lots of candidate additive genetic markers and genes had been detected by using genome-widely association study (GWAS). However, restricted to computing power and practical tool, the interactive effect of markers and genes were not revealed clearly. And utilization of these interactive markers is difficult in the breeding and prediction, such as genome selection (GS).

RESULTS

Through the Power-FDR curve, the GbyE algorithm can detect more significant genetic loci at different levels of genetic correlation and heritability, especially at low heritability levels. The additive effect of GbyE exhibits high significance on certain chromosomes, while the interactive effect detects more significant sites on other chromosomes, which were not detected in the first two parts. In prediction accuracy testing, in most cases of heritability and genetic correlation, the majority of prediction accuracy of GbyE is significantly higher than that of the mean method, regardless of whether the rrBLUP model or BGLR model is used for statistics. The GbyE algorithm improves the prediction accuracy of the three Bayesian models BRR, BayesA, and BayesLASSO using information from genetic by environmental interaction (G × E) and increases the prediction accuracy by 9.4%, 9.1%, and 11%, respectively, relative to the Mean value method. The GbyE algorithm is significantly superior to the mean method in the absence of a single environment, regardless of the combination of heritability and genetic correlation, especially in the case of high genetic correlation and heritability.

CONCLUSIONS

Therefore, this study constructed a new genotype design model program (GbyE) for GWAS and GS using Kronecker product. which was able to clearly estimate the additive and interactive effects separately. The results showed that GbyE can provide higher statistical power for the GWAS and more prediction accuracy of the GS models. In addition, GbyE gives varying degrees of improvement of prediction accuracy in three Bayesian models (BRR, BayesA, and BayesCpi). Whatever the phenotype were missed in the single environment or multiple environments, the GbyE also makes better prediction for inference population set. This study helps us understand the interactive relationship between genomic and environment in the complex traits. The GbyE source code is available at the GitHub website ( https://github.com/liu-xinrui/GbyE ).

Decoding biological age from face photographs using deep learning

Because humans age at different rates, a person's physical appearance may yield insights into their biological age and physiological health more reliably than their chronological age. In medicine, however, appearance is incorporated into medical judgments in a subjective and non-standardized fashion. In this study, we developed and validated FaceAge, a deep learning system to estimate biological age from easily obtainable and low-cost face photographs. FaceAge was trained on data from 58,851 healthy individuals, and clinical utility was evaluated on data from 6,196 patients with cancer diagnoses from two institutions in the United States and The Netherlands. To assess the prognostic relevance of FaceAge estimation, we performed Kaplan Meier survival analysis. To test a relevant clinical application of FaceAge, we assessed the performance of FaceAge in end-of-life patients with metastatic cancer who received palliative treatment by incorporating FaceAge into clinical prediction models. We found that, on average, cancer patients look older than their chronological age, and looking older is correlated with worse overall survival. FaceAge demonstrated significant independent prognostic performance in a range of cancer types and stages. We found that FaceAge can improve physicians' survival predictions in incurable patients receiving palliative treatments, highlighting the clinical utility of the algorithm to support end-of-life decision-making. FaceAge was also significantly associated with molecular mechanisms of senescence through gene analysis, while age was not. These findings may extend to diseases beyond cancer, motivating using deep learning algorithms to translate a patient's visual appearance into objective, quantitative, and clinically useful measures.

Genetic association-based functional analysis detects HOGA1 as a potential gene involved in fat accumulation

Although there are a number of discoveries from genome-wide association studies (GWAS) for obesity, it has not been successful in linking GWAS results to biology. We sought to discover causal genes for obesity by conducting functional studies on genes detected from genetic association analysis. Gene-based association analysis of 917 individual exome sequences showed that HOGA1 attains exome-wide significance (p-value &lt; 2.7 × 10–6) for body mass index (BMI). The mRNA expression of HOGA1 is significantly increased in human adipose tissues from obese individuals in the Genotype-Tissue Expression (GTEx) dataset, which supports the genetic association of HOGA1 with BMI. Functional analyses employing cell- and animal model-based approaches were performed to gain insights into the functional relevance of Hoga1 in obesity. Adipogenesis was retarded when Hoga1 was knocked down by siRNA treatment in a mouse 3T3-L1 cell line and a similar inhibitory effect was confirmed in mice with down-regulated Hoga1. Hoga1 antisense oligonucleotide (ASO) treatment reduced body weight, blood lipid level, blood glucose, and adipocyte size in high-fat diet-induced mice. In addition, several lipogenic genes including Srebf1, Scd1, Lp1, and Acaca were down-regulated, while lipolytic genes Cpt1l, Ppara, and Ucp1 were up-regulated. Taken together, HOGA1 is a potential causal gene for obesity as it plays a role in excess body fat development.

Revisiting the Wald Test in Small Case-Control Studies With a Skewed Covariate

The Wald test is routinely used in case-control studies to test for association between a covariate and disease. However, when the evidence for association is high, the Wald test tends to inflate small P values as a result of the Hauck-Donner effect (HDE). Here, we investigate the HDE in the context of genetic burden, both with and without additional covariates. First, we examine the burden-based P values in the absence of association using whole-exome sequence data from 1000 Genomes Project reference samples (n = 54) and selected preterm infants with neonatal complications (n = 74). Our careful analysis of the burden-based P values shows that the HDE is present and that the cause of the HDE in this setting is likely a natural extension of the well-known cause of the HDE in 2 × 2 contingency tables. Second, in a reanalysis of real data, we find that the permutation test provides increased power over the Wald, Firth, and likelihood ratio tests, which agrees with our intuition since the permutation test is valid for any sample size and since it does not suffer from the HDE. Therefore, we propose a powerful and computationally efficient permutation-based approach for the analysis and reanalysis of small case-control association studies.

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.

A Comparison of Methods for Gene-Based Testing That Account for Linkage Disequilibrium

Controlling the type I error rate while retaining sufficient power is a major concern in genome-wide association studies, which nowadays often examine more than a million single-nucleotide polymorphisms (SNPs) simultaneously. Methods such as the Bonferroni correction can lead to a considerable decrease in power due to the large number of tests conducted. Shifting the focus to higher functional structures (e.g., genes) can reduce the loss of power. This can be accomplished via the combination of p-values of SNPs that belong to the same structural unit to test their joint null hypothesis. However, standard methods for this purpose (e.g., Fisher’s method) do not account for the dependence among the tests due to linkage disequilibrium (LD). In this paper, we review various adjustments to methods for combining p-values that take LD information explicitly into consideration and evaluate their performance in a simulation study based on data from the HapMap project. The results illustrate the importance of incorporating LD information into the methods for controlling the type I error rate at the desired level. Furthermore, some methods are more successful in controlling the type I error rate than others. Among them, Brown’s method was the most robust technique with respect to the characteristics of the genes and outperformed the Bonferroni method in terms of power in many scenarios. Examining the genetic factors of a phenotype of interest at the gene-rather than SNP-level can provide researchers benefits in terms of the power of the study. While doing so, one should be careful to account for LD in SNPs belonging to the same gene, for which Brown’s method seems the most robust technique.

A computationally efficient clustering linear combination approach to jointly analyze multiple phenotypes for GWAS

There has been an increasing interest in joint analysis of multiple phenotypes in genome-wide association studies (GWAS) because jointly analyzing multiple phenotypes may increase statistical power to detect genetic variants associated with complex diseases or traits. Recently, many statistical methods have been developed for joint analysis of multiple phenotypes in genetic association studies, including the Clustering Linear Combination (CLC) method. The CLC method works particularly well with phenotypes that have natural groupings, but due to the unknown number of clusters for a given data, the final test statistic of CLC method is the minimum p-value among all p-values of the CLC test statistics obtained from each possible number of clusters. Therefore, a simulation procedure needs to be used to evaluate the p-value of the final test statistic. This makes the CLC method computationally demanding. We develop a new method called computationally efficient CLC (ceCLC) to test the association between multiple phenotypes and a genetic variant. Instead of using the minimum p-value as the test statistic in the CLC method, ceCLC uses the Cauchy combination test to combine all p-values of the CLC test statistics obtained from each possible number of clusters. The test statistic of ceCLC approximately follows a standard Cauchy distribution, so the p-value can be obtained from the cumulative density function without the need for the simulation procedure. Through extensive simulation studies and application on the COPDGene data, the results demonstrate that the type I error rates of ceCLC are effectively controlled in different simulation settings and ceCLC either outperforms all other methods or has statistical power that is very close to the most powerful method with which it has been compared.

Sweet Taste Preference: Relationships with Other Tastes, Liking for Sugary Foods and Exploratory Genome-Wide Association Analysis in Subjects with Metabolic Syndrome

Taste perception and its association with nutrition and related diseases (type 2 diabetes, obesity, metabolic syndrome, cardiovascular, etc.) are emerging fields of biomedicine. There is currently great interest in investigating the environmental and genetic factors that influence sweet taste and sugary food preferences for personalized nutrition. Our aims were: (1) to carry out an integrated analysis of the influence of sweet taste preference (both in isolation and in the context of other tastes) on the preference for sugary foods and its modulation by type 2 diabetes status; (2) as well as to explore new genetic factors associated with sweet taste preference. We studied 425 elderly white European subjects with metabolic syndrome and analyzed taste preference, taste perception, sugary-foods liking, biochemical and genetic markers. We found that type 2 diabetic subjects (38%) have a small, but statistically higher preference for sweet taste (p = 0.021) than non-diabetic subjects. No statistically significant differences (p > 0.05) in preferences for the other tastes (bitter, salty, sour or umami) were detected. For taste perception, type 2 diabetic subjects have a slightly lower perception of all tastes (p = 0.026 for the combined “total taste score”), bitter taste being statistically lower (p = 0.023). We also carried out a principal component analysis (PCA), to identify latent variables related to preferences for the five tastes. We identified two factors with eigenvalues >1. Factor 2 was the one with the highest correlation with sweet taste preference. Sweet taste preference was strongly associated with a liking for sugary foods. In the exploratory SNP-based genome-wide association study (GWAS), we identified some SNPs associated with sweet taste preference, both at the suggestive and at the genome-wide level, especially a lead SNP in the PTPRN2 (Protein Tyrosine Phosphatase Receptor Type N2) gene, whose minor allele was associated with a lower sweet taste preference. The PTPRN2 gene was also a top-ranked gene obtained in the gene-based exploratory GWAS analysis. In conclusion, sweet taste preference was strongly associated with sugary food liking in this population. Our exploratory GWAS identified an interesting candidate gene related with sweet taste preference, but more studies in other populations are required for personalized nutrition.

Different genes may be involved in distal and local sensitisation: a genome-wide gene-based association study and meta-analysis

Background

Neuropathic pain symptoms and signs of increased pain sensitization in osteoarthritis (OA) patients may explain persistent pain after total joint replacement (TJR). Therefore, identifying genetic markers associated with pain sensitization and neuropathic‐like pain phenotypes could be clinically important in identifying targets for early intervention.

Methods

We performed a genome‐wide gene‐based association study (GWGAS) using pressure pain detection thresholds (PPTs) from distal pain‐free sites (anterior tibia), a measure of distal sensitization, and from proximal pain‐affected sites (lateral joint line), a measure of local sensitization, in 320 knee OA participants from the Knee Pain and related health in the Community (KPIC) cohort. We next performed gene‐based fixed‐effects meta‐analysis of PPTs and a neuropathic‐like pain phenotype using genome‐wide association study (GWAS) data from KPIC and from an independent cohort of 613 post‐TJR participants, respectively.

Results

The most significant genes associated with distal and local sensitization were OR5B3 and BRDT, respectively. We also found previously identified neuropathic pain‐associated genes—KCNA1, MTOR, ADORA1 and SCN3B—associated with PPT at the anterior tibia and an inflammatory pain gene—PTAFR—associated with PPT at the lateral joint line. Meta‐analysis results of anterior tibia and neuropathic‐like pain phenotypes revealed genes associated with bone morphogenesis, neuro‐inflammation, obesity, type 2 diabetes, cardiovascular disease and cognitive function.

Conclusions

Overall, our results suggest that different biological processes might be involved in distal and local sensitization, and common genetic mechanisms might be implicated in distal sensitization and neuropathic‐like pain. Future studies are needed to replicate these findings.

Significance

To the best of our knowledge, this is the first GWAS for pain sensitization and the first gene‐based meta‐analysis of pain sensitization and neuropathic‐like pain. Higher pain sensitization and neuropathic pain symptoms are associated with persistent pain after surgery hence, identifying genetic biomarkers and molecular pathways associated with these traits is clinically relevant.

Advances in Genetics and Epigenetic Alterations in Alzheimer's Disease: A Notion for Therapeutic Treatment

Alzheimer's disease (AD) is a disabling neurodegenerative disorder that leads to long-term functional and cognitive impairment and greatly reduces life expectancy. Early genetic studies focused on tracking variations in genome-wide DNA sequences discovered several polymorphisms and novel susceptibility genes associated with AD. However, despite the numerous risk factors already identified, there is still no fully satisfactory explanation for the mechanisms underlying the onset of the disease. Also, as with other complex human diseases, the causes of low heritability are unclear. Epigenetic mechanisms, in which changes in gene expression do not depend on changes in genotype, have attracted considerable attention in recent years and are key to understanding the processes that influence age-related changes and various neurological diseases. With the recent use of massive sequencing techniques, methods for studying epigenome variations in AD have also evolved tremendously, allowing the discovery of differentially expressed disease traits under different conditions and experimental settings. This is important for understanding disease development and for unlocking new potential AD therapies. In this work, we outline the genomic and epigenomic components involved in the initiation and development of AD and identify potentially effective therapeutic targets for disease control.

Identification of novel genes for triple-negative breast cancer with semiparametric gene-based analysis

Triple-negative breast cancer (TNBC) is generally considered an aggressive breast cancer subtype associated with poor prognostic outcomes. Up to now, the molecular and cellular mechanisms underlying TNBC pathology have not been fully understood. In this manuscript, we propose a novel semiparametric model with kernel for gene-based analysis with a breast cancer GWAS data. The software of SPMGBA (semiparametric method for gene-based analysis) in MATLAB is available at GitHub (https://github.com/zliu3/SPMGBA). Genetic signatures associated with breast cancer are discovered. We further validate the prognostic power of the identified genes with a large cohort of expression data from the European Genome-Phenome Archive, and discover that SEL1L is associated with the overall survival of TNBC with the p-value of .0002. We conclude that gene SEL1L is down-regulated in TNBC and the expression of SEL1L is positively associated with patient survival.

Multivariate Analysis Identifies Eight Novel Loci Associated with Meat Productivity Traits in Sheep

Despite their economic value, sheep remain relatively poorly studied animals in terms of the number of known loci and genes associated with commercially important traits. This gap in our knowledge can be filled in by performing new genome-wide association studies (GWAS) or by re-analyzing previously documented data using novel powerful statistical methods. This study is focused on the search for new loci associated with meat productivity and carcass traits in sheep. With a multivariate approach applied to publicly available GWAS results, we identified eight novel loci associated with the meat productivity and carcass traits in sheep. Using an in silico follow-up approach, we prioritized 13 genes in these loci. One of eight novel loci near the FAM3C and WNT16 genes has been replicated in an independent sample of Russian sheep populations (N = 108). The novel loci were added to our regularly updated database increasing the number of known loci to more than 140.

Integrating Genome-Wide Association Analysis With Transcriptome Sequencing to Identify Candidate Genes Related to Blooming Time in Prunus mume

Prunus mume is one of the most important woody perennials for edible and ornamental use. Despite a substantial variation in the flowering phenology among the P. mume germplasm resources, the genetic control for flowering time remains to be elucidated. In this study, we examined five blooming time-related traits of 235 P. mume landraces for 2 years. Based on the phenotypic data, we performed genome-wide association studies, which included a combination of marker- and gene-based association tests, and identified 1,445 candidate genes that are consistently linked with flowering time across multiple years. Furthermore, we assessed the global transcriptome change of floral buds from the two P. mume cultivars exhibiting contrasting bloom dates and detected 617 associated genes that were differentially expressed during the flowering process. By integrating a co-expression network analysis, we screened out 191 gene candidates of conserved transcriptional pattern during blooming across cultivars. Finally, we validated the temporal expression profiles of these candidates and highlighted their putative roles in regulating floral bud break and blooming time in P. mume. Our findings are important to expand the understanding of flowering time control in woody perennials and will boost the molecular breeding of novel varieties in P. mume.

Genome-wide pathway-based quantitative multiple phenotypes analysis

For complex diseases, genome-wide pathway association studies have become increasingly promising. Currently, however, pathway-based association analysis mainly focus on a single phenotype, which may insufficient to describe the complex diseases and physiological processes. This work proposes a combination model to evaluate the association between a pathway and multiple phenotypes and to reduce the run time based on asymptotic results. For a single phenotype, we propose a semi-supervised maximum kernel-based U-statistics (mSKU) method to assess the pathway-based association analysis. For multiple phenotypes, we propose the fisher combination function with dependent phenotypes (FC) to transform the p-values between the pathway and each marginal phenotype individually to achieve pathway-based multiple phenotypes analysis. With real data from the Alzheimer Disease Neuroimaging Initiative (ADNI) study and Human Liver Cohort (HLC) study, the FC-mSKU method allows us to specify which pathways are specific to a single phenotype or contribute to common genetic constructions of multiple phenotypes. If we only focus on single-phenotype tests, we may miss some findings for etiology studies. Through extensive simulation studies, the FC-mSKU method demonstrates its advantages compared with its counterparts.

Multi-trait analysis of rare-variant association summary statistics using MTAR

Integrating association evidence across multiple traits can improve the power of gene discovery and reveal pleiotropy. Most multi-trait analysis methods focus on individual common variants in genome-wide association studies. Here, we introduce multi-trait analysis of rare-variant associations (MTAR), a framework for joint analysis of association summary statistics between multiple rare variants and different traits. MTAR achieves substantial power gain by leveraging the genome-wide genetic correlation measure to inform the degree of gene-level effect heterogeneity across traits. We apply MTAR to rare-variant summary statistics for three lipid traits in the Global Lipids Genetics Consortium. 99 genome-wide significant genes were identified in the single-trait-based tests, and MTAR increases this to 139. Among the 11 novel lipid-associated genes discovered by MTAR, 7 are replicated in an independent UK Biobank GWAS analysis. Our study demonstrates that MTAR is substantially more powerful than single-trait-based tests and highlights the value of MTAR for novel gene discovery.

MetaPhat: Detecting and Decomposing Multivariate Associations From Univariate Genome-Wide Association Statistics

Background

Multivariate testing tools that integrate multiple genome-wide association studies (GWAS) have become important as the number of phenotypes gathered from study cohorts and biobanks has increased. While these tools have been shown to boost statistical power considerably over univariate tests, an important remaining challenge is to interpret which traits are driving the multivariate association and which traits are just passengers with minor contributions to the genotype-phenotypes association statistic.

Results

We introduce MetaPhat, a novel bioinformatics tool to conduct GWAS of multiple correlated traits using univariate GWAS results and to decompose multivariate associations into sets of central traits based on intuitive trace plots that visualize Bayesian Information Criterion (BIC) and P-value statistics of multivariate association models. We validate MetaPhat with Global Lipids Genetics Consortium GWAS results, and we apply MetaPhat to univariate GWAS results for 21 heritable and correlated polyunsaturated lipid species from 2,045 Finnish samples, detecting seven independent loci associated with a cluster of lipid species. In most cases, we are able to decompose these multivariate associations to only three to five central traits out of all 21 traits included in the analyses. We release MetaPhat as an open source tool written in Python with built-in support for multi-processing, quality control, clumping and intuitive visualizations using the R software.

Conclusion

MetaPhat efficiently decomposes associations between multivariate phenotypes and genetic variants into smaller sets of central traits and improves the interpretation and specificity of genome-phenome associations. MetaPhat is freely available under the MIT license at: https://sourceforge.net/projects/meta-pheno-association-tracer.

Information Provided to Consumers about Direct-to-Consumer Nutrigenetic Testing

BACKGROUND

Nutrigenetic tests are often considered to be less serious compared to other health-related genetic tests, although they share similar ethical concerns. Nutrigenetic tests are mainly available through direct-to-consumer genetic testing (DTC GT) and increasing in popularity.

OBJECTIVE

To analyze the contents of nutrigenetic DTC GT websites with respect to the adequacy of the information provided to support a well-informed decision of purchasing the tests.

METHODS

The websites of DTC GT companies selling nutrigenetic tests that could be ordered online without involving any healthcare professional, available in English, marketing tests in Europe, the USA, Australia, or Canada, and accessible from Finland were included in the study (n = 38). Quantitative and qualitative content analyses of the websites were carried out with the help of a codebook.

RESULTS

Of the 38 websites, 8 included a clearly identifiable and easy-to-find information section about genetics. The quality and contents of these sections were often insufficient and/or misleading. Fourteen websites had specific sections discussing the risks related to GT, and on 13 signed informed consent was requested for GT. Furthermore, only 2 of the companies offered any kind of pretest consultation and 13 offered mostly separately charged posttest consultation. The complex structure of the websites made it difficult to find all key information, with many important aspects buried in legal documents, which were challenging to comprehend even for a professional.

CONCLUSION

The structure of the websites and the amount and quality of the content therein do not support a well-informed decision.