A model of developmental evolution: selection, pleiotropy and compensation

Temporal genomics in Hawaiian crickets reveals compensatory intragenomic coadaptation during adaptive evolution

Theory predicts that compensatory genetic changes reduce negative indirect effects of selected variants during adaptive evolution, but evidence is scarce. Here, we test this in a wild population of Hawaiian crickets using temporal genomics and a high-quality chromosome-level cricket genome. In this population, a mutation, flatwing, silences males and rapidly spread due to an acoustically-orienting parasitoid. Our sampling spanned a social transition during which flatwing fixed and the population went silent. We find long-range linkage disequilibrium around the putative flatwing locus was maintained over time, and hitchhiking genes had functions related to negative flatwing-associated effects. We develop a combinatorial enrichment approach using transcriptome data to test for compensatory, intragenomic coevolution. Temporal changes in genomic selection were distributed genome-wide and functionally associated with the population's transition to silence, particularly behavioural responses to silent environments. Our results demonstrate how 'adaptation begets adaptation'; changes to the sociogenetic environment accompanying rapid trait evolution can generate selection provoking further, compensatory adaptation.

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.

The cliff edge model of the evolution of schizophrenia: Mathematical, epidemiological, and genetic evidence

How has schizophrenia, a condition that significantly reduces an individual's evolutionary fitness, remained common across generations and cultures? Numerous theories about the evolution of schizophrenia have been proposed, most of which are not consistent with modern epidemiological and genetic evidence. Here, we briefly review this evidence and explore the cliff edge model of schizophrenia. It suggests that schizophrenia is the extreme manifestation of a polygenic trait or a combination of traits that, within a normal range of variation, confer cognitive, linguistic, and/or social advantages. Only beyond a certain threshold, these traits precipitate the onset of schizophrenia and reduce fitness. We provide the first mathematical model of this qualitative concept and show that it requires only very weak positive selection of the underlying trait(s) to explain today's schizophrenia prevalence. This prediction, along with expectations about the effect size of schizophrenia risk alleles, are surprisingly well matched by empirical evidence. The cliff edge model predicts a dynamic change of selection of risk alleles, which explains the contradictory findings of evolutionary genetic studies.

Brief Communication: The Predictable Network Topology of Evolutionary Genomic Constraint

Large-scale comparative genomics studies offer valuable resources for understanding both functional and evolutionary rate constraints. It is suggested that constraint aligns with the topology of genomic networks, increasing toward the center, with intermediate nodes combining relaxed constraint with higher contributions to the phenotype due to pleiotropy. However, this pattern has yet to be demonstrated in vertebrates. This study shows that constraint intensifies toward the network's center in placental mammals. Genes with rate changes associated with emergence of hibernation cluster mostly toward intermediate positions, with higher constraint in faster-evolving genes, which is indicative of a "sweet spot" for adaptation. If this trend holds universally, network node metrics could predict high-constraint regions even in clades lacking empirical constraint data.

A mathematical framework for evo-devo dynamics

Natural selection acts on phenotypes constructed over development, which raises the question of how development affects evolution. Classic evolutionary theory indicates that development affects evolution by modulating the genetic covariation upon which selection acts, thus affecting genetic constraints. However, whether genetic constraints are relative, thus diverting adaptation from the direction of steepest fitness ascent, or absolute, thus blocking adaptation in certain directions, remains uncertain. This limits understanding of long-term evolution of developmentally constructed phenotypes. Here we formulate a general, tractable mathematical framework that integrates age progression, explicit development (i.e., the construction of the phenotype across life subject to developmental constraints), and evolutionary dynamics, thus describing the evolutionary and developmental (evo-devo) dynamics. The framework yields simple equations that can be arranged in a layered structure that we call the evo-devo process, whereby five core elementary components generate all equations including those mechanistically describing genetic covariation and the evo-devo dynamics. The framework recovers evolutionary dynamic equations in gradient form and describes the evolution of genetic covariation from the evolution of genotype, phenotype, environment, and mutational covariation. This shows that genotypic and phenotypic evolution must be followed simultaneously to yield a dynamically sufficient description of long-term phenotypic evolution in gradient form, such that evolution described as the climbing of a fitness landscape occurs in "geno-phenotype" space. Genetic constraints in geno-phenotype space are necessarily absolute because the phenotype is related to the genotype by development. Thus, the long-term evolutionary dynamics of developed phenotypes is strongly non-standard: (1) evolutionary equilibria are either absent or infinite in number and depend on genetic covariation and hence on development; (2) developmental constraints determine the admissible evolutionary path and hence which evolutionary equilibria are admissible; and (3) evolutionary outcomes occur at admissible evolutionary equilibria, which do not generally occur at fitness landscape peaks in geno-phenotype space, but at peaks in the admissible evolutionary path where "total genotypic selection" vanishes if exogenous plastic response vanishes and mutational variation exists in all directions of genotype space. Hence, selection and development jointly define the evolutionary outcomes if absolute mutational constraints and exogenous plastic response are absent, rather than the outcomes being defined only by selection. Moreover, our framework provides formulas for the sensitivities of a recurrence and an alternative method to dynamic optimization (i.e., dynamic programming or optimal control) to identify evolutionary outcomes in models with developmentally dynamic traits. These results show that development has major evolutionary effects.

Rapid and biased evolution of canalization during adaptive divergence revealed by dominance in gene expression variability during Arctic charr early development

Adaptive evolution may be influenced by canalization, the buffering of developmental processes from environmental and genetic perturbations, but how this occurs is poorly understood. Here, we explore how gene expression variability evolves in diverging and hybridizing populations, by focusing on the Arctic charr (Salvelinus alpinus) of Thingvallavatn, a classic case of divergence between feeding habitats. We report distinct profiles of gene expression variance for both coding RNAs and microRNAs between the offspring of two contrasting morphs (benthic/limnetic) and their hybrids reared in common conditions and sampled at two key points of cranial development. Gene expression variance in the hybrids is substantially affected by maternal effects, and many genes show biased expression variance toward the limnetic morph. This suggests that canalization, as inferred by gene expression variance, can rapidly diverge in sympatry through multiple gene pathways, which are associated with dominance patterns possibly biasing evolutionary trajectories and mitigating the effects of hybridization on adaptive evolution.

Phenotypic heterogeneity in human genetic diseases: ultrasensitivity-mediated threshold effects as a unifying molecular mechanism

Phenotypic heterogeneity is very common in genetic systems and in human diseases and has important consequences for disease diagnosis and treatment. In addition to the many genetic and non-genetic (e.g., epigenetic, environmental) factors reported to account for part of the heterogeneity, we stress the importance of stochastic fluctuation and regulatory network topology in contributing to phenotypic heterogeneity. We argue that a threshold effect is a unifying principle to explain the phenomenon; that ultrasensitivity is the molecular mechanism for this threshold effect; and discuss the three conditions for phenotypic heterogeneity to occur. We suggest that threshold effects occur not only at the cellular level, but also at the organ level. We stress the importance of context-dependence and its relationship to pleiotropy and edgetic mutations. Based on this model, we provide practical strategies to study human genetic diseases. By understanding the network mechanism for ultrasensitivity and identifying the critical factor, we may manipulate the weak spot to gently nudge the system from an ultrasensitive state to a stable non-disease state. Our analysis provides a new insight into the prevention and treatment of genetic diseases.

