Meta-analysis identifies gene-by-environment interactions as demonstrated in a study of 4,965 mice

GenomeMUSter mouse genetic variation service enables multitrait, multipopulation data integration and analysis

Hundreds of inbred mouse strains and intercross populations have been used to characterize the function of genetic variants that contribute to disease. Thousands of disease-relevant traits have been characterized in mice and made publicly available. New strains and populations including consomics, the collaborative cross, expanded BXD, and inbred wild-derived strains add to existing complex disease mouse models, mapping populations, and sensitized backgrounds for engineered mutations. The genome sequences of inbred strains, along with dense genotypes from others, enable integrated analysis of trait-variant associations across populations, but these analyses are hampered by the sparsity of genotypes available. Moreover, the data are not readily interoperable with other resources. To address these limitations, we created a uniformly dense variant resource by harmonizing multiple data sets. Missing genotypes were imputed using the Viterbi algorithm with a data-driven technique that incorporates local phylogenetic information, an approach that is extendable to other model organisms. The result is a web- and programmatically accessible data service called GenomeMUSter, comprising single-nucleotide variants covering 657 strains at 106.8 million segregating sites. Interoperation with phenotype databases, analytic tools, and other resources enable a wealth of applications, including multitrait, multipopulation meta-analysis. We show this in cross-species comparisons of type 2 diabetes and substance use disorder meta-analyses, leveraging mouse data to characterize the likely role of human variant effects in disease. Other applications include refinement of mapped loci and prioritization of strain backgrounds for disease modeling to further unlock extant mouse diversity for genetic and genomic studies in health and disease.

Identification of genetic drivers of plasma lipoprotein size in the Diversity Outbred mouse population

Despite great progress in understanding lipoprotein physiology, there is still much to be learned about the genetic drivers of lipoprotein abundance, composition, and function. We used ion mobility spectrometry to survey 16 plasma lipoprotein subfractions in 500 Diversity Outbred mice maintained on a Western-style diet. We identified 21 quantitative trait loci (QTL) affecting lipoprotein abundance. To refine the QTL and link them to disease risk in humans, we asked if the human homologs of genes located at each QTL were associated with lipid traits in human genome-wide association studies. Integration of mouse QTL with human genome-wide association studies yielded candidate gene drivers for 18 of the 21 QTL. This approach enabled us to nominate the gene encoding the neutral ceramidase, Asah2, as a novel candidate driver at a QTL on chromosome 19 for large HDL particles (HDL-2b). To experimentally validate Asah2, we surveyed lipoproteins in Asah2-/- mice. Compared to wild-type mice, female Asah2-/- mice showed an increase in several lipoproteins, including HDL. Our results provide insights into the genetic regulation of circulating lipoproteins, as well as mechanisms by which lipoprotein subfractions may affect cardiovascular disease risk in humans.

The potential of integrating human and mouse discovery platforms to advance our understanding of cardiometabolic diseases

Cardiometabolic diseases encompass a range of interrelated conditions that arise from underlying metabolic perturbations precipitated by genetic, environmental, and lifestyle factors. While obesity, dyslipidaemia, smoking, and insulin resistance are major risk factors for cardiometabolic diseases, individuals still present in the absence of such traditional risk factors, making it difficult to determine those at greatest risk of disease. Thus, it is crucial to elucidate the genetic, environmental, and molecular underpinnings to better understand, diagnose, and treat cardiometabolic diseases. Much of this information can be garnered using systems genetics, which takes population-based approaches to investigate how genetic variance contributes to complex traits. Despite the important advances made by human genome-wide association studies (GWAS) in this space, corroboration of these findings has been hampered by limitations including the inability to control environmental influence, limited access to pertinent metabolic tissues, and often, poor classification of diseases or phenotypes. A complementary approach to human GWAS is the utilisation of model systems such as genetically diverse mouse panels to study natural genetic and phenotypic variation in a controlled environment. Here, we review mouse genetic reference panels and the opportunities they provide for the study of cardiometabolic diseases and related traits. We discuss how the post-GWAS era has prompted a shift in focus from discovery of novel genetic variants to understanding gene function. Finally, we highlight key advantages and challenges of integrating complementary genetic and multi-omics data from human and mouse populations to advance biological discovery.

Role of gut microbe-derived metabolites in cardiometabolic diseases: Systems based approach

BACKGROUND

The gut microbiome influences host physiology and cardiometabolic diseases by interacting directly with intestinal cells or by producing molecules that enter the host circulation. Given the large number of microbial species present in the gut and the numerous factors that influence gut bacterial composition, it has been challenging to understand the underlying biological mechanisms that modulate risk of cardiometabolic disease.

SCOPE OF THE REVIEW

Here we discuss a systems-based approach that involves simultaneously examining individuals in populations for gut microbiome composition, molecular traits using "omics" technologies, such as circulating metabolites quantified by mass spectrometry, and clinical traits. We summarize findings from landmark studies using this approach and discuss future applications.

MAJOR CONCLUSIONS

Population-based integrative approaches have identified a large number of microbe-derived or microbe-modified metabolites that are associated with cardiometabolic traits. The knowledge gained from these studies provide new opportunities for understanding the mechanisms involved in gut microbiome-host interactions and may have potentially important implications for developing novel therapeutic approaches.

Cross-ancestry genome-wide meta-analysis of 61,047 cases and 947,237 controls identifies new susceptibility loci contributing to lung cancer

