Detection of Allelic Frequency Differences between the Sexes in Humans: A Signature of Sexually Antagonistic Selection

The maintenance of genetic polymorphism underlying sexually antagonistic traits

Selection often favors different trait values in males and females, leading to genetic conflicts between the sexes when traits have a shared genetic basis. Such sexual antagonism has been proposed to maintain genetic polymorphism. However, this notion is based on insights from population genetic models of single loci with fixed fitness effects. It is thus unclear how readily polymorphism emerges from sex-specific selection acting on continuous traits, where fitness effects arise from the genotype-phenotype map and the fitness landscape. Here, we model the evolution of a continuous trait that has a shared genetic basis but different optima in males and females, considering a wide variety of genetic architectures and fitness landscapes. For autosomal loci, the long-term maintenance of polymorphism requires strong conflict between males and females that generates uncharacteristic sex-specific fitness patterns. Instead, more plausible sex-specific fitness landscapes typically generate stabilizing selection leading to an evolutionarily stable state that consists of a single homozygous genotype. Except for sites tightly linked to the sex-determining region, our results indicate that genetic variation due to sexual antagonism should arise only rarely and often be transient, making these signatures challenging to detect in genomic data.

Evolutionary divergence between homologous X-Y chromosome genes shapes sex-biased biology

Sex chromosomes are a fundamental aspect of sex-biased biology, but the extent to which homologous X-Y gene pairs ('the gametologs') contribute to sex-biased phenotypes remains hotly debated. Although these genes tend to exhibit large sex differences in expression throughout the body (XX females can express both X members, and XY males can express one X and one Y member), there is conflicting evidence regarding the degree of functional divergence between the X and Y members. Here we develop and apply co-expression fingerprint analysis to characterize functional divergence between the X and Y members of 17 gametolog gene pairs across >40 human tissues. Gametolog pairs exhibit functional divergence between the sexes that is driven by divergence between the X versus Y members (assayed in males), and this within-pair divergence is greatest among pairs with evolutionarily distant X and Y members. These patterns reflect that X versus Y gametologs show coordinated patterns of asymmetric coupling with large sets of autosomal genes, which are enriched for functional pathways and gene sets implicated in sex-biased biology and disease. Our findings suggest that the X versus Y gametologs have diverged in function and prioritize specific gametolog pairs for future targeted experimental studies.

The battle of the sexes in humans is highly polygenic

Sex-differential selection (SDS), which occurs when the fitness effects of alleles differ between males and females, can have profound impacts on the maintenance of genetic variation, disease risk, and other key aspects of natural populations. Because the sexes mix their autosomal genomes each generation, quantifying SDS is not possible using conventional population genetic approaches. Here, we introduce a method that exploits subtle sex differences in haplotype frequencies resulting from SDS acting in the current generation. Using data from 300K individuals in the UK Biobank, we estimate the strength of SDS throughout the genome. While only a handful of loci under SDS are individually significant, we uncover highly polygenic signals of genome-wide SDS for both viability and fecundity. Selection coefficients of [Formula: see text] may be typical. Despite its ubiquity, SDS may impose a mortality load of less than 1%. An interesting life-history tradeoff emerges: Alleles that increase viability more strongly in females than males tend to increase fecundity more strongly in males than in females. Finally, we find marginal evidence of SDS on fecundity acting on alleles affecting arm fat-free mass. Taken together, our findings connect the long-standing evidence of SDS acting on human phenotypes with its impact on the genome.

Sexually discordant selection is associated with trait-specific morphological changes and a complex genomic response

Sexes often have differing fitness optima, potentially generating intra-locus sexual conflict, as each sex bears a genetic "load" of alleles beneficial to the other sex. One strategy to evaluate conflict in the genome is to artificially select populations discordantly against established sexual dimorphism (SD), reintroducing attenuated conflict. We investigate a long-term artificial selection experiment reversing sexual size dimorphism in Drosophila melanogaster during ~350 generations of sexually discordant selection. We explore morphological and genomic changes to identify loci under selection between the sexes in discordantly and concordantly size-selected treatments. Despite substantial changes to overall size, concordant selection maintained ancestral SD. However, discordant selection altered size dimorphism in a trait-specific manner. We observe multiple possible soft selective sweeps in the genome, with size-related genes showing signs of selection. Patterns of genomic differentiation between the sexes within lineages identified potential sites maintained by sexual conflict. One discordant selected lineage shows a pattern of elevated genomic differentiation between males and females on chromosome 3L, consistent with the maintenance of sexual conflict. Our results suggest visible signs of conflict and differentially segregating alleles between the sexes due to discordant selection.

On the resolution of sexual conflict over shared traits

Anisogamy, different-sized male and female gametes, sits at the heart of sexual selection and conflict between the sexes. Sperm producers (males) and egg producers (females) of the same species generally share most, if not all, of the same genome, but selection frequently favours different trait values in each sex for traits common to both. The extent to which this conflict might be resolved, and the potential mechanisms by which this can occur, have been widely debated. Here, we summarize recent findings and emphasize that once the sexes evolve, sexual selection is ongoing, and therefore new conflict is always possible. In addition, sexual conflict is largely a multivariate problem, involving trait combinations underpinned by networks of interconnected genes. Although these complexities can hinder conflict resolution, they also provide multiple possible routes to decouple male and female phenotypes and permit sex-specific evolution. Finally, we highlight difficulty in the study of sexual conflict over shared traits and promising directions for future research.

Unveiling recent and ongoing adaptive selection in human populations

Genome-wide scans for signals of selection have become a routine part of the analysis of population genomic variation datasets and have resulted in compelling evidence of selection during recent human evolution. This Essay spotlights methodological innovations that have enabled the detection of selection over very recent timescales, even in contemporary human populations. By harnessing large-scale genomic and phenotypic datasets, these new methods use different strategies to uncover connections between genotype, phenotype, and fitness. This Essay outlines the rationale and key findings of each strategy, discusses challenges in interpretation, and describes opportunities to improve detection and understanding of ongoing selection in human populations.

