Sex: differences in mutation, recombination, selection, gene flow, and genetic drift

Identifying the sex chromosome and sex determination genes in the cattle tick, Rhipicephalus (Boophilus) microplus

Rhipicephalus (Boophilus) microplus is globally one of the most economically important ectoparasites of cattle costing the agriculture industry billions of dollars annually. Resistance to chemical control measures has prompted the development of novel methods of control. Recent advancements in genetic control measures for human and other animal vectors have utilized sex determination research to manipulate sex ratios, which have shown promising results in mosquitoes namely Aedes aegypti and Anopheles stephensi. Here, we use R. (B.) microplus as a model to provide foundational research to allow similar avenues of investigation in ticks using R. (B.) microplus as a model. Karyotypes for R. (B.) microplus show an XX:XO sex determining system with the largest chromosome being the sex chromosome. Using flow cytometric methods, the size of the sex chromosome was estimated at 526.91 Mb. All measures to identify the sex chromosome within the cattle tick genome assembly associated sex chromosomal characteristics to two chromosomes. This discrepancy between the assembly and karyotypes of the tick led to generating a new genome assembly with a single adult male specimen. The two chromosomes in question aligned with a single scaffold within the new genome that had a length of 513.29 Mb and was the first time the sex chromosome was identified in an Ixodid genome assembly.

Role-reversed polyandry is associated with faster fast-Z in shorebirds

In birds, males are homogametic and carry two copies of the Z chromosome ('ZZ'), while females are heterogametic and exhibit a 'ZW' genotype. The Z chromosome evolves at a faster rate than similarly sized autosomes, a phenomenon termed 'fast-Z evolution'. This is thought to be caused by two independent processes-greater Z chromosome genetic drift owing to a reduced effective population size, and stronger Z chromosome positive selection owing to the exposure of partially recessive alleles to selection. Here, we investigate the relative contributions of these processes by considering the effect of role-reversed polyandry on fast-Z in shorebirds, a paraphyletic group of wading birds that exhibit unusually diverse mating systems. We find stronger fast-Z effects under role-reversed polyandry, which is consistent with particularly strong selection on polyandrous females driving the fixation of recessive beneficial alleles. This result contrasts with previous research in birds, which has tended to implicate a primary role of genetic drift in driving fast-Z variation. We suggest that this discrepancy can be interpreted in two ways-stronger sexual selection acting on polyandrous females overwhelms an otherwise central role of genetic drift, and/or sexual antagonism is also contributing significantly to fast-Z and is exacerbated in sexually dimorphic species.

How much does the unguarded X contribute to sex differences in life span?

Females and males often have markedly different mortality rates and life spans, but it is unclear why these forms of sexual dimorphism evolve. The unguarded X hypothesis contends that dimorphic life spans arise from sex differences in X or Z chromosome copy number (i.e., one copy in the "heterogametic" sex; two copies in the "homogametic" sex), which leads to a disproportionate expression of deleterious mutations by the heterogametic sex (e.g., mammalian males; avian females). Although data on adult sex ratios and sex-specific longevity are consistent with predictions of the unguarded X hypothesis, direct experimental evidence remains scant, and alternative explanations are difficult to rule out. Using a simple population genetic model, we show that the unguarded X effect on sex differential mortality is a function of several reasonably well-studied evolutionary parameters, including the proportion of the genome that is sex linked, the genomic deleterious mutation rate, the mean dominance of deleterious mutations, the relative rates of mutation and strengths of selection in each sex, and the average effect of mutations on survival and longevity relative to their effects on fitness. We review published estimates of these parameters, parameterize our model with them, and show that unguarded X effects are too small to explain observed sex differences in life span across species. For example, sex differences in mean life span are known to often exceed 20% (e.g., in mammals), whereas our parameterized models predict unguarded X effects of a few percent (e.g., 1-3% in Drosophila and mammals). Indeed, these predicted unguarded X effects fall below statistical thresholds of detectability in most experiments, potentially explaining why direct tests of the hypothesis have generated little support for it. Our results suggest that evolution of sexually dimorphic life spans is predominantly attributable to other mechanisms, potentially including "toxic Y" effects and sexual dimorphism for optimal investment in survival versus reproduction.

The rates of introgression and barriers to genetic exchange between hybridizing species: sex chromosomes vs autosomes

Interspecific crossing experiments have shown that sex chromosomes play a major role in reproductive isolation between many pairs of species. However, their ability to act as reproductive barriers, which hamper interspecific genetic exchange, has rarely been evaluated quantitatively compared to Autosomes. This genome-wide limitation of gene flow is essential for understanding the complete separation of species, and thus speciation. Here, we develop a mainland-island model of secondary contact between hybridizing species of an XY (or ZW) sexual system. We obtain theoretical predictions for the frequency of introgressed alleles, and the strength of the barrier to neutral gene flow for the two types of chromosomes carrying multiple interspecific barrier loci. Theoretical predictions are obtained for scenarios where introgressed alleles are rare. We show that the same analytical expressions apply for sex chromosomes and autosomes, but with different sex-averaged effective parameters. The specific features of sex chromosomes (hemizygosity and absence of recombination in the heterogametic sex) lead to reduced levels of introgression on the X (or Z) compared to autosomes. This effect can be enhanced by certain types of sex-biased forces, but it remains overall small (except when alleles causing incompatibilities are recessive). We discuss these predictions in the light of empirical data comprising model-based tests of introgression and cline surveys in various biological systems.

Support for faster and more adaptive Z chromosome evolution in two divergent lepidopteran lineages

