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

Comprehending the dynamism of B chromosomes in their journey towards becoming unselfish

Investigated for more than a century now, B chromosomes (Bs) research has come a long way from Bs being considered parasitic or neutral to becoming unselfish and bringing benefits to their hosts. B chromosomes exist as accessory chromosomes along with the standard A chromosomes (As) across eukaryotic taxa. Represented singly or in multiple copies, B chromosomes are largely heterochromatic but also contain euchromatic and organellar segments. Although B chromosomes are derived entities, they follow their species-specific evolutionary pattern. B chromosomes fail to pair with the standard chromosomes during meiosis and vary in their number, size, composition and structure across taxa and ensure their successful transmission through non-mendelian mechanisms like mitotic, pre-meiotic, meiotic or post-meiotic drives, unique non-disjunction, self-pairing or even imparting benefits to the host when they lack drive. B chromosomes have been associated with cellular processes like sex determination, pathogenicity, resistance to pathogens, phenotypic effects, and differential gene expression. With the advancements in B-omics research, novel insights have been gleaned on their functions, some of which have been associated with the regulation of gene expression of A chromosomes through increased expression of miRNAs or differential expression of transposable elements located on them. The next-generation sequencing and emerging technologies will further likely unravel the cellular, molecular and functional behaviour of these enigmatic entities. Amidst the extensive fluidity shown by B chromosomes in their structural and functional attributes, we perceive that the existence and survival of B chromosomes in the populations most likely seem to be a trade-off between the drive efficiency and adaptive significance versus their adverse effects on reproduction.

Widespread chromosomal rearrangements preceded genetic divergence in a monitor lizard, Varanus acanthurus (Varanidae)

Chromosomal rearrangements are often associated with local adaptation and speciation because they suppress recombination, and as a result, rearrangements have been implicated in disrupting gene flow. Although there is strong evidence to suggest that chromosome rearrangements are a factor in genetic isolation of divergent populations, the underlying mechanism remains elusive. Here, we applied an integrative cytogenetics and genomics approach testing whether chromosomal rearrangements are the initial process, or a consequence, of population divergence in the dwarf goanna, Varanus acanthurus. Specifically, we tested whether chromosome rearrangements are indicators of genetic barriers that can be used to identify divergent populations by looking at gene flow within and between populations with rearrangements. We found that gene flow was present between individuals with chromosome rearrangements within populations, but there was no gene flow between populations that had similar chromosome rearrangements. Moreover, we identified a correlation between reduced genetic variation in populations with a higher frequency of homozygous submetacentric individuals. These findings suggest that chromosomal rearrangements were widespread prior to divergence, and because we found populations with higher frequencies of submetacentric chromosomes were associated with lower genetic diversity, this could indicate that polymorphisms within populations are early indicators of genetic drift.

Allele-specific cis-regulatory methylation of the gene for vasoactive intestinal peptide in white-throated sparrows

White-throated sparrows (Zonotrichia albicollis) offer a unique opportunity to connect genotype with behavioral phenotype. In this species, a rearrangement of the second chromosome is linked with territorial aggression; birds with a copy of this "supergene" rearrangement are more aggressive than those without it. The supergene has captured the gene VIP, which encodes vasoactive intestinal peptide, a neuromodulator that drives aggression in other songbirds. In white-throated sparrows, VIP expression is higher in the anterior hypothalamus of birds with the supergene than those without it, and expression of VIP in this region predicts the level of territorial aggression regardless of genotype. Here, we aimed to identify epigenetic mechanisms that could contribute to differential expression of VIP both in breeding adults, which exhibit morph differences in territorial aggression, and in nestlings, before territorial behavior develops. We extracted and bisulfite-converted DNA from samples of the hypothalamus in wild-caught adults and nestlings and used high-throughput sequencing to measure DNA methylation of a region upstream of the VIP start site. We found that the allele inside the supergene was less methylated than the alternative allele in both adults and nestlings. The differential methylation was attributed primarily to CpG sites that were shared between the alleles, not to polymorphic sites, which suggests that epigenetic regulation is occurring independently of the genetic differentiation within the supergene. This work represents an initial step toward understanding how epigenetic differentiation inside chromosomal inversions leads to the development of alternative behavioral phenotypes.

Dynamic molecular evolution of a supergene with suppressed recombination in white-throated sparrows

In white-throated sparrows, two alternative morphs differing in plumage and behavior segregate with a large chromosomal rearrangement. As with sex chromosomes such as the mammalian Y, the rearranged version of chromosome two (ZAL2^m) is in a near-constant state of heterozygosity, offering opportunities to investigate both degenerative and selective processes during the early evolutionary stages of 'supergenes.' Here, we generated, synthesized, and analyzed extensive genome-scale data to better understand the forces shaping the evolution of the ZAL2 and ZAL2^m chromosomes in this species. We found that features of ZAL2^m are consistent with substantially reduced recombination and low levels of degeneration. We also found evidence that selective sweeps took place both on ZAL2^m and its standard counterpart, ZAL2, after the rearrangement event. Signatures of positive selection were associated with allelic bias in gene expression, suggesting that antagonistic selection has operated on gene regulation. Finally, we discovered a region exhibiting long-range haplotypes inside the rearrangement on ZAL2^m. These haplotypes appear to have been maintained by balancing selection, retaining genetic diversity within the supergene. Together, our analyses illuminate mechanisms contributing to the evolution of a young chromosomal polymorphism, revealing complex selective processes acting concurrently with genetic degeneration to drive the evolution of supergenes.

Supergenes on steroids

At the birth of supergenes, the genomic landscape is dramatically re-organized leading to pronounced differences in phenotypes and increased intrasexual diversity. Two of the best-studied supergenes in vertebrates are arguably the inversion polymorphisms on chromosomes 2 and 11 in the white-throated sparrow (Zonotrichia albicollis) and the ruff (Calidris pugnax), respectively. In both species, regions of suppressed recombination determine plumage coloration and social behavioural phenotypes. Despite the apparent lack of gene overlap between these two supergenes, in both cases the alternative phenotypes seem to be driven largely by alterations in steroid hormone pathways. Here, we explore the interplay between genomic architecture and steroid-related genes. Due to the highly pleiotropic effects of steroid-related genes and their universal involvement in social behaviour and transcriptomic regulation, processes favouring their linkage are likely to have substantial effects on the evolution of behavioural phenotypes, individual fitness, and life-history strategies. We propose that inversion-related differentiation and regulatory changes in steroid-related genes lie at the core of phenotypic differentiation in both of these interesting species. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Functional genomics of supergene-controlled behavior in the white-throated sparrow

Supergenes are regions of suppressed recombination that may span hundreds of genes and can control variation in key ecological phenotypes. Since genetic analysis is made impossible by the absence of recombination between genes, it has been difficult to establish how individual genes within these regions contribute to supergene-controlled phenotypes. The white-throated sparrow is a classic example in which a supergene controls behavioral differences as well as distinct coloration that determines mate choice. A landmark study now demonstrates that differences between supergene variants in the promoter sequences of a hormone receptor gene change its expression and control changes in behavior. To unambiguously establish the link between genotype and phenotype, the authors used antisense oligonucleotides to alter the level of gene expression in a focal brain region targeted through a cannula. The study showcases a powerful approach to the functional genomic manipulation of a wild vertebrate species.

Sparrows and supergenes: Ecological epigenetics in action

Despite the promise of ecological epigenetics, there remain few cases that clearly link epigenetic variation in wild animal populations to evolutionary change. In this issue of Molecular Ecology, Sun et al. provide such an example in white-throated sparrows-a fascinating system in which a large chromosomal rearrangement generates a "supergene" polymorphism linked to plumage colour, aggression and parenting behaviour. By combining whole genome bisulphite sequencing with RNA-sequencing and chromatin accessibility data, they show that the two alleles of this chromosomal polymorphism also exhibit substantial differences in DNA methylation levels, with implications for gene expression and transposable element activity. Their results provide a compelling case study for how genetic and epigenetic evolution proceed in concert. They also demonstrate the importance of integrating multiple types of genomic information to understand how gene regulation evolves, providing a model for future work in nonmodel species.

Gene Expression Modification by an Autosomal Inversion Associated With Three Male Mating Morphs

Chromosomal inversions are structural rearrangements that frequently provide genomic substrate for phenotypic diversity. In the ruff Philomachus pugnax, three distinct male reproductive morphs (Independents, Satellites and Faeders) are genetically determined by a 4.5 Mb autosomal inversion. Here we test how this stable inversion polymorphism affects gene expression in males during the lekking season. Gene expression may be altered through disruptions at the breakpoints and the accumulation of mutations due to suppressed recombination. We used quantitative PCR to measure expression of 11 candidate inversion genes across three different tissues (liver, adrenal glands and gonads) and tested for allelic imbalance in four inversion genes across 12 males of all three morphs (8 Independents, 2 Satellites, 2 Faeders). We quantified transcripts of CENPN, an essential gene disrupted by the inversion at the proximal breakpoint, at different exons distributed near and across the breakpoint region. Consistent with dosage dependent gene expression for the breakpoint gene CENPN, we found that expression in Independents was broadly similar for transcripts segments from inside and outside the inversion regions, whereas for Satellites and Faeders, transcript segments outside of the inversion showed at least twofold higher expression than those spanning over the breakpoint. Within the inversion, observed expression differences for inversion males across all four genes with allele-specific primers were consistent with allelic imbalance. We further analyzed gonadal expression of two inversion genes, HSD17B2 and SDR42E1, along with 12 non-inversion genes related to steroid metabolism and signaling in 25 males (13 Independents, 7 Satellites, 5 Faeders). Although we did not find clear morph differentiation for many individual genes, all three morphs could be separated based on gene expression differences when using linear discriminant analysis (LDA), regardless of genomic location (i.e., inside or outside of the inversion). This was robust to the removal of genes with the highest loadings. Pairwise correlations in the expression of genes showed significant correlations for 9–18 pairs of genes within morphs. However, between morphs, we only found a single association between genes SDR42E1 and AROM for Independents and Satellites. Our results suggest complex and wide-ranging changes in gene expression caused by structural variants.

The Genomic Architecture and Evolutionary Fates of Supergenes

Supergenes are genomic regions containing sets of tightly linked loci that control multi-trait phenotypic polymorphisms under balancing selection. Recent advances in genomics have uncovered significant variation in both the genomic architecture as well as the mode of origin of supergenes across diverse organismal systems. Although the role of genomic architecture for the origin of supergenes has been much discussed, differences in the genomic architecture also subsequently affect the evolutionary trajectory of supergenes and the rate of degeneration of supergene haplotypes. In this review, we synthesize recent genomic work and historical models of supergene evolution, highlighting how the genomic architecture of supergenes affects their evolutionary fate. We discuss how recent findings on classic supergenes involved in governing ant colony social form, mimicry in butterflies, and heterostyly in flowering plants relate to theoretical expectations. Furthermore, we use forward simulations to demonstrate that differences in genomic architecture affect the degeneration of supergenes. Finally, we discuss implications of the evolution of supergene haplotypes for the long-term fate of balanced polymorphisms governed by supergenes.

Genome-wide variation in DNA methylation linked to developmental stage and chromosomal suppression of recombination in white-throated sparrows

Much of our knowledge on regulatory impacts of DNA methylation has come from laboratory‐bred model organisms, which may not exhibit the full extent of variation found in wild populations. Here, we investigated naturally‐occurring variation in DNA methylation in a wild avian species, the white‐throated sparrow (Zonotrichia albicollis). This species offers exceptional opportunities for studying the link between genetic differentiation and phenotypic traits because of a nonrecombining chromosome pair linked to both plumage and behavioural phenotypes. Using novel single‐nucleotide resolution methylation maps and gene expression data, we show that DNA methylation and the expression of DNA methyltransferases are significantly higher in adults than in nestlings. Genes for which DNA methylation varied between nestlings and adults were implicated in development and cell differentiation and were located throughout the genome. In contrast, differential methylation between plumage morphs was concentrated in the nonrecombining chromosome pair. Interestingly, a large number of CpGs on the nonrecombining chromosome, localized to transposable elements, have undergone dramatic loss of DNA methylation since the split of the ZAL2 and ZAL2^m chromosomes. Changes in methylation predicted changes in gene expression for both chromosomes. In summary, we demonstrate changes in genome‐wide DNA methylation that are associated with development and with specific functional categories of genes in white‐throated sparrows. Moreover, we observe substantial DNA methylation reprogramming associated with the suppression of recombination, with implications for genome integrity and gene expression divergence. These results offer an unprecedented view of ongoing epigenetic reprogramming in a wild population.

see also the Perspective by Jordan A. Anderson and Jenny Tung.

Behavioral Genetics: Dissecting a Supergene to Understand Behavior

New research shows how alternative 'supergene' alleles of Estrogen Receptor 1 are differentially expressed in specific brain nuclei causing aggressive behavior in the white-throated sparrow.

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.

A supergene-linked estrogen receptor drives alternative phenotypes in a polymorphic songbird

Behavioral evolution relies on genetic changes, yet few behaviors can be traced to specific genetic sequences in vertebrates. Here we provide experimental evidence showing that differentiation of a single gene has contributed to the evolution of divergent behavioral phenotypes in the white-throated sparrow, a common backyard songbird. In this species, a series of chromosomal inversions has formed a supergene that segregates with an aggressive phenotype. The supergene has captured ESR1, the gene that encodes estrogen receptor α (ERα); as a result, this gene is accumulating changes that now distinguish the supergene allele from the standard allele. Our results show that in birds of the more aggressive phenotype, ERα knockdown caused a phenotypic change to that of the less aggressive phenotype. We next showed that in a free-living population, aggression is predicted by allelic imbalance favoring the supergene allele. Finally, we identified cis-regulatory features, both genetic and epigenetic, that explain the allelic imbalance. This work provides a rare illustration of how genotypic divergence has led to behavioral phenotypic divergence in a vertebrate.

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.

Genomic architecture and evolutionary antagonism drive allelic expression bias in the social supergene of red fire ants

Supergene regions maintain alleles of multiple genes in tight linkage through suppressed recombination. Despite their importance in determining complex phenotypes, our empirical understanding of early supergene evolution is limited. Here we focus on the young ‘social’ supergene of fire ants, a powerful system for disentangling the effects of evolutionary antagonism and suppressed recombination. We hypothesize that gene degeneration and social antagonism shaped the evolution of the fire ant supergene, resulting in distinct patterns of gene expression. We test these ideas by identifying allelic differences between supergene variants, characterizing allelic expression across populations, castes and body parts, and contrasting allelic expression biases with differences in expression between social forms. We find strong signatures of gene degeneration and gene-specific dosage compensation. On this background, a small portion of the genes has the signature of adaptive responses to evolutionary antagonism between social forms.

Red fire ants (Solenopsis invicta) are native to South America, but the species has spread to North America, Australia and New Zealand where it can be an invasive pest. A reason for this species’ invasiveness types of colonies : one with a single egg-laying queen and another with several queens. However, it is not possible to simply add more queens to a colony with one queen. Instead, the number of queens in a colony is controlled genetically, by a chromosome known as the ‘social chromosome’.

Like many other animals, red fire ants are diploid: their cells have two copies of each chromosome, which can carry two different versions of each gene. The social chromosome is no different, and it comes in two variants, SB and Sb. Each ant can therefore have either two SB chromosomes, leading to a colony with a single queen; or one SB chromosome and one Sb chromosome, leading to a colony with multiple queens. Ants with two copies of the Sb variant die when they are young, so the Sb version is inherited in a similar way to how the Y chromosome is passed on in humans. However, the social chromosome in red fire ants appeared less than one million years ago, making it much younger than the human Y chromosome, which is 180 million years old. This makes the social chromosome a good candidate for examining the early evolution of special chromosome variants that are only inherited.

How differences between the SB and the Sb chromosomes are evolving is an open question, however. Perhaps each version of the social chromosome has been optimised through natural selection to one colony type. Another suggestion is that the Sb chromosome has degenerated over time because its genes cannot be ‘reshuffled’ as they would be on normal chromosomes.

Martinez-Ruiz et al. compared genetic variants on the SB and Sb chromosomes, along with their expression in different types of ant colonies. The analysis showed that the Sb variant is in fact breaking down because of the lack of gene shuffling. This loss is compensated by intact copies of the same genes found on the SB variant, which explains why ants with the Sb variant can only survive if they also carry the SB version. Only a handful of genes on the social chromosomes appear to have been optimised by natural selection. Therefore Martinez-Ruiz et al. concluded the differences between the two chromosomes that lead to different colony types are collateral effects of Sb’s inability to reshuffle its genes.

This work reveals how a special chromosome similar to the Y chromosome in humans evolved. It also shows how multiple complex evolutionary forces can shape a species’ genetic makeup and social forms.

Vasoactive intestinal peptide as a mediator of the effects of a supergene on social behaviour

Supergenes, or linked groups of alleles that are inherited together, present excellent opportunities to understand gene-behaviour relationships. In white-throated sparrows (Zonotrichia albicollis), a supergene on the second chromosome associates with a more aggressive and less parental phenotype. This supergene includes the gene for vasoactive intestinal peptide (VIP), a neuropeptide known to play a causal role in both aggression and parental behaviour. Here, using a free-living population, we compared the levels of VIP mRNA between birds with and without the supergene. We focused on the anterior hypothalamus and infundibular region, two brain regions containing VIP neurons known to play a causal role in aggression and parental behaviour, respectively. First, we show that the supergene enhances VIP expression in the anterior hypothalamus and that expression positively predicts vocal aggression independently of genotype in both sexes. Next, we show that the supergene reduces VIP expression in the infundibular region, which suggests reduced secretion of prolactin, a pro-parental hormone. Thus, the patterns of VIP expression in these two regions are consistent with the enhanced aggression and reduced parental behaviour of birds with the supergene allele. Our results illustrate mechanisms by which elements of genomic architecture, such as supergenes, can contribute to the evolution of alternative behavioural phenotypes.

Regional epigenetic differentiation of the Z Chromosome between sexes in a female heterogametic system

In male heterogametic systems, the X Chromosome is epigenetically differentiated between males and females, to facilitate dosage compensation. For example, the X Chromosome in female mammals is largely inactivated. Relative to well-studied male heterogametic systems, the extent of epigenetic differentiation between male and female Z Chromosomes in female heterogametic species, which often lack complete dosage compensation, is poorly understood. Here, we examined the chromosomal DNA methylation landscapes of male and female Z Chromosomes in two distantly related avian species, namely chicken and white-throated sparrow. We show that, in contrast to the pattern in mammals, male and female Z Chromosomes in these species exhibit highly similar patterns of DNA methylation, which is consistent with weak or absent dosage compensation. We further demonstrate that the epigenetic differences between male and female chicken Z Chromosomes are localized to a few regions, including a previously identified male hypermethylated region 1 (MHM1; CGNC: 80601). We discovered a novel region with elevated male-to-female methylation ratios on the chicken Z Chromosome (male hypermethylated region 2 [MHM2]; CGNC: 80602). The MHM1 and MHM2, despite little sequence similarity between them, bear similar molecular features that are likely associated with their functions. We present evidence consistent with female hypomethylation of MHMs and up-regulation of nearby genes. Therefore, despite little methylation differentiation between sexes, extremely localized DNA methylation differences between male and female chicken Z Chromosomes have evolved and affect expression of nearby regions. Our findings offer new insights into epigenetic regulation of gene expression between sexes in female heterogametic systems.

A chromosomal inversion predicts the expression of sex steroid-related genes in a species with alternative behavioral phenotypes

In white-throated sparrows, a chromosomal rearrangement has led to alternative phenotypes that differ in sex steroid-dependent behaviors. The rearrangement has captured the genes estrogen receptor alpha and 5-alpha reductase, making these genes strong candidates for mediating the behavioral phenotypes. We report here that of the two genes, expression of estrogen receptor alpha mRNA differs between the morphs and predicts behavior to a much greater extent than does expression of 5-alpha reductase mRNA. Differentiation of estrogen receptor alpha, therefore, is likely more important for the behavioral phenotypes. We also found that in some brain regions, the degree to which testosterone treatment affects the expression of steroid-related genes depends strongly on morph. A large morph difference in estrogen receptor alpha mRNA expression in the amygdala appears to be independent of plasma testosterone; this difference persists during the non-breeding season and is detectable in nestlings at post-hatch day seven. The latter result suggests a substrate for organizational effects of hormones during development.

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.

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.

Comparative analyses of sperm DNA methylomes among human, mouse and cattle provide insights into epigenomic evolution and complex traits

Sperm DNA methylation is crucial for fertility and viability of offspring but epigenome evolution in mammals is largely understudied. By comparing sperm DNA methylomes and large-scale genome-wide association study (GWAS) signals between human and cattle, we aimed to examine the DNA methylome evolution and its associations with complex phenotypes in mammals. Our analysis revealed that genes with conserved non-methylated promoters (e.g., ANKS1A and WNT7A) among human and cattle were involved in common system and embryo development, and enriched for GWAS signals of body conformation traits in both species, while genes with conserved hypermethylated promoters (e.g., TCAP and CD80) were engaged in immune responses and highlighted by immune-related traits. On the other hand, genes with human-specific hypomethylated promoters (e.g., FOXP2 and HYDIN) were engaged in neuron system development and enriched for GWAS signals of brain-related traits, while genes with cattle-specific hypomethylated promoters (e.g., LDHB and DGAT2) mainly participated in lipid storage and metabolism. We validated our findings using sperm-retained nucleosome, preimplantation transcriptome, and adult tissue transcriptome data, as well as sequence evolutionary features, including motif binding sites, mutation rates, recombination rates and evolution signatures. In conclusion, our results demonstrate important roles of epigenome evolution in shaping the genetic architecture underlying complex phenotypes, hence enhance signal prioritization in GWAS and provide valuable information for human neurological disorders and livestock genetic improvement.