Gene functional trade-offs and the evolution of pleiotropy

Pleiotropy alleviates the fitness costs associated with resource allocation trade-offs in immune signalling networks

Many genes and signalling pathways within plant and animal taxa drive the expression of multiple organismal traits. This form of genetic pleiotropy instigates trade-offs among life-history traits if a mutation in the pleiotropic gene improves the fitness contribution of one trait at the expense of another. Whether or not pleiotropy gives rise to conflict among traits, however, likely depends on the resource costs and timing of trait deployment during organismal development. To investigate factors that could influence the evolutionary maintenance of pleiotropy in gene networks, we developed an agent-based model of co-evolution between parasites and hosts. Hosts comprise signalling networks that must faithfully complete a developmental programme while also defending against parasites, and trait signalling networks could be independent or share a pleiotropic component as they evolved to improve host fitness. We found that hosts with independent developmental and immune networks were significantly more fit than hosts with pleiotropic networks when traits were deployed asynchronously during development. When host genotypes directly competed against each other, however, pleiotropic hosts were victorious regardless of trait synchrony because the pleiotropic networks were more robust to parasite manipulation, potentially explaining the abundance of pleiotropy in immune systems despite its contribution to life history trade-offs.

From molecular variations to behavioral adaptations: Unveiling adaptive epistasis in primate oxytocin system

OBJECTIVE

Our primary objective was to investigate the variability of oxytocin (OT) and the GAMEN binding motif within the LNPEP oxytocinase in primates.

MATERIALS AND METHODS

We sequenced the LNPEP segment encompassing the GAMEN motif in 34 Platyrrhini species, with 21 of them also sequenced for the OT gene. Our dataset was supplemented with primate sequences of LNPEP, OT, and the oxytocin receptor (OTR) sourced from public databases. Evolutionary analysis and coevolution predictions were made followed by the macroevolution analysis of relevant amino acids associated with phenotypic traits, such as mating systems, parental care, and litter size. To account for phylogenetic structure, we utilized two distinct statistical tests. Additionally, we calculated binding energies focusing on the interaction between Callithtrix jacchus VAMEN and Pro8OT.

RESULTS

We identified two novel motifs (AAMEN and VAMEN), challenging the current knowledge of motif conservation in placental mammals. Coevolution analysis demonstrated a correlation between GAMEN, AAMEN, and VAMEN and their corresponding OTs and OTRs. Callithrix jacchus exhibited a higher binding energy between VAMEN and Pro8OT than orthologous molecules found in humans (GAMEN and Leu8OT).

DISCUSSION

The coevolution of AAMEN and VAMEN with their corresponding OTs and OTRs suggests a functional relationship that could have contributed to specific reproductive and adaptive behaviors, including paternal care, social monogamy, and twin births, prominent traits in Cebidae species, such as marmosets and tamarins. Our findings underscore the coevolution of taxon-specific amino acids among the three studied molecules, shedding light on the oxytocinergic system as an adaptive epistatic repertoire in primates.

Evolutionarily new genes in humans with disease phenotypes reveal functional enrichment patterns shaped by adaptive innovation and sexual selection

New genes (or young genes) are genetic novelties pivotal in mammalian evolution. Their phenotypic impacts and evolutionary pattern over time, however, remain elusive in humans due to the technical and ethical complexities in functional studies. By combining human gene age dating and Mendelian disease phenotyping, our research reveals a gradual increase in disease gene proportions with gene age. Logistic regression modeling indicates that this increase could be related to longer protein lengths and higher burdens of deleterious de novo germline variants (DNVs) for older genes. We also find a steady integration of new genes with biomedical phenotypes into the human genome over macroevolutionary timescales (~0.07% per million years). Despite this stable pace, we observe distinct patterns in phenotypic enrichment, pleiotropy, and selective pressures across gene ages. Notably, young genes show significant enrichment in diseases related to the male reproductive system, indicating strong sexual selection. Young genes also exhibit disease-related functions in tissues and systems potentially linked to human phenotypic innovations, such as increased brain size, musculoskeletal phenotypes, and color vision. We further reveal a logistic growth pattern of pleiotropy over evolutionary time, indicating a diminishing marginal growth of new functions for older genes due to intensifying selective constraints over time. We propose a "pleiotropy-barrier" model that delineates higher potentials of phenotypic innovation for young genes than for older genes, a process subject to natural selection. Our study demonstrates that evolutionary new genes are critical in influencing human reproductive evolution and adaptive phenotypic innovations driven by sexual and natural selection, with low pleiotropy as a selective advantage.

Asymmetric and parallel subgenome selection co-shape common carp domestication

BACKGROUND

The common carp (Cyprinus carpio) might best represent the domesticated allopolyploid animals. Although subgenome divergence which is well-known to be a key to allopolyploid domestication has been comprehensively characterized in common carps, the link between genetic architecture underlying agronomic traits and subgenome divergence is unknown in the selective breeding of common carps globally.

RESULTS

We utilized a comprehensive SNP dataset in 13 representative common carp strains worldwide to detect genome-wide genetic variations associated with scale reduction, vibrant skin color, and high growth rate in common carp domestication. We identified numerous novel candidate genes underlie the three agronomically most desirable traits in domesticated common carps, providing potential molecular targets for future genetic improvement in the selective breeding of common carps. We found that independently selective breeding of the same agronomic trait (e.g., fast growing) in common carp domestication could result from completely different genetic variations, indicating the potential advantage of allopolyploid in domestication. We observed that candidate genes associated with scale reduction, vibrant skin color, and/or high growth rate are repeatedly enriched in the immune system, suggesting that domestication of common carps was often accompanied by the disease resistance improvement.

CONCLUSIONS

In common carp domestication, asymmetric subgenome selection is prevalent, while parallel subgenome selection occurs in selective breeding of common carps. This observation is not due to asymmetric gene retention/loss between subgenomes but might be better explained by reduced pleiotropy through transposable element-mediated expression divergence between ohnologs. Our results demonstrate that domestication benefits from polyploidy not only in plants but also in animals.

Evolutionarily new genes in humans with disease phenotypes reveal functional enrichment patterns shaped by adaptive innovation and sexual selection