The rateof divergence for Z or X chromosomes is usually observed to be greater than autosomes, but the proposed evolutionary causes for this pattern vary, as do empirical results from diverse taxa. Even among moths and butterflies (Lepidoptera), which generally share a single-origin Z chromosome, the handful of available studies give mixed support for faster or more adaptive evolution of the Z chromosome, depending on the species assayed. Here, we examine the molecular evolution of Z chromosomes in two additional lepidopteran species: the Carolina sphinx moth and the monarch butterfly, the latter of which possesses a recent chromosomal fusion yielding a segment of newly Z-linked DNA. We find evidence for both faster and more adaptive Z chromosome evolution in both species, although this effect is strongest in the neo-Z portion of the monarch sex chromosome. The neo-Z is less male-biased than expected of a Z chromosome, and unbiased and female-biased genes drive the signal for adaptive evolution here. Together these results suggest that male-biased gene accumulation and haploid selection have opposing effects on long-term rates of adaptation and may help explain the discrepancies in previous findings as well as the repeated evolution of neo-sex chromosomes in Lepidoptera.

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.

Congruent geographic variation in saccular otolith shape across multiple species of African cichlids

The otoliths of teleost fishes exhibit a great deal of inter- and intra-species shape variation. The ecomorphology of the saccular otolith is often studied by comparing its shape across species and populations inhabiting a range of environments. However, formal tests are often lacking to examine how closely variation in otolith shape follows the genetic drift of a neutral trait. Here, we examine patterns of saccular otolith shape variation in four species of African cichlid fishes, each sampled from three field sites. All four species showed the greatest level of otolith shape variation along two principal component axes, one pertaining to otolith height and another to the prominence of an anterior notch. Fish collected from the same site possessed similarities in saccular otolith shape relative to fish from other sites, and these ‘site-difference’ signatures were consistent across species and observable in both sexes. Sex-differences in saccular otolith shape differed in magnitude from site to site. Population differences in saccular otolith shape did not covary with neutral genetic differentiation between those populations. Otolith height, in particular, displayed large site similarities across species, weak correlation with neutral genetic variation, and strong sex differences, collectively suggesting that otolith shape represents a selectively non-neutral trait.

Life History Effects on Neutral Diversity Levels of Autosomes and Sex Chromosomes

Understanding the determinants of neutral diversity patterns on autosomes and sex chromosomes provides a bedrock for the interpretation of population genetic data; in particular, differences between the two informs our understanding of sex-specific demographic and mutation processes. While sex-specific age-structure and variation in reproductive success have long been known to affect neutral diversity, theoretical descriptions of these effects were complicated and lacking in generality, stymying attempts to relate diversity patterns of species with their life history. Here, we derive general yet simple expressions for these effects. In particular, we show that life history effects on X-to-autosome ratios of pairwise diversity levels (X:A diversity ratios) depend only on the male-to-female ratios of mutation rates, generation times, and reproductive variances. Our results reveal that changing the male-to-female ratio of generation times has opposite effects on X:A ratios of diversity and divergence. They also explain how sex-specific life histories modulate the response of X:A diversity ratios to changes in population size. More generally, they clarify that sex-specific life history-generation times in particular-should have marked effects on X:A diversity ratios in many taxa and enable further investigation of these effects.

The Population Genetics of Evolutionary Rescue in Diploids: X Chromosomal versus Autosomal Rescue

Most population genetic theory assumes that populations adapt to an environmental change without a change in population size. However, environmental changes might be so severe that populations decline in size and, without adaptation, become extinct. This "evolutionary rescue" scenario differs from traditional models of adaptation in that rescue involves a race between adaptation and extinction. While most previous work has usually focused on models of evolutionary rescue in haploids, here we consider diploids. In many species, diploidy introduces a novel feature into adaptation: adaptive evolution might occur either on sex chromosomes or on autosomes. Previous studies of nonrescue adaptation revealed that the relative rates of adaptation on the X chromosome versus autosomes depend on the dominance of beneficial mutations, reflecting differences in effective population size and the efficacy of selection. Here, we extend these results to evolutionary rescue and find that, given equal-sized chromosomes, there is greater parameter space in which the X is more likely to contribute to adaptation than the autosomes relative to standard nonrescue models. We also discuss how subtle effects of dominance can increase the chance of evolutionary rescue in diploids when absolute heterozygote fitness is close to 1. These effects do not arise in standard nonrescue models.

Sex Differences in the Recombination Landscape

Sex differences in overall recombination rates are well known, but little theoretical or empirical attention has been given to how and why sexes differ in their recombination landscapes: the patterns of recombination along chromosomes. In the first scientific review of this phenomenon, we find that recombination is biased toward telomeres in males and more uniformly distributed in females in most vertebrates and many other eukaryotes. Notable exceptions to this pattern exist, however. Fine-scale recombination patterns also frequently differ between males and females. The molecular mechanisms responsible for sex differences remain unclear, but chromatin landscapes play a role. Why these sex differences evolve also is unclear. Hypotheses suggest that they may result from sexually antagonistic selection acting on coding genes and their regulatory elements, meiotic drive in females, selection during the haploid phase of the life cycle, selection against aneuploidy, or mechanistic constraints. No single hypothesis, however, can adequately explain the evolution of sex differences in all cases. Sex-specific recombination landscapes have important consequences for population differentiation and sex chromosome evolution.

Disentangling the intertwined roles of mutation, selection and drift in the mitochondrial genome

Understanding and quantifying the rates of change in the mitochondrial genome is a major component of many areas of biological inquiry, from phylogenetics to human health. A critical parameter in understanding rates of change is estimating the mitochondrial mutation rate (mtDNA MR). Although the first direct estimates of mtDNA MRs were reported almost 20 years ago, the number of estimates has not grown markedly since that time. This is largely owing to the challenges associated with time- and labour-intensive mutation accumulation (MA) experiments. But even MA experiments do not solve a major problem with estimating mtDNA MRs-the challenge of disentangling the role of mutation from other evolutionary forces acting within the cell. Now that it is widely understood that any newly generated mutant allele in the mitochondria will initially be at very low frequency (1/N, where N is the number of mtDNA molecules in the cell), the importance of understanding the effective population size (Ne) of the mtDNA and the size of genetic bottlenecks during gametogenesis and development has come into the spotlight. In addition to these factors regulating the role of genetic drift, advances in our understanding of mitochondrial replication and turnover allow us to more easily envision how natural selection within the cell might favour or purge mutations in multi-copy organellar genomes. Here, we review the unique features of the mitochondrial genome that pose a challenge for accurate MR estimation and discuss ways to overcome those challenges. Estimates of mtDNA MRs remain one of the most widely used parameters in biology, thus accurate quantification and a deeper understanding of how and why they may vary within and between individuals, populations and species is an important goal. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

