Supergenes and complex phenotypes

No mate preference associated with the supergene controlling social organization in Alpine silver ants

Disassortative mating is a powerful mechanism stabilizing polymorphisms at sex chromosomes and other supergenes. The Alpine silver ant, Formica selysi, has two forms of social organization-single-queen and multiple-queen colonies-determined by alternate haplotypes at a large supergene. Here, we explore whether mate preference contributes to the maintenance of the genetic polymorphism at the social supergene. With mate choice experiments, we found that females and males mated randomly with respect to social form. Moreover, queens were able to produce offspring irrespective of whether they had mated with a male from the same or the alternative social form. Yet, females originating from single-queen colonies were more fertile, suggesting that they may be more successful at independent colony founding. We conclude that the pattern of asymmetric assortative mating documented from mature F. selysi colonies in the field is not caused by mate preferences or major genetic incompatibilities between social forms. More generally, we found no evidence that disassortative mate preference contributes to the maintenance of polymorphism at this supergene controlling ant social organization.

Assessing the Architecture of Drosophila mojavensis Locomotor Evolution with Bulk Segregant Analysis

Behavior is frequently predicted to be especially important for evolution in novel environments. If these predictions are accurate, there might be particular patterns of genetic architecture associated with recently diverged behaviors. Specifically, it has been predicted that behaviors linked to population divergence should be underpinned by a few genes of relatively large effect, compared to architectures of intrapopulation behavioral variation, which is considered to be highly polygenic. More mapping studies of behavioral variation between recently diverged populations are needed to continue assessing the generality of these predictions. Here, we used a bulk segregant mapping approach to dissect the genetic architecture of a locomotor trait that has evolved between two populations of the cactophilic fly Drosophila mojavensis We created an F8 mapping population of 1,500 individuals from advanced intercross lines and sequenced the 10% of individuals with the highest and lowest levels of locomotor activity. Using three alternative statistical approaches, we found strong evidence for two relatively large-effect QTL that is localized in a region homologous to a region of densely packed behavior loci in Drosophila melanogaster, suggesting that clustering of behavior genes may display relatively deep evolutionary conservation. Broadly, our data are most consistent with a polygenic architecture, though with several loci explaining a high proportion of variation in comparison to similar behavioral traits. We further note the presence of several antagonistic QTL linked to locomotion and discuss these results in light of theories regarding behavioral evolution and the effect size and direction of QTL for diverging traits in general.

Genetics and evidence for balancing selection of a sex-linked colour polymorphism in a songbird

Colour polymorphisms play a key role in sexual selection and speciation, yet the mechanisms that generate and maintain them are not fully understood. Here, we use genomic and transcriptomic tools to identify the precise genetic architecture and evolutionary history of a sex-linked colour polymorphism in the Gouldian finch Erythrura gouldiae that is also accompanied by remarkable differences in behaviour and physiology. We find that differences in colour are associated with an ~72-kbp region of the Z chromosome in a putative regulatory region for follistatin, an antagonist of the TGF-β superfamily genes. The region is highly differentiated between morphs, unlike the rest of the genome, yet we find no evidence that an inversion is involved in maintaining the distinct haplotypes. Coalescent simulations confirm that there is elevated nucleotide diversity and an excess of intermediate frequency alleles at this locus. We conclude that this pleiotropic colour polymorphism is most probably maintained by balancing selection.

GOOGA: A platform to synthesize mapping experiments and identify genomic structural diversity

Understanding genomic structural variation such as inversions and translocations is a key challenge in evolutionary genetics. We develop a novel statistical approach to comparative genetic mapping to detect large-scale structural mutations from low-level sequencing data. The procedure, called Genome Order Optimization by Genetic Algorithm (GOOGA), couples a Hidden Markov Model with a Genetic Algorithm to analyze data from genetic mapping populations. We demonstrate the method using both simulated data (calibrated from experiments on Drosophila melanogaster) and real data from five distinct crosses within the flowering plant genus Mimulus. Application of GOOGA to the Mimulus data corrects numerous errors (misplaced sequences) in the M. guttatus reference genome and confirms or detects eight large inversions polymorphic within the species complex. Finally, we show how this method can be applied in genomic scans to improve the accuracy and resolution of Quantitative Trait Locus (QTL) mapping.

Regulatory changes in pterin and carotenoid genes underlie balanced color polymorphisms in the wall lizard

Reptiles use pterin and carotenoid pigments to produce yellow, orange, and red colors. These conspicuous colors serve a diversity of signaling functions, but their molecular basis remains unresolved. Here, we show that the genomes of sympatric color morphs of the European common wall lizard (Podarcis muralis), which differ in orange and yellow pigmentation and in their ecology and behavior, are virtually undifferentiated. Genetic differences are restricted to two small regulatory regions near genes associated with pterin [sepiapterin reductase (SPR)] and carotenoid [beta-carotene oxygenase 2 (BCO2)] metabolism, demonstrating that a core gene in the housekeeping pathway of pterin biosynthesis has been coopted for bright coloration in reptiles and indicating that these loci exert pleiotropic effects on other aspects of physiology. Pigmentation differences are explained by extremely divergent alleles, and haplotype analysis revealed abundant transspecific allele sharing with other lacertids exhibiting color polymorphisms. The evolution of these conspicuous color ornaments is the result of ancient genetic variation and cross-species hybridization.

Regulatory changes in pterin and carotenoid genes underlie balanced color polymorphisms in the wall lizard

Reptiles show an amazing color diversity based on variation in melanins, carotenoids, and pterins. This study reveals genes controlling differences between three color morphs (white, orange, and yellow) in the common wall lizard. Orange pigmentation, due to high levels of orange/red pterins in skin, is caused by genetic changes in the sepiapterin reductase gene. Yellow skin, showing high levels of yellow carotenoids, is controlled by the beta-carotene oxygenase 2 locus. Thus, the color polymorphism in the common wall lizard is associated with changes in two small regions of the genome containing genes with crucial roles in pterin and carotenoid metabolism. These genes are likely to have pleiotropic effects on behavior and other traits associated with the different color morphs.

