Strategies For Improving Epistasis Detection And Replication

Genome-wide association studies (GWAS) have been extensively critiqued for their perceived inability to adequately elucidate the genetic underpinnings of complex disease. Of particular concern is “missing heritability,” or the difference between the total estimated heritability of a phenotype and that explained by GWAS-identified loci. There are numerous proposed explanations for this missing heritability, but a frequently ignored and potentially vastly informative alternative explanation is the ubiquity of epistasis underlying complex phenotypes. Given our understanding of how biomolecules interact in networks and pathways, it is not unreasonable to conclude that the effect of variation at individual genetic loci may non-additively depend on and should be analyzed in the context of their interacting partners. It has been recognized for over a century that deviation from expected Mendelian proportions can be explained by the interaction of multiple loci, and the epistatic underpinnings of phenotypes in model organisms have been extensively experimentally quantified. Therefore, the dearth of inspiring single locus GWAS hits for complex human phenotypes (and the inconsistent replication of these between populations) should not be surprising, as one might expect the joint effect of multiple perturbations to interacting partners within a functional biological module to be more important than individual main effects. Current methods for analyzing data from GWAS are not well-equipped to detect epistasis or replicate significant interactions. The multiple testing burden associated with testing each pairwise interaction quickly becomes nearly insurmountable with increasing numbers of loci. Statistical and machine learning approaches that have worked well for other types of high-dimensional data are appealing and may be useful for detecting epistasis, but potentially require tweaks to function appropriately. Biological knowledge may also be leveraged to guide the search for epistasis candidates, but requires context-appropriate application (as, for example, two loci with significant main effects may not have a significant interaction, and vice versa). Rather than renouncing GWAS and the wealth of associated data that has been accumulated as a failure, I propose the development of new techniques and incorporation of diverse data sources to analyze GWAS data in an epistasis-centric framework

Evidence for the contribution of COMT gene Val158/108Met polymorphism (rs4680) to working memory training-related prefrontal plasticity.

BackgroundGenetic factors have been suggested to affect the efficacy of working memory training. However, few studies have attempted to identify the relevant genes.MethodsIn this study, we first performed a randomized controlled trial (RCT) to identify brain regions that were specifically affected by working memory training. Sixty undergraduate students were randomly assigned to either the adaptive training group (N = 30) or the active control group (N = 30). Both groups were trained for 20 sessions during 4 weeks and received fMRI scans before and after the training. Afterward, we combined the data from the 30 participants in the RCT study who received adaptive training with data from 71 additional participants who also received the same adaptive training but were not part of the RCT study (total N = 101) to test the contribution of the COMT Val158/108Met polymorphism to the interindividual difference in the training effect within the identified brain regions.ResultsIn the RCT study, we found that the adaptive training significantly decreased brain activation in the left prefrontal cortex (TFCE-FWE corrected p = .030). In the genetic study, we found that compared with the Val allele homozygotes, the Met allele carriers' brain activation decreased more after the training at the left prefrontal cortex (TFCE-FWE corrected p = .025).ConclusionsThis study provided evidence for the neural effect of a visual-spatial span training and suggested that genetic factors such as the COMT Val158/108Met polymorphism may have to be considered in future studies of such training

A multivariant recall-by-genotype study of the metabolomic signature of BMI

OBJECTIVE: This study estimated the effect of BMI on circulating metabolites in young adults using a recall‐by‐genotype study design. METHODS: A recall‐by‐genotype study was implemented in the Avon Longitudinal Study of Parents and Children. Samples from 756 participants were selected for untargeted metabolomics analysis based on low versus high genetic liability for higher BMI defined by a genetic risk score (GRS). Regression analyses were performed to investigate associations between BMI GRS group and relative abundance of 973 metabolites. RESULTS: After correction for multiple testing, 29 metabolites were associated with BMI GRS group. Bilirubin was among the most strongly associated metabolites, with reduced levels measured in individuals in the high‐BMI GRS group (β = −0.32, 95% CI: −0.46 to −0.18, Benjamini‐Hochberg adjusted p = 0.005). This study observed associations between BMI GRS group and the levels of several potentially diet‐related metabolites, including hippurate, which had lower mean abundance in individuals in the high‐BMI GRS group (β = −0.29, 95% CI: −0.44 to −0.15, Benjamini‐Hochberg adjusted p = 0.008). CONCLUSIONS: Together with existing literature, these results suggest that a genetic predisposition to higher BMI captures differences in metabolism leading to adiposity gain. In the absence of prospective data, separating these effects from the downstream consequences of weight gain is challenging

Prevalence of sexual dimorphism in mammalian phenotypic traits

The role of sex in biomedical studies has often been overlooked, despite evidence of sexually dimorphic effects in some biological studies. Here, we used high-throughput phenotype data from 14,250 wildtype and 40,192 mutant mice (representing 2,186 knockout lines), analysed for up to 234 traits, and found a large proportion of mammalian traits both in wildtype and mutants are influenced by sex. This result has implications for interpreting disease phenotypes in animal models and humans