Local adaptation and the evolution of inversions on sex chromosomes and autosomes

Spatially varying selection with gene flow can favour the evolution of inversions that bind locally adapted alleles together, facilitate local adaptation and ultimately drive genomic divergence between species. Several studies have shown that the rates of spread and establishment of new inversions capturing locally adaptive alleles depend on a suite of evolutionary factors, including the strength of selection for local adaptation, rates of gene flow and recombination, and the deleterious mutation load carried by inversions. Because the balance of these factors is expected to differ between X (or Z) chromosomes and autosomes, opportunities for inversion evolution are likely to systematically differ between these genomic regions, though such scenarios have not been formally modelled. Here, we consider the evolutionary dynamics of X-linked and autosomal inversions in populations evolving at a balance between migration and local selection. We identify three factors that lead to asymmetric rates of X-linked and autosome inversion establishment: (1) sex-biased migration, (2) dominance of locally adapted alleles and (3) chromosome-specific deleterious mutation loads. This theory predicts an elevated rate of fixation, and depressed opportunities for polymorphism, for X-linked inversions. Our survey of data on the genomic distribution of polymorphic and fixed inversions supports both theoretical predictions.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

Linking local adaptation with the evolution of sex differences

Many conspicuous forms of evolutionary diversity occur within species. Two prominent examples include evolutionary divergence between populations differentially adapted to their local environments (local adaptation), and divergence between females and males in response to sex differences in selection (sexual dimorphism sensu lato). These two forms of diversity have inspired vibrant research programmes, yet these fields have largely developed in isolation from one another. Nevertheless, conceptual parallels between these research traditions are striking. Opportunities for local adaptation strike a balance between local selection, which promotes divergence, and gene flow-via dispersal and interbreeding between populations-which constrains it. Sex differences are similarly constrained by fundamental features of inheritance that mimic gene flow. Offspring of each sex inherit genes from same-sex and opposite-sex parents, leading to gene flow between each differentially selected half of the population, and raising the question of how sex differences arise and are maintained. This special issue synthesizes and extends emerging research at the interface between the research traditions of local adaptation and sex differences. Each field can promote understanding of the other, and interactions between local adaptation and sex differences can generate new empirical predictions about the evolutionary consequences of selection that varies across space, time, and between the sexes.This article is part of the theme issue 'Linking local adaptation with the evolution of sex differences'.

Evidence for Faster X Chromosome Evolution in Spiders

In species with chromosomal sex determination, X chromosomes are predicted to evolve faster than autosomes because of positive selection on recessive alleles or weak purifying selection. We investigated X chromosome evolution in Stegodyphus spiders that differ in mating system, sex ratio, and population dynamics. We assigned scaffolds to X chromosomes and autosomes using a novel method based on flow cytometry of sperm cells and reduced representation sequencing. We estimated coding substitution patterns (dN/dS) in a subsocial outcrossing species (S. africanus) and its social inbreeding and female-biased sister species (S. mimosarum), and found evidence for faster-X evolution in both species. X chromosome-to-autosome diversity (piX/piA) ratios were estimated in multiple populations. The average piX/piA estimates of S. africanus (0.57 [95% CI: 0.55-0.60]) was lower than the neutral expectation of 0.75, consistent with more hitchhiking events on X-linked loci and/or a lower X chromosome mutation rate, and we provide evidence in support of both. The social species S. mimosarum has a significantly higher piX/piA ratio (0.72 [95% CI: 0.65-0.79]) in agreement with its female-biased sex ratio. Stegodyphus mimosarum also have different piX/piA estimates among populations, which we interpret as evidence for recurrent founder events. Simulations show that recurrent founder events are expected to decrease the piX/piA estimates in S. mimosarum, thus underestimating the true effect of female-biased sex ratios. Finally, we found lower synonymous divergence on X chromosomes in both species, and the male-to-female substitution ratio to be higher than 1, indicating a higher mutation rate in males.

Complex interactions between local adaptation, phenotypic plasticity and sex affect vulnerability to warming in a widespread marine copepod

Predicting the response of populations to climate change requires an understanding of how various factors affect thermal performance. Genetic differentiation is well known to affect thermal performance, but the effects of sex and developmental phenotypic plasticity often go uncharacterized. We used common garden experiments to test for effects of local adaptation, developmental phenotypic plasticity and individual sex on thermal performance of the ubiquitous copepod, Acartia tonsa (Calanoida, Crustacea) from two populations strongly differing in thermal regimes (Florida and Connecticut, USA). Females had higher thermal tolerance than males in both populations, while the Florida population had higher thermal tolerance compared with the Connecticut population. An effect of developmental phenotypic plasticity on thermal tolerance was observed only in the Connecticut population. Our results show clearly that thermal performance is affected by complex interactions of the three tested variables. Ignoring sex-specific differences in thermal performance may result in a severe underestimation of population-level impacts of warming because of population decline due to sperm limitation. Furthermore, despite having a higher thermal tolerance, low-latitude populations may be more vulnerable to warming as they lack the ability to respond to increases in temperature through phenotypic plasticity.

Where to Cross Over? Defining Crossover Sites in Plants

It is believed that recombination in meiosis serves to reshuffle genetic material from both parents to increase genetic variation in the progeny. At the same time, the number of crossovers is usually kept at a very low level. As a consequence, many organisms need to make the best possible use from the one or two crossovers that occur per chromosome in meiosis. From this perspective, the decision of where to allocate rare crossover events becomes an important issue, especially in self-pollinating plant species, which experience limited variation due to inbreeding. However, the freedom in crossover allocation is significantly limited by other, genetic and non-genetic factors, including chromatin structure. Here we summarize recent progress in our understanding of those processes with a special emphasis on plant genomes. First, we focus on factors which influence the distribution of recombination initiation sites and discuss their effects at both, the single hotspot level and at the chromosome scale. We also briefly explain the aspects of hotspot evolution and their regulation. Next, we analyze how recombination initiation sites translate into the development of crossovers and their location. Moreover, we provide an overview of the sequence polymorphism impact on crossover formation and chromosomal distribution.

