Extent and context dependence of pleiotropy revealed by high-throughput single-cell phenotyping

Low-dimensional genotype-fitness mapping across divergent environments suggests a limiting functions model of fitness

A central goal in evolutionary biology is to be able to predict the effect of a genetic mutation on fitness. This is a major challenge because fitness depends both on phenotypic changes due to the mutation, and how these phenotypes map onto fitness in a particular environment. Genotype, phenotype, and environment spaces are all extremely complex, rendering bottom-up prediction unlikely. Here we show, using a large collection of adaptive yeast mutants, that fitness across a set of lab environments can be well-captured by top-down, low-dimensional linear models that generate abstract genotype-phenotype-fitness maps. We find that these maps are low-dimensional not only in the environment where the adaptive mutants evolved, but also in more divergent environments. We further find that the genotype-phenotype-fitness spaces implied by these maps overlap only partially across environments. We argue that these patterns are consistent with a "limiting functions" model of fitness, whereby only a small number of limiting functions can be modified to affect fitness in any given environment. The pleiotropic side-effects on non-limiting functions are effectively hidden from natural selection locally, but can be revealed globally. These results combine to emphasize the importance of environmental context in genotype-phenotype-fitness mapping, and have implications for the predictability and trajectory of evolution in complex environments.

Disentangling variational bias: the roles of development, mutation, and selection

The extraordinary diversity and adaptive fit of organisms to their environment depends fundamentally on the availability of variation. While most population genetic frameworks assume that random mutations produce isotropic phenotypic variation, the distribution of variation available to natural selection is more restricted, as the distribution of phenotypic variation is affected by a range of factors in developmental systems. Here, we revisit the concept of developmental bias - the observation that the generation of phenotypic variation is biased due to the structure, character, composition, or dynamics of the developmental system - and argue that a more rigorous investigation into the role of developmental bias in the genotype-to-phenotype map will produce fundamental insights into evolutionary processes, with potentially important consequences on the relation between micro- and macro-evolution. We discuss the hierarchical relationships between different types of variational biases, including mutation bias and developmental bias, and their roles in shaping the realized phenotypic space. Furthermore, we highlight the challenges in studying variational bias and propose potential approaches to identify developmental bias using modern tools.

Distinguishing mutants that resist drugs via different mechanisms by examining fitness tradeoffs

There is growing interest in designing multidrug therapies that leverage tradeoffs to combat resistance. Tradeoffs are common in evolution and occur when, for example, resistance to one drug results in sensitivity to another. Major questions remain about the extent to which tradeoffs are reliable, specifically, whether the mutants that provide resistance to a given drug all suffer similar tradeoffs. This question is difficult because the drug-resistant mutants observed in the clinic, and even those evolved in controlled laboratory settings, are often biased towards those that provide large fitness benefits. Thus, the mutations (and mechanisms) that provide drug resistance may be more diverse than current data suggests. Here, we perform evolution experiments utilizing lineage-tracking to capture a fuller spectrum of mutations that give yeast cells a fitness advantage in fluconazole, a common antifungal drug. We then quantify fitness tradeoffs for each of 774 evolved mutants across 12 environments, finding these mutants group into classes with characteristically different tradeoffs. Their unique tradeoffs may imply that each group of mutants affects fitness through different underlying mechanisms. Some of the groupings we find are surprising. For example, we find some mutants that resist single drugs do not resist their combination, while others do. And some mutants to the same gene have different tradeoffs than others. These findings, on one hand, demonstrate the difficulty in relying on consistent or intuitive tradeoffs when designing multidrug treatments. On the other hand, by demonstrating that hundreds of adaptive mutations can be reduced to a few groups with characteristic tradeoffs, our findings may yet empower multidrug strategies that leverage tradeoffs to combat resistance. More generally speaking, by grouping mutants that likely affect fitness through similar underlying mechanisms, our work guides efforts to map the phenotypic effects of mutation.

FLT4 causes developmental disorders of the cardiovascular and lymphovascular systems via pleiotropic molecular mechanisms

