Supergenes and complex phenotypes

Ecology shapes epistasis in a genotype-phenotype-fitness map for stick insect colour

Genetic interactions such as epistasis are widespread in nature and can shape evolutionary dynamics. Epistasis occurs due to nonlinearity in biological systems, which can arise via cellular processes that convert genotype to phenotype and via selective processes that connect phenotype to fitness. Few studies in nature have connected genotype to phenotype to fitness for multiple potentially interacting genetic variants. Thus, the causes of epistasis in the wild remain poorly understood. Here, we show that epistasis for fitness is an emergent and predictable property of nonlinear selective processes. We do so by measuring the genetic basis of cryptic colouration and survival in a field experiment with stick insects. We find that colouration shows a largely additive genetic basis but with some effects of epistasis that enhance differentiation between colour morphs. In terms of fitness, different combinations of loci affecting colouration confer high survival in one host-plant treatment. Specifically, nonlinear correlational selection for specific combinations of colour traits in this treatment drives the emergence of pairwise and higher-order epistasis for fitness at loci underlying colour. In turn, this results in a rugged fitness landscape for genotypes. In contrast, fitness epistasis was dampened in another treatment, where selection was weaker. Patterns of epistasis that are shaped by ecologically based selection could be common and central to understanding fitness landscapes, the dynamics of evolution and potentially other complex systems.

Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal, poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab-adapted ecotypes that exhibit partial reproductive isolation from wave-adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co-association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes-consistent with models of adaptation with gene flow-but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co-association of loci during ecological speciation with ongoing gene flow.

Inside the supergene of the bird with four sexes

The white-throated sparrow (Zonotrichia albicollis) offers unique opportunities to understand the adaptive value of supergenes, particularly their role in alternative phenotypes. In this species, alternative plumage morphs segregate with a nonrecombining segment of chromosome 2, which has been called a 'supergene'. The species mates disassortatively with respect to the supergene; that is, each breeding pair consists of one individual with it and one without it. This species has therefore been called the "bird with four sexes". The supergene segregates with a behavioral phenotype; birds with it are more aggressive and less parental than birds without it. Here, we review our efforts to identify the genes inside the supergene that are responsible for the behavioral polymorphism. The gene ESR1, which encodes estrogen receptor α, differs between the morphs and predicts both territorial and parental behavior. Variation in the regulatory regions of ESR1 causes an imbalance in expression of the two alleles, and the degree to which this imbalance favors the supergene allele predicts territorial singing. In heterozygotes, knockdown of ESR1 causes a phenotypic switch, from more aggressive to less aggressive. We recently showed that another gene important for social behavior, vasoactive intestinal peptide (VIP), is differentially expressed between the morphs and predicts territorial singing. We hypothesize that ESR1 and VIP contribute to behavior in a coordinated way and could represent co-adapted alleles. Because the supergene contains more than 1000 individual genes, this species provides rich possibilities for discovering alleles that work together to mediate life-history trade-offs and maximize the fitness of alternative complex phenotypes.

Identification of transcription factors controlling floral morphology in wild Petunia species with contrasting pollination syndromes

Adaptation to different pollinators is an important driver of speciation in the angiosperms. Genetic approaches such as QTL mapping have been successfully used to identify the underlying speciation genes. However, these methods are often limited by widespread suppression of recombination due to divergence between species. While the mutations that caused the interspecific differences in floral color and scent have been elucidated in a variety of plant genera, the genes that are responsible for morphological differences remain mostly unknown. Differences in floral organ length determine the pollination efficiency of hawkmoths and hummingbirds, and therefore the genes that control these differences are potential speciation genes. Identifying such genes is challenging, especially in non-model species and when studying complex traits for which little prior genetic and biochemical knowledge is available. Here we combine transcriptomics with detailed growth analysis to identify candidate transcription factors underlying interspecific variation in the styles of Petunia flowers. Starting from a set of 2284 genes, stepwise filtering for expression in styles, differential expression between species, correlation with growth-related traits, allele-specific expression in interspecific hybrids, and/or high-impact polymorphisms resulted in a set of 43 candidate speciation genes. Validation by virus-induced gene silencing identified two MYB transcription factors, EOBI and EOBII, that were previously shown to regulate floral scent emission, a trait associated with pollination by hawkmoths.

Supergene evolution via stepwise duplications and neofunctionalization of a floral-organ identity gene

Heterostyly represents a fascinating adaptation to promote outbreeding in plants that evolved multiple times independently. While l-morph individuals form flowers with long styles, short anthers, and small pollen grains, S-morph individuals have flowers with short styles, long anthers, and large pollen grains. The difference between the morphs is controlled by an S-locus "supergene" consisting of several distinct genes that determine different traits of the syndrome and are held together, because recombination between them is suppressed. In Primula, the S locus is a roughly 300-kb hemizygous region containing five predicted genes. However, with one exception, their roles remain unclear, as does the evolutionary buildup of the S locus. Here we demonstrate that the MADS-box GLOBOSA2 (GLO2) gene at the S locus determines anther position. In Primula forbesii S-morph plants, GLO2 promotes growth by cell expansion in the fused tube of petals and stamen filaments beneath the anther insertion point; by contrast, neither pollen size nor male incompatibility is affected by GLO2 activity. The paralogue GLO1, from which GLO2 arose by duplication, has maintained the ancestral B-class function in specifying petal and stamen identity, indicating that GLO2 underwent neofunctionalization, likely at the level of the encoded protein. Genetic mapping and phylogenetic analysis indicate that the duplications giving rise to the style-length-determining gene CYP734A50 and to GLO2 occurred sequentially, with the CYP734A50 duplication likely the first. Together these results provide the most detailed insight into the assembly of a plant supergene yet and have important implications for the evolution of heterostyly.

Host-parasite co-evolution and its genomic signature

Studies in diverse biological systems have indicated that host-parasite co-evolution is responsible for the extraordinary genetic diversity seen in some genomic regions, such as major histocompatibility (MHC) genes in jawed vertebrates and resistance genes in plants. This diversity is believed to evolve under balancing selection on hosts by parasites. However, the mechanisms that link the genomic signatures in these regions to the underlying co-evolutionary process are only slowly emerging. We still lack a clear picture of the co-evolutionary concepts and of the genetic basis of the co-evolving phenotypic traits in the interacting antagonists. Emerging genomic tools that provide new options for identifying underlying genes will contribute to a fuller understanding of the co-evolutionary process.

