Genotypic Context and Epistasis in Individuals and Populations

Unraveling schizophrenia's genetic complexity through advanced causal inference and chromatin 3D conformation

Schizophrenia is a polygenic complex disease with a heritability as high as 80 %, yet the mechanism of polygenic interaction in its pathogenesis remains unclear. Studying the interaction and regulation of schizophrenia susceptibility genes is crucial for unraveling the pathogenesis of schizophrenia and developing antipsychotic drugs. Therefore, we developed a bioinformatics method named GRACI (Gene Regulation Analysis based on Causal Inference) based on the principles of information theory, a causal inference model, and high order chromatin 3D conformation. GRACI captures the interaction and regulatory relationships between schizophrenia susceptibility genes by analyzing genotyping data. Two datasets, comprising 1459 and 2065 samples respectively, were analyzed, and the gene networks from both datasets were constructed. GRACI showcased superior accuracy when compared to widely adopted methods for detecting gene-gene interactions and intergenic regulation. This alignment was further substantiated by its correlation with chromatin high-order conformation patterns. Using GRACI, we identified three potential genes-KCNN3, KCNH1, and KCND3-that are directly associated with schizophrenia pathogenesis. Furthermore, the results of GRACI on the standalone dataset illustrated the method's applicability to other complex diseases. GRACI download: https://github.com/liuliangjie19/GRACI.

Environment by environment interactions (ExE) differ across genetic backgrounds (ExExG)

While the terms "gene-by-gene interaction" (GxG) and "gene-by-environment interaction" (GxE) are commonplace within the fields of quantitative and evolutionary genetics, "environment-by-environment interaction" (ExE) is a term used less often. In this study, we find that environment-by-environment interactions are a meaningful driver of phenotypes, and that they differ across different genotypes (suggestive of ExExG). To reach this conclusion, we analyzed a large dataset of roughly 1,000 mutant yeast strains with varying degrees of resistance to different antifungal drugs. We show that the effectiveness of a drug combination, relative to single drugs, often varies across different drug resistant mutants. Even mutants that differ by only a single nucleotide change can have dramatically different drug x drug (ExE) interactions. We also introduce a new framework that better predicts the direction and magnitude of ExE interactions for some mutants. Studying how ExE interactions change across genotypes (ExExG) is not only important when modeling the evolution of pathogenic microbes, but also for broader efforts to understand the cell biology underlying these interactions and to resolve the source of phenotypic variance across populations. The relevance of ExExG interactions have been largely omitted from canon in evolutionary and population genetics, but these fields and others stand to benefit from perspectives that highlight how interactions between external forces craft the complex behavior of living systems.

Tissue-specific metabolic enzyme levels covary with whole-animal metabolic rates and life-history loci via epistatic effects

Metabolic rates, including standard (SMR) and maximum (MMR) metabolic rate have often been linked with life-history strategies. Variation in context- and tissue-level metabolism underlying SMR and MMR may thus provide a physiological basis for life-history variation. This raises a hypothesis that tissue-specific metabolism covaries with whole-animal metabolic rates and is genetically linked to life history. In Atlantic salmon (Salmo salar), variation in two loci, vgll3 and six6, affects life history via age-at-maturity as well as MMR. Here, using individuals with known SMR and MMR with different vgll3 and six6 genotype combinations, we measured proxies of mitochondrial density and anaerobic metabolism, i.e. maximal activities of the mitochondrial citrate synthase (CS) and lactate dehydrogenase (LDH) enzymes, in four tissues (heart, intestine, liver, white muscle) across low- and high-food regimes. We found enzymatic activities were related to metabolic rates, mainly SMR, in the intestine and heart. Individual loci were not associated with the enzymatic activities, but we found epistatic effects and genotype-by-environment interactions in CS activity in the heart and epistasis in LDH activity in the intestine. These effects suggest that mitochondrial density and anaerobic capacity in the heart and intestine may partly mediate variation in metabolic rates and life history via age-at-maturity. This article is part of the theme issue 'The evolutionary significance of variation in metabolic rates'.

Genetic interactions of schizophrenia using gene-based statistical epistasis exclusively identify nervous system-related pathways and key hub genes

Background: The relationship between genotype and phenotype is governed by numerous genetic interactions (GIs), and the mapping of GI networks is of interest for two main reasons: 1) By modelling biological robustness, GIs provide a powerful opportunity to infer compensatory biological mechanisms via the identification of functional relationships between genes, which is of interest for biological discovery and translational research. Biological systems have evolved to compensate for genetic (i.e., variations and mutations) and environmental (i.e., drug efficacy) perturbations by exploiting compensatory relationships between genes, pathways and biological processes; 2) GI facilitates the identification of the direction (alleviating or aggravating interactions) and magnitude of epistatic interactions that influence the phenotypic outcome. The generation of GIs for human diseases is impossible using experimental biology approaches such as systematic deletion analysis. Moreover, the generation of disease-specific GIs has never been undertaken in humans. Methods: We used our Indian schizophrenia case-control (case-816, controls-900) genetic dataset to implement the workflow. Standard GWAS sample quality control procedure was followed. We used the imputed genetic data to increase the SNP coverage to analyse epistatic interactions across the genome comprehensively. Using the odds ratio (OR), we identified the GIs that increase or decrease the risk of a disease phenotype. The SNP-based epistatic results were transformed into gene-based epistatic results. Results: We have developed a novel approach by conducting gene-based statistical epistatic analysis using an Indian schizophrenia case-control genetic dataset and transforming these results to infer GIs that increase the risk of schizophrenia. There were ∼9.5 million GIs with a p-value ≤ 1 × 10-5. Approximately 4.8 million GIs showed an increased risk (OR > 1.0), while ∼4.75 million GIs had a decreased risk (OR <1.0) for schizophrenia. Conclusion: Unlike model organisms, this approach is specifically viable in humans due to the availability of abundant disease-specific genome-wide genotype datasets. The study exclusively identified brain/nervous system-related processes, affirming the findings. This computational approach fills a critical gap by generating practically non-existent heritable disease-specific human GIs from human genetic data. These novel datasets can train innovative deep-learning models, potentially surpassing the limitations of conventional GWAS.

The genetics of gaits in Icelandic horses goes beyond DMRT3, with RELN and STAU2 identified as two new candidate genes

BACKGROUND

In domesticated animals, many important traits are complex and regulated by a large number of genes, genetic interactions, and environmental influences. The ability of Icelandic horses to perform the gait 'pace' is largely influenced by a single mutation in the DMRT3 gene, but genetic modifiers likely exist. The aim of this study was to identify novel genetic factors that influence pacing ability and quality of the gait through a genome-wide association study (GWAS) and correlate new findings to previously identified quantitative trait loci (QTL) and mutations.

RESULTS

Three hundred and seventy-two Icelandic horses were genotyped with the 670 K+ Axiom Equine Genotyping Array, of which 362 had gait scores from breeding field tests. A GWAS revealed several SNPs on Equus caballus chromosomes (ECA) 4, 9, and 20 that were associated (p < 1.0 × 10-5) with the breeding field test score for pace. The two novel QTL on ECA4 and 9 were located within the RELN and STAU2 genes, respectively, which have previously been associated with locomotor behavior in mice. Haplotypes were identified and the most frequent one for each of these two QTL had a large favorable effect on pace score. The second most frequent haplotype for the RELN gene was positively correlated with scores for tölt, trot, gallop, and canter. Similarly, the second most frequent haplotype for the STAU2 gene had favorable effects on scores for trot and gallop. Different genotype ratios of the haplotypes in the RELN and STAU2 genes were also observed in groups of horses with different levels of pacing ability. Furthermore, interactions (p < 0.05) were detected for the QTL in the RELN and STAU2 genes with the DMRT3 gene. The novel QTL on ECA4, 9, and 20, along with the effects of the DMRT3 variant, were estimated to account jointly for 27.4% of the phenotypic variance of the gait pace.