Evolutionary Change in Gut Specification in Caenorhabditis Centers on the GATA Factor ELT-3 in an Example of Developmental System Drift

Cells in a developing animal embryo become specified by the activation of cell-type-specific gene regulatory networks. The network that specifies the gut in the nematode Caenorhabditis elegans has been the subject of study for more than two decades. In this network, the maternal factors SKN-1/Nrf and POP-1/TCF activate a zygotic GATA factor cascade consisting of the regulators MED-1,2 → END-1,3 → ELT-2,7, leading to the specification of the gut in early embryos. Paradoxically, the MED, END, and ELT-7 regulators are present only in species closely related to C. elegans, raising the question of how the gut can be specified without them. Recent work found that ELT-3, a GATA factor without an endodermal role in C. elegans, acts in a simpler ELT-3 → ELT-2 network to specify gut in more distant species. The simpler ELT-3 → ELT-2 network may thus represent an ancestral pathway. In this review, we describe the elucidation of the gut specification network in C. elegans and related species and propose a model by which the more complex network might have formed. Because the evolution of this network occurred without a change in phenotype, it is an example of the phenomenon of Developmental System Drift.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

Genome-Wide Association Mapping for Yield and Yield-Related Traits in Rice (Oryza Sativa L.) Using SNPs Markers

