Genetic architecture of complex traits: large phenotypic effects and pervasive epistasis

Pleiotropy, epistasis and the genetic architecture of quantitative traits

Pleiotropy (whereby one genetic polymorphism affects multiple traits) and epistasis (whereby non-linear interactions between genetic polymorphisms affect the same trait) are fundamental aspects of the genetic architecture of quantitative traits. Recent advances in the ability to characterize the effects of polymorphic variants on molecular and organismal phenotypes in human and model organism populations have revealed the prevalence of pleiotropy and unexpected shared molecular genetic bases among quantitative traits, including diseases. By contrast, epistasis is common between polymorphic loci associated with quantitative traits in model organisms, such that alleles at one locus have different effects in different genetic backgrounds, but is rarely observed for human quantitative traits and common diseases. Here, we review the concepts and recent inferences about pleiotropy and epistasis, and discuss factors that contribute to similarities and differences between the genetic architecture of quantitative traits in model organisms and humans.

Factorizing polygenic epistasis improves prediction and uncovers biological pathways in complex traits

Epistasis is central in many domains of biology, but it has not yet been proven useful for understanding the etiology of complex traits. This is partly because complex-trait epistasis involves polygenic interactions that are poorly captured in current models. To address this gap, we developed a model called Epistasis Factor Analysis (EFA). EFA assumes that polygenic epistasis can be factorized into interactions between a few epistasis factors (EFs), which represent latent polygenic components of the observed complex trait. The statistical goals of EFA are to improve polygenic prediction and to increase power to detect epistasis, while the biological goal is to unravel genetic effects into more-homogeneous units. We mathematically characterize EFA and use simulations to show that EFA outperforms current epistasis models when its assumptions approximately hold. Applied to predicting yeast growth rates, EFA outperforms the additive model for several traits with large epistasis heritability and uniformly outperforms the standard epistasis model. We replicate these prediction improvements in a second dataset. We then apply EFA to four previously characterized traits in the UK Biobank and find statistically significant epistasis in all four, including two that are robust to scale transformation. Moreover, we find that the inferred EFs partly recover pre-defined biological pathways for two of the traits. Our results demonstrate that more realistic models can identify biologically and statistically meaningful epistasis in complex traits, indicating that epistasis has potential for precision medicine and characterizing the biology underlying GWAS results.

The genotype-phenotype relationship and evolutionary genetics in the light of the Metabolic Control Analysis

Metabolic control analysis has long been used as a systemic model of the genotype-phenotype (GP) relationship. By considering kinetic parameters and enzyme concentrations as reflecting the genotype level and metabolic fluxes or pools as phenotypes related to fitness, MCA has given a biological basis to the relationship between these two levels. The non-linear and concave relationship between enzymes and fluxes can account for common genetic effects that reductionist approaches have been powerless to explain, such as the dominance of active alleles over less active alleles, the various types of epistasis and heterosis, and reveals the structural links between these genetic effects. The summation property of the flux control coefficients accounts for the L-shaped distribution of Quantitative Trait Locus (QTL) effects, irrespective of other possible causes. Metabolic models of response to selection results in evolutionary scenarios that are markedly different from those derived from the classical infinitesimal model of quantitative genetics. In particular, evolution towards selective neutrality appears to be a consequence of the diminishing return of the flux-enzyme relationship. In this paper, we survey the historical and recent achievements of MCA in genetics, quantitative genetics and evolution, focusing on epistasis and the evolution of flux in relation to enzyme concentrations.

Leveraging the genetic correlation between traits improves the detection of epistasis in genome-wide association studies

Epistasis, commonly defined as the interaction between genetic loci, is known to play an important role in the phenotypic variation of complex traits. As a result, many statistical methods have been developed to identify genetic variants that are involved in epistasis, and nearly all of these approaches carry out this task by focusing on analyzing one trait at a time. Previous studies have shown that jointly modeling multiple phenotypes can often dramatically increase statistical power for association mapping. In this study, we present the "multivariate MArginal ePIstasis Test" (mvMAPIT)-a multioutcome generalization of a recently proposed epistatic detection method which seeks to detect marginal epistasis or the combined pairwise interaction effects between a given variant and all other variants. By searching for marginal epistatic effects, one can identify genetic variants that are involved in epistasis without the need to identify the exact partners with which the variants interact-thus, potentially alleviating much of the statistical and computational burden associated with conventional explicit search-based methods. Our proposed mvMAPIT builds upon this strategy by taking advantage of correlation structure between traits to improve the identification of variants involved in epistasis. We formulate mvMAPIT as a multivariate linear mixed model and develop a multitrait variance component estimation algorithm for efficient parameter inference and P-value computation. Together with reasonable model approximations, our proposed approach is scalable to moderately sized genome-wide association studies. With simulations, we illustrate the benefits of mvMAPIT over univariate (or single-trait) epistatic mapping strategies. We also apply mvMAPIT framework to protein sequence data from two broadly neutralizing anti-influenza antibodies and approximately 2,000 heterogeneous stock of mice from the Wellcome Trust Centre for Human Genetics. The mvMAPIT R package can be downloaded at https://github.com/lcrawlab/mvMAPIT.

Unveiling the Mysteries of Non-Mendelian Heredity in Plant Breeding

Mendelian heredity is the cornerstone of plant breeding and has been used to develop new varieties of plants since the 19th century. However, there are several breeding cases, such as cytoplasmic inheritance, methylation, epigenetics, hybrid vigor, and loss of heterozygosity (LOH), where Mendelian heredity is not applicable, known as non-Mendelian heredity. This type of inheritance can be influenced by several factors besides the genetic architecture of the plant and its breeding potential. Therefore, exploring various non-Mendelian heredity mechanisms, their prevalence in plants, and the implications for plant breeding is of paramount importance to accelerate the pace of crop improvement. In this review, we examine the current understanding of non-Mendelian heredity in plants, including the mechanisms, inheritance patterns, and applications in plant breeding, provide an overview of the various forms of non-Mendelian inheritance (including epigenetic inheritance, cytoplasmic inheritance, hybrid vigor, and LOH), explore insight into the implications of non-Mendelian heredity in plant breeding, and the potential it holds for future research.

Allograft inflammatory factor-1-like is a situational regulator of leptin levels, hyperphagia, and obesity

Mouse models enable the study of genetic factors affecting the complex pathophysiology of metabolic disorders. Here, we identify reductions in leptin levels, food intake, and obesity due to high-fat diet, accompanied by increased leptin sensitivity, in mice that harbor the E2a-Cre transgene within Obrq2, an obesity quantitative trait locus (QTL) that includes the leptin gene. Interestingly, loss of allograft inflammatory factor-1-like (AIF1L) protein in these transgenic mice leads to similar leptin sensitivity, yet marked reversal of the obesity phenotype, with accelerated weight gain and increased food intake. Transgenic mice lacking AIF1L also have low circulating leptin, which suggests that benefits of enhanced leptin sensitivity are lost with further impairment of leptin expression due to loss of AIF1L. Together, our results identify AIF1L as a genetic modifier of Obrq2 and leptin that affects leptin levels, food intake, and obesity during the metabolic stress imposed by HFD.

Quantitative trait loci, G×E and G×G for glycemic traits: response to metformin and placebo in the Diabetes Prevention Program (DPP)