Prevalence of sexual dimorphism in mammalian phenotypic traits.

The role of sex in biomedical studies has often been overlooked, despite evidence of sexually dimorphic effects in some biological studies. Here, we used high-throughput phenotype data from 14,250 wildtype and 40,192 mutant mice (representing 2,186 knockout lines), analysed for up to 234 traits, and found a large proportion of mammalian traits both in wildtype and mutants are influenced by sex. This result has implications for interpreting disease phenotypes in animal models and humans

A cautionary note on the impact of protocol changes for Genome-Wide Association SNP x SNP Interaction studies: an example on ankylosing spondylitis

Genome-wide association interaction (GWAI) studies have increased in popularity. Yet to date, no standard protocol exists. In practice, any GWAI workflow involves making choices about quality control strategy, SNP filtering, linkage disequilibrium (LD) pruning, analytic tool to model or to test for genetic interactions. Each of these can have an impact on the final epistasis findings and may affect their reproducibility in follow-up analyses. Choosing an analytic tool is not straightforward, as different such tools exist and current understanding about their performance is based on often very particular simulation settings. In the present study, we wish to create awareness for the impact of (minor) changes in a GWAI analysis protocol can have on final epistasis findings. In particular, we investigate the influence of marker selection and marker prioritization strategies, LD pruning and the choice of epistasis detection analytics on study results, giving rise to 8 GWAI protocols. Discussions are made in the context of the ankylosing spondylitis (AS) data obtained via the Wellcome Trust Case Control Consortium (WTCCC2). As expected, the largest impact on AS epistasis findings is caused by the choice of marker selection criterion, followed by marker coding and LD pruning. In MB-MDR, co-dominant coding of main effects is more robust to the effects of LD pruning than additive coding. We were able to reproduce previously reported epistasis involvement of HLA-B and ERAP1 in AS pathology. In addition, our results suggest involvement of MAGI3 and PARK2, responsible for cell adhesion and cellular trafficking. Gene Ontology (GO) biological function enrichment analysis across the 8 considered GWAI protocols also suggested that AS could be associated to the Central Nervous System (CNS) malfunctions, specifically, in nerve impulse propagation and in neurotransmitters metabolic processes

The emerging landscape of health research based on biobanks linked to electronic health records: Existing resources, statistical challenges, and potential opportunities

Peer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/154448/1/sim8445_am.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/154448/2/sim8445.pd

Computational genomics in the era of precision medicine: Applications to variant analysis and gene therapy

Rapid methodological advances in statistical and computational genomics have enabled researchers to better identify and interpret both rare and common variants responsible for complex human diseases. As we continue to see an expansion of these advances in the field, it is now imperative for researchers to understand the resources and methodologies available for various data types and study designs. In this review, we provide an overview of recent methods for identifying rare and common variants and understanding their roles in disease etiology. Additionally, we discuss the strategy, challenge, and promise of gene therapy. As computational and statistical approaches continue to improve, we will have an opportunity to translate human genetic findings into personalized health care

A guide for understanding and designing Mendelian randomization studies in the musculoskeletal field

Mendelian randomization (MR) is an increasingly popular component of an epidemiologist's toolkit, used to provide evidence of a causal effect of one trait (an exposure, eg, body mass index [BMI]) on an outcome trait or disease (eg, osteoarthritis). Identifying these effects is important for understanding disease etiology and potentially identifying targets for therapeutic intervention. MR uses genetic variants as instrumental variables for the exposure, which should not be influenced by the outcome or confounding variables, overcoming key limitations of traditional epidemiological analyses. For MR to generate a valid estimate of effect, key assumptions must be met. In recent years, there has been a rapid rise in MR methods that aim to test, or are robust to violations of, these assumptions. In this review, we provide an overview of MR for a non‐expert audience, including an explanation of these key assumptions and how they are often tested, to aid a better reading and understanding of the MR literature. We highlight some of these new methods and how they can be useful for specific methodological challenges in the musculoskeletal field, including for conditions or traits that share underlying biological pathways, such as bone and joint disease. © 2022 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research

Genetic Diversity of the Cestode Echinococcus multilocularis in Red Foxes at a Continental Scale in Europe

Echinococcus multilocularis is a tapeworm of the red fox, which represents a considerable health threat to respectively infected humans. Main endemic areas are located in China, Siberia, and central Europe. Alarmed by an emerging or reemerging situation in Europe, the question of how the parasite gets spatially and temporally spread and transmitted becomes essential to prepare appropriate control programs. The question was tackled by using genetic data on a large sample size of E. multilocularis adult stage tapeworms, combined with geographical site location data input. The historically documented endemic area, represented by the northern Alpine arch, was shown to harbour the highest genetic richness and diversity, as compared to surrounding areas in northern and eastern Europe. The spatial and temporal spread of different E. multilocularis genotypes in Europe seems to be ruled by a founder event, linked to exportation of parasites from the central core to newly identified (western and eastern) areas or subregions, where these parasites could subsequently disseminate under geographical separation from the original foci