Contrasting Patterns of Genomic Diversity Reveal Accelerated Genetic Drift but Reduced Directional Selection on X-Chromosome in Wild and Domestic Sheep Species

Analyses of genomic diversity along the X chromosome and of its correlation with autosomal diversity can facilitate understanding of evolutionary forces in shaping sex-linked genomic architecture. Strong selective sweeps and accelerated genetic drift on the X-chromosome have been inferred in primates and other model species, but no such insight has yet been gained in domestic animals compared with their wild relatives. Here, we analyzed X-chromosome variability in a large ovine data set, including a BeadChip array for 943 ewes from the world’s sheep populations and 110 whole genomes of wild and domestic sheep. Analyzing whole-genome sequences, we observed a substantially reduced X-to-autosome diversity ratio (∼0.6) compared with the value expected under a neutral model (0.75). In particular, one large X-linked segment (43.05–79.25 Mb) was found to show extremely low diversity, most likely due to a high density of coding genes, featuring highly conserved regions. In general, we observed higher nucleotide diversity on the autosomes, but a flat diversity gradient in X-linked segments, as a function of increasing distance from the nearest genes, leading to a decreased X: autosome (X/A) diversity ratio and contrasting to the positive correlation detected in primates and other model animals. Our evidence suggests that accelerated genetic drift but reduced directional selection on X chromosome, as well as sex-biased demographic events, explain low X-chromosome diversity in sheep species. The distinct patterns of X-linked and X/A diversity we observed between Middle Eastern and non-Middle Eastern sheep populations can be explained by multiple migrations, selection, and admixture during the domestic sheep’s recent postdomestication demographic expansion, coupled with natural selection for adaptation to new environments. In addition, we identify important novel genes involved in abnormal behavioral phenotypes, metabolism, and immunity, under selection on the sheep X-chromosome.

Quantifying the relative contributions of the X chromosome, autosomes, and mitochondrial genome to local adaptation

During local adaptation with gene flow, some regions of the genome are inherently more responsive to selection than others. Recent theory predicts that X-linked genes should disproportionately contribute to local adaptation relative to other genomic regions, yet this prediction remains to be tested. We carried out a multigeneration crossing scheme, using two cline-end populations of Drosophila melanogaster, to estimate the relative contributions of the X chromosome, autosomes, and mitochondrial genome to divergence in four traits involved in local adaptation (wing size, resistance to heat, desiccation, and starvation stresses). We found that the mitochondrial genome and autosomes contributed significantly to clinal divergence in three of the four traits. In contrast, the X made no significant contribution to divergence in these traits. Given the small size of the mitochondrial genome, our results indicate that it plays a surprisingly large role in clinal adaptation. In contrast, the X, which represents roughly 20% of the Drosophila genome, contributes negligibly-a pattern that conflicts with theoretical predictions. These patterns reinforce recent work implying a central role of mitochondria in climatic adaptation, and suggest that different genomic regions may play fundamentally different roles in processes of divergence with gene flow.

Evaluating genomic signatures of "the large X-effect" during complex speciation

The ubiquity of the "two rules of speciation"-Haldane's rule and the large X-effect-implies a general, special role for sex chromosomes in the evolution of intrinsic postzygotic reproductive isolation. The recent proliferation of genome-scale analyses has revealed two further general observations: (a) complex speciation involving some form of gene flow is not uncommon, and (b) sex chromosomes in male- and in female-heterogametic taxa tend to show elevated differentiation relative to autosomes. Together, these observations are consistent with speciation histories in which population genetic differentiation at autosomal loci is reduced by gene flow while natural selection against hybrid incompatibilities renders sex chromosomes relatively refractory to gene flow. Here, I summarize multilocus population genetic and population genomic evidence for greater differentiation on the X (or Z) vs. the autosomes and consider the possible causes. I review common population genetic circumstances involving no selection and/or no interspecific gene flow that are nevertheless expected to elevate differentiation on sex chromosomes relative to autosomes. I then review theory for why large X-effects exist for hybrid incompatibilities and, more generally, for loci mediating local adaptation. The observed levels of sex chromosome vs. autosomal differentiation, in many cases, appear consistent with simple explanations requiring neither large X-effects nor gene flow. Discerning signatures of large X-effects during complex speciation will therefore require analyses that go beyond chromosome-scale summaries of population genetic differentiation, explicitly test for differential introgression, and/or integrate experimental genetic data.

Variation in recombination frequency and distribution across eukaryotes: patterns and processes

Recombination, the exchange of DNA between maternal and paternal chromosomes during meiosis, is an essential feature of sexual reproduction in nearly all multicellular organisms. While the role of recombination in the evolution of sex has received theoretical and empirical attention, less is known about how recombination rate itself evolves and what influence this has on evolutionary processes within sexually reproducing organisms. Here, we explore the patterns of, and processes governing recombination in eukaryotes. We summarize patterns of variation, integrating current knowledge with an analysis of linkage map data in 353 organisms. We then discuss proximate and ultimate processes governing recombination rate variation and consider how these influence evolutionary processes. Genome-wide recombination rates (cM/Mb) can vary more than tenfold across eukaryotes, and there is large variation in the distribution of recombination events across closely related taxa, populations and individuals. We discuss how variation in rate and distribution relates to genome architecture, genetic and epigenetic mechanisms, sex, environmental perturbations and variable selective pressures. There has been great progress in determining the molecular mechanisms governing recombination, and with the continued development of new modelling and empirical approaches, there is now also great opportunity to further our understanding of how and why recombination rate varies.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

IBD Sharing between Africans, Neandertals, and Denisovans