CONCLUSIONS

Our findings provide valuable information about the genetic architecture of pace beyond the contribution of the DMRT3 gene and indicate genetic interactions that contribute to the complexity of this trait. Further investigation is needed to fully understand the underlying genetic factors and interactions.

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.

Idiosyncratic and dose-dependent epistasis drives variation in tomato fruit size

Epistasis between genes is traditionally studied with mutations that eliminate protein activity, but most natural genetic variation is in cis-regulatory DNA and influences gene expression and function quantitatively. In this study, we used natural and engineered cis-regulatory alleles in a plant stem-cell circuit to systematically evaluate epistatic relationships controlling tomato fruit size. Combining a promoter allelic series with two other loci, we collected over 30,000 phenotypic data points from 46 genotypes to quantify how allele strength transforms epistasis. We revealed a saturating dose-dependent relationship but also allele-specific idiosyncratic interactions, including between alleles driving a step change in fruit size during domestication. Our approach and findings expose an underexplored dimension of epistasis, in which cis-regulatory allelic diversity within gene regulatory networks elicits nonlinear, unpredictable interactions that shape phenotypes.

Role of genetic architecture in phenotypic plasticity

Phenotypic plasticity, the ability of an organism to display different phenotypes across environments, is widespread in nature. Plasticity aids survival in novel environments. Herein, we review studies from yeast that allow us to start uncovering the genetic architecture of phenotypic plasticity. Genetic variants and their interactions impact the phenotype in different environments, and distinct environments modulate the impact of genetic variants and their interactions on the phenotype. Because of this, certain hidden genetic variation is expressed in specific genetic and environmental backgrounds. A better understanding of the genetic mechanisms of phenotypic plasticity will help to determine short- and long-term responses to selection and how wide variation in disease manifestation occurs in human populations.

Genomic prediction in hybrid breeding: I. Optimizing the training set design

Training sets produced by maximizing the number of parent lines, each involved in one cross, had the highest prediction accuracy for H0 hybrids, but lowest for H1 and H2 hybrids. Genomic prediction holds great promise for hybrid breeding but optimum composition of the training set (TS) as determined by the number of parents (nTS) and crosses per parent (c) has received little attention. Our objective was to examine prediction accuracy ([Formula: see text]) of GCA for lines used as parents of the TS (I1 lines) or not (I0 lines), and H0, H1 and H2 hybrids, comprising crosses of type I0 × I0, I1 × I0 and I1 × I1, respectively, as function of nTS and c. In the theory, we developed estimates for [Formula: see text] of GBLUPs for hybrids: (i)[Formula: see text] based on the expected prediction accuracy, and (ii) [Formula: see text] based on [Formula: see text] of GBLUPs of GCA and SCA effects. In the simulation part, hybrid populations were generated using molecular data from two experimental maize data sets. Additive and dominance effects of QTL borrowed from literature were used to simulate six scenarios of traits differing in the proportion (τSCA = 1%, 6%, 22%) of SCA variance in σG2 and heritability (h2 = 0.4, 0.8). Values of [Formula: see text] and [Formula: see text] closely agreed with [Formula: see text] for hybrids. For given size NTS = nTS × c of TS, [Formula: see text] of H0 hybrids and GCA of I0 lines was highest for c = 1. Conversely, for GCA of I1 lines and H1 and H2 hybrids, c = 1 yielded lowest [Formula: see text] with concordant results across all scenarios for both data sets. In view of these opposite trends, the optimum choice of c for maximizing selection response across all types of hybrids depends on the size and resources of the breeding program.

Higher-order epistasis shapes natural variation in germ stem cell niche activity

To study how natural allelic variation explains quantitative developmental system variation, we characterized natural differences in germ stem cell niche activity, measured as progenitor zone (PZ) size, between two Caenorhabditis elegans isolates. Linkage mapping yielded candidate loci on chromosomes II and V, and we found that the isolate with a smaller PZ size harbours a 148 bp promoter deletion in the Notch ligand, lag-2/Delta, a central signal promoting germ stem cell fate. As predicted, introducing this deletion into the isolate with a large PZ resulted in a smaller PZ size. Unexpectedly, restoring the deleted ancestral sequence in the isolate with a smaller PZ did not increase-but instead further reduced-PZ size. These seemingly contradictory phenotypic effects are explained by epistatic interactions between the lag-2/Delta promoter, the chromosome II locus, and additional background loci. These results provide first insights into the quantitative genetic architecture regulating an animal stem cell system.

Widespread epistasis among beneficial genetic variants revealed by high-throughput genome editing

The phenotypic effect of any genetic variant can be altered by variation at other genomic loci. Known as epistasis, these genetic interactions shape the genotype-phenotype map of every species, yet their origins remain poorly understood. To investigate this, we employed high-throughput genome editing to measure the fitness effects of 1,826 naturally polymorphic variants in four strains of Saccharomyces cerevisiae. About 31% of variants affect fitness, of which 24% have strain-specific fitness effects indicative of epistasis. We found that beneficial variants are more likely to exhibit genetic interactions and that these interactions can be mediated by specific traits such as flocculation ability. This work suggests that adaptive evolution will often involve trade-offs where a variant is only beneficial in some genetic backgrounds, potentially explaining why many beneficial variants remain polymorphic. In sum, we provide a framework to understand the factors influencing epistasis with single-nucleotide resolution, revealing widespread epistasis among beneficial variants.

Large-scale analysis of de novo mutations identifies risk genes for female infertility characterized by oocyte and early embryo defects

BACKGROUND

Oocyte maturation arrest and early embryonic arrest are important reproductive phenotypes resulting in female infertility and cause the recurrent failure of assisted reproductive technology (ART). However, the genetic etiologies of these female infertility-related phenotypes are poorly understood. Previous studies have mainly focused on inherited mutations based on large pedigrees or consanguineous patients. However, the role of de novo mutations (DNMs) in these phenotypes remains to be elucidated.

RESULTS

To decipher the role of DNMs in ART failure and female infertility with oocyte and embryo defects, we explore the landscape of DNMs in 473 infertile parent-child trios and identify a set of 481 confident DNMs distributed in 474 genes. Gene ontology analysis reveals that the identified genes with DNMs are enriched in signaling pathways associated with female reproductive processes such as meiosis, embryonic development, and reproductive structure development. We perform functional assays on the effects of DNMs in a representative gene Tubulin Alpha 4a (TUBA4A), which shows the most significant enrichment of DNMs in the infertile parent-child trios. DNMs in TUBA4A disrupt the normal assembly of the microtubule network in HeLa cells, and microinjection of DNM TUBA4A cRNAs causes abnormalities in mouse oocyte maturation or embryo development, suggesting the pathogenic role of these DNMs in TUBA4A.

CONCLUSIONS

Our findings suggest novel genetic insights that DNMs contribute to female infertility with oocyte and embryo defects. This study also provides potential genetic markers and facilitates the genetic diagnosis of recurrent ART failure and female infertility.

Intragenic compensation through the lens of deep mutational scanning