The roles of sexual selection and sexual conflict in shaping patterns of genome and transcriptome variation

Sexual dimorphism is one of the most prevalent, and often the most extreme, examples of phenotypic variation within species, and arises primarily from genomic variation that is shared between females and males. Many sexual dimorphisms arise through sex differences in gene expression, and sex-biased expression is one way that a single, shared genome can generate multiple, distinct phenotypes. Although many sexual dimorphisms are expected to result from sexual selection, and many studies have invoked the possible role of sexual selection to explain sex-specific traits, the role of sexual selection in the evolution of sexually dimorphic gene expression remains difficult to differentiate from other forms of sex-specific selection. In this Review, we propose a holistic framework for the study of sex-specific selection and transcriptome evolution. We advocate for a comparative approach, across tissues, developmental stages and species, which incorporates an understanding of the molecular mechanisms, including genomic variation and structure, governing gene expression. Such an approach is expected to yield substantial insights into the evolution of genetic variation and have important applications in a variety of fields, including ecology, evolution and behaviour.

Evolutionary dynamics of pseudoautosomal region 1 in humans and great apes

BACKGROUND

The pseudoautosomal region 1 (PAR1) is a 2.7 Mb telomeric region of human sex chromosomes. PAR1 has a crucial role in ensuring proper segregation of sex chromosomes during male meiosis, exposing it to extreme recombination and mutation processes. We investigate PAR1 evolution using population genomic datasets of extant humans, eight populations of great apes, and two archaic human genome sequences.

RESULTS

We find that PAR1 is fast evolving and closer to evolutionary nucleotide equilibrium than autosomal telomeres. We detect a difference between substitution patterns and extant diversity in PAR1, mainly driven by the conflict between strong mutation and recombination-associated fixation bias at CpG sites. We detect excess C-to-G mutations in PAR1 of all great apes, specific to the mutagenic effect of male recombination. Despite recent evidence for Y chromosome introgression from humans into Neanderthals, we find that the Neanderthal PAR1 retained similarity to the Denisovan sequence. We find differences between substitution spectra of these archaics suggesting rapid evolution of PAR1 in recent hominin history. Frequency analysis of alleles segregating in females and males provided no evidence for recent sexual antagonism in this region. We study repeat content and double-strand break hotspot regions in PAR1 and find that they may play roles in ensuring the obligate X-Y recombination event during male meiosis.

CONCLUSIONS

Our study provides an unprecedented quantification of population genetic forces governing PAR1 biology across extant and extinct hominids. PAR1 evolutionary dynamics are predominantly governed by recombination processes with a strong impact on mutation patterns across all species.

Polygenic signals of sex differences in selection in humans from the UK Biobank

Sex differences in the fitness effects of genetic variants can influence the rate of adaptation and the maintenance of genetic variation. For example, "sexually antagonistic" (SA) variants, which are beneficial for one sex and harmful for the other, can both constrain adaptation and increase genetic variability for fitness components such as survival, fertility, and disease susceptibility. However, detecting variants with sex-differential fitness effects is difficult, requiring genome sequences and fitness measurements from large numbers of individuals. Here, we develop new theory for studying sex-differential selection across a complete life cycle and test our models with genotypic and reproductive success data from approximately 250,000 UK Biobank individuals. We uncover polygenic signals of sex-differential selection affecting survival, reproductive success, and overall fitness, with signals of sex-differential reproductive selection reflecting a combination of SA polymorphisms and sexually concordant polymorphisms in which the strength of selection differs between the sexes. Moreover, these signals hold up to rigorous controls that minimise the contributions of potential confounders, including sequence mapping errors, population structure, and ascertainment bias. Functional analyses reveal that sex-differentiated sites are enriched in phenotype-altering genomic regions, including coding regions and loci affecting a range of quantitative traits. Population genetic analyses show that sex-differentiated sites exhibit evolutionary histories dominated by genetic drift and/or transient balancing selection, but not long-term balancing selection, which is consistent with theoretical predictions of effectively weak SA balancing selection in historically small populations. Overall, our results are consistent with polygenic sex-differential-including SA-selection in humans. Evidence for sex-differential selection is particularly strong for variants affecting reproductive success, in which the potential contributions of nonrandom sampling to signals of sex differentiation can be excluded.

Genome-wide interaction analysis identified low-frequency variants with sex disparity in lung cancer risk

Differences by sex in lung cancer incidence and mortality have been reported which cannot be fully explained by sex differences in smoking behavior, implying existence of genetic and molecular basis for sex disparity in lung cancer development. However, the information about sex dimorphism in lung cancer risk is quite limited despite the great success in lung cancer association studies. By adopting a stringent two-stage analysis strategy, we performed a genome-wide gene-sex interaction analysis using genotypes from a lung cancer cohort including ~ 47 000 individuals with European ancestry. Three low-frequency variants (minor allele frequency < 0.05), rs17662871 [odds ratio (OR) = 0.71, P = 4.29×10-8); rs79942605 (OR = 2.17, P = 2.81×10-8) and rs208908 (OR = 0.70, P = 4.54×10-8) were identified with different risk effect of lung cancer between men and women. Further expression quantitative trait loci and functional annotation analysis suggested rs208908 affects lung cancer risk through differential regulation of Coxsackie virus and adenovirus receptor gene expression in lung tissues between men and women. Our study is one of the first studies to provide novel insights about the genetic and molecular basis for sex disparity in lung cancer development.

Genetic and genomic analysis of oxygen consumption in mice