Simple inheritance, complex regulation: Supergene-mediated fire ant queen polymorphism

The fire ant Solenopsis invicta exists in two alternate social forms: monogyne nests contain a single reproductive queen and polygyne nests contain multiple reproductive queens. This colony-level social polymorphism corresponds with individual differences in queen physiology, queen dispersal patterns and worker discrimination behaviours, all evidently regulated by an inversion-based supergene that spans more than 13 Mb of a "social chromosome," contains over 400 protein-coding genes and rarely undergoes recombination. The specific mechanisms by which this supergene influences expression of the many distinctive features that characterize the alternate forms remain almost wholly unknown. To advance our understanding of these mechanisms, we explore the effects of social chromosome genotype and natal colony social form on gene expression in queens sampled as they embarked on nuptial flights, using RNA-sequencing of brains and ovaries. We observe a large effect of natal social form, that is, of the social/developmental environment, on gene expression profiles, with similarly substantial effects of genotype, including: (a) supergene-associated gene upregulation, (b) allele-specific expression and (c) pronounced extra-supergene trans-regulatory effects. These findings, along with observed spatial variation in differential and allele-specific expression within the supergene region, highlight the complex gene regulatory landscape that emerged following divergence of the inversion-mediated Sb haplotype from its homologue, which presumably largely retained the ancestral gene order. The distinctive supergene-associated gene expression trajectories we document at the onset of a queen's reproductive life expand the known record of relevant molecular correlates of a complex social polymorphism and point to putative genetic factors underpinning the alternate social syndromes.

Mimicry genes reduce pre-adult survival rate in Papilio polytes: A possible new mechanism for maintaining female-limited polymorphism in Batesian mimicry

Batesian mimicry, in which harmless organisms resemble unpalatable or harmful species, is a well-studied adaptation for predation avoidance. The females of some Batesian mimic species comprise mimetic and nonmimetic individuals. Mimetic females of such polymorphic species clearly have a selective advantage due to decreased predation pressure, but the selective forces that maintain nonmimetic females in a population remain unclear. In the swallowtail butterfly, Papilio polytes, female polymorphism is controlled by the H (mimetic) and h (nonmimetic) alleles at a single autosomal locus. Here, we examined whether the dominant H allele has a deleterious effect on the pre-adult survival rate (egg-to-adult emergence rate). We repeated an assortative mating-like treatment-that is breeding of males and females whose mothers had the same phenotype (mimetic or nonmimetic)-for three consecutive generations, while avoiding inbreeding. Results showed that pre-adult survival rate decreased over generations only in lines derived from mothers with the mimetic phenotype (hereafter, mimetic-assorted lines). This lowered survival was due to an increased mortality at the final instar larval stage and the pupal stages. Interestingly, the pre-adult mortality in the mimetic-assorted lines seemed to be associated with a male-biased sex ratio at adult emergence. These results suggest that the dominant H allele displays a mildly deleterious effect that is expressed more strongly in females and homozygous individuals than in heterozygous individuals. We propose that this cost of mimicry in larval and pupal stages contributes to the maintenance of female-limited polymorphism in P. polytes.

Consequences of Single-Locus and Tightly Linked Genomic Architectures for Evolutionary Responses to Environmental Change

Genetic and genomic architectures of traits under selection are key factors influencing evolutionary responses. Yet, knowledge of their impacts has been limited by a widespread assumption that most traits are controlled by unlinked polygenic architectures. Recent advances in genome sequencing and eco-evolutionary modeling are unlocking the potential for integrating genomic information into predictions of population responses to environmental change. Using eco-evolutionary simulations, we demonstrate that hypothetical single-locus control of a life history trait produces highly variable and unpredictable harvesting-induced evolution relative to the classically applied multilocus model. Single-locus control of complex traits is thought to be uncommon, yet blocks of linked genes, such as those associated with some types of structural genomic variation, have emerged as taxonomically widespread phenomena. Inheritance of linked architectures resembles that of single loci, thus enabling single-locus-like modeling of polygenic adaptation. Yet, the number of loci, their effect sizes, and the degree of linkage among them all occur along a continuum. We review how linked architectures are often associated, directly or indirectly, with traits expected to be under selection from anthropogenic stressors and are likely to play a large role in adaptation to environmental disturbance. We suggest using single-locus models to explore evolutionary extremes and uncertainties when the trait architecture is unknown, refining parameters as genomic information becomes available, and explicitly incorporating linkage among loci when possible. By overestimating the complexity (e.g., number of independent loci) of the genomic architecture of traits under selection, we risk underestimating the complexity (e.g., nonlinearity) of their evolutionary dynamics.

Structural genomic variation leads to genetic differentiation in Lake Tanganyika's sardines

Identifying patterns in genetic structure and the genetic basis of ecological adaptation is a core goal of evolutionary biology and can inform the management and conservation of species that are vulnerable to population declines exacerbated by climate change. We used reduced-representation genomic sequencing methods to gain a better understanding of genetic structure among and within populations of Lake Tanganyika's two sardine species, Limnothrissa miodon and Stolothrissa tanganicae. Samples of these ecologically and economically important species were collected across the length of Lake Tanganyika, as well as from nearby Lake Kivu, where L. miodon was introduced in 1959. Our results reveal differentiation within both S. tanganicae and L. miodon that is not explained by geography. Instead, this genetic differentiation is due to the presence of large sex-specific regions in the genomes of both species, but involving different polymorphic sites in each species. Our results therefore indicate rapidly evolving XY sex determination in the two species. Additionally, we found evidence of a large chromosomal rearrangement in L. miodon, creating two homokaryotypes and one heterokaryotype. We found all karyotypes throughout Lake Tanganyika, but the frequencies vary along a north-south gradient and differ substantially in the introduced Lake Kivu population. We do not find evidence for significant isolation by distance, even over the hundreds of kilometres covered by our sampling, but we do find shallow population structure.

Large-scale mutation in the evolution of a gene complex for cryptic coloration