A significant fraction of mutations in proteins are deleterious and result in adverse consequences for protein function, stability, or interaction with other molecules. Intragenic compensation is a specific case of positive epistasis when a neutral missense mutation cancels effect of a deleterious mutation in the same protein. Permissive compensatory mutations facilitate protein evolution, since without them all sequences would be extremely conserved. Understanding compensatory mechanisms is an important scientific challenge at the intersection of protein biophysics and evolution. In human genetics, intragenic compensatory interactions are important since they may result in variable penetrance of pathogenic mutations or fixation of pathogenic human alleles in orthologous proteins from related species. The latter phenomenon complicates computational and clinical inference of an allele's pathogenicity. Deep mutational scanning is a relatively new technique that enables experimental studies of functional effects of thousands of mutations in proteins. We review the important aspects of the field and discuss existing limitations of current datasets. We reviewed ten published DMS datasets with quantified functional effects of single and double mutations and described rates and patterns of intragenic compensation in eight of them.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01005-w.

Loss-of-function mutation survey revealed that genes with background-dependent fitness are rare and functionally related in yeast

In different individuals, the same mutation can lead to different phenotypes due to genetic background effects. This is commonly observed in various systems, including many human diseases. While isolated examples of such background effects have been observed, a systematic view across a large number of individuals is still lacking. Here, we surveyed genetic background effects associated with gene loss-of-function mutations across a population of natural isolates of the yeast Saccharomyces cerevisiae. We found that ∼15% of genes can display a background-dependent fitness change. Genes related to mitochondrial functions are significantly overrepresented, and showed reversed patterns of fitness gain or loss with genes involved in transcription and chromatin remodeling as well as in nuclear–cytoplasmic transport, suggesting a potential functional rewiring.

In natural populations, the same mutation can lead to different phenotypic outcomes due to the genetic variation that exists among individuals. Such genetic background effects are commonly observed, including in the context of many human diseases. However, systematic characterization of these effects at the species level is still lacking to date. Here, we sought to comprehensively survey background-dependent traits associated with gene loss-of-function (LoF) mutations in 39 natural isolates of Saccharomyces cerevisiae using a transposon saturation strategy. By analyzing the modeled fitness variability of a total of 4,469 genes, we found that 15% of them, when impacted by a LoF mutation, exhibited a significant gain- or loss-of-fitness phenotype in certain natural isolates compared with the reference strain S288C. Out of these 632 genes with predicted background-dependent fitness effects, around 2/3 impact multiple backgrounds with a gradient of predicted fitness change while 1/3 are specific to a single genetic background. Genes related to mitochondrial function are significantly overrepresented in the set of genes showing a continuous variation and display a potential functional rewiring with other genes involved in transcription and chromatin remodeling as well as in nuclear–cytoplasmic transport. Such rewiring effects are likely modulated by both the genetic background and the environment. While background-specific cases are rare and span diverse cellular processes, they can be functionally related at the individual level. All genes with background-dependent fitness effects tend to have an intermediate connectivity in the global genetic interaction network and have shown relaxed selection pressure at the population level, highlighting their potential evolutionary characteristics.

A Random Forest-Based Genome-Wide Scan Reveals Fertility-Related Candidate Genes and Potential Inter-Chromosomal Epistatic Regions Associated With Age at First Calving in Nellore Cattle

This study aimed to perform a genome-wide association analysis (GWAS) using the Random Forest (RF) approach for scanning candidate genes for age at first calving (AFC) in Nellore cattle. Additionally, potential epistatic effects were investigated using linear mixed models with pairwise interactions between all markers with high importance scores within the tree ensemble non-linear structure. Data from Nellore cattle were used, including records of animals born between 1984 and 2015 and raised in commercial herds located in different regions of Brazil. The estimated breeding values (EBV) were computed and used as the response variable in the genomic analyses. After quality control, the remaining number of animals and SNPs considered were 3,174 and 360,130, respectively. Five independent RF analyses were carried out, considering different initialization seeds. The importance score of each SNP was averaged across the independent RF analyses to rank the markers according to their predictive relevance. A total of 117 SNPs associated with AFC were identified, which spanned 10 autosomes (2, 3, 5, 10, 11, 17, 18, 21, 24, and 25). In total, 23 non-overlapping genomic regions embedded 262 candidate genes for AFC. Enrichment analysis and previous evidence in the literature revealed that many candidate genes annotated close to the lead SNPs have key roles in fertility, including embryo pre-implantation and development, embryonic viability, male germinal cell maturation, and pheromone recognition. Furthermore, some genomic regions previously associated with fertility and growth traits in Nellore cattle were also detected in the present study, reinforcing the effectiveness of RF for pre-screening candidate regions associated with complex traits. Complementary analyses revealed that many SNPs top-ranked in the RF-based GWAS did not present a strong marginal linear effect but are potentially involved in epistatic hotspots between genomic regions in different autosomes, remarkably in the BTAs 3, 5, 11, and 21. The reported results are expected to enhance the understanding of genetic mechanisms involved in the biological regulation of AFC in this cattle breed.

Three legs of the missing heritability problem

The so-called 'missing heritability problem' is often characterized by behavior geneticists as a numerical discrepancy between alternative kinds of heritability. For example, while 'traditional heritability' derived from twin and family studies indicates that approximately ∼50% of variation in intelligence is attributable to genetics, 'SNP heritability' derived from genome-wide association studies indicates that only ∼10% of variation in intelligence is attributable to genetics. This 40% gap in variance accounted for by alternative kinds of heritability is frequently referred to as what's "missing." Philosophers have picked up on this reading, suggesting that "dissolving" the missing heritability problem is merely a matter of closing the numerical gap between traditional and molecular kinds of heritability. We argue that this framing of the problem undervalues the severity of the many challenges to scientific understanding of the "heritability" of human behavior. On our view, resolving the numerical discrepancies between alternative kinds of heritability will do little to advance scientific explanation and understanding of behavior genetics. Thus, we propose a new conceptual framework of the missing heritability problem that comprises three independent methodological and explanatory challenges: the numerical gap, the prediction gap, and the mechanism gap.

An Epistatic Network Describes oppA and glgB as Relevant Genes for Mycobacterium tuberculosis

Mycobacterium tuberculosis is an acid-fast bacterium that causes tuberculosis worldwide. The role of epistatic interactions among different loci of the M. tuberculosis genome under selective pressure may be crucial for understanding the disease and the molecular basis of antibiotic resistance acquisition. Here, we analyzed polymorphic loci interactions by applying a model-free method for epistasis detection, SpydrPick, on a pan–genome-wide alignment created from a set of 254 complete reference genomes. By means of the analysis of an epistatic network created with the detected epistatic interactions, we found that glgB (α-1,4-glucan branching enzyme) and oppA (oligopeptide-binding protein) are putative targets of co-selection in M. tuberculosis as they were associated in the network with M. tuberculosis genes related to virulence, pathogenesis, transport system modulators of the immune response, and antibiotic resistance. In addition, our work unveiled potential pharmacological applications for genotypic antibiotic resistance inherent to the mutations of glgB and oppA as they epistatically interact with fprA and embC, two genes recently included as antibiotic-resistant genes in the catalog of the World Health Organization. Our findings showed that this approach allows the identification of relevant epistatic interactions that may lead to a better understanding of M. tuberculosis by deciphering the complex interactions of molecules involved in its metabolism, virulence, and pathogenesis and that may be applied to different bacterial populations.

Complex fitness landscape shapes variation in a hyperpolymorphic species