The complex genetic architecture of type-2-diabetes (T2D) includes gene-by-environment (G×E) and gene-by-gene (G×G) interactions. To identify G×E and G×G, we screened markers for patterns indicative of interactions (relationship loci [rQTL] and variance heterogeneity loci [vQTL]). rQTL exist when the correlation between multiple traits varies by genotype and vQTL occur when the variance of a trait differs by genotype (potentially flagging G×G and G×E). In the metformin and placebo arms of the DPP (n = 1762) we screened 280,965 exomic and intergenic SNPs, for rQTL and vQTL patterns in association with year one changes from baseline in glycemia and related traits (insulinogenic index [IGI], insulin sensitivity index [ISI], fasting glucose and fasting insulin). Significant (p < 1.8 × 10^-7) rQTL and vQTL generated a priori hypotheses of individual G×E tests for a SNP × metformin treatment interaction and secondarily for G×G screens. Several rQTL and vQTL identified led to 6 nominally significant (p < 0.05) metformin treatment × SNP interactions (4 for IGI, one insulin, and one glucose) and 12G×G interactions (all IGI) that exceeded experiment-wide significance (p < 4.1 × 10^-9). Some loci are directly associated with incident diabetes, and others are rQTL and modify a trait's relationship with diabetes (2 diabetes/glucose, 2 diabetes/insulin, 1 diabetes/IGI). rs3197999, an ISI/insulin rQTL, is a possible gene damaging missense mutation in MST1, is associated with ulcerative colitis, sclerosing cholangitis, Crohn's disease, BMI and coronary artery disease. This study demonstrates evidence for context-dependent effects (G×G & G×E) and the complexity of these T2D-related traits.

Biomarker interaction selection and disease detection based on multivariate gain ratio

Background

Disease detection is an important aspect of biotherapy. With the development of biotechnology and computer technology, there are many methods to detect disease based on single biomarker. However, biomarker does not influence disease alone in some cases. It’s the interaction between biomarkers that determines disease status. The existing influence measure I-score is used to evaluate the importance of interaction in determining disease status, but there is a deviation about the number of variables in interaction when applying I-score. To solve the problem, we propose a new influence measure Multivariate Gain Ratio (MGR) based on Gain Ratio (GR) of single-variate, which provides us with multivariate combination called interaction.

Results

We propose a preprocessing verification algorithm based on partial predictor variables to select an appropriate preprocessing method. In this paper, an algorithm for selecting key interactions of biomarkers and applying key interactions to construct a disease detection model is provided. MGR is more credible than I-score in the case of interaction containing small number of variables. Our method behaves better with average accuracy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$93.13\%$$\end{document}93.13% than I-score of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$91.73\%$$\end{document}91.73% in Breast Cancer Wisconsin (Diagnostic) Dataset. Compared to the classification results \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$89.80\%$$\end{document}89.80% based on all predictor variables, MGR identifies the true main biomarkers and realizes the dimension reduction. In Leukemia Dataset, the experiment results show the effectiveness of MGR with the accuracy of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$97.32\%$$\end{document}97.32% compared to I-score with accuracy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$89.11\%$$\end{document}89.11%. The results can be explained by the nature of MGR and I-score mentioned above because every key interaction contains a small number of variables in Leukemia Dataset.

Conclusions

MGR is effective for selecting important biomarkers and biomarker interactions even in high-dimension feature space in which the interaction could contain more than two biomarkers. The prediction ability of interactions selected by MGR is better than I-score in the case of interaction containing small number of variables. MGR is generally applicable to various types of biomarker datasets including cell nuclei, gene, SNPs and protein datasets.

The interplay of additivity, dominance, and epistasis on fitness in a diploid yeast cross

In diploid species, genetic loci can show additive, dominance, and epistatic effects. To characterize the contributions of these different types of genetic effects to heritable traits, we use a double barcoding system to generate and phenotype a panel of ~200,000 diploid yeast strains that can be partitioned into hundreds of interrelated families. This experiment enables the detection of thousands of epistatic loci, many whose effects vary across families. Here, we show traits are largely specified by a small number of hub loci with major additive and dominance effects, and pervasive epistasis. Genetic background commonly influences both the additive and dominance effects of loci, with multiple modifiers typically involved. The most prominent dominance modifier in our data is the mating locus, which has no effect on its own. Our findings show that the interplay between additivity, dominance, and epistasis underlies a complex genotype-to-phenotype map in diploids.

Genetic basis of a spontaneous mutation’s expressivity

Genetic background often influences the phenotypic consequences of mutations, resulting in variable expressivity. How standing genetic variants collectively cause this phenomenon is not fully understood. Here, we comprehensively identify loci in a budding yeast cross that impact the growth of individuals carrying a spontaneous missense mutation in the nuclear-encoded mitochondrial ribosomal gene MRP20. Initial results suggested that a single large effect locus influences the mutation’s expressivity, with one allele causing inviability in mutants. However, further experiments revealed this simplicity was an illusion. In fact, many additional loci shape the mutation’s expressivity, collectively leading to a wide spectrum of mutational responses. These results exemplify how complex combinations of alleles can produce a diversity of qualitative and quantitative responses to the same mutation.

Natural genetic variation as a tool for discovery in Caenorhabditis nematodes

Over the last 20 years, studies of Caenorhabditis elegans natural diversity have demonstrated the power of quantitative genetic approaches to reveal the evolutionary, ecological, and genetic factors that shape traits. These studies complement the use of the laboratory-adapted strain N2 and enable additional discoveries not possible using only one genetic background. In this chapter, we describe how to perform quantitative genetic studies in Caenorhabditis, with an emphasis on C. elegans. These approaches use correlations between genotype and phenotype across populations of genetically diverse individuals to discover the genetic causes of phenotypic variation. We present methods that use linkage, near-isogenic lines, association, and bulk-segregant mapping, and we describe the advantages and disadvantages of each approach. The power of C. elegans quantitative genetic mapping is best shown in the ability to connect phenotypic differences to specific genes and variants. We will present methods to narrow genomic regions to candidate genes and then tests to identify the gene or variant involved in a quantitative trait. The same features that make C. elegans a preeminent experimental model animal contribute to its exceptional value as a tool to understand natural phenotypic variation.

Genetic prediction of complex traits with polygenic scores: a statistical review

Accurate genetic prediction of complex traits can facilitate disease screening, improve early intervention, and aid in the development of personalized medicine. Genetic prediction of complex traits requires the development of statistical methods that can properly model polygenic architecture and construct a polygenic score (PGS). We present a comprehensive review of 46 methods for PGS construction. We connect the majority of these methods through a multiple linear regression framework which can be instrumental for understanding their prediction performance for traits with distinct genetic architectures. We discuss the practical considerations of PGS analysis as well as challenges and future directions of PGS method development. We hope our review serves as a useful reference both for statistical geneticists who develop PGS methods and for data analysts who perform PGS analysis.

Capturing the rapidly evolving study of adaptation