Rice (Oryza sativa L.) is a staple food for more than 50% of the world's population. Rice cultivar improvement is critical in order to feed the world's growing population. Improving yield is one of the main aims of rice breeders. However, yield is a complex quantitative trait controlled by many genes. The presence of genetic diversity is the key factor to improve the yield hence, the presence of diversity in any germplasm is important for yield improvement. In the current study, the rice germplasm was collected from Pakistan and the United States of America and a panel of 100 diverse genotypes was utilized to identify important yield and yield-related traits. For this, a genome-wide association study (GWAS) was performed to identify the genetic loci related to yield. The GWAS on the diverse germplasm will lead to the identification of new genes which can be utilized in the breeding program for improvement of yield. For this reason, firstly, the germplasm was phenotypically evaluated in two growing seasons for yield and yield-related traits. The analysis of variance results showed significant differences among traits which showed the presence of diversity in the current germplasm. Secondly, the germplasm was also genotypically evaluated using 10K SNP. Genetic structure analysis showed the presence of four groups which showed that enough genetic diversity was present in the rice germplasm to be used for association mapping analysis. The results of GWAS identified 201 significant marker trait associations (MTAs. 16 MTAs were identified for plant height, 49 for days to flowering, three for days to maturity, four for tillers per plant, four for panicle length, eight for grains per panicle, 20 unfilled grains per panicle, 81 for seed setting %, four for thousand-grain weight, five for yield per plot and seven for yield per hectare. Apart from this, some pleiotropic loci were also identified. The results showed that panicle length (PL) and thousand-grain weight (TGW) were controlled by a pleiotropic locus OsGRb23906 on chromosome 1 at 10,116,371 cM. The loci OsGRb25803 and OsGRb15974 on chromosomes 4 and 8 at the position of 14,321,111 cM and 6,205,816 cM respectively, showed pleiotropic effects for seed setting % (SS) and unfilled grain per panicle (UG/P). A locus OsGRb09180 on chromosome 4 at 19,850,601 cM was significantly linked with SS and yield/ha. Furthermore, gene annotation was performed, and results indicated that the 190 candidate genes or QTLs that closely linked with studied traits. These candidate genes and novel significant markers could be useful in marker-assisted gene selection and QTL pyramiding to improve rice yield and the selection of potential parents, recombinants and MTAs which could be used in rice breeding programs to develop high-yielding rice varieties for sustainable food security.

Genetic characterization of Macaca arctoides: A highlight of key genes and pathways

When compared to the approximately 22 other macaque species, Macaca arctoides has many unique phenotypes. These traits fall into various phenotypic categories, including genitalia, coloration, mating, and olfactory traits. Here we used a previously identified whole genome set of 690 outlier genes to look for possible genetic explanations of these unique traits. Of these, 279 genes were annotated miRNAs, which are non-coding. Patterns within the remaining outliers in coding genes were investigated using GO (n = 370) and String (n = 383) analysis, which showed many interconnected immune-related genes. Further, we compared the outliers to candidate pathways associated with M. arcotides' unique phenotypes, revealing 10/690 outlier genes that overlapped these four pathways: hedgehog signaling, WNT signaling, olfactory, and melanogenesis. Of these, genes in all pathways except olfactory had higher F_ST values than the rest of the genes in the genome based on permutation tests. Overall, our results point to many genes each having a small impact on phenotype, working in tandem to cause large systemic changes. Additionally, these results may indicate pleiotropy. This seems to be especially true with the development and coloration of M. arctoides. Our results highlight that development, melanogenesis, immune function, and miRNAs may be heavily involved in M. arctoides' evolutionary history.

Evidence for a single, ancient origin of a genus-wide alternative life history strategy

Understanding the evolutionary origins and factors maintaining alternative life history strategies (ALHS) within species is a major goal of evolutionary research. While alternative alleles causing discrete ALHS are expected to purge or fix over time, one-third of the ~90 species of Colias butterflies are polymorphic for a female-limited ALHS called Alba. Whether Alba arose once, evolved in parallel, or has been exchanged among taxa is currently unknown. Using comparative genome-wide association study (GWAS) and population genomic analyses, we placed the genetic basis of Alba in time-calibrated phylogenomic framework, revealing that Alba evolved once near the base of the genus and has been subsequently maintained via introgression and balancing selection. CRISPR-Cas9 mutagenesis was then used to verify a putative cis-regulatory region of Alba, which we identified using phylogenetic foot printing. We hypothesize that this cis-regulatory region acts as a modular enhancer for the induction of the Alba ALHS, which has likely facilitated its long evolutionary persistence.

Multiple solutions at the genomic level in response to selective breeding for high locomotor activity

Replicate lines under uniform selection often evolve in different ways. Previously, analyses using whole-genome sequence data for individual mice (Mus musculus) from 4 replicate High Runner lines and 4 nonselected control lines demonstrated genomic regions that have responded consistently to selection for voluntary wheel-running behavior. Here, we ask whether the High Runner lines have evolved differently from each other, even though they reached selection limits at similar levels. We focus on 1 High Runner line (HR3) that became fixed for a mutation at a gene of major effect (Myh4Minimsc) that, in the homozygous condition, causes a 50% reduction in hindlimb muscle mass and many pleiotropic effects. We excluded HR3 from SNP analyses and identified 19 regions not consistently identified in analyses with all 4 lines. Repeating analyses while dropping each of the other High Runner lines identified 12, 8, and 6 such regions. (Of these 45 regions, 37 were unique.) These results suggest that each High Runner line indeed responded to selection somewhat uniquely, but also that HR3 is the most distinct. We then applied 2 additional analytical approaches when dropping HR3 only (based on haplotypes and nonstatistical tests involving fixation patterns). All 3 approaches identified 7 new regions (as compared with analyses using all 4 High Runner lines) that include genes associated with activity levels, dopamine signaling, hippocampus morphology, heart size, and body size, all of which differ between High Runner and control lines. Our results illustrate how multiple solutions and "private" alleles can obscure general signatures of selection involving "public" alleles.

Allopatric divergence of cooperators confers cheating resistance and limits effects of a defector mutation

BACKGROUND

Social defectors may meet diverse cooperators. Genotype-by-genotype interactions may constrain the ranges of cooperators upon which particular defectors can cheat, limiting cheater spread. Upon starvation, the soil bacterium Myxococcus xanthus cooperatively develops into spore-bearing fruiting bodies, using a complex regulatory network and several intercellular signals. Some strains (cheaters) are unable to sporulate effectively in pure culture due to mutations that reduce signal production but can exploit and outcompete cooperators within mixed groups.

RESULTS

In this study, interactions between a cheater disrupted at the signaling gene csgA and allopatrically diversified cooperators reveal a very small cheating range. Expectedly, the cheater failed to cheat on all natural-isolate cooperators owing to non-cheater-specific antagonisms. Surprisingly, some lab-evolved cooperators had already exited the csgA mutant's cheating range after accumulating fewer than 20 mutations and without experiencing cheating during evolution. Cooperators might also diversify in the potential for a mutation to reduce expression of a cooperative trait or generate a cheating phenotype. A new csgA mutation constructed in several highly diverged cooperators generated diverse sporulation phenotypes, ranging from a complete defect to no defect, indicating that genetic backgrounds can limit the set of genomes in which a mutation creates a defector.

CONCLUSIONS

Our results demonstrate that natural populations may feature geographic mosaics of cooperators that have diversified in their susceptibility to particular cheaters, limiting defectors' cheating ranges and preventing them from spreading. This diversification may also lead to variation in the phenotypes generated by any given cooperation-gene mutation, further decreasing the chance of a cheater emerging which threatens the persistence of cooperation in the system.

Unusual occurrence of domestication syndrome amongst African mole-rats: Is the naked mole-rat a domestic animal?

The Naked mole-rat (NMR) is becoming a prominent model organism due to its peculiar traits, such as eusociality, extreme longevity, cancer resistance, and reduced pain sensitivity. It belongs to the African mole-rats (AMR), a family of subterranean rodents that includes solitary, cooperative breeding and eusocial species. We identified and quantified the domestication syndrome (DS) across AMR, a set of morphological and behavioural traits significantly more common and pronounced amongst domesticated animals than in their wild counterparts. Surprisingly, the NMR shows apparent DS traits when compared to the solitary AMR. Animals can self-domesticate when a reduction of the fear response is naturally selected, such as in islands with no predators, or to improve the group’s harmony in cooperative breeding species. The DS may be caused by alterations in the physiology of the neural crest cells (NCC), a transient population of cells that generate a full range of tissues during development. The NCC contribute to organs responsible for transmitting the fear response and various other tissues, including craniofacial bones. Therefore, mutations affecting the NCC can manifest as behavioural and morphological alterations in many structures across the body, as seen in neurocristopathies. We observed that all social AMRs are chisel-tooth diggers, an adaption to hard soils that requires the flattening of the skull. We hypothesise that chisel-tooth digging could impose a selective pressure on the NCC that triggered the DS’s appearance, possibly facilitating the evolution of sociality. Finally, we discuss how DS traits are neutral or beneficial for the subterranean niche, strategies to test this hypothesis and report well-studied mutations in the NMR that are associated with the NCC physiology or with the control of the fear response. In conclusion, we argue that many of the NMR’s unconventional traits are compatible with the DS and provide a hypothesis about its origins. Our model proposes a novel avenue to enhance the understanding of the extraordinary biology of the NMR.

Comparative analysis of wild-type accessions reveals novel determinants of Arabidopsis seed longevity

Understanding the genetic factors involved in seed longevity is of paramount importance in agricultural and ecological contexts. The polygenic nature of this trait suggests that many of them remain undiscovered. Here, we exploited the contrasting seed longevity found amongst Arabidopsis thaliana accessions to further understand this phenomenon. Concentrations of glutathione were higher in longer-lived than shorter-lived accessions, supporting that redox poise plays a prominent role in seed longevity. However, high seed permeability, normally associated with shorter longevity, is also present in long-lived accessions. Dry seed transcriptome analysis indicated that the contribution to longevity of stored messenger RNA (mRNAs) is complex, including mainly accession-specific mechanisms. The detrimental effect on longevity caused by other factors may be counterbalanced by higher levels of specific mRNAs stored in dry seeds, for instance those of heat-shock proteins. Indeed, loss-of-function mutant analysis demonstrated that heat-shock factors HSF1A and 1B contributed to longevity. Furthermore, mutants of the stress-granule zinc-finger protein TZF9 or the spliceosome subunits MOS4 or MAC3A/MAC3B, extended seed longevity, positioning RNA as a novel player in the regulation of seed viability. mRNAs of proteins with putative relevance to longevity were also abundant in shorter-lived accessions, reinforcing the idea that resistance to ageing is determined by multiple factors.

Thirty years of geometric morphometrics: Achievements, challenges, and the ongoing quest for biological meaningfulness

The foundations of geometric morphometrics were worked out about 30 years ago and have continually been refined and extended. What has remained as a central thrust and source of debate in the morphometrics community is the shared goal of meaningful biological inference through a tight connection between biological theory, measurement, multivariate biostatistics, and geometry. Here we review the building blocks of modern geometric morphometrics: the representation of organismal geometry by landmarks and semilandmarks, the computation of shape or form variables via superimposition, the visualization of statistical results as actual shapes or forms, the decomposition of shape variation into symmetric and asymmetric components and into different spatial scales, the interpretation of various geometries in shape or form space, and models of the association between shape or form and other variables, such as environmental, genetic, or behavioral data. We focus on recent developments and current methodological challenges, especially those arising from the increasing number of landmarks and semilandmarks, and emphasize the importance of thorough exploratory multivariate analyses rather than single scalar summary statistics. We outline promising directions for further research and for the evaluation of new developments, such as "landmark-free" approaches. To illustrate these methods, we analyze three-dimensional human face shape based on data from the Avon Longitudinal Study of Parents and Children (ALSPAC).

Gene loss and compensatory evolution promotes the emergence of morphological novelties in budding yeast

Deleterious mutations are generally considered to be irrelevant for morphological evolution. However, they could be compensated by conditionally beneficial mutations, thereby providing access to new adaptive paths. Here we use high-dimensional phenotyping of laboratory-evolved budding yeast lineages to demonstrate that new cellular morphologies emerge exceptionally rapidly as a by-product of gene loss and subsequent compensatory evolution. Unexpectedly, the capacities for invasive growth, multicellular aggregation and biofilm formation also spontaneously evolve in response to gene loss. These multicellular phenotypes can be achieved by diverse mutational routes and without reactivating the canonical regulatory pathways. These ecologically and clinically relevant traits originate as pleiotropic side effects of compensatory evolution and have no obvious utility in the laboratory environment. The extent of morphological diversity in the evolved lineages is comparable to that of natural yeast isolates with diverse genetic backgrounds and lifestyles. Finally, we show that both the initial gene loss and subsequent compensatory mutations contribute to new morphologies, with their synergistic effects underlying specific morphological changes. We conclude that compensatory evolution is a previously unrecognized source of morphological diversity and phenotypic novelties.

System drift and speciation

Even if a species' phenotype does not change over evolutionary time, the underlying mechanism may change, as distinct molecular pathways can realize identical phenotypes. Here we use linear system theory to explore the consequences of this idea, describing how a gene network underlying a conserved phenotype evolves, as the genetic drift of small changes to these molecular pathways causes a population to explore the set of mechanisms with identical phenotypes. To do this, we model an organism's internal state as a linear system of differential equations for which the environment provides input and the phenotype is the output, in which context there exists an exact characterization of the set of all mechanisms that give the same input-output relationship. This characterization implies that selectively neutral directions in genotype space should be common and that the evolutionary exploration of these distinct but equivalent mechanisms can lead to the reproductive incompatibility of independently evolving populations. This evolutionary exploration, or system drift, is expected to proceed at a rate proportional to the amount of intrapopulation genetic variation divided by the effective population size ( Ne$N_e$ ). At biologically reasonable parameter values this could lead to substantial interpopulation incompatibility, and thus speciation, on a time scale of Ne$N_e$ generations. This model also naturally predicts Haldane's rule, thus providing a concrete explanation of why heterogametic hybrids tend to be disrupted more often than homogametes during the early stages of speciation.

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.

Differences in Germination of ACCase-Resistant Biotypes Containing Isoleucine-1781-Leucine Mutation and Susceptible Biotypes of Wild Oat (Avena sterilis ssp. ludoviciana)

Herbicide resistance can affect seed germination and the optimal conditions required for seed germination, which in turn may impose a fitness cost in resistant populations. Winter wild oat [Avena sterilis L. ssp. ludoviciana (Durieu) Gillet and Magne] is a serious weed in cereal fields. In this study, the molecular basis of resistance to an ACCase herbicide, clodinafop-propargyl, in four A. ludoviciana biotypes was assessed. Germination differences between susceptible (S) and ACCase-resistant biotypes (WR₁, WR₂, WR₃, WR₄) and the effect of Isoleucine-1781-Leucine mutation on germination were also investigated through germination models. The results indicated that WR₁ and WR₄ were very highly resistant (RI > 214.22) to clodinafop-propargyl-contained Isoleucine to Leucine amino acid substitution. However, Isoleucine-1781-Leucine mutation was not detected in other very highly resistant biotypes. Germination studies indicated that resistant biotypes (in particular WR₁ and WR₄) had higher base water potentials than the susceptible one. This shows that resistant biotypes need more soil water to initiate their germination. However, the hydrotime constant for germination was higher in resistant biotypes than in the susceptible one in most cases, showing faster germination in susceptible biotypes. ACCase-resistant biotypes containing the Isoleucine-1781-Leucine mutation had lower seed weight but used more seed reserve to produce seedlings. Hence, integrated management practices such as stale seedbed and implementing it at the right time could be used to take advantage of the differential soil water requirement and relatively late germination characteristics of ACCase-resistant biotypes.

Functional genomics of abiotic environmental adaptation in lacertid lizards and other vertebrates

Understanding the genomic basis of adaptation to different abiotic environments is important in the context of climate change and resulting short-term environmental fluctuations. Using functional and comparative genomics approaches, we here investigated whether signatures of genomic adaptation to a set of environmental parameters are concentrated in specific subsets of genes and functions in lacertid lizards and other vertebrates. We first identify 200 genes with signatures of positive diversifying selection from transcriptomes of 24 species of lacertid lizards and demonstrate their involvement in physiological and morphological adaptations to climate. To understand how functionally similar these genes are to previously predicted candidate functions for climate adaptation and to compare them with other vertebrate species, we then performed a meta-analysis of 1,100 genes under selection obtained from -omics studies in vertebrate species adapted to different abiotic factors. We found that the vertebrate gene set formed a tightly connected interactome, which was to 23% enriched in previously predicted functions of adaptation to climate, and to a large part (18%) involved in organismal stress response. We found a much higher degree of identical genes being repeatedly selected among different animal groups (43.6%), and of functional similarity and post-translational modifications than expected by chance, and no clear functional division between genes used for ectotherm and endotherm physiological strategies. In total, 171 out of 200 genes of Lacertidae were part of this network. These results highlight an important role of a comparatively small set of genes and their functions in environmental adaptation and narrow the set of candidate pathways and markers to be used in future research on adaptation and stress response related to climate change.

Favorable pleiotropic loci for fiber yield and quality in upland cotton (Gossypium hirsutum)

Upland cotton (Gossypium hirsutum L.) is an important economic crop for renewable textile fibers. However, the simultaneous improvement of yield and fiber quality in cotton is difficult as the linkage drag. Compared with breaking the linkage drag, identification of the favorable pleiotropic loci on the genome level by genome-wide association study (GWAS) provides a new way to improve the yield and fiber quality simultaneously. In our study restriction-site-associated DNA sequencing (RAD-seq) was used to genotype 316 cotton accessions. Eight major traits in three categories including yield, fiber quality and maturation were investigated in nine environments (3 sites × 3 years). 231 SNPs associated with these eight traits (− log₁₀(P) > 5.27) were identified, located in 27 genomic regions respectively by linkage disequilibrium analysis. Further analysis showed that four genomic regions (the region 1, 6, 8 and 23) held favorable pleiotropic loci and 6 candidate genes were identified. Through genotyping, 14 elite accessions carrying the favorable loci on four pleiotropic regions were identified. These favorable pleiotropic loci and elite genotypes identified in this study will be utilized to improve the yield and fiber quality simultaneously in future cotton breeding.

Molecular evolutionary insights from PRLR in mammals

Prolactin (PRL) is a pleiotropic neurohormone secreted by the mammalian pituitary gland into the blood, thus reaching many tissues and organs beyond the brain. PRL binds to its receptor, PRLR, eliciting a molecular signaling cascade. This system modulates essential mammalian behaviors and promotes notable modifications in the reproductive female tissues and organs. Here, we explore how the intracellular domain of PRLR (PRLR-ICD) modulates the expression of the PRLR gene. Despite differences in the reproductive strategies between eutherian and metatherian mammals, there is no clear distinction between PRLR-ICD functional motifs. However, we found selection signatures that showed differences between groups, with many conserved functional elements strongly maintained through purifying selection across the class Mammalia. We observed a few residues under relaxed selection, the levels of which were more pronounced in Eutheria and particularly striking in primates (Simiiformes), which could represent a pre-adaptive genetic element protected from purifying selection. Alternative, new motifs, such as YLDP (318-321) and others with residues Y283 and Y290, may already be functional. These motifs would have been co-opted in primates as part of a complex genetic repertoire related to some derived adaptive phenotypes, but these changes would have no impact on the primordial functions that characterize the mammals as a whole and that are related to the PRL-PRLR system.

Morphometric Variation at Different Spatial Scales: Coordination and Compensation in the Emergence of Organismal Form

It is a classic aim of quantitative and evolutionary biology to infer genetic architecture and potential evolutionary responses to selection from the variance–covariance structure of measured traits. But a meaningful genetic or developmental interpretation of raw covariances is difficult, and classic concepts of morphological integration do not directly apply to modern morphometric data. Here, we present a new morphometric strategy based on the comparison of morphological variation across different spatial scales. If anatomical elements vary completely independently, then their variance accumulates at larger scales or for structures composed of multiple elements: morphological variance would be a power function of spatial scale. Deviations from this pattern of “variational self-similarity” (serving as a null model of completely uncoordinated growth) indicate genetic or developmental coregulation of anatomical components. We present biometric strategies and R scripts for identifying patterns of coordination and compensation in the size and shape of composite anatomical structures. In an application to human cranial variation, we found that coordinated variation and positive correlations are prevalent for the size of cranial components, whereas their shape was dominated by compensatory variation, leading to strong canalization of cranial shape at larger scales. We propose that mechanically induced bone formation and remodeling are key mechanisms underlying compensatory variation in cranial shape. Such epigenetic coordination and compensation of growth are indispensable for stable, canalized development and may also foster the evolvability of complex anatomical structures by preserving spatial and functional integrity during genetic responses to selection.[Cranial shape; developmental canalization; evolvability; morphological integration; morphometrics; phenotypic variation; self-similarity.]

Speciation and the developmental alarm clock

New species arise as the genomes of populations diverge. The developmental 'alarm clock' of speciation sounds off when sufficient divergence in genetic control of development leads hybrid individuals to infertility or inviability, the world awoken to the dawn of new species with intrinsic post-zygotic reproductive isolation. Some developmental stages will be more prone to hybrid dysfunction due to how molecular evolution interacts with the ontogenetic timing of gene expression. Considering the ontogeny of hybrid incompatibilities provides a profitable connection between 'evo-devo' and speciation genetics to better link macroevolutionary pattern, microevolutionary process, and molecular mechanisms. Here, we explore speciation alongside development, emphasizing their mutual dependence on genetic network features, fitness landscapes, and developmental system drift. We assess models for how ontogenetic timing of reproductive isolation can be predictable. Experiments and theory within this synthetic perspective can help identify new rules of speciation as well as rules in the molecular evolution of development.

Experimental Hybridization Studies Suggest That Pleiotropic Alleles Commonly Underlie Adaptive Divergence between Natural Populations

The alleles used for adaptation can pleiotropically affect traits under stabilizing selection. The fixation of alleles with deleterious pleiotropic side effects causes compensatory alleles to be favored by selection. Such compensatory alleles might segregate in interpopulation hybrids, resulting in segregation variance for traits where parents have indistinguishable phenotypes. If adaptation typically involves pleiotropy and compensation, then the segregation variance for traits under stabilizing selection is expected to increase with the magnitude of adaptive phenotypic divergence between parents. This prediction has not been tested empirically, and I gathered data from experimental hybridization studies to evaluate it. I found that pairs of parents that are more phenotypically divergent beget hybrids with more segregation variance in traits for which the parents are statistically indistinguishable. This result suggests that adaptive divergence between pairs of natural populations proceeds via pleiotropy and compensation and that deleterious transgressive segregation variance accumulates systematically as populations diverge.

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.

Evolution of the Mammalian Ear: An Evolvability Hypothesis

Encapsulated within the temporal bone and comprising the smallest elements of the vertebrate skeleton, the ear is key to multiple senses: balance, posture control, gaze stabilization, and hearing. The transformation of the primary jaw joint into the mammalian ear ossicles is one of the most iconic transitions in vertebrate evolution, but the drivers of this complex evolutionary trajectory are not fully understood. We propose a novel hypothesis: The incorporation of the bones of the primary jaw joint into the middle ear has considerably increased the genetic, regulatory, and developmental complexity of the mammalian ear. This increase in the number of genetic and developmental factors may, in turn, have increased the evolutionary degrees of freedom for independent adaptations of the different functional ear units. The simpler ear anatomy in birds and reptiles may be less susceptible to developmental instabilities and disorders than in mammals but also more constrained in its evolution. Despite the tight spatial entanglement of functional ear components, the increased "evolvability" of the mammalian ear may have contributed to the evolutionary success and adaptive diversification of mammals in the vast diversity of ecological and behavioral niches observable today. A brief literature review revealed supporting evidence for this hypothesis.

Aligning functional network constraint to evolutionary outcomes

BACKGROUND

Functional constraint through genomic architecture is suggested to be an important dimension of genome evolution, but quantitative evidence for this idea is rare. In this contribution, existing evidence and discussions on genomic architecture as constraint for convergent evolution, rapid adaptation, and genic adaptation are summarized into alternative, testable hypotheses. Network architecture statistics from protein-protein interaction networks are then used to calculate differences in evolutionary outcomes on the example of genomic evolution in yeast, and the results are used to evaluate statistical support for these longstanding hypotheses.

RESULTS

A discriminant function analysis lent statistical support to classifying the yeast interactome into hub, intermediate and peripheral nodes based on network neighborhood connectivity, betweenness centrality, and average shortest path length. Quantitative support for the existence of genomic architecture as a mechanistic basis for evolutionary constraint is then revealed through utilizing these statistical parameters of the protein-protein interaction network in combination with estimators of protein evolution.

CONCLUSIONS

As functional genetic networks are becoming increasingly available, it will now be possible to evaluate functional genetic network constraint against variables describing complex phenotypes and environments, for better understanding of commonly observed deterministic patterns of evolution in non-model organisms. The hypothesis framework and methodological approach outlined herein may help to quantify the extrinsic versus intrinsic dimensions of evolutionary constraint, and result in a better understanding of how fast, effectively, or deterministically organisms adapt.

Rapid Evolution of Gained Essential Developmental Functions of a Young Gene via Interactions with Other Essential Genes

New genes are of recent origin and only present in a subset of species in a phylogeny. Accumulated evidence suggests that new genes, like old genes that are conserved across species, can also take on important functions and be essential for the survival and reproductive success of organisms. Although there are detailed analyses of the mechanisms underlying new genes' gaining fertility functions, how new genes rapidly become essential for viability remains unclear. We focused on a young retro-duplicated gene (CG7804, which we named Cocoon) in Drosophila that originated between 4 and 10 Ma. We found that, unlike its evolutionarily conserved parental gene, Cocoon has evolved under positive selection and accumulated many amino acid differences at functional sites from the parental gene. Despite its young age, Cocoon is essential for the survival of Drosophila melanogaster at multiple developmental stages, including the critical embryonic stage, and its expression is essential in different tissues from those of its parental gene. Functional genomic analyses found that Cocoon acquired unique DNA-binding sites and has a contrasting effect on gene expression to that of its parental gene. Importantly, Cocoon binding predominantly locates at genes that have other essential functions and/or have multiple gene-gene interactions, suggesting that Cocoon acquired novel essential function to survival through forming interactions that have large impacts on the gene interaction network. Our study is an important step toward deciphering the evolutionary trajectory by which new genes functionally diverge from parental genes and become essential.

Evolutionarily novel genes are expressed in transgenic fish tumors and their orthologs are involved in development of progressive traits in humans

Abstract

Earlier we suggested a new hypothesis of the possible evolutionary role of hereditary tumors (Kozlov, Evolution by tumor Neofunctionalization, 2014), and described a new class of genes – tumor specifically expressed, evolutionarily novel (TSEEN) genes - that are predicted by this hypothesis (Kozlov, Infect Agents Cancer 11:34, 2016). In this paper we studied evolutionarily novel genes expressed in fish tumors after regression, as a model of evolving organs. As evolutionarily novel genes may not yet have organismal functions, we studied the acquisition of new gene functions by comparing fish evolutionarily novel genes with their human orthologs. We found that many genes involved in development of progressive traits in humans (lung, mammary gland, placenta, ventricular septum, etc.) originated in fish and are expressed in fish tumors and tumors after regression. These findings support a possible evolutionary role of hereditary tumors, and in particular the hypothesis of evolution by tumor neofunctionalization.

Research highlights

Earlier we described a new class of genes that are tumor-specifically expressed and evolutionarily novel (TSEEN). As the functions of TSEEN genes are often uncertain, we decided to study TSEEN genes of fishes so that we could trace the appearance of their new functions in higher vertebrates. We found that many human genes which are involved in development of progressive traits (placenta development, mammary gland and lung development etc.,) originated in fishes and are expressed in fish tumors.

Developmental Systems Drift and the Drivers of Sex Chromosome Evolution

Phenotypic invariance—the outcome of purifying selection—is a hallmark of biological importance. However, invariant phenotypes might be controlled by diverged genetic systems in different species. Here, we explore how an important and invariant phenotype—the development of sexually differentiated individuals—is controlled in over two dozen species in the frog family Pipidae. We uncovered evidence in different species for 1) an ancestral W chromosome that is not found in many females and is found in some males, 2) independent losses and 3) autosomal segregation of this W chromosome, 4) changes in male versus female heterogamy, and 5) substantial variation among species in recombination suppression on sex chromosomes. We further provide evidence of, and evolutionary context for, the origins of at least seven distinct systems for regulating sex determination among three closely related genera. These systems are distinct in their genomic locations, evolutionary origins, and/or male versus female heterogamy. Our findings demonstrate that the developmental control of sexual differentiation changed via loss, sidelining, and empowerment of a mechanistically influential gene, and offer insights into novel factors that impinge on the diverse evolutionary fates of sex chromosomes.

Non-consumptive effects of predation: does perceived risk strengthen the genetic integration of behaviour and morphology in stickleback?

Predators can shape genetic correlations in prey by altering prey perception of risk. We manipulated perceived risk to test whether such non-consumptive effects tightened behavioural trait correlations in wild-caught stickleback from high- compared to low-risk environments due to genetic variation in plasticity. We expected tighter genetic correlations within perceived risk treatments than across them, and tighter genetic correlations in high-risk than in low-risk treatments. We identified genetic variation in plasticity, with genetic correlations between boldness, sociality, and antipredator morphology, as expected, being tighter within treatments than across them, for both of two populations. By contrast, genetic correlations did not tighten with exposure to risk. Tighter phenotypic correlations in wild stickleback may thus arise because predators induce correlational selection on environmental components of these traits, or because predators tighten residual correlations by causing environmental heterogeneity that is controlled in the laboratory. Our study places phenotypic integration firmly into an ecological context.

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.

Compensatory trans-regulatory alleles minimizing variation in TDH3 expression are common within Saccharomyces cerevisiae

Heritable variation in gene expression is common within species. Much of this variation is due to genetic differences outside of the gene with altered expression and is trans-acting. This trans-regulatory variation is often polygenic, with individual variants typically having small effects, making the genetic architecture and evolution of trans-regulatory variation challenging to study. Consequently, key questions about trans-regulatory variation remain, including the variability of trans-regulatory variation within a species, how selection affects trans-regulatory variation, and how trans-regulatory variants are distributed throughout the genome and within a species. To address these questions, we isolated and measured trans-regulatory differences affecting TDH3 promoter activity among 56 strains of Saccharomyces cerevisiae, finding that trans-regulatory backgrounds varied approximately twofold in their effects on TDH3 promoter activity. Comparing this variation to neutral models of trans-regulatory evolution based on empirical measures of mutational effects revealed that despite this variability in the effects of trans-regulatory backgrounds, stabilizing selection has constrained trans-regulatory differences within this species. Using a powerful quantitative trait locus mapping method, we identified ∼100 trans-acting expression quantitative trait locus in each of three crosses to a common reference strain, indicating that regulatory variation is more polygenic than previous studies have suggested. Loci altering expression were located throughout the genome, and many loci were strain specific. This distribution and prevalence of alleles is consistent with recent theories about the genetic architecture of complex traits. In all mapping experiments, the nonreference strain alleles increased and decreased TDH3 promoter activity with similar frequencies, suggesting that stabilizing selection maintained many trans-acting variants with opposing effects. This variation may provide the raw material for compensatory evolution and larger scale regulatory rewiring observed in developmental systems drift among species.

Investigating the evolution and development of biological complexity under the framework of epigenetics

Biological complexity is a key component of evolvability, yet its study has been hampered by a focus on evolutionary trends of complexification and inconsistent definitions. Here, we demonstrate the utility of bringing complexity into the framework of epigenetics to better investigate its utility as a concept in evolutionary biology. We first analyze the existing metrics of complexity and explore the link between complexity and adaptation. Although recently developed metrics allow for a unified framework, they omit developmental mechanisms. We argue that a better approach to the empirical study of complexity and its evolution includes developmental mechanisms. We then consider epigenetic mechanisms and their role in shaping developmental and evolutionary trajectories, as well as the development and organization of complexity. We argue that epigenetics itself could have emerged from complexity because of a need to self‐regulate. Finally, we explore hybridization complexes and hybrid organisms as potential models for studying the association between epigenetics and complexity. Our goal is not to explain trends in biological complexity but to help develop and elucidate novel questions in the investigation of biological complexity and its evolution.

This manuscript argues that biological complexity is better understood under the framework of epigenetics and that the epigenetic interactions emerge from the self‐regulation of complex systems. Hybrids are offered as models to study these properties.

The making of elaborate petals in Nigella through developmental repatterning

Elaborate petals are present in many flowering plants lineages and have greatly promoted the success and evolutionary radiation of these groups. How elaborate petals are made, however, remains largely unclear. Petals of Nigella (Ranunculaceae) have long been recognized as elaborate and can thus be an excellent model for the study of petal elaboration. Here, by conducting detailed morphological, micromorphological, anatomical, developmental and evolutionary studies on the petals of Nigella species, we explored the processes, general patterns and underlying mechanisms of petal elaboration. We found that petals of Nigella are highly complex, and the complexity can be reflected at various levels. We also found that evolutionary elaboration of the Nigella petals is a gradual process, involving not only modifications of pre-existing structures but also de novo origination of new characters. Further investigations indicated that the elaboration and diversification of Nigella petals were accomplished by modifying the ancestral trajectory of petal development, a process known as developmental repatterning. Our results not only provide new insights into the development and evolution of elaborate petals, but also highlight the necessity of conducting multiple-level investigations for understanding the processes, patterns and underlying mechanisms of plant evolution.

Regulation and evolution of muscle development in tunicates

For more than a century, studies on tunicate muscle formation have revealed many principles of cell fate specification, gene regulation, morphogenesis, and evolution. Here, we review the key studies that have probed the development of all the various muscle cell types in a wide variety of tunicate species. We seize this occasion to explore the implications and questions raised by these findings in the broader context of muscle evolution in chordates.

CCR5-∆32 is deleterious in the homozygous state in humans

We use the genotyping and death register information of 409,693 individuals of British ancestry to investigate fitness effects of the CCR5-∆32 mutation. We estimate a 21% increase in the all-cause mortality rate in individuals who are homozygous for the ∆32 allele. A deleterious effect of the ∆32/∆32 mutation is also independently supported by a significant deviation from the Hardy-Weinberg equilibrium (HWE) due to a deficiency of ∆32/∆32 individuals at the time of recruitment.

Sex-specific dominance reversal of genetic variation for fitness

The maintenance of genetic variance in fitness represents one of the most longstanding enigmas in evolutionary biology. Sexually antagonistic (SA) selection may contribute substantially to maintaining genetic variance in fitness by maintaining alternative alleles with opposite fitness effects in the two sexes. This is especially likely if such SA loci exhibit sex-specific dominance reversal (SSDR)-wherein the allele that benefits a given sex is also dominant in that sex-which would generate balancing selection and maintain stable SA polymorphisms for fitness. However, direct empirical tests of SSDR for fitness are currently lacking. Here, we performed a full diallel cross among isogenic strains derived from a natural population of the seed beetle Callosobruchus maculatus that is known to exhibit SA genetic variance in fitness. We measured sex-specific competitive lifetime reproductive success (i.e., fitness) in >500 sex-by-genotype F1 combinations and found that segregating genetic variation in fitness exhibited pronounced contributions from dominance variance and sex-specific dominance variance. A closer inspection of the nature of dominance variance revealed that the fixed allelic variation captured within each strain tended to be dominant in one sex but recessive in the other, revealing genome-wide SSDR for SA polymorphisms underlying fitness. Our findings suggest that SA balancing selection could play an underappreciated role in maintaining fitness variance in natural populations.

Heterologous Hsp90 promotes phenotypic diversity through network evolution

Biological processes in living cells are often carried out by gene networks in which signals and reactions are integrated through network hubs. Despite their functional importance, it remains unclear to what extent network hubs are evolvable and how alterations impact long-term evolution. We investigated these issues using heat shock protein 90 (Hsp90), a central hub of proteostasis networks. When native Hsp90 in Saccharomyces cerevisiae cells was replaced by the ortholog from hypersaline-tolerant Yarrowia lipolytica that diverged from S. cerevisiae about 270 million years ago, the cells exhibited improved growth in hypersaline environments but compromised growth in others, indicating functional divergence in Hsp90 between the two yeasts. Laboratory evolution shows that evolved Y. lipolytica-HSP90–carrying S. cerevisiae cells exhibit a wider range of phenotypic variation than cells carrying native Hsp90. Identified beneficial mutations are involved in multiple pathways and are often pleiotropic. Our results show that cells adapt to a heterologous Hsp90 by modifying different subnetworks, facilitating the evolution of phenotypic diversity inaccessible to wild-type cells.

Biological processes in living cells are often carried out by gene networks. Hubs are highly connected network components important for integrating signal inputs and generating responsive functional outputs. Heat shock protein 90 (Hsp90), a versatile hub in the protein homeostasis network, is a molecular chaperone essential for cell viability in all tested eukaryotic cells. In yeast, about a quarter of the expressed proteins are profoundly influenced when Hsp90 activity is reduced. Despite its pivotal role, we found that the function of Hsp90 has diverged between two yeast species, Yarrowia lipolytica and Saccharomyces cerevisiae, which split about 270 million years ago. To understand the impacts and adaptive strategies in cells with an altered network hub, we conducted laboratory evolution experiments using a S. cerevisiae strain in which native Hsp90 is replaced by its counterpart in Y. lipolytica. We observed different fitness gain or loss under various stress conditions in individual evolved clones, suggesting that cells adapted via different evolutionary paths. Genome sequencing and mutation reconstitution experiments show that beneficial mutations occurred in multiple Hsp90-related pathways that interact with each other. Our results show that a perturbed network allows cells to evolve a broader range of phenotypic diversity unavailable to wild-type cells.

Revisiting the evolutionary origins of obesity: lazy versus peppy-thrifty genotype hypothesis

The recent global obesity epidemic is attributed to major societal and environmental changes, such as excessive energy intake and sedentary lifestyle. However, exposure to 'obesogenic' environments does not necessarily result in obesity at the individual level, as 40-75% of body mass index variation in population is attributed to genetic differences. The thrifty genotype theory posits that genetic variants promoting efficient food sequestering and optimal deposition of fat during periods of food abundance were evolutionarily advantageous for the early hunter-gatherer and were positively selected. However, the thrifty genotype is likely too simplistic and fails to provide a justification for the complex distribution of obesity predisposing gene variants and for the broad range of body mass index observed in diverse ethnic groups. This review proposes that gene pleiotropy may better account for the variability in the distribution of obesity susceptibility alleles across modern populations. We outline the lazy-thrifty versus peppy-thrifty genotype hypothesis and detail the body of evidence in the literature in support of this novel concept. Future population genetics and mathematical modelling studies that account for pleiotropy may further improve our understanding of the evolutionary origins of the current obesity epidemic.

The emerging structure of the Extended Evolutionary Synthesis: where does Evo-Devo fit in?

The Extended Evolutionary Synthesis (EES) debate is gaining ground in contemporary evolutionary biology. In parallel, a number of philosophical standpoints have emerged in an attempt to clarify what exactly is represented by the EES. For Massimo Pigliucci, we are in the wake of the newest instantiation of a persisting Kuhnian paradigm; in contrast, Telmo Pievani has contended that the transition to an EES could be best represented as a progressive reformation of a prior Lakatosian scientific research program, with the extension of its Neo-Darwinian core and the addition of a brand-new protective belt of assumptions and auxiliary hypotheses. Here, we argue that those philosophical vantage points are not the only ways to interpret what current proposals to 'extend' the Modern Synthesis-derived 'standard evolutionary theory' (SET) entail in terms of theoretical change in evolutionary biology. We specifically propose the image of the emergent EES as a vast network of models and interweaved representations that, instantiated in diverse practices, are connected and related in multiple ways. Under that assumption, the EES could be articulated around a paraconsistent network of evolutionary theories (including some elements of the SET), as well as models, practices and representation systems of contemporary evolutionary biology, with edges and nodes that change their position and centrality as a consequence of the co-construction and stabilization of facts and historical discussions revolving around the epistemic goals of this area of the life sciences. We then critically examine the purported structure of the EES-published by Laland and collaborators in 2015-in light of our own network-based proposal. Finally, we consider which epistemic units of Evo-Devo are present or still missing from the EES, in preparation for further analyses of the topic of explanatory integration in this conceptual framework.

Stress-induced evolution of herbicide resistance and related pleiotropic effects

Herbicide-resistant weeds, especially those with resistance to multiple herbicides, represent a growing worldwide threat to agriculture and food security. Natural selection for resistant genotypes may act on standing genetic variation, or on a genetic and physiological background that is fundamentally altered because of stress responses to sublethal herbicide exposure. Stress-induced changes include DNA mutations, epigenetic alterations, transcriptional remodeling, and protein modifications, all of which can lead to herbicide resistance and a wide range of pleiotropic effects. Resistance selected in this manner is termed systemic acquired herbicide resistance, and the associated pleiotropic effects are manifested as a suite of constitutive transcriptional and post-translational changes related to biotic and abiotic stress adaptation, representing the evolutionary signature of selection. This phenotype is being investigated in two multiple herbicide-resistant populations of the hexaploid, self-pollinating weedy monocot Avena fatua that display such changes as well as constitutive reductions in certain heat shock proteins and their transcripts, which are well known as global regulators of diverse stress adaptation pathways. Herbicide-resistant populations of most weedy plant species exhibit pleiotropic effects, and their association with resistance genes presents a fertile area of investigation. This review proposes that more detailed studies of resistant A. fatua and other species through the lens of plant evolution under stress will inform improved resistant weed prevention and management strategies. © 2018 Society of Chemical Industry.

Developmental system drift in motor ganglion patterning between distantly related tunicates

BACKGROUND

The larval nervous system of the solitary tunicate Ciona is a simple model for the study of chordate neurodevelopment. The development and connectivity of the Ciona motor ganglion have been studied in fine detail, but how this important structure develops in other tunicates is not well known.

METHODS AND RESULTS

By comparing gene expression patterns in the developing MG of the distantly related tunicate Molgula occidentalis, we found that its patterning is highly conserved compared to the Ciona MG. MG neuronal subtypes in Molgula were specified in the exact same positions as in Ciona, though the timing of subtype-specific gene expression onset was slightly shifted to begin earlier, relative to mitotic exit and differentiation. In transgenic Molgula embryos electroporated with Dmbx reporter plasmids, we were also able to characterize the morphology of the lone pair of descending decussating neurons (ddNs) in Molgula, revealing the same unique contralateral projection seen in Ciona ddNs and their putative vertebrate homologs the Mauthner cells. Although Dmbx expression labels the ddNs in both species, cross-species transgenic assays revealed significant changes to the regulatory logic underlying Dmbx transcription. We found that Dmbx cis-regulatory DNAs from Ciona can drive highly specific reporter gene expression in Molgula ddNs, but Molgula sequences are not active in Ciona ddNs.

CONCLUSIONS

This acute divergence in the molecular mechanisms that underlie otherwise functionally conserved cis-regulatory DNAs supports the recently proposed idea that the extreme genetic plasticity observed in tunicates may be attributed to the extreme rigidity of the spatial organization of their embryonic cell lineages.

Gene duplicates resolving sexual conflict rapidly evolved essential gametogenesis functions

Males and females have different fitness optima but share the vast majority of their genomes, causing an inherent genetic conflict between the two sexes that must be resolved to achieve maximal population fitness. We show that two tandem duplicate genes found specifically in Drosophila melanogaster are sexually antagonistic, but rapidly evolved sex-specific functions and expression patterns that mitigate their antagonistic effects. We use copy-specific knockouts and rescue experiments to show that Apollo (Apl) is essential for male fertility but detrimental to female fertility, in addition to its important role in development, while Artemis (Arts) is essential for female fertility but detrimental to male fertility. Further analyses show that Apl and Arts have essential roles in spermatogenesis and oogenesis. These duplicates formed ~200,000 years ago, underwent a strong selective sweep and lost most expression in the antagonized sex. These data provide direct evidence that gene duplication allowed rapid mitigation of sexual conflict by allowing Apl and Arts to evolve essential sex-specific reproductive functions and complementary expression in male and female gonads.

A damped oscillator imposes temporal order on posterior gap gene expression in Drosophila

Insects determine their body segments in two different ways. Short-germband insects, such as the flour beetle Tribolium castaneum, use a molecular clock to establish segments sequentially. In contrast, long-germband insects, such as the vinegar fly Drosophila melanogaster, determine all segments simultaneously through a hierarchical cascade of gene regulation. Gap genes constitute the first layer of the Drosophila segmentation gene hierarchy, downstream of maternal gradients such as that of Caudal (Cad). We use data-driven mathematical modelling and phase space analysis to show that shifting gap domains in the posterior half of the Drosophila embryo are an emergent property of a robust damped oscillator mechanism, suggesting that the regulatory dynamics underlying long- and short-germband segmentation are much more similar than previously thought. In Tribolium, Cad has been proposed to modulate the frequency of the segmentation oscillator. Surprisingly, our simulations and experiments show that the shift rate of posterior gap domains is independent of maternal Cad levels in Drosophila. Our results suggest a novel evolutionary scenario for the short- to long-germband transition and help explain why this transition occurred convergently multiple times during the radiation of the holometabolan insects.

Different insect species exhibit one of two distinct modes of determining their body segments (known as segmentation) during development: they either use a molecular oscillator to position segments sequentially, or they generate segments simultaneously through a hierarchical gene-regulatory cascade. The sequential mode is ancestral, while the simultaneous mode has been derived from it independently several times during evolution. In this paper, we present evidence suggesting that simultaneous segmentation also involves an oscillator in the posterior end of the embryo of the vinegar fly, Drosophila melanogaster. This surprising result indicates that both modes of segment determination are much more similar than previously thought. Such similarity provides an important step towards our understanding of the frequent evolutionary transitions observed between sequential and simultaneous segmentation.