AIMS

Rare, deleterious genetic variants in FLT4 are associated with Tetralogy of Fallot (TOF), the most common cyanotic congenital heart disease. The distinct genetic variants in FLT4 are also an established cause of Milroy disease, the most prevalent form of primary hereditary lymphoedema. The phenotypic features of these two conditions are non-overlapping, implying pleiotropic cellular mechanisms during development.

METHODS AND RESULTS

In this study, we show that FLT4 variants identified in patients with TOF, when expressed in primary human endothelial cells, cause aggregation of FLT4 protein in the perinuclear endoplasmic reticulum, activating proteostatic and metabolic signalling, whereas lymphoedema-associated FLT4 variants and wild-type (WT) FLT4 do not. FLT4 TOF variants display characteristic gene expression profiles in key developmental signalling pathways, revealing a role for FLT4 in cardiogenesis distinct from its role in lymphatic development. Inhibition of proteostatic signalling abrogates these effects, identifying potential avenues for therapeutic intervention. Depletion of flt4 in zebrafish caused cardiac phenotypes of reduced heart size and altered heart looping. These phenotypes were rescued with coinjection of WT human FLT4 mRNA, but incompletely or not at all by mRNA harbouring FLT4 TOF variants.

CONCLUSION

Taken together, we identify a pathogenic mechanism for FLT4 variants predisposing to TOF that is distinct from the known dominant negative mechanism of Milroy-causative variants. FLT4 variants give rise to conditions of the two circulatory subdivisions of the vascular system via distinct developmental pleiotropic molecular mechanisms.

Distinguishing mutants that resist drugs via different mechanisms by examining fitness tradeoffs

There is growing interest in designing multidrug therapies that leverage tradeoffs to combat resistance. Tradeoffs are common in evolution and occur when, for example, resistance to one drug results in sensitivity to another. Major questions remain about the extent to which tradeoffs are reliable, specifically, whether the mutants that provide resistance to a given drug all suffer similar tradeoffs. This question is difficult because the drug-resistant mutants observed in the clinic, and even those evolved in controlled laboratory settings, are often biased towards those that provide large fitness benefits. Thus, the mutations (and mechanisms) that provide drug resistance may be more diverse than current data suggests. Here, we perform evolution experiments utilizing lineage-tracking to capture a fuller spectrum of mutations that give yeast cells a fitness advantage in fluconazole, a common antifungal drug. We then quantify fitness tradeoffs for each of 774 evolved mutants across 12 environments, finding these mutants group into 6 classes with characteristically different tradeoffs. Their unique tradeoffs may imply that each group of mutants affects fitness through different underlying mechanisms. Some of the groupings we find are surprising. For example, we find some mutants that resist single drugs do not resist their combination, while others do. And some mutants to the same gene have different tradeoffs than others. These findings, on one hand, demonstrate the difficulty in relying on consistent or intuitive tradeoffs when designing multidrug treatments. On the other hand, by demonstrating that hundreds of adaptive mutations can be reduced to a few groups with characteristic tradeoffs, our findings may yet empower multidrug strategies that leverage tradeoffs to combat resistance. More generally speaking, by grouping mutants that likely affect fitness through similar underlying mechanisms, our work guides efforts to map the phenotypic effects of mutation.

Sex-specific associations of the maturation locus vgll3 with exploratory behavior and boldness in Atlantic salmon juveniles

Studies linking genetics, behavior and life history in any species are rare. In Atlantic salmon (Salmo salar), age at maturity is a key life-history trait and associates strongly with the vgll3 locus, whereby the vgll3*E allele is linked with younger age at maturity, and higher body condition than the vgll3*L allele. However, the relationship between this genetic variation and behaviors like boldness and exploration which may impact growth and reproductive strategies is poorly understood. The pace-of-life syndrome (POLS) framework provides predictions, whereby heightened exploratory behavior and boldness are predicted in individuals with the early maturation-associated vgll3 genotype (EE). Here, we tested these predictions by investigating the relationship between vgll3 genotypes and exploration and boldness behaviors in 129 juveniles using the novel environment and novel object trials. Our results indicated that contrary to POLS predictions, vgll3*LL fish were bolder and more explorative, suggesting a genotype-level syndrome including several behaviors. Interestingly, clear sex differences were observed in the latency to move in a new environment, with vgll3*EE males, but not females, taking longer to move than their vgll3*LL counterparts. Our results provide further empirical support for recent calls to consider more nuanced explanations than the pace of life theory for integrating behavior into life-history theory.