To identify new susceptibility loci to lung cancer among diverse populations, we performed cross-ancestry genome-wide association studies in European, East Asian and African populations and discovered five loci that have not been previously reported. We replicated 26 signals and identified 10 new lead associations from previously reported loci. Rare-variant associations tended to be specific to populations, but even common-variant associations influencing smoking behavior, such as those with CHRNA5 and CYP2A6, showed population specificity. Fine-mapping and expression quantitative trait locus colocalization nominated several candidate variants and susceptibility genes such as IRF4 and FUBP1. DNA damage assays of prioritized genes in lung fibroblasts indicated that a subset of these genes, including the pleiotropic gene IRF4, potentially exert effects by promoting endogenous DNA damage.

Genetic Background Matters: Population-Based Studies in Model Organisms for Translational Research

We are all similar but a bit different. These differences are partially due to variations in our genomes and are related to the heterogeneity of symptoms and responses to treatments that patients exhibit. Most animal studies are performed in one single strain with one manipulation. However, due to the lack of variability, therapies are not always reproducible when treatments are translated to humans. Panels of already sequenced organisms are valuable tools for mimicking human phenotypic heterogeneities and gene mapping. This review summarizes the current knowledge of mouse, fly, and yeast panels with insightful applications for translational research.

Analysis of independent cohorts of outbred CFW mice reveals novel loci for behavioral and physiological traits and identifies factors determining reproducibility