Research on the genomics of adaptation is rapidly changing. In the last few decades, progress in this area has been driven by methodological advances, not only in the way increasingly large amounts of molecular data are generated (e.g. with high-throughput sequencing), but also in the way these data are analysed. This includes a growing appreciation and quantitative treatment of covariation among units within the same data type (e.g. genes) or across data types (e.g. genes and phenotypes). The development and adoption of more and more integrative tools have resulted in richer and more interesting empirical work. This special issue - comprising methodological, empirical, and review papers - aims to capture a 'snapshot' of this rapidly evolving field. We discuss in particular three important themes in the study of adaptation: the genetic architecture of adaptive variation, protein-coding and regulatory changes, and parallel evolution. We highlight how more traditional key themes in the study of genetic architecture (e.g. the number of loci underlying adaptive traits and the distribution of their effects) are now being complemented by other factors (e.g. how patterns of linkage and number of loci interact to affect the ability to adapt). Similarly, apart from addressing the relative importance of protein-coding and regulatory changes, we now have the tools to look in-depth at specific types of regulatory variation to gain a clearer picture of regulatory networks. Finally, parallel evolution has always been central to the study of adaptation, but now we are often able to address the question of whether - and to what extent - parallelism at the organismal or phenotypic level is matched by parallelism at the genetic level. Perhaps most importantly, we can now determine what mechanisms are driving parallelism (or lack thereof) across levels of biological organization. All these recent methodological developments open up new directions for future studies of adaptive changes across traits, levels of biological organization, demographic contexts and time scales.

A model and test for coordinated polygenic epistasis in complex traits

Interactions between genetic variants-epistasis-is pervasive in model systems and can profoundly impact evolutionary adaption, population disease dynamics, genetic mapping, and precision medicine efforts. In this work, we develop a model for structured polygenic epistasis, called coordinated epistasis (CE), and prove that several recent theories of genetic architecture fall under the formal umbrella of CE. Unlike standard epistasis models that assume epistasis and main effects are independent, CE captures systematic correlations between epistasis and main effects that result from pathway-level epistasis, on balance skewing the penetrance of genetic effects. To test for the existence of CE, we propose the even-odd (EO) test and prove it is calibrated in a range of realistic biological models. Applying the EO test in the UK Biobank, we find evidence of CE in 18 of 26 traits spanning disease, anthropometric, and blood categories. Finally, we extend the EO test to tissue-specific enrichment and identify several plausible tissue-trait pairs. Overall, CE is a dimension of genetic architecture that can capture structured, systemic forms of epistasis in complex human traits.

Hypothetical novel simulations to explain the evolutionary survival of the hypo-reproductive extreme tail in the complex human diversity

Hierarchical structures which lie hidden between human complex conditions and reproductivity cannot be simple, and trends of each population component does not necessarily pertain to evolutionary theories. As an illustration, the fitness of individuals with heritable extreme conditions can be low across continuing generations in observational data. Autism and schizophrenia are characterized by such evolutionary paradox of survival and hypo-reproductivity in the complex human diversity. Theoretical mechanisms for the observational fact were evaluated using a simple formula which was established to simulate stochastic epistasis-mediated phenotypic diversity. The survival of the hypo-reproductive extreme tail could be imitated just by the predominant presence of stochastic epistasis mechanism, suggesting that stochastic epistasis might be a genetic prerequisite for the evolutionary paradox. As supplemental cofactors of stochastic epistasis, a random link of the extreme tail to both un- and hyper-reproductivity and group assortative mating were shown to be effective for the paradox. Especially, the mixed localization of un- and hyper-reproductivity in the tail of a generational population evidently induced the continuous survival of outliers and extremes. These hypothetical considerations and mathematical simulations may suggest the significance of stochastic epistasis as the essential genetic background of complex human diversity.

Competitive history shapes rapid evolution in a seasonal climate

Eco-evolutionary dynamics will play a critical role in determining species' fates as climatic conditions change. Unfortunately, we have little understanding of how rapid evolutionary responses to climate play out when species are embedded in the competitive communities that they inhabit in nature. We tested the effects of rapid evolution in response to interspecific competition on subsequent ecological and evolutionary trajectories in a seasonally changing climate using a field-based evolution experiment with Drosophila melanogaster Populations of D. melanogaster were either exposed, or not exposed, to interspecific competition with an invasive competitor, Zaprionus indianus, over the summer. We then quantified these populations' ecological trajectories (abundances) and evolutionary trajectories (heritable phenotypic change) when exposed to a cooling fall climate. We found that competition with Z. indianus in the summer affected the subsequent evolutionary trajectory of D. melanogaster populations in the fall, after all interspecific competition had ceased. Specifically, flies with a history of interspecific competition evolved under fall conditions to be larger and have lower cold fecundity and faster development than flies without a history of interspecific competition. Surprisingly, this divergent fall evolutionary trajectory occurred in the absence of any detectible effect of the summer competitive environment on phenotypic evolution over the summer or population dynamics in the fall. This study demonstrates that competitive interactions can leave a legacy that shapes evolutionary responses to climate even after competition has ceased, and more broadly, that evolution in response to one selective pressure can fundamentally alter evolution in response to subsequent agents of selection.

The Combined Analysis of Pleiotropy and Epistasis (CAPE)

Epistasis, or gene-gene interaction, contributes substantially to trait variation in organisms ranging from yeast to humans, and modeling epistasis directly is critical to understanding the genotype-phenotype map. However, inference of genetic interactions is challenging compared to inference of individual allele effects due to low statistical power. Furthermore, genetic interactions can appear inconsistent across different quantitative traits, presenting a challenge for the interpretation of detected interactions. Here we present a method called the Combined Analysis of Pleiotropy and Epistasis (CAPE) that combines information across multiple quantitative traits to infer directed epistatic interactions. By combining information across multiple traits, CAPE not only increases power to detect genetic interactions but also interprets these interactions across traits to identify a single interaction that is consistent across all observed data. This method generates informative, interpretable interaction networks that explain how variants interact with each other to influence groups of related traits. This method could potentially be used to link genetic variants to gene expression, physiological endophenotypes, and higher-level disease traits.

The complex role of genetic background in shaping the effects of spontaneous and induced mutations

Spontaneous and induced mutations frequently show different phenotypic effects across genetically distinct individuals. It is generally appreciated that these background effects mainly result from genetic interactions between the mutations and segregating loci. However, the architectures and molecular bases of these genetic interactions are not well understood. Recent work in a number of model organisms has tried to advance knowledge of background effects both by using large‐scale screens to find mutations that exhibit this phenomenon and by identifying the specific loci that are involved. Here, we review this body of research, emphasizing in particular the insights it provides into both the prevalence of background effects across different mutations and the mechanisms that cause these background effects.

A large fraction of mutations show different effects in distinct individuals.

These background effects are mainly caused by epistasis with segregating loci.

Mapping studies show a diversity of genetic architectures can be involved.

Genetically complex changes in gene expression are often, but not always, causative.

A Novel Mapping Strategy Utilizing Mouse Chromosome Substitution Strains Identifies Multiple Epistatic Interactions That Regulate Complex Traits

The genetic contribution of additive vs. non-additive (epistatic) effects in the regulation of complex traits is unclear. While genome-wide association studies typically ignore gene-gene interactions, in part because of the lack of statistical power for detecting them, mouse chromosome substitution strains (CSSs) represent an alternate approach for detecting epistasis given their limited allelic variation. Therefore, we utilized CSSs to identify and map both additive and epistatic loci that regulate a range of hematologic- and metabolism-related traits, as well as hepatic gene expression. Quantitative trait loci (QTL) were identified using a CSS-based backcross strategy involving the segregation of variants on the A/J-derived substituted chromosomes 4 and 6 on an otherwise C57BL/6J genetic background. In the liver transcriptomes of offspring from this cross, we identified and mapped additive QTL regulating the hepatic expression of 768 genes, and epistatic QTL pairs for 519 genes. Similarly, we identified additive QTL for fat pad weight, platelets, and the percentage of granulocytes in blood, as well as epistatic QTL pairs controlling the percentage of lymphocytes in blood and red cell distribution width. The variance attributed to the epistatic QTL pairs was approximately equal to that of the additive QTL; however, the SNPs in the epistatic QTL pairs that accounted for the largest variances were undetected in our single locus association analyses. These findings highlight the need to account for epistasis in association studies, and more broadly demonstrate the importance of identifying genetic interactions to understand the complete genetic architecture of complex traits.