We estimated genetic parameters for oxygen consumption (OC), OC per metabolic body weight (OCMBW) and body weight at three through 8 weeks of age in divergently selected mice populations, with an animal model considering maternal genetic, common litter environmental and cytoplasmic inheritance effects. Cytoplasmic inheritance was considered based on maternal lineage information. With respect to OC, estimated direct heritability was moderate (0.32) and the estimated proportion of the variance of cytoplasmic inheritance effects to the phenotypic variance was very low (0.01), implying that causal genes for OC could be located on autosomes. To assess this hypothesis, we attempted to identify possible candidate causal genes through selective signature detection with the results of pooled whole-genome resequencing using pooled DNA samples from high and low OC mice. We made a list of possible candidate causal genes for OC, including those relating to electron transport chain and ATP-binding proteins (Ndufa12, Sdhc, Atp10b, etc.), Prr16 encoding Largen protein, Cry1 encoding a key component of the circadian core oscillator and so on. The results, although careful interpretation must be required, could contribute to elucidate the genetic mechanism of OC, an indicator for maintenance energy requirement, and therefore feed efficiency.

Detection of sexually antagonistic transmission distortions in trio datasets

Sexual dimorphisms are widespread in animals and plants, for morphological as well as physiological traits. Understanding the genetic basis of sexual dimorphism and its evolution is crucial for understanding biological differences between the sexes. Genetic variants with sex‐antagonistic effects on fitness are expected to segregate in populations at the early phases of sexual dimorphism emergence. Detecting such variants is notoriously difficult, and the few genome‐scan methods employed so far have limited power and little specificity. Here, we propose a new framework to detect a signature of sexually antagonistic (SA) selection. We rely on trio datasets where sex‐biased transmission distortions can be directly tracked from parents to offspring, and identify signals of SA transmission distortions in genomic regions. We report the genomic location of six candidate regions detected in human populations as potentially under sexually antagonist selection. We find an enrichment of genes associated with embryonic development within these regions. Last, we highlight two candidate regions for SA selection in humans.

An unbiased test reveals no enrichment of sexually antagonistic polymorphisms on the human X chromosome

Mutations with beneficial effects in one sex can have deleterious effects in the other. Such 'sexually antagonistic' (SA) variants contribute to variation in life-history traits and overall fitness, yet their genomic distribution is poorly resolved. Theory predicts that SA variants could be enriched on the X chromosome or autosomes, yet current empirical tests face two formidable challenges: (i) identifying SA selection in genomic data is difficult; and (ii) metrics of SA variation show persistent biases towards the X, even when SA variants are randomly distributed across the genome. Here, we present an unbiased test of the theory that SA variants are enriched on the X. We first develop models for reproductive FST-a metric for quantifying sex-differential (including SA) effects of genetic variants on lifetime reproductive success-that control for X-linked biases. Comparing data from approximately 250 000 UK Biobank individuals to our models, we find FST elevations consistent with both X-linked and autosomal SA polymorphisms affecting reproductive success in humans. However, the extent of FST elevations does not differ from a model in which SA polymorphisms are randomly distributed across the genome. We argue that the polygenic nature of SA variation, along with sex asymmetries in SA effects, might render X-linked enrichment of SA polymorphisms unlikely.

Gene duplication to the Y chromosome in Trinidadian Guppies

Differences in allele frequencies at autosomal genes between males and females in a population can result from two scenarios. First, unresolved sexual conflict over survival can produce allelic differentiation between the sexes. However, given the substantial mortality costs required to produce allelic differences between males and females at each generation, it remains unclear how many loci within the genome experience significant sexual conflict over survival. Alternatively, recent studies have shown that similarity between autosomal and Y sequences can create perceived allelic differences between the sexes. However, Y duplications are most likely in species with large non-recombining regions, in part because they simply represent larger targets for duplications. We assessed the genomes of 120 wild-caught guppies, which experience extensive predation- and pathogen-induced mortality and have a relatively small ancestral Y chromosome. We identified seven autosomal genes that show allelic differences between male and female adults. Five of these genes show clear evidence of whole or partial gene duplication between the Y chromosome and the autosomes. The remaining two genes show evidence of partial homology to the Y. Overall, our findings suggest that the guppy genome experiences a very low level of unresolved sexual conflict over survival, and instead the Y chromosome, despite its small ancestral size and recent origin, may nonetheless accumulate genes with male-specific functions.

Experimental evolution reveals sex-specific dominance for surviving bacterial infection in laboratory populations of Drosophila melanogaster

Males and females are subjected to distinct kinds of selection pressures, often leading to the evolution of sex‐specific genetic architecture, an example being sex‐specific dominance. Sex‐specific dominance reversals (SSDRs), where alleles at sexually antagonistic loci are at least partially dominant in the sex they benefit, have been documented in Atlantic salmon, rainbow trout, and seed beetles. Another interesting feature of many sexually reproducing organisms is the asymmetric inheritance pattern of X chromosomes, which often leads to distinct evolutionary outcomes on X chromosomes compared to autosomes. Examples include the higher efficacy of sexually concordant selection on X chromosomes, and X chromosomes being more conducive to the maintenance of sexually antagonistic polymorphisms under certain conditions. Immunocompetence is a trait that has been extensively investigated for sexual dimorphism with growing evidence for sex‐specific or sexually antagonistic variation. X chromosomes have been shown to harbor substantial immunity‐related genetic variation in the fruit fly, Drosophila melanogaster. Here, using interpopulation crosses and cytogenetic cloning, we investigated sex‐specific dominance and the role of the X chromosome in improved postinfection survivorship of laboratory populations of D. melanogaster selected against pathogenic challenge by Pseudomonas entomophila. We could not detect any contribution of the X chromosome to the evolved immunocompetence of our selected populations, as well as to within‐population variation in immunocompetence. However, we found strong evidence of sex‐specific dominance related to surviving bacterial infection. Our results indicate that alleles that confer a survival advantage to the selected populations are, on average, partially dominant in females but partially recessive in males. This could also imply an SSDR for overall fitness, given the putative evidence for sexually antagonistic selection affecting immunocompetence in Drosophila melanogaster. We also highlight sex‐specific dominance as a potential mechanism of sex differences in immunocompetence, with population‐level sex differences primarily driven by sex differences in heterozygotes.

Diversity of Modes of Reproduction and Sex Determination Systems in Invertebrates, and the Putative Contribution of Genetic Conflict