Interbreeding between ancestors of humans and other hominins outside of Africa has been studied intensively, while their common history within Africa still lacks proper attention. However, shedding light on human evolution in this time period about which little is known, is essential for understanding subsequent events outside of Africa. We investigate the genetic relationships of humans, Neandertals, and Denisovans by identifying very short DNA segments in the 1000 Genomes Phase 3 data that these hominins share identical by descent (IBD). By focusing on low frequency and rare variants, we identify very short IBD segments with high confidence. These segments reveal events from a very distant past because shorter IBD segments are presumably older than longer ones. We extracted two types of very old IBD segments that are not only shared among humans, but also with Neandertals and/or Denisovans. The first type contains longer segments that are found primarily in Asians and Europeans where more segments are found in South Asians than in East Asians for both Neandertal and Denisovan. These longer segments indicate complex admixture events outside of Africa. The second type consists of shorter segments that are shared mainly by Africans and therefore may indicate events involving ancestors of humans and other ancient hominins within Africa. Our results from the autosomes are further supported by an analysis of chromosome X, on which segments that are shared by Africans and match the Neandertal and/or Denisovan genome were even more prominent. Our results indicate that interbreeding with other hominins was a common feature of human evolution starting already long before ancestors of modern humans left Africa.

Evolutionary mysteries in meiosis

Meiosis is a key event of sexual life cycles in eukaryotes. Its mechanistic details have been uncovered in several model organisms, and most of its essential features have received various and often contradictory evolutionary interpretations. In this perspective, we present an overview of these often 'weird' features. We discuss the origin of meiosis (origin of ploidy reduction and recombination, two-step meiosis), its secondary modifications (in polyploids or asexuals, inverted meiosis), its importance in punctuating life cycles (meiotic arrests, epigenetic resetting, meiotic asymmetry, meiotic fairness) and features associated with recombination (disjunction constraints, heterochiasmy, crossover interference and hotspots). We present the various evolutionary scenarios and selective pressures that have been proposed to account for these features, and we highlight that their evolutionary significance often remains largely mysterious. Resolving these mysteries will likely provide decisive steps towards understanding why sex and recombination are found in the majority of eukaryotes.This article is part of the themed issue 'Weird sex: the underappreciated diversity of sexual reproduction'.

The impact of mating systems and dispersal on fine-scale genetic structure at maternally, paternally and biparentally inherited markers

For decades, studies have focused on how dispersal and mating systems influence genetic structure across populations or social groups. However, we still lack a thorough understanding of how these processes and their interaction shape spatial genetic patterns over a finer scale (tens-hundreds of metres). Using uniparentally inherited markers may help answer these questions, yet their potential has not been fully explored. Here, we use individual-level simulations to investigate the effects of dispersal and mating system on fine-scale genetic structure at autosomal, mitochondrial and Y chromosome markers. Using genetic spatial autocorrelation analysis, we found that dispersal was the major driver of fine-scale genetic structure across maternally, paternally and biparentally inherited markers. However, when dispersal was restricted (mean distance = 100 m), variation in mating behaviour created strong differences in the comparative level of structure detected at maternally and paternally inherited markers. Promiscuity reduced spatial genetic structure at Y chromosome loci (relative to monogamy), whereas structure increased under polygyny. In contrast, mitochondrial and autosomal markers were robust to differences in the specific mating system, although genetic structure increased across all markers when reproductive success was skewed towards fewer individuals. Comparing males and females at Y chromosome vs. mitochondrial markers, respectively, revealed that some mating systems can generate similar patterns to those expected under sex-biased dispersal. This demonstrates the need for caution when inferring ecological and behavioural processes from genetic results. Comparing patterns between the sexes, across a range of marker types, may help us tease apart the processes shaping fine-scale genetic structure.

Regulatory and evolutionary signatures of sex-biased genes on both the X chromosome and the autosomes

Background

Sex is an important but understudied factor in the genetics of human diseases. Analyses using a combination of gene expression data, ENCODE data, and evolutionary data of sex-biased gene expression in human tissues can give insight into the regulatory and evolutionary forces acting on sex-biased genes.

Methods

In this study, we analyzed the differentially expressed genes between males and females. On the X chromosome, we used a novel method and investigated the status of genes that escape X-chromosome inactivation (escape genes), taking into account the clonality of lymphoblastoid cell lines (LCLs). To investigate the regulation of sex-biased differentially expressed genes (sDEG), we conducted pathway and transcription factor enrichment analyses on the sDEGs, as well as analyses on the genomic distribution of sDEGs. Evolutionary analyses were also conducted on both sDEGs and escape genes.

Results

Genome-wide, we characterized differential gene expression between sexes in 462 RNA-seq samples and identified 587 sex-biased genes, or 3.2% of the genes surveyed. On the X chromosome, sDEGs were distributed in evolutionary strata in a similar pattern as escape genes. We found a trend of negative correlation between the gene expression breadth and nonsynonymous over synonymous mutation (dN/dS) ratios, showing a possible pleiotropic constraint on evolution of genes. Genome-wide, nine transcription factors were found enriched in binding to the regions surrounding the transcription start sites of female-biased genes. Many pathways and protein domains were enriched in sex-biased genes, some of which hint at sex-biased physiological processes.

Conclusions

These findings lend insight into the regulatory and evolutionary forces shaping sex-biased gene expression and their involvement in the physiological and pathological processes in human health and diseases.

Electronic supplementary material

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

Genomics of parallel adaptation at two timescales in Drosophila