Natural variation in the sequestosome-related gene, sqst-5, underlies zinc homeostasis in Caenorhabditis elegans

Zinc is an essential trace element that acts as a co-factor for many enzymes and transcription factors required for cellular growth and development. Altering intracellular zinc levels can produce dramatic effects ranging from cell proliferation to cell death. To avoid such fates, cells have evolved mechanisms to handle both an excess and a deficiency of zinc. Zinc homeostasis is largely maintained via zinc transporters, permeable channels, and other zinc-binding proteins. Variation in these proteins might affect their ability to interact with zinc, leading to either increased sensitivity or resistance to natural zinc fluctuations in the environment. We can leverage the power of the roundworm nematode Caenorhabditis elegans as a tractable metazoan model for quantitative genetics to identify genes that could underlie variation in responses to zinc. We found that the laboratory-adapted strain (N2) is resistant and a natural isolate from Hawaii (CB4856) is sensitive to micromolar amounts of exogenous zinc supplementation. Using a panel of recombinant inbred lines, we identified two large-effect quantitative trait loci (QTL) on the left arm of chromosome III and the center of chromosome V that are associated with zinc responses. We validated and refined both QTL using near-isogenic lines (NILs) and identified a naturally occurring deletion in sqst-5, a sequestosome-related gene, that is associated with resistance to high exogenous zinc. We found that this deletion is relatively common across strains within the species and that variation in sqst-5 is associated with zinc resistance. Our results offer a possible mechanism for how organisms can respond to naturally high levels of zinc in the environment and how zinc homeostasis varies among individuals.

Dissecting the eQTL Micro-Architecture in Caenorhabditis elegans

The study of expression quantitative trait loci (eQTL) using natural variation in inbred populations has yielded detailed information about the transcriptional regulation of complex traits. Studies on eQTL using recombinant inbred lines (RILs) led to insights on cis and trans regulatory loci of transcript abundance. However, determining the underlying causal polymorphic genes or variants is difficult, but ultimately essential for the understanding of regulatory networks of complex traits. This requires insight into whether associated loci are single eQTL or a combination of closely linked eQTL, and how this QTL micro-architecture depends on the environment. We addressed these questions by testing for independent replication of previously mapped eQTL in Caenorhabditis elegans using new data from introgression lines (ILs). Both populations indicate that the overall heritability of gene expression, number, and position of eQTL differed among environments. Across environments we were able to replicate 70% of the cis- and 40% of the trans-eQTL using the ILs. Testing eight different simulation models, we suggest that additive effects explain up to 60-93% of RIL/IL heritability for all three environments. Closely linked eQTL explained up to 40% of RIL/IL heritability in the control environment whereas only 7% in the heat-stress and recovery environments. In conclusion, we show that reproducibility of eQTL was higher for cis vs. trans eQTL and that the environment affects the eQTL micro-architecture.

Flexible Mixture Model Approaches That Accommodate Footprint Size Variability for Robust Detection of Balancing Selection

Long-term balancing selection typically leaves narrow footprints of increased genetic diversity, and therefore most detection approaches only achieve optimal performances when sufficiently small genomic regions (i.e., windows) are examined. Such methods are sensitive to window sizes and suffer substantial losses in power when windows are large. Here, we employ mixture models to construct a set of five composite likelihood ratio test statistics, which we collectively term B statistics. These statistics are agnostic to window sizes and can operate on diverse forms of input data. Through simulations, we show that they exhibit comparable power to the best-performing current methods, and retain substantially high power regardless of window sizes. They also display considerable robustness to high mutation rates and uneven recombination landscapes, as well as an array of other common confounding scenarios. Moreover, we applied a specific version of the B statistics, termed B2, to a human population-genomic data set and recovered many top candidates from prior studies, including the then-uncharacterized STPG2 and CCDC169-SOHLH2, both of which are related to gamete functions. We further applied B2 on a bonobo population-genomic data set. In addition to the MHC-DQ genes, we uncovered several novel candidate genes, such as KLRD1, involved in viral defense, and SCN9A, associated with pain perception. Finally, we show that our methods can be extended to account for multiallelic balancing selection and integrated the set of statistics into open-source software named BalLeRMix for future applications by the scientific community.

The Effects of Quantitative Trait Architecture on Detection Power in Short-Term Artificial Selection Experiments

Evolve and resequence (E&R) experiments, in which artificial selection is imposed on organisms in a controlled environment, are becoming an increasingly accessible tool for studying the genetic basis of adaptation. Previous work has assessed how different experimental design parameters affect the power to detect the quantitative trait loci (QTL) that underlie adaptive responses in such experiments, but so far there has been little exploration of how this power varies with the genetic architecture of the evolving traits. In this study, we use forward simulation to build a more realistic model of an E&R experiment in which a quantitative polygenic trait experiences a short, but strong, episode of truncation selection. We study the expected power for QTL detection in such an experiment and how this power is influenced by different aspects of trait architecture, including the number of QTL affecting the trait, their starting frequencies, effect sizes, clustering along a chromosome, dominance, and epistasis patterns. We show that all of these parameters can affect allele frequency dynamics at the QTL and linked loci in complex and often unintuitive ways, and thus influence our power to detect them. One consequence of this is that existing detection methods based on models of independent selective sweeps at individual QTL often have lower detection power than a simple measurement of allele frequency differences before and after selection. Our findings highlight the importance of taking trait architecture into account when designing and interpreting studies of molecular adaptation with temporal data. We provide a customizable modeling framework that will enable researchers to easily simulate E&R experiments with different trait architectures and parameters tuned to their specific study system, allowing for assessment of expected detection power and optimization of experimental design.

Hidden genetic variance contributes to increase the short-term adaptive potential of selfing populations

Standing genetic variation is considered a major contributor to the adaptive potential of species. The low heritable genetic variation observed in self-fertilizing populations has led to the hypothesis that species with this mating system would be less likely to adapt. However, a non-negligible amount of cryptic genetic variation for polygenic traits, accumulated through negative linkage disequilibrium, could prove to be an important source of standing variation in self-fertilizing species. To test this hypothesis, we simulated populations under stabilizing selection subjected to an environmental change. We demonstrate that, when the mutation rate is high (but realistic), selfing populations are better able to store genetic variance than outcrossing populations through genetic associations, notably due to the reduced effective recombination rate associated with predominant selfing. Following an environmental shift, this diversity can be partially remobilized, which increases the additive variance and adaptive potential of predominantly (but not completely) selfing populations. In such conditions, despite initially lower observed genetic variance, selfing populations adapt as readily as outcrossing ones within a few generations. For low mutation rates, purifying selection impedes the storage of diversity through genetic associations, in which case, as previously predicted, the lower genetic variance of selfing populations results in lower adaptability compared to their outcrossing counterparts. The population size and the mutation rate are the main parameters to consider, as they are the best predictors of the amount of stored diversity in selfing populations. Our results and their impact on our knowledge of adaptation under high selfing rates are discussed.

Genetic dissection of the gene coexpression network underlying photosynthesis in Populus