About eight million animal species are estimated to live on Earth, and all except those belonging to one subphylum are invertebrates. Invertebrates are incredibly diverse in their morphologies, life histories, and in the range of the ecological niches that they occupy. A great variety of modes of reproduction and sex determination systems is also observed among them, and their mosaic-distribution across the phylogeny shows that transitions between them occur frequently and rapidly. Genetic conflict in its various forms is a long-standing theory to explain what drives those evolutionary transitions. Here, we review (1) the different modes of reproduction among invertebrate species, highlighting sexual reproduction as the probable ancestral state; (2) the paradoxical diversity of sex determination systems; (3) the different types of genetic conflicts that could drive the evolution of such different systems.

Identification of genomic regions that exhibit sexual dimorphism for size and muscularity in cattle

Sexual dimorphism, the phenomenon whereby males and females of the same species are distinctive in some aspect of appearance or size, has previously been documented in cattle for traits such as growth rate and carcass merit using a quantitative genetics approach. No previous study in cattle has attempted to document sexual dimorphism at a genome level; therefore, the objective of the present study was to determine whether genomic regions associated with size and muscularity in cattle exhibited signs of sexual dimorphism. Analyses were undertaken on 10 linear-type traits that describe the muscular and skeletal characteristics of both males and females of five beef cattle breeds: 1,444 Angus (AA), 6,433 Charolais (CH), 1,129 Hereford, 8,745 Limousin (LM), and 1,698 Simmental. Genome-wide association analyses were undertaken using imputed whole-genome sequence data for each sex separately by breed. For each single-nucleotide polymorphism (SNP) that was segregating in both sexes, the difference between the allele substitution effect sizes for each sex, in each breed separately, was calculated. Suggestively (P ≤ 1 × 10-5) sexually dimorphic SNPs that were segregating in both males and females were detected for all traits in all breeds, although the location of these SNPs differed by both trait and breed. Significantly (P ≤ 1 × 10-8) dimorphic SNPs were detected in just three traits in the AA, seven traits in the CH, and three traits in the LM. The vast majority of all segregating autosomal SNPs (86% in AA to 94% in LM) had the same minor allele in both males and females. Differences (P ≤ 0.05) in allele frequencies between the sexes were observed for between 36% (LM) and 66% (AA) of the total autosomal SNPs that were segregating in both sexes. Dimorphic SNPs were located within a number of genes related to muscularity and/or size including the NAB1, COL5A2, and IWS1 genes on BTA2 that are located close to, and thought to be co-inherited with, the MSTN gene. Overall, sexual dimorphism exists in cattle at the genome level, but it is not consistent by either trait or breed.

An intronic DHCR7 genetic polymorphism associates with vitamin D serum level and incidence of acute coronary syndrome

INTRODUCTION

Vitamin D deficiency has been linked to cardiovascular pathologies including acute coronary syndrome (ACS). Polymorphisms in vitamin D associated genes have been confounding to vitamin D serum levels and pathological predispositions. 7-hydrocholesterol is a common precursor in cholesterol and vitamin D synthesis. DHCR7/NADSYN1 genetic locus expresses 7-hydrocholesterol reductase (DHCR7), an enzyme that recruits 7-hydrocholesterol in cholesterol biosynthesis, and NAD synthetase 1 (NADSYN1), which participates in the hydroxylation of 25 hydroxyvitamin D.

AIM

This study aims to correlate two polymorphisms in the DHCR7/NADSYN1 genetic locus with levels of circulatory vitamin D and the presentation of ACS in an Egyptian population.

METHODS

In a case control study, 189 ACS patients and 106 healthy control subjects were genotyped for SNPs rs11606033 of the DHCR7 gene and rs2276360 of the NADSYN1 gene using the amplification-refractory mutation system (ARMS). The levels of 25(OH)D2 and 25(OH)D3 were measured using an in-house developed and validated ultra-performance liquid chromatography tandem mass spectrometry (UPLC-MS/MS) based protocol.

RESULTS

ACS patients have significantly lower levels of circulating vitamin D in comparison to healthy controls. Allele A of the DHCR7 polymorphism was found to correlate with serum vitamin D deficiency and incidence of ACS classes: NSTEMI, STEMI and unstable angina, when compared to allele G. On the other hand, the NADSYN1 polymorphism rs2276360 correlated with serum 25(OH)D3 deficiency. Yet, no significant correlation was found with incidences of ACS.

CONCLUSION

We conclude that rs11606033, which is an intronic SNP between exon 4 and exon 5 of the DHCR7 gene, influences vitamin D serum abundance and more importantly ACS incidence.

Dispersal Alters the Nature and Scope of Sexually Antagonistic Variation

AbstractIntralocus sexual conflict, or sexual antagonism, occurs when alleles have opposing fitness effects in the two sexes. Previous theory suggests that sexual antagonism is a driver of genetic variation by generating balancing selection. However, most of these studies assume that populations are well mixed, neglecting the effects of spatial subdivision. Here, we use mathematical modeling to show that limited dispersal changes evolution at sexually antagonistic autosomal and X-linked loci as a result of inbreeding and sex-specific kin competition. We find that if the sexes disperse at different rates, kin competition within the philopatric sex biases intralocus conflict in favor of the more dispersive sex. Furthermore, kin competition diminishes the strength of balancing selection relative to genetic drift, reducing genetic variation in small subdivided populations. Meanwhile, by decreasing heterozygosity, inbreeding reduces the scope for sexually antagonistic polymorphism due to nonadditive allelic effects, and this occurs to a greater extent on the X chromosome than autosomes. Overall, our results indicate that spatial structure is a relevant factor in predicting where sexually antagonistic alleles might be observed. We suggest that sex-specific dispersal ecology and demography can contribute to interspecific and intragenomic variation in sexual antagonism.

Evolutionary dynamics of the human pseudoautosomal regions