New genes (or young genes) are structural novelties pivotal in mammalian evolution. Their phenotypic impact on humans, however, remains elusive due to the technical and ethical complexities in functional studies. Through combining gene age dating with Mendelian disease phenotyping, our research reveals that new genes associated with disease phenotypes steadily integrate into the human genome at a rate of ~ 0.07% every million years over macroevolutionary timescales. Despite this stable pace, we observe distinct patterns in phenotypic enrichment, pleiotropy, and selective pressures between young and old genes. Notably, young genes show significant enrichment in the male reproductive system, indicating strong sexual selection. Young genes also exhibit functions in tissues and systems potentially linked to human phenotypic innovations, such as increased brain size, bipedal locomotion, and color vision. Our findings further reveal increasing levels of pleiotropy over evolutionary time, which accompanies stronger selective constraints. We propose a "pleiotropy-barrier" model that delineates different potentials for phenotypic innovation between young and older genes subject to natural selection. Our study demonstrates that evolutionary new genes are critical in influencing human reproductive evolution and adaptive phenotypic innovations driven by sexual and natural selection, with low pleiotropy as a selective advantage.

Pleiotropy alleviates the fitness costs associated with resource allocation trade-offs in immune signaling networks

Many genes and signaling pathways within plant and animal taxa drive the expression of multiple organismal traits. This form of genetic pleiotropy instigates trade-offs among life-history traits if a mutation in the pleiotropic gene improves the fitness contribution of one trait at the expense of another. Whether or not pleiotropy gives rise to conflict among traits, however, likely depends on the resource costs and timing of trait deployment during organismal development. To investigate factors that could influence the evolutionary maintenance of pleiotropy in gene networks, we developed an agent-based model of co-evolution between parasites and hosts. Hosts comprise signaling networks that must faithfully complete a developmental program while also defending against parasites, and trait signaling networks could be independent or share a pleiotropic component as they evolved to improve host fitness. We found that hosts with independent developmental and immune networks were significantly more fit than hosts with pleiotropic networks when traits were deployed asynchronously during development. When host genotypes directly competed against each other, however, pleiotropic hosts were victorious regardless of trait synchrony because the pleiotropic networks were more robust to parasite manipulation, potentially explaining the abundance of pleiotropy in immune systems despite its contribution to life history trade-offs.

Bivariate GWAS reveals pleiotropic regions among feed efficiency and beef quality-related traits in Nelore cattle

Feed-efficient cattle selection is among the most leading solutions to reduce cost for beef cattle production. However, technical difficulties in measuring feed efficiency traits had limited the application in livestock. Here, we performed a Bivariate Genome-Wide Association Study (Bi-GWAS) and presented candidate biological mechanisms underlying the association between feed efficiency and meat quality traits in a half-sibling design with 353 Nelore steers derived from 34 unrelated sires. A total of 13 Quantitative Trait Loci (QTL) were found explaining part of the phenotypic variations. An important transcription factor of adipogenesis in cattle, the TAL1 (rs133408775) gene located on BTA3 was associated with intramuscular fat and average daily gain (IMF-ADG), and a region located on BTA20, close to CD180 and MAST4 genes, both related to fat accumulation. We observed a low positive genetic correlation between IMF-ADG (r = 0.30 ± 0.0686), indicating that it may respond to selection in the same direction. Our findings contributed to clarifying the pleiotropic modulation of the complex traits, indicating new QTLs for bovine genetic improvement.

Three laws of teleonometrics

We define teleonometrics as the theoretical and empirical study of teleonomy. We propose three laws for teleonometrics. The first law describes the hierarchical organization of teleonomic functions across biological levels from genes to individuals. According to this law, the number of goal-directed functions increases from individuals (one goal, maximizing inclusive fitness) to intermediate levels and to genes and alleles (myriad time-, space- and context-dependent goals, depending upon degrees and patterns of pleiotropy). The second law describes the operation of teleonomic functions under trade-offs, coadaptations and negative and positive pleiotropies, which are universal in biological systems. According to this law, the functions of an allele, gene or trait are described and defined by patterns of antagonistic (trading off) and compatible (coadapted) functions. The third law of teleonometrics is that the major transitions in evolution are driven by the origins of novel, emergent goals associated with functional changes and by the breaking and reshaping of trade-offs, especially by mechanisms involving increases in resources or time, and new divisions of labour or function. We illustrate the application of these laws using data from three empirical vignettes, which help to show the usefulness of teleonometric viewpoints for understanding the interfaces between function, trade-offs and dysfunctions manifest as disease.

Pleiotropic constraints promote the evolution of cooperation in cellular groups

The evolution of cooperation in cellular groups is threatened by lineages of cheaters that proliferate at the expense of the group. These cell lineages occur within microbial communities, and multicellular organisms in the form of tumours and cancer. In contrast to an earlier study, here we show how the evolution of pleiotropic genetic architectures—which link the expression of cooperative and private traits—can protect against cheater lineages and allow cooperation to evolve. We develop an age-structured model of cellular groups and show that cooperation breaks down more slowly within groups that tie expression to a private trait than in groups that do not. We then show that this results in group selection for pleiotropy, which strongly promotes cooperation by limiting the emergence of cheater lineages. These results predict that pleiotropy will rapidly evolve, so long as groups persist long enough for cheater lineages to threaten cooperation. Our results hold when pleiotropic links can be undermined by mutations, when pleiotropy is itself costly, and in mixed-genotype groups such as those that occur in microbes. Finally, we consider features of multicellular organisms—a germ line and delayed reproductive maturity—and show that pleiotropy is again predicted to be important for maintaining cooperation. The study of cancer in multicellular organisms provides the best evidence for pleiotropic constraints, where abberant cell proliferation is linked to apoptosis, senescence, and terminal differentiation. Alongside development from a single cell, we propose that the evolution of pleiotropic constraints has been critical for cooperation in many cellular groups.

The evolution of cooperation in cellular groups is threatened by lineages of cheaters that proliferate at the expense of the group. In this study, an age-structured model of cellular groups shows that pleiotropy promotes the evolution of cooperation and may have been important for the origins of multicellularity.

Trade-Offs (and Constraints) in Organismal Biology