Photosynthesis is a key reaction that ultimately generates the carbohydrates needed to form woody tissues in trees. However, the genetic regulatory network of protein-encoding genes (PEGs) and regulatory noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), underlying the photosynthetic pathway is unknown. Here, we integrated data from coexpression analysis, association studies (additive, dominance and epistasis), and expression quantitative trait nucleotide (eQTN) mapping to dissect the causal variants and genetic interaction network underlying photosynthesis in Populus. We initially used 30 PEGs, 6 miRNAs and 12 lncRNAs to construct a coexpression network based on the tissue-specific gene expression profiles of 15 Populus samples. Then, we performed association studies using a natural population of 435 unrelated Populus tomentosa individuals, and identified 72 significant associations (P ≤ 0.001, q ≤ 0.05) with diverse additive and dominance patterns underlying photosynthesis-related traits. Analysis of epistasis and eQTNs revealed that the complex genetic interactions in the coexpression network contribute to phenotypes at various levels. Finally, we demonstrated that heterologously expressing the most highly linked gene (PtoPsbX1) in this network significantly improved photosynthesis in Arabidopsis thaliana, pointing to the functional role of PtoPsbX1 in the photosynthetic pathway. This study provides an integrated strategy for dissecting a complex genetic interaction network, which should accelerate marker-assisted breeding efforts to genetically improve woody plants.

Evolution of Epistatic Networks and the Genetic Basis of Innate Behaviors

Instinctive behaviors are genetically programmed behaviors that occur independent of experience. How genetic programs that give rise to the manifestation of such behaviors evolve remains an unresolved question. I propose that evolution of species-specific innate behaviors is accomplished through progressive modifications of pre-existing genetic networks composed of allelic variants. I hypothesize that changes in frequencies of one or more constituent allelic variants within the network leads to changes in gene network connectivity and the emergence of a reorganized network that can support the emergence of a novel behavioral phenotype and becomes stabilized when key allelic variants are driven to fixation.

Tightly linked antagonistic-effect loci underlie polygenic phenotypic variation in C. elegans

Recent work has provided strong empirical support for the classic polygenic model for trait variation. Population-based findings suggest that most regions of genome harbor variation affecting most traits. Here, we use the approach of experimental genetics to show that, indeed, most genomic regions carry variants with detectable effects on growth and reproduction in Caenorhabditis elegans populations sensitized by nickel stress. Nine of 15 adjacent intervals on the X chromosome, each encompassing ∼0.001 of the genome, have significant effects when tested individually in near-isogenic lines (NILs). These intervals have effects that are similar in magnitude to those of genome-wide significant loci that we mapped in a panel of recombinant inbred advanced intercross lines (RIAILs). If NIL-like effects were randomly distributed across the genome, the RIAILs would exhibit phenotypic variance that far exceeds the observed variance. However, the NIL intervals are arranged in a pattern that significantly reduces phenotypic variance relative to a random arrangement; adjacent intervals antagonize one another, cancelling each other's effects. Contrary to the expectation of small additive effects, our findings point to large-effect variants whose effects are masked by epistasis or linkage disequilibrium between alleles of opposing effect.

Dissecting Epistatic QTL for Blood Pressure in Rats: Congenic Strains versus Heterogeneous Stocks, a Reality Check

Advances in molecular genetics have provided well-defined physical genetic maps and large numbers of genetic markers for both model organisms and humans. It is now possible to gain a fundamental understanding of the genetic architecture underlying quantitative traits, of which blood pressure (BP) is an important example. This review emphasizes analytical techniques and results obtained using the Dahl salt-sensitive (S) rat as a model of hypertension by presenting results in detail for three specific chromosomal regions harboring genetic elements of increasing complexity controlling BP. These results highlight the critical importance of genetic interactions (epistasis) on BP at all levels of structure, intragenic, intergenic, intrachromosomal, interchromosomal, and across whole genomes. In two of the three examples presented, specific DNA structural variations leading to biochemical, physiological, and pathological mechanisms are well defined. This proves the usefulness of the techniques involving interval mapping followed by substitution mapping using congenic strains. These classic techniques are compared to newer approaches using sophisticated statistical analysis on various segregating or outbred model-organism populations, which in some cases are uniquely useful in demonstrating the existence of higher-order interactions. It is speculated that hypertension as an outlier quantitative phenotype is dependent on higher-order genetic interactions. The obstacle to the identification of genetic elements and the biochemical/physiological mechanisms involved in higher-order interactions is not theoretical or technical but the lack of future resources to finish the job of identifying the individual genetic elements underlying the quantitative trait loci for BP and ascertaining their molecular functions. © 2019 American Physiological Society. Compr Physiol 9:1305-1337, 2019.

Genetically Distinct Behavioral Modules Underlie Natural Variation in Thermal Performance Curves

Thermal reaction norms pervade organismal traits as stereotyped responses to temperature, a fundamental environmental input into sensory and physiological systems. Locomotory behavior represents an especially plastic read-out of animal response, with its dynamic dependence on environmental stimuli presenting a challenge for analysis and for understanding the genomic architecture of heritable variation. Here we characterize behavioral reaction norms as thermal performance curves for the nematode Caenorhabditis briggsae, using a collection of 23 wild isolate genotypes and 153 recombinant inbred lines to quantify the extent of genetic and plastic variation in locomotory behavior to temperature changes. By reducing the dimensionality of the multivariate phenotypic response with a function-valued trait framework, we identified genetically distinct behavioral modules that contribute to the heritable variation in the emergent overall behavioral thermal performance curve. Quantitative trait locus mapping isolated regions on Chromosome II associated with locomotory activity at benign temperatures and Chromosome V loci related to distinct aspects of sensitivity to high temperatures, with each quantitative trait locus explaining up to 28% of trait variation. These findings highlight how behavioral responses to environmental inputs as thermal reaction norms can evolve through independent changes to genetically distinct modular components of such complex phenotypes.

Low Citrate Synthase Activity Is Associated with Glucose Intolerance and Lipotoxicity

Citrate synthase (CS) is a key mitochondrial enzyme. The aim of this study was to test the hypothesis that low CS activity impairs the metabolic health of mice fed a high fat diet (HFD) and promotes palmitate-induced lipotoxicity in muscle cells. C57BL/6J (B6) mice and congenic B6.A-(rs3676616-D10Utsw1)/KjnB6 (B6.A), a strain which carries the A/J allele of CS on the B6 strain background, were fed HFD (45% kcal from fat) for 12 weeks. C2C12 mouse muscle cells were used to investigate effects of CS knockdown on cell viability and signalling after incubation in 0.8 mM palmitate. CS activity, but not that of β-hydroxyacyl-coenzyme-A dehydrogenase was lower in the gastrocnemius muscle and heart of B6.A mice compared to B6 mice (P < 0.001). During HFD feeding, glucose tolerance of mice decreased progressively and to a greater extent in B6.A females compared to B6 females, with males showing a similar trend. Body weight and fat gain did not differ between B6.A and B6 mice. After an 18 h incubation in 0.8 mM palmitate C2C12 muscle cells with ∼50% shRNA mediated reduction in CS activity showed lower (P < 0.001) viability and increased (P < 0.001) levels of cleaved caspase-3 compared to the scramble shRNA treated C2C12 cells. A/J strain variant of CS is associated with low enzyme activity and impaired metabolic health. This could be due to impaired lipid metabolism in muscle cells.

Decoupling the Variances of Heterosis and Inbreeding Effects Is Evidenced in Yeast's Life-History and Proteomic Traits