Two interesting unanswered questions are the extent to which both the broad patterns and genetic details of adaptive divergence are repeatable across species, and the timescales over which parallel adaptation may be observed. Drosophila melanogaster is a key model system for population and evolutionary genomics. Findings from genetics and genomics suggest that recent adaptation to latitudinal environmental variation (on the timescale of hundreds or thousands of years) associated with Out-of-Africa colonization plays an important role in maintaining biological variation in the species. Additionally, studies of interspecific differences between D. melanogaster and its sister species D. simulans have revealed that a substantial proportion of proteins and amino acid residues exhibit adaptive divergence on a roughly few million years long timescale. Here we use population genomic approaches to attack the problem of parallelism between D. melanogaster and a highly diverged conger, D. hydei, on two timescales. D. hydei, a member of the repleta group of Drosophila, is similar to D. melanogaster, in that it too appears to be a recently cosmopolitan species and recent colonizer of high latitude environments. We observed parallelism both for genes exhibiting latitudinal allele frequency differentiation within species and for genes exhibiting recurrent adaptive protein divergence between species. Greater parallelism was observed for long-term adaptive protein evolution and this parallelism includes not only the specific genes/proteins that exhibit adaptive evolution, but extends even to the magnitudes of the selective effects on interspecific protein differences. Thus, despite the roughly 50 million years of time separating D. melanogaster and D. hydei, and despite their considerably divergent biology, they exhibit substantial parallelism, suggesting the existence of a fundamental predictability of adaptive evolution in the genus.

Both local adaptation on short timescales and the long-term accumulation of adaptive differences between species have recently been investigated using comparative genomic and population genomic approaches in several species. However, the repeatability of adaptive evolution at the genetic level is poorly understood. Here we attack this problem by comparing patterns of long and short-term adaptation in Drosophila melanogaster to patterns of adaptation on two timescales in a highly diverged congener, Drosophila hydei. We found, despite the fact that these species diverged from a common ancestor roughly 50 million years ago, the population genomics of latitudinal allele frequency differentiation shows that there is a substantial shared set of genes likely playing a role in the short term adaptive divergence of populations in both species. Analyses of longer-term adaptive protein divergence for the D. hydei-D. mojavensis and D. melanogaster-D. simulans clades reveal a striking level of parallel adaptation. This parallelism includes not only the specific genes/proteins that exhibit adaptive evolution, but extends even to the magnitudes of the selective effects on interspecific protein differences.

The Relative Contributions of the X Chromosome and Autosomes to Local Adaptation

Models of sex chromosome and autosome evolution yield key predictions about the genomic basis of adaptive divergence, and such models have been important in guiding empirical research in comparative genomics and studies of speciation. In addition to the adaptive differentiation that occurs between species over time, selection also favors genetic divergence across geographic space, with subpopulations of single species evolving conspicuous differences in traits involved in adaptation to local environmental conditions. The potential contribution of sex chromosomes (the X or Z) to local adaptation remains unclear, as we currently lack theory that directly links spatial variation in selection to local adaptation of X-linked and autosomal genes. Here, we develop population genetic models that explicitly consider the effects of genetic dominance, effective population size, and sex-specific migration and selection on the relative contributions of X-linked and autosomal genes to local adaptation. We show that X-linked genes should nearly always disproportionately contribute to local adaptation in the presence of gene flow. We also show that considerations of dominance and effective population size-which play pivotal roles in the theory of faster-X adaptation between species-have surprisingly little influence on the relative contribution of the X chromosome to local adaptation. Instead, sex-biased migration is the primary mediator of the strength of spatial large-X effects. Our results yield novel predictions about the role of sex chromosomes in local adaptation. We outline empirical approaches in evolutionary quantitative genetics and genomics that could build upon this new theory.

Genomic signatures of sex-biased demography: progress and prospects

Sex-biased demographic events have played a crucial role in shaping human history. Many of these processes affect genetic variation and can therefore leave detectable signatures in the genome because autosomal, X-linked, Y-linked, and mitochondrial DNA inheritance differ between sexes. Here, we discuss how sex-biased processes shape patterns of genetic diversity across the genome, review recent genomic evidence for sex-biased demography in modern human populations, and suggest directions for future research.

Accumulation of Deleterious Mutations Near Sexually Antagonistic Genes

Mutation generates a steady supply of genetic variation that, while occasionally useful for adaptation, is more often deleterious for fitness. Recent research has emphasized that the fitness effects of mutations often differ between the sexes, leading to important evolutionary consequences for the maintenance of genetic variation and long-term population viability. Some forms of sex-specific selection-i.e., stronger purifying selection in males than females-can help purge a population's load of female-harming mutations and promote population growth. Other scenarios-e.g., sexually antagonistic selection, in which mutations that harm females are beneficial for males-inflate genetic loads and potentially dampen population viability. Evolutionary processes of sexual antagonism and purifying selection are likely to impact the evolutionary dynamics of different loci within a genome, yet theory has mostly ignored the potential for interactions between such loci to jointly shape the evolutionary genetic basis of female and male fitness variation. Here, we show that sexually antagonistic selection at a locus tends to elevate the frequencies of deleterious alleles at tightly linked loci that evolve under purifying selection. Moreover, haplotypes that segregate for different sexually antagonistic alleles accumulate different types of deleterious mutations. Haplotypes that carry female-benefit sexually antagonistic alleles preferentially accumulate mutations that are primarily male harming, whereas male-benefit haplotypes accumulate mutations that are primarily female harming. The theory predicts that sexually antagonistic selection should shape the genomic organization of genetic variation that differentially impacts female and male fitness, and contribute to sexual dimorphism in the genetic basis of fitness variation.

Towards cancer-aware life-history modelling

Studies of body size evolution, and life-history theory in general, are conducted without taking into account cancer as a factor that can end an organism's reproductive lifespan. This reflects a tacit assumption that predation, parasitism and starvation are of overriding importance in the wild. We argue here that even if deaths directly attributable to cancer are a rarity in studies of natural populations, it remains incorrect to infer that cancer has not been of importance in shaping observed life histories. We present first steps towards a cancer-aware life-history theory, by quantifying the decrease in the length of the expected reproductively active lifespan that follows from an attempt to grow larger than conspecific competitors. If all else is equal, a larger organism is more likely to develop cancer, but, importantly, many factors are unlikely to be equal. Variations in extrinsic mortality as well as in the pace of life—larger organisms are often near the slow end of the fast–slow life-history continuum—can make realized cancer incidences more equal across species than what would be observed in the absence of adaptive responses to cancer risk (alleviating the so-called Peto's paradox). We also discuss reasons why patterns across species can differ from within-species predictions. Even if natural selection diminishes cancer susceptibility differences between species, within-species differences can remain. In many sexually dimorphic cases, we predict males to be more cancer-prone than females, forming an understudied component of sexual conflict.