AbstractTrade-offs and constraints are inherent to life, and studies of these phenomena play a central role in both organismal and evolutionary biology. Trade-offs can be defined, categorized, and studied in at least six, not mutually exclusive, ways. (1) Allocation constraints are caused by a limited resource (e.g., energy, time, space, essential nutrients), such that increasing allocation to one component necessarily requires a decrease in another (if only two components are involved, this is referred to as the Y-model, e.g., energy devoted to size versus number of offspring). (2) Functional conflicts occur when features that enhance performance of one task decrease performance of another (e.g., relative lengths of in-levers and out-levers, force-velocity trade-offs related to muscle fiber type composition). (3) Shared biochemical pathways, often involving integrator molecules (e.g., hormones, neurotransmitters, transcription factors), can simultaneously affect multiple traits, with some effects being beneficial for one or more components of Darwinian fitness (e.g., survival, age at first reproduction, fecundity) and others detrimental. (4) Antagonistic pleiotropy describes genetic variants that increase one component of fitness (or a lower-level trait) while simultaneously decreasing another. (5) Ecological circumstances (or selective regime) may impose trade-offs, such as when foraging behavior increases energy availability yet also decreases survival. (6) Sexual selection may lead to the elaboration of (usually male) secondary sexual characters that improve mating success but handicap survival and/or impose energetic costs that reduce other fitness components. Empirical studies of trade-offs often search for negative correlations between two traits that are the expected outcomes of the trade-offs, but this will generally be inadequate if more than two traits are involved and especially for complex physiological networks of interacting traits. Moreover, trade-offs often occur only in populations that are experiencing harsh environmental conditions or energetic challenges at the extremes of phenotypic distributions, such as among individuals or species that have exceptional athletic abilities. Trade-offs may be (partially) circumvented through various compensatory mechanisms, depending on the timescale involved, ranging from acute to evolutionary. Going forward, a pluralistic view of trade-offs and constraints, combined with integrative analyses that cross levels of biological organization and traditional boundaries among disciplines, will enhance the study of evolutionary organismal biology.

Pleiotropy or linkage? Their relative contributions to the genetic correlation of quantitative traits and detection by multitrait GWA studies

Genetic correlations between traits may cause correlated responses to selection. Previous models described the conditions under which genetic correlations are expected to be maintained. Selection, mutation, and migration are all proposed to affect genetic correlations, regardless of whether the underlying genetic architecture consists of pleiotropic or tightly linked loci affecting the traits. Here, we investigate the conditions under which pleiotropy and linkage have different effects on the genetic correlations between traits by explicitly modeling multiple genetic architectures to look at the effects of selection strength, degree of correlational selection, mutation rate, mutational variance, recombination rate, and migration rate. We show that at mutation-selection(-migration) balance, mutation rates differentially affect the equilibrium levels of genetic correlation when architectures are composed of pairs of physically linked loci compared to architectures of pleiotropic loci. Even when there is perfect linkage (no recombination within pairs of linked loci), a lower genetic correlation is maintained than with pleiotropy, with a lower mutation rate leading to a larger decrease. These results imply that the detection of causal loci in multitrait association studies will be affected by the type of underlying architectures, whereby pleiotropic variants are more likely to be underlying multiple detected associations. We also confirm that tighter linkage between nonpleiotropic causal loci maintains higher genetic correlations at the traits and leads to a greater proportion of false positives in association analyses.

The relative impact of evolving pleiotropy and mutational correlation on trait divergence

Both pleiotropic connectivity and mutational correlations can restrict the decoupling of traits under divergent selection, but it is unknown which is more important in trait evolution. To address this question, we create a model that permits within-population variation in both pleiotropic connectivity and mutational correlation, and compare their relative importance to trait evolution. Specifically, we developed an individual-based stochastic model where mutations can affect whether a locus affects a trait and the extent of mutational correlations in a population. We find that traits can decouple whether there is evolution in pleiotropic connectivity or mutational correlation, but when both can evolve, then evolution in pleiotropic connectivity is more likely to allow for decoupling to occur. The most common genotype found in this case is characterized by having one locus that maintains connectivity to all traits and another that loses connectivity to the traits under stabilizing selection (subfunctionalization). This genotype is favored because it allows the subfunctionalized locus to accumulate greater effect size alleles, contributing to increasingly divergent trait values in the traits under divergent selection without changing the trait values of the other traits (genetic modularization). These results provide evidence that partial subfunctionalization of pleiotropic loci may be a common mechanism of trait decoupling under regimes of corridor selection.

Oxygen Plasma Substrate and Specific Nanopattern Promote Early Differentiation of HepaRG Progenitors

Fully differentiated HepaRG™ cells are the hepatic cell line of choice for in vitro study in toxicology and drug trials. They are derived from a hepatoblast-like progenitor (HepaRG-P) that differentiates into a coculture of hepatocyte-like and cholangiocyte-like cells. This process that requires 2 weeks of proliferation followed by 2 weeks of differentiation using dimethyl sulfoxide (DMSO) can be time consuming and costly. Identifying a method to accelerate HepaRG-Ps toward a mature lineage would save both time and money. The ability to do this in the absence of DMSO would remove the possibility of confounding toxicology results caused by DMSO induction of CYP pathways. It has been shown that tissue culture substrates play an important role in the development and maturity of a cell line, and this is particularly important for progenitor cells, which retain some form of plasticity. Oxygen plasma treatment is used extensively to modify cell culture substrates. There is also evidence that patterned rather than planar surfaces have a positive effect on proliferation and differentiation. In this study, we compared the effect of standard tissue culture plastic (TCP), oxygen plasma coated (OPC), and nanopatterned substrates (NPS) on early differentiation and function of HepaRG-P cells. Since NPS were OPC we initially compared the effect of TCP and OPC to enable comparison between all three culture surfaces using OPC as control to asses if patterning further enhanced early differentiation and functionality. The results show that HepaRG-P's grown on OPC substrate exhibited earlier differentiation, proliferation, and function compared with TCP. Culturing HepaRG-P's on OPC with the addition of NPS did not confer any additional advantage. In conclusion, OPC surface appeared to enhance hepatic differentiation and functionality and could replace traditional methods of differentiating HepaRG-P cells into fully differentiated and functional HepaRGs earlier than standard methods. Impact statement We show significantly earlier differentiation and function of HepaRG progenitor cells when grown in dimethyl sulfoxide-free medium on oxygen plasma substrates versus standard tissue culture plastic. Further investigation showed that nanopatterning of oxygen plasma substrates did not confer any additional advantage over smooth oxygen plasma, although one pattern (DSQ120) showed comparable early differentiation and function.