Reptiles use pterin and carotenoid pigments to produce yellow, orange, and red colors. These conspicuous colors serve a diversity of signaling functions, but their molecular basis remains unresolved. Here, we show that the genomes of sympatric color morphs of the European common wall lizard (Podarcis muralis), which differ in orange and yellow pigmentation and in their ecology and behavior, are virtually undifferentiated. Genetic differences are restricted to two small regulatory regions near genes associated with pterin [sepiapterin reductase (SPR)] and carotenoid [beta-carotene oxygenase 2 (BCO2)] metabolism, demonstrating that a core gene in the housekeeping pathway of pterin biosynthesis has been coopted for bright coloration in reptiles and indicating that these loci exert pleiotropic effects on other aspects of physiology. Pigmentation differences are explained by extremely divergent alleles, and haplotype analysis revealed abundant transspecific allele sharing with other lacertids exhibiting color polymorphisms. The evolution of these conspicuous color ornaments is the result of ancient genetic variation and cross-species hybridization.

A selfish genetic element linked to increased lifespan impacts metabolism in female house mice

Gene drive systems can lead to the evolution of traits that further enhance the transmission of the driving element. In gene drive, one allele is transmitted to offspring at a higher frequency than the homologous allele. This has a range of consequences, which generally include a reduction in fitness of the carrier of the driving allele, making such systems “selfish”. The t haplotype is one such driver, found in house mice. It is linked to a reduction in litter size in matings among heterozygous animals, but also to increased lifespan in wild females that carry it. Here, we tested whether carrying the t haplotype was associated with altered resting metabolic rate (RMR). We show that females carrying the t haplotype decrease RMR as they increase in size, compared to wildtype females or males of either genotype. Our study elucidates a plausible mechanism by which a selfish genetic element increases lifespan.

Unexpected patterns of segregation distortion at a selfish supergene in the fire ant Solenopsis invicta

BACKGROUND

The Sb supergene in the fire ant Solenopsis invicta determines the form of colony social organization, with colonies whose inhabitants bear the element containing multiple reproductive queens and colonies lacking it containing only a single queen. Several features of this supergene - including suppressed recombination, presence of deleterious mutations, association with a large centromere, and "green-beard" behavior - suggest that it may be a selfish genetic element that engages in transmission ratio distortion (TRD), defined as significant departures in progeny allele frequencies from Mendelian inheritance ratios. We tested this possibility by surveying segregation ratios in embryo progenies of 101 queens of the "polygyne" social form (3512 embryos) using three supergene-linked markers and twelve markers outside the supergene.

RESULTS

Significant departures from Mendelian ratios were observed at the supergene loci in 3-5 times more progenies than expected in the absence of TRD and than found, on average, among non-supergene loci. Also, supergene loci displayed the greatest mean deviations from Mendelian ratios among all study loci, although these typically were modest. A surprising feature of the observed inter-progeny variation in TRD was that significant deviations involved not only excesses of supergene alleles but also similarly frequent excesses of the alternate alleles on the homologous chromosome. As expected given the common occurrence of such "drive reversal" in this system, alleles associated with the supergene gain no consistent transmission advantage over their alternate alleles at the population level. Finally, we observed low levels of recombination and incomplete gametic disequilibrium across the supergene, including between adjacent markers within a single inversion.

CONCLUSIONS

Our data confirm the prediction that the Sb supergene is a selfish genetic element capable of biasing its own transmission during reproduction, yet counterselection for suppressor loci evidently has produced an evolutionary stalemate in TRD between the variant homologous haplotypes on the "social chromosome". Evidence implicates prezygotic segregation distortion as responsible for the TRD we document, with "true" meiotic drive the most likely mechanism. Low levels of recombination and incomplete gametic disequilibrium across the supergene suggest that selection does not preserve a single uniform supergene haplotype responsible for inducing polygyny.

Behavior-dependent cis regulation reveals genes and pathways associated with bower building in cichlid fishes

Many behaviors are associated with heritable genetic variation [Kendler and Greenspan (2006) Am J Psychiatry 163:1683-1694]. Genetic mapping has revealed genomic regions or, in a few cases, specific genes explaining part of this variation [Bendesky and Bargmann (2011) Nat Rev Gen 12:809-820]. However, the genetic basis of behavioral evolution remains unclear. Here we investigate the evolution of an innate extended phenotype, bower building, among cichlid fishes of Lake Malawi. Males build bowers of two types, pits or castles, to attract females for mating. We performed comparative genome-wide analyses of 20 bower-building species and found that these phenotypes have evolved multiple times with thousands of genetic variants strongly associated with this behavior, suggesting a polygenic architecture. Remarkably, F₁ hybrids of a pit-digging and a castle-building species perform sequential construction of first a pit and then a castle bower. Analysis of brain gene expression in these hybrids showed that genes near behavior-associated variants display behavior-dependent allele-specific expression with preferential expression of the pit-digging species allele during pit digging and of the castle-building species allele during castle building. These genes are highly enriched for functions related to neurodevelopment and neural plasticity. Our results suggest that natural behaviors are associated with complex genetic architectures that alter behavior via cis-regulatory differences whose effects on gene expression are specific to the behavior itself.

Modularity of genes involved in local adaptation to climate despite physical linkage

BACKGROUND

Linkage among genes experiencing different selection pressures can make natural selection less efficient. Theory predicts that when local adaptation is driven by complex and non-covarying stresses, increased linkage is favored for alleles with similar pleiotropic effects, with increased recombination favored among alleles with contrasting pleiotropic effects. Here, we introduce a framework to test these predictions with a co-association network analysis, which clusters loci based on differing associations. We use this framework to study the genetic architecture of local adaptation to climate in lodgepole pine, Pinus contorta, based on associations with environments.

RESULTS

We identify many clusters of candidate genes and SNPs associated with distinct environments, including aspects of aridity and freezing, and discover low recombination rates among some candidate genes in different clusters. Only a few genes contain SNPs with effects on more than one distinct aspect of climate. There is limited correspondence between co-association networks and gene regulatory networks. We further show how associations with environmental principal components can lead to misinterpretation. Finally, simulations illustrate both benefits and caveats of co-association networks.

CONCLUSIONS

Our results support the prediction that different selection pressures favor the evolution of distinct groups of genes, each associating with a different aspect of climate. But our results went against the prediction that loci experiencing different sources of selection would have high recombination among them. These results give new insight into evolutionary debates about the extent of modularity, pleiotropy, and linkage in the evolution of genetic architectures.