It is natural to assume that patterns of genetic variation in hyperpolymorphic species can reveal large-scale properties of the fitness landscape that are hard to detect by studying species with ordinary levels of genetic variation. Here, we study such patterns in a fungus Schizophyllum commune, the most polymorphic species known. Throughout the genome, short-range linkage disequilibrium (LD) caused by attraction of minor alleles is higher between pairs of nonsynonymous than of synonymous variants. This effect is especially pronounced for pairs of sites that are located within the same gene, especially if a large fraction of the gene is covered by haploblocks, genome segments where the gene pool consists of two highly divergent haplotypes, which is a signature of balancing selection. Haploblocks are usually shorter than 1000 nucleotides, and collectively cover about 10% of the S. commune genome. LD tends to be substantially higher for pairs of nonsynonymous variants encoding amino acids that interact within the protein. There is a substantial correlation between LDs at the same pairs of nonsynonymous mutations in the USA and the Russian populations. These patterns indicate that selection in S. commune involves positive epistasis due to compensatory interactions between nonsynonymous alleles. When less polymorphic species are studied, analogous patterns can be detected only through interspecific comparisons.

Changes to DNA known as mutations may alter how the proteins and other components of a cell work, and thus play an important role in allowing living things to evolve new traits and abilities over many generations. Whether a mutation is beneficial or harmful may differ depending on the genetic background of the individual – that is, depending on other mutations present in other positions within the same gene – due to a phenomenon called epistasis.

Epistasis is known to affect how various species accumulate differences in their DNA compared to each other over time. For example, a mutation that is rare in humans and known to cause disease may be widespread in other primates because its negative effect is canceled out by another mutation that is standard for these species but absent in humans. However, it remains unclear whether epistasis plays a significant part in shaping genetic differences between individuals of the same species.

A type of fungus known as Schizophyllum commune lives on rotting wood and is found across the world. It is one of the most genetically diverse species currently known, so there is a higher chance of pairs of compensatory mutations occurring and persisting for a long time in S. commune than in most other species, providing a unique opportunity to study epistasis.

Here, Stolyarova et al. studied two distinct populations of S. commune, one from the USA and one from Russia. The team found that – unlike in humans, flies and other less genetically diverse species – epistasis maintains combinations of mutations in S. commune that individually would be harmful to the fungus but together compensate for each other. For example, pairs of mutations affecting specific molecules known as amino acids – the building blocks of proteins – that physically interact with each other tended to be found together in the same individuals.

One potential downside of having pairs of compensatory mutations in the genome is that when the organism reproduces, the process of making sex cells may split up these pairs so that harmful mutations are inherited without their partner mutations. Thus, epistasis may have helped shape the way S. commune and other genetically diverse species have evolved.

Prediction performance of linear models and gradient boosting machine on complex phenotypes in outbred mice

We compared the performance of linear (GBLUP, BayesB, and elastic net) methods to a nonparametric tree-based ensemble (gradient boosting machine) method for genomic prediction of complex traits in mice. The dataset used contained genotypes for 50,112 SNP markers and phenotypes for 835 animals from 6 generations. Traits analyzed were bone mineral density, body weight at 10, 15, and 20 weeks, fat percentage, circulating cholesterol, glucose, insulin, triglycerides, and urine creatinine. The youngest generation was used as a validation subset, and predictions were based on all older generations. Model performance was evaluated by comparing predictions for animals in the validation subset against their adjusted phenotypes. Linear models outperformed gradient boosting machine for 7 out of 10 traits. For bone mineral density, cholesterol, and glucose, the gradient boosting machine model showed better prediction accuracy and lower relative root mean squared error than the linear models. Interestingly, for these 3 traits, there is evidence of a relevant portion of phenotypic variance being explained by epistatic effects. Using a subset of top markers selected from a gradient boosting machine model helped for some of the traits to improve the accuracy of prediction when these were fitted into linear and gradient boosting machine models. Our results indicate that gradient boosting machine is more strongly affected by data size and decreased connectedness between reference and validation sets than the linear models. Although the linear models outperformed gradient boosting machine for the polygenic traits, our results suggest that gradient boosting machine is a competitive method to predict complex traits with assumed epistatic effects.

Global urban environmental change drives adaptation in white clover

Urbanization transforms environments in ways that alter biological evolution. We examined whether urban environmental change drives parallel evolution by sampling 110,019 white clover plants from 6169 populations in 160 cities globally. Plants were assayed for a Mendelian antiherbivore defense that also affects tolerance to abiotic stressors. Urban-rural gradients were associated with the evolution of clines in defense in 47% of cities throughout the world. Variation in the strength of clines was explained by environmental changes in drought stress and vegetation cover that varied among cities. Sequencing 2074 genomes from 26 cities revealed that the evolution of urban-rural clines was best explained by adaptive evolution, but the degree of parallel adaptation varied among cities. Our results demonstrate that urbanization leads to adaptation at a global scale.

Systemic lupus erythematosus as a genetic disease

Systemic lupus erythematosus is the prototypical systemic autoimmune disease, as it is characterized both by protean multi-organ system manifestations and by the uniform presence of pathogenic autoantibodies directed against components of the nucleus. Prior to the modern genetic era, the diverse clinical manifestations of SLE suggested to many that SLE patients were unlikely to share a common genetic risk basis. However, modern genetic studies have revealed that SLE usually arises when an environmental exposure occurs in an individual with a collection of genetic risk variants passing a liability threshold. Here, we summarize the current state of the field aimed at: (1) understanding the genetic architecture of this complex disease, (2) synthesizing how this genetic risk architecture impacts cellular and molecular disease pathophysiology, (3) providing illustrative examples that highlight the rich complexity of the pathobiology of this prototypical autoimmune disease and (4) communicating this complex etiopathogenesis to patients.

Epistatic genetic interactions govern morphogenesis during sexual reproduction and infection in a global human fungal pathogen

Cellular development is orchestrated by evolutionarily conserved signaling pathways, which are often pleiotropic and involve intra- and interpathway epistatic interactions that form intricate, complex regulatory networks. Cryptococcus species are a group of closely related human fungal pathogens that grow as yeasts yet transition to hyphae during sexual reproduction. Additionally, during infection they can form large, polyploid titan cells that evade immunity and develop drug resistance. Multiple known signaling pathways regulate cellular development, yet how these are coordinated and interact with genetic variation is less well understood. Here, we conducted quantitative trait locus (QTL) analyses of a mapping population generated by sexual reproduction of two parents, only one of which is unisexually fertile. We observed transgressive segregation of the unisexual phenotype among progeny, as well as a large-cell phenotype under mating-inducing conditions. These large-cell progeny were found to produce titan cells both in vitro and in infected animals. Two major QTLs and corresponding quantitative trait genes (QTGs) were identified: RIC8 (encoding a guanine-exchange factor) and CNC06490 (encoding a putative Rho-GTPase activator), both involved in G protein signaling. The two QTGs interact epistatically with each other and with the mating-type locus in phenotypic determination. These findings provide insights into the complex genetics of morphogenesis during unisexual reproduction and pathogenic titan cell formation and illustrate how QTL analysis can be applied to identify epistasis between genes. This study shows that phenotypic outcomes are influenced by the genetic background upon which mutations arise, implicating dynamic, complex genotype-to-phenotype landscapes in fungal pathogens and beyond.

Mitochondria as environments for the nuclear genome in Drosophila: mitonuclear G×G×E