Application of unimodal probability distribution models for morphological phenotyping of budding yeast

Morphological phenotyping of the budding yeast Saccharomyces cerevisiae has helped to greatly clarify the functions of genes and increase our understanding of cellular functional networks. It is necessary to understand cell morphology and perform quantitative morphological analysis (QMA) but assigning precise values to morphological phenotypes has been challenging. We recently developed the Unimodal Morphological Data image analysis pipeline for this purpose. All true values can be estimated theoretically by applying an appropriate probability distribution if the distribution of experimental values follows a unimodal pattern. This reliable pipeline allows several downstream analyses, including detection of subtle morphological differences, selection of mutant strains with similar morphology, clustering based on morphology, and study of morphological diversity. In addition to basic research, morphological analyses of yeast cells can also be used in applied research to monitor breeding and fermentation processes and control the fermentation activity of yeast cells.

Joint analysis of multiple phenotypes for extremely unbalanced case-control association studies using multi-layer network

MOTIVATION

Genome-wide association studies is an essential tool for analyzing associations between phenotypes and single nucleotide polymorphisms (SNPs). Most of binary phenotypes in large biobanks are extremely unbalanced, which leads to inflated type I error rates for many widely used association tests for joint analysis of multiple phenotypes. In this article, we first propose a novel method to construct a Multi-Layer Network (MLN) using individuals with at least one case status among all phenotypes. Then, we introduce a computationally efficient community detection method to group phenotypes into disjoint clusters based on the MLN. Finally, we propose a novel approach, MLN with Omnibus (MLN-O), to jointly analyse the association between phenotypes and a SNP. MLN-O uses the score test to test the association of each merged phenotype in a cluster and a SNP, then uses the Omnibus test to obtain an overall test statistic to test the association between all phenotypes and a SNP.

RESULTS

We conduct extensive simulation studies to reveal that the proposed approach can control type I error rates and is more powerful than some existing methods. Meanwhile, we apply the proposed method to a real data set in the UK Biobank. Using phenotypes in Chapter XIII (Diseases of the musculoskeletal system and connective tissue) in the UK Biobank, we find that MLN-O identifies more significant SNPs than other methods we compare with.

AVAILABILITY AND IMPLEMENTATION

https://github.com/Hongjing-Xie/Multi-Layer-Network-with-Omnibus-MLN-O.

Patterns and evolutionary consequences of pleiotropy

Pleiotropy refers to the phenomenon of one gene or one mutation affecting multiple phenotypic traits. While the concept of pleiotropy is as old as Mendelian genetics, functional genomics has finally allowed the first glimpses of the extent of pleiotropy for a large fraction of genes in a genome. After describing conceptual and operational difficulties in quantifying pleiotropy and the pros and cons of various methods for measuring pleiotropy, I review empirical data on pleiotropy, which generally show an L-shaped distribution of the degree of pleiotropy (i.e., the number of traits affected) with most genes having low pleiotropy. I then review the current understanding of the molecular basis of pleiotropy. The rest of the review discusses evolutionary consequences of pleiotropy, focusing on advances in topics including the cost of complexity, regulatory vs. coding evolution, environmental pleiotropy and adaptation, evolution of ageing and other seemingly harmful traits, and evolutionary resolution of pleiotropy.

A pleiotropic-epistatic entangelement model of drug response