The types of mutations affecting adaptation in the wild are only beginning to be understood. In particular, whether structural changes shape adaptation by suppressing recombination or by creating new mutations is unresolved. Here, we show that multiple linked but recombining loci underlie cryptic color morphs of Timema chumash stick insects. In a related species, these loci are found in a region of suppressed recombination, forming a supergene. However, in seven species of Timema, we found that a megabase-size "supermutation" has deleted color loci in green morphs. Moreover, we found that balancing selection likely contributes more to maintaining this mutation than does introgression. Our results show how suppressed recombination and large-scale mutation can help to package gene complexes into discrete units of diversity such as morphs, ecotypes, or species.

Adaptation via pleiotropy and linkage: Association mapping reveals a complex genetic architecture within the stickleback Eda locus

Genomic mapping of the loci associated with phenotypic evolution has revealed genomic "hotspots," or regions of the genome that control multiple phenotypic traits. This clustering of loci has important implications for the speed and maintenance of adaptation and could be due to pleiotropic effects of a single mutation or tight genetic linkage of multiple causative mutations affecting different traits. The threespine stickleback (Gasterosteus aculeatus) is a powerful model for the study of adaptive evolution because the marine ecotype has repeatedly adapted to freshwater environments across the northern hemisphere in the last 12,000 years. Freshwater ecotypes have repeatedly fixed a 16 kilobase haplotype on chromosome IV that contains Ectodysplasin (Eda), a gene known to affect multiple traits, including defensive armor plates, lateral line sensory hair cells, and schooling behavior. Many additional traits have previously been mapped to a larger region of chromosome IV that encompasses the Eda freshwater haplotype. To identify which of these traits specifically map to this adaptive haplotype, we made crosses of rare marine fish heterozygous for the freshwater haplotype in an otherwise marine genetic background. Further, we performed fine-scale association mapping in a fully interbreeding, polymorphic population of freshwater stickleback to disentangle the effects of pleiotropy and linkage on the phenotypes affected by this haplotype. Although we find evidence that linked mutations have small effects on a few phenotypes, a small 1.4-kb region within the first intron of Eda has large effects on three phenotypic traits: lateral plate count, and both the number and patterning of the posterior lateral line neuromasts. Thus, the Eda haplotype is a hotspot of adaptation in stickleback due to both a small, pleiotropic region affecting multiple traits as well as multiple linked mutations affecting additional traits.

Maternal effect killing by a supergene controlling ant social organization

Supergenes underlie striking polymorphisms in nature, yet the evolutionary mechanisms by which they arise and persist remain enigmatic. These clusters of linked loci can spread in populations because they captured coadapted alleles or by selfishly distorting the laws of Mendelian inheritance. Here, we show that the supergene haplotype associated with multiple-queen colonies in Alpine silver ants is a maternal effect killer. All eggs from heterozygous queens failed to hatch when they did not inherit this haplotype. Hence, the haplotype specific to multiple-queen colonies is a selfish genetic element that enhances its own transmission by causing developmental arrest of progeny that do not carry it. At the population level, such transmission ratio distortion favors the spread of multiple-queen colonies, to the detriment of the alternative haplotype associated with single-queen colonies. Hence, selfish gene drive by one haplotype will impact the evolutionary dynamics of alternative forms of colony social organization. This killer hidden in a social supergene shows that large nonrecombining genomic regions are prone to cause multifarious effects across levels of biological organization.

QTL for horticulturally important traits associated with pleiotropic andromonoecy and carpel number loci, and a paracentric inversion in cucumber

The legendary cucumber inbred line WI2757 possesses a rare combination of resistances against nine pathogens, which is an important germplasm for cucumber breeding. However, WI2757 flowers late and does not perform well under field conditions. The genetic basis for horticulturally important traits other than disease resistances in WI2757 is largely unknown. In this study, we conducted QTL mapping using F₂ and recombinant inbred line (RIL) populations from the WI2757 × True Lemon cross that were segregating for multiple traits. Phenotypic data were collected in replicated field trials across multiple years for seven traits including fruit carpel number (CN) and sex expression. A high-density SNP-based genetic map was developed with genotyping by sequencing of the RIL population, which revealed a region on chromosome 1 with strong recombination suppression. The reduced recombination in this region was due to a ~ 10-Mbp paracentric inversion in WI2757 that was confirmed with additional segregation and cytological (FISH) analyses. Thirty-six QTL were detected for flowering time, fruit length (FL), fruit diameter (FD), fruit shape (LD), fruit number (FN), CN, and powdery mildew resistance. Five moderate- or major-effect QTL for FL, FD, LD, and FN inside the inversion are likely the pleiotropic effects of the andromonoecy (m), or the cn locus. The major-effect flowering time QTL ft1.1 was also mapped inside the inversion, which seems to be different from the previously assigned delayed flowering in WI2757. Implications of these findings on the use of WI2757 in cucumber breeding are discussed.

An alternative route of bacterial infection associated with a novel resistance locus in the Daphnia-Pasteuria host-parasite system

To understand the mechanisms of antagonistic coevolution, it is crucial to identify the genetics of parasite resistance. In the Daphnia magna–Pasteuria ramosa host–parasite system, the most important step of the infection process is the one in which P. ramosa spores attach to the host’s foregut. A matching-allele model (MAM) describes the host–parasite genetic interactions underlying attachment success. Here we describe a new P. ramosa genotype, P15, which, unlike previously studied genotypes, attaches to the host’s hindgut, not to its foregut. Host resistance to P15 attachment shows great diversity across natural populations. In contrast to P. ramosa genotypes that use foregut attachment, P15 shows some quantitative variation in attachment success and does not always lead to successful infections, suggesting that hindgut attachment represents a less-efficient infection mechanism than foregut attachment. Using a Quantitative Trait Locus (QTL) approach, we detect two significant QTLs in the host genome: one that co-localizes with the previously described D. magna PR locus of resistance to foregut attachment, and a second, major QTL located in an unlinked genomic region. We find no evidence of epistasis. Fine mapping reveals a genomic region, the D locus, of ~13 kb. The discovery of a second P. ramosa attachment site and of a novel host-resistance locus increases the complexity of this system, with implications for both for the coevolutionary dynamics (e.g., Red Queen and the role of recombination), and for the evolution and epidemiology of the infection process.