Multicopy gene family evolution on primate Y chromosomes

Background

The primate Y chromosome is distinguished by a lack of inter-chromosomal recombination along most of its length, extensive gene loss, and a prevalence of repetitive elements. A group of genes on the male-specific portion of the Y chromosome known as the “ampliconic genes” are present in multiple copies that are sometimes part of palindromes, and that undergo a form of intra-chromosomal recombination called gene conversion, wherein the nucleotides of one copy are homogenized by those of another. With the aim of further understanding gene family evolution of these genes, we collected nucleotide sequence and gene copy number information for several species of papionin monkey. We then tested for evidence of gene conversion, and developed a novel statistical framework to evaluate alternative models of gene family evolution using our data combined with other information from a human, a chimpanzee, and a rhesus macaque.

Results

Our results (i) recovered evidence for several novel examples of gene conversion in papionin monkeys and indicate that (ii) ampliconic gene families evolve faster than autosomal gene families and than single-copy genes on the Y chromosome and that (iii) Y-linked singleton and autosomal gene families evolved faster in humans and chimps than they do in the other Old World Monkey lineages we studied.

Conclusions

Rapid evolution of ampliconic genes cannot be attributed solely to residence on the Y chromosome, nor to variation between primate lineages in the rate of gene family evolution. Instead other factors, such as natural selection and gene conversion, appear to play a role in driving temporal and genomic evolutionary heterogeneity in primate gene families.

Electronic supplementary material

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

Construction of Ultradense Linkage Maps with Lep-MAP2: Stickleback F2 Recombinant Crosses as an Example

High-density linkage maps are important tools for genome biology and evolutionary genetics by quantifying the extent of recombination, linkage disequilibrium, and chromosomal rearrangements across chromosomes, sexes, and populations. They provide one of the best ways to validate and refine de novo genome assemblies, with the power to identify errors in assemblies increasing with marker density. However, assembly of high-density linkage maps is still challenging due to software limitations. We describe Lep-MAP2, a software for ultradense genome-wide linkage map construction. Lep-MAP2 can handle various family structures and can account for achiasmatic meiosis to gain linkage map accuracy. Simulations show that Lep-MAP2 outperforms other available mapping software both in computational efficiency and accuracy. When applied to two large F₂-generation recombinant crosses between two nine-spined stickleback (Pungitius pungitius) populations, it produced two high-density (∼6 markers/cM) linkage maps containing 18,691 and 20,054 single nucleotide polymorphisms. The two maps showed a high degree of synteny, but female maps were 1.5–2 times longer than male maps in all linkage groups, suggesting genome-wide recombination suppression in males. Comparison with the genome sequence of the three-spined stickleback (Gasterosteus aculeatus) revealed a high degree of interspecific synteny with a low frequency (<5%) of interchromosomal rearrangements. However, a fairly large (ca. 10 Mb) translocation from autosome to sex chromosome was detected in both maps. These results illustrate the utility and novel features of Lep-MAP2 in assembling high-density linkage maps, and their usefulness in revealing evolutionarily interesting properties of genomes, such as strong genome-wide sex bias in recombination rates.

Z chromosome divergence, polymorphism and relative effective population size in a genus of lekking birds

Sex chromosomes contribute disproportionately to species boundaries as they diverge faster than autosomes and often have reduced diversity. Their hemizygous nature contributes to faster divergence and reduced diversity, as do some types of selection. In birds, other factors (mating system and bottlenecks) can further decrease the effective population size of Z-linked loci and accelerate divergence (Fast-Z). We assessed Z-linked divergence and effective population sizes for two polygynous sage-grouse species and compared them to estimates from birds with various mating systems. We found lower diversity and higher FST for Z-linked loci than for autosomes, as expected. The π(Z)/π(A) ratio was 0.38 in Centrocercus minimus, 0.48 in Centrocercus urophasianus and 0.59 in a diverged, parapatric population of C. urophasianus, a broad range given the mating system among these groups is presumably equivalent. The full data set had unequal males and females across groups, so we compared an equally balanced reduced set of C. minimus and individuals pooled from both C. urophasianus subgroups recovering similar estimates: 0.54 for C. urophasianus and 0.38 for C. minimus. We provide further evidence that N(eZ)/N(eA) in birds is often lower than expected under random mating or monogamy. The lower ratio in C. minimus could be a consequence of stronger selection or drift acting on Z loci during speciation, as this species differs strongly from C. urophasianus in sexually selected characters with minimal mitochondrial divergence. As C. minimus also exhibited lower genomic diversity, it is possible that a more severe demographic history may contribute to its lower ratio.

The effect of temperature on the male and female recombination landscape of barley

Barley (Hordeum vulgare) is a crop of global significance. However, a third of the genes of barley are largely inaccessible to conventional breeding programmes as crossovers are localised to the ends of the chromosomes. This work examines whether crossovers can be shifted to more proximal regions simply by elevating growth temperature. We utilised a genome-wide marker set for linkage analysis combined with cytological mapping of crossover events to examine the recombination landscape of plants grown at different temperatures. We found that barley shows heterochiasmy, that is, differences between female and male recombination frequencies. In addition, we found that elevated temperature significantly changes patterns of recombination in male meiosis only, with a repositioning of Class I crossovers determined by cytological mapping of HvMLH3 foci. We show that the length of synaptonemal complexes in male meiocytes increases in response to temperature. The results demonstrate that the distribution of crossover events are malleable and can be shifted to proximal regions by altering the growth temperature. The shift in recombination is the result of altering the distribution of Class I crossovers, but the higher recombination at elevated temperatures is potentially not the result of an increase in Class I events.

Beyond 2/3 and 1/3: The Complex Signatures of Sex-Biased Admixture on the X Chromosome