Heterosis (hybrid vigor) and inbreeding depression, commonly considered as corollary phenomena, could nevertheless be decoupled under certain assumptions according to theoretical population genetics works. To explore this issue on real data, we analyzed the components of genetic variation in a population derived from a half-diallel cross between strains from Saccharomyces cerevisiae and S. uvarum, two related yeast species involved in alcoholic fermentation. A large number of phenotypic traits, either molecular (coming from quantitative proteomics) or related to fermentation and life history, were measured during alcoholic fermentation. Because the parental strains were included in the design, we were able to distinguish between inbreeding effects, which measure phenotypic differences between inbred and hybrids, and heterosis, which measures phenotypic differences between a specific hybrid and the other hybrids sharing a common parent. The sources of phenotypic variation differed depending on the temperature, indicating the predominance of genotype-by-environment interactions. Decomposing the total genetic variance into variances of additive (intra- and interspecific) effects, of inbreeding effects, and of heterosis (intra- and interspecific) effects, we showed that the distribution of variance components defined clear-cut groups of proteins and traits. Moreover, it was possible to cluster fermentation and life-history traits into most proteomic groups. Within groups, we observed positive, negative, or null correlations between the variances of heterosis and inbreeding effects. To our knowledge, such a decoupling had never been experimentally demonstrated. This result suggests that, despite a common evolutionary history of individuals within a species, the different types of traits have been subject to different selective pressures.

Large-effect loci affect survival in Tasmanian devils (Sarcophilus harrisii) infected with a transmissible cancer

Identifying the genetic architecture of complex phenotypes is a central goal of modern biology, particularly for disease-related traits. Genome-wide association methods are a classical approach for identifying the genomic basis of variation in disease phenotypes, but such analyses are particularly challenging in natural populations due to sample size difficulties. Extensive mark-recapture data, strong linkage disequilibrium and a lethal transmissible cancer make the Tasmanian devil (Sarcophilus harrisii) an ideal model for such an association study. We used a RAD-capture approach to genotype 624 devils at ~16,000 loci and then used association analyses to assess the heritability of three cancer-related phenotypes: infection case-control (where cases were infected devils and controls were devils that were never infected), age of first infection and survival following infection. The SNP array explained much of the phenotypic variance for female survival (>80%) and female case-control (>61%). We found that a few large-effect SNPs explained much of the variance for female survival (~5 SNPs explained >61% of the total variance), whereas more SNPs (~56) of smaller effect explained less of the variance for female case-control (~23% of the total variance). By contrast, these same SNPs did not account for a significant proportion of phenotypic variance in males, suggesting that the genetic bases of these traits and/or selection differ across sexes. Loci involved with cell adhesion and cell-cycle regulation underlay trait variation, suggesting that the devil immune system is rapidly evolving to recognize and potentially suppress cancer growth through these pathways. Overall, our study provided necessary data for genomics-based conservation and management in Tasmanian devils.

Synergistic interaction between APOE and family history of Alzheimer's disease on cerebral amyloid deposition and glucose metabolism

Background

Recently, the field of gene-gene or gene-environment interaction research appears to have gained growing interest, although it is seldom investigated in Alzheimer’s disease (AD). Hence, the current study aims to investigate interaction effects of the key genetic and environmental risks—the apolipoprotein ε4 allele (APOE4) and family history of late-onset AD (FH)—on AD-related brain changes in cognitively normal (CN) middle-aged and older adults.

Methods

[¹¹C] Pittsburg compound-B (PiB) positron emission tomography (PET) imaging as well as [¹⁸F] fluoro-2-deoxyglucose (FDG) PET that were simultaneously taken with T1-weighted magnetic resonance imaging (MRI) were obtained from 268 CNs from the Korean Brain Aging Study for Early Diagnosis and Prediction of AD (KBASE). Composite standardized uptake value ratios were obtained from PiB-PET and FDG-PET images in the AD signature regions of interests (ROIs) and analyzed. Voxel-wise analyses were also performed to examine detailed regional changes not captured by the ROI analyses.

Results

A significant synergistic interaction effect was found between the APOE4 and FH on amyloid-beta (Aβ) deposition in the AD signature ROIs as well as other regions. Synergistic interaction effects on cerebral glucose metabolism were observed in the regions not captured by the AD signature ROIs, particularly in the medial temporal regions.

Conclusions

Strong synergistic effects of APOE4 and FH on Aβ deposition and cerebral glucose metabolism in CN adults indicate possible gene-to-gene or gene-to-environment interactions that are crucial for pathogenesis of AD involving Aβ. Other unspecified risk factors—genes and/or environmental—that are captured by the positive FH status might either coexpress or interact with APOE4 to alter AD-related brain changes in CN. Healthy people with both FH and APOE4 need more attention for AD prevention.

Electronic supplementary material

The online version of this article (10.1186/s13195-018-0411-x) contains supplementary material, which is available to authorized users.

WISH-R- a fast and efficient tool for construction of epistatic networks for complex traits and diseases

Background

Genetic epistasis is an often-overlooked area in the study of the genomics of complex traits. Genome-wide association studies are a useful tool for revealing potential causal genetic variants, but in this context, epistasis is generally ignored. Data complexity and interpretation issues make it difficult to process and interpret epistasis. As the number of interaction grows exponentially with the number of variants, computational limitation is a bottleneck. Gene Network based strategies have been successful in integrating biological data and identifying relevant hub genes and pathways related to complex traits. In this study, epistatic interactions and network-based analysis are combined in the Weighted Interaction SNP hub (WISH) method and implemented in an efficient and easy to use R package.

Results

The WISH R package (WISH-R) was developed to calculate epistatic interactions on a genome-wide level based on genomic data. It is easy to use and install, and works on regular genomic data. The package filters data based on linkage disequilibrium and calculates epistatic interaction coefficients between SNP pairs based on a parallelized efficient linear model and generalized linear model implementations. Normalized epistatic coefficients are analyzed in a network framework, alleviating multiple testing issues and integrating biological signal to identify modules and pathways related to complex traits. Functions for visualizing results and testing runtimes are also provided.

Conclusion

The WISH-R package is an efficient implementation for analyzing genome-wide epistasis for complex diseases and traits. It includes methods and strategies for analyzing epistasis from initial data filtering until final data interpretation. WISH offers a new way to analyze genomic data by combining epistasis and network based analysis in one method and provides options for visualizations. This alleviates many of the existing hurdles in the analysis of genomic interactions.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2291-2) contains supplementary material, which is available to authorized users.

Genetic background effects in quantitative genetics: gene-by-system interactions

Proper cell function depends on networks of proteins that interact physically and functionally to carry out physiological processes. Thus, it seems logical that the impact of sequence variation in one protein could be significantly influenced by genetic variants at other loci in a genome. Nonetheless, the importance of such genetic interactions, known as epistasis, in explaining phenotypic variation remains a matter of debate in genetics. Recent work from our lab revealed that genes implicated from an association study of toxin tolerance in Saccharomyces cerevisiae show extensive interactions with the genetic background: most implicated genes, regardless of allele, are important for toxin tolerance in only one of two tested strains. The prevalence of background effects in our study adds to other reports of widespread genetic-background interactions in model organisms. We suggest that these effects represent many-way interactions with myriad features of the cellular system that vary across classes of individuals. Such gene-by-system interactions may influence diverse traits and require new modeling approaches to accurately represent genotype-phenotype relationships across individuals.

Fine mapping of epistatic genetic determinants of blood pressure on rat chromosome 5

BACKGROUND