Mitochondria evolved from a union of microbial cells belonging to distinct lineages that were likely anaerobic. The evolution of eukaryotes required a massive reorganization of the 2 genomes and eventual adaptation to aerobic environments. The nutrients and oxygen that sustain eukaryotic metabolism today are processed in mitochondria through coordinated expression of 37 mitochondrial genes and over 1000 nuclear genes. This puts mitochondria at the nexus of gene-by-gene (G×G) and gene-by-environment (G×E) interactions that sustain life. Here we use a Drosophila model of mitonuclear genetic interactions to explore the notion that mitochondria are environments for the nuclear genome, and vice versa. We construct factorial combinations of mtDNA and nuclear chromosomes to test for epistatic interactions (G×G), and expose these mitonuclear genotypes to altered dietary environments to examine G×E interactions. We use development time and genome-wide RNAseq analyses to assess the relative contributions of mtDNA, nuclear chromosomes, and environmental effects on these traits (mitonuclear G×G×E). We show that the nuclear transcriptional response to alternative mitochondrial "environments" (G×G) has significant overlap with the transcriptional response of mitonuclear genotypes to altered dietary environments. These analyses point to specific transcription factors (e.g., giant) that mediated these interactions, and identified coexpressed modules of genes that may account for the overlap in differentially expressed genes. Roughly 20% of the transcriptome includes G×G genes that are concordant with G×E genes, suggesting that mitonuclear interactions are part of an organism's environment.

The mutation effect reaction norm (mu-rn) highlights environmentally dependent mutation effects and epistatic interactions

Since the modern synthesis, the fitness effects of mutations and epistasis have been central yet provocative concepts in evolutionary and population genetics. Studies of how the interactions between parcels of genetic information can change as a function of environmental context have added a layer of complexity to these discussions. Here I introduce the "mutation effect reaction norm" (Mu-RN), a new instrument through which one can analyze the phenotypic consequences of mutations and interactions across environmental contexts. It embodies the fusion of measurements of genetic interactions with the reaction norm, a classic depiction of the performance of genotypes across environments. I demonstrate the utility of the Mu-RN through the signature of a "compensatory ratchet" mutation that undermines reverse evolution of antimicrobial resistance. More broadly, I argue that the mutation effect reaction norm may help us resolve the dynamism and unpredictability of evolution, with implications for theoretical biology, genetic modification technology, and public health. This article is protected by copyright. All rights reserved.

Gene-Environment Interactions for Cardiovascular Disease

PURPOSE OF REVIEW

We provide an overview of recent findings with respect to gene-environment (GxE) interactions for cardiovascular disease (CVD) risk and discuss future opportunities for advancing the field.

RECENT FINDINGS

Over the last several years, GxE interactions for CVD have mostly been identified for smoking and coronary artery disease (CAD) or related risk factors. By comparison, there is more limited evidence for GxE interactions between CVD outcomes and other exposures, such as physical activity, air pollution, diet, and sex. The establishment of large consortia and population-based cohorts, in combination with new computational tools and mouse genetics platforms, can potentially overcome some of the limitations that have hindered human GxE interaction studies and reveal additional association signals for CVD-related traits. The identification of novel GxE interactions is likely to provide a better understanding of the pathogenesis and genetic liability of CVD, with significant implications for healthy lifestyles and therapeutic strategies.

Decoding 'Unnecessary Complexity': A Law of Complexity and a Concept of Hidden Variation Behind "Missing Heritability" in Precision Medicine

The high hopes for the Human Genome Project and personalized medicine were not met because the relationship between genotypes and phenotypes turned out to be more complex than expected. In a previous study we laid the foundation of a theory of complexity and showed that because of the blind nature of evolution, and molecular and historical contingency, cells have accumulated unnecessary complexity, complexity beyond what is necessary and sufficient to describe an organism. Here we provide empirical evidence and show that unnecessary complexity has become integrated into the genome in the form of redundancy and is relevant to molecular evolution of phenotypic complexity. Unnecessary complexity creates uncertainty between molecular and phenotypic complexity, such that phenotypic complexity (CP) is higher than molecular complexity (CM), which is higher than DNA complexity (CD). The qualitative inequality in complexity is based on the following hierarchy: CP > CM > CD. This law-like relationship holds true for all complex traits, including complex diseases. We present a hypothesis of two types of variation, namely open and closed (hidden) systems, show that hidden variation provides a hitherto undiscovered "third source" of phenotypic variation, beside genotype and environment, and argue that "missing heritability" for some complex diseases is likely to be a case of "diluted heritability". There is a need for radically new ways of thinking about the principles of genotype-phenotype relationship. Understanding how cells use hidden, pathway variation to respond to stress can shed light on why two individuals who share the same risk factors may not develop the same disease, or how cancer cells escape death.

Inferring multilayer interactome networks shaping phenotypic plasticity and evolution

Phenotypic plasticity represents a capacity by which the organism changes its phenotypes in response to environmental stimuli. Despite its pivotal role in adaptive evolution, how phenotypic plasticity is genetically controlled remains elusive. Here, we develop a unified framework for coalescing all single nucleotide polymorphisms (SNPs) from a genome-wide association study (GWAS) into a quantitative graph. This framework integrates functional genetic mapping, evolutionary game theory, and predator-prey theory to decompose the net genetic effect of each SNP into its independent and dependent components. The independent effect arises from the intrinsic capacity of a SNP, only expressed when it is in isolation, whereas the dependent effect results from the extrinsic influence of other SNPs. The dependent effect is conceptually beyond the traditional definition of epistasis by not only characterizing the strength of epistasis but also capturing the bi-causality of epistasis and the sign of the causality. We implement functional clustering and variable selection to infer multilayer, sparse, and multiplex interactome networks from any dimension of genetic data. We design and conduct two GWAS experiments using Staphylococcus aureus, aimed to test the genetic mechanisms underlying the phenotypic plasticity of this species to vancomycin exposure and Escherichia coli coexistence. We reconstruct the two most comprehensive genetic networks for abiotic and biotic phenotypic plasticity. Pathway analysis shows that SNP-SNP epistasis for phenotypic plasticity can be annotated to protein-protein interactions through coding genes. Our model can unveil the regulatory mechanisms of significant loci and excavate missing heritability from some insignificant loci. Our multilayer genetic networks provide a systems tool for dissecting environment-induced evolution.

Human genetic analyses of organelles highlight the nucleus in age-related trait heritability

Most age-related human diseases are accompanied by a decline in cellular organelle integrity, including impaired lysosomal proteostasis and defective mitochondrial oxidative phosphorylation. An open question, however, is the degree to which inherited variation in or near genes encoding each organelle contributes to age-related disease pathogenesis. Here, we evaluate if genetic loci encoding organelle proteomes confer greater-than-expected age-related disease risk. As mitochondrial dysfunction is a 'hallmark' of aging, we begin by assessing nuclear and mitochondrial DNA loci near genes encoding the mitochondrial proteome and surprisingly observe a lack of enrichment across 24 age-related traits. Within nine other organelles, we find no enrichment with one exception: the nucleus, where enrichment emanates from nuclear transcription factors. In agreement, we find that genes encoding several organelles tend to be 'haplosufficient,' while we observe strong purifying selection against heterozygous protein-truncating variants impacting the nucleus. Our work identifies common variation near transcription factors as having outsize influence on age-related trait risk, motivating future efforts to determine if and how this inherited variation then contributes to observed age-related organelle deterioration.

The Genetic Architecture of a Congenital Heart Defect Is Related to Its Fitness Cost