The genomics of coloration provides insights into adaptive evolution

Coloration is an easily quantifiable visual trait that has proven to be a highly tractable system for genetic analysis and for studying adaptive evolution. The application of genomic approaches to evolutionary studies of coloration is providing new insight into the genetic architectures underlying colour traits, including the importance of large-effect mutations and supergenes, the role of development in shaping genetic variation and the origins of adaptive variation, which often involves adaptive introgression. Improved knowledge of the genetic basis of traits can facilitate field studies of natural selection and sexual selection, making it possible for strong selection and its influence on the genome to be demonstrated in wild populations.

Mechanisms for Color Convergence in a Mimetic Radiation of Poison Frogs

In animals, bright colors often evolve to mimic other species when a resemblance is selectively favored. Understanding the proximate mechanisms underlying such color mimicry can give insights into how mimicry evolves-for example, whether color convergence evolves from a shared set of mechanisms or through the evolution of novel color production mechanisms. We studied color production mechanisms in poison frogs (Dendrobatidae), focusing on the mimicry complex of Ranitomeya imitator. Using reflectance spectrometry, skin pigment analysis, electron microscopy, and color modeling, we found that the bright colors of these frogs, both within and outside the mimicry complex, are largely structural and produced by iridophores but that color production depends crucially on interactions with pigments. Color variation and mimicry are regulated predominantly by iridophore platelet thickness and, to a lesser extent, concentration of the red pteridine pigment drosopterin. Compared with each of the four morphs of model species that it resembles, R. imitator displays greater variation in both structural and pigmentary mechanisms, which may have facilitated phenotypic divergence in this species. Analyses of nonmimetic dendrobatids in other genera demonstrate that these mechanisms are widespread within the family and that poison frogs share a complex physiological "color palette" that can produce diverse and highly reflective colors.

Retracted: The evolutionary history of colour polymorphism in Ischnura damselflies