Because drug response is multifactorial, graph models are uniquely powerful for comprehending its genetic architecture. We deconstruct drug response into many different and interdependent sub-traits, with each sub-trait controlled by multiple genes that act and interact in a complicated manner. The outcome of drug response is the consequence of multileveled intertwined interactions between pleiotropic effects and epistatic effects. Here, we propose a general statistical physics framework to chart the 3D geometric network that codes how epistasis pleiotropically influences a complete set of sub-traits to shape body-drug interactions. This model can dissect the topological architecture of epistatically induced pleiotropic networks (EiPN) and pleiotropically influenced epistatic networks (PiEN). We analyze and interpret the practical implications of the pleiotropic-epistatic entanglement model for pharmacogenomic studies.

Detecting natural selection in trait-trait coevolution

No phenotypic trait evolves independently of all other traits, but the cause of trait-trait coevolution is poorly understood. While the coevolution could arise simply from pleiotropic mutations that simultaneously affect the traits concerned, it could also result from multivariate natural selection favoring certain trait relationships. To gain a general mechanistic understanding of trait-trait coevolution, we examine the evolution of 220 cell morphology traits across 16 natural strains of the yeast Saccharomyces cerevisiae and the evolution of 24 wing morphology traits across 110 fly species of the family Drosophilidae, along with the variations of these traits among gene deletion or mutation accumulation lines (a.k.a. mutants). For numerous trait pairs, the phenotypic correlation among evolutionary lineages differs significantly from that among mutants. Specifically, we find hundreds of cases where the evolutionary correlation between traits is strengthened or reversed relative to the mutational correlation, which, according to our population genetic simulation, is likely caused by multivariate selection. Furthermore, we detect selection for enhanced modularity of the yeast traits analyzed. Together, these results demonstrate that trait-trait coevolution is shaped by natural selection and suggest that the pleiotropic structure of mutation is not optimal. Because the morphological traits analyzed here are chosen largely because of their measurability and thereby are not expected to be biased with regard to natural selection, our conclusion is likely general.

Extreme Sensitivity of Fitness to Environmental Conditions: Lessons from #1BigBatch

The phrase "survival of the fittest" has become an iconic descriptor of how natural selection works. And yet, precisely measuring fitness, even for single-celled microbial populations growing in controlled laboratory conditions, remains a challenge. While numerous methods exist to perform these measurements, including recently developed methods utilizing DNA barcodes, all methods are limited in their precision to differentiate strains with small fitness differences. In this study, we rule out some major sources of imprecision, but still find that fitness measurements vary substantially from replicate to replicate. Our data suggest that very subtle and difficult to avoid environmental differences between replicates create systematic variation across fitness measurements. We conclude by discussing how fitness measurements should be interpreted given their extreme environment dependence. This work was inspired by the scientific community who followed us and gave us tips as we live tweeted a high-replicate fitness measurement experiment at #1BigBatch.

In Silico Pleiotropy Analysis in KEGG Signaling Networks Using a Boolean Network Model

Pleiotropy, which refers to the ability of different mutations on the same gene to cause different pathological effects in human genetic diseases, is important in understanding system-level biological diseases. Although some biological experiments have been proposed, still little is known about pleiotropy on gene–gene dynamics, since most previous studies have been based on correlation analysis. Therefore, a new perspective is needed to investigate pleiotropy in terms of gene–gene dynamical characteristics. To quantify pleiotropy in terms of network dynamics, we propose a measure called in silico Pleiotropic Scores (sPS), which represents how much a gene is affected against a pair of different types of mutations on a Boolean network model. We found that our model can identify more candidate pleiotropic genes that are not known to be pleiotropic than the experimental database. In addition, we found that many types of functionally important genes tend to have higher sPS values than other genes; in other words, they are more pleiotropic. We investigated the relations of sPS with the structural properties in the signaling network and found that there are highly positive relations to degree, feedback loops, and centrality measures. This implies that the structural characteristics are principles to identify new pleiotropic genes. Finally, we found some biological evidence showing that sPS analysis is relevant to the real pleiotropic data and can be considered a novel candidate for pleiotropic gene research. Taken together, our results can be used to understand the dynamics pleiotropic characteristics in complex biological systems in terms of gene–phenotype relations.