In newborns, severe congenital heart defects are rarer than mild ones. This epidemiological relationship between heart defect severity and incidence lacks explanation. Here, an analysis of ~10,000 Nkx2-5^+/− mice from two inbred strain crosses illustrates the fundamental role of epistasis. Modifier genes raise or lower the risk of specific defects via pairwise (G×G_Nkx) and higher-order (G×G×G_Nkx) interactions with Nkx2-5. Their effect sizes correlate with the severity of a defect. The risk loci for mild, atrial septal defects exert predominantly small G×G_Nkx effects, while the loci for severe, atrioventricular septal defects exert large G×G_Nkx and G×G×G_Nkx effects. The loci for moderately severe ventricular septal defects have intermediate effects. Interestingly, G×G×G_Nkx effects are three times more likely to suppress risk when the genotypes at the first two loci are from the same rather than different parental inbred strains. This suggests the genetic coadaptation of interacting G×G×G_Nkx loci, a phenomenon that Dobzhansky first described in Drosophila. Thus, epistasis plays dual roles in the pathogenesis of congenital heart disease and the robustness of cardiac development. The empirical results suggest a relationship between the fitness cost and genetic architecture of a disease phenotype and a means for phenotypic robustness to have evolved.

Relevance of Higher-Order Epistasis in Drug Resistance

We studied five chemically distinct but related 1,3,5-triazine antifolates with regard to their effects on growth of a set of mutants in dihydrofolate reductase. The mutants comprise a combinatorially complete data set of all 16 possible combinations of four amino acid replacements associated with resistance to pyrimethamine in the malaria parasite Plasmodium falciparum. Pyrimethamine was a mainstay medication for malaria for many years, and it is still in use in intermittent treatment during pregnancy or as a partner drug in artemisinin combination therapy. Our goal was to investigate the extent to which the alleles yield similar adaptive topographies and patterns of epistasis across chemically related drugs. We find that the adaptive topographies are indeed similar with the same or closely related alleles being fixed in computer simulations of stepwise evolution. For all but one of the drugs the topography features at least one suboptimal fitness peak. Our data are consistent with earlier results indicating that third order and higher epistatic interactions appear to contribute only modestly to the overall adaptive topography, and they are largely conserved. In regard to drug development, our data suggest that higher-order interactions are likely to be of little value as an advisory tool in the choice of lead compounds.

Modeling genome-wide by environment interactions through omnigenic interactome networks

How genes interact with the environment to shape phenotypic variation and evolution is a fundamental question intriguing to biologists from various fields. Existing linear models built on single genes are inadequate to reveal the complexity of genotype-environment (G-E) interactions. Here, we develop a conceptual model for mechanistically dissecting G-E interplay by integrating previously disconnected theories and methods. Under this integration, evolutionary game theory, developmental modularity theory, and a variable selection method allow us to reconstruct environment-induced, maximally informative, sparse, and casual multilayer genetic networks. We design and conduct two mapping experiments by using a desert-adapted tree species to validate the biological application of the model proposed. The model identifies previously uncharacterized molecular mechanisms that mediate trees' response to saline stress. Our model provides a tool to comprehend the genetic architecture of trait variation and evolution and trace the information flow of each gene toward phenotypes within omnigenic networks.

Diversifying maize genomic selection models

Genomic selection (GS) is one of the most powerful tools available for maize breeding. Its use of genome-wide marker data to estimate breeding values translates to increased genetic gains with fewer breeding cycles. In this review, we cover the history of GS and highlight particular milestones during its adaptation to maize breeding. We discuss how GS can be applied to developing superior maize inbreds and hybrids. Additionally, we characterize refinements in GS models that could enable the encapsulation of non-additive genetic effects, genotype by environment interactions, and multiple levels of the biological hierarchy, all of which could ultimately result in more accurate predictions of breeding values. Finally, we suggest the stages in a maize breeding program where it would be beneficial to apply GS. Given the current sophistication of high-throughput phenotypic, genotypic, and other -omic level data currently available to the maize community, now is the time to explore the implications of their incorporation into GS models and thus ensure that genetic gains are being achieved as quickly and efficiently as possible.

The Origin of Additive Genetic Variance Driven by Positive Selection

Fisher's fundamental theorem of natural selection predicts no additive variance of fitness in a natural population. Consistently, studies in a variety of wild populations show virtually no narrow-sense heritability (h2) for traits important to fitness. However, counterexamples are occasionally reported, calling for a deeper understanding on the evolution of additive variance. In this study, we propose adaptive divergence followed by population admixture as a source of the additive genetic variance of evolutionarily important traits. We experimentally tested the hypothesis by examining a panel of ∼1,000 yeast segregants produced by a hybrid of two yeast strains that experienced adaptive divergence. We measured >400 yeast cell morphological traits and found a strong positive correlation between h2 and evolutionary importance. Because adaptive divergence followed by population admixture could happen constantly, particularly in species with wide geographic distribution and strong migratory capacity (e.g., humans), the finding reconciles the observation of abundant additive variances in evolutionarily important traits with Fisher's fundamental theorem of natural selection. Importantly, the revealed role of positive selection in promoting rather than depleting additive variance suggests a simple explanation for why additive genetic variance can be dominant in a population despite the ubiquitous between-gene epistasis observed in functional assays.

Robust modeling of additive and nonadditive variation with intuitive inclusion of expert knowledge

We propose a novel Bayesian approach that robustifies genomic modeling by leveraging expert knowledge (EK) through prior distributions. The central component is the hierarchical decomposition of phenotypic variation into additive and nonadditive genetic variation, which leads to an intuitive model parameterization that can be visualized as a tree. The edges of the tree represent ratios of variances, for example broad-sense heritability, which are quantities for which EK is natural to exist. Penalized complexity priors are defined for all edges of the tree in a bottom-up procedure that respects the model structure and incorporates EK through all levels. We investigate models with different sources of variation and compare the performance of different priors implementing varying amounts of EK in the context of plant breeding. A simulation study shows that the proposed priors implementing EK improve the robustness of genomic modeling and the selection of the genetically best individuals in a breeding program. We observe this improvement in both variety selection on genetic values and parent selection on additive values; the variety selection benefited the most. In a real case study, EK increases phenotype prediction accuracy for cases in which the standard maximum likelihood approach did not find optimal estimates for the variance components. Finally, we discuss the importance of EK priors for genomic modeling and breeding, and point to future research areas of easy-to-use and parsimonious priors in genomic modeling.

Robust modeling of additive and nonadditive variation with intuitive inclusion of expert knowledge

We propose a novel Bayesian approach that robustifies genomic modeling by leveraging expert knowledge (EK) through prior distributions. The central component is the hierarchical decomposition of phenotypic variation into additive and nonadditive genetic variation, which leads to an intuitive model parameterization that can be visualized as a tree. The edges of the tree represent ratios of variances, for example broad-sense heritability, which are quantities for which EK is natural to exist. Penalized complexity priors are defined for all edges of the tree in a bottom-up procedure that respects the model structure and incorporates EK through all levels. We investigate models with different sources of variation and compare the performance of different priors implementing varying amounts of EK in the context of plant breeding. A simulation study shows that the proposed priors implementing EK improve the robustness of genomic modeling and the selection of the genetically best individuals in a breeding program. We observe this improvement in both variety selection on genetic values and parent selection on additive values; the variety selection benefited the most. In a real case study, EK increases phenotype prediction accuracy for cases in which the standard maximum likelihood approach did not find optimal estimates for the variance components. Finally, we discuss the importance of EK priors for genomic modeling and breeding, and point to future research areas of easy-to-use and parsimonious priors in genomic modeling.

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.

