The role of microbial symbionts in governing social traits of their hosts is an exciting and developing research area. Just like symbionts influence host reproductive behaviour and can cause mating incompatibilities to promote symbiont transmission through host populations (Engelstadter and Hurst 2009; Correa and Ballard 2016; Johnson and Foster 2018) (see also discussion on conflict resolution in Johnsen and Foster 2018), microbial symbionts could enhance transmission by promoting the social behaviour of their hosts (Ezenwa et al. 2012; Lewin-Epstein et al. 2017; Gurevich et al. 2020). Here I apply the term ‘symbiosis’ in the broad sense, following De Bary 1879 as “the living together of two differently named organisms“ independent of effects on the organisms involved (De Bary 1879), i.e. the biological interaction between the host and its symbionts may include mutualism, parasitism and commensalism.
So far, we have relative few studies that explore the role of symbionts in promoting social behaviours such as parental care. Clearly, disentangling cause and effect when assessing the functional significance of symbiotic relationships in general is extremely challenging, and perhaps even more caution is needed when assessing the role of symbionts in the evolution of parental care, given the high fitness benefits to the offspring of receiving care. An interesting study on the symbiotic relationship between termites and their eukaryotic gut symbionts proposes a role of gut flagellates in the origin of subsocial behaviour (extended offspring care) in the termites through proctodeal trophallaxis (i.e. anus-to-mouth feeding), driven by mutualistic beneficial interactions (Nalepa 2020). Van Meyel et al. (2021) hypothesized a role of gut symbionts in promoting maternal care behaviour in the European earwig, and set out to test this idea in a carefully executed experimental study. They used a broad-spectrum antibiotic treatment to alter gut microbiota in mothers and examined its effect on maternal care provisioning. While the antibiotic treatment altered the gut microbiome, no effect on pre- or post-hatching maternal care was detected. The authors also investigated a broad range of physiological and reproductive traits measured over a major part of a female’s lifetime, and detected no effect of microbiome alteration on these traits. The study therefore does not provide evidence for a direct role of the gut microbiome in shaping offspring care in this population of European earwigs.
Within populations, earwigs show inter-individual variation in the expression of maternal care (Meunier et al. 2012; Ratz et al. 2016), and there is evidence that genetic and environmental factors contribute to this this variation (Meunier and Kolliker 2012; Kramer et al. 2017). The study by Van Meyel et al. (2021) is the first to analyse microbiome composition of the European earwig, and they study host-symbiont associations in a single population. A next step could be to explore among population variation in the gut microbiome, to achieve a better understanding on host-microbiome variation and dynamics in wild populations. Depending on the nature of host-symbiont associations across populations, new perspectives on their functional significance may arise (Hird 2017; Johnson and Foster 2018). It is therefore too early to conclusively confirm or reject the role of microbial symbionts in the expression of parental care in this system.

References

Correa, C. C., and Ballard, J. W. O. (2016). Wolbachia associations with insects: winning or losing against a master manipulator. Frontiers in Ecology and Evolution, 3, 153. doi: https://doi.org/10.3389/fevo.2015.00153

De Bary, A. (1879). Die Erscheinung der Symbiose. Verlag von Karl J. Trubner, Strassburg.

Engelstädter, J., and Hurst, G. D. (2009). The ecology and evolution of microbes that manipulate host reproduction. Annual Review of Ecology, Evolution, and Systematics, 40, 127-149. doi: https://doi.org/10.1146/annurev.ecolsys.110308.120206

Ezenwa, V. O., Gerardo, N. M., Inouye, D. W., Medina, M., and Xavier, J. B. (2012). Animal behavior and the microbiome. Science, 338(6104), 198-199. doi: https://doi.org/10.1126/science.1227412

Gurevich, Y., Lewin-Epstein, O., and Hadany, L. (2020). The evolution of paternal care: a role for microbes?. Philosophical Transactions of the Royal Society B, 375(1808), 20190599. doi: https://doi.org/10.1098/rstb.2019.0599

Hird, S. M. (2017). Evolutionary biology needs wild microbiomes. Frontiers in microbiology, 8, 725. doi: https://doi.org/10.3389/fmicb.2017.00725

Johnson, K. V. A., and Foster, K. R. (2018). Why does the microbiome affect behaviour?. Nature reviews microbiology, 16(10), 647-655. doi: https://doi.org/10.1038/s41579-018-0014-3