Maintenance of quantitative genetic variance in complex, multi-trait phenotypes: The contribution of rare, large effect variants in two Drosophila species

The interaction of evolutionary processes to determine quantitative genetic variation has implications for contemporary and future phenotypic evolution, as well as for our ability to detect causal genetic variants. While theoretical studies have provided robust predictions to discriminate among competing models, empirical assessment of these has been limited. In particular, theory highlights the importance of pleiotropy in resolving observations of selection and mutation, but empirical investigations have typically been limited to few traits. Here, we applied high-dimensional Bayesian Sparse Factor Genetic modeling to gene expression datasets in 2 species, Drosophila melanogaster and Drosophila serrata, to explore the distributions of genetic variance across high-dimensional phenotypic space. Surprisingly, most of the heritable trait covariation was due to few lines (genotypes) with extreme [>3 interquartile ranges (IQR) from the median] values. Intriguingly, while genotypes extreme for a multivariate factor also tended to have a higher proportion of individual traits that were extreme, we also observed genotypes that were extreme for multivariate factors but not for any individual trait. We observed other consistent differences between heritable multivariate factors with outlier lines vs those factors without extreme values, including differences in gene functions. We use these observations to identify further data required to advance our understanding of the evolutionary dynamics and nature of standing genetic variation for quantitative traits.

Challenges and potential solutions for studying the genetic and phenotypic architecture of adaptation in microbes

All organisms are defined by the makeup of their DNA. Over billions of years, the structure and information contained in that DNA, often referred to as genetic architecture, have been honed by a multitude of evolutionary processes. Mutations that cause genetic elements to change in a way that results in beneficial phenotypic change are more likely to survive and propagate through the population in a process known as adaptation. Recent work reveals that the genetic targets of adaptation are varied and can change with genetic background. Further, seemingly similar adaptive mutations, even within the same gene, can have diverse and unpredictable effects on phenotype. These challenges represent major obstacles in predicting adaptation and evolution. In this review, we cover these concepts in detail and identify three emerging synergistic solutions: higher-throughput evolution experiments combined with updated genotype-phenotype mapping strategies and physiological models. Our review largely focuses on recent literature in yeast, and the field seems to be on the cusp of a new era with regard to studying the predictability of evolution.

The acquisitive-conservative axis of leaf trait variation emerges even in homogeneous environments

BACKGROUND AND AIMS

The acquisitive-conservative axis of plant ecological strategies results in a pattern of leaf trait covariation that captures the balance between leaf construction costs and plant growth potential. Studies evaluating trait covariation within species are scarcer, and have mostly dealt with variation in response to environmental gradients. Little work has been published on intraspecific patterns of leaf trait covariation in the absence of strong environmental variation.

METHODS

We analysed covariation of four leaf functional traits [specific leaf area (SLA) leaf dry matter content (LDMC), force to tear (Ft) and leaf nitrogen content (Nm)] in six Poaceae and four Fabaceae species common in the dry Chaco forest of Central Argentina, growing in the field and in a common garden. We compared intraspecific covariation patterns (slopes, correlation and effect size) of leaf functional traits with global interspecific covariation patterns. Additionally, we checked for possible climatic and edaphic factors that could affect the intraspecific covariation pattern.

KEY RESULTS

We found negative correlations for the LDMC-SLA, Ft-SLA, LDMC-Nm and Ft-Nm trait pairs. This intraspecific covariation pattern found both in the field and in the common garden and not explained by climatic or edaphic variation in the field follows the expected acquisitive-conservative axis. At the same time, we found quantitative differences in slopes among different species, and between these intraspecific patterns and the interspecific ones. Many of these differences seem to be idiosyncratic, but some appear consistent among species (e.g. all the intraspecific LDMC-SLA and LDMC-Nm slopes tend to be shallower than the global pattern).

CONCLUSIONS