Every breath you take: assessing metabolic costs of toxin resistance in garter snakes (Thamnophis)

Trait specialization often comes at the expense of original trait function, potentially causing evolutionary tradeoffs that may render specialist populations vulnerable to extinction. However, many specialized adaptations evolve repeatedly, suggesting selection favors specialization in specific environments. Some garter snake (Thamnophis) populations possess specialized mutations in voltage-gated sodium channels that allow them to consume Pacific newts (Taricha) defended by a highly potent neurotoxin (tetrodotoxin). These mutations, however, also decrease protein and muscle function, suggesting garter snakes may suffer evolutionary tradeoffs. We measured a key physiological process, standard metabolic rate (SMR), to investigate whether specialized adaptations in toxin-resistant garter snakes affect baseline energy expenditure. In snakes, skeletal muscles influence metabolism and power ventilation, so inefficiencies of sodium channels in these muscles might impact whole-animal energy expenditure. Further, because sodium channels are membrane-bound proteins, inefficiencies of channel kinetics and performance might be exacerbated at suboptimal temperatures. We measured SMR in 2 species, Thamnophis atratus and Thamnophis sirtalis, that independently evolved tetrodotoxin resistance through unique mutations, providing replicate experiments with distinct underlying genetics and potential physiological costs. Despite our expectations, neither resistance phenotype nor sodium channel genotype affected metabolism and resistant snakes did not perform worse under suboptimal body temperature. Instead, T. atratus and T. sirtalis show nearly identical rates of mass-adjusted energy expenditure at both temperatures, despite differing eco-morphologies, life histories, and distant phylogenetic positions. These findings suggest SMR may be a conserved feature of Thamnophis, and that any organismal tradeoffs may be compensated to retain whole-animal function.

Enhancer pleiotropy, gene expression, and the architecture of human enhancer-gene interactions

Enhancers are often studied as noncoding regulatory elements that modulate the precise spatiotemporal expression of genes in a highly tissue-specific manner. This paradigm has been challenged by recent evidence of individual enhancers acting in multiple tissues or developmental contexts. However, the frequency of these enhancers with high degrees of “pleiotropy” out of all putative enhancers is not well understood. Consequently, it is unclear how the variation of enhancer pleiotropy corresponds to the variation in expression breadth of target genes. Here, we use multi-tissue chromatin maps from diverse human tissues to investigate the enhancer–gene interaction architecture while accounting for 1) the distribution of enhancer pleiotropy, 2) the variations of regulatory links from enhancers to target genes, and 3) the expression breadth of target genes. We show that most enhancers are tissue-specific and that highly pleiotropy enhancers account for <1% of all putative regulatory sequences in the human genome. Notably, several genomic features are indicative of increasing enhancer pleiotropy, including longer sequence length, greater number of links to genes, increasing abundance and diversity of encoded transcription factor motifs, and stronger evolutionary conservation. Intriguingly, the number of enhancers per gene remains remarkably consistent for all genes (∼14). However, enhancer pleiotropy does not directly translate to the expression breadth of target genes. We further present a series of Gaussian Mixture Models to represent this organization architecture. Consequently, we demonstrate that a modest trend of more pleiotropic enhancers targeting more broadly expressed genes can generate the observed diversity of expression breadths in the human genome.

Stressful development: integrating endoderm development, stress, and longevity

In addition to performing digestion and nutrient absorption, the intestine serves as one of the first barriers to the external environment, crucial for protecting the host from environmental toxins, pathogenic invaders, and other stress inducers. The gene regulatory network (GRN) governing embryonic development of the endoderm and subsequent differentiation and maintenance of the intestine has been well-documented in C. elegans. A key regulatory input that initiates activation of the embryonic GRN for endoderm and mesoderm in this animal is the maternally provided SKN-1 transcription factor, an ortholog of the vertebrate Nrf1 and 2, which, like C. elegans SKN-1, perform conserved regulatory roles in mediating a variety of stress responses across metazoan phylogeny. Other key regulatory factors in early gut development also participate in stress response as well as in innate immunity and aging and longevity. In this review, we discuss the intersection between genetic nodes that mediate endoderm/intestine differentiation and regulation of stress and homeostasis. We also consider how direct signaling from the intestine to the germline, in some cases involving SKN-1, facilitates heritable epigenetic changes, allowing transmission of adaptive stress responses across multiple generations. These connections between regulation of endoderm/intestine development and stress response mechanisms suggest that varying selective pressure exerted on the stress response pathways may influence the architecture of the endoderm GRN, thereby leading to genetic and epigenetic variation in early embryonic GRN regulatory events.

Marker Selection in Multivariate Genomic Prediction Improves Accuracy of Low Heritability Traits

Multivariate analysis using mixed models allows for the exploration of genetic correlations between traits. Additionally, the transition to a genomic based approach is simplified by substituting classic pedigrees with a marker-based relationship matrix. It also enables the investigation of correlated responses to selection, trait integration and modularity in different kinds of populations. This study investigated a strategy for the construction of a marker-based relationship matrix that prioritized markers using Partial Least Squares. The efficiency of this strategy was found to depend on the correlation structure between investigated traits. In terms of accuracy, we found no benefit of this strategy compared with the all-marker-based multivariate model for the primary trait of diameter at breast height (DBH) in a radiata pine (Pinus radiata) population, possibly due to the presence of strong and well-estimated correlation with other highly heritable traits. Conversely, we did see benefit in a shining gum (Eucalyptus nitens) population, where the primary trait had low or only moderate genetic correlation with other low/moderately heritable traits. Marker selection in multivariate analysis can therefore be an efficient strategy to improve prediction accuracy for low heritability traits due to improved precision in poorly estimated low/moderate genetic correlations. Additionally, our study identified the genetic diversity as a factor contributing to the efficiency of marker selection in multivariate approaches due to higher precision of genetic correlation estimates.

Trade-offs, Pleiotropy, and Shared Molecular Pathways: A Unified View of Constraints on Adaptation