Kramer et al. (2017). When earwig mothers do not care to share: parent–offspring competition and the evolution of family life. Functional Ecology, 31(11), 2098-2107. doi: https://doi.org/10.1111/1365-2435.12915

Lewin-Epstein, O., Aharonov, R., and Hadany, L. (2017). Microbes can help explain the evolution of host altruism. Nature communications, 8(1), 1-7. doi: https://doi.org/10.1038/ncomms14040

Meunier, J., and Kölliker, M. (2012). Parental antagonism and parent–offspring co-adaptation interact to shape family life. Proceedings of the Royal Society B: Biological Sciences, 279(1744), 3981-3988. doi: https://doi.org/10.1098/rspb.2012.1416

Meunier, J., Wong, J. W., Gómez, Y., Kuttler, S., Röllin, L., Stucki, D., and Kölliker, M. (2012). One clutch or two clutches? Fitness correlates of coexisting alternative female life-histories in the European earwig. Evolutionary Ecology, 26(3), 669-682. doi: https://doi.org/10.1007/s10682-011-9510-x

Nalepa, C. A. (2020). Origin of mutualism between termites and flagellated gut protists: transition from horizontal to vertical transmission. Frontiers in Ecology and Evolution, 8, 14. doi: https://doi.org/10.3389/fevo.2020.00014

Ratz, T., Kramer, J., Veuille, M., and Meunier, J. (2016). The population determines whether and how life-history traits vary between reproductive events in an insect with maternal care. Oecologia, 182(2), 443-452. doi: https://doi.org/10.1007/s00442-016-3685-3

Van Meyel, S., Devers, S., Dupont, S., Dedeine, F. and Meunier, J. (2021) Alteration of gut microbiota with a broad-spectrum antibiotic does not impair maternal care in the European earwig. bioRxiv, 2020.10.08.331363. ver. 5 peer-reviewed and recommended by PCI Evol Biol. https://doi.org/10.1101/2020.10.08.331363

Despite a vast array of ways in which organisms can reproduce (Bell, 1982), most animals engage in sexual reproduction (Otto & Lenormand, 2002). A fascinating alternative to sex is parthenogenesis, where offspring are produced asexually from a gamete, typically the egg, without receiving genetic material from another gamete (Simon, Delmotte, Rispe, & Crease, 2003). One of the long-standing questions in the field is why parthenogenesis is not more widespread, given the costs associated with sex (Otto & Lenormand, 2002).  Natural populations of most species appear to be reproducing either sexually or parthenogenetically, even if a species can employ both reproductive modes (Larose et al 2023). Larose et al (2023) highlight the conundrum in this pattern, as organisms that are capable of employing parthenogenesis facultatively would be able to gain the benefits of both modes of reproduction. Why then, is facultative parthenogenesis not more common?

Larose et al (2023) propose that constraints on being efficient in both sexual and asexual reproduction could cause a trade-off between reproductive modes that favours an obligate strategy of either sex or no sex. This would provide an explanation for why facultative parthenogenesis is rare. Timema stick insects provide an excellent system to investigate reproductive strategies, as some species have parthenogenetic females, while other species are sexual, and they show repeated transitions from sexual reproduction to obligate parthenogenesis (Schwander & Crespi, 2009). The authors performed comprehensive and complementary studies in a recently discovered species T. douglasi, in which populations show both modes of reproduction, with some populations consisting only of females and others showing equal proportions of males and females. The sex ratio varied significantly, with the proportion of females ranging between 43-100% across 29 populations. These populations form a monophyletic clade with clustering into three genetic lineages and only a few cases of admixture. Females from all populations were capable of producing unfertilized eggs, but the hatching success varied hugely among populations and lineages (3-100%). Parthenogenetically produced offspring were homozygous, showing that parthenogenesis causes a complete loss of heterozygosity in a single generation. After producing eggs as virgins, females were mated to assess the capacity to also reproduce sexually, and fertilization increased the hatching success of eggs in two lineages. In one lineage, in which the hatching success of unfertilized eggs is similar to that of other sexually reproducing Timema species, fertilization reduced egg-hatching success, indicating a trade-off between reproductive modes with parthenogenetic reproduction performing best. Approximately 58% of the offspring produced after mating were fertilized, demonstrating the capacity of females to reproduce parthenogenetically also after mating has occurred, however with huge variation among individuals.