Our study indicates that the acquisitive-conservative leaf functional trait covariation pattern occurs at the intraspecific level even in the absence of relevant environmental variation in the field. This suggests a high degree of variation-covariation in leaf functional traits not driven by environmental variables.

Shared Genetic Etiology between Cortical Brain Morphology and Tobacco, Alcohol, and Cannabis Use

Genome-wide association studies (GWAS) have identified genetic variants associated with brain morphology and substance use behaviors (SUB). However, the genetic overlap between brain structure and SUB has not been well characterized. We leveraged GWAS summary data of 71 brain imaging measures and alcohol, tobacco, and cannabis use to investigate their genetic overlap using linkage disequilibrium score regression. We used genomic structural equation modeling to model a "common SUB genetic factor" and investigated its genetic overlap with brain structure. Furthermore, we estimated SUB polygenic risk scores (PRS) and examined whether they predicted brain imaging traits using the Adolescent Behavior and Cognitive Development (ABCD) study. We identified 8 significant negative genetic correlations, including between (1) alcoholic drinks per week and average cortical thickness, and (2) intracranial volume with age of smoking initiation. We observed 5 positive genetic correlations, including those between (1) insula surface area and lifetime cannabis use, and (2) the common SUB genetic factor and pericalcarine surface area. SUB PRS were associated with brain structure variation in ABCD. Our findings highlight a shared genetic etiology between cortical brain morphology and SUB and suggest that genetic variants associated with SUB may be causally related to brain structure differences.

High-throughput platform for yeast morphological profiling predicts the targets of bioactive compounds

Morphological profiling is an omics-based approach for predicting intracellular targets of chemical compounds in which the dose-dependent morphological changes induced by the compound are systematically compared to the morphological changes in gene-deleted cells. In this study, we developed a reliable high-throughput (HT) platform for yeast morphological profiling using drug-hypersensitive strains to minimize compound use, HT microscopy to speed up data generation and analysis, and a generalized linear model to predict targets with high reliability. We first conducted a proof-of-concept study using six compounds with known targets: bortezomib, hydroxyurea, methyl methanesulfonate, benomyl, tunicamycin, and echinocandin B. Then we applied our platform to predict the mechanism of action of a novel diferulate-derived compound, poacidiene. Morphological profiling of poacidiene implied that it affects the DNA damage response, which genetic analysis confirmed. Furthermore, we found that poacidiene inhibits the growth of phytopathogenic fungi, implying applications as an effective antifungal agent. Thus, our platform is a new whole-cell target prediction tool for drug discovery.

Transcriptomics and Functional Analysis of Copper Stress Response in the Sulfate-Reducing Bacterium Desulfovibrio alaskensis G20

Copper (Cu) is an essential micronutrient required as a co-factor in the catalytic center of many enzymes. However, excess Cu can generate pleiotropic effects in the microbial cell. In addition, leaching of Cu from pipelines results in elevated Cu concentration in the environment, which is of public health concern. Sulfate-reducing bacteria (SRB) have been demonstrated to grow in toxic levels of Cu. However, reports on Cu toxicity towards SRB have primarily focused on the degree of toxicity and subsequent elimination. Here, Cu(II) stress-related effects on a model SRB, Desulfovibrio alaskensis G20, is reported. Cu(II) stress effects were assessed as alterations in the transcriptome through RNA-Seq at varying Cu(II) concentrations (5 µM and 15 µM). In the pairwise comparison of control vs. 5 µM Cu(II), 61.43% of genes were downregulated, and 38.57% were upregulated. In control vs. 15 µM Cu(II), 49.51% of genes were downregulated, and 50.5% were upregulated. The results indicated that the expression of inorganic ion transporters and translation machinery was massively modulated. Moreover, changes in the expression of critical biological processes such as DNA transcription and signal transduction were observed at high Cu(II) concentrations. These results will help us better understand the Cu(II) stress-response mechanism and provide avenues for future research.

A Multilayer Interactome Network Constructed in a Forest Poplar Population Mediates the Pleiotropic Control of Complex Traits