Combining samples for genetic association is standard practice in human genetic analysis of complex traits, but is rarely undertaken in rodent genetics. Here, using 23 phenotypes and genotypes from two independent laboratories, we obtained a sample size of 3076 commercially available outbred mice and identified 70 loci, more than double the number of loci identified in the component studies. Fine-mapping in the combined sample reduced the number of likely causal variants, with a median reduction in set size of 51%, and indicated novel gene associations, including Pnpo, Ttll6, and GM11545 with bone mineral density, and Psmb9 with weight. However, replication at a nominal threshold of 0.05 between the two component studies was low, with less than one-third of loci identified in one study replicated in the second. In addition to overestimates in the effect size in the discovery sample (Winner's Curse), we also found that heterogeneity between studies explained the poor replication, but the contribution of these two factors varied among traits. Leveraging these observations, we integrated information about replication rates, study-specific heterogeneity, and Winner's Curse corrected estimates of power to assign variants to one of four confidence levels. Our approach addresses concerns about reproducibility and demonstrates how to obtain robust results from mapping complex traits in any genome-wide association study.

Estimation of genetic variance contributed by a quantitative trait locus: correcting the bias associated with significance tests

The Beavis effect in quantitative trait locus (QTL) mapping describes a phenomenon that the estimated effect size of a statistically significant QTL (measured by the QTL variance) is greater than the true effect size of the QTL if the sample size is not sufficiently large. This is a typical example of the Winners' curse applied to molecular quantitative genetics. Theoretical evaluation and correction for the Winners' curse have been studied for interval mapping. However, similar technologies have not been available for current models of QTL mapping and genome-wide association studies where a polygene is often included in the linear mixed models to control the genetic background effect. In this study, we developed the theory of the Beavis effect in a linear mixed model using a truncated noncentral Chi-square distribution. We equated the observed Wald test statistic of a significant QTL to the expectation of a truncated noncentral Chi-square distribution to obtain a bias-corrected estimate of the QTL variance. The results are validated from replicated Monte Carlo simulation experiments. We applied the new method to the grain width (GW) trait of a rice population consisting of 524 homozygous varieties with over 300 k single nucleotide polymorphism markers. Two loci were identified and the estimated QTL heritability were corrected for the Beavis effect. Bias correction for the larger QTL on chromosome 5 (GW5) with an estimated heritability of 12% did not change the QTL heritability due to the extremely large test score and estimated QTL effect. The smaller QTL on chromosome 9 (GW9) had an estimated QTL heritability of 9% reduced to 6% after the bias-correction.

Integrative Analyses Reveal Tstd1 as a Potential Modulator of HDL Cholesterol and Mitochondrial Function in Mice

High-density lipoprotein (HDL) cholesterol levels are closely associated with human health and diseases. To identify genes modulating plasma HDL levels, we integrated HDL measurements and multi-omics data collected from diverse mouse cohorts and combined a list of systems genetics methods, including quantitative trait loci (QTL) mapping analysis, mediation analysis, transcriptome-wide association analysis (TWAS), and correlation analysis. We confirmed a significant and conserved QTL for plasma HDL on chromosome 1 and identified that Tstd1 liver transcript correlates with plasma HDL in several independent mouse cohorts, suggesting Tstd1 may be a potential modulator of plasma HDL levels. Correlation analysis using over 70 transcriptomics datasets in humans and mice revealed consistent correlations between Tstd1 and genes known to be involved in cholesterol and HDL regulation. Consistent with strong enrichment in gene sets related to cholesterol and lipoproteins in the liver, mouse strains with high Tstd1 exhibited higher plasma levels of HDL, total cholesterol and other lipid markers. GeneBridge using large-scale expression datasets identified conserved and positive associations between TSTD1/Tstd1 and mitochondrial pathways, as well as cholesterol and lipid pathways in human, mouse and rat. In summary, we identified Tstd1 as a new modulator of plasma HDL and mitochondrial function through integrative systems analyses, and proposed a new mechanism of HDL modulation and a potential therapeutic target for relevant diseases. This study highlights the value of such integrative approaches in revealing molecular mechanisms of complex traits or diseases.

Genome-wide association analyses identify genotype-by-environment interactions of growth traits in Simmental cattle

Understanding genotype-by-environment interactions (G × E) is crucial to understand environmental adaptation in mammals and improve the sustainability of agricultural production. Here, we present an extensive study investigating the interaction of genome-wide SNP markers with a vast assortment of environmental variables and searching for SNPs controlling phenotypic variance (vQTL) using a large beef cattle dataset. We showed that G × E contribute 10.1%, 3.8%, and 2.8% of the phenotypic variance of birth weight, weaning weight, and yearling weight, respectively. G × E genome-wide association analysis (GWAA) detected a large number of G × E loci affecting growth traits, which the traditional GWAA did not detect, showing that functional loci may have non-additive genetic effects regardless of differences in genotypic means. Further, variance-heterogeneity GWAA detected loci enriched with G × E effects without requiring prior knowledge of the interacting environmental factors. Functional annotation and pathway analysis of G × E genes revealed biological mechanisms by which cattle respond to changes in their environment, such as neurotransmitter activity, hypoxia-induced processes, keratinization, hormone, thermogenic and immune pathways. We unraveled the relevance and complexity of the genetic basis of G × E underlying growth traits, providing new insights into how different environmental conditions interact with specific genes influencing adaptation and productivity in beef cattle and potentially across mammals.

Genome-wide association study and Mendelian randomization analysis provide insights for improving rice yield potential

Rice yield per plant has a complex genetic architecture, which is mainly determined by its three component traits: the number of grains per panicle (GPP), kilo-grain weight (KGW), and tillers per plant (TP). Exploring ideotype breeding based on selection for genetically less complex component traits is an alternative route for further improving rice production. To understand the genetic basis of the relationship between rice yield and component traits, we investigated the four traits of two rice hybrid populations (575 + 1495 F₁) in different environments and conducted meta-analyses of genome-wide association study (meta-GWAS). In total, 3589 significant loci for three components traits were detected, while only 3 loci for yield were detected. It indicated that rice yield is mainly controlled by minor-effect loci and hardly to be identified. Selecting quantitative trait locus/gene affected component traits to further enhance yield is recommended. Mendelian randomization design is adopted to investigate the genetic effects of loci on yield through component traits and estimate the genetic relationship between rice yield and its component traits by these loci. The loci for GPP or TP mainly had a positive genetic effect on yield, but the loci for KGW with different direction effects (positive effect or negative effect). Additionally, TP (Beta = 1.865) has a greater effect on yield than KGW (Beta = 1.016) and GPP (Beta = 0.086). Five significant loci for component traits that had an indirect effect on yield were identified. Pyramiding superior alleles of the five loci revealed improved yield. A combination of direct and indirect effects may better contribute to the yield potential of rice. Our findings provided a rationale for using component traits as indirect indices to enhanced rice yield, which will be helpful for further understanding the genetic basis of yield and provide valuable information for improving rice yield potential.

Coupled mixed model for joint genetic analysis of complex disorders with two independently collected data sets

BACKGROUND

In the last decade, Genome-wide Association studies (GWASs) have contributed to decoding the human genome by uncovering many genetic variations associated with various diseases. Many follow-up investigations involve joint analysis of multiple independently generated GWAS data sets. While most of the computational approaches developed for joint analysis are based on summary statistics, the joint analysis based on individual-level data with consideration of confounding factors remains to be a challenge.

RESULTS

In this study, we propose a method, called Coupled Mixed Model (CMM), that enables a joint GWAS analysis on two independently collected sets of GWAS data with different phenotypes. The CMM method does not require the data sets to have the same phenotypes as it aims to infer the unknown phenotypes using a set of multivariate sparse mixed models. Moreover, CMM addresses the confounding variables due to population stratification, family structures, and cryptic relatedness, as well as those arising during data collection such as batch effects that frequently appear in joint genetic studies. We evaluate the performance of CMM using simulation experiments. In real data analysis, we illustrate the utility of CMM by an application to evaluating common genetic associations for Alzheimer's disease and substance use disorder using datasets independently collected for the two complex human disorders. Comparison of the results with those from previous experiments and analyses supports the utility of our method and provides new insights into the diseases. The software is available at https://github.com/HaohanWang/CMM .

Meta-analysis of SNP-environment interaction with heterogeneity for overlapping data

Meta-analysis is a popular method used in genome-wide association studies, by which the results of multiple studies are combined to identify associations. This process generates heterogeneity. Recently, we proposed a random effect model meta-regression method (MR) to study the effect of single nucleotide polymorphism (SNP)-environment interactions. This method takes heterogeneity into account and produces high power. We also proposed a fixed effect model overlapping MR in which the overlapping data is taken into account. In the present study, a random effect model overlapping MR that simultaneously considers heterogeneity and overlapping data is proposed. This method is based on the random effect model MR and the fixed effect model overlapping MR. A new way of solving the logarithm of the determinant of covariance matrices in likelihood functions is also provided. Tests for the likelihood ratio statistic of the SNP-environment interaction effect and the SNP and SNP-environment joint effects are given. In our simulations, null distributions and type I error rates were proposed to verify the suitability of our method, and powers were applied to evaluate the superiority of our method. Our findings indicate that this method is effective in cases of overlapping data with a high heterogeneity.

A genome-wide association and meta-analysis reveal regions associated with seed size in cowpea [Vigna unguiculata (L.) Walp]

This paper combined GWAS, meta-analysis and sequence homology comparison with common bean to identify regions associated with seed size variation in domesticated cowpea. Seed size is an important trait for yield and commercial value in dry-grain cowpea. Seed size varies widely among different cowpea accessions, and the genetic basis of such variation is not yet well understood. To better decipher the genetic basis of seed size, a genome-wide association study (GWAS) and meta-analysis were conducted on a panel of 368 cowpea diverse accessions from 51 countries. Four traits, including seed weight, length, width and density were evaluated across three locations. Using 51,128 single nucleotide polymorphisms covering the cowpea genome, 17 loci were identified for these traits. One locus was common to weight, width and length, suggesting pleiotropy. By integrating synteny-based analysis with common bean, six candidate genes (Vigun05g036000, Vigun05g039600, Vigun05g204200, Vigun08g217000, Vigun11g187000, and Vigun11g191300) which are implicated in multiple functional categories related to seed size such as endosperm development, embryo development, and cell elongation were identified. These results suggest that a combination of GWAS meta-analysis with synteny comparison in a related plant is an efficient approach to identify candidate gene (s) for complex traits in cowpea. The identified loci and candidate genes provide useful information for improving cowpea varieties and for molecular investigation of seed size.

Association of Behavior With Noise-Induced Hearing Loss Among Attendees of an Outdoor Music Festival: A Secondary Analysis of a Randomized Clinical Trial

Importance

To date, factors associated with noise-induced hearing loss at music festivals have not yet been analyzed in a single comprehensive data set. In addition, little is known about the hearing loss-associated behavior of music festival attendees.

Objectives

To assess which factors are associated with the occurrence of a temporary threshold shift (TTS) after music exposure and to investigate the behavior of music festival attendees.

Design, Setting, and Participants

This prospective post hoc analysis gathered data from a randomized, single-blind clinical trial conducted on September 5, 2015, at an outdoor music festival in Amsterdam, the Netherlands. Adult volunteers with normal hearing were recruited via social media from August 26 through September 3, 2015. Intention to use earplugs was an exclusion criterion. Of 86 volunteers assessed, 51 were included. This post hoc analysis was performed from October 3, 2016, through February 27, 2017.

Interventions

Music festival visit for 4.5 hours.

Main Outcomes and Measures

The primary outcome was a TTS on a standard audiogram for the frequencies 3.0- and 4.0-kHz. Multivariable linear regression was performed to determine which factors are associated with a TTS. A questionnaire on behavior, hearing, and tinnitus was distributed to the participants before and after the festival visit.

Results

A total of 51 participants were included (18 men [35%] and 33 women [65%]) with a mean (SD) age of 27 (6) years. Mean (SD) threshold change across 3.0 and 4.0 kHz was 5.4 (5.7) dB for the right ear and 4.0 (6.1) dB for the left ear. Earplug use (absolute difference in the left ear, -6.0 dB [95% CI, -8.7 to -3.2 dB]; in the right ear, -6.4 dB [95% CI, -8.8 to -4.1 dB]), quantity of alcohol use (absolute difference per unit in the left ear, 1.1 dB [95% CI, 0.5 to 1.7 dB]; in the right ear, 0.7 dB [95% CI, 0.1 to 1.4 dB]), drug use (absolute difference in the right ear, 6.0 dB [95% CI, 0.9 to 11.1 dB]), and male sex (absolute difference in the right ear, 4.1 dB [95% CI, 0.3 to 5.9 dB]) were independently associated with hearing loss, with earplug use being the most important factor. Unprotected participants reported significantly worse subjective hearing performance and tinnitus after the festival visit than did participants using earplugs (Cramer V, 0.62 [95% CI, 0.47-0.79] and 0.39 [95% CI, 0.16-0.62], respectively). In the earplug group, the perceived loudness (r = -0.72; 95% CI, -1.00 to -0.43) and appreciation (r = 0.53; 95% CI, 0.29 to 0.78) of music and speech perception (r = 0.21; 95% CI, 0.09 to 0.35) were correlated with the duration of earplug use.

Conclusions and Relevance

The present study identified nonuse of earplugs, use of alcohol and drugs, and male sex as associated with a TTS at an outdoor music festival. Physicians should consider these factors to raise awareness about the combined risk of attending music festivals without using earplugs while consuming alcohol and/or drugs. The intention to use earplugs was correlated with the loudness and appreciation of music with earplugs, which may advocate for the use of personalized earplugs.

Trial Registration

trialregister.nl Identifier: NTR5401.

Molecular Genetic Analysis Subdivided by Adversity Exposure Suggests Etiologic Heterogeneity in Major Depression

OBJECTIVE

The extent to which major depression is the outcome of a single biological mechanism or represents a final common pathway of multiple disease processes remains uncertain. Genetic approaches can potentially identify etiologic heterogeneity in major depression by classifying patients on the basis of their experience of major adverse events.

METHOD

Data are from the China, Oxford, and VCU Experimental Research on Genetic Epidemiology (CONVERGE) project, a study of Han Chinese women with recurrent major depression aimed at identifying genetic risk factors for major depression in a rigorously ascertained cohort carefully assessed for key environmental risk factors (N=9,599). To detect etiologic heterogeneity, genome-wide association studies, heritability analyses, and gene-by-environment interaction analyses were performed.

RESULTS

Genome-wide association studies stratified by exposure to adversity revealed three novel loci associated with major depression only in study participants with no history of adversity. Significant gene-by-environment interactions were seen between adversity and genotype at all three loci, and 13.2% of major depression liability can be attributed to genome-wide interaction with adversity exposure. The genetic risk in major depression for participants who reported major adverse life events (27%) was partially shared with that in participants who did not (73%; genetic correlation=+0.64). Together with results from simulation studies, these findings suggest etiologic heterogeneity within major depression as a function of environmental exposures.

CONCLUSIONS

The genetic contributions to major depression may differ between women with and those without major adverse life events. These results have implications for the molecular dissection of major depression and other complex psychiatric and biomedical diseases.

An Association Mapping Framework To Account for Potential Sex Difference in Genetic Architectures

Over the past few years, genome-wide association studies have identified many trait-associated loci that have different effects on females and males, which increased attention to the genetic architecture differences between the sexes. The between-sex differences in genetic architectures can cause a variety of phenomena such as differences in the effect sizes at trait-associated loci, differences in the magnitudes of polygenic background effects, and differences in the phenotypic variances. However, current association testing approaches for dealing with sex, such as including sex as a covariate, cannot fully account for these phenomena and can be suboptimal in statistical power. We present a novel association mapping framework, MetaSex, that can comprehensively account for the genetic architecture differences between the sexes. Through simulations and applications to real data, we show that our framework has superior performance than previous approaches in association mapping.

Increasing the power of meta-analysis of genome-wide association studies to detect heterogeneous effects

Motivation

Meta-analysis is essential to combine the results of genome-wide association studies (GWASs). Recent large-scale meta-analyses have combined studies of different ethnicities, environments and even studies of different related phenotypes. These differences between studies can manifest as effect size heterogeneity. We previously developed a modified random effects model (RE2) that can achieve higher power to detect heterogeneous effects than the commonly used fixed effects model (FE). However, RE2 cannot perform meta-analysis of correlated statistics, which are found in recent research designs, and the identified variants often overlap with those found by FE.

Results

Here, we propose RE2C, which increases the power of RE2 in two ways. First, we generalized the likelihood model to account for correlations of statistics to achieve optimal power, using an optimization technique based on spectral decomposition for efficient parameter estimation. Second, we designed a novel statistic to focus on the heterogeneous effects that FE cannot detect, thereby, increasing the power to identify new associations. We developed an efficient and accurate p-value approximation procedure using analytical decomposition of the statistic. In simulations, RE2C achieved a dramatic increase in power compared with the decoupling approach (71% vs. 21%) when the statistics were correlated. Even when the statistics are uncorrelated, RE2C achieves a modest increase in power. Applications to real genetic data supported the utility of RE2C. RE2C is highly efficient and can meta-analyze one hundred GWASs in one day.

Availability and implementation

The software is freely available at http://software.buhmhan.com/RE2C.

Supplementary information

Supplementary data are available at Bioinformatics online.

Applying meta-analysis to genotype-tissue expression data from multiple tissues to identify eQTLs and increase the number of eGenes

Motivation

There is recent interest in using gene expression data to contextualize findings from traditional genome-wide association studies (GWAS). Conditioned on a tissue, expression quantitative trait loci (eQTLs) are genetic variants associated with gene expression, and eGenes are genes whose expression levels are associated with genetic variants. eQTLs and eGenes provide great supporting evidence for GWAS hits and important insights into the regulatory pathways involved in many diseases. When a significant variant or a candidate gene identified by GWAS is also an eQTL or eGene, there is strong evidence to further study this variant or gene. Multi-tissue gene expression datasets like the Gene Tissue Expression (GTEx) data are used to find eQTLs and eGenes. Unfortunately, these datasets often have small sample sizes in some tissues. For this reason, there have been many meta-analysis methods designed to combine gene expression data across many tissues to increase power for finding eQTLs and eGenes. However, these existing techniques are not scalable to datasets containing many tissues, like the GTEx data. Furthermore, these methods ignore a biological insight that the same variant may be associated with the same gene across similar tissues.

Results

We introduce a meta-analysis model that addresses these problems in existing methods. We focus on the problem of finding eGenes in gene expression data from many tissues, and show that our model is better than other types of meta-analyses.

Availability and Implementation

Source code is at https://github.com/datduong/RECOV.

Supplementary information

Supplementary data are available at Bioinformatics online.

An ancestry-based approach for detecting interactions

BACKGROUND

Epistasis and gene-environment interactions are known to contribute significantly to variation of complex phenotypes in model organisms. However, their identification in human association studies remains challenging for myriad reasons. In the case of epistatic interactions, the large number of potential interacting sets of genes presents computational, multiple hypothesis correction, and other statistical power issues. In the case of gene-environment interactions, the lack of consistently measured environmental covariates in most disease studies precludes searching for interactions and creates difficulties for replicating studies.

RESULTS

In this work, we develop a new statistical approach to address these issues that leverages genetic ancestry, defined as the proportion of ancestry derived from each ancestral population (e.g., the fraction of European/African ancestry in African Americans), in admixed populations. We applied our method to gene expression and methylation data from African American and Latino admixed individuals, respectively, identifying nine interactions that were significant at P<5×10-8. We show that two of the interactions in methylation data replicate, and the remaining six are significantly enriched for low P-values (P<1.8×10-6).

CONCLUSION

We show that genetic ancestry can be a useful proxy for unknown and unmeasured covariates in the search for interaction effects. These results have important implications for our understanding of the genetic architecture of complex traits.

Genetics of metabolic syndrome: potential clues from wild-derived inbred mouse strains

The metabolic syndrome (MetS) is a complex constellation of metabolic abnormalities including obesity, abnormal glucose metabolism, dyslipidemia, and elevated blood pressure that together substantially increase risk for cardiovascular disease and Type 2 diabetes. Both genetic and environmental factors contribute to the development of MetS, but this process is still far from understood. Human studies have revealed only part of the underlying basis. Studies in mice offer many strengths that can complement human studies to help elucidate the etiology and pathophysiology of MetS. Here we review the ways mice can contribute to MetS research. In particular, we focus on the information that can be obtained from studies of the inbred strains, with specific focus on the phenotypes of the wild-derived inbred strains. These are newly derived inbred strains that were created from wild-caught mice. They contain substantial genetic variation that is not present in the classical inbred strains, have phenotypes of relevance for MetS, and various mouse strain resources have been created to facilitate the mining of this new genetic variation. Thus studies using wild-derived inbred strains hold great promise for increasing our understanding of MetS.

Adiposity and Genetic Factors in Relation to Triglycerides and Triglyceride-Rich Lipoproteins in the Women's Genome Health Study

BACKGROUND

Previous results from Scandinavian cohorts have shown that obesity accentuates the effects of common genetic susceptibility variants on increased triglycerides (TG). Whether such interactions are present in the US population and further selective for particular TG-rich lipoprotein subfractions is unknown.

METHODS

We examined these questions using body mass index (BMI) and waist circumference (WC) among women of European ancestry from the Women's Genome Health Study (WGHS) (n = 21840 for BMI; n = 19313 for WC). A weighted genetic risk score (TG-wGRS) based on 40 published TG-associated single-nucleotide polymorphisms was calculated using published effect estimates.

RESULTS

Comparing overweight (BMI ≥ 25 kg/m²) and normal weight (BMI < 25 kg/m²) WGHS women, each unit increase of TG-wGRS was associated with TG increases of 1.013% and 1.011%, respectively, and this differential association was significant (P_interaction = 0.014). Metaanalyses combining results for WGHS BMI with the 4 Scandinavian cohorts (INTER99, HEALTH2006, GLACIER, MDC) (total n = 40026) yielded a more significant interaction (P_interaction = 0.001). Similarly, we observed differential association of the TG-wGRS with TG (P_interaction = 0.006) in strata of WC (<80 cm vs ≥80 cm). Metaanalysis with 2 additional cohorts reporting WC (INTER99 and HEALTH2006) (total n = 27834) was significant with consistent effects (P_interaction = 0.006). We also observed highly significant interactions of the TG-wGRS across the strata of BMI with very large, medium, and small TG-rich lipoprotein subfractions measured by nuclear magnetic resonance spectroscopy (all P_interactions < 0.0001). The differential effects were strongest for very large TG-rich lipoprotein.

CONCLUSIONS

Our results support the original findings and suggest that obese individuals may be more susceptible to aggregated genetic risk associated with common TG-raising alleles, with effects accentuated in the large TG-rich lipoprotein subfraction.

New Rodent Population Models May Inform Human Health Risk Assessment and Identification of Genetic Susceptibility to Environmental Exposures

BACKGROUND

This paper provides an introduction for environmental health scientists to emerging population-based rodent resources. Mouse reference populations provide an opportunity to model environmental exposures and gene-environment interactions in human disease and to inform human health risk assessment.

OBJECTIVES

This review will describe several mouse populations for toxicity assessment, including older models such as the Mouse Diversity Panel (MDP), and newer models that include the Collaborative Cross (CC) and Diversity Outbred (DO) models.

METHODS

This review will outline the features of the MDP, CC, and DO mouse models and will discuss published case studies investigating the use of these mouse population resources in each step of the risk assessment paradigm.

DISCUSSION

These unique resources have the potential to be powerful tools for generating hypotheses related to gene-environment interplay in human disease, performing controlled exposure studies to understand the differential responses in humans for susceptibility or resistance to environmental exposures, and identifying gene variants that influence sensitivity to toxicity and disease states.

CONCLUSIONS

These new resources offer substantial advances to classical toxicity testing paradigms by including genetically sensitive individuals that may inform toxicity risks for sensitive subpopulations. Both in vivo and complementary in vitro resources provide platforms with which to reduce uncertainty by providing population-level data around biological variability. https://doi.org/10.1289/EHP1274.

The Genetic Architecture of Noise-Induced Hearing Loss: Evidence for a Gene-by-Environment Interaction

The discovery of environmentally specific genetic effects is crucial to the understanding of complex traits, such as susceptibility to noise-induced hearing loss (NIHL). We describe the first genome-wide association study (GWAS) for NIHL in a large and well-characterized population of inbred mouse strains, known as the Hybrid Mouse Diversity Panel (HMDP). We recorded auditory brainstem response (ABR) thresholds both pre and post 2-hr exposure to 10-kHz octave band noise at 108 dB sound pressure level in 5–6-wk-old female mice from the HMDP (4–5 mice/strain). From the observation that NIHL susceptibility varied among the strains, we performed a GWAS with correction for population structure and mapped a locus on chromosome 6 that was statistically significantly associated with two adjacent frequencies. We then used a “genetical genomics” approach that included the analysis of cochlear eQTLs to identify candidate genes within the GWAS QTL. In order to validate the gene-by-environment interaction, we compared the effects of the postnoise exposure locus with that from the same unexposed strains. The most significant SNP at chromosome 6 (rs37517079) was associated with noise susceptibility, but was not significant at the same frequencies in our unexposed study. These findings demonstrate that the genetic architecture of NIHL is distinct from that of unexposed hearing levels and provide strong evidence for gene-by-environment interactions in NIHL.

ForestPMPlot: A Flexible Tool for Visualizing Heterogeneity Between Studies in Meta-analysis

Meta-analysis has become a popular tool for genetic association studies to combine different genetic studies. A key challenge in meta-analysis is heterogeneity, or the differences in effect sizes between studies. Heterogeneity complicates the interpretation of meta-analyses. In this paper, we describe ForestPMPlot, a flexible visualization tool for analyzing studies included in a meta-analysis. The main feature of the tool is visualizing the differences in the effect sizes of the studies to understand why the studies exhibit heterogeneity for a particular phenotype and locus pair under different conditions. We show the application of this tool to interpret a meta-analysis of 17 mouse studies, and to interpret a multi-tissue eQTL study.

The Hybrid Mouse Diversity Panel: a resource for systems genetics analyses of metabolic and cardiovascular traits

The Hybrid Mouse Diversity Panel (HMDP) is a collection of approximately 100 well-characterized inbred strains of mice that can be used to analyze the genetic and environmental factors underlying complex traits. While not nearly as powerful for mapping genetic loci contributing to the traits as human genome-wide association studies, it has some important advantages. First, environmental factors can be controlled. Second, relevant tissues are accessible for global molecular phenotyping. Finally, because inbred strains are renewable, results from separate studies can be integrated. Thus far, the HMDP has been studied for traits relevant to obesity, diabetes, atherosclerosis, osteoporosis, heart failure, immune regulation, fatty liver disease, and host-gut microbiota interactions. High-throughput technologies have been used to examine the genomes, epigenomes, transcriptomes, proteomes, metabolomes, and microbiomes of the mice under various environmental conditions. All of the published data are available and can be readily used to formulate hypotheses about genes, pathways and interactions.

Collaborative Cross and Diversity Outbred data resources in the Mouse Phenome Database

The Mouse Phenome Database was originally conceived as a platform for the integration of phenotype data collected on a defined collection of 40 inbred mouse strains--the "phenome panel." This model provided an impetus for community data sharing, and integration was readily achieved through the reproducible genotypes of the phenome panel strains. Advances in the development of mouse populations lead to an expanded role of the Mouse Phenome Database to encompass new strain panels and inbred strain crosses. The recent introduction of the Collaborative Cross and Diversity Outbred mice, which share an extensive pool of genetic variation from eight founder inbred strains, presents new opportunities and challenges for community data resources. A wide variety of molecular and clinical phenotypes are being collected across genotypes, tissues, ages, environmental exposures, interventions, and treatments. The Mouse Phenome Database provides a framework for retrieval, integration, analysis, and display of these data, enabling them to be evaluated in the context of existing data from standard inbred strains. Primary data in the Mouse Phenome Database are supported by extensive metadata on protocols and procedures. These are centrally curated to ensure accuracy and reproducibility and to provide data in consistent formats. The Mouse Phenome Database represents an established and growing community data resource for mouse phenotype data and encourages submissions from new mouse resources, enabling investigators to integrate existing data into their studies of the phenotypic consequences of genetic variation.

GeneWeaver: finding consilience in heterogeneous cross-species functional genomics data

A persistent challenge lies in the interpretation of consensus and discord from functional genomics experimentation. Harmonizing and analyzing this data will enable investigators to discover relations of many genes to many diseases, and from many phenotypes and experimental paradigms to many diseases through their genomic substrates. The GeneWeaver.org system provides a platform for cross-species integration and interrogation of heterogeneous curated and experimentally derived functional genomics data. GeneWeaver enables researchers to store, share, analyze, and compare results of their own genome-wide functional genomics experiments in an environment containing rich companion data obtained from major curated repositories, including the Mouse Genome Database and other model organism databases, along with derived data from highly specialized resources, publications, and user submissions. The data, largely consisting of gene sets and putative biological networks, are mapped onto one another through gene identifiers and homology across species. A versatile suite of interactive tools enables investigators to perform a variety of set analysis operations to find consilience among these often noisy experimental results. Fast algorithms enable real-time analysis of large queries. Specific applications include prioritizing candidate genes for quantitative trait loci, identifying biologically valid mouse models and phenotypic assays for human disease, finding the common biological substrates of related diseases, classifying experiments and the biological concepts they represent from empirical data, and applying patterns of genomic evidence to implicate novel genes in disease. These results illustrate an alternative to strict emphasis on replicability, whereby researchers classify experimental results to identify the conditions that lead to their similarity.

Contrasting genetic architectures in different mouse reference populations used for studying complex traits

Quantitative trait loci (QTLs) are being used to study genetic networks, protein functions, and systems properties that underlie phenotypic variation and disease risk in humans, model organisms, agricultural species, and natural populations. The challenges are many, beginning with the seemingly simple tasks of mapping QTLs and identifying their underlying genetic determinants. Various specialized resources have been developed to study complex traits in many model organisms. In the mouse, remarkably different pictures of genetic architectures are emerging. Chromosome Substitution Strains (CSSs) reveal many QTLs, large phenotypic effects, pervasive epistasis, and readily identified genetic variants. In contrast, other resources as well as genome-wide association studies (GWAS) in humans and other species reveal genetic architectures dominated with a relatively modest number of QTLs that have small individual and combined phenotypic effects. These contrasting architectures are the result of intrinsic differences in the study designs underlying different resources. The CSSs examine context-dependent phenotypic effects independently among individual genotypes, whereas with GWAS and other mouse resources, the average effect of each QTL is assessed among many individuals with heterogeneous genetic backgrounds. We argue that variation of genetic architectures among individuals is as important as population averages. Each of these important resources has particular merits and specific applications for these individual and population perspectives. Collectively, these resources together with high-throughput genotyping, sequencing and genetic engineering technologies, and information repositories highlight the power of the mouse for genetic, functional, and systems studies of complex traits and disease models.

Gene-environment interaction in dyslipidemia

PURPOSE OF REVIEW

Recent genome-wide association studies have identified numerous common genetic variants associated with plasma lipid traits and have provided new insights into the regulation of lipoprotein metabolism including the identification of novel biological processes. These findings add to a body of existing data on dietary and environmental factors affecting plasma lipids. Here we explore how interactions between genetic risk factors and other phenotypes may explain some of the missing heritability of plasma lipid traits.

RECENT FINDINGS

Recent studies have identified true statistical interaction between several environmental and genetic risk factors and their effects on plasma lipid fractions. These include interactions between behaviors such as smoking or exercise as well as specific dietary nutrients and the effect size of specific genetic variants on plasma lipid traits risk and modifying effects of measures of adiposity on the cumulative impact of a number of common genetic variants on each of plasma triglycerides and HDL cholesterol.

SUMMARY

Interactions between genetic risk factors and clinical phenotypes may account for some of the unexplained heritability of plasma lipid traits. Recent studies provide biological insight into specific genetic associations and may aid in the identification of dyslipidemic patients for whom specific lifestyle interventions are likely to be most effective.

Trans-ethnic genome-wide association studies: advantages and challenges of mapping in diverse populations

Genome-wide association studies (GWASs) are the method most often used by geneticists to interrogate the human genome, and they provide a cost-effective way to identify the genetic variants underpinning complex traits and diseases. Most initial GWASs have focused on genetically homogeneous cohorts from European populations given the limited availability of ethnic minority samples and so as to limit population stratification effects. Transethnic studies have been invaluable in explaining the heritability of common quantitative traits, such as height, and in examining the genetic architecture of complex diseases, such as type 2 diabetes. They provide an opportunity for large-scale signal replication in independent populations and for cross-population meta-analyses to boost statistical power. In addition, transethnic GWASs enable prioritization of candidate genes, fine-mapping of functional variants, and potentially identification of SNPs associated with disease risk in admixed populations, by taking advantage of natural differences in genomic linkage disequilibrium across ethnically diverse populations. Recent efforts to assess the biological function of variants identified by GWAS have highlighted the need for large-scale replication, meta-analyses and fine-mapping across worldwide populations of ethnically diverse genetic ancestries. Here, we review recent advances and new approaches that are important to consider when performing, designing or interpreting transethnic GWASs, and we highlight existing challenges, such as the limited ability to handle heterogeneity in linkage disequilibrium across populations and limitations in dissecting complex architectures, such as those found in recently admixed populations.

Genome-wide association study for age-related hearing loss (AHL) in the mouse: a meta-analysis

Age-related hearing loss (AHL) is characterized by a symmetric sensorineural hearing loss primarily in high frequencies and individuals have different levels of susceptibility to AHL. Heritability studies have shown that the sources of this variance are both genetic and environmental, with approximately half of the variance attributable to hereditary factors as reported by Huag and Tang (Eur Arch Otorhinolaryngol 267(8):1179-1191, 2010). Only a limited number of large-scale association studies for AHL have been undertaken in humans, to date. An alternate and complementary approach to these human studies is through the use of mouse models. Advantages of mouse models include that the environment can be more carefully controlled, measurements can be replicated in genetically identical animals, and the proportion of the variability explained by genetic variation is increased. Complex traits in mouse strains have been shown to have higher heritability and genetic loci often have stronger effects on the trait compared to humans. Motivated by these advantages, we have performed the first genome-wide association study of its kind in the mouse by combining several data sets in a meta-analysis to identify loci associated with age-related hearing loss. We identified five genome-wide significant loci (<10(-6)). One of these loci confirmed a previously identified locus (ahl8) on distal chromosome 11 and greatly narrowed the candidate region. Specifically, the most significant associated SNP is located 450 kb upstream of Fscn2. These data confirm the utility of this approach and provide new high-resolution mapping information about variation within the mouse genome associated with hearing loss.

No boundaries: genomes, organisms, and ecological interactions responsible for divergence and reproductive isolation

Revealing the genetic basis of traits that cause reproductive isolation, particularly premating or sexual isolation, usually involves the same challenges as most attempts at genotype-phenotype mapping and so requires knowledge of how these traits are expressed in different individuals, populations, and environments, particularly under natural conditions. Genetic dissection of speciation phenotypes thus requires understanding of the internal and external contexts in which underlying genetic elements are expressed. Gene expression is a product of complex interacting factors internal and external to the organism including developmental programs, the genetic background including nuclear-cytotype interactions, epistatic relationships, interactions among individuals or social effects, stochasticity, and prevailing variation in ecological conditions. Understanding of genomic divergence associated with reproductive isolation will be facilitated by functional expression analysis of annotated genomes in organisms with well-studied evolutionary histories, phylogenetic affinities, and known patterns of ecological variation throughout their life cycles. I review progress and prospects for understanding the pervasive role of host plant use on genetic and phenotypic expression of reproductive isolating mechanisms in cactophilic Drosophila mojavensis and suggest how this system can be used as a model for revealing the genetic basis for species formation in organisms where speciation phenotypes are under the joint influences of genetic and environmental factors.