Investigating Mitonuclear Genetic Interactions Through Machine Learning: A Case Study on Cold Adaptation Genes in Human Populations From Different European Climate Regions

Cold climates represent one of the major environmental challenges that anatomically modern humans faced during their dispersal out of Africa. The related adaptive traits have been achieved by modulation of thermogenesis and thermoregulation processes where nuclear (nuc) and mitochondrial (mt) genes play a major role. In human populations, mitonuclear genetic interactions are the result of both the peculiar genetic history of each human group and the different environments they have long occupied. This study aims to investigate mitonuclear genetic interactions by considering all the mitochondrial genes and 28 nuclear genes involved in brown adipose tissue metabolism, which have been previously hypothesized to be crucial for cold adaptation. For this purpose, we focused on three human populations (i.e., Finnish, British, and Central Italian people) of European ancestry from different biogeographical and climatic areas, and we used a machine learning approach to identify relevant nucDNA–mtDNA interactions that characterized each population. The obtained results are twofold: (i) at the methodological level, we demonstrated that a machine learning approach is able to detect patterns of genetic structure among human groups from different latitudes both at single genes and by considering combinations of mtDNA and nucDNA loci; (ii) at the biological level, the analysis identified population-specific nuclear genes and variants that likely play a relevant biological role in association with a mitochondrial gene (such as the “obesity gene” FTO in Finnish people). Further studies are needed to fully elucidate the evolutionary dynamics (e.g., migration, admixture, and/or local adaptation) that shaped these nucDNA–mtDNA interactions and their functional role.

Short-term heritable variation overwhelms 200 generations of mutational variance for metabolic traits in Caenorhabditis elegans

Metabolic disorders have a large heritable component, and have increased markedly in human populations over the past few generations. Genome-wide association studies of metabolic traits typically find a substantial unexplained fraction of total heritability, suggesting an important role of spontaneous mutation. An alternative explanation is that epigenetic effects contribute significantly to the heritable variation. Here, we report a study designed to quantify the cumulative effects of spontaneous mutation on adenosine metabolism in the nematode Caenorhabditis elegans, including both the activity and concentration of two metabolic enzymes and the standing pools of their associated metabolites. The only prior studies on the effects of mutation on metabolic enzyme activity, in Drosophila melanogaster, found that total enzyme activity presents a mutational target similar to that of morphological and life-history traits. However, those studies were not designed to account for short-term heritable effects. We find that the short-term heritable variance for most traits is of similar magnitude as the variance among MA lines. This result suggests that the potential heritable effects of epigenetic variation in metabolic disease warrant additional scrutiny.

Major Impacts of Widespread Structural Variation on Gene Expression and Crop Improvement in Tomato

Structural variants (SVs) underlie important crop improvement and domestication traits. However, resolving the extent, diversity, and quantitative impact of SVs has been challenging. We used long-read nanopore sequencing to capture 238,490 SVs in 100 diverse tomato lines. This panSV genome, along with 14 new reference assemblies, revealed large-scale intermixing of diverse genotypes, as well as thousands of SVs intersecting genes and cis-regulatory regions. Hundreds of SV-gene pairs exhibit subtle and significant expression changes, which could broadly influence quantitative trait variation. By combining quantitative genetics with genome editing, we show how multiple SVs that changed gene dosage and expression levels modified fruit flavor, size, and production. In the last example, higher order epistasis among four SVs affecting three related transcription factors allowed introduction of an important harvesting trait in modern tomato. Our findings highlight the underexplored role of SVs in genotype-to-phenotype relationships and their widespread importance and utility in crop improvement.

Homeostatic plasticity fails at the intersection of autism-gene mutations and a novel class of common genetic modifiers

We identify a set of common phenotypic modifiers that interact with five independent autism gene orthologs (RIMS1, CHD8, CHD2, WDFY3, ASH1L) causing a common failure of presynaptic homeostatic plasticity (PHP) in Drosophila. Heterozygous null mutations in each autism gene are demonstrated to have normal baseline neurotransmission and PHP. However, PHP is sensitized and rendered prone to failure. A subsequent electrophysiology-based genetic screen identifies the first known heterozygous mutations that commonly genetically interact with multiple ASD gene orthologs, causing PHP to fail. Two phenotypic modifiers identified in the screen, PDPK1 and PPP2R5D, are characterized. Finally, transcriptomic, ultrastructural and electrophysiological analyses define one mechanism by which PHP fails; an unexpected, maladaptive up-regulation of CREG, a conserved, neuronally expressed, stress response gene and a novel repressor of PHP. Thus, we define a novel genetic landscape by which diverse, unrelated autism risk genes may converge to commonly affect the robustness of synaptic transmission.

Genetic control of non-genetic inheritance in mammals: state-of-the-art and perspectives

Thought to be directly and uniquely dependent from genotypes, the ontogeny of individual phenotypes is much more complicated. Individual genetics, environmental exposures, and their interaction are the three main determinants of individual's phenotype. This picture has been further complicated a decade ago when the Lamarckian theory of acquired inheritance has been rekindled with the discovery of epigenetic inheritance, according to which acquired phenotypes can be transmitted through fertilization and affect phenotypes across generations. The results of Genome-Wide Association Studies have also highlighted a big degree of missing heritability in genetics and have provided hints that not only acquired phenotypes, but also individual's genotypes affect phenotypes intergenerationally through indirect genetic effects. Here, we review available examples of indirect genetic effects in mammals, what is known of the underlying molecular mechanisms and their potential impact for our understanding of missing heritability, phenotypic variation. and individual disease risk.

Systems genetics applications in metabolism research

The common forms of metabolic diseases are highly complex, involving hundreds of genes, environmental and lifestyle factors, age-related changes, sex differences and gut-microbiome interactions. Systems genetics is a population-based approach to address this complexity. In contrast to commonly used 'reductionist' approaches, such as gain or loss of function, that examine one element at a time, systems genetics uses high-throughput 'omics' technologies to quantitatively assess the many molecular differences among individuals in a population and then to relate these to physiologic functions or disease states. Unlike genome-wide association studies, systems genetics seeks to go beyond the identification of disease-causing genes to understand higher-order interactions at the molecular level. The purpose of this review is to introduce the systems genetics applications in the areas of metabolic and cardiovascular disease. Here, we explain how large clinical and omics-level data and databases from both human and animal populations are available to help researchers place genes in the context of pathways and networks and formulate hypotheses that can then be experimentally examined. We provide lists of such databases and examples of the integration of reductionist and systems genetics data. Among the important applications emerging is the development of improved nutritional and pharmacological strategies to address the rise of metabolic diseases.

Transcriptomic and presence/absence variation in the barley genome assessed from multi-tissue mRNA sequencing and their power to predict phenotypic traits

Background

Barley is the world’s fourth most cultivated cereal and is an important crop model for genetic studies. One layer of genomic information that remains poorly explored in barley is presence/absence variation (PAV), which has been suggested to contribute to phenotypic variation of agronomic importance in various crops.

Results

An mRNA sequencing approach was used to study genomic PAV and transcriptomic variation in 23 spring barley inbreds. 1502 new genes identified here were physically absent from the Morex reference sequence, and 11,523 previously unannotated genes were not expressed in Morex. The procedure applied to detect expression PAV revealed that more than 50% of all genes of our data set are not expressed in all inbreds. Interestingly, expression PAV were not in strong linkage disequilibrium with neighboring sequence variants (SV), and therefore provided an additional layer of genetic information. Optimal combinations of expression PAV, SV, and gene abundance data could enhance the prediction accuracy of predicting three different agronomic traits.