The effects of genes on physiological and biochemical processes are interrelated and interdependent; it is common for genes to express pleiotropic control of complex traits. However, the study of gene expression and participating pathways in vivo at the whole-genome level is challenging. Here, we develop a coupled regulatory interaction differential equation to assess overall and independent genetic effects on trait growth. Based on evolutionary game theory and developmental modularity theory, we constructed multilayer, omnigenic networks of bidirectional, weighted, and positive or negative epistatic interactions using a forest poplar tree mapping population, which were organized into metagalactic, intergalactic, and local interstellar networks that describe layers of structure between modules, submodules, and individual single nucleotide polymorphisms, respectively. These multilayer interactomes enable the exploration of complex interactions between genes, and the analysis of not only differential expression of quantitative trait loci but also previously uncharacterized determinant SNPs, which are negatively regulated by other SNPs, based on the deconstruction of genetic effects to their component parts. Our research framework provides a tool to comprehend the pleiotropic control of complex traits and explores the inherent directional connections between genes in the structure of omnigenic networks.

Dynamics and variability in the pleiotropic effects of adaptation in laboratory budding yeast populations

Evolutionary adaptation to a constant environment is driven by the accumulation of mutations which can have a range of unrealized pleiotropic effects in other environments. These pleiotropic consequences of adaptation can influence the emergence of specialists or generalists, and are critical for evolution in temporally or spatially fluctuating environments. While many experiments have examined the pleiotropic effects of adaptation at a snapshot in time, very few have observed the dynamics by which these effects emerge and evolve. Here, we propagated hundreds of diploid and haploid laboratory budding yeast populations in each of three environments, and then assayed their fitness in multiple environments over 1000 generations of evolution. We find that replicate populations evolved in the same condition share common patterns of pleiotropic effects across other environments, which emerge within the first several hundred generations of evolution. However, we also find dynamic and environment-specific variability within these trends: variability in pleiotropic effects tends to increase over time, with the extent of variability depending on the evolution environment. These results suggest shifting and overlapping contributions of chance and contingency to the pleiotropic effects of adaptation, which could influence evolutionary trajectories in complex environments that fluctuate across space and time.

Pleiotropy data resource as a primer for investigating co-morbidities/multi-morbidities and their role in disease

Most current biomedical and protein research focuses only on a small proportion of genes, which results in a lost opportunity to identify new gene-disease associations and explore new opportunities for therapeutic intervention. The International Mouse Phenotyping Consortium (IMPC) focuses on elucidating gene function at scale for poorly characterized and/or under-studied genes. A key component of the IMPC initiative is the implementation of a broad phenotyping pipeline, which is facilitating the discovery of pleiotropy. Characterizing pleiotropy is essential to identify gene-disease associations, and it is of particular importance when elucidating the genetic causes of syndromic disorders. Here we show how the IMPC is effectively uncovering pleiotropy and how the new mouse models and gene function hypotheses generated by the IMPC are increasing our understanding of the mammalian genome, forming the basis of new research and identifying new gene-disease associations.

Fitness variation across subtle environmental perturbations reveals local modularity and global pleiotropy of adaptation

Building a genotype-phenotype-fitness map of adaptation is a central goal in evolutionary biology. It is difficult even when adaptive mutations are known because it is hard to enumerate which phenotypes make these mutations adaptive. We address this problem by first quantifying how the fitness of hundreds of adaptive yeast mutants responds to subtle environmental shifts. We then model the number of phenotypes these mutations collectively influence by decomposing these patterns of fitness variation. We find that a small number of inferred phenotypes can predict fitness of the adaptive mutations near their original glucose-limited evolution condition. Importantly, inferred phenotypes that matter little to fitness at or near the evolution condition can matter strongly in distant environments. This suggests that adaptive mutations are locally modular - affecting a small number of phenotypes that matter to fitness in the environment where they evolved - yet globally pleiotropic - affecting additional phenotypes that may reduce or improve fitness in new environments.