Environment-dependent epistasis increases phenotypic diversity in gene regulatory networks

Using supervised machine-learning approaches to understand abiotic stress tolerance and design resilient crops

Abiotic stresses such as drought, heat, cold, salinity and flooding significantly impact plant growth, development and productivity. As the planet has warmed, these abiotic stresses have increased in frequency and intensity, affecting the global food supply and making it imperative to develop stress-resilient crops. In the past 20 years, the development of omics technologies has contributed to the growth of datasets for plants grown under a wide range of abiotic environments. Integration of these rapidly growing data using machine-learning (ML) approaches can complement existing breeding efforts by providing insights into the mechanisms underlying plant responses to stressful conditions, which can be used to guide the design of resilient crops. In this review, we introduce ML approaches and provide examples of how researchers use these approaches to predict molecular activities, gene functions and genotype responses under stressful conditions. Finally, we consider the potential and challenges of using such approaches to enable the design of crops that are better suited to a changing environment.This article is part of the theme issue 'Crops under stress: can we mitigate the impacts of climate change on agriculture and launch the 'Resilience Revolution'?'.

A direct experimental test of Ohno's hypothesis

Gene duplication drives evolution by providing raw material for proteins with novel functions. An influential hypothesis by Ohno (1970) posits that gene duplication helps genes tolerate new mutations and thus facilitates the evolution of new phenotypes. Competing hypotheses argue that deleterious mutations will usually inactivate gene duplicates too rapidly for Ohno's hypothesis to work. We experimentally tested Ohno's hypothesis by evolving one or exactly two copies of a gene encoding a fluorescent protein in Escherichia coli through several rounds of mutation and selection. We analyzed the genotypic and phenotypic evolutionary dynamics of the evolving populations through high-throughput DNA sequencing, biochemical assays, and engineering of selected variants. In support of Ohno's hypothesis, populations carrying two gene copies displayed higher mutational robustness than those carrying a single gene copy. Consequently, the double-copy populations experienced relaxed purifying selection, evolved higher phenotypic and genetic diversity, carried more mutations and accumulated combinations of key beneficial mutations earlier. However, their phenotypic evolution was not accelerated, possibly because one gene copy rapidly became inactivated by deleterious mutations. Our work provides an experimental platform to test models of evolution by gene duplication, and it supports alternatives to Ohno's hypothesis that point to the importance of gene dosage.

Global regulators enable bacterial adaptation to a phenotypic trade-off

Cellular fitness depends on multiple phenotypes that must be balanced during evolutionary adaptation. For instance, coordinating growth and motility is critical for microbial colonization and cancer invasiveness. In bacteria, these phenotypes are controlled by local regulators that target single operons, as well as by global regulators that impact hundreds of genes. However, how the different levels of regulation interact during evolution is unclear. Here, we measured in Escherichia coli how CRISPR-mediated knockdowns of global and local transcription factors impact growth and motility in three environments. We found that local regulators mostly modulate motility, whereas global regulators jointly modulate growth and motility. Simulated evolutionary trajectories indicate that local regulators are typically altered first to improve motility before global regulators adjust growth and motility following their trade-off. These findings highlight the role of pleiotropic regulators in the adaptation of multiple phenotypes.

Data-driven discovery of the interplay between genetic and environmental factors in bacterial growth

A complex interplay of genetic and environmental factors influences bacterial growth. Understanding these interactions is crucial for insights into complex living systems. This study employs a data-driven approach to uncover the principles governing bacterial growth changes due to genetic and environmental variation. A pilot survey is conducted across 115 Escherichia coli strains and 135 synthetic media comprising 45 chemicals, generating 13,944 growth profiles. Machine learning analyzes this dataset to predict the chemicals' priorities for bacterial growth. The primary gene-chemical networks are structured hierarchically, with glucose playing a pivotal role. Offset in bacterial growth changes is frequently observed across 1,445,840 combinations of strains and media, with its magnitude correlating to individual alterations in strains or media. This counterbalance in the gene-chemical interplay is supposed to be a general feature beneficial for bacterial population growth.

Understanding developmental system drift

Developmental system drift (DSD) occurs when the genetic basis for homologous traits diverges over time despite conservation of the phenotype. In this Review, we examine the key ideas, evidence and open problems arising from studies of DSD. Recent work suggests that DSD may be pervasive, having been detected across a range of different organisms and developmental processes. Although developmental research remains heavily reliant on model organisms, extrapolation of findings to non-model organisms can be error-prone if the lineages have undergone DSD. We suggest how existing data and modelling approaches may be used to detect DSD and estimate its frequency. More direct study of DSD, we propose, can inform null hypotheses for how much genetic divergence to expect on the basis of phylogenetic distance, while also contributing to principles of gene regulatory evolution.

Patterns of Fitness and Gene Expression Epistasis Generated by Beneficial Mutations in the rho and rpoB Genes of Escherichia coli during High-Temperature Adaptation

Epistasis is caused by genetic interactions among mutations that affect fitness. To characterize properties and potential mechanisms of epistasis, we engineered eight double mutants that combined mutations from the rho and rpoB genes of Escherichia coli. The two genes encode essential functions for transcription, and the mutations in each gene were chosen because they were beneficial for adaptation to thermal stress (42.2 °C). The double mutants exhibited patterns of fitness epistasis that included diminishing returns epistasis at 42.2 °C, stronger diminishing returns between mutations with larger beneficial effects and both negative and positive (sign) epistasis across environments (20.0 °C and 37.0 °C). By assessing gene expression between single and double mutants, we detected hundreds of genes with gene expression epistasis. Previous work postulated that highly connected hub genes in coexpression networks have low epistasis, but we found the opposite: hub genes had high epistasis values in both coexpression and protein-protein interaction networks. We hypothesized that elevated epistasis in hub genes reflected that they were enriched for targets of Rho termination but that was not the case. Altogether, gene expression and coexpression analyses revealed that thermal adaptation occurred in modules, through modulation of ribonucleotide biosynthetic processes and ribosome assembly, the attenuation of expression in genes related to heat shock and stress responses, and with an overall trend toward restoring gene expression toward the unstressed state.