The concept of trade-offs permeates our thinking about adaptive evolution because they are exhibited at every level of biological organization, from molecular and cellular processes to organismal and ecological functions. Trade-offs inevitably arise because different traits do not occur in isolation, but instead are imbedded within complex, integrated systems that make up whole organisms. The genetic and mechanistic underpinning of trade-offs can be found in the pleiotropic nodes that occur in the biological pathways shared between traits. Yet, often trade-offs are only understood as statistical correlations, limiting the ability to evaluate the interplay between how selection and constraint interact during adaptive evolution. Here, we first review the classic paradigms in which physiologists and evolutionary biologists have studied trade-offs and highlight the ways in which network and molecular pathway approaches unify these paradigms. We discuss how these approaches allow researchers to evaluate why trade-offs arise and how selection can act to overcome trait correlations and evolutionary constraints. We argue that understanding how the conserved molecular pathways are shared between different traits and functions provides a conceptual framework for evolutionary biologists, physiologists, and molecular biologists to meaningfully work together toward the goal of understanding why correlations and trade-offs occur between traits. We briefly highlight the melanocortin system and the hormonal control of osmoregulation as two case studies where an understanding of shared molecular pathways reveals why trade-offs occur between seemingly unrelated traits. While we recognize that applying such approaches poses challenges and limitations particularly in the context of natural populations, we advocate for the view that focusing on the biological pathways responsible for trade-offs provides a unified conceptual context accessible to a broad range of integrative biologists.

Computational Modeling of Anthocyanin Pathway Evolution: Biases, Hotspots, and Trade-offs

The alteration of metabolic pathways is a common mechanism underlying the evolution of new phenotypes. Flower color is a striking example of the importance of metabolic evolution in a complex phenotype, wherein shifts in the activity of the underlying pathway lead to a wide range of pigments. Although experimental work has identified common classes of mutations responsible for transitions among colors, we lack a unifying model that relates pathway function and activity to the evolution of distinct pigment phenotypes. One challenge in creating such a model is the branching structure of pigment pathways, which may lead to evolutionary trade-offs due to competition for shared substrates. In order to predict the effects of shifts in enzyme function and activity on pigment production, we created a simple kinetic model of a major plant pigmentation pathway: the anthocyanin pathway. This model describes the production of the three classes of blue, purple, and red anthocyanin pigments, and accordingly, includes multiple branches and substrate competition. We first studied the general behavior of this model using a naïve set of parameters. We then stochastically evolved the pathway toward a defined optimum and analyzed the patterns of fixed mutations. This approach allowed us to quantify the probability density of trajectories through pathway state space and identify the types and number of changes. Finally, we examined whether our simulated results qualitatively align with experimental observations, i.e., the predominance of mutations which change color by altering the function of branching genes in the pathway. These analyses provide a theoretical framework that can be used to predict the consequences of new mutations in terms of both pigment phenotypes and pleiotropic effects.

Pleiotropy in developmental regulation by flowering-pathway genes: is it an evolutionary constraint?

Pleiotropy occurs when one gene influences more than one trait, contributing to genetic correlations among traits. Consequently, it is considered a constraint on the evolution of adaptive phenotypes because of potential antagonistic selection on correlated traits, or, alternatively, preservation of functional trait combinations. Such evolutionary constraints may be mitigated by the evolution of different functions of pleiotropic genes in their regulation of different traits. Arabidopsis thaliana flowering-time genes, and the pathways in which they operate, are among the most thoroughly studied regarding molecular functions, phenotypic effects, and adaptive significance. Many of them show strong pleiotropic effects. Here, we review examples of pleiotropy of flowering-time genes and highlight those that also influence seed germination. Some genes appear to operate in the same genetic pathways when regulating both traits, whereas others show diversity of function in their regulation, either interacting with the same genetic partners but in different ways or potentially interacting with different partners. We discuss how functional diversification of pleiotropic genes in the regulation of different traits across the life cycle may mitigate evolutionary constraints of pleiotropy, permitting traits to respond more independently to environmental cues, and how it may even contribute to the evolutionary divergence of gene function across taxa.

Genome wide identification and comparative analysis of glutathione transferases (GST) family genes in Brassica napus

Glutathione transferases (GSTs) are multifunctional enzymes that play important roles in plant development and responses to biotic and abiotic stress. However, a systematic analysis of GST family members in Brassica napus has not yet been reported. In this study, we identified 179 full-length GST genes in B. napus, 44.2% of which are clustered on various chromosomes. In addition, we identified 141 duplicated GST gene pairs in B. napus. Molecular evolutionary analysis showed that speciation and whole-genome triplication played important roles in the divergence of the B. napus GST duplicated genes. Transcriptome analysis of 21 tissues at different developmental stages showed that 47.6% of duplicated GST gene pairs have divergent expression patterns, perhaps due to structural divergence. We constructed a GST gene coexpression network with genes encoding various transcription factors (NAC, MYB, WRKY and bZIP) and identified six modules, including genes expressed during late seed development (after 40 days; BnGSTU19, BnGSTU20 and BnGSTZ1) and in the seed coat (BnGSTF6 and BnGSTF12), stamen and anther (BnGSTF8), root and stem (BnGSTU21), leaves and funiculus, as well as during the late stage of pericarp development (after 40 days; BnGSTU12 and BnGSTF2) and in the radicle during seed germination (BnGSTF14, BnGSTU1, BnGSTU28, and BnGSTZ1). These findings lay the foundation for elucidating the roles of GSTs in B. napus.

Symbiont evolution during the free-living phase can improve host colonization

For micro-organisms cycling between free-living and host-associated stages, where reproduction occurs in both of these lifestyles, an interesting inquiry is whether evolution during the free-living stage can be positively pleiotropic to microbial fitness in a host environment. To address this topic, the squid host Euprymna tasmanica and the marine bioluminescent bacterium Vibrio fischeri were utilized. Microbial ecological diversification in static liquid microcosms was used to simulate symbiont evolution during the free-living stage. Thirteen genetically distinct V. fischeri strains from a broad diversity of ecological sources (e.g. squid light organs, fish light organs and seawater) were examined to see if the results were reproducible in many different genetic settings. Genetic backgrounds that are closely related can be predisposed to considerable differences in how they respond to similar selection pressures. For all strains examined, new mutations with striking and facilitating effects on host colonization arose quickly during microbial evolution in the free-living stage, regardless of the ecological context under consideration for a strain’s genetic background. Microbial evolution outside a host environment promoted host range expansion, improved host colonization for a micro-organism, and diminished the negative correlation between biofilm formation and motility.