This wonderful and meticulously performed study produces strong and complementary evidence for facultative parthenogenesis in T. douglasi populations. The study shows large variation in how reproductive mode is employed, supporting the existence of a trade-off between sexual and parthenogenetic reproduction. This might be an example of an ongoing transition from sexual to asexual reproduction, which indicates that obligate parthenogenesis may derive via transient facultative parthenogenesis. 

REFERENCES

Bell, G. (1982) The Masterpiece of Nature: The Evolution and Genetics of Sexuality. University of California Press. 635 p.

Otto, S. P., & Lenormand, T. (2002). Resolving the paradox of sex and recombination. Nature Reviews Genetics, 3(4), 252-261. https://doi.org/10.1038/nrg761

Schwander, T., & Crespi, B. J. (2009). Multiple direct transitions from sexual reproduction to apomictic parthenogenesis in Timema stick insects. Evolution, 63(1), 84-103. 
https://doi.org/10.1111/j.1558-5646.2008.00524.x

Simon, J.-C., Delmotte, F., Rispe, C., & Crease, T. (2003). Phylogenetic relationships between parthenogens and their sexual relatives: the possible routes to parthenogenesis in animals. Biological Journal of the Linnean Society, 79(1), 151-163. https://doi.org/10.1046/j.1095-8312.2003.00175.x

Larose, C., Lavanchy,  G., Freitas, S., Parker, D.J., Schwander, T. (2023) Facultative parthenogenesis: a transient state in transitions between sex and obligate asexuality in stick insects? bioRxiv, 2022.03.25.485836, ver. 4 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.03.25.485836

The advent of genome-wide association studies (GWAS) has been a great promise for our understanding of the connection between genotype and phenotype. Today, the NHGRI-EBI GWAS catalog contains 251,401 associations from 4,961 studies (1). This wealth of studies has also generated interest to use the summary statistics beyond the few top hits in order to make predictions for individuals without known phenotype, e.g. to predict polygenic risk scores or to study polygenic selection by comparing different groups. For instance, polygenic selection acting on the most studied polygenic trait, height, has been subject to multiple studies during the past decade (e.g. 2–6). They detected north-south gradients in Europe which were consistent with expectations. However, their GWAS summary statistics were based on the GIANT consortium data set, a meta-analysis of GWAS conducted in different European cohorts (7,8). The availability of large data sets with less stratification such as the UK Biobank (9) has led to a re-evaluation of those results. The nature of the GIANT consortium data set was realized to represent a potential problem for studies of polygenic adaptation which led several of the authors of the original articles to caution against the interpretations of polygenic selection on height (10,11). This was a great example on how the scientific community assessed their own earlier results in a critical way as more data became available. At the same time it left the question whether there is detectable polygenic selection separating populations more open than ever.

Generally, recent years have seen several articles critically assessing the portability of GWAS results and risk score predictions to other populations (12–14). Refoyo-Martínez et al. (15) are now presenting a systematic assessment on the robustness of cross-population signatures of polygenic adaptation in humans. They compiled GWAS results for complex traits which have been studied in more than one cohort and then use allele frequencies from the 1000 Genomes Project data (16) set to detect signals of polygenic score overdispersion. As the source for the allele frequencies is kept the same across all tests, differences between the signals must be caused by the underlying GWAS. The results are concerning as the level of overdispersion largely depends on the choice of GWAS cohort. Cohorts with homogenous ancestries show little to no overdispersion compared to cohorts of mixed ancestries such as meta-analyses. It appears that the meta-analyses fail to fully account for stratification in their data sets.

The authors based most of their analyses on the heavily studied trait height. Additionally, they use educational attainment (measured as the number of school years of an individual) as an example. This choice was due to the potential over- or misinterpretation of results by the media, the general public and by far right hate groups. Such traits are potentially confounded by unaccounted cultural and socio-economic factors. Showing that previous results about polygenic selection on educational attainment are not robust is an important result that needs to be communicated well. This forms a great example for everyone working in human genomics. We need to be aware that our results can sometimes be misinterpreted. And we need to make an effort to write our papers and communicate our results in a way that is honest about the limitations of our research and that prevents the misuse of our results by hate groups.

This article represents an important contribution to the field. It is cruicial to be aware of potential methodological biases and technical artifacts. Future studies of polygenic adaptation need to be cautious with their interpretations of polygenic score overdispersion. A recommendation would be to use GWAS results obtained in homogenous cohorts. But even if different biobank-scale cohorts of homogeneous ancestry are employed, there will always be some remaining risk of unaccounted stratification. These conclusions may seem sobering but they are part of the scientific process. We need additional controls and new, different methods than polygenic score overdispersion for assessing polygenic selection. Last year also saw the presentation of a novel approach using sequence data and GWAS summary statistics to detect directional selection on a polygenic trait (17). This new method appears to be robust to bias stemming from stratification in the GWAS cohort as well as other confounding factors. Such new developments show light at the end of the tunnel for the use of GWAS summary statistics in the study of polygenic adaptation.