Identification and validation of reference genes for qPCR in the terrestrial gastropod Cepaea nemoralis

Identifying and understanding mechanisms that generate phenotypic diversity is a fundamental goal of evolutionary biology. With a diversity of pigmented shell morphotypes governed by Mendelian patterns of inheritance, the common grove snail Cepaea nemoralis (Linnaeus, 1758) has been a model for evolutionary biologists and population geneticists for decades. However, the genetic mechanisms by which C. nemoralis generates this pigmented shell diversity remain unknown. An important first step in investigating this pigmentation pattern is to establish a set of validated reference genes for differential gene expression assays. Here we have evaluated eleven candidate genes for reverse transcription quantitative polymerase chain reaction (qPCR) in C. nemoralis. Five of these were housekeeping genes traditionally employed as qPCR reference genes in other species, while six alternative genes were selected de novo from C. nemoralis transcriptome data based on the stability of their expression levels. We tested all eleven candidates for expression stability in four sub-adult tissues of C. nemoralis: pigmented mantle, unpigmented mantle, head and foot. We find that two commonly employed housekeeping genes (alpha tubulin, glyceraldehyde 3-phosphate dehydrogenase) are unsuitable for use as qPCR reference genes in C. nemoralis. The traditional housekeeping gene UBIquitin on the other hand performed very well. Additionally, an RNA-directed DNA polymerase (RNAP), a Potassium Channel Protein (KCHP) and a Prenylated Rab acceptor protein 1 (PRAP), identified de novo from transcriptomic data, were the most stably expressed genes in different tissue combinations. We also tested expression stability over two seasons and found that, although other genes are more stable within a single season, beta actin (BACT) and elongation factor 1 alpha (EF1α) were the most reliable reference genes across seasons.

Keystone Genes

The keystone species concept is used in ecology to describe individual species with disproportionately large effects on their communities. We extend this idea to the level of genes with disproportionately large effects on ecological processes. Such 'keystone genes' (KGs) would underlie traits involved in species interactions or causing critical biotic and/or abiotic changes that influence emergent community and ecosystem properties. We propose a general framework for how KGs could be identified, while keeping KGs under the umbrella of 'ecologically important genes' (EIGs) that also include categories such as 'foundation genes', 'ecosystem engineering genes', and more. Although likely rare, KGs and other EIGs could dominate certain ecological processes; thus, their discovery and study are relevant for understanding eco-evolutionary dynamics.

Asymmetric assortative mating and queen polyandry are linked to a supergene controlling ant social organization

Nonrecombining genomic variants underlie spectacular social polymorphisms, from bird mating systems to ant social organization. Because these "social supergenes" affect multiple phenotypic traits linked to survival and reproduction, explaining their persistence remains a substantial challenge. Here, we investigate how large nonrecombining genomic variants relate to colony social organization, mating system and dispersal in the Alpine silver ant, Formica selysi. The species has colonies headed by a single queen (monogynous) and colonies headed by multiple queens (polygynous). We confirmed that a supergene with alternate haplotypes-Sm and Sp-underlies this polymorphism in social structure: Females from mature monogynous colonies had the Sm/Sm genotype, while those from polygynous colonies were Sm/Sp and Sp/Sp. Queens heading monogynous colonies were exclusively mated with Sm males. In contrast, queens heading polygynous colonies were mated with Sp males and Sm males. Sm males, which are only produced by monogynous colonies, accounted for 22.9% of the matings with queens from mature polygynous colonies. This asymmetry between social forms in the degree of assortative mating generates unidirectional male-mediated gene flow from the monogynous to the polygynous social form. Biased gene flow was confirmed by a significantly higher number of private alleles in the polygynous social form. Moreover, heterozygous queens were three times as likely as homozygous queens to be multiply mated. This study reveals that the supergene variants jointly affect social organization and multiple components of the mating system that alter the transmission of the variants and thus influence the dynamics of the system.

The Genetics of a Behavioral Speciation Phenotype in an Island System

Mating behavior divergence can make significant contributions to reproductive isolation and speciation in various biogeographic contexts. However, whether the genetic architecture underlying mating behavior divergence is related to the biogeographic history and the tempo and mode of speciation remains poorly understood. Here, we use quantitative trait locus (QTL) mapping to infer the number, distribution, and effect size of mating song rhythm variations in the crickets Laupala eukolea and Laupala cerasina, which occur on different islands (Maui and Hawaii). We then compare these results with a similar study of an independently evolving species pair that diverged within the same island. Finally, we annotate the L. cerasina transcriptome and test whether the QTL fall in functionally enriched genomic regions. We document a polygenic architecture behind the song rhythm divergence in the inter-island species pair that is remarkably similar to that previously found for an intra-island species pair in the same genus. Importantly, the QTL regions were significantly enriched for potential homologs of the genes involved in pathways that may be modulating the cricket song rhythm. These clusters of loci could constrain the spatial genomic distribution of the genetic variation underlying the cricket song variation and harbor several candidate genes that merit further study.

Micro- and macroevolutionary adaptation through repeated loss of a complete metabolic pathway

There is growing evidence for the convergent evolution of physically linked gene clusters encoding chemical defense pathways. Metabolic clusters are proposed to evolve because they ensure co-inheritance of all required genes where the defense is favored, and prevent inheritance of toxic partial pathways where it is not. This hypothesis rests on the assumption that clusters evolve in species where selection favors intraspecific polymorphism for the defense; however, they have not been examined in polymorphic species. We examined metabolic cluster evolution in relation to an adaptive polymorphism for cyanogenic glucoside (CNglc) production in clover. Using 163 accessions, we performed CNglc assays, BAC sequencing, Southern hybridizations and molecular evolutionary analyses. We find that the CNglc pathway forms a 138-kb cluster in white clover, and that the adaptive polymorphism occurs through presence/absence of the complete cluster. Component genes are orthologous to those in the distantly related legume Lotus japonicus. These findings provide empirical support for the co-inheritance hypothesis, and they indicate that adaptive CNglc variation in white clover evolves through recurrent deletions of the entire pathway. They further indicate that the shared ancestor of many important legume crops was likely cyanogenic and that this defense was lost repeatedly over the last 50 Myr.

An inversion supergene in Drosophila underpins latitudinal clines in survival traits