The distribution of fitness effects during adaptive walks using a simple genetic network

The tempo and mode of adaptation depends on the availability of beneficial alleles. Genetic interactions arising from gene networks can restrict this availability. However, the extent to which networks affect adaptation remains largely unknown. Current models of evolution consider additive genotype-phenotype relationships while often ignoring the contribution of gene interactions to phenotypic variance. In this study, we model a quantitative trait as the product of a simple gene regulatory network, the negative autoregulation motif. Using forward-time genetic simulations, we measure adaptive walks towards a phenotypic optimum in both additive and network models. A key expectation from adaptive walk theory is that the distribution of fitness effects of new beneficial mutations is exponential. We found that both models instead harbored distributions with fewer large-effect beneficial alleles than expected. The network model also had a complex and bimodal distribution of fitness effects among all mutations, with a considerable density at deleterious selection coefficients. This behavior is reminiscent of the cost of complexity, where correlations among traits constrain adaptation. Our results suggest that the interactions emerging from genetic networks can generate complex and multimodal distributions of fitness effects.

Competition-driven eco-evolutionary feedback reshapes bacteriophage lambda's fitness landscape and enables speciation

A major challenge in evolutionary biology is explaining how populations navigate rugged fitness landscapes without getting trapped on local optima. One idea illustrated by adaptive dynamics theory is that as populations adapt, their newly enhanced capacities to exploit resources alter fitness payoffs and restructure the landscape in ways that promote speciation by opening new adaptive pathways. While there have been indirect tests of this theory, to our knowledge none have measured how fitness landscapes deform during adaptation, or test whether these shifts promote diversification. Here, we achieve this by studying bacteriophage [Formula: see text], a virus that readily speciates into co-existing receptor specialists under controlled laboratory conditions. We use a high-throughput gene editing-phenotyping technology to measure [Formula: see text]'s fitness landscape in the presence of different evolved-[Formula: see text] competitors and find that the fitness effects of individual mutations, and their epistatic interactions, depend on the competitor. Using these empirical data, we simulate [Formula: see text]'s evolution on an unchanging landscape and one that recapitulates how the landscape deforms during evolution. [Formula: see text] heterogeneity only evolves in the shifting landscape regime. This study provides a test of adaptive dynamics, and, more broadly, shows how fitness landscapes dynamically change during adaptation, potentiating phenomena like speciation by opening new adaptive pathways.

Competition-driven eco-evolutionary feedback reshapes bacteriophage lambda's fitness landscape and enables speciation

A major challenge in evolutionary biology is explaining how populations navigate rugged fitness landscapes without getting trapped on local optima. One idea illustrated by adaptive dynamics theory is that as populations adapt, their newly enhanced capacities to exploit resources alter fitness payoffs and restructure the landscape in ways that promote speciation by opening new adaptive pathways. While there have been indirect tests of this theory, none have measured how fitness landscapes deform during adaptation, or test whether these shifts promote diversification. Here, we achieve this by studying bacteriophage λ, a virus that readily speciates into co-existing receptor specialists under controlled laboratory conditions. We used a high-throughput gene editing-phenotyping technology to measure λ's fitness landscape in the presence of different evolved-λ competitors and found that the fitness effects of individual mutations, and their epistatic interactions, depend on the competitor. Using these empirical data, we simulated λ's evolution on an unchanging landscape and one that recapitulates how the landscape deforms during evolution. λ heterogeneity only evolved in the shifting landscape regime. This study provides a test of adaptive dynamics, and, more broadly, shows how fitness landscapes dynamically change during adaptation, potentiating phenomena like speciation by opening new adaptive pathways.

Improving Genetic Association Studies with a Novel Methodology that Unveils the Hidden Complexity of All-Cause Heart Failure

Motivation

Genome-Wide Association Studies (GWAS) commonly assume phenotypic and genetic homogeneity that is not present in complex conditions. We designed Transformative Regression Analysis of Combined Effects (TRACE), a GWAS methodology that better accounts for clinical phenotype heterogeneity and identifies gene-by-environment (GxE) interactions. We demonstrated with UK Biobank (UKB) data that TRACE increased the variance explained in All-Cause Heart Failure (AHF) via the discovery of novel single nucleotide polymorphism (SNP) and SNP-by-environment (i.e. GxE) interaction associations. First, we transformed 312 AHF-related ICD10 codes (including AHF) into continuous low-dimensional features (i.e., latent phenotypes) for a more nuanced disease representation. Then, we ran a standard GWAS on our latent phenotypes to discover main effects and identified GxE interactions with target encoding. Genes near associated SNPs subsequently underwent enrichment analysis to explore potential functional mechanisms underlying associations. Latent phenotypes were regressed against their SNP hits and the estimated latent phenotype values were used to measure the amount of AHF variance explained.

Results

Our method identified over 100 main GWAS effects that were consistent with prior studies and hundreds of novel gene-by-smoking interactions, which collectively accounted for approximately 10% of AHF variance. This represents an improvement over traditional GWAS whose results account for a negligible proportion of AHF variance. Enrichment analyses suggested that hundreds of miRNAs mediated the SNP effect on various AHF-related biological pathways. The TRACE framework can be applied to decode the genetics of other complex diseases.

Availability

All code is available at https://github.com/EpistasisLab/latent_phenotype_project.