The above article from Journal of Evolutionary Biology, published online on 24 May 2018 in Wiley Online Library (http://wileyonlinelibrary.com), has been retracted on the request of the authors and with the agreement of the Journal's Editor in Chief Wolf Blanckenhorn and John Wiley & Sons, following disagreement on potential corrections to the article after publication. The decision to retract followed significant issues with the methods and analyses of the manuscript that were originally not uncovered during peer-review, but which were subsequently brought to the Journal's attention following publication of the Article on Early View. [Correction added on 2 July 2021, after first online publication: retraction statement has been modified.].

Butterfly Mimicry Polymorphisms Highlight Phylogenetic Limits of Gene Reuse in the Evolution of Diverse Adaptations

Some genes have repeatedly been found to control diverse adaptations in a wide variety of organisms. Such gene reuse reveals not only the diversity of phenotypes these unique genes control but also the composition of developmental gene networks and the genetic routes available to and taken by organisms during adaptation. However, the causes of gene reuse remain unclear. A small number of large-effect Mendelian loci control a huge diversity of mimetic butterfly wing color patterns, but reasons for their reuse are difficult to identify because the genetic basis of mimicry has primarily been studied in two systems with correlated factors: female-limited Batesian mimicry in Papilio swallowtails (Papilionidae) and non-sex-limited Müllerian mimicry in Heliconius longwings (Nymphalidae). Here, we break the correlation between phylogenetic relationship and sex-limited mimicry by identifying loci controlling female-limited mimicry polymorphism Hypolimnas misippus (Nymphalidae) and non-sex-limited mimicry polymorphism in Papilio clytia (Papilionidae). The Papilio clytia polymorphism is controlled by the genome region containing the gene cortex, the classic P supergene in Heliconius numata, and loci controlling color pattern variation across Lepidoptera. In contrast, female-limited mimicry polymorphism in Hypolimnas misippus is associated with a locus not previously implicated in color patterning. Thus, although many species repeatedly converged on cortex and its neighboring genes over 120 My of evolution of diverse color patterns, female-limited mimicry polymorphisms each evolved using a different gene. Our results support conclusions that gene reuse occurs mainly within ∼10 My and highlight the puzzling diversity of genes controlling seemingly complex female-limited mimicry polymorphisms.

Negative frequency dependent selection contributes to the maintenance of a global polymorphism in mitochondrial DNA

BACKGROUND

Understanding the forces that maintain diversity across a range of scales is at the very heart of biology. Frequency-dependent processes are generally recognized as the most central process for the maintenance of ecological diversity. The same is, however, not generally true for genetic diversity. Negative frequency dependent selection, where rare genotypes have an advantage, is often regarded as a relatively weak force in maintaining genetic variation in life history traits because recombination disassociates alleles across many genes. Yet, many regions of the genome show low rates of recombination and genetic variation in such regions (i.e., supergenes) may in theory be upheld by frequency dependent selection.

RESULTS

We studied what is essentially a ubiquitous life history supergene (i.e., mitochondrial DNA) in the fruit fly Drosophila subobscura, showing sympatric polymorphism with two main mtDNA genotypes co-occurring in populations world-wide. Using an experimental evolution approach involving manipulations of genotype starting frequencies, we show that negative frequency dependent selection indeed acts to maintain genetic variation in this region. Moreover, the strength of selection was affected by food resource conditions.

CONCLUSIONS

Our work provides novel experimental support for the view that balancing selection through negative frequency dependency acts to maintain genetic variation in life history genes. We suggest that the emergence of negative frequency dependent selection on mtDNA is symptomatic of the fundamental link between ecological processes related to resource use and the maintenance of genetic variation.

Balancing selection via life-history trade-offs maintains an inversion polymorphism in a seaweed fly

How natural diversity is maintained is an evolutionary puzzle. Genetic variation can be eroded by drift and directional selection but some polymorphisms persist for long time periods, implicating a role for balancing selection. Here, we investigate the maintenance of a chromosomal inversion polymorphism in the seaweed fly Coelopa frigida. Using experimental evolution and quantifying fitness, we show that the inversion underlies a life-history trade-off, whereby each haplotype has opposing effects on larval survival and adult reproduction. Numerical simulations confirm that such antagonistic pleiotropy can maintain polymorphism. Our results also highlight the importance of sex-specific effects, dominance and environmental heterogeneity, whose interaction enhances the maintenance of polymorphism through antagonistic pleiotropy. Overall, our findings directly demonstrate how overdominance and sexual antagonism can emerge from a life-history trade-off, inviting reconsideration of antagonistic pleiotropy as a key part of multi-headed balancing selection processes that enable the persistence of genetic variation.

Few studies empirically pinpoint how balanced polymorphisms are maintained. “Mérot et al”. identify an inversion polymorphism that is maintained in seaweed fly populations because of antagonistic pleiotropy that mediates a classic life history tradeoff between larval survival and adult reproduction.

Evolution of a supergene that regulates a trans-species social polymorphism

Supergenes are clusters of linked genetic loci that jointly affect the expression of complex phenotypes, such as social organization. Little is known about the origin and evolution of these intriguing genomic elements. Here we analyse whole-genome sequences of males from native populations of six fire ant species and show that variation in social organization is under the control of a novel supergene haplotype (termed Sb), which evolved by sequential incorporation of three inversions spanning half of a 'social chromosome'. Two of the inversions interrupt protein-coding genes, resulting in the increased expression of one gene and modest truncation in the primary protein structure of another. All six socially polymorphic species studied harbour the same three inversions, with the single origin of the supergene in their common ancestor inferred by phylogenomic analyses to have occurred half a million years ago. The persistence of Sb along with the ancestral SB haplotype through multiple speciation events provides a striking example of a functionally important trans-species social polymorphism presumably maintained by balancing selection. We found that while recombination between the Sb and SB haplotypes is severely restricted in all species, a low level of gene flux between the haplotypes has occurred following the appearance of the inversions, potentially mitigating the evolutionary degeneration expected at genomic regions that cannot freely recombine. These results provide a detailed picture of the structural genomic innovations involved in the formation of a supergene controlling a complex social phenotype.

An Ancient and Eroded Social Supergene Is Widespread across Formica Ants

Supergenes, clusters of tightly linked genes, play a key role in the evolution of complex adaptive variation [1, 2]. Although supergenes have been identified in many species, we lack an understanding of their origin, evolution, and persistence [3]. Here, we uncover 20-40 Ma of evolutionary history of a supergene associated with polymorphic social organization in Formica ants [4]. We show that five Formica species exhibit homologous divergent haplotypes spanning 11 Mbp on chromosome 3. Despite the supergene's size, only 142 single nucleotide polymorphisms (SNPs) consistently distinguish alternative supergene haplotypes across all five species. These conserved trans-species SNPs are localized in a small number of disjunct clusters distributed across the supergene. This unexpected pattern of divergence indicates that the Formica supergene does not follow standard models of sex chromosome evolution, in which distinct evolutionary strata reflect an expanding region of suppressed recombination [5]. We propose an alternative "eroded strata model" in which clusters of conserved trans-species SNPs represent functionally important areas maintained by selection in the face of rare recombination between ancestral haplotypes. The comparison of whole-genome sequences across 10 additional Formica species reveals that the most conserved region of the supergene contains a transcription factor essential for motor neuron development in Drosophila [6]. The discovery that a very small portion of this large and ancient supergene harbors conserved trans-species SNPs linked to colony social organization suggests that the ancestral haplotypes have been eroded by recombination, with selection preserving differentiation at one or a few genes generating alternative social organization.

Genome of an iconic Australian bird: High-quality assembly and linkage map of the superb fairy-wren (Malurus cyaneus)

The superb fairy-wren, Malurus cyaneus, is one of the most iconic Australian passerine species. This species belongs to an endemic Australasian clade, Meliphagides, which diversified early in the evolution of the oscine passerines. Today, the oscine passerines comprise almost half of all avian species diversity. Despite the rapid increase of available bird genome assemblies, this part of the avian tree has not yet been represented by a high-quality reference. To rectify that, we present the first high-quality genome assembly of a Meliphagides representative: the superb fairy-wren. We combined Illumina shotgun and mate-pair sequences, PacBio long-reads, and a genetic linkage map from an intensively sampled pedigree of a wild population to generate this genome assembly. Of the final assembled 1.07-Gb genome, 975 Mb (90.4%) was anchored onto 25 pseudochromosomes resulting in a final superscaffold N50 of 68.11 Mb. This high-quality bird genome assembly is one of only a handful which is also accompanied by a genetic map and recombination landscape. In comparison to other pedigree-based bird genetic maps, we find that the fairy-wren genetic map more closely resembles those of Taeniopygia guttata and Parus major maps, unlike the Ficedula albicollis map which more closely resembles that of Gallus gallus. Lastly, we also provide a predictive gene and repeat annotation of the genome assembly. This new high-quality, annotated genome assembly will be an invaluable resource not only regarding the superb fairy-wren species and relatives but also broadly across the avian tree by providing a novel reference point for comparative genomic analyses.

A transposable element insertion is associated with an alternative life history strategy

Tradeoffs affect resource allocation during development and result in fitness consequences that drive the evolution of life history strategies. Yet despite their importance, we know little about the mechanisms underlying life history tradeoffs. Many species of Colias butterflies exhibit an alternative life history strategy (ALHS) where females divert resources from wing pigment synthesis to reproductive and somatic development. Due to this reallocation, a wing color polymorphism is associated with the ALHS: either yellow/orange or white. Here we map the locus associated with this ALHS in Colias crocea to a transposable element insertion located downstream of the Colias homolog of BarH-1, a homeobox transcription factor. Using CRISPR/Cas9 gene editing, antibody staining, and electron microscopy we find white-specific expression of BarH-1 suppresses the formation of pigment granules in wing scales and gives rise to white wing color. Lipid and transcriptome analyses reveal physiological differences associated with the ALHS. Together, these findings characterize a mechanism for a female-limited ALHS.

Tradeoffs are central to life history theory and evolutionary biology, yet almost nothing is known about their mechanistic basis. Here the authors characterize one such mechanism and find a transposable element insertion is associated with the switch between alternative life history strategies.

Disruptive selection of shell colour in land snails: a mark–recapture study of Euhadra peliomphala simodae

Many theoretical studies have suggested that disruptive selection plays an important role in phenotypic divergence, but few studies have determined the action of disruptive selection on phenotypic divergence via field studies. This study investigated the effect of disruptive selection on shell colour polymorphism in the Japanese land snail Euhadra peliomphala simodae to determine whether extreme phenotypes of snail shell colour are favoured over intermediate phenotypes. We conducted field surveys on an oceanic island with black, yellow and intermediate-coloured E. p. simodae snails. We captured and marked ~1800 individual snails and monitored their survival over 18 months. We quantified shell colours against images and examined the frequency distribution of shell colour variation. The variation exhibited a bimodal distribution with a far lower frequency of intermediate-coloured snails than of black or yellow snails. The population sizes of the three snail groups fluctuated synchronously with the changing seasons. Bayesian estimates showed lower survival rates for juvenile intermediate-coloured snails than for juvenile black and yellow snails, implying there was disruptive selection associated with shell colour. We suggest this disruptive selection may have resulted in the evolutionary divergence of the snail’s shell colour within the lineage having high shell colour variation.

The role of structural genomic variants in population differentiation and ecotype formation in Timema cristinae walking sticks

Theory predicts that structural genomic variants such as inversions can promote adaptive diversification and speciation. Despite increasing empirical evidence that adaptive divergence can be triggered by one or a few large inversions, the degree to which widespread genomic regions under divergent selection are associated with structural variants remains unclear. Here we test for an association between structural variants and genomic regions that underlie parallel host-plant-associated ecotype formation in Timema cristinae stick insects. Using mate-pair resequencing of 20 new whole genomes we find that moderately sized structural variants such as inversions, deletions and duplications are widespread across the genome, being retained as standing variation within and among populations. Using 160 previously published, standard-orientation whole genome sequences we find little to no evidence that the DNA sequences within inversions exhibit accentuated differentiation between ecotypes. In contrast, a formerly described large region of reduced recombination that harbours genes controlling colour-pattern exhibits evidence for accentuated differentiation between ecotypes, which is consistent with differences in the frequency of colour-pattern morphs between host-associated ecotypes. Our results suggest that some types of structural variants (e.g., large inversions) are more likely to underlie adaptive divergence than others, and that structural variants are not required for subtle yet genome-wide genetic differentiation with gene flow.

An Ultra High-Density Arabidopsis thaliana Crossover Map That Refines the Influences of Structural Variation and Epigenetic Features

Many environmental, genetic, and epigenetic factors are known to affect the frequency and positioning of meiotic crossovers (COs). Suppression of COs by large, cytologically visible inversions and translocations has long been recognized, but relatively little is known about how smaller structural variants (SVs) affect COs. To examine fine-scale determinants of the CO landscape, including SVs, we used a rapid, cost-effective method for high-throughput sequencing to generate a precise map of >17,000 COs between the Col-0 and Ler-0 accessions of Arabidopsis thaliana COs were generally suppressed in regions with SVs, but this effect did not depend on the size of the variant region, and was only marginally affected by the variant type. CO suppression did not extend far beyond the SV borders and CO rates were slightly elevated in the flanking regions. Disease resistance gene clusters, which often exist as SVs, exhibited high CO rates at some loci, but there was a tendency toward depressed CO rates at loci where large structural differences exist between the two parents. Our high-density map also revealed in fine detail how CO positioning relates to genetic (DNA motifs) and epigenetic (chromatin structure) features of the genome. We conclude that suppression of COs occurs over a narrow region spanning large- and small-scale SVs, representing an influence on the CO landscape in addition to sequence and epigenetic variation along chromosomes.

A Developmental Switch Generating Phenotypic Plasticity Is Part of a Conserved Multi-gene Locus

Switching between alternative complex phenotypes is often regulated by "supergenes," polymorphic clusters of linked genes such as in butterfly mimicry. In contrast, phenotypic plasticity results in alternative complex phenotypes controlled by environmental influences rather than polymorphisms. Here, we show that the developmental switch gene regulating predatory versus non-predatory mouth-form plasticity in the nematode Pristionchus pacificus is part of a multi-gene locus containing two sulfatases and two α-N-acetylglucosaminidases (nag). We provide functional characterization of all four genes, using CRISPR-Cas9-based reverse genetics, and show that nag genes and the previously identified eud-1/sulfatase have opposing influences. Members of the multi-gene locus show non-overlapping neuronal expression and epistatic relationships. The locus architecture is conserved in the entire genus Pristionchus. Interestingly, divergence between paralogs is counteracted by gene conversion, as inferred from phylogenies and genotypes of CRISPR-Cas9-induced mutants. Thus, we found that physical linkage accompanies regulatory linkage between switch genes controlling plasticity in P. pacificus.

A chromosome-level assembly of the Atlantic herring genome-detection of a supergene and other signals of selection

The Atlantic herring is a model species for exploring the genetic basis for ecological adaptation, due to its huge population size and extremely low genetic differentiation at selectively neutral loci. However, such studies have so far been hampered because of a highly fragmented genome assembly. Here, we deliver a chromosome-level genome assembly based on a hybrid approach combining a de novo Pacific Biosciences (PacBio) assembly with Hi-C-supported scaffolding. The assembly comprises 26 autosomes with sizes ranging from 12.4 to 33.1 Mb and a total size, in chromosomes, of 726 Mb, which has been corroborated by a high-resolution linkage map. A comparison between the herring genome assembly with other high-quality assemblies from bony fishes revealed few inter-chromosomal but frequent intra-chromosomal rearrangements. The improved assembly facilitates analysis of previously intractable large-scale structural variation, allowing, for example, the detection of a 7.8-Mb inversion on Chromosome 12 underlying ecological adaptation. This supergene shows strong genetic differentiation between populations. The chromosome-based assembly also markedly improves the interpretation of previously detected signals of selection, allowing us to reveal hundreds of independent loci associated with ecological adaptation.

Genetic coupling of signal and preference facilitates sexual isolation during rapid speciation

The divergence of sexual signals is ultimately a coevolutionary process: while signals and preferences diverge between lineages, they must remain coordinated within lineages for matings to occur. Divergence in sexual signals makes a major contribution to evolving species barriers. Therefore, the genetic architecture underlying signal-preference coevolution is essential to understanding speciation but remains largely unknown. In Laupala crickets where male song pulse rate and female pulse rate preferences have coevolved repeatedly and rapidly, we tested two contrasting hypotheses for the genetic architecture underlying signal-preference coevolution: linkage disequilibrium between unlinked loci and genetic coupling (linkage disequilibrium resulting from pleiotropy of a shared locus or tight physical linkage). Through selective introgression and quantitative trait locus (QTL) fine mapping, we estimated the location of QTL underlying interspecific variation in both female preference and male pulse rate from the same mapping populations. Remarkably, map estimates of the pulse rate and preference loci are as close as 0.06 cM apart, the strongest evidence to date for genetic coupling between signal and preference loci. As the second pair of colocalizing signal and preference loci in the Laupala genome, our finding supports an intriguing pattern, pointing to a major role for genetic coupling in the quantitative evolution of a reproductive barrier and rapid speciation in Laupala. Owing to its effect on suppressing recombination, a coupled, quantitative genetic architecture offers a powerful and parsimonious genetic mechanism for signal-preference coevolution and the establishment of positive genetic covariance on which the Fisherian runaway process of sexual selection relies.

Parent and offspring genotypes influence gene expression in early life

Parents can have profound effects on offspring fitness. Little, however, is known about the mechanisms through which parental genetic variation influences offspring physiology in natural systems. White-throated sparrows (Zonotrichia albicollis, WTSP) exist in two genetic morphs, tan and white, controlled by a large polymorphic supergene. Morphs mate disassortatively, resulting in two pair types: tan male × white female (T × W) pairs, which provide biparental care and white male × tan female (W × T) pairs, which provide female-biased care. To investigate how parental composition impacts offspring, we performed RNA-seq on whole blood of WTSP nestlings sampled from nests of both pair types. Parental pair type had a large effect on nestling gene expression, with 881 genes differentially expressed (DE) and seven correlated gene coexpression modules. The DE genes and modules expressed at higher levels in W × T nests with female-biased parental care function in metabolism and stress-related pathways resulting from the overrepresentation of proteolysis and stress-response genes (e.g., SOD2, NR3C1). These results show that parental genotypes and/or associated behaviours influence nestling physiology, and highlight avenues of further research investigating the ultimate implications for the maintenance of this polymorphism. Nestlings also exhibited morph-specific gene expression, with 92 differentially expressed genes, comprising immunity genes and genes encompassed by the supergene. Remarkably, we identified the same regulatory hub genes in these blood-derived expression networks as were previously identified in adult WTSP brains (EPM2A, BPNT1, TAF5L). These hub genes were located within the supergene, highlighting the importance of this gene complex in structuring regulatory networks across diverse tissues.

Unravelling the genes forming the wing pattern supergene in the polymorphic butterfly Heliconius numata

Background

Unravelling the genetic basis of polymorphic characters is central to our understanding of the origins and diversification of living organisms. Recently, supergenes have been implicated in a wide range of complex polymorphisms, from adaptive colouration in butterflies and fish to reproductive strategies in birds and plants. The concept of a supergene is now a hot topic in biology, and identification of its functional elements is needed to shed light on the evolution of highly divergent adaptive traits. Here, we apply different gene expression analyses to study the supergene P that controls polymorphism of mimetic wing colour patterns in the neotropical butterfly Heliconius numata.

Results

We performed de novo transcriptome assembly and differential expression analyses using high-throughput Illumina RNA sequencing on developing wing discs of different H. numata morphs. Within the P interval, 30 and 17 of the 191 transcripts were expressed differentially in prepupae and day-1 pupae, respectively. Among these is the gene cortex, known to play a role in wing pattern formation in Heliconius and other Lepidoptera. Our in situ hybridization experiments confirmed the relationship between cortex expression and adult wing patterns.

Conclusions

This study found the majority of genes in the P interval to be expressed in the developing wing discs during the critical stages of colour pattern formation, and detect drastic changes in expression patterns in multiple genes associated with structural variants. The patterns of expression of cortex only partially recapitulate the variation in adult phenotype, suggesting that the remaining phenotypic variation could be controlled by other genes within the P interval. Although functional studies on cortex are now needed to determine its exact developmental role, our results are in accordance with the classical supergene hypothesis, whereby several genes inherited together due to tight linkage control a major developmental switch.

Electronic supplementary material

The online version of this article (10.1186/s13227-019-0129-2) contains supplementary material, which is available to authorized users.

Ecological factors and morphological traits are associated with repeated genomic differentiation between lake and stream stickleback

The repeated evolution of similar phenotypes in independent populations (i.e. parallel or convergent evolution) provides an opportunity to identify genetic and ecological factors that influence the process of adaptation. Threespine stickleback fish ( Gasterosteus aculeatus) are an excellent model for such studies, as they have repeatedly adapted to divergent habitats across the Northern hemisphere. Here, we use genomic, ecological and morphological data from 16 independent pairs of stickleback populations adapted to divergent lake and stream habitats. We combine a population genomic approach to identify regions of the genome that are likely under selection in these divergent habitats with an association mapping approach to identify regions of the genome that underlie variation in ecological factors and morphological traits. Over 37% of genomic windows are repeatedly differentiated across lake-stream pairs. Similarly, many genomic windows are associated with variation in abiotic factors, diet items and morphological phenotypes. Both the highly differentiated windows and candidate trait windows are non-randomly distributed across the genome and show some overlap. However, the overlap is not significant on a genome-wide scale. Together, our data suggest that adaptation to divergent food resources and predation regimes are drivers of differentiation in lake-stream stickleback, but that additional ecological factors are also important. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Molecular Characterisation of a Supergene Conditioning Super-High Vitamin C in Kiwifruit Hybrids

During analysis of kiwifruit derived from hybrids between the high vitamin C (ascorbic acid; AsA) species Actinidia eriantha and A. chinensis, we observed bimodal segregation of fruit AsA concentration suggesting major gene segregation. To test this hypothesis, we performed whole-genome sequencing on pools of hybrid genotypes with either high or low AsA fruit. Pool-GWAS (genome-wide association study) revealed a single Quantitative Trait Locus (QTL) spanning more than 5 Mbp on chromosome 26, which we denote as qAsA26.1. A co-dominant PCR marker was used to validate this association in four diploid (A. chinensis × A. eriantha) × A. chinensis backcross families, showing that the A. eriantha allele at this locus increases fruit AsA levels by 250 mg/100 g fresh weight. Inspection of genome composition and recombination in other A. chinensis genetic maps confirmed that the qAsA26.1 region bears hallmarks of suppressed recombination. The molecular fingerprint of this locus was examined in leaves of backcross validation families by RNA sequencing (RNASEQ). This confirmed strong allelic expression bias across this region as well as differential expression of transcripts on other chromosomes. This evidence suggests that the region harbouring qAsA26.1 constitutes a supergene, which may condition multiple pleiotropic effects on metabolism.

The mimetic wing pattern of Papilio polytes butterflies is regulated by a doublesex-orchestrated gene network

The swallowtail butterfly Papilio polytes is sexually dimorphic and exhibits female-limited Batesian mimicry. This species also has two female forms, a non-mimetic form with male-like wing patterns, and a mimetic form resembling an unpalatable model, Pachliopta aristolochiae. The mimicry locus H constitutes a dimorphic Mendelian 'supergene', including a transcription factor gene doublesex (dsx). However, how the mimetic-type dsx (dsx-H) orchestrates the downstream gene network and causes the mimetic traits remains unclear. Here we performed RNA-seq-based gene screening and found that Wnt1 and Wnt6 are up-regulated by dsx-H during the early pupal stage and are involved in the red/white pigmentation and patterning of mimetic female wings. In contrast, a homeobox gene abdominal-A is repressed by dsx-H and involved in the non-mimetic colouration pattern. These findings suggest that dual regulation by dsx-H, induction of mimetic gene networks and repression of non-mimetic gene networks, is essential for the switch from non-mimetic to mimetic pattern in mimetic female wings.

Supergenes

Daniella Black and David M. Shuker introduce supergenes - clusters of genes that control complex traits and are inherited together.

Degenerative Expansion of a Young Supergene

Long-term suppression of recombination ultimately leads to gene loss, as demonstrated by the depauperate Y and W chromosomes of long-established pairs of XY and ZW chromosomes. The young social supergene of the Solenopsis invicta red fire ant provides a powerful system to examine the effects of suppressed recombination over a shorter timescale. The two variants of this supergene are carried by a pair of heteromorphic chromosomes, referred to as the social B and social b (SB and Sb) chromosomes. The Sb variant of this supergene changes colony social organization and has an inheritance pattern similar to a Y or W chromosome because it is unable to recombine. We used high-resolution optical mapping, k-mer distribution analysis, and quantification of repetitive elements on haploid ants carrying alternate variants of this young supergene region. We find that instead of shrinking, the Sb variant of the supergene has increased in length by more than 30%. Surprisingly, only a portion of this length increase is due to consistent increases in the frequency of particular classes of repetitive elements. Instead, haplotypes of this supergene variant differ dramatically in the amounts of other repetitive elements, indicating that the accumulation of repetitive elements is a heterogeneous and dynamic process. This is the first comprehensive demonstration of degenerative expansion in an animal and shows that it occurs through nonlinear processes during the early evolution of a region of suppressed recombination.

The Origin of a New Sex Chromosome by Introgression between Two Stickleback Fishes

Introgression is increasingly recognized as a source of genetic diversity that fuels adaptation. Its role in the evolution of sex chromosomes, however, is not well known. Here, we confirm the hypothesis that the Y chromosome in the ninespine stickleback, Pungitius pungitius, was established by introgression from the Amur stickleback, P. sinensis. Using whole genome resequencing, we identified a large region of Chr 12 in P. pungitius that is diverged between males and females. Within but not outside of this region, several lines of evidence show that the Y chromosome of P. pungitius shares a most recent common ancestor not with the X chromosome, but with the homologous chromosome in P. sinensis. Accumulation of repetitive elements and gene expression changes on the new Y are consistent with a young sex chromosome in early stages of degeneration, but other hallmarks of Y chromosomes have not yet appeared. Our findings indicate that porous species boundaries can trigger rapid sex chromosome evolution.

Consequences of Asexuality in Natural Populations: Insights from Stick Insects

Recombination is a fundamental process with significant impacts on genome evolution. Predicted consequences of the loss of recombination include a reduced effectiveness of selection, changes in the amount of neutral polymorphisms segregating in populations, and an arrest of GC-biased gene conversion. Although these consequences are empirically well documented for nonrecombining genome portions, it remains largely unknown if they extend to the whole genome scale in asexual organisms. We identify the consequences of asexuality using de novo transcriptomes of five independently derived, obligately asexual lineages of stick insects, and their sexual sister-species. We find strong evidence for higher rates of deleterious mutation accumulation, lower levels of segregating polymorphisms and arrested GC-biased gene conversion in asexuals as compared with sexuals. Taken together, our study conclusively shows that predicted consequences of genome evolution under asexuality can indeed be found in natural populations.

The Effect of Neutral Recombination Variation on Genome Scans for Selection

Recently, there has been an increasing interest in identifying the role that regions of low recombination or inversion play in adaptation of species to local environments. Many examples of groups of adapted genes located within inversions are arising in the literature, in part inspired by theory that predicts the evolution of these so-called “supergenes.” We still, however, have a poor understanding of how genomic heterogeneity, such as varying rates of recombination, may confound signals of selection. Here, I evaluate the effect of neutral inversions and recombination variation on genome scans for selection, including tests for selective sweeps, differentiation outlier tests, and association tests. There is considerable variation among methods in their performance, with some methods being unaffected and some showing elevated false positive signals within a neutral inversion or region of low recombination. In some cases the false positive signal can be dampened or removed, if it is possible to use a quasi-independent set of SNPs to parameterize the model before performing the test. These results will be helpful to those seeking to understand the importance of regions of low recombination in adaptation.

Coevolution of Genome Architecture and Social Behavior

Although social behavior can have a strong genetic component, it can also result in selection on genome structure and function, thereby influencing the evolution of the genome itself. Here we explore the bidirectional links between social behavior and genome architecture by considering variation in social and/or mating behavior among populations (social polymorphisms) and across closely related species. We propose that social behavior can influence genome architecture via associated demographic changes due to social living. We establish guidelines to exploit emerging whole-genome sequences using analytical approaches that examine genome structure and function at different levels (regulatory vs structural variation) from the perspective of both molecular biology and population genetics in an ecological context.

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.