Chromosomal inversions often contribute to local adaptation across latitudinal clines, but the underlying selective mechanisms remain poorly understood. We and others have previously shown that a clinal inversion polymorphism in Drosophila melanogaster, In(3R)Payne, underpins body size clines along the North American and Australian east coasts. Here, we ask whether this polymorphism also contributes to clinal variation in other fitness-related traits, namely survival traits (lifespan, survival upon starvation and survival upon cold shock). We generated homokaryon lines, either carrying the inverted or standard chromosomal arrangement, isolated from populations approximating the endpoints of the North American cline (Florida, Maine) and phenotyped the flies at two growth temperatures (18 °C, 25 °C). Across both temperatures, high-latitude flies from Maine lived longer and were more stress resistant than low-latitude flies from Florida, as previously observed. Interestingly, we find that this latitudinal pattern is partly explained by the clinal distribution of the In(3R)P polymorphism, which is at ~ 50% frequency in Florida but absent in Maine: inverted karyotypes tended to be shorter-lived and less stress resistant than uninverted karyotypes. We also detected an interaction between karyotype and temperature on survival traits. As In(3R)P influences body size and multiple survival traits, it can be viewed as a 'supergene', a cluster of tightly linked loci affecting multiple complex phenotypes. We conjecture that the inversion cline is maintained by fitness trade-offs and balancing selection across geography; elucidating the mechanisms whereby this inversion affects alternative, locally adapted phenotypes across the cline is an important task for future work.

Toward Understanding the Repeated Occurrence of Associations between Melanin-Based Coloration and Multiple Phenotypes

Melanin is the most widespread pigment in organisms. Melanin-based coloration has been repeatedly observed to be associated with the same traits and in the same direction in different vertebrate and insect species. However, whether any factors that are common to different taxa account for the repeated evolution of melanin-phenotype associations remains unclear. We propose to approach this question from the perspective of convergent and parallel evolution to clarify to what extent different species have evolved the same associations owing to a shared genetic basis and being subjected to similar selective pressures. Our current understanding of the genetic basis of melanin-phenotype associations allows for both convergent and parallel evolution, but this understanding is still limited. Further research is needed to clarify the generality and interdependencies of the different proposed mechanisms (supergenes, pleiotropy based on hormones, or neural crest cells). The general ecological scenarios whereby melanin-based coloration is under selection-protection from ultraviolet radiation, thermoregulation in cold environments, or as a signal of social status-offer a good opportunity to study how melanin-phenotype associations evolve. Reviewing these scenarios shows that some traits associated with melanin-based coloration might be selected together with coloration by also favoring adaptation but that other associated traits might impede adaptation, which may be indicative of genetic constraints. We therefore encourage further research on the relative roles that selection and genetic constraints play in shaping multiple melanin-phenotype associations. Placed into a phylogenetic context, this will help clarify to what extent these associations result from convergent or parallel evolutionary processes and why melanin-phenotype associations are so common across the tree of life.

Multiple large inversions and breakpoint rewiring of gene expression in the evolution of the fire ant social supergene

Supergenes consist of co-adapted loci that segregate together and are associated with adaptive traits. In the fire ant Solenopsis invicta, two 'social' supergene variants regulate differences in colony queen number and other traits. Suppressed recombination in this system is maintained, in part, by a greater than 9 Mb inversion, but the supergene is larger. Has the supergene in S. invicta undergone multiple large inversions? The initial gene content of the inverted allele of a supergene would be the same as that of the wild-type allele. So, how did the inversion increase in frequency? To address these questions, we cloned one extreme breakpoint in the fire ant supergene. In doing so, we found a second large (greater than 800 Kb) rearrangement. Furthermore, we determined the temporal order of the two big inversions based on the translocation pattern of a third small fragment. Because the S. invicta supergene lacks evolutionary strata, our finding of multiple inversions may support an introgression model of the supergene. Finally, we showed that one of the inversions swapped the promoter of a breakpoint-adjacent gene, which might have conferred a selective advantage relative to the non-inverted allele. Our findings provide a rare example of gene alterations arising directly from an inversion event.

Parallel evolution of Batesian mimicry supergene in two Papilio butterflies, P. polytes and P. memnon

Batesian mimicry protects animals from predators when mimics resemble distasteful models. The female-limited Batesian mimicry in Papilio butterflies is controlled by a supergene locus switching mimetic and nonmimetic forms. In Papilio polytes, recent studies revealed that a highly diversified region (HDR) containing doublesex (dsx-HDR) constitutes the supergene with dimorphic alleles and is likely maintained by a chromosomal inversion. In the closely related Papilio memnon, which exhibits a similar mimicry polymorphism, we performed whole-genome sequence analyses in 11 butterflies, which revealed a nearly identical dsx-HDR containing three genes (dsx, Nach-like, and UXT) with dimorphic sequences strictly associated with the mimetic/nonmimetic phenotypes. In addition, expression of these genes, except that of Nach-like in female hind wings, showed differences correlated with phenotype. The dimorphic dsx-HDR in P. memnon is maintained without a chromosomal inversion, suggesting that a separate mechanism causes and maintains allelic divergence in these genes. More abundant accumulation of transposable elements and repetitive sequences in the dsx-HDR than in other genomic regions may contribute to the suppression of chromosomal recombination. Gene trees for Dsx, Nach-like, and UXT indicated that mimetic alleles evolved independently in the two Papilio species. These results suggest that the genomic region involving the above three genes has repeatedly diverged so that two allelic sequences of this region function as developmental switches for mimicry polymorphism in the two Papilio species. The supergene structures revealed here suggest that independent evolutionary processes with different genetic mechanisms have led to parallel evolution of similar female-limited polymorphisms underlying Batesian mimicry in Papilio butterflies.

Genetic Dissection of a Supergene Implicates Tfap2a in Craniofacial Evolution of Threespine Sticklebacks