Recombination between the X and Y human sex chromosomes is limited to the two pseudoautosomal regions (PARs) that present quite distinct evolutionary origins. Despite the crucial importance for male meiosis, genetic diversity patterns and evolutionary dynamics of these regions are poorly understood. In the present study, we analyzed and compared the genetic diversity of the PAR regions using publicly available genomic sequences encompassing both PAR1 and PAR2. Comparisons were performed through allele diversities, linkage disequilibrium status and recombination frequencies within and between X and Y chromosomes. In agreement with previous studies, we confirmed the role of PAR1 as a male-specific recombination hotspot, but also observed similar characteristic patterns of diversity in both regions although male recombination occurs at PAR2 to a much lower extent (at least one recombination event at PAR1 and in ≈1% in normal male meioses at PAR2). Furthermore, we demonstrate that both PARs harbor significantly different allele frequencies between X and Y chromosomes, which could support that recombination is not sufficient to homogenize the pseudoautosomal gene pool or is counterbalanced by other evolutionary forces. Nevertheless, the observed patterns of diversity are not entirely explainable by sexually antagonistic selection. A better understanding of such processes requires new data from intergenerational transmission studies of PARs, which would be decisive on the elucidation of PARs evolution and their role in male-driven heterosomal aneuploidies.

Evaluating human autosomal loci for sexually antagonistic viability selection in two large biobanks

Sex and sexual differentiation are pervasive across the tree of life. Because females and males often have substantially different functional requirements, we expect selection to differ between the sexes. Recent studies in diverse species, including humans, suggest that sexually antagonistic viability selection creates allele frequency differences between the sexes at many different loci. However, theory and population-level simulations indicate that sex-specific differences in viability would need to be very large to produce and maintain reported levels of between-sex allelic differentiation. We address this contradiction between theoretical predictions and empirical observations by evaluating evidence for sexually antagonistic viability selection on autosomal loci in humans using the largest cohort to date (UK Biobank, n = 487,999) along with a second large, independent cohort (BioVU, n = 93,864). We performed association tests between genetically ascertained sex and autosomal loci. Although we found dozens of genome-wide significant associations, none replicated across cohorts. Moreover, closer inspection revealed that all associations are likely due to cross-hybridization with sex chromosome regions during genotyping. We report loci with potential for mis-hybridization found on commonly used genotyping platforms that should be carefully considered in future genetic studies of sex-specific differences. Despite being well powered to detect allele frequency differences of up to 0.8% between the sexes, we do not detect clear evidence for this signature of sexually antagonistic viability selection on autosomal variation. These findings suggest a lack of strong ongoing sexually antagonistic viability selection acting on single locus autosomal variation in humans.

Disagreement in FST estimators: A case study from sex chromosomes

Sewall Wright developed FST for describing population differentiation and it has since been extended to many novel applications, including the detection of homomorphic sex chromosomes. However, there has been confusion regarding the expected estimate of FST for a fixed difference between the X- and Y-chromosome when comparing males and females. Here, we attempt to resolve this confusion by contrasting two common FST estimators and explain why they yield different estimates when applied to the case of sex chromosomes. We show that this difference is true for many allele frequencies, but the situation characterized by fixed differences between the X- and Y-chromosome is among the most extreme. To avoid additional confusion, we recommend that all authors using FST clearly state which estimator of FST their work uses.

Is the X chromosome a hot spot for sexually antagonistic polymorphisms? Biases in current empirical tests of classical theory

Females and males carry nearly identical genomes, which can constrain the evolution of sexual dimorphism and generate conditions that are favourable for maintaining sexually antagonistic (SA) polymorphisms, in which alleles beneficial for one sex are deleterious for the other. An influential theoretical prediction, by Rice (Rice 1984 Evolution38, 735-742), is that the X chromosome should be a 'hot spot' (i.e. enriched) for SA polymorphisms. While important caveats to Rice's theoretical prediction have since been highlighted (e.g. by Fry (2010) Evolution64, 1510-1516), several empirical studies appear to support it. Here, we show that current tests of Rice's theory-most of which are based on quantitative genetic measures of fitness (co)variance-are frequently biased towards detecting X-linked effects. We show that X-linked genes tend to contribute disproportionately to quantitative genetic patterns of SA fitness variation whether or not the X is enriched for SA polymorphisms. Population genomic approaches for detecting SA loci, including genome-wide association study of fitness and analyses of intersexual FST, are similarly biased towards detecting X-linked effects. In the light of our models, we critically re-evaluate empirical evidence for Rice's theory and discuss prospects for empirically testing it.

The search for sexually antagonistic genes: Practical insights from studies of local adaptation and statistical genomics

Sexually antagonistic (SA) genetic variation-in which alleles favored in one sex are disfavored in the other-is predicted to be common and has been documented in several animal and plant populations, yet we currently know little about its pervasiveness among species or its population genetic basis. Recent applications of genomics in studies of SA genetic variation have highlighted considerable methodological challenges to the identification and characterization of SA genes, raising questions about the feasibility of genomic approaches for inferring SA selection. The related fields of local adaptation and statistical genomics have previously dealt with similar challenges, and lessons from these disciplines can therefore help overcome current difficulties in applying genomics to study SA genetic variation. Here, we integrate theoretical and analytical concepts from local adaptation and statistical genomics research-including F ST and F IS statistics, genome-wide association studies, pedigree analyses, reciprocal transplant studies, and evolve-and-resequence experiments-to evaluate methods for identifying SA genes and genome-wide signals of SA genetic variation. We begin by developing theoretical models for between-sex F ST and F IS, including explicit null distributions for each statistic, and using them to critically evaluate putative multilocus signals of sex-specific selection in previously published datasets. We then highlight new statistics that address some of the limitations of F ST and F IS, along with applications of more direct approaches for characterizing SA genetic variation, which incorporate explicit fitness measurements. We finish by presenting practical guidelines for the validation and evolutionary analysis of candidate SA genes and discussing promising empirical systems for future work.

Gene-based analyses of the maternal genome implicate maternal effect genes as risk factors for conotruncal heart defects