Engineered Functional Redundancy Relaxes Selective Constraints upon Endogenous Genes in Viral RNA Genomes

Functional redundancy, understood as the functional overlap of different genes, is a double-edge sword. At the one side, it is thought to serve as a robustness mechanism that buffers the deleterious effect of mutations hitting one of the redundant copies, thus resulting in pseudogenization. At the other side, it is considered as a source of genetic and functional innovation. In any case, genetically redundant genes are expected to show an acceleration in the rate of molecular evolution. Here, we tackle the role of functional redundancy in viral RNA genomes. To this end, we have evaluated the rates of compensatory evolution for deleterious mutations affecting an essential function, the suppression of RNA silencing plant defense, of tobacco etch potyvirus (TEV). TEV genotypes containing deleterious mutations in presence/absence of engineered functional redundancy were evolved and the pattern of fitness and pathogenicity recovery evaluated. Genetically redundant genotypes suffered less from the effect of deleterious mutations and showed relatively minor changes in fitness and pathogenicity. By contrast, nongenetically redundant genotypes had very low fitness and pathogenicity at the beginning of the evolution experiment that were fully recovered by the end. At the molecular level, the outcome depended on the combination of the actual mutations being compensated and the presence/absence of functional redundancy. Reversions to wild-type alleles were the norm in the nonredundant genotypes while redundant ones either did not fix any mutation at all or showed a higher nonsynonymous mutational load.

Homeotic transformations reflect departure from the mammalian 'rule of seven' cervical vertebrae in sloths: inferences on the Hox code and morphological modularity of the mammalian neck

Background

Sloths are one of only two exceptions to the mammalian ‘rule of seven’ vertebrae in the neck. As a striking case of breaking the evolutionary constraint, the explanation for the exceptional number of cervical vertebrae in sloths is still under debate. Two diverging hypotheses, both ultimately linked to the low metabolic rate of sloths, have been proposed: hypothesis 1 involves morphological transformation of vertebrae due to changes in the Hox gene expression pattern and hypothesis 2 assumes that the Hox gene expression pattern is not altered and the identity of the vertebrae is not changed. Direct evidence supporting either hypothesis would involve knowledge of the vertebral Hox code in sloths, but the realization of such studies is extremely limited. Here, on the basis of the previously established correlation between anterior Hox gene expression and the quantifiable vertebral shape, we present the morphological regionalization of the neck in three different species of sloths with aberrant cervical count providing indirect insight into the vertebral Hox code.

Results

Shape differences within the cervical vertebral column suggest a mouse-like Hox code in the neck of sloths. We infer an anterior shift of HoxC-6 expression in association with the first thoracic vertebra in short-necked sloths with decreased cervical count, and a posterior shift of HoxC-5 and HoxC-6 expression in long-necked sloths with increased cervical count.

Conclusion

Although only future developmental analyses in non-model organisms, such as sloths, will yield direct evidence for the evolutionary mechanism responsible for the aberrant number of cervical vertebrae, our observations lend support to hypothesis 1 indicating that the number of modules is retained but their boundaries are displaced. Our approach based on quantified morphological differences also provides a reliable basis for further research including fossil taxa such as extinct ‘ground sloths’ in order to trace the pattern and the underlying genetic mechanisms in the evolution of the vertebral column in mammals.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1202-5) contains supplementary material, which is available to authorized users.

Gene Regulatory Enhancers with Evolutionarily Conserved Activity Are More Pleiotropic than Those with Species-Specific Activity

Studies of regulatory activity and gene expression have revealed an intriguing dichotomy: There is substantial turnover in the regulatory activity of orthologous sequences between species; however, the expression level of orthologous genes is largely conserved. Understanding how distal regulatory elements, for example, enhancers, evolve and function is critical, as alterations in gene expression levels can drive the development of both complex disease and functional divergence between species. In this study, we investigated determinants of the conservation of regulatory enhancer activity for orthologous sequences across mammalian evolution. Using liver enhancers identified from genome-wide histone modification profiles in ten diverse mammalian species, we compared orthologous sequences that exhibited regulatory activity in all species (conserved-activity enhancers) to shared sequences active only in a single species (species-specific-activity enhancers). Conserved-activity enhancers have greater regulatory potential than species-specific-activity enhancers, as quantified by both the density and diversity of transcription factor binding motifs. Consistent with their greater regulatory potential, conserved-activity enhancers have greater regulatory activity in humans than species-specific-activity enhancers: They are active across more cellular contexts, and they regulate more genes than species-specific-activity enhancers. Furthermore, the genes regulated by conserved-activity enhancers are expressed in more tissues and are less tolerant of loss-of-function mutations than those targeted by species-specific-activity enhancers. These consistent results across various stages of gene regulation demonstrate that conserved-activity enhancers are more pleiotropic than their species-specific-activity counterparts. This suggests that pleiotropy is associated with the conservation of regulatory across mammalian evolution.

Parallel evolution of the POQR prolyl oligo peptidase gene conferring plant quantitative disease resistance

Plant pathogens with a broad host range are able to infect plant lineages that diverged over 100 million years ago. They exert similar and recurring constraints on the evolution of unrelated plant populations. Plants generally respond with quantitative disease resistance (QDR), a form of immunity relying on complex genetic determinants. In most cases, the molecular determinants of QDR and how they evolve is unknown. Here we identify in Arabidopsis thaliana a gene mediating QDR against Sclerotinia sclerotiorum, agent of the white mold disease, and provide evidence of its convergent evolution in multiple plant species. Using genome wide association mapping in A. thaliana, we associated the gene encoding the POQR prolyl-oligopeptidase with QDR against S. sclerotiorum. Loss of this gene compromised QDR against S. sclerotiorum but not against a bacterial pathogen. Natural diversity analysis associated POQR sequence with QDR. Remarkably, the same amino acid changes occurred after independent duplications of POQR in ancestors of multiple plant species, including A. thaliana and tomato. Genome-scale expression analyses revealed that parallel divergence in gene expression upon S. sclerotiorum infection is a frequent pattern in genes, such as POQR, that duplicated both in A. thaliana and tomato. Our study identifies a previously uncharacterized gene mediating QDR against S. sclerotiorum. It shows that some QDR determinants are conserved in distantly related plants and have emerged through the repeated use of similar genetic polymorphisms at different evolutionary time scales.