In nature, multiple adaptive phenotypes often coevolve and can be controlled by tightly linked genetic loci known as supergenes. Dissecting the genetic basis of these linked phenotypes is a major challenge in evolutionary genetics. Multiple freshwater populations of threespine stickleback fish (Gasterosteus aculeatus) have convergently evolved two constructive craniofacial traits, longer branchial bones and increased pharyngeal tooth number, likely as adaptations to dietary differences between marine and freshwater environments. Prior QTL mapping showed that both traits are partially controlled by overlapping genomic regions on chromosome 21 and that a regulatory change in Bmp6 likely underlies the tooth number QTL. Here, we mapped the branchial bone length QTL to a 155 kb, eight-gene interval tightly linked to, but excluding the coding regions of Bmp6 and containing the candidate gene Tfap2a Further recombinant mapping revealed this bone length QTL is separable into at least two loci. During embryonic and larval development, Tfap2a was expressed in the branchial bone primordia, where allele specific expression assays revealed the freshwater allele of Tfap2a was expressed at lower levels relative to the marine allele in hybrid fish. Induced loss-of-function mutations in Tfap2a revealed an essential role in stickleback craniofacial development and show that bone length is sensitive to Tfap2a dosage in heterozygotes. Combined, these results suggest that closely linked but genetically separable changes in Bmp6 and Tfap2a contribute to a supergene underlying evolved skeletal gain in multiple freshwater stickleback populations.

Negative-Assortative Mating in the White-Throated Sparrow

Nonrandom mating based on phenotype has been observed in a number of organisms, but a very high proportion of these examples are of assortative mating. The strongest example of negative-assortative mating is for white-striped versus tan-striped crown in the white-throated sparrow, where about 98% of the observed pairings (mated pairs or social pairs) are between mates with different phenotypes and the correlation between mating types is -0.964. Although nonrandom mating has been explored theoretically for decades, these models have generally not focused on specific well-documented examples. Here we have developed a model to investigate the dynamics and equilibrium of this iconic example. The observed pattern of mating appears to be the result of 96% negative-assortative mating and a 17% advantage of W (white) male × T (tan) female matings compared to the reciprocal T male × W female matings. The equilibrium heterozygosity given these values is 0.500, very close to the 0.501 observed in our large sample of pairings, and this heterozygosity has been maintained for the 29 years from 1988 to 2016. In addition, the equilibrium frequency of 2m inversion determining the white-striped phenotype has been maintained at a frequency very close to its equilibrium frequency of 0.25. Overall, this model demonstrates how evolutionary genetic models can be used to understand negative-assortative mating.

Darwin endures, despite disparagement

Evolution lies at the heart of the life sciences, and Charles Darwin is a towering historical figure within evolutionary science. One testimony to his lasting influence is that declaring Darwin to have been wrong all along remains a provocative way to command attention. The present paper discusses various strands of 'Darwin was wrong' partisans and their divergent views and motives: some are looking to Darwin to justify or condemn the political ideologies that they support or reject; others are concerned with the corrupting influence that the bleak cosmic outlook of evolution is deemed to exert on the moral or religious rectitude of impressionable minds, or regard Darwinism as a direct assault on religion; philosophers question the very coherence of the entire enterprise; and certain biologists aspire to go down in history as even greater than Darwin. It is sobering to reflect that this diverse group is united only by their poor grasp of Darwin's theory of natural selection.

Rapid regulatory evolution of a nonrecombining autosome linked to divergent behavioral phenotypes

In the white-throated sparrow (Zonotrichia albicollis), the second chromosome bears a striking resemblance to sex chromosomes. First, within each breeding pair of birds, one bird is homozygous for the standard arrangement of the chromosome (ZAL2/ZAL2) and its mate is heterozygous for a different version (ZAL2/ZAL2^m). Second, recombination is profoundly suppressed between the two versions, leading to genetic differentiation between them. Third, the ZAL2^m version is linked with phenotypic traits, such as bright plumage, high aggression, and low parental behavior, which are usually associated with males. These similarities to sex chromosomes suggest that the evolutionary mechanisms that shape sex chromosomes, in particular genetic degeneration of the heterogametic version due to the suppression of recombination, are likely important in this system as well. Here, we investigated patterns of protein sequence evolution and gene expression evolution between the ZAL2 and ZAL2^m chromosomes by whole-genome sequencing and transcriptome analyses. Patterns of protein evolution exhibited only weak signals of genetic degeneration, and few genes harbored signatures of positive selection. We found substantial evidence of transcriptome evolution, such as significant expression divergence between ZAL2 and ZAL2^m alleles and signatures of dosage compensation for highly expressed genes. These results suggest that, early in the evolution of heteromorphic chromosomes, gene expression divergence and dosage compensation can prevail before large-scale genetic degeneration. Our results show further that suppression of recombination between heteromorphic chromosomes can lead to the evolution of alternative (sex-like) behavioral phenotypes before substantial genetic degeneration.

A short story gets longer: recent insights into the molecular basis of heterostyly

Heterostyly is a fascinating adaptation to promote outbreeding and a classical paradigm of botany. In the most common type of heterostyly, plants either form flowers with long styles and short stamens, or short styles and long stamens. This reciprocal organ positioning reduces pollen wastage and promotes cross-pollination, thus increasing male fitness. In addition, in many heterostylous species selfing and the generation of unfit progeny due to inbreeding depression is limited by a self-incompatibility system, thus promoting female fitness. The two floral forms are genetically determined by the S locus as a complex supergene, namely a chromosomal region containing several individual genes that control the different traits, such as style or stamen length, and are held together by very tight linkage due to suppressed recombination. Recent molecular-genetic studies in several systems, including Turnera, Fagopyrum, Linum, and Primula have begun to identify and characterize the causal heterostyly genes residing at the S locus. An emerging theme from several families is that the dominant S haplotype represents a hemizygous region not present on the recessive s haplotype. This provides an explanation for the suppressed recombination and suggests a scenario for the chromosomal evolution of the S locus. In this review, we discuss the results from recent molecular-genetic analyses in light of the classical models on the genetics and evolution of heterostyly.

Social polymorphism is favoured by the co-evolution of dispersal with social behaviour

Dispersal determines gene flow among groups in a population and so plays a major role in many ecological and evolutionary processes. As gene flow shapes kin structure, dispersal is important to the evolution of social behaviours that influence reproduction within groups. Conversely, dispersal depends on kin structure and social behaviour. Dispersal and social behaviour therefore co-evolve, but the nature and consequences of this interplay are not well understood. Here, we show that it readily leads to the emergence of two social morphs: a sessile, benevolent morph expressed by individuals who tend to increase the reproduction of others within their group relative to their own; and a dispersive, self-serving morph expressed by individuals who tend to increase their own reproduction. This social polymorphism arises due to a positive linkage between the loci responsible for dispersal and social behaviour, leading to benevolent individuals preferentially interacting with relatives and self-serving individuals with non-relatives. We find that this linkage is favoured under a large spectrum of conditions, suggesting that associations between dispersal and other social traits should be common in nature. In line with this prediction, dispersers across a wide range of organisms have been reported to differ in their social tendencies from non-dispersers.