Congenital heart defects (CHDs) affect approximately 1% of newborns. Epidemiological studies have identified several genetically-mediated maternal phenotypes (e.g., pregestational diabetes, chronic hypertension) that are associated with the risk of CHDs in offspring. However, the role of the maternal genome in determining CHD risk has not been defined. We present findings from gene-level, genome-wide studies that link CHDs to maternal effect genes as well as to maternal genes related to hypertension and proteostasis. Maternal effect genes, which provide the mRNAs and proteins in the oocyte that guide early embryonic development before zygotic gene activation, have not previously been implicated in CHD risk. Our findings support a role for and suggest new pathways by which the maternal genome may contribute to the development of CHDs in offspring.

Recessive Z-linked lethals and the retention of haplotype diversity in a captive butterfly population

Sex chromosomes are predicted to harbour elevated levels of sexually antagonistic variation due to asymmetries in the heritability of recessive traits in the homogametic versus heterogametic sex. This evolutionary dynamic may manifest as high recessive load specifically affecting the homogametic sex, and the retention of haplotype diversity in small populations. We tested the hypothesis that the Z chromosome in the butterfly Bicyclus anynana carries a high inbred load for male fertility and viability. Homozygosity of Z chromosome blocks was produced by daughter-father backcrosses, and inferred from marker loci positioned via a linkage map. Male sterility was, in general, unrelated to homozygosity in any region of the Z, but there was an extreme deficit of homozygous males within a 2 cM interval in all families. In contrast, no corresponding skew in Z genotype was detected in their (hemizygous) sisters. The same pattern was observed in historically inbred lines, indicating a high frequency of recessive lethals in the ancestral population. Allele-frequency changes between 1993 and 2006 (70 generations at N_e ~ 160) show that, despite the loss of many haplotypes, diversity was retained significantly above the neutral expectation. Effective overdominance in the lethal region can account for this effect locally but not in other parts of the chromosome, that are also associated with persistent linkage disequilibrium. These unexpected patterns suggest the operation of other factors, such as epistatic selection, recombination suppression, assortative mating and meiotic drive. Our results highlight the role of balancing selection in maintaining the inbred load and linked genetic diversity.

Signature of sexual conflict is actually conflict resolved

There has been substantial interest of late in using population genetic methods to study sexual conflict, where an allele increases the fitness of one sex at some cost to the other (Mank, 2017). Population genomic scans for sexual conflict offer an important advance given the difficulties of identifying antagonistic alleles from more traditional methods, and could greatly increase our understanding of the extent and loci of sexual conflict. This is particularly true for studies in natural populations, for which obtaining accurate fitness measurements for each sex can be challenging. In this issue of Molecular Ecology, Bissegger, Laurentino, Roesti, and Berner (2019) present a cautionary tale about how to interpret these population genomic data.

Sex matters: Otolith shape and genomic variation in deacon rockfish (Sebastes diaconus)

Little is known about intraspecific variation within the deacon rockfish (Sebastes diaconus), a recently described species found in the northeast Pacific Ocean. We investigated population structure among fish sampled from two nearshore reefs (Siletz Reef and Seal Rock) and one offshore site (Stonewall Bank) within a <50-km2 area off the Oregon coast. Fish from the three sample sites exhibited small but statistically significant differences based on genetic variation at >15,000 neutral loci, whether analyzed independently or classified into nearshore and offshore groups. Male and females were readily distinguished using genetic data and 92 outlier loci were associated with sex, potentially indicating differential selection between males and females. Morphometric results indicated that there was significant secondary sexual dimorphism in otolith shape, but further sampling is required to disentangle potential confounding influence of age. This study is the first step toward understanding intraspecific variation within the deacon rockfish and the potential management implications. Since differentiation among the three sample sites was small, we consider the results to be suggestive of a single stock. However, future studies should evaluate how the stock is affected by differences in sex, age, and gene flow between the nearshore and offshore environments.

The genomic footprint of sexual conflict

Genes with sex-biased expression show a number of unique properties and this has been seen as evidence for conflicting selection pressures in males and females, forming a genetic ‘tug-of-war’ between the sexes. However, we lack studies of taxa where an understanding of conflicting phenotypic selection in the sexes has been linked with studies of genomic signatures of sexual conflict. Here, we provide such a link. We used an insect where sexual conflict is unusually well understood, the seed beetle Callosobruchus maculatus, to test for molecular genetic signals of sexual conflict across genes with varying degrees of sex-bias in expression. We sequenced, assembled and annotated its genome and performed population resequencing of three divergent populations. Sex-biased genes showed increased levels of genetic diversity and bore a remarkably clear footprint of relaxed purifying selection. Yet, segregating genetic variation was also affected by balancing selection in weakly female-biased genes, while male-biased genes showed signs of overall purifying selection. Female-biased genes contributed disproportionally to shared polymorphism across populations, while male-biased genes, male seminal fluid protein genes and sex-linked genes did not. Genes showing genomic signatures consistent with sexual conflict generally matched life-history phenotypes known to experience sexually antagonistic selection in this species. Our results highlight metabolic and reproductive processes, confirming the key role of general life-history traits in sexual conflict.

Limits to Genomic Divergence Under Sexually Antagonistic Selection

Since the autosomal genome is shared between the sexes, sex-specific fitness optima present an evolutionary challenge. While sexually antagonistic selection might favor different alleles within females and males, segregation randomly reassorts alleles at autosomal loci between sexes each generation. This process of homogenization during transmission thus prevents between-sex allelic divergence generated by sexually antagonistic selection from accumulating across multiple generations. However, recent empirical studies have reported high male-female F_ST statistics. Here, we use a population genetic model to evaluate whether these observations could plausibly be produced by sexually antagonistic selection. To do this, we use both a single-locus model with nonrandom mate choice, and individual-based simulations to study the relationship between strength of selection, degree of between-sex divergence, and the associated genetic load. We show that selection must be exceptionally strong to create measurable divergence between the sexes and that the decrease in population fitness due to this process is correspondingly high. Individual-based simulations with selection genome-wide recapitulate these patterns and indicate that small sample sizes and sampling variance can easily generate substantial male-female divergence. We therefore conclude that caution should be taken when interpreting autosomal allelic differentiation between the sexes.