References

1. Buniello A, MacArthur JAL, Cerezo M, Harris LW, Hayhurst J, Malangone C, et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Research. 2019 Jan 8;47(D1):D1005–12. doi: https://doi.org/10.1093/nar/gky1120

2. Turchin MC, Chiang CW, Palmer CD, Sankararaman S, Reich D, Hirschhorn JN. Evidence of widespread selection on standing variation in Europe at height-associated SNPs. Nature Genetics. 2012 Sep;44(9):1015–9. doi: https://doi.org/10.1038/ng.2368

3. Berg JJ, Coop G. A Population Genetic Signal of Polygenic Adaptation. PLOS Genetics. 2014 Aug 7;10(8):e1004412. doi: https://doi.org/10.1371/journal.pgen.1004412

4. Robinson MR, Hemani G, Medina-Gomez C, Mezzavilla M, Esko T, Shakhbazov K, et al. Population genetic differentiation of height and body mass index across Europe. Nature Genetics. 2015 Nov;47(11):1357–62. doi: https://doi.org/10.1038/ng.3401

5. Mathieson I, Lazaridis I, Rohland N, Mallick S, Patterson N, Roodenberg SA, et al. Genome-wide patterns of selection in 230 ancient Eurasians. Nature. 2015 Dec;528(7583):499–503. doi: https://doi.org/10.1038/nature16152

6. Racimo F, Berg JJ, Pickrell JK. Detecting polygenic adaptation in admixture graphs. Genetics. 2018. Arp;208(4):1565–1584. doi: https://doi.org/10.1534/genetics.117.300489

7. Lango Allen H, Estrada K, Lettre G, Berndt SI, Weedon MN, Rivadeneira F, et al. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature. 2010 Oct;467(7317):832–8. doi: https://doi.org/10.1038/nature09410

8. Wood AR, Esko T, Yang J, Vedantam S, Pers TH, Gustafsson S, et al. Defining the role of common variation in the genomic and biological architecture of adult human height. Nat Genet. 2014 Nov;46(11):1173–86. doi: https://doi.org/10.1038/ng.3097

9. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018 Oct;562(7726):203–9. doi: https://doi.org/10.1038/s41586-018-0579-z

10. Berg JJ, Harpak A, Sinnott-Armstrong N, Joergensen AM, Mostafavi H, Field Y, et al. Reduced signal for polygenic adaptation of height in UK Biobank. eLife. 2019 Mar 21;8:e39725. doi: https://doi.org/10.7554/eLife.39725

11. Sohail M, Maier RM, Ganna A, Bloemendal A, Martin AR, Turchin MC, et al. Polygenic adaptation on height is overestimated due to uncorrected stratification in genome-wide association studies. eLife. 2019 Mar 21;8:e39702. doi: https://doi.org/10.7554/eLife.39702

12. Martin AR, Kanai M, Kamatani Y, Okada Y, Neale BM, Daly MJ. Clinical use of current polygenic risk scores may exacerbate health disparities. Nature Genetics. 2019 Apr;51(4):584–91. doi: https://doi.org/10.1038/s41588-019-0379-x

13. Bitarello BD, Mathieson I. Polygenic Scores for Height in Admixed Populations. G3: Genes, Genomes, Genetics. 2020 Nov 1;10(11):4027–36. doi: https://doi.org/10.1534/g3.120.401658

14. Uricchio LH, Kitano HC, Gusev A, Zaitlen NA. An evolutionary compass for detecting signals of polygenic selection and mutational bias. Evolution Letters. 2019;3(1):69–79. doi: https://doi.org/10.1002/evl3.97

15. Refoyo-Martínez A, Liu S, Jørgensen AM, Jin X, Albrechtsen A, Martin AR, Racimo F. How robust are cross-population signatures of polygenic adaptation in humans? bioRxiv, 2021, 2020.07.13.200030, version 5 peer-reviewed and recommended by Peer community in Evolutionary Biology. doi: https://doi.org/10.1101/2020.07.13.200030

16. Auton A, Abecasis GR, Altshuler DM, Durbin RM, Abecasis GR, Bentley DR, et al. A global reference for human genetic variation. Nature. 2015 Sep 30;526(7571):68–74. doi: https://doi.org/10.1038/nature15393

17. Stern AJ, Speidel L, Zaitlen NA, Nielsen R. Disentangling selection on genetically correlated polygenic traits using whole-genome genealogies. bioRxiv. 2020 May 8;2020.05.07.083402. doi: https://doi.org/10.1101/2020.05.07.083402

Much research effort has been extended toward developing systems for managing soil inhabiting insect pests of crops with entomopathogenic nematodes as biocontrol agents. Although small plot or laboratory experiments may suggest a particular insect pest is vulnerable to management in this way, it is often difficult to scale-up nematode production for application at the field- and farm scale to make such a tactic viable. Part of the problem is that entomopathogenic nematode strains must be propagated by serial passage in vivo, because storage by freezing decreases fitness. At the same time, serial propagation results in loss of virulence (ability to infect) over generations in the laboratory, a phenomenon called attenuation.

To probe the underlying reasons for development of attenuation, as a prerequisite to designing strategies to mitigate it, Shapiro-Ilan and Raymond [1] turned to evolutionary theory of social conflict as a possible explanatory framework. Virulence of entomopathogenic nematodes depends on a combination of virulence factors, like various proteases, secreted by both the nematode and symbiotic bacteria to overcome host defenses. Attenuation is characterized in part by a reduced production of these factors. Invasion of a host involves simultaneous attack by a group of nematodes ("cooperators"), which together neutralize host defenses enough to allow individuals to successfully invade. "Cheaters" in the invading population can avoid the metabolic costs of producing virulence factors while reaping the benefits of infecting the host made vulnerable by the cooperators in the population. The authors hypothesize that an increase in frequency of cheaters may contribute to attenuation of virulence during serial propagation in the laboratory. The evolutionary dynamics of cheater frequency in a population have been explored in many contexts as part of kin selection theory. Cheaters can increase in a population by outcompeting cooperators in a host if overall relatedness within the invading population is low. Conversely, frequency of altruism, or costly cooperation, increases in a population if relatedness is high, which is enhanced by low effective dispersal. However, a population that is too isolated can suffer from inbreeding effects, and competition will occur mainly among relatives, which decreases the fitness benefits of altruism.

Shapiro-Ilan and Raymond [1] tested changes in virulence and reproductive output in a serially propagated entomopathogenic nematode, Heterorhabditis floridensis. They compared lines of high or low relatedness, manipulated via multiplicity of infection (MOI) rates (where a low dose of nematodes gives high relatedness and a high dose gives low relatedness); and under global or local competition, manipulated by pooling populations emerging from all or only two host cadavers per generation, respectively. As predicted, treatments of high relatedness (low MOI) and global competition had the greatest level of reproduction, while all lines of low relatedness (high MOI) evolved decreased reproduction and decreased virulence, which led to extinction. The key finding was that lines in the high relatedness (low MOI) and low (local) competition treatment exhibited the most stable virulence through the 12 generations tested. Thus, to minimize attenuation of virulence while maintaining fitness of recently isolated entomopathogenic nematodes, the authors recommend insect hosts be inoculated with low doses of nematodes from inocula pools from as few cadavers as possible.

The application of evolutionary theory, with a clever experimental design, to an important problem in pest management makes this paper particularly noteworthy.

Reference

[1] Shapiro-Ilan D, Raymond B. 2016. Limiting opportunities for cheating stabilizes virulence in insect parasitic nematodes. Evolutionary Applications 9:462-470. doi: 10.1111/eva.12348

Two back-to-back papers published earlier this year in Nature Genetics provide compelling evidence for the control of a male reproductive polymorphism in a wading bird by a "supergene", a cluster of tightly linked genes [1-2]. The bird in question, the ruff (Philomachus pugnax), has a rather unusual reproductive system that consists of three distinct types of males ("reproductive morphs"): aggressive "independents" who represent the majority of males; a smaller fraction of non-territorial "satellites" who are submissive towards "independents"; and "faeders" who mimic females and are rare. Previous work has shown that the male morphs differ in major aspects of mating and aggression behavior, plumage coloration and body size, and that – intriguingly – this complex multi-trait polymorphism is apparently controlled by a single autosomal Mendelian locus with three alleles [3]. To uncover the genetic control of this polymorphism two independent teams, led by Terry Burke [1] and Leif Andersson [2], have set out to analyze the genomes of male ruffs. Using a combination of genomics and genetics, both groups managed to pin down the supergene locus and map it to a non-recombining, 4.5 Mb large inversion which arose 3.8 million years ago. While "independents" are homozygous for the ancestral uninverted sequence, "satellites" and "faeders" carry evolutionarily divergent, dominant alternative haplotypes of the inversion. Thus, as in several other notable cases, for example the supergene control of disassortative mating, aggressiveness and plumage color in white-throated sparrows [4], of mimicry in Heliconius and Papilio butterflies [5-6], or of social structure in ants [7], an inversion – behaving as a single "locus" – underpins the mechanistic basis of the supergene. More generally, and beyond inversions, a growing number of studies now shows that selection can favor the evolution of suppressed recombination, thereby leading to the emergence of clusters of tightly linked loci which can then control – presumably due to polygenic gene action – a suite of complex phenotypes [8-10]. A largely unresolved question in this field concerns the identity of the causative alleles and loci within a given supergene. Recent progress on this question has been made for example in Papilio polytes butterflies where a mimicry supergene has been found to involve – surprisingly – only a single but large gene: multiple mimicry alleles in the doublesex gene are maintained in strong linkage disequilibrium via an inversion. It will clearly be of great interest to see future examples of such a fine-scale genetic dissection of supergenes. In conclusion, we were impressed by the data and analyses of Küpper et al. [1] and Lamichhaney et al. [2]: both papers beautifully illustrate how genomics and evolutionary ecology can be combined to make new, exciting discoveries. Both papers will appeal to readers with an interest in supergenes, inversions, the interplay of selection and recombination, or the genetic control of complex phenotypes.

References

[1] Küpper C, Stocks M, Risse JE, dos Remedios N, Farrell LL, McRae SB, Morgan TC, Karlionova N, Pinchuk P, Verkuil YI, et al. 2016. A supergene determines highly divergent male reproductive morphs in the ruff. Nature Genetics 48:79-83. doi: 10.1038/ng.3443

[2] Lamichhaney S, Fan G, Widemo F, Gunnarsson U, Thalmann DS, Hoeppner MP, Kerje S, Gustafson U, Shi C, Zhang H, et al. 2016. Structural genomic changes underlie alternative reproductive strategies in the ruff (Philomachus pugnax). Nature Genetics 48:84-88. doi: 10.1038/ng.3430

[3] Lank DB, Smith CM, Hanotte O, Burke T, Cooke F. 1995. Genetic polymorphism for alternative mating behaviour in lekking male ruff Philomachus pugnax. Nature 378:59-62. doi: 10.1038/378059a0

[4] Tuttle Elaina M, Bergland Alan O, Korody Marisa L, Brewer Michael S, Newhouse Daniel J, Minx P, Stager M, Betuel A, Cheviron Zachary A, Warren Wesley C, et al. 2016. Divergence and Functional Degradation of a Sex Chromosome-like Supergene. Current Biology 26:344-350. doi: 10.1016/j.cub.2015.11.069

[5] Joron M, Frezal L, Jones RT, Chamberlain NL, Lee SF, Haag CR, Whibley A, Becuwe M, Baxter SW, Ferguson L, et al. 2011. Chromosomal rearrangements maintain a polymorphic supergene controlling butterfly mimicry. Nature 477:203-206. doi: 10.1038/nature10341

[6] Kunte K, Zhang W, Tenger-Trolander A, Palmer DH, Martin A, Reed RD, Mullen SP, Kronforst MR. 2014. doublesex is a mimicry supergene. Nature 507:229-232. doi: 10.1038/nature13112

[7] Wang J, Wurm Y, Nipitwattanaphon M, Riba-Grognuz O, Huang Y-C, Shoemaker D, Keller L. 2013. A Y-like social chromosome causes alternative colony organization in fire ants. Nature 493:664-668. doi: 10.1038/nature11832

[8] Thompson MJ, Jiggins CD. 2014. Supergenes and their role in evolution. Heredity 113:1-8. doi: 10.1038/hdy.2014.20

[9] Schwander T, Libbrecht R, Keller L. 2014. Supergenes and Complex Phenotypes. Current Biology 24:R288-R294. doi: 10.1016/j.cub.2014.01.056

[10] Charlesworth D. 2015. The status of supergenes in the 21st century: recombination suppression in Batesian mimicry and sex chromosomes and other complex adaptations. Evolutionary Applications 9:74-90. doi: 10.1111/eva.12291