The genetic and molecular architecture of phenotypic diversity in sticklebacks

A major goal of evolutionary biology is to identify the genotypes and phenotypes that underlie adaptation to divergent environments. Stickleback fish, including the threespine stickleback (Gasterosteus aculeatus) and the ninespine stickleback (Pungitius pungitius), have been at the forefront of research to uncover the genetic and molecular architecture that underlies phenotypic diversity and adaptation. A wealth of quantitative trait locus (QTL) mapping studies in sticklebacks have provided insight into long-standing questions about the distribution of effect sizes during adaptation as well as the role of genetic linkage in facilitating adaptation. These QTL mapping studies have also provided a basis for the identification of the genes that underlie phenotypic diversity. These data have revealed that mutations in regulatory elements play an important role in the evolution of phenotypic diversity in sticklebacks. Genetic and molecular studies in sticklebacks have also led to new insights on the genetic basis of repeated evolution and suggest that the same loci are involved about half of the time when the same phenotypes evolve independently. When the same locus is involved, selection on standing variation and repeated mutation of the same genes have both contributed to the evolution of similar phenotypes in independent populations.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Complex History and Differentiation Patterns of the t-Haplotype, a Mouse Meiotic Driver

The t-haplotype, a mouse meiotic driver found on chromosome 17, has been a model for autosomal segregation distortion for close to a century, but several questions remain regarding its biology and evolutionary history. A recently published set of population genomics resources for wild mice includes several individuals heterozygous for the t-haplotype, which we use to characterize this selfish element at the genomic and transcriptomic level. Our results show that large sections of the t-haplotype have been replaced by standard homologous sequences, possibly due to occasional events of recombination, and that this complicates the inference of its history. As expected for a long genomic segment of very low recombination, the t-haplotype carries an excess of fixed nonsynonymous mutations compared to the standard chromosome. This excess is stronger for regions that have not undergone recent recombination, suggesting that occasional gene flow between the t and the standard chromosome may provide a mechanism to regenerate coding sequences that have accumulated deleterious mutations. Finally, we find that t-complex genes with altered expression largely overlap with deleted or amplified regions, and that carrying a t-haplotype alters the testis expression of genes outside of the t-complex, providing new leads into the pathways involved in the biology of this segregation distorter.

Tracing the origin and evolution of supergene mimicry in butterflies

Supergene mimicry is a striking phenomenon but we know little about the evolution of this trait in any species. Here, by studying genomes of butterflies from a recent radiation in which supergene mimicry has been isolated to the gene doublesex, we show that sexually dimorphic mimicry and female-limited polymorphism are evolutionarily related as a result of ancient balancing selection combined with independent origins of similar morphs in different lineages and secondary loss of polymorphism in other lineages. Evolutionary loss of polymorphism appears to have resulted from an interaction between natural selection and genetic drift. Furthermore, molecular evolution of the supergene is dominated not by adaptive protein evolution or balancing selection, but by extensive hitchhiking of linked variants on the mimetic dsx haplotype that occurred at the origin of mimicry. Our results suggest that chance events have played important and possibly opposing roles throughout the history of this classic example of adaptation.

Large-scale suppression of recombination predates genomic rearrangements in Neurospora tetrasperma

A common feature of eukaryote genomes is large chromosomal regions where recombination is absent or strongly reduced, but the factors that cause this reduction are not well understood. Genomic rearrangements have often been implicated, but they may also be a consequence of recombination suppression rather than a cause. In this study, we generate eight high-quality genomic data sets of the filamentous ascomycete Neurospora tetrasperma, a fungus that lacks recombination over most of its largest chromosome. The genomes surprisingly reveal collinearity of the non-recombining regions and although large inversions are enriched in these regions, we conclude these inversions to be derived and not the cause of the suppression. To our knowledge, this is the first time that non-recombining, genic regions as large as 86% of a full chromosome (or 8 Mbp), are shown to be collinear. These findings are of significant interest for our understanding of the evolution of sex chromosomes and other supergene complexes.

Chromosome Synapsis and Recombination in Male Hybrids between Two Chromosome Races of the Common Shrew (Sorex araneus L., Soricidae, Eulipotyphla)

Hybrid zones between chromosome races of the common shrew (Sorex araneus) provide exceptional models to study the potential role of chromosome rearrangements in the initial steps of speciation. The Novosibirsk and Tomsk races differ by a series of Robertsonian fusions with monobrachial homology. They form a narrow hybrid zone and generate hybrids with both simple (chain of three chromosomes) and complex (chain of eight or nine) synaptic configurations. Using immunolocalisation of the meiotic proteins, we examined chromosome pairing and recombination in males from the hybrid zone. Homozygotes and simple heterozygotes for Robertsonian fusions showed a low frequency of synaptic aberrations (<10%). The carriers of complex synaptic configurations showed multiple pairing abnormalities, which might lead to reduced fertility. The recombination frequency in the proximal regions of most chromosomes of all karyotypes was much lower than in the other regions. The strong suppression of recombination in the pericentromeric regions and co-segregation of race specific chromosomes involved in the long chains would be expected to lead to linkage disequilibrium between genes located there. Genic differentiation, together with the high frequency of pairing aberrations in male carriers of the long chains, might contribute to maintenance of the narrow hybrid zone.

Significant Synteny and Colocalization of Ecologically Relevant Quantitative Trait Loci Within and Across Species of Salmonid Fishes