Previously, using linkage analysis and substitution mapping, two closely-linked interactive blood pressure quantitative trait loci (QTLs), BP QTL1 and BP QTL2, were located within a 13.96 Mb region from 117894038 to 131853815 bp (RGSC 3.4 version) on rat chromosome 5 (RNO5). This was done by using a series of congenic strains consisting of genomic segments of the Dahl salt-sensitive (S) rat substituted with that of the normotensive Lewis (LEW) rat. The interactive nature of the two loci was further confirmed by the construction and characterization of a panel of S.LEW bicongenic strains and corresponding S.LEW monocongenic strains, which provided definitive evidence of epistasis (genetic interaction) between BP QTL1 (7.77 Mb) and BP QTL2 (4.18 Mb). The purpose of this work was to further map these interacting QTLs.

METHOD

A new panel of seven new S.LEW bicongenic strains was constructed and characterized for BP.

RESULTS

The data obtained from these new strains further resolved BP QTL1 from 7.77 to 2.93 Mb. Further, BP QTL2 was traceable as not being a single QTL, but a composite of at least three QTLs, LEW alleles at two of which located within 2.26 Mb and 175 kb lowered BP but the third one located within 1.31 Mb increased BP.

CONCLUSION

Lack of coding variation within any of the regions further mapped within the previous QTL2 suggests noncoding variation as likely responsible for the observed epistasis.

Dynamic and epistatic QTL mapping reveals the complex genetic architecture of waterlogging tolerance in chrysanthemum

37 unconditional QTLs, 51 conditional QTLs and considerable epistatic QTLs were detected for waterlogging tolerance, and six favourable combinations were selected accelerating the possible application of MAS in chrysanthemum breeding. Chrysanthemum is seriously impacted by soil waterlogging. To determine the genetic characteristics of waterlogging tolerance (WAT) in chrysanthemum, a population of 162 F₁ lines was used to construct a genetic map to identify the dynamic and epistatic quantitative trait loci (QTLs) for four WAT traits: wilting index (WI), dead leaf ratio (DLR), chlorosis score (Score) and membership function value of waterlogging (MFVW). The h _B² for the WAT traits ranged from 0.49 to 0.64, and transgressive segregation was observed in both directions. A total of 37 unconditional consensus QTLs with 5.81-18.21% phenotypic variation explanation (PVE) and 51 conditional consensus QTLs with 5.90-24.56% PVE were detected. Interestingly, three unconditional consensus QTLs were consistently identified across different stages, whereas no conditional consensus QTLs were consistently expressed. In addition, considerable epistatic QTLs, all with PVE values ranging from 0.01 to 8.87%, were detected by a joint analysis of WAT phenotypes. These results illustrated that the QTLs (genes) controlling WAT were environmentally dependent and selectively expressed at different times and indicated that both additive and epistatic effects underlie the inheritance of WAT in chrysanthemum. The findings of the current study provide insights into the complex genetic architecture of WAT, and the identification of favourable alleles represents an important step towards the application of molecular marker-assisted selection (MAS) and QTL pyramiding in chrysanthemum WAT breeding programmes.

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.

C57BL/6J mice as a polygenic developmental model of diet-induced obesity

Susceptibility to obesity changes during the course of life. We utilized the C57BL/6J (B6) and 129S mouse as a genetic model for variation in diet‐induced obesity to define the adiposity phenotypes from birth to maturity at 8 weeks‐of‐age. From birth to 8 weeks‐of‐age, both male and female 129S mice had significantly higher fat mass and adiposity index than B6 mice, although they were not obese. After 8 weeks‐of‐age, B6 had greater adiposity/obesity than 129S mice in response to a high fat (HF). We sought to determine the mechanism activating the fat accumulation in B6 mice at 8‐weeks‐of‐age. We used microarray analysis of gene expression during development of inguinal fat to show that molecular networks of lipogenesis were maximally expressed at 8 weeks‐of‐age. In addition, the DNA methylation analysis of the Sfrp5 promoter and binding of acetylated histones to Sfrp5 and Acly promoter regions showed that major differences in the expression of genes of lipogenesis and chromatin structure occur during development. Differences in lipogenesis networks could account for the strain‐dependent differences in adiposity up to 8 weeks‐of‐age; however, changes in the expression of genes in these networks were not associated with the susceptibility to DIO in B6 male mice beyond 8 weeks‐of‐age.

Polygenicity and Epistasis Underlie Fitness-Proximal Traits in the Caenorhabditis elegans Multiparental Experimental Evolution (CeMEE) Panel

Understanding the genetic basis of complex traits remains a major challenge in biology. Polygenicity, phenotypic plasticity, and epistasis contribute to phenotypic variance in ways that are rarely clear. This uncertainty can be problematic for estimating heritability, for predicting individual phenotypes from genomic data, and for parameterizing models of phenotypic evolution. Here, we report an advanced recombinant inbred line (RIL) quantitative trait locus mapping panel for the hermaphroditic nematode Caenorhabditis elegans, the C. elegans multiparental experimental evolution (CeMEE) panel. The CeMEE panel, comprising 507 RILs at present, was created by hybridization of 16 wild isolates, experimental evolution for 140-190 generations, and inbreeding by selfing for 13-16 generations. The panel contains 22% of single-nucleotide polymorphisms known to segregate in natural populations, and complements existing C. elegans mapping resources by providing fine resolution and high nucleotide diversity across > 95% of the genome. We apply it to study the genetic basis of two fitness components, fertility and hermaphrodite body size at time of reproduction, with high broad-sense heritability in the CeMEE. While simulations show that we should detect common alleles with additive effects as small as 5%, at gene-level resolution, the genetic architectures of these traits do not feature such alleles. We instead find that a significant fraction of trait variance, approaching 40% for fertility, can be explained by sign epistasis with main effects below the detection limit. In congruence, phenotype prediction from genomic similarity, while generally poor ([Formula: see text]), requires modeling epistasis for optimal accuracy, with most variance attributed to the rapidly evolving chromosome arms.

Genetic Dissection of Trabecular Bone Structure with Mouse Intersubspecific Consomic Strains

Trabecular bone structure has an important influence on bone strength, but little is known about its genetic regulation. To elucidate the genetic factor(s) regulating trabecular bone structure, we compared the trabecular bone structures of two genetically remote mouse strains, C57BL/6J and Japanese wild mouse-derived MSM/Ms. Phenotyping by X-ray micro-CT revealed that MSM/Ms has structurally more fragile trabecular bone than C57BL/6J. Toward identification of genetic determinants for the difference in fragility of trabecular bone between the two mouse strains, we employed phenotype screening of consomic mouse strains in which each C57BL/6J chromosome is substituted by its counterpart from MSM/Ms. The results showed that many chromosomes affect trabecular bone structure, and that the consomic strain B6-Chr15^MSM, carrying MSM/Ms-derived chromosome 15 (Chr15), has the lowest values for the parameters BV/TV, Tb.N, and Conn.D, and the highest values for the parameters Tb.Sp and SMI. Subsequent phenotyping of subconsomic strains for Chr15 mapped four novel trabecular bone structure-related QTL (Tbsq1-4) on mouse Chr15. These results collectively indicate that genetic regulation of trabecular bone structure is highly complex, and that even in the single Chr15, the combined action of the four Tbsqs controls the fragility of trabecular bone. Given that Tbsq4 is syntenic to human Chr 12q12-13.3, where several bone-related SNPs are assigned, further study of Tbsq4 should facilitate our understanding of the genetic regulation of bone formation in humans.