Conclusions

Our results highlight the advantage of mRNA sequencing for genomic prediction over other technologies, as it allows extracting multiple layers of genomic data from a single sequencing experiment. Finally, we propose low coverage mRNA sequencing based characterization of breeding material harvested as seedlings in petri dishes as a powerful and cost efficient approach to replace current single nucleotide polymorphism (SNP) based characterizations.

Using tree-based methods for detection of gene-gene interactions in the presence of a polygenic signal: simulation study with application to educational attainment in the Generation Scotland Cohort Study

Motivation

The genomic architecture of human complex diseases is thought to be attributable to single markers, polygenic components and epistatic components. No study has examined the ability of tree-based methods to detect epistasis in the presence of a polygenic signal. We sought to apply decision tree-based methods, C5.0 and logic regression, to detect epistasis under several simulated conditions, varying strength of interaction and linkage disequilibrium (LD) structure. We then applied the same methods to the phenotype of educational attainment in a large population cohort.

Results

LD pruning improved the power and reduced the type I error. C5.0 had a conservative type I error rate whereas logic regression had a type I error rate that exceeded 5%. Despite the more conservative type I error, C5.0 was observed to have higher power than logic regression across several conditions. In the presence of a polygenic signal, power was generally reduced. Applying both methods on educational attainment in a large population cohort yielded numerous interacting SNPs; notably a SNP in RCAN3 which is associated with reading and spelling and a SNP in NPAS3, a neurodevelopmental gene.

Availability and implementation

All methods used are implemented and freely available in R.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genetics of Human Longevity Within an Eco-Evolutionary Nature-Nurture Framework

Human longevity is a complex trait, and to disentangle its basis has a great theoretical and practical consequences for biomedicine. The genetics of human longevity is still poorly understood despite several investigations that used different strategies and protocols. Here, we argue that such rather disappointing harvest is largely because of the extraordinary complexity of the longevity phenotype in humans. The capability to reach the extreme decades of human lifespan seems to be the result of an intriguing mixture of gene-environment interactions. Accordingly, the genetics of human longevity is here described as a highly context-dependent phenomenon, within a new integrated, ecological, and evolutionary perspective, and is presented as a dynamic process, both historically and individually. The available literature has been scrutinized within this perspective, paying particular attention to factors (sex, individual biography, family, population ancestry, social structure, economic status, and education, among others) that have been relatively neglected. The strength and limitations of the most powerful and used tools, such as genome-wide association study and whole-genome sequencing, have been discussed, focusing on prominently emerged genes and regions, such as apolipoprotein E, Forkhead box O3, interleukin 6, insulin-like growth factor-1, chromosome 9p21, 5q33.3, and somatic mutations among others. The major results of this approach suggest that (1) the genetics of longevity is highly population specific; (2) small-effect alleles, pleiotropy, and the complex allele timing likely play a major role; (3) genetic risk factors are age specific and need to be integrated in the light of the geroscience perspective; (4) a close relationship between genetics of longevity and genetics of age-related diseases (especially cardiovascular diseases) do exist. Finally, the urgent need of a global approach to the largely unexplored interactions between the 3 genetics of human body, that is, nuclear, mitochondrial, and microbiomes, is stressed. We surmise that the comprehensive approach here presented will help in increasing the above-mentioned harvest.

Discovering genetic interactions bridging pathways in genome-wide association studies

Genetic interactions have been reported to underlie phenotypes in a variety of systems, but the extent to which they contribute to complex disease in humans remains unclear. In principle, genome-wide association studies (GWAS) provide a platform for detecting genetic interactions, but existing methods for identifying them from GWAS data tend to focus on testing individual locus pairs, which undermines statistical power. Importantly, a global genetic network mapped for a model eukaryotic organism revealed that genetic interactions often connect genes between compensatory functional modules in a highly coherent manner. Taking advantage of this expected structure, we developed a computational approach called BridGE that identifies pathways connected by genetic interactions from GWAS data. Applying BridGE broadly, we discover significant interactions in Parkinson's disease, schizophrenia, hypertension, prostate cancer, breast cancer, and type 2 diabetes. Our novel approach provides a general framework for mapping complex genetic networks underlying human disease from genome-wide genotype data.

The modulation of HBsAg level by sI126T is affected by additional amino acid substitutions in the S region of HBV

The hepatitis B surface antigen (HBsAg) is a vital serum marker for hepatitis B virus (HBV) infection. Amino acid (AA) substitutions in small hepatitis B surface protein (SHBs) are known to affect HBsAg level. However, how the genetic backbones of SHBs sequences would affect the roles of a specific AA substitution on HBsAg level remains unclear. In this study, we found that sI126 had a very high substitution detection rate of 17.54% (40/228) in untreated chronic hepatitis B cohort with subgenotype C2 HBV infection. Among different substitution types at sI126, the sI126T (N = 28) was found to be associated with significantly lower serum HBsAg level. Clone sequencing revealed that sI126T-harboring SHBs sequences had varied genetic backbones with zero to nine additional AA substitutions. Thus, we constructed 24 HBsAg expression plasmids harboring sI126T without (plasmid 1, P1) or with (P2-P24) additional AA substitution(s) and studied them in the HepG2 cells. The HBsAg levels were determined by both ELISA and Western blot. In vitro experiments showed that P1 significantly reduced HBsAg level and its secretion (p < .05), however, P2-P24 showed various extracellular and intracellular HBsAg levels. No significant differences were detected among the HBsAg mRNA levels of nine representative mutant plasmids. Our findings suggest that the modulation of HBsAg level by sI126T is affected by additional AA substitution(s) in the S region of HBV. The effects of AA combination substitutions in SHBs sequences on HBsAg levels are worthwhile for more attentions in terms of HBV biology and its clinical application.

Lexical Landscapes as large in silico data for examining advanced properties of fitness landscapes

In silico approaches have served a central role in the development of evolutionary theory for generations. This especially applies to the concept of the fitness landscape, one of the most important abstractions in evolutionary genetics, and one which has benefited from the presence of large empirical data sets only in the last decade or so. In this study, we propose a method that allows us to generate enormous data sets that walk the line between in silico and empirical: word usage frequencies as catalogued by the Google ngram corpora. These data can be codified or analogized in terms of a multidimensional empirical fitness landscape towards the examination of advanced concepts-adaptive landscape by environment interactions, clonal competition, higher-order epistasis and countless others. We argue that the greater Lexical Landscapes approach can serve as a platform that offers an astronomical number of fitness landscapes for exploration (at least) or theoretical formalism (potentially) in evolutionary biology.

Network approaches to systems biology analysis of complex disease: integrative methods for multi-omics data

In the past decade, significant progress has been made in complex disease research across multiple omics layers from genome, transcriptome and proteome to metabolome. There is an increasing awareness of the importance of biological interconnections, and much success has been achieved using systems biology approaches. However, because of the typical focus on one single omics layer at a time, existing systems biology findings explain only a modest portion of complex disease. Recent advances in multi-omics data collection and sharing present us new opportunities for studying complex diseases in a more comprehensive fashion, and yet simultaneously create new challenges considering the unprecedented data dimensionality and diversity. Here, our goal is to review extant and emerging network approaches that can be applied across multiple biological layers to facilitate a more comprehensive and integrative multilayered omics analysis of complex diseases.