The organization of functional regions within genomes has important implications for evolutionary potential. Considerable research effort has gone toward identifying the genomic basis of phenotypic traits of interest through quantitative trait loci (QTL) analyses. Less research has assessed the arrangement of QTL in the genome within and across species. To investigate the distribution, extent of colocalization, and the synteny of QTL for ecologically relevant traits, we used a comparative genomic mapping approach within and across a range of salmonid species. We compiled 943 QTL from all available species [lake whitefish (Coregonus clupeaformis), coho salmon (Oncorhynchus kisutch), rainbow trout (O. mykiss), Chinook salmon (O. tshawytscha), Atlantic salmon (Salmo salar), and Arctic charr (Salvelinus alpinus)]. We developed a novel analytical framework for mapping and testing the distribution of these QTL. We found no correlation between QTL density and gene density at the chromosome level but did at the fine-scale. Two chromosomes were significantly enriched for QTL. We found multiple synteny blocks for morphological, life history, and physiological traits across species, but only morphology and physiology had significantly more than expected. Two or three pairs of traits were significantly colocalized in three species (lake whitefish, coho salmon, and rainbow trout). Colocalization and fine-scale synteny suggest genetic linkage between traits within species and a conserved genetic basis across species. However, this pattern was weak overall, with colocalization and synteny being relatively rare. These findings advance our understanding of the role of genomic organization in the renowned ecological and phenotypic variability of salmonid fishes.

Candidate genes for shell colour polymorphism in Cepaea nemoralis

The characteristic ground colour and banding patterns on shells of the land snail Cepaea nemoralis form a classic study system for genetics and adaptation as it varies widely between individuals. We use RNAseq analysis to identify candidate genes underlying this polymorphism. We sequenced cDNA from the foot and the mantle (the shell-producing tissue) of four individuals of two phenotypes and produced a de novo transcriptome of 147,397 contigs. Differential expression analysis identified a set of 1,961 transcripts that were upregulated in mantle tissue. Sequence variant analysis resulted in a set of 2,592 transcripts with single nucleotide polymorphisms (SNPs) that differed consistently between the phenotypes. Inspection of the overlap between the differential expression analysis and SNP analysis yielded a set of 197 candidate transcripts, of which 38 were annotated. Four of these transcripts are thought to be involved in production of the shell’s nacreous layer. Comparison with morph-associated Restriction-site Associated DNA (RAD)-tags from a published study yielded eight transcripts that were annotated as metallothionein, a protein that is thought to inhibit the production of melanin in melanocytes. These results thus provide an excellent starting point for the elucidation of the genetic regulation of the Cepaea nemoralis shell colour polymorphism.

Fire ant social chromosomes: Differences in number, sequence and expression of odorant binding proteins

Variation in social behavior is common yet our knowledge of the mechanisms underpinning its evolution is limited. The fire ant Solenopsis invicta provides a textbook example of a Mendelian element controlling social organization: alternate alleles of a genetic element first identified as encoding an odorant binding protein (OBP) named Gp-9 determine whether a colony accepts one or multiple queens. The potential roles of such a protein in perceiving olfactory cues and evidence of positive selection on its amino acid sequence made it an appealing candidate gene. However, we recently showed that recombination is suppressed between Gp-9 and hundreds of other genes as part of a >19 Mb supergene-like region carried by a pair of social chromosomes. This finding raises the need to reassess the potential role of Gp-9. We identify 23 OBPs in the fire ant genome assembly, including nine located in the region of suppressed recombination with Gp-9. For six of these, the alleles carried by the two variants of the supergene-like region differ in protein-coding sequence and thus likely in function, with Gp-9 showing the strongest evidence of positive selection. We identify an additional OBP specific to the Sb variant of the region. Finally, we find that 14 OBPs are differentially expressed between single- and multiple-queen colonies. These results are consistent with multiple OBPs playing a role in determining social structure.

Predator Perspective Drives Geographic Variation in Frequency-Dependent Polymorphism

Color polymorphism in natural populations can manifest as a striking patchwork of phenotypes in space, with neighboring populations characterized by dramatic differences in morph composition. These geographic mosaics can be challenging to explain in the absence of localized selection because they are unlikely to result from simple isolation-by-distance or clinal variation in selective regimes. To identify processes that can lead to the formation of geographic mosaics, we developed a simulation-based model to explore the influence of predator perspective, selection, migration, and genetic linkage of color loci on allele frequencies in polymorphic populations over space and time. Using simulated populations inspired by the biology of Heliconius longwing butterflies, Cepaea land snails, Oophaga poison frogs, and Sonora ground snakes, we found that the relative sizes of predator and prey home ranges can produce large differences in morph composition between neighboring populations under both positive and negative frequency-dependent selection. We also demonstrated the importance of the interaction of predator perspective with the type of frequency dependence and localized directional selection across migration and selection intensities. Our results show that regional-scale predation can promote the formation of phenotypic mosaics in prey species, without the need to invoke spatial variation in selective regimes. We suggest that predator behavior can play an important and underappreciated role in the formation and maintenance of geographic mosaics in polymorphic species.

Genetics and Evolution of Social Behavior in Insects

The study of insect social behavior has offered tremendous insight into the molecular mechanisms mediating behavioral and phenotypic plasticity. Genomic applications to the study of eusocial insect species, in particular, have led to several hypotheses for the processes underlying the molecular evolution of behavior. Advances in understanding the genetic control of social organization have also been made, suggesting an important role for supergenes in the evolution of divergent behavioral phenotypes. Intensive study of social phenotypes across species has revealed that behavior and caste are controlled by an interaction between genetic and environmentally mediated effects and, further, that gene expression and regulation mediate plastic responses to environmental signals. However, several key methodological flaws that are hindering progress in the study of insect social behavior remain. After reviewing the current state of knowledge, we outline ongoing challenges in experimental design that remain to be overcome in order to advance the field.

The fire ant social chromosome supergene variant Sb shows low diversity but high divergence from SB

Variation in social behaviour is common, yet little is known about the genetic architectures underpinning its evolution. A rare exception is in the fire ant Solenopsis invicta: Alternative variants of a supergene region determine whether a colony will have exactly one or up to dozens of queens. The two variants of this region are carried by a pair of 'social chromosomes', SB and Sb, which resemble a pair of sex chromosomes. Recombination is suppressed between the two chromosomes in the supergene region. While the X-like SB can recombine with itself in SB/SB queens, recombination is effectively absent in the Y-like Sb because Sb/Sb queens die before reproducing. Here, we analyse whole-genome sequences of eight haploid SB males and eight haploid Sb males. We find extensive SB-Sb differentiation throughout the >19-Mb-long supergene region. We find no evidence of 'evolutionary strata' with different levels of divergence comparable to those reported in several sex chromosomes. A high proportion of substitutions between the SB and Sb haplotypes are nonsynonymous, suggesting inefficacy of purifying selection in Sb sequences, similar to that for Y-linked sequences in XY systems. Finally, we show that the Sb haplotype of the supergene region has 635-fold less nucleotide diversity than the rest of the genome. We discuss how this reduction could be due to a recent selective sweep affecting Sb specifically or associated with a population bottleneck during the invasion of North America by the sampled population.