Sex-biased demography, in which parameters governing migration and population size differ between females and males, has been studied through comparisons of X chromosomes, which are inherited sex-specifically, and autosomes, which are not. A common form of sex bias in humans is sex-biased admixture, in which at least one of the source populations differs in its proportions of females and males contributing to an admixed population. Studies of sex-biased admixture often examine the mean ancestry for markers on the X chromosome in relation to the autosomes. A simple framework noting that in a population with equally many females and males, two-thirds of X chromosomes appear in females, suggests that the mean X-chromosomal admixture fraction is a linear combination of female and male admixture parameters, with coefficients 2/3 and 1/3, respectively. Extending a mechanistic admixture model to accommodate the X chromosome, we demonstrate that this prediction is not generally true in admixture models, although it holds in the limit for an admixture process occurring as a single event. For a model with constant ongoing admixture, we determine the mean X-chromosomal admixture, comparing admixture on female and male X chromosomes to corresponding autosomal values. Surprisingly, in reanalyzing African-American genetic data to estimate sex-specific contributions from African and European sources, we find that the range of contributions compatible with the excess African ancestry on the X chromosome compared to autosomes has a wide spread, permitting scenarios either without male-biased contributions from Europe or without female-biased contributions from Africa.

Evolution: The Island of Misfit Mammoths

The genomes of two woolly mammoths have been sequenced. One of the last survivors had reduced genetic diversity. Although divergent in their mitochondrial genomes, the mammoths had similar nuclear genomes, a finding germane to elephant conservation.

Does sex-biased dispersal account for the lack of geographic and host-associated differentiation in introduced populations of an aphid parasitoid?

Host recognition and use in female parasitoids strongly relies on host fidelity, a plastic behavior which can significantly restrict the host preferences of parasitoids, thus reducing the gene flow between parasitoid populations attacking different insect hosts. However, the effect of migrant males on the genetic differentiation of populations has been frequently ignored in parasitoids, despite its known impact on gene flow between populations. Hence, we studied the extent of gene flow mediated by female and male parasitoids by assessing sibship relationships among parasitoids within and between populations, and its impact on geographic and host-associated differentiation in the aphid parasitoid Aphidius ervi. We report evidences of a high gene flow among parasitoid populations on different aphid hosts and geographic locations. The high gene flow among parasitoid populations was found to be largely male mediated, suggested by significant differences in the distribution of full-sib and paternal half-sib dyads of parasitoid populations.

Reduced representation genome sequencing suggests low diversity on the sex chromosomes of tonkean macaque monkeys

In species with separate sexes, social systems can differ in the relative variances of male versus female reproductive success. Papionin monkeys (macaques, mangabeys, mandrills, drills, baboons, and geladas) exhibit hallmarks of a high variance in male reproductive success, including a female-biased adult sex ratio and prominent sexual dimorphism. To explore the potential genomic consequences of such sex differences, we used a reduced representation genome sequencing approach to quantifying polymorphism at sites on autosomes and sex chromosomes of the tonkean macaque (Macaca tonkeana), a species endemic to the Indonesian island of Sulawesi. The ratio of nucleotide diversity of the X chromosome to that of the autosomes was less than the value (0.75) expected with a 1:1 sex ratio and no sex differences in the variance in reproductive success. However, the significance of this difference was dependent on which outgroup was used to standardize diversity levels. Using a new model that includes the effects of varying population size, sex differences in mutation rate between the autosomes and X chromosome, and GC-biased gene conversion (gBGC) or selection on GC content, we found that the maximum-likelihood estimate of the ratio of effective population size of the X chromosome to that of the autosomes was 0.68, which did not differ significantly from 0.75. We also found evidence for 1) a higher level of purifying selection on genic than nongenic regions, 2) gBGC or natural selection favoring increased GC content, 3) a dynamic demography characterized by population growth and contraction, 4) a higher mutation rate in males than females, and 5) a very low polymorphism level on the Y chromosome. These findings shed light on the population genomic consequences of sex differences in the variance in reproductive success, which appear to be modest in the tonkean macaque; they also suggest the occurrence of hitchhiking on the Y chromosome.

Phylogeography of the Lutzomyia gomezi (Diptera: Phlebotominae) on the Panama Isthmus

BACKGROUND

Lutzomyia gomezi (Nitzulescu, 1931) is one of the main Leishmania (Vianna) panamensis vectors in Panama, and despite its medical significance, there are no population genetic studies regarding this species. In this study, we used the sequences of the mitochondrial gene cytochrome b/start of NADH1 and the nuclear elongation gene α-1 in order to analyze genetic variation and phylogeographic structure of the Lu. gomezi populations.

METHODS

A total of 86 Lu. gomezi individuals were captured in 38 locations where cutaneous leishmaniasis occurred. DNA was extracted with phenol/chloroform methods and amplification of genes was performed using PCR primers for mitochondrial and nuclear markers.

RESULTS

We found a total of 37 and 26 haplotypes of mitochondrial and nuclear genes, high haplotype diversity (h) for all three populations were detected with both molecular markers. Nucleotide diversity (π) was estimated to be high for all three populations with the mitochondrial marker, which was opposite to the estimate with the nuclear marker. In the AMOVA Φst recorded moderate (mitochondrial) and small (nuclear) population structure with statistical significance among populations. The analysis of the fixation index (Fst) used to measure the differentiation of populations showed that with the exception of the population located in the region of Bocas del Toro, the other populations presented with minor genetic differentiation. The median-Joining network of the mitochondrial marker reveled three clusters and recorded four haplotypes exclusively of localities sampled from Western Panama, demonstrating strong divergence. We found demographic population expansion with Fu´s Fs neutrality test. In the analysis mismatch distribution was observed as a bimodal curve.

CONCLUSION

Lu. gomezi is a species with higher genetic pool or variability and mild population structure, due to possible capacity migration and local adaptation to environmental changes or colonization potential. Thus, knowledge of the genetic population and evolutionary history is useful to understand the implications of different population genetic structures for cutaneous leishmaniasis epidemiology.