The integrative biology of genetic dominance

An explanation for the prevalence of XY over ZW sex determination in species derived from hermaphroditism

The many independent transitions from hermaphroditism to separate sexes (dioecy) in flowering plants and some animal clades must often have involved the emergence of a heterogametic sex-determining locus, the basis of XY and ZW sex determination (i.e., male and female heterogamety). Current estimates indicate that XY sex determination is much more frequent than ZW, but the reasons for this asymmetry are unclear. One proposition is that separate sexes evolve through the invasion of sterility mutations at closely linked loci, in which case XY sex determination evolves if the initial male sterility mutation is fully recessive. Alternatively, dioecy may evolve via the gradual divergence of male and female phenotypes, but the genetic basis of such divergence and its connection to XY and ZW systems remain poorly understood. Using mathematical modeling, we show how dioecy with XY or ZW sex determination can emerge from the joint evolution of resource allocation to male and female function with its genetic architecture. Our model reveals that whether XY or ZW sex determination evolves depends on the trade-off between allocation to male and female function, and on the mating system of the ancestral hermaphrodites, with selection for female specialization or inbreeding avoidance both favoring XY sex determination. Together, our results cast light on an important but poorly understood path from hermaphroditism to dioecy, and provide an adaptive hypothesis for the preponderance of XY systems. Beyond sex and sex determination, our model shows how ecology can influence the way selection shapes the genetic architecture of polymorphic traits.

Consequences of partially recessive deleterious genetic variation for the evolution of inversions suppressing recombination between sex chromosomes1

The evolution of suppressed recombination between sex chromosomes is widely hypothesized to be driven by sexually antagonistic selection (SA), where tighter linkage between the sex-determining gene(s) and nearby SA loci is favored when it couples male-beneficial alleles to the proto-Y chromosome, and female-beneficial alleles to the proto-X. Although difficult to test empirically, the SA selection hypothesis overshadows several alternatives, including an incomplete but often-repeated "sheltering" hypothesis which suggests that expansion of the sex-linked region (SLR) reduces the homozygous expression of deleterious mutations at selected loci. Here, we use population genetic models to evaluate the consequences of partially recessive deleterious mutational variation for the evolution of otherwise neutral chromosomal inversions expanding the SLR on proto-Y chromosomes. Both autosomal and SLR-expanding inversions face a race against time: lightly-loaded inversions are initially beneficial, but eventually become deleterious as they accumulate new mutations, after which their chances of fixing become negligible. In contrast, initially unloaded inversions eventually become neutral as their deleterious load reaches the same equilibrium as non-inverted haplotypes. Despite the differences in inheritance and indirect selection, SLR-expanding inversions exhibit similar evolutionary dynamics to autosomal inversions over many biologically plausible parameter conditions. Differences emerge when the population average mutation load is quite high; in this case large autosomal inversions that are lucky enough to be mutation-free can rise to intermediate to high frequencies where selection in homozygotes becomes important (Y-linked inversions never appear as homozygous karyotypes); conditions requiring either high mutation rates, highly recessive deleterious mutations, weak selection, or a combination thereof.

Dominance and multi-locus interaction

Dominance is usually considered a constant value that describes the relative difference in fitness or phenotype between heterozygotes and the average of homozygotes at a focal polymorphic locus. However, the observed dominance can vary with the genetic background of the focal locus. Here, alleles at other loci modify the observed phenotype through position effects or dominance modifiers that are sometimes associated with pathogen resistance, lineage, sex, or mating type. Theoretical models have illustrated how variable dominance appears in the context of multi-locus interaction (epistasis). Here, we review empirical evidence for variable dominance and how the observed patterns may be captured by proposed epistatic models. We highlight how integrating epistasis and dominance is crucial for comprehensively understanding adaptation and speciation.

The effect of divergent and parallel selection on the genomic landscape of divergence

While the role of selection in divergence along the speciation continuum is theoretically well understood, defining specific signatures of selection in the genomic landscape of divergence is empirically challenging. Modelling approaches can provide insight into the potential role of selection on the emergence of a heterogenous genomic landscape of divergence. Here, we extend and apply an individual-based approach that simulates the phenotypic and genotypic distributions of two populations under a variety of selection regimes, genotype-phenotype maps, modes of migration, and genotype-environment interactions. We show that genomic islands of high differentiation and genomic valleys of similarity may respectively form under divergent and parallel selection between populations. For both types of between-population selection, negative and positive frequency-dependent selection within populations generated genomic islands of higher magnitude and genomic valleys of similarity, respectively. Divergence rates decreased under strong dominance with divergent selection, as well as in models including genotype-environment interactions under parallel selection. For both divergent and parallel selection models, divergence rate was higher under an intermittent migration regime between populations, in contrast to a constant level of migration across generations, despite an equal number of total migrants. We highlight that interpreting a particular evolutionary history from an observed genomic pattern must be done cautiously, as similar patterns may be obtained from different combinations of evolutionary processes. Modelling approaches such as ours provide an opportunity to narrow the potential routes that generate the genomic patterns of specific evolutionary histories.

Does the definition of a novel environment affect the ability to detect cryptic genetic variation?

Anthropogenic change exposes populations to environments that have been rare or entirely absent from their evolutionary past. Such novel environments are hypothesized to release cryptic genetic variation, a hidden store of variance that can fuel evolution. However, support for this hypothesis is mixed. One possible reason is a lack of clarity in what is meant by 'novel environment', an umbrella term encompassing conditions with potentially contrasting effects on the exposure or concealment of cryptic variation. Here, we use a meta-analysis approach to investigate changes in the total genetic variance of multivariate traits in ancestral versus novel environments. To determine whether the definition of a novel environment could explain the mixed support for a release of cryptic genetic variation, we compared absolute novel environments, those not represented in a population's evolutionary past, to extreme novel environments, those involving frequency or magnitude changes to environments present in a population's ancestry. Despite sufficient statistical power, we detected no broad-scale pattern of increased genetic variance in novel environments, and finding the type of novel environment did not explain any significant variation in effect sizes. When effect sizes were partitioned by experimental design, we found increased genetic variation in studies based on broad-sense measures of variance, and decreased variation in narrow-sense studies, in support of previous research. Therefore, the source of genetic variance, not the definition of a novel environment, was key to understanding environment-dependant genetic variation, highlighting non-additive genetic variance as an important component of cryptic genetic variation and avenue for future research.

The genotype-phenotype relationship and evolutionary genetics in the light of the Metabolic Control Analysis

Metabolic control analysis has long been used as a systemic model of the genotype-phenotype (GP) relationship. By considering kinetic parameters and enzyme concentrations as reflecting the genotype level and metabolic fluxes or pools as phenotypes related to fitness, MCA has given a biological basis to the relationship between these two levels. The non-linear and concave relationship between enzymes and fluxes can account for common genetic effects that reductionist approaches have been powerless to explain, such as the dominance of active alleles over less active alleles, the various types of epistasis and heterosis, and reveals the structural links between these genetic effects. The summation property of the flux control coefficients accounts for the L-shaped distribution of Quantitative Trait Locus (QTL) effects, irrespective of other possible causes. Metabolic models of response to selection results in evolutionary scenarios that are markedly different from those derived from the classical infinitesimal model of quantitative genetics. In particular, evolution towards selective neutrality appears to be a consequence of the diminishing return of the flux-enzyme relationship. In this paper, we survey the historical and recent achievements of MCA in genetics, quantitative genetics and evolution, focusing on epistasis and the evolution of flux in relation to enzyme concentrations.

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.

Analysis of genetic dominance in the UK Biobank

Classical statistical genetics theory defines dominance as any deviation from a purely additive, or dosage, effect of a genotype on a trait, which is known as the dominance deviation. Dominance is well documented in plant and animal breeding. Outside of rare monogenic traits, however, evidence in humans is limited. We systematically examined common genetic variation across 1060 traits in a large population cohort (UK Biobank, N = 361,194 samples analyzed) for evidence of dominance effects. We then developed a computationally efficient method to rapidly assess the aggregate contribution of dominance deviations to heritability. Lastly, observing that dominance associations are inherently less correlated between sites at a genomic locus than their additive counterparts, we explored whether they may be leveraged to identify causal variants more confidently.

Using inbreeding to test the contribution of non-additive genetic effects to additive genetic variance: a case study in Drosophila serrata

Additive genetic variance, V_A, is the key parameter for predicting adaptive and neutral phenotypic evolution. Changes in demography (e.g. increased close-relative inbreeding) can alter V_A, but how they do so depends on the (typically unknown) gene action and allele frequencies across many loci. For example, V_A increases proportionally with the inbreeding coefficient when allelic effects are additive, but smaller (or larger) increases can occur when allele frequencies are unequal at causal loci with dominance effects. Here, we describe an experimental approach to assess the potential for dominance effects to deflate V_A under inbreeding. Applying a powerful paired pedigree design in Drosophila serrata, we measured 11 wing traits on half-sibling families bred via either random or sibling mating, differing only in homozygosity (not allele frequency). Despite close inbreeding and substantial power to detect small V_A, we detected no deviation from the expected additive effect of inbreeding on genetic (co)variances. Our results suggest the average dominance coefficient is very small relative to the additive effect, or that allele frequencies are relatively equal at loci affecting wing traits. We outline the further opportunities for this paired pedigree approach to reveal the characteristics of V_A, providing insight into historical selection and future evolutionary potential.

Allele frequencies and selection coefficients in locally adapted populations

Understanding the role of natural selection in driving evolutionary change requires accurate estimates of the strength of selection acting at the genetic level in the wild. This is challenging to achieve but may be easier in the case of populations in migration-selection balance. When two populations are at equilibrium under migration-selection balance, there exist loci whose alleles are selected different ways in the two populations. Such loci can be identified from genome sequencing by their high values of F_ST. This raises the question of what is the strength of selection on locally-adaptive alleles. To answer this question we analyse a 1-locus 2-allele model of a population distributed between two niches. We show by simulation of selected cases that the outputs from finite-population models are essentially the same as those from deterministic infinite-population models. We then derive theory for the infinite-population model showing the dependence of selection coefficients on equilibrium allele frequencies, migration rates, dominance and relative population sizes in the two niches. An Excel spreadsheet is provided for the calculation of selection coefficients and their approximate standard errors from observed values of population parameters. We illustrate our results with a worked example, with graphs showing the dependence of selection coefficients on equilibrium allele frequencies, and graphs showing how F_ST depends on the selection coefficients acting on the alleles at a locus. Given the extent of recent progress in ecological genomics, we hope our methods may help those studying migration-selection balance to quantify the advantages conferred by adaptive genes.

Dominance mechanisms in supergene alleles controlling butterfly wing pattern variation: insights from gene expression in Heliconius numata

Loci under balancing selection, where multiple alleles are maintained, offer a relevant opportunity to investigate the role of natural selection in shaping genetic dominance: the high frequency of heterozygotes at these loci has been shown to enable the evolution of dominance among alleles. In the butterfly Heliconius numata, mimetic wing color variations are controlled by an inversion polymorphism of a circa 2 Mb genomic region (supergene P), with strong dominance between sympatric alleles. To test how differences in dominance observed on wing patterns correlate with variations in expression levels throughout the supergene region, we sequenced the complete transcriptome of heterozygotes at the prepupal stage and compared it to corresponding homozygotes. By defining dominance based on non-overlapping ranges of transcript expression between genotypes, we found contrasting patterns of dominance between the supergene and the rest of the genome; the patterns of transcript expression in the heterozygotes were more similar to the expression observed in the dominant homozygotes in the supergene region. Dominance also differed among the three subinversions of the supergene, suggesting possible epistatic interactions among their gene contents underlying dominance evolution. We found the expression pattern of the melanization gene cortex located in the P-region to predict wing pattern phenotype in the heterozygote. We also identify new candidate genes that are potentially involved in mimetic color pattern variations highlighting the relevance of transcriptomic analyses in heterozygotes to pinpoint candidate genes in non-recombining regions.

An exploration into the conversion of dominance to additive genetic variance in contrasting environments

PREMISE

The evolutionary response of a trait to environmental change depends upon the level of additive genetic variance. It has been long argued that sustained selection will tend to deplete additive genetic variance as favored alleles approach fixation. Non-additive genetic variance, due to interactions among alleles within and between loci, does not immediately contribute to an evolutionary response. However, shifts in the allele frequencies within and between interacting loci may convert non-additive variance into additive variance. Here we consider the possibility that an environmental shift may alter allelic interactions in ways that convert dominance into additive genetic variance.

METHODS

We grew a pedigreed population of Brassica rapa in greenhouse and field conditions. The field conditions mimicked agricultural conditions from which the base population was drawn, while the greenhouse featured benign conditions. We used Bayesian models to estimate the additive, dominance, and maternal components of quantitative genetic variance. We also estimated genetic correlations across environments using parental breeding values.

RESULTS

Although the additive genetic variance was elevated in the greenhouse condition, no consistent pattens emerged that would indicate a conversion of dominance variance. The unusually low genetic variance and broad confidence intervals for the variance estimates obtained through this analysis preclude definitive interpretations.

CONCLUSIONS

Further studies are needed to determine whether between-environment changes in additive genetic variance can be traced to conversion of dominance variance.

Sexual imprinting leads to speciation in locally adapted populations

Sexual imprinting is widespread in birds and other species but its existence requires explanation. Our results suggest that sexual imprinting leads to speciation in locally-adapted populations if a neutral mating cue-e.g., novel plumage coloration-arises through mutation. Importantly, the mating cue locus is not linked to adaptation loci. Local adaptation is a necessary precursor to speciation and occurs when evolution results in stable genetic polymorphisms with one allele predominating in some areas while others predominate elsewhere. Here we use a deterministic two-niche population genetic model to map the set of migration and selection rates for which polymorphic evolutionary outcomes, i.e., local adaptations, can occur. Approximate equations for the boundaries of the set of polymorphic evolutionary outcomes were derived by Bulmer (American Naturalist, 106, 254, 1972), but our results, obtained by deterministic simulation of the evolutionary process, show that one of Bulmer's equations is inaccurate except when the level of dominance is 0.5, and fails if one of the alleles is dominant. Having an accurate map of the set of migration and selection rates for which polymorphic evolutionary outcomes can occur, we then show using the model of Sibly et al. (Ecology and Evolution, 9, 13506, 2019) that local adaptation in all analyzed cases leads to speciation if a new neutral mating cue arises by mutation. We finish by considering how genome sequencing makes possible testing our model and its predictions.

Selection on modifiers of genetic architecture under migration load

Gene flow between populations adapting to differing local environmental conditions might be costly because individuals can disperse to habitats where their survival is low or because they can reproduce with locally maladapted individuals. The amount by which the mean relative population fitness is kept below one creates an opportunity for modifiers of the genetic architecture to spread due to selection. Prior work that separately considered modifiers changing dispersal, recombination rates, or altering dominance or epistasis, has typically focused on the direction of selection rather than its absolute magnitude. We here develop methods to determine the strength of selection on modifiers of the genetic architecture, including modifiers of the dispersal rate, in populations that have previously evolved local adaptation. We consider scenarios with up to five loci contributing to local adaptation and derive a new model for the deterministic spread of modifiers. We find that selection for modifiers of epistasis and dominance is stronger than selection for decreased recombination, and that selection for partial reductions in recombination are extremely weak, regardless of the number of loci contributing to local adaptation. The spread of modifiers that reduce dispersal depends on the number of loci, epistasis and extent of local adaptation in the ancestral population. We identify a novel effect, that modifiers of dominance are more strongly selected when they are unlinked to the locus that they modify. These findings help explain population differentiation and reproductive isolation and provide a benchmark to compare selection on modifiers under finite population sizes and demographic stochasticity.

When populations of a species are spread over different habitats the populations can adapt to their local conditions, provided dispersal between habitats is low enough. Natural selection allows the populations to maintain local adaptation, but dispersal and gene flow create a cost called the migration load. The migration load measures how much fitness is lost because of dispersal between different habitats, and also creates an opportunity for selection to act on the arrangement and interaction between genes that are involved in local adaptation. Modifier genes can spread in these linked populations and cause functional, local adaptation genes, to become more closely linked on a chromosome, or change the way that these genes are expressed so that the locally adapted gene copy becomes dominant. We modeled this process and found that selection on modifiers that create tighter linkage between locally adapted genes is generally weak, and modifiers that cause gene interactions are more strongly selected. Even after these gene interactions have begun to evolve, further selection for increased gene interaction is still strong. Our results show that populations are more likely to adapt to local conditions by evolving new gene interactions than by evolving tightly linked gene clusters.

Intrinsic emergence and modulation of sex-specific dominance reversals in threshold traits

Sex-specific dominance reversals (SSDRs) in fitness-related traits, where heterozygotes' phenotypes resemble those of alternative homozygotes in females versus males, can simultaneously maintain genetic variation in fitness and resolve sexual conflict and thereby shape key evolutionary outcomes. However, the full implications of SSDRs will depend on how they arise and the resulting potential for evolutionary, ecological and environmental modulation. Recent field and laboratory studies have demonstrated SSDRs in threshold(-like) traits with dichotomous or competitive phenotypic outcomes, implying that such traits could promote the emergence of SSDRs. However, such possibilities have not been explicitly examined. I show how phenotypic SSDRs can readily emerge in threshold traits given genetic architectures involving large-effect loci alongside sexual dimorphism in the mean and variance in polygenic liability. I also show how multilocus SSDRs can arise in line-cross experiments, especially given competitive reproductive systems that generate nonlinear fitness outcomes. SSDRs can consequently emerge in threshold(-like) traits as functions of sexual antagonism, sexual dimorphism and reproductive systems, even with purely additive underlying genetic effects. Accordingly, I identify theoretical and empirical advances that are now required to discern the basis and occurrence of SSDRs in nature, probe forms of (co-)evolutionary, ecological and environmental modulation, and evaluate net impacts on sexual conflict.

Consequences of partially recessive deleterious genetic variation for the evolution of inversions suppressing recombination between sex chromosomes

The evolution of suppressed recombination between sex chromosomes is widely hypothesized to be driven by sexually antagonistic selection (SA), where tighter linkage between the sex-determining gene(s) and nearby SA loci is favored when it couples male-beneficial alleles to the proto-Y chromosome, and female-beneficial alleles to the proto-X. Despite limited empirical evidence, the SA selection hypothesis overshadows several alternatives, including an incomplete but often-repeated "sheltering hypothesis" that suggests that expansion of the sex-linked region (SLR) reduces homozygous expression of partially recessive deleterious mutations at selected loci. Here, we use population genetic models to evaluate the consequences of deleterious mutational variation for the evolution of neutral chromosomal inversions expanding the SLR on proto-Y chromosomes. We find that SLR-expanding inversions face a race against time: lightly loaded inversions are initially beneficial, but eventually become deleterious as they accumulate new mutations, and must fix before this window of opportunity closes. The outcome of this race is strongly influenced by inversion size, the mutation rate, and the dominance coefficient of deleterious mutations. Yet, small inversions have elevated fixation probabilities relative to neutral expectations for biologically plausible parameter values. Our results demonstrate that deleterious genetic variation can plausibly drive recombination suppression in small steps and would be most consistent with empirical patterns of small evolutionary strata or gradual recombination arrest.