Constraints on the evolution of function‐valued traits: a study of growth in Tribolium castaneum

Growth trajectories often impact individual fitness. They are continuous by nature and so are amenable to analysis using a function‐valued (FV) trait framework to reveal their underlying genetic architecture. Previous studies have found high levels of standing additive genetic (co)variance for growth trajectories despite the expectation that growth should be responding to frequent strong directional selection. In this study, the FV framework is used to estimate the additive genetic covariance function for growth trajectories in larval Tribolium castaneum to address questions about standing additive genetic (co)variance and possible evolutionary constraints on growth and to predict responses to four plausible selection regimes. Results show that additive genetic (co)variance is high at the early ages, but decreases towards later ages in the larval period. A selection gradient function of the same size and in the same direction of the first eigenfunction of the G‐function should give the maximal response. However, evolutionary constraints may be acting to keep this maximal response from being realized, through either conflicting effects on survivability and fecundity of larger body size, few evolutionary directions having sufficient additive variance for a response, genetic trade‐offs with other traits or physiological regulatory mechanisms. More light may be shed on these constraints through the development of more sophisticated statistical approaches and implementation of additional empirical studies to explicitly test for specific types of constraints.     

Gene flow does not prevent personality and morphological differentiation between two blue tit populations

Understanding the causes and consequences of population phenotypic divergence is a central goal in ecology and evolution. Phenotypic divergence among populations can result from genetic divergence, phenotypic plasticity or a combination of the two. However, few studies have deciphered these mechanisms for populations geographically close and connected by gene flow, especially in the case of personality traits. In this study, we used a common garden experiment to explore the genetic basis of the phenotypic divergence observed between two blue tit (Cyanistes caeruleus) populations inhabiting contrasting habitats separated by 25 km, for two personality traits (exploration speed and handling aggression), one physiological trait (heart rate during restraint) and two morphological traits (tarsus length and body mass). Blue tit nestlings were removed from their population and raised in a common garden for up to 5 years. We then compared adult phenotypes between the two populations, as well as trait‐specific Q_st and F_st. Our results revealed differences between populations similar to those found in the wild, suggesting a genetic divergence for all traits. Q_st–F_st comparisons revealed that the trait divergences likely result from dissimilar selection patterns rather than from genetic drift. Our study is one of the first to report a Q_st–F_st comparison for personality traits and adds to the growing body of evidence that population genetic divergence is possible at a small scale for a variety of traits including behavioural traits.     

Phenotypic variation and covariation indicate high evolvability of acoustic communication in crickets

Studying the genetic architecture of sexual traits provides insight into the rate and direction at which traits can respond to selection. Traits associated with few loci and limited genetic and phenotypic constraints tend to evolve at high rates typically observed for secondary sexual characters. Here, we examined the genetic architecture of song traits and female song preferences in the field crickets Gryllus rubens and Gryllus texensis. Song and preference data were collected from both species and interspecific F1 and F2 hybrids. We first analysed phenotypic variation to examine interspecific differentiation and trait distributions in parental and hybrid generations. Then, the relative contribution of additive and additive‐dominance variation was estimated. Finally, phenotypic variance–covariance ( P ) matrices were estimated to evaluate the multivariate phenotype available for selection. Song traits and preferences had unimodal trait distributions, and hybrid offspring were intermediate with respect to the parents. We uncovered additive and dominance variation in song traits and preferences. For two song traits, we found evidence for X‐linked inheritance. On the one hand, the observed genetic architecture does not suggest rapid divergence, although sex linkage may have allowed for somewhat higher evolutionary rates. On the other hand, P matrices revealed that multivariate variation in song traits aligned with major dimensions in song preferences, suggesting a strong selection response. We also found strong covariance between the main traits that are sexually selected and traits that are not directly selected by females, providing an explanation for the striking multivariate divergence in male calling songs despite limited divergence in female preferences.     

Explosive radiation and spatial expansion across the cold environments of the Old World in an avian family