Phenotypic sexual dimorphism is associated with genomic signatures of resolved sexual conflict

Intralocus sexual conflict, where an allele benefits one sex at the expense of the other, has an important role in shaping genetic diversity of populations through balancing selection. However, the potential for mating systems to exert balancing selection through sexual conflict on the genome remains unclear. Furthermore, the nature and potential for resolution of sexual conflict across the genome has been hotly debated. To address this, we analysed de novo transcriptomes from six avian species, chosen to reflect the full range of sexual dimorphism and mating systems. Our analyses combine expression and population genomic statistics across reproductive and somatic tissue, with measures of sperm competition and promiscuity. Our results reveal that balancing selection is weakest in the gonad, consistent with the resolution of sexual conflict and evolutionary theory that phenotypic sex differences are associated with lower levels of ongoing conflict. We also demonstrate a clear link between variation in sexual conflict and levels of genetic variation across phylogenetic space in a comparative framework. Our observations suggest that this conflict is short-lived, and is resolved via the decoupling of male and female gene expression patterns, with important implications for the role of sexual selection in adaptive potential and role of dimorphism in facilitating sex-specific fitness optima.

Genome-wide sexually antagonistic variants reveal long-standing constraints on sexual dimorphism in fruit flies

The evolution of sexual dimorphism is constrained by a shared genome, leading to ‘sexual antagonism’, in which different alleles at given loci are favoured by selection in males and females. Despite its wide taxonomic incidence, we know little about the identity, genomic location, and evolutionary dynamics of antagonistic genetic variants. To address these deficits, we use sex-specific fitness data from 202 fully sequenced hemiclonal Drosophila melanogaster fly lines to perform a genome-wide association study (GWAS) of sexual antagonism. We identify approximately 230 chromosomal clusters of candidate antagonistic single nucleotide polymorphisms (SNPs). In contradiction to classic theory, we find no clear evidence that the X chromosome is a hot spot for sexually antagonistic variation. Characterising antagonistic SNPs functionally, we find a large excess of missense variants but little enrichment in terms of gene function. We also assess the evolutionary persistence of antagonistic variants by examining extant polymorphism in wild D. melanogaster populations and closely related species. Remarkably, antagonistic variants are associated with multiple signatures of balancing selection across the D. melanogaster distribution range and in their sister species D. simulans, indicating widespread and evolutionarily persistent (about 1 million years) genomic constraints on the evolution of sexual dimorphism. Based on our results, we propose that antagonistic variation accumulates because of constraints on the resolution of sexual conflict over protein coding sequences, thus contributing to the long-term maintenance of heritable fitness variation.

This study characterises antagonistic loci across the genome of the fruit fly Drosophila melanogaster, finding them to be preferentially associated with variation in coding sequences and to be selectively maintained across worldwide populations of D. melanogaster, and even its sister species D. simulans.

Sexual conflict and the Trivers-Willard hypothesis: Females prefer daughters and males prefer sons

Because parental care is expected to depend on the fitness returns generated by each unit of investment, it should be sensitive to both offspring condition and parental ability to invest. The Trivers-Willard Hypothesis (TWH) predicts that parents who are in good condition will bias investment towards sons, while parents who are in poor condition will bias investment towards daughters because high-quality sons are expected to out-reproduce high quality daughters, while low-quality daughters are expected to out-reproduce low quality sons. We report results from an online experiment testing the Trivers-Willard effect by measuring implicit and explicit psychological preferences and behaviorally implied preferences for sons or daughters both as a function of their social and economic status and in the aftermath of a priming task designed to make participants feel wealthy or poor. We find only limited support for predictions derived from the TWH and instead find that women have strong preferences for girls and men have preferences for boys.

The Evolutionary Consequences of Selection at the Haploid Gametic Stage

As an immediate consequence of sexual reproduction, biphasic life cycles with alternating diploid and haploid phases are a common characteristic of sexually reproducing eukaryotes. Much of our focus in evolutionary biology has been directed toward dynamics in diploid or haploid populations, but we rarely consider selection occurring during both phases when studying evolutionary processes. One of the reasons for this apparent omission is the fact that many flowering plants and metazoans are predominantly diploid with a very short haploid gametic phase. While this gametic phase may be short, it can play a crucial role in fundamental processes including the rate of adaptation, the load of mutation, and the evolution of features such as recombination. In addition, if selection acts in different directions between the two phases, a genetic conflict will occur, impacting the maintenance of genetic variation. Here we provide an overview of theoretical and empirical studies investigating the importance of selection at the haploid gametic phase in predominantly diploid organisms and discuss future directions to improve our understanding of the underlying dynamics and the general implications of haploid selection.

Male-biased gene expression resolves sexual conflict through the evolution of sex-specific genetic architecture

Many genes are subject to contradictory selection pressures in males and females, and balancing selection resulting from sexual conflict has the potential to substantially increase standing genetic diversity in populations and thereby act as an important force in adaptation. However, the underlying causes of sexual conflict, and the potential for resolution, remains hotly debated. Using transcriptome‐resequencing data from male and female guppies, we use a novel approach, combining patterns of genetic diversity and intersexual divergence in allele frequency, to distinguish the different scenarios that give rise to sexual conflict, and how this conflict may be resolved through regulatory evolution. We show that reproductive fitness is the main source of sexual conflict, and this is resolved via the evolution of male‐biased expression. Furthermore, resolution of sexual conflict produces significant differences in genetic architecture between males and females, which in turn lead to specific alleles influencing sex‐specific viability. Together, our findings suggest an important role for sexual conflict in shaping broad patterns of genome diversity, and show that regulatory evolution is a rapid and efficient route to the resolution of conflict.

Inferring sex-specific demographic history from SNP data