Highly localized divergence within supergenes in Atlantic cod (Gadus morhua) within the Gulf of Maine

BACKGROUND

Atlantic cod (Gadus morhua), is known to vary genetically across the North Atlantic, Greenland, and Newfoundland. This genetic variation occurs both spatially and temporally through decades of heavy fishing, and is concentrated in three linkage disequilibrium blocks, previously defined by pedigreed linkage mapping analysis. Variation within these genomic regions is correlated with both seawater temperature and behavioral ecotype. The full extent and nature of these linkage groups is important information for interpreting cod genetic structure as a tool for future fisheries management.

RESULTS

We conducted whole genome sequencing for 31 individual cod from three sub-populations in the Gulf of Maine. Across the genome, we found 3,390,654 intermediate to high frequency Single Nucleotide Polymorphisms (SNPs). We show that pairwise linkage analysis among these SNPs is a powerful tool to detect linkage disequilibrium clusters by recovering the three previously detected linkage groups and identifying the 1031 genes contained therein. Across these genes, we found significant population differentiation among spawning groups in the Gulf of Maine and between Georges Bank and Gulf of Maine. Coordinated divergence among these genes and their differentiation at both short and long spatial scales suggests that they are acting as linked supergenes in local adaptation of cod populations.

CONCLUSIONS

Differentiation between SNPs in linkage disequilibrium blocks is the major signal of genetic differentiation between all groups tested within the Gulf of Maine. Our data provide a map of genes contained in these blocks, allowing an enhanced search for neutral genetic structure for demographic inference and fisheries modeling. Patterns of selection and the history of populations may be possible to identify in cod using this description of linkage disequilibrium blocks and future data sets to robustly separate neutral and selected genetic markers.

Coat of many colours-DNA reveals polymorphism of mantle patterns and colouration in Caribbean Cyphoma Röding, 1798 (Gastropoda, Ovulidae)

The iconic gastropod genus Cyphoma is commonly observed in the Caribbean, where it lives in association with various octocorallian hosts. Each species in the genus Cyphoma has a unique, characteristic mantle pattern and colouration, which separates the valid taxa. Because of its abundance and recognisability Cyphoma gibbosum has been used as a model organism in several studies concerning allelochemicals, reef degradation, and physical defence mechanisms. Molecular analyses based on four molecular markers (COI, 16S, H3 and 28S) for three Cyphoma species (C. gibbosum, C. mcgintyi, C. signatum) and an unidentified black morph, collected from three localities in the Caribbean, show that they represent morphological varieties of a single, genetically homogeneous species. This outcome is in agreement with previous anatomical studies. As a result C. mcgintyi and C. signatum are synonymised with C. gibbosum, which is a key result for future work using C. gibbosum as a model organism. The striking morphological differences in mantle pattern and colouration are hypothesised to be the result of one of three possible scenarios: rapid divergence, supergenes (including balanced polymorphism), or incipient speciation.

Genetic architecture and balancing selection: the life and death of differentiated variants

Balancing selection describes any form of natural selection, which results in the persistence of multiple variants of a trait at intermediate frequencies within populations. By offering up a snapshot of multiple co-occurring functional variants and their interactions, systems under balancing selection can reveal the evolutionary mechanisms favouring the emergence and persistence of adaptive variation in natural populations. We here focus on the mechanisms by which several functional variants for a given trait can arise, a process typically requiring multiple epistatic mutations. We highlight how balancing selection can favour specific features in the genetic architecture and review the evolutionary and molecular mechanisms shaping this architecture. First, balancing selection affects the number of loci underlying differentiated traits and their respective effects. Control by one or few loci favours the persistence of differentiated functional variants by limiting intergenic recombination, or its impact, and may sometimes lead to the evolution of supergenes. Chromosomal rearrangements, particularly inversions, preventing adaptive combinations from being dissociated are increasingly being noted as features of such systems. Similarly, due to the frequency of heterozygotes maintained by balancing selection, dominance may be a key property of adaptive variants. High heterozygosity and limited recombination also influence associated genetic load, as linked recessive deleterious mutations may be sheltered. The capture of deleterious elements in a locus under balancing selection may reinforce polymorphism by further promoting heterozygotes. Finally, according to recent genomewide scans, balanced polymorphism might be more pervasive than generally thought. We stress the need for both functional and ecological studies to characterize the evolutionary mechanisms operating in these systems.

Integrative Approaches for Studying Mitochondrial and Nuclear Genome Co-evolution in Oxidative Phosphorylation

In animals, interactions among gene products of mitochondrial and nuclear genomes (mitonuclear interactions) are of profound fitness, evolutionary, and ecological significance. Most fundamentally, the oxidative phosphorylation (OXPHOS) complexes responsible for cellular bioenergetics are formed by the direct interactions of 13 mitochondrial-encoded and ∼80 nuclear-encoded protein subunits in most animals. It is expected that organisms will develop genomic architecture that facilitates co-adaptation of these mitonuclear interactions and enhances biochemical efficiency of OXPHOS complexes. In this perspective, we present principles and approaches to understanding the co-evolution of these interactions, with a novel focus on how genomic architecture might facilitate it. We advocate that recent interdisciplinary advances assist in the consolidation of links between genotype and phenotype. For example, advances in genomics allow us to unravel signatures of selection in mitochondrial and nuclear OXPHOS genes at population-relevant scales, while newly published complete atomic-resolution structures of the OXPHOS machinery enable more robust predictions of how these genes interact epistatically and co-evolutionarily. We use three case studies to show how integrative approaches have improved the understanding of mitonuclear interactions in OXPHOS, namely those driving high-altitude adaptation in bar-headed geese, allopatric population divergence in Tigriopus californicus copepods, and the genome architecture of nuclear genes coding for mitochondrial functions in the eastern yellow robin.