Our objective was to elucidate the biogeography and speciation patterns in an entire avian family, which shows a complex pattern of overlapping and nonoverlapping geographical distributions, and much variation in plumage, but less in size and structure. We estimated the phylogeny and divergence times for all of the world's species of Prunella based on multiple genetic loci, and analyzed morphometric divergence and biogeographical history. The common ancestor of Prunella was present in the Sino‐Himalayan Mountains or these mountains and Central Asia–Mongolia more than 9 million years ago (mya), but a burst of speciations took place during the mid‐Pliocene to early Pleistocene. The relationships among the six primary lineages resulting from that differentiation are unresolved, probably because of the rapid radiation. A general increase in sympatry with increasing time since divergence is evident. With one exception, species in clades younger than c. 3.7 my are allopatric. Species that are widely sympatric, including the most recently diverged (2.4 mya) sympatric sisters, are generally more divergent in size/structure than allo‐/parapatric close relatives. The distributional pattern and inferred ages suggest divergence in allopatry and substantial waiting time until secondary contact, likely due to competitive exclusion. All sympatrically breeding species are ecologically segregated, as suggested by differences in size/structure and habitat. Colonizations of new areas were facilitated during glacial periods, followed by fragmentation during interglacials—contrary to the usual view that glacial periods resulted mainly in fragmentations.     

Divergent selection and drift shape the genomes of two avian sister species spanning a saline–freshwater ecotone

The role of species divergence due to ecologically based divergent selection—or ecological speciation—in generating and maintaining biodiversity is a central question in evolutionary biology. Comparison of the genomes of phylogenetically related taxa spanning a selective habitat gradient enables discovery of divergent signatures of selection and thereby provides valuable insight into the role of divergent ecological selection in speciation. Tidal marsh ecosystems provide tractable opportunities for studying organisms' adaptations to selective pressures that underlie ecological divergence. Sharp environmental gradients across the saline–freshwater ecotone within tidal marshes present extreme adaptive challenges to terrestrial vertebrates. Here, we sequence 20 whole genomes of two avian sister species endemic to tidal marshes—the saltmarsh sparrow (Ammospiza caudacutus) and Nelson's sparrow (A. nelsoni)—to evaluate the influence of selective and demographic processes in shaping genome‐wide patterns of divergence. Genome‐wide divergence between these two recently diverged sister species was notably high (genome‐wide F_ST = 0.32). Against a background of high genome‐wide divergence, regions of elevated divergence were widespread throughout the genome, as opposed to focused within islands of differentiation. These patterns may be the result of genetic drift resulting from past tidal march colonization events in conjunction with divergent selection to different environments. We identified several candidate genes that exhibited elevated divergence between saltmarsh and Nelson's sparrows, including genes linked to osmotic regulation, circadian rhythm, and plumage melanism—all putative candidates linked to adaptation to tidal marsh environments. These findings provide new insights into the roles of divergent selection and genetic drift in generating and maintaining biodiversity.     

Geographic variation in advertisement calls of a Microhylid frog – testing the role of drift and ecology

Acoustic signals for mating are important traits that could drive population differentiation and speciation. Ecology may play a role in acoustic divergence through direct selection (e.g., local adaptation to abiotic environment), constraint of correlated traits (e.g., acoustic traits linked to another trait under selection), and/or interspecific competition (e.g., character displacement). However, genetic drift alone can also drive acoustic divergence. It is not always easy to differentiate the role of ecology versus drift in acoustic divergence. In this study, we tested the role of ecology and drift in shaping geographic variation in the advertisement calls of Microhyla fissipes. We examined three predictions based on ecological processes: (1) the correlation between temperature and call properties across M. fissipes populations; (2) the correlation between call properties and body size across M. fissipes populations; and (3) reproductive character displacement (RCD) in call properties between M. fissipes populations that are sympatric with and allopatric to a congener M. heymonsi. To test genetic drift, we examined correlations among call divergence, geographic distance, and genetic distance across M. fissipes populations. We recorded the advertisement calls from 11 populations of M. fissipes in Taiwan, five of which are sympatrically distributed with M. heymonsi. We found geographic variation in both temporal and spectral properties of the advertisement calls of M. fissipes. However, the call properties were not correlated with local temperature or the callers' body size. Furthermore, we did not detect RCD. By contrast, call divergence, geographic distance, and genetic distance between M. fissipes populations were all positively correlated. The comparisons between phenotypic Q_st (P_st) and F_st values did not show significant differences, suggesting a role of drift. We concluded that genetic drift, rather than ecological processes, is the more likely driver for the geographic variation in the advertisement calls of M. fissipes.     