Widespread epistasis regulates glucose homeostasis and gene expression

The relative contributions of additive versus non-additive interactions in the regulation of complex traits remains controversial. This may be in part because large-scale epistasis has traditionally been difficult to detect in complex, multi-cellular organisms. We hypothesized that it would be easier to detect interactions using mouse chromosome substitution strains that simultaneously incorporate allelic variation in many genes on a controlled genetic background. Analyzing metabolic traits and gene expression levels in the offspring of a series of crosses between mouse chromosome substitution strains demonstrated that inter-chromosomal epistasis was a dominant feature of these complex traits. Epistasis typically accounted for a larger proportion of the heritable effects than those due solely to additive effects. These epistatic interactions typically resulted in trait values returning to the levels of the parental CSS host strain. Due to the large epistatic effects, analyses that did not account for interactions consistently underestimated the true effect sizes due to allelic variation or failed to detect the loci controlling trait variation. These studies demonstrate that epistatic interactions are a common feature of complex traits and thus identifying these interactions is key to understanding their genetic regulation.

Most complex traits and diseases are regulated by the combined influence of multiple genetic variants. However, it remains controversial whether these genetic variants independently influence complex traits, and therefore the impact of each variant could be simply added together (additivity), or whether the variants work together to influence trait variation, in which case the combined impact of multiple variants would differ from the summed impact of each individual variant (epistasis). In this study in mice, we discovered that the genetic regulation of blood sugar levels and gene expression in the liver were predominantly controlled by non-additive interactions, whereas body weight was predominantly controlled by additive interactions. Remarkably, the expression level of nearly 25% of all genes in the liver was controlled by non-additive interactions. The non-additive interactions typically acted to return trait values to the levels detected in control mice, thus contributing to a reduction in trait variation. We also demonstrated that not accounting for non-additive interactions significantly underestimated the phenotypic effect of a genetic variant on a particular genetic background, suggesting that many previously identified risk loci may have significantly larger effects on disease susceptibility in a subset of individuals. These studies highlight the importance of understanding interactions between genetic variants to better understand disease risk and personalize clinical care.

Comparative Analyses of Phenotypic Trait Covariation within and among Populations

Many morphological, behavioral, physiological, and life-history traits covary across the biological scales of individuals, populations, and species. However, the processes that cause traits to covary also change over these scales, challenging our ability to use patterns of trait covariance to infer process. Trait relationships are also widely assumed to have generic functional relationships with similar evolutionary potentials, and even though many different trait relationships are now identified, there is little appreciation that these may influence trait covariation and evolution in unique ways. We use a trait-performance-fitness framework to classify and organize trait relationships into three general classes, address which ones more likely generate trait covariation among individuals in a population, and review how selection shapes phenotypic covariation. We generate predictions about how trait covariance changes within and among populations as a result of trait relationships and in response to selection and consider how these can be tested with comparative data. Careful comparisons of covariation patterns can narrow the set of hypothesized processes that cause trait covariation when the form of the trait relationship and how it responds to selection yield clear predictions about patterns of trait covariation. We discuss the opportunities and limitations of comparative approaches to evaluate hypotheses about the evolutionary causes and consequences of trait covariation and highlight the importance of evaluating patterns within populations replicated in the same and in different selective environments. Explicit hypotheses about trait relationships are key to generating effective predictions about phenotype and its evolution using covariance data.

From Peas to Disease: Modifier Genes, Network Resilience, and the Genetics of Health

Phenotypes are rarely consistent across genetic backgrounds and environments, but instead vary in many ways depending on allelic variants, unlinked genes, epigenetic factors, and environmental exposures. In the extreme, individuals carrying the same causal DNA sequence variant but on different backgrounds can be classified as having distinct conditions. Similarly, some individuals that carry disease alleles are nevertheless healthy despite affected family members in the same environment. These genetic background effects often result from the action of so-called "modifier genes" that modulate the phenotypic manifestation of target genes in an epistatic manner. While complicating the prospects for gene discovery and the feasibility of mechanistic studies, such effects are opportunities to gain a deeper understanding of gene interaction networks that provide organismal form and function as well as resilience to perturbation. Here, we review the principles of modifier genetics and assess progress in studies of modifier genes and their targets in both simple and complex traits. We propose that modifier effects emerge from gene interaction networks whose structure and function vary with genetic background and argue that these effects can be exploited as safe and effective ways to prevent, stabilize, and reverse disease and dysfunction.

Epistatic Networks Jointly Influence Phenotypes Related to Metabolic Disease and Gene Expression in Diversity Outbred Mice

In this study, Tyler et al. analyzed the complex genetic architecture of metabolic disease-related traits using the Diversity Outbred mouse population

Genetic studies of multidimensional phenotypes can potentially link genetic variation, gene expression, and physiological data to create multi-scale models of complex traits. The challenge of reducing these data to specific hypotheses has become increasingly acute with the advent of genome-scale data resources. Multi-parent populations derived from model organisms provide a resource for developing methods to understand this complexity. In this study, we simultaneously modeled body composition, serum biomarkers, and liver transcript abundances from 474 Diversity Outbred mice. This population contained both sexes and two dietary cohorts. Transcript data were reduced to functional gene modules with weighted gene coexpression network analysis (WGCNA), which were used as summary phenotypes representing enriched biological processes. These module phenotypes were jointly analyzed with body composition and serum biomarkers in a combined analysis of pleiotropy and epistasis (CAPE), which inferred networks of epistatic interactions between quantitative trait loci that affect one or more traits. This network frequently mapped interactions between alleles of different ancestries, providing evidence of both genetic synergy and redundancy between haplotypes. Furthermore, a number of loci interacted with sex and diet to yield sex-specific genetic effects and alleles that potentially protect individuals from the effects of a high-fat diet. Although the epistatic interactions explained small amounts of trait variance, the combination of directional interactions, allelic specificity, and high genomic resolution provided context to generate hypotheses for the roles of specific genes in complex traits. Our approach moves beyond the cataloging of single loci to infer genetic networks that map genetic etiology by simultaneously modeling all phenotypes.

Detecting epistasis with the marginal epistasis test in genetic mapping studies of quantitative traits

Epistasis, commonly defined as the interaction between multiple genes, is an important genetic component underlying phenotypic variation. Many statistical methods have been developed to model and identify epistatic interactions between genetic variants. However, because of the large combinatorial search space of interactions, most epistasis mapping methods face enormous computational challenges and often suffer from low statistical power due to multiple test correction. Here, we present a novel, alternative strategy for mapping epistasis: instead of directly identifying individual pairwise or higher-order interactions, we focus on mapping variants that have non-zero marginal epistatic effects-the combined pairwise interaction effects between a given variant and all other variants. By testing marginal epistatic effects, we can identify candidate variants that are involved in epistasis without the need to identify the exact partners with which the variants interact, thus potentially alleviating much of the statistical and computational burden associated with standard epistatic mapping procedures. Our method is based on a variance component model, and relies on a recently developed variance component estimation method for efficient parameter inference and p-value computation. We refer to our method as the "MArginal ePIstasis Test", or MAPIT. With simulations, we show how MAPIT can be used to estimate and test marginal epistatic effects, produce calibrated test statistics under the null, and facilitate the detection of pairwise epistatic interactions. We further illustrate the benefits of MAPIT in a QTL mapping study by analyzing the gene expression data of over 400 individuals from the GEUVADIS consortium.