Protein Moonlighting Revealed by Noncatalytic Phenotypes of Yeast Enzymes

A single gene can partake in several biological processes, and therefore gene deletions can lead to different-sometimes unexpected-phenotypes. However, it is not always clear whether such pleiotropy reflects the loss of a unique molecular activity involved in different processes or the loss of a multifunctional protein. Here, using Saccharomyces cerevisiae metabolism as a model, we systematically test the null hypothesis that enzyme phenotypes depend on a single annotated molecular function, namely their catalysis. We screened a set of carefully selected genes by quantifying the contribution of catalysis to gene deletion phenotypes under different environmental conditions. While most phenotypes were explained by loss of catalysis, slow growth was readily rescued by a catalytically inactive protein in about one-third of the enzymes tested. Such noncatalytic phenotypes were frequent in the Alt1 and Bat2 transaminases and in the isoleucine/valine biosynthetic enzymes Ilv1 and Ilv2, suggesting novel "moonlighting" activities in these proteins. Furthermore, differential genetic interaction profiles of gene deletion and catalytic mutants indicated that ILV1 is functionally associated with regulatory processes, specifically to chromatin modification. Our systematic study shows that gene loss phenotypes and their genetic interactions are frequently not driven by the loss of an annotated catalytic function, underscoring the moonlighting nature of cellular metabolism.

What affects the predictability of evolutionary constraints using a G-matrix? The relative effects of modular pleiotropy and mutational correlation

Phenotypic traits do not always respond to selection independently from each other and often show correlated responses to selection. The structure of a genotype-phenotype map (GP map) determines trait covariation, which involves variation in the degree and strength of the pleiotropic effects of the underlying genes. It is still unclear, and debated, how much of that structure can be deduced from variational properties of quantitative traits that are inferred from their genetic (co) variance matrix (G-matrix). Here we aim to clarify how the extent of pleiotropy and the correlation among the pleiotropic effects of mutations differentially affect the structure of a G-matrix and our ability to detect genetic constraints from its eigen decomposition. We show that the eigenvectors of a G-matrix can be predictive of evolutionary constraints when they map to underlying pleiotropic modules with correlated mutational effects. Without mutational correlation, evolutionary constraints caused by the fitness costs associated with increased pleiotropy are harder to infer from evolutionary metrics based on a G-matrix's geometric properties because uncorrelated pleiotropic effects do not affect traits' genetic correlations. Correlational selection induces much weaker modular partitioning of traits' genetic correlations in absence then in presence of underlying modular pleiotropy.

Epistasis and Pleiotropy Affect the Modularity of the Genotype-Phenotype Map of Cross-Resistance in HIV-1

The genotype–phenotype (GP) map is a central concept in evolutionary biology as it describes the mapping of molecular genetic variation onto phenotypic trait variation. Our understanding of that mapping remains partial, especially when trying to link functional clustering of pleiotropic gene effects with patterns of phenotypic trait co-variation. Only on rare occasions have studies been able to fully explore that link and tend to show poor correspondence between modular structures within the GP map and among phenotypes. By dissecting the structure of the GP map of the replicative capacity of HIV-1 in 15 drug environments, we provide a detailed view of that mapping from mutational pleiotropic variation to phenotypic co-variation, including epistatic effects of a set of amino-acid substitutions in the reverse transcriptase and protease genes. We show that epistasis increases the pleiotropic degree of single mutations and provides modularity to the GP map of drug resistance in HIV-1. Moreover, modules of epistatic pleiotropic effects within the GP map match the phenotypic modules of correlated replicative capacity among drug classes. Epistasis thus increases the evolvability of cross-resistance in HIV by providing more drug- and class-specific pleiotropic profiles to the main effects of the mutations. We discuss the implications for the evolution of cross-resistance in HIV.

Evolution of a transcriptional regulator from a transmembrane nucleoporin

Franks et al. identify a widely expressed variant of the transmembrane nucleoporin Pom121 (named sPom121, for “soluble Pom121”) that arose by genomic rearrangement before the divergence of hominoids. Instead of localizing to the NPC, sPom121 colocalizes and interacts with nucleoplasmic Nup98, a previously identified transcriptional regulator, at gene promoters to control transcription of its target genes in human cells.

Nuclear pore complexes (NPCs) emerged as nuclear transport channels in eukaryotic cells ∼1.5 billion years ago. While the primary role of NPCs is to regulate nucleo–cytoplasmic transport, recent research suggests that certain NPC proteins have additionally acquired the role of affecting gene expression at the nuclear periphery and in the nucleoplasm in metazoans. Here we identify a widely expressed variant of the transmembrane nucleoporin (Nup) Pom121 (named sPom121, for “soluble Pom121”) that arose by genomic rearrangement before the divergence of hominoids. sPom121 lacks the nuclear membrane-anchoring domain and thus does not localize to the NPC. Instead, sPom121 colocalizes and interacts with nucleoplasmic Nup98, a previously identified transcriptional regulator, at gene promoters to control transcription of its target genes in human cells. Interestingly, sPom121 transcripts appear independently in several mammalian species, suggesting convergent innovation of Nup-mediated transcription regulation during mammalian evolution. Our findings implicate alternate transcription initiation as a mechanism to increase the functional diversity of NPC components.

Convergent evolution of mechanically optimal locomotion in aquatic invertebrates and vertebrates

Examples of animals evolving similar traits despite the absence of that trait in the last common ancestor, such as the wing and camera-type lens eye in vertebrates and invertebrates, are called cases of convergent evolution. Instances of convergent evolution of locomotory patterns that quantitatively agree with the mechanically optimal solution are very rare. Here, we show that, with respect to a very diverse group of aquatic animals, a mechanically optimal method of swimming with elongated fins has evolved independently at least eight times in both vertebrate and invertebrate swimmers across three different phyla. Specifically, if we take the length of an undulation along an animal's fin during swimming and divide it by the mean amplitude of undulations along the fin length, the result is consistently around twenty. We call this value the optimal specific wavelength (OSW). We show that the OSW maximizes the force generated by the body, which also maximizes swimming speed. We hypothesize a mechanical basis for this optimality and suggest reasons for its repeated emergence through evolution.

Integrated phenotypes: understanding trait covariation in plants and animals

Integration and modularity refer to the patterns and processes of trait interaction and independence. Both terms have complex histories with respect to both conceptualization and quantification, resulting in a plethora of integration indices in use. We review briefly the divergent definitions, uses and measures of integration and modularity and make conceptual links to allometry. We also discuss how integration and modularity might evolve. Although integration is generally thought to be generated and maintained by correlational selection, theoretical considerations suggest the relationship is not straightforward. We caution here against uncontrolled comparisons of indices across studies. In the absence of controls for trait number, dimensionality, homology, development and function, it is difficult, or even impossible, to compare integration indices across organisms or traits. We suggest that care be invested in relating measurement to underlying theory or hypotheses, and that summative, theory-free descriptors of integration generally be avoided. The papers that follow in this Theme Issue illustrate the diversity of approaches to studying integration and modularity, highlighting strengths and pitfalls that await researchers investigating integration in plants and animals.

The evolution of phenotypic correlations and "developmental memory"

Development introduces structured correlations among traits that may constrain or bias the distribution of phenotypes produced. Moreover, when suitable heritable variation exists, natural selection may alter such constraints and correlations, affecting the phenotypic variation available to subsequent selection. However, exactly how the distribution of phenotypes produced by complex developmental systems can be shaped by past selective environments is poorly understood. Here we investigate the evolution of a network of recurrent nonlinear ontogenetic interactions, such as a gene regulation network, in various selective scenarios. We find that evolved networks of this type can exhibit several phenomena that are familiar in cognitive learning systems. These include formation of a distributed associative memory that can "store" and "recall" multiple phenotypes that have been selected in the past, recreate complete adult phenotypic patterns accurately from partial or corrupted embryonic phenotypes, and "generalize" (by exploiting evolved developmental modules) to produce new combinations of phenotypic features. We show that these surprising behaviors follow from an equivalence between the action of natural selection on phenotypic correlations and associative learning, well-understood in the context of neural networks. This helps to explain how development facilitates the evolution of high-fitness phenotypes and how this ability changes over evolutionary time.

Flower color as a model system for studies of plant evo-devo

Even though pigmentation traits have had substantial impacts on the field of animal evolutionary developmental biology, they have played only relatively minor roles in plant evo-devo. This is surprising given the often direct connection between flower color and fitness variation mediated through the effects of pollinators. At the same time, ecological and evolutionary genetic studies have utilized the molecular resources available for the anthocyanin pathway to generate several examples of the molecular basis of putatively adaptive transitions in flower color. Despite this opportunity to synthesize experimental approaches in ecology, evolution, and developmental biology, the investigation of many fundamental questions in evo-devo using this powerful model is only at its earliest stages. For example, a long-standing question is whether predictable genetic changes accompany the repeated evolution of a trait. Due to the conserved nature of the biochemical and regulatory control of anthocyanin biosynthesis, it has become possible to determine whether, and under what circumstances, different types of mutations responsible for flower color variation are preferentially targeted by natural selection. In addition, because plants use anthocyanin and related compounds in vegetative tissue for other important physiological functions, the identification of naturally occurring transitions from unpigmented to pigmented flowers provides the opportunity to examine the mechanisms by which regulatory networks are co-opted into new developmental domains. Here, we review what is known about the ecological and molecular basis of anthocyanic flower color transitions in natural systems, focusing on the evolutionary and developmental features involved. In doing so, we provide suggestions for future work on this trait and suggest that there is still much to be learned from the evolutionary development of flower color transitions in nature.

Sex-differential selection and the evolution of X inactivation strategies

X inactivation—the transcriptional silencing of one X chromosome copy per female somatic cell—is universal among therian mammals, yet the choice of which X to silence exhibits considerable variation among species. X inactivation strategies can range from strict paternally inherited X inactivation (PXI), which renders females haploid for all maternally inherited alleles, to unbiased random X inactivation (RXI), which equalizes expression of maternally and paternally inherited alleles in each female tissue. However, the underlying evolutionary processes that might account for this observed diversity of X inactivation strategies remain unclear. We present a theoretical population genetic analysis of X inactivation evolution and specifically consider how conditions of dominance, linkage, recombination, and sex-differential selection each influence evolutionary trajectories of X inactivation. The results indicate that a single, critical interaction between allelic dominance and sex-differential selection can select for a broad and continuous range of X inactivation strategies, including unequal rates of inactivation between maternally and paternally inherited X chromosomes. RXI is favored over complete PXI as long as alleles deleterious to female fitness are sufficiently recessive, and the criteria for RXI evolution is considerably more restrictive when fitness variation is sexually antagonistic (i.e., alleles deleterious to females are beneficial to males) relative to variation that is deleterious to both sexes. Evolutionary transitions from PXI to RXI also generally increase mean relative female fitness at the expense of decreased male fitness. These results provide a theoretical framework for predicting and interpreting the evolution of chromosome-wide expression of X-linked genes and lead to several useful predictions that could motivate future studies of allele-specific gene expression variation.

With the exception of its most primitive members, mammal species practice X inactivation, where one copy of each X chromosome pair is silenced in each cell of the female body. The particular copy of the X that is silenced nevertheless shows considerable variability among species, and the evolutionary causes for this variability remain unclear. Here, we show that X inactivation strategies are likely to evolve in response to the sex-differential fitness properties of X-linked genetic variation. Genetic variation with similar effects on male and female fitness will generally favor the evolution of random X inactivation, potentially including preferential inactivation of the maternally inherited X chromosome. Variation with opposing fitness effects in each sex (“sexually antagonistic” variation, which includes mutations that both decrease female fitness and enhance male fitness) selects for preferential or complete inactivation of the paternally inherited X. Paternally biased X inactivation patterns appear to be common in nature, which suggests that sexually antagonistic genetic variation might be an important factor underlying the evolution of X inactivation. The theory provides a conceptual framework for understanding the evolution of X inactivation strategies and generates several novel predictions that may soon be tested with modern genome sequencing technologies.