The relative female and male contributions to demography are of great importance to better understand the history and dynamics of populations. While earlier studies relied on uniparental markers to investigate sex-specific questions, the increasing amount of sequence data now enables us to take advantage of tens to hundreds of thousands of independent loci from autosomes and the X chromosome. Here, we develop a novel method to estimate effective sex ratios or ESR (defined as the female proportion of the effective population) from allele count data for each branch of a rooted tree topology that summarizes the history of the populations of interest. Our method relies on Kimura’s time-dependent diffusion approximation for genetic drift, and is based on a hierarchical Bayesian model to integrate over the allele frequencies along the branches. We show via simulations that parameters are inferred robustly, even under scenarios that violate some of the model assumptions. Analyzing bovine SNP data, we infer a strongly female-biased ESR in both dairy and beef cattle, as expected from the underlying breeding scheme. Conversely, we observe a strongly male-biased ESR in early domestication times, consistent with an easier taming and management of cows, and/or introgression from wild auroch males, that would both cause a relative increase in male effective population size. In humans, analyzing a subsample of non-African populations, we find a male-biased ESR in Oceanians that may reflect complex marriage patterns in Aboriginal Australians. Because our approach relies on allele count data, it may be applied on a wide range of species.

The history of populations and their social organization is often intricate due to breeding structures, migration patterns or population bottlenecks. Estimation of the female proportion of the effective population (sex ratio) is therefore important to better understand this underlying social structure and dynamics. This question has been mainly investigated so far by comparing genetic variation of mitochondrial DNA and the Y chromosome, two uniparentally inherited markers that reflect the demographic history of females and males, respectively. To overcome the intrinsic limitations of these genetic markers, and to take advantage of the increasing amount of sequence data, we propose a new approach that uses large numbers of independent polymorphisms from autosomes and the X chromosome to estimate sex ratios, throughout the history of populations. This method allows us to confirm a strongly female-biased sex ratio in modern dairy and beef cattle breeds. Yet, we find a strongly male-biased sex ratio during domestication times, consistent with an easier taming and management of cows, and/or introgression from wild auroch males. Analyzing human data from a sample of non-African populations, we find a male bias in Oceanians, possibly indicating complex marriage patterns among Aboriginal Australian groups.

Alternative patterns of sex chromosome differentiation in Aedes aegypti (L)

Background

Some populations of West African Aedes aegypti, the dengue and zika vector, are reproductively incompatible; our earlier study showed that divergence and rearrangements of genes on chromosome 1, which bears the sex locus (M), may be involved. We also previously described a proposed cryptic subspecies SenAae (PK10, Senegal) that had many more high inter-sex F_ST genes on chromosome 1 than did Ae.aegypti aegypti (Aaa, Pai Lom, Thailand). The current work more thoroughly explores the significance of those findings.

Results

Intersex standardized variance (F_ST) of single nucleotide polymorphisms (SNPs) was characterized from genomic exome capture libraries of both sexes in representative natural populations of Aaa and SenAae. Our goal was to identify SNPs that varied in frequency between males and females, and most were expected to occur on chromosome 1. Use of the assembled AaegL4 reference alleviated the previous problem of unmapped genes. Because the M locus gene nix was not captured and not present in AaegL4, the male-determining locus, per se, was not explored. Sex-associated genes were those with F_ST values ≥ 0.100 and/or with increased expected heterozygosity (H_exp, one-sided T-test, p < 0.05) in males. There were 85 genes common to both collections with high inter-sex F_ST values; all genes but one were located on chromosome 1. Aaa showed the expected cluster of high inter-sex F_ST genes proximal to the M locus, whereas SenAae had inter-sex F_ST genes along the length of chromosome 1. In addition, the Aaa M-locus proximal region showed increased H_exp levels in males, whereas SenAae did not. In SenAae, chromosomal rearrangements and subsequent suppressed recombination may have accelerated X-Y differentiation.

Conclusions

The evidence presented here is consistent with differential evolution of proto-Y chromosomes in Aaa and SenAae.

Electronic supplementary material

The online version of this article (doi: 10.1186/s12864-017-4348-4) contains supplementary material, which is available to authorized users.

Genomic Signatures of Sexual Conflict

Sexual conflict is a specific class of intergenomic conflict that describes the reciprocal sex-specific fitness costs generated by antagonistic reproductive interactions. The potential for sexual conflict is an inherent property of having a shared genome between the sexes and, therefore, is an extreme form of an environment-dependent fitness effect. In this way, many of the predictions from environment-dependent selection can be used to formulate expected patterns of genome evolution under sexual conflict. However, the pleiotropic and transmission constraints inherent to having alleles move across sex-specific backgrounds from generation to generation further modulate the anticipated signatures of selection. We outline methods for detecting candidate sexual conflict loci both across and within populations. Additionally, we consider the ability of genome scans to identify sexually antagonistic loci by modeling allele frequency changes within males and females due to a single generation of selection. In particular, we highlight the need to integrate genotype, phenotype, and functional information to truly distinguish sexual conflict from other forms of sexual differentiation.

Sex-Specific Selection and Sex-Biased Gene Expression in Humans and Flies

Sexual dimorphism results from sex-biased gene expression, which evolves when selection acts differently on males and females. While there is an intimate connection between sex-biased gene expression and sex-specific selection, few empirical studies have studied this relationship directly. Here we compare the two on a genome-wide scale in humans and flies. We find a distinctive "Twin Peaks" pattern in humans that relates the strength of sex-specific selection, quantified by genetic divergence between male and female adults at autosomal loci, to the degree of sex-biased expression. Genes with intermediate degrees of sex-biased expression show evidence of ongoing sex-specific selection, while genes with either little or completely sex-biased expression do not. This pattern apparently results from differential viability selection in males and females acting in the current generation. The Twin Peaks pattern is also found in Drosophila using a different measure of sex-specific selection acting on fertility. We develop a simple model that successfully recapitulates the Twin Peaks. Our results suggest that many genes with intermediate sex-biased expression experience ongoing sex-specific selection in humans and flies.