Are habitat fragmentation,local adaptation and isolation‐by‐distance driving population divergence in wild rice Oryza rufipogon?

Habitat fragmentation weakens the connection between populations and is accompanied with isolation by distance (IBD) and local adaptation (isolation by adaptation, IBA), both leading to genetic divergence between populations. To understand the evolutionary potential of a population and to formulate proper conservation strategies, information on the roles of IBD and IBA in driving population divergence is critical. The putative ancestor of Asian cultivated rice (Oryza sativa) is endangered in China due to habitat loss and fragmentation. We investigated the genetic variation in 11 Chinese Oryza rufipogon populations using 79 microsatellite loci to infer the effects of habitat fragmentation, IBD and IBA on genetic structure. Historical and current gene flows were found to be rare (m_h = 0.0002–0.0013, m_c = 0.007–0.029), indicating IBD and resulting in a high level of population divergence (F_ST = 0.343). High within‐population genetic variation (H_E = 0.377–0.515), relatively large effective population sizes (N_e = 96–158), absence of bottlenecks and limited gene flow were found, demonstrating little impact of recent habitat fragmentation on these populations. Eleven gene‐linked microsatellite loci were identified as outliers, indicating local adaptation. Hierarchical AMOVA and partial Mantel tests indicated that population divergence of Chinese O. rufipogon was significantly correlated with environmental factors, especially habitat temperature. Common garden trials detected a significant adaptive population divergence associated with latitude. Collectively, these findings imply that IBD due to historical rather than recent fragmentation, followed by local adaptation, has driven population divergence in O. rufipogon.     

Conservation of genetic diversity hotspots of the high‐valued relic yellowhorn (Xanthoceras sorbifolium) considering climate change predictions

Genetic structure and major climate factors may contribute to the distribution of genetic diversity of a highly valued oil tree species Xanthoceras sorbifolium (yellowhorn). Long‐term over utilization along with climate change is affecting the viability of yellowhorn wild populations. To preserve the species known and unknown valuable gene pools, the identification of genetic diversity “hotspots” is a prerequisite for their consideration as in situ conservation high priority. Chloroplast DNA (cpDNA) diversity was high among 38 natural populations (H_d = 0.717, K = 4.616, Tajmas’ D = ?0.22) and characterized by high genetic divergence (F_ST = 0.765) and relatively low gene flow (N_m = 0.03), indicating populations isolation reflecting the species’ habitat fragmentation and inbreeding depression. Six out of the studied 38 populations are defined as genetic diversity “hotspots.” The number and geographic direction of cpDNA mutation steps supported the species southwest to northeast migration history. Climatic factors such as extreme minimum temperature over 30 years indicated that the identified genetic “hotspots” are expected to experience 5°C temperature increase in next following 50 years. The results identified vulnerable genetic diversity “hotspots” and provided fundamental information for the species’ future conservation and breeding activities under the anticipated climate change. More specifically, the role of breeding as a component of a gene resource management strategy aimed at fulfilling both utilization and conservation goals.     

Testing for a genetic response to sexual selection in a wild Drosophila population

In accordance with the consensus that sexual selection is responsible for the rapid evolution of display traits on macroevolutionary scales, microevolutionary studies suggest sexual selection is a widespread and often strong form of directional selection in nature. However, empirical evidence for the contemporary evolution of sexually selected traits via sexual rather than natural selection remains weak. In this study, we used a novel application of quantitative genetic breeding designs to test for a genetic response to sexual selection on eight chemical display traits from a field population of the fly, Drosophila serrata. Using our quantitative genetic approach, we were able to detect a genetically based difference in means between groups of males descended from fathers who had either successfully sired offspring or were randomly collected from the same wild population for one of these display traits, the diene (Z,Z)‐5,9‐C_27 : 2. Our experimental results, in combination with previous laboratory studies on this system, suggest that both natural and sexual selection may be influencing the evolutionary trajectories of these traits in nature, limiting the capacity for a contemporary evolutionary response.     

Genetic dissection of the gene coexpression network underlying photosynthesis in Populus

Photosynthesis is a key reaction that ultimately generates the carbohydrates needed to form woody tissues in trees. However, the genetic regulatory network of protein‐encoding genes (PEGs) and regulatory noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs), underlying the photosynthetic pathway is unknown. Here, we integrated data from coexpression analysis, association studies (additive, dominance and epistasis), and expression quantitative trait nucleotide (eQTN) mapping to dissect the causal variants and genetic interaction network underlying photosynthesis in Populus. We initially used 30 PEGs, 6 miRNAs and 12 lncRNAs to construct a coexpression network based on the tissue‐specific gene expression profiles of 15 Populus samples. Then, we performed association studies using a natural population of 435 unrelated Populus tomentosa individuals, and identified 72 significant associations (P ≤ 0.001, q ≤ 0.05) with diverse additive and dominance patterns underlying photosynthesis‐related traits. Analysis of epistasis and eQTNs revealed that the complex genetic interactions in the coexpression network contribute to phenotypes at various levels. Finally, we demonstrated that heterologously expressing the most highly linked gene (PtoPsbX1) in this network significantly improved photosynthesis in Arabidopsis thaliana, pointing to the functional role of PtoPsbX1 in the photosynthetic pathway. This study provides an integrated strategy for dissecting a complex genetic interaction network, which should accelerate marker‐assisted breeding efforts to genetically improve woody plants.     

Divergence in mating signals correlates with genetic distance and behavioural responses to playback

Animals use acoustic signals to defend resources against rivals and attract breeding partners. As with many biological traits, acoustic signals may reflect ancestry; closely related species often produce more similar signals than do distantly related species. Whether this similarity in acoustic signals is biologically relevant to animals is poorly understood. We conducted a playback experiment to measure the physical and vocal responses of male songbirds to the songs of both conspecific and allopatric‐congeneric animals that varied in their acoustic and genetic similarity. Our subjects were territorial males of four species of neotropical Troglodytes wrens: Brown‐throated Wrens (Troglodytes brunneicollis), Cozumel Wrens (T. beani), Clarion Wrens (T. tanneri) and Socorro Wrens (T. sissonii). Our results indicate that birds respond to playback of both conspecific and allopatric‐congeneric animals; that acoustic differences increase with genetic distance; and that genetic divergence predicts the strength of behavioural responses to playback, after removing the effects of acoustic similarity between subjects’ songs and playback stimuli. Collectively, these results demonstrate that the most distantly related species have the most divergent songs; that male wrens perceive divergence in fine structural characteristics of songs; and that perceptual differences between species reflect evolutionary history. This study offers novel insight into the importance of acoustic divergence of learned signals and receiver responses in species recognition.     

Bud phenology and growth are subject to divergent selection across a latitudinal gradient in Populus angustifolia and impact adaptation across the distributional range and associated arthropods

Temperate forest tree species that span large geographical areas and climatic gradients often have high levels of genetic variation. Such species are ideal for testing how neutral demographic factors and climate‐driven selection structure genetic variation within species, and how this genetic variation can affect ecological communities. Here, we quantified genetic variation in vegetative phenology and growth traits in narrowleaf cottonwood, Populus angustifolia, using three common gardens planted with genotypes originating from source populations spanning the species' range along the Rocky Mountains of North America (ca. 1700 km). We present three main findings. First, we found strong evidence of divergent selection (Q_ST > F_ST) on fall phenology (bud set) with adaptive consequences for frost avoidance. We also found evidence for selection on bud flush duration, tree height, and basal diameter, resulting in population differentiation. Second, we found strong associations with climate variables that were strongly correlated with latitude of origin. More strongly differentiated traits also showed stronger climate correlations, which emphasizes the role that climate has played in divergent selection throughout the range. We found population × garden interaction effects; for some traits, this accounted for more of the variance than either factor alone. Tree height was influenced by the difference in climate of the source and garden locations and declined with increasing transfer distance. Third, growth traits were correlated with dependent arthropod community diversity metrics. Synthesis. Overall, we conclude that climate has influenced genetic variation and structure in phenology and growth traits and leads to local adaptation in P. angustifolia, which can then impact dependent arthropod species. Importantly, relocation of genotypes far northward or southward often resulted in poor growth, likely due to a phenological mismatch with photoperiod, the proximate cue for fall growth cessation. Genotypes moved too far southward suffer from early growth cessation, whereas those moved too far northward are prone to fall frost and winter dieback. In the face of current and forecasted climate change, habitat restoration, forestry, and tree breeding efforts should utilize these findings to better match latitudinal and climatic source environments with management locations for optimal future outcomes.     

Genome divergence and diversification within a geographic mosaic of coevolution

Despite substantial interest in coevolution's role in diversification, examples of coevolution contributing to speciation have been elusive. Here, we build upon past studies that have shown both coevolution between South Hills crossbills and lodgepole pine (Pinus contorta), and high levels of reproductive isolation between South Hills crossbills and other ecotypes in the North American red crossbill (Loxia curvirostra) complex. We used genotyping by sequencing to generate population genomic data and applied phylogenetic and population genetic analyses to characterize the genetic structure within and among nine of the ecotypes. Although genome‐wide divergence was slight between ecotypes (F_ST = 0.011–0.035), we found evidence of relative genetic differentiation (as measured by F_ST) between and genetic cohesiveness within many of them. As expected for nomadic and opportunistic breeders, we detected no evidence of isolation by distance. The one sedentary ecotype, the South Hills crossbill, was genetically most distinct because of elevated divergence at a small number of loci rather than pronounced overall genome‐wide divergence. These findings suggest that mechanisms related to recent local coevolution between South Hills crossbills and lodgepole pine (e.g. strong resource‐based density dependence limiting gene flow) have been associated with genome divergence in the face of gene flow. Our results further characterize a striking example of coevolution driving speciation within perhaps as little as 6000 years.     

Intronic variation at the CHD1‐Z gene in Black‐tailed Godwits Limosa limosa limosa: correlations with fitness components revisited

Recently, Schroeder et al. (2010, Ibis 152: 368–377) suggested that intronic variation in the CHD1‐Z gene of Black‐tailed Godwits breeding in southwest Friesland, The Netherlands, correlated with fitness components. Here we re‐examine this surprising result using an expanded dataset (2088 birds sampled from 2004 to 2010 vs. 284 birds from 2004 to 2007). We find that the presence of the Z* allele (9% of the birds) is not associated with breeding habitat type, egg size, adult survival, adult body mass or adult body condition. The results presented here, when used in synergy with the previously reported results by Schroeder et al., suggest that there might be a tendency towards female adults with the Z* allele laying earlier clutches than adult females without the Z* allele. The occurrence of the Z* allele was also associated with a higher chick body mass and return rate. Chicks with the Z* allele that had hatched early in the breeding season were heavier at birth than chicks without the Z* allele and chicks with the Z* allele that had hatched late. Collectively, the results suggest that variation in the CHD1‐Z gene may indeed have arisen as a byproduct of selection acting on females during the egg fase and on chicks during the rearing stages of the reproductive cycle.     

Genetic diversity and phylogeographic structure of Bactrian camels shown by mitochondrial sequence variations

The Bactrian camel includes various domestic (Camelus bactrianus) and wild (Camelus ferus) breeds that are important for transportation and for their nutritional value. However, there is a lack of extensive information on their genetic diversity and phylogeographic structure. Here, we studied these parameters by examining an 809‐bp mtDNA fragment from 113 individuals, representing 11 domestic breeds, one wild breed and two hybrid individuals. We found 15 different haplotypes, and the phylogenetic analysis suggests that domestic and wild Bactrian camels have two distinct lineages. The analysis of molecular variance placed most of the genetic variance (90.14%, P < 0.01) between wild and domestic camel lineages, suggesting that domestic and wild Bactrian camel do not have the same maternal origin. The analysis of domestic Bactrian camels from different geographical locations found there was no significant genetic divergence in China, Russia and Mongolia. This suggests a strong gene flow due to wide movement of domestic Bactrian camels.     

Development of an integrated 200K SNP genotyping array and application for genetic mapping,genome assembly improvement and genome wide association studies in pear (Pyrus)

Pear (Pyrus; 2n = 34), the third most important temperate fruit crop, has great nutritional and economic value. Despite the availability of many genomic resources in pear, it is challenging to genotype novel germplasm resources and breeding progeny in a timely and cost‐effective manner. Genotyping arrays can provide fast, efficient and high‐throughput genetic characterization of diverse germplasm, genetic mapping and breeding populations. We present here 200K AXIOM^® PyrSNP, a large‐scale single nucleotide polymorphism (SNP) genotyping array to facilitate genotyping of Pyrus species. A diverse panel of 113 re‐sequenced pear genotypes was used to discover SNPs to promote increased adoption of the array. A set of 188 diverse accessions and an F₁ population of 98 individuals from ‘Cuiguan’ × ‘Starkrimson’ was genotyped with the array to assess its effectiveness. A large majority of SNPs (166 335 or 83%) are of high quality. The high density and uniform distribution of the array SNPs facilitated prediction of centromeric regions on 17 pear chromosomes, and significantly improved the genome assembly from 75.5% to 81.4% based on genetic mapping. Identification of a gene associated with flowering time and candidate genes linked to size of fruit core via genome wide association studies showed the usefulness of the array in pear genetic research. The newly developed high‐density SNP array presents an important tool for rapid and high‐throughput genotyping in pear for genetic map construction, QTL identification and genomic selection.     

Analysis of the Inheritance, Selection and Evolution of Growth Trajectories

We present methods for estimating the parameters of inheritance and selection that appear in a quantitative genetic model for the evolution growth trajectories and other "infinite-dimensional" traits that we recently introduced. Two methods for estimating the additive genetic covariance function are developed, a "full" model that fully fits the data and a "reduced" model that generates a smoothed estimate consistent with the sampling errors in the data. By decomposing the covariance function into its eigenvalues and eigenfunctions, it is possible to identify potential evolutionary changes in the population's mean growth trajectory for which there is (and those for which there is not) genetic variation. Algorithms for estimating these quantities, their confidence intervals, and for testing hypotheses about them are developed. These techniques are illustrated by an analysis of early growth in mice. Compatible methods for estimating the selection gradient function acting on growth trajectories in natural or domesticated populations are presented. We show how the estimates for the additive genetic covariance function and the selection gradient function can be used to predict the evolutionary change in a population's mean growth trajectory.     

Assessment of genetic diversity and population structure in cultured Australian Pacific oysters

The recent development of Pacific oyster (Crassostrea gigas) SNP genotyping arrays has allowed detailed characterisation of genetic diversity and population structure within and between oyster populations. It also raises the potential of harnessing genomic selection for genetic improvement in oyster breeding programmes. The aim of this study was to characterise a breeding population of Australian oysters through genotyping and analysis of 18 027 SNPs, followed by comparison with genotypes of oyster sampled from Europe and Asia. This revealed that the Australian populations had similar population diversity (H_E) to oysters from New Zealand, the British Isles, France and Japan. Population divergence was assessed using PCA of genetic distance and revealed that Australian oysters were distinct from all other populations tested. Australian Pacific oysters originate from planned introductions sourced from three Japanese populations. Approximately 95% of these introductions were from geographically, and potentially genetically, distinct populations from the Nagasaki oysters assessed in this study. Finally, in preparation for the application of genomic selection in oyster breeding programmes, the strength of LD was evaluated and subsets of loci were tested for their ability to accurately infer relationships. Weak LD was observed on average; however, SNP subsets were shown to accurately reconstitute a genomic relationship matrix constructed using all loci. This suggests that low‐density SNP panels may have utility in the Australian population tested, and the findings represent an important first step towards the design and implementation of genomic approaches for applied breeding in Pacific oysters.     

Body size divergence in nine‐spined sticklebacks: disentangling additive genetic and maternal effects

Interpopulation differences in body size are of common occurrence in vertebrates, but the relative importance of genetic, maternal, and environmental effects as causes of observed differentiation have seldom been assessed in the wild. Gigantism in pond nine‐spined sticklebacks (Pungitius pungitius Linnaeus, 1758) has been repeatedly observed, but the quantitative genetic basis of population divergence in size has remained unstudied. We conducted a common garden experiment – using ‘pure’ and reciprocal crosses between two populations (‘giant’ pond versus ‘normal’ marine) – to test for the relative importance of additive genetic, non‐additive genetic, and maternal effects on body size after 11 months of growth in the laboratory. We found that body size difference between the two populations in laboratory conditions owed mainly to additive genetic effects, and only to a minor degree to maternal effects. Furthermore, the weak maternal effects were seen only in the offspring of ‘giant’ mothers, and appeared to be mediated through differences in egg size. Thus, the results suggest that gigantism in pond populations of P. pungitius is based on the effects of additively acting genes, rather than to direct environmental induction, or maternal or non‐additive gene action. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 107 , 521–528.