On the effective size of populations with separate sexes, with particular reference to sex-linked genes

Mitonuclear Interactions and the Origin of Macaque Societies

In most eukaryotes, aerobic respiration requires interactions between autosomally encoded genes (Ninteract genes) and mitochondrial DNA, RNA, and protein. In species where females are philopatric, contrasting distributions of genetic variation in mitochondrial and nuclear genomes create variation in mitonuclear interactions that may be subject to natural selection. To test this expectation, we turned to a group with extreme female philopatry: the macaque monkeys. We examined four genomic data sets from (1) wild caught and (2) captive populations of rhesus macaque, which is the most widely distributed nonhuman primate, and (3) the stump-tailed macaque and (4) a subspecies of longtail macaque, both of whose mitochondrial DNA is introgressed from a highly diverged ancestor. We identified atypically long runs of homozygosity, low polymorphism, high differentiation, and/or rapid protein evolution associated with Ninteract genes compared with non-Ninteract genes. These metrics suggest a subset of Ninteract genes were independently subject to atypically pervasive natural selection in multiple species. These findings suggest that natural selection on mitonuclear interactions could have influenced several aspects of macaque societies including species diversity, ecological breadth, female-biased adult sex ratio and demography, sexual dimorphism, and mitonuclear phylogenomics.

X chromosomes show relaxed selection and complete somatic dosage compensation across Timema stick insect species

Sex chromosomes have evolved repeatedly across the tree of life. As they are present in different copy numbers in males and females, they are expected to experience different selection pressures than the autosomes, with consequences including a faster rate of evolution, increased accumulation of sexually antagonistic alleles and the evolution of dosage compensation. Whether these consequences are general or linked to idiosyncrasies of specific taxa is not clear as relatively few taxa have been studied thus far. Here, we use whole-genome sequencing to identify and characterize the evolution of the X chromosome in five species of Timema stick insects with XX:X0 sex determination. The X chromosome had a similar size (approximately 12% of the genome) and gene content across all five species, suggesting that the X chromosome originated prior to the diversification of the genus. Genes on the X showed evidence of relaxed selection (elevated dN/dS) and a slower evolutionary rate (dN + dS) than genes on the autosomes, likely due to sex-biased mutation rates. Genes on the X also showed almost complete dosage compensation in somatic tissues (heads and legs), but dosage compensation was absent in the reproductive tracts. Contrary to prediction, sex-biased genes showed little enrichment on the X, suggesting that the advantage X-linkage provides to the accumulation of sexually antagonistic alleles is weak. Overall, we found the consequences of X-linkage on gene sequences and expression to be similar across Timema species, showing the characteristics of the X chromosome are surprisingly consistent over 30 million years of evolution.

Consequences of coupled barriers to gene flow for the build-up of genomic differentiation

Theory predicts that when different barriers to gene flow become coincident, their joint effects enhance reproductive isolation and genomic divergence beyond their individual effects, but empirical tests of this "coupling" hypothesis are rare. Here, we analyze patterns of gene exchange among populations of European corn borer moths that vary in the number of acting barriers, allowing for comparisons of genomic variation when barrier traits or loci are in coincident or independent states. We find that divergence is mainly restricted to barrier loci when populations differ by a single barrier, whereas the coincidence of temporal and behavioral barriers is associated with divergence of two chromosomes harboring barrier loci. Furthermore, differentiation at temporal barrier loci increases in the presence of behavioral divergence and differentiation at behavioral barrier loci increases in the presence of temporal divergence. Our results demonstrate how the joint action of coincident barrier effects leads to levels of genomic differentiation that far exceed those of single barriers acting alone, consistent with theory arguing that coupling allows indirect selection to combine with direct selection and thereby lead to a stronger overall barrier to gene flow. Thus, the state of barriers-independent or coupled-strongly influences the accumulation of genomic differentiation.

A chromosome-level genome assembly of Cairina moschata and comparative genomic analyses

Background

The Muscovy duck (Cairina moschata) is an economically important duck species, with favourable growth and carcass composition parameters in comparison to other ducks. However, limited genomic resources for Muscovy duck hinder our understanding of its evolution and genetic diversity.

Results

We combined linked-reads sequencing technology and reference-guided methods for de novo genome assembly. The final draft assembly was 1.12 Gbp with 29 autosomes, one sex chromosome and 4,583 unlocalized scaffolds with an N50 size of 77.35 Mb. Based on universal single-copy orthologues (BUSCO), the draft genome assembly completeness was estimated to be 93.30 %. Genome annotation identified 15,580 genes, with 15,537 (99.72 %) genes annotated in public databases. We conducted comparative genomic analyses and found that species-specific and rapidly expanding gene families (compared to other birds) in Muscovy duck are mainly involved in Calcium signaling, Adrenergic signaling in cardiomyocytes, and GnRH signaling pathways. In comparison to the common domestic duck (Anas platyrhynchos), we identified 104 genes exhibiting strong signals of adaptive evolution (Ka/Ks > 1). Most of these genes were associated with immune defence pathways (e.g. IFNAR1 and TLR5). This is indicative of the existence of differences in the immune responses between the two species. Additionally, we combined divergence and polymorphism data to demonstrate the “faster-Z effect” of chromosome evolution.

Conclusions

The chromosome-level genome assembly of Muscovy duck and comparative genomic analyses provide valuable resources for future molecular ecology studies, as well as the evolutionary arms race between the host and influenza viruses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07897-4.

Sex-linked genetic diversity and differentiation in a globally distributed avian species complex

Sex chromosomes often bear distinct patterns of genetic variation due to unique patterns of inheritance and demography. The processes of mutation, recombination, genetic drift and selection also influence rates of evolution on sex chromosomes differently than autosomes. Measuring such differences provides information about how these processes shape genomic variation and their roles in the origin of species. To test hypotheses and predictions about patterns of autosomal and sex-linked genomic diversity and differentiation, we measured population genetic statistics within and between populations and subspecies of the barn swallow (Hirundo rustica) and performed explicit comparisons between autosomal and Z-linked genomic regions. We first tested for evidence of low Z-linked genetic diversity and high Z-linked population differentiation relative to autosomes, then for evidence that the Z chromosome bears greater ancestry information due to faster lineage sorting. Finally, we investigated geographical clines across hybrid zones for evidence that the Z chromosome is resistant to introgression due to selection against hybrids. We found evidence that the barn swallow mating system, demographic history and linked selection each contribute to low Z-linked diversity and high Z-linked differentiation. While incomplete lineage sorting is rampant across the genome, our results indicate faster sorting of ancestral polymorphism on the Z. Finally, hybrid zone analyses indicate barriers to introgression on the Z chromosome, suggesting that sex-linked traits are important in reproductive isolation, especially in migratory divide regions. Our study highlights how selection, gene flow and demography shape sex-linked genetic diversity and underlines the relevance of the Z chromosome in speciation.

How Good Are Predictions of the Effects of Selective Sweeps on Levels of Neutral Diversity?

Selective sweeps are thought to play a significant role in shaping patterns of variability across genomes; accurate predictions of their effects are, therefore, important for understanding these patterns. A commonly used model of selective sweeps assumes that alleles sampled at the end of a sweep, and that fail to recombine with wild-type haplotypes during the sweep, coalesce instantaneously, leading to a simple expression for sweep effects on diversity. It is shown here that there can be a significant probability that a pair of alleles sampled at the end of a sweep coalesce during the sweep before a recombination event can occur, reducing their expected coalescent time below that given by the simple approximation. Expressions are derived for the expected reductions in pairwise neutral diversities caused by both single and recurrent sweeps in the presence of such within-sweep coalescence, although the effects of multiple recombination events during a sweep are only treated heuristically. The accuracies of the resulting expressions were checked against the results of simulations. For even moderate ratios of the recombination rate to the selection coefficient, the simple approximation can be substantially inaccurate. The selection model used here can be applied to favorable mutations with arbitrary dominance coefficients, to sex-linked loci with sex-specific selection coefficients, and to inbreeding populations. Using the results from this model, the expected differences between the levels of variability on X chromosomes and autosomes with selection at linked sites are discussed, and compared with data on a population of Drosophila melanogaster.

Female-specific resource limitation does not make the opportunity for selection more female biased

Competition for limiting resources and stress can magnify variance in fitness and therefore selection. But even in a common environment, the strength of selection can differ across the sexes, as their fitness is often limited by different factors. Indeed, most taxa show stronger selection in males, a bias often ascribed to intense competition for access to mating partners. This sex bias could reverberate on many aspects of evolution, from speed of adaptation to genome evolution. It is unclear, however, whether stronger opportunity for selection in males is a pattern robust to sex-specific stress or resource limitation. We test this in the model species Callosobruchus maculatus by comparing female and male opportunity for selection (i) with and without limitation of quality oviposition sites, and (ii) under delayed age at oviposition. Decreasing the abundance of the resource key to females or increasing their reproductive age was challenging, as shown by a reduction in mean fitness, but opportunity for selection remained stronger in males across all treatments, and even more so when oviposition sites were limiting. This suggests that males remain the more variable sex independent of context, and that the opportunity for selection through males is indirectly affected by female-specific resource limitation.

Demographic History of the Human Commensal Drosophila melanogaster

The cohabitation of Drosophila melanogaster with humans is nearly ubiquitous. Though it has been well established that this fly species originated in sub-Saharan Africa, and only recently has spread globally, many details of its swift expansion remain unclear. Elucidating the demographic history of D. melanogaster provides a unique opportunity to investigate how human movement might have impacted patterns of genetic diversity in a commensal species, as well as providing neutral null models for studies aimed at identifying genomic signatures of local adaptation. Here, we use whole-genome data from five populations (Africa, North America, Europe, Central Asia, and the South Pacific) to carry out demographic inferences, with particular attention to the inclusion of migration and admixture. We demonstrate the importance of these parameters for model fitting and show that how previous estimates of divergence times are likely to be significantly underestimated as a result of not including them. Finally, we discuss how human movement along early shipping routes might have shaped the present-day population structure of D. melanogaster.

Maximum avoidance of inbreeding in haplodiploid populations

Maximum avoidance of inbreeding (MAI) is a mating system, in which mates are as distantly related as possible. Although theoretical aspects and applications of MAI in diploid populations have been studied by many researchers, extension of MAI to haplodiploid populations is an unresolved problem. In this paper, this problem is addressed, and the following conclusions are derived. For a haplodiploid population with a Fibonacci number of females, a set of mating systems (one cycle MAI-hd) to avoid inbreeding to the maximum after one cycle practice of the set can be defined. But unlike MAI in diploid populations, repetition of one cycle MAI-hd cannot be MAI in the global range of generations. Numerical comparison with random mating and circular half-sib mating shows that as in diploid populations, repetition of one cycle MAI-hd in haplodiploid populations attains a lower inbreeding coefficient in early generations at the expense of a higher asymptotic rate of inbreeding.

Disentangling the Causes for Faster-X Evolution in Aphids

The faster evolution of X chromosomes has been documented in several species, and results from the increased efficiency of selection on recessive alleles in hemizygous males and/or from increased drift due to the smaller effective population size of X chromosomes. Aphids are excellent models for evaluating the importance of selection in faster-X evolution because their peculiar life cycle and unusual inheritance of sex chromosomes should generally lead to equivalent effective population sizes for X and autosomes. Because we lack a high-density genetic map for the pea aphid, whose complete genome has been sequenced, we first assigned its entire genome to the X or autosomes based on ratios of sequencing depth in males (X0) to females (XX). Then, we computed nonsynonymous to synonymous substitutions ratios (dN/dS) for the pea aphid gene set and found faster evolution of X-linked genes. Our analyses of substitution rates, together with polymorphism and expression data, showed that relaxed selection is likely to be the greatest contributor to faster-X because a large fraction of X-linked genes are expressed at low rates and thus escape selection. Yet, a minor role for positive selection is also suggested by the difference between substitution rates for X and autosomes for male-biased genes (but not for asexual female-biased genes) and by lower Tajima’s D for X-linked compared with autosomal genes with highly male-biased expression patterns. This study highlights the relevance of organisms displaying alternative chromosomal inheritance to the understanding of forces shaping genome evolution.

Genetic Diversity on the Sex Chromosomes

Levels and patterns of genetic diversity can provide insights into a population’s history. In species with sex chromosomes, differences between genomic regions with unique inheritance patterns can be used to distinguish between different sets of possible demographic and selective events. This review introduces the differences in population history for sex chromosomes and autosomes, provides the expectations for genetic diversity across the genome under different evolutionary scenarios, and gives an introductory description for how deviations in these expectations are calculated and can be interpreted. Predominantly, diversity on the sex chromosomes has been used to explore and address three research areas: 1) Mating patterns and sex-biased variance in reproductive success, 2) signatures of selection, and 3) evidence for modes of speciation and introgression. After introducing the theory, this review catalogs recent studies of genetic diversity on the sex chromosomes across species within the major research areas that sex chromosomes are typically applied to, arguing that there are broad similarities not only between male-heterogametic (XX/XY) and female-heterogametic (ZZ/ZW) sex determination systems but also any mating system with reduced recombination in a sex-determining region. Further, general patterns of reduced diversity in nonrecombining regions are shared across plants and animals. There are unique patterns across populations with vastly different patterns of mating and speciation, but these do not tend to cluster by taxa or sex determination system.

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.

Cultural hitchhiking and competition between patrilineal kin groups explain the post-Neolithic Y-chromosome bottleneck

In human populations, changes in genetic variation are driven not only by genetic processes, but can also arise from cultural or social changes. An abrupt population bottleneck specific to human males has been inferred across several Old World (Africa, Europe, Asia) populations 5000–7000 BP. Here, bringing together anthropological theory, recent population genomic studies and mathematical models, we propose a sociocultural hypothesis, involving the formation of patrilineal kin groups and intergroup competition among these groups. Our analysis shows that this sociocultural hypothesis can explain the inference of a population bottleneck. We also show that our hypothesis is consistent with current findings from the archaeogenetics of Old World Eurasia, and is important for conceptions of cultural and social evolution in prehistory.

A population bottleneck 5000-7000 years ago in human males, but not females, has been inferred across several African, European and Asian populations. Here, Zeng and colleagues synthesize theory and data to suggest that competition among patrilineal kin groups produced the bottleneck pattern.

Sex chromosomes and speciation in birds and other ZW systems

Theory and empirical patterns suggest a disproportionate role for sex chromosomes in evolution and speciation. Focusing on ZW sex determination (females ZW, males ZZ; the system in birds, many snakes, and lepidopterans), I review how evolutionary dynamics are expected to differ between the Z, W and the autosomes, discuss how these differences may lead to a greater role of the sex chromosomes in speciation and use data from birds to compare relative evolutionary rates of sex chromosomes and autosomes. Neutral mutations, partially or completely recessive beneficial mutations, and deleterious mutations under many conditions are expected to accumulate faster on the Z than on autosomes. Sexually antagonistic polymorphisms are expected to arise on the Z, raising the possibility of the spread of preference alleles. The faster accumulation of many types of mutations and the potential for complex evolutionary dynamics of sexually antagonistic traits and preferences contribute to a role for the Z chromosome in speciation. A quantitative comparison among a wide variety of bird species shows that the Z tends to have less within-population diversity and greater between-species differentiation than the autosomes, likely due to both adaptive evolution and a greater rate of fixation of deleterious alleles. The W chromosome also shows strong potential to be involved in speciation, in part because of its co-inheritance with the mitochondrial genome. While theory and empirical evidence suggest a disproportionate role for sex chromosomes in speciation, the importance of sex chromosomes is moderated by their small size compared to the whole genome.

Faster-X evolution: Theory and evidence from Drosophila

A faster rate of adaptive evolution of X-linked genes compared with autosomal genes can be caused by the fixation of recessive or partially recessive advantageous mutations, due to the full expression of X-linked mutations in hemizygous males. Other processes, including recombination rate and mutation rate differences between X chromosomes and autosomes, may also cause faster evolution of X-linked genes. We review population genetics theory concerning the expected relative values of variability and rates of evolution of X-linked and autosomal DNA sequences. The theoretical predictions are compared with data from population genomic studies of several species of Drosophila. We conclude that there is evidence for adaptive faster-X evolution of several classes of functionally significant nucleotides. We also find evidence for potential differences in mutation rates between X-linked and autosomal genes, due to differences in mutational bias towards GC to AT mutations. Many aspects of the data are consistent with the male hemizygosity model, although not all possible confounding factors can be excluded.

Migration-Selection Balance Drives Genetic Differentiation in Genes Associated with High-Altitude Function in the Speckled Teal (Anas flavirostris) in the Andes

Local adaptation frequently occurs across populations as a result of migration-selection balance between divergent selective pressures and gene flow associated with life in heterogeneous landscapes. Studying the effects of selection and gene flow on the adaptation process can be achieved in systems that have recently colonized extreme environments. This study utilizes an endemic South American duck species, the speckled teal (Anas flavirostris), which has both high- and low-altitude populations. High-altitude speckled teal (A. f. oxyptera) are locally adapted to the Andean environment and mostly allopatric from low-altitude birds (A. f. flavirostris); however, there is occasional gene flow across altitudinal gradients. In this study, we used next-generation sequencing to explore genetic patterns associated with high-altitude adaptation in speckled teal populations, as well as the extent to which the balance between selection and migration have affected genetic architecture. We identified a set of loci with allele frequencies strongly correlated with altitude, including those involved in the insulin-like signaling pathway, bone morphogenesis, oxidative phosphorylation, responders to hypoxia-induced DNA damage, and feedback loops to the hypoxia-inducible factor pathway. These same outlier loci were found to have depressed gene flow estimates, as well as being highly concentrated on the Z-chromosome. Our results suggest a multifactorial response to life at high altitudes through an array of interconnected pathways that are likely under positive selection and whose genetic components seem to be providing an effective genomic barrier to interbreeding, potentially functioning as an avenue for population divergence and speciation.

Hill-Robertson Interference Reduces Genetic Diversity on a Young Plant Y-Chromosome

X and Y chromosomes differ in effective population size (N_e ), rates of recombination, and exposure to natural selection, all of which can affect patterns of genetic diversity. On Y chromosomes with suppressed recombination, natural selection is expected to eliminate linked neutral variation, and lower the N_e of Y compared to X chromosomes or autosomes. However, female-biased sex ratios and high variance in male reproductive success can also reduce Y-linked N_e , making it difficult to infer the causes of low Y-diversity. Here, we investigate the factors affecting levels of polymorphism during sex chromosome evolution in the dioecious plant Rumexhastatulus (Polygonaceae). Strikingly, we find that neutral diversity for genes on the Y chromosome is, on average, 2.1% of the value for their X-linked homologs, corresponding to a chromosome-wide reduction of 93% compared to the standard neutral expectation. We demonstrate that the magnitude of this diversity loss is inconsistent with reduced male N_e caused by neutral processes. Instead, using forward simulations and estimates of the distribution of deleterious fitness effects, we show that Y chromosome diversity loss can be explained by purifying selection acting in aggregate over a large number of genetically linked sites. Simulations also suggest that our observed level of Y-diversity is consistent with the joint action of purifying and positive selection, but only for models in which there were fewer constrained sites than we empirically estimated. Given the relatively recent origin of R. hastatulus sex chromosomes, our results imply that Y-chromosome degeneration in the early stages may be largely driven by selective interference rather than by neutral genetic drift of silenced Y-linked genes.

Detecting past male-mediated expansions using the Y chromosome

Males and females display biological differences that lead to a higher variance of offspring number in males, and this is frequently exacerbated in human societies by mating practices, and possibly by past socio-cultural circumstances. This implies that the genetic record might contain the imprint of past male-mediated expansions, which can be investigated by analysing the male-specific region of the Y chromosome (MSY). Here, we review studies that have used MSY data to infer such expansions. Sets of short-tandem repeats define haplotypes of very low average frequencies, but in a few cases, high-frequency haplotypes are observed, forming the core of descent clusters. Estimates of the ages of such clusters, together with geographical information, have been used to propose powerful historical founders, including Genghis Khan, although without direct supporting evidence. Resequencing of multi-megabase segments of MSY has allowed the construction of detailed phylogenies in which branch lengths are proportional to time, leading to the identification of lineage expansions in the last few millennia as well as the more distant past. Comparisons with maternally-inherited mitochondrial DNA sequence data allow the male specificity of some of these expansions to be demonstrated. These include expansions in Europe in the last ~5000 years that may be associated with a cultural shift during the Bronze Age, as well as expansions elsewhere in the world for which explanations from archaeological evidence are not yet clear.

Fixation Probability in a Haploid-Diploid Population

Classical population genetic theory generally assumes either a fully haploid or fully diploid life cycle. However, many organisms exhibit more complex life cycles, with both free-living haploid and diploid stages. Here we ask what the probability of fixation is for selected alleles in organisms with haploid-diploid life cycles. We develop a genetic model that considers the population dynamics using both the Moran model and Wright-Fisher model. Applying a branching process approximation, we obtain an accurate fixation probability assuming that the population is large and the net effect of the mutation is beneficial. We also find the diffusion approximation for the fixation probability, which is accurate even in small populations and for deleterious alleles, as long as selection is weak. These fixation probabilities from branching process and diffusion approximations are similar when selection is weak for beneficial mutations that are not fully recessive. In many cases, particularly when one phase predominates, the fixation probability differs substantially for haploid-diploid organisms compared to either fully haploid or diploid species.

Prediction and estimation of effective population size

Effective population size (Ne) is a key parameter in population genetics. It has important applications in evolutionary biology, conservation genetics and plant and animal breeding, because it measures the rates of genetic drift and inbreeding and affects the efficacy of systematic evolutionary forces, such as mutation, selection and migration. We review the developments in predictive equations and estimation methodologies of effective size. In the prediction part, we focus on the equations for populations with different modes of reproduction, for populations under selection for unlinked or linked loci and for the specific applications to conservation genetics. In the estimation part, we focus on methods developed for estimating the current or recent effective size from molecular marker or sequence data. We discuss some underdeveloped areas in predicting and estimating Ne for future research.

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.

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.

Positive Selection Underlies Faster-Z Evolution of Gene Expression in Birds

The elevated rate of evolution for genes on sex chromosomes compared with autosomes (Fast-X or Fast-Z evolution) can result either from positive selection in the heterogametic sex or from nonadaptive consequences of reduced relative effective population size. Recent work in birds suggests that Fast-Z of coding sequence is primarily due to relaxed purifying selection resulting from reduced relative effective population size. However, gene sequence and gene expression are often subject to distinct evolutionary pressures; therefore, we tested for Fast-Z in gene expression using next-generation RNA-sequencing data from multiple avian species. Similar to studies of Fast-Z in coding sequence, we recover clear signatures of Fast-Z in gene expression; however, in contrast to coding sequence, our data indicate that Fast-Z in expression is due to positive selection acting primarily in females. In the soma, where gene expression is highly correlated between the sexes, we detected Fast-Z in both sexes, although at a higher rate in females, suggesting that many positively selected expression changes in females are also expressed in males. In the gonad, where intersexual correlations in expression are much lower, we detected Fast-Z for female gene expression, but crucially, not males. This suggests that a large amount of expression variation is sex-specific in its effects within the gonad. Taken together, our results indicate that Fast-Z evolution of gene expression is the product of positive selection acting on recessive beneficial alleles in the heterogametic sex. More broadly, our analysis suggests that the adaptive potential of Z chromosome gene expression may be much greater than that of gene sequence, results which have important implications for the role of sex chromosomes in speciation and sexual selection.

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.

Assessing effective population size, coancestry and inbreeding effects on litter size using the pedigree and SNP data in closed lines of the Iberian pig breed

The complete pedigree of two closed Iberian pig lines (Gamito and Torbiscal), with 798 and 4077 reproducers, has been used to measure the evolution of coancestry (f) and inbreeding (F) for autosomal and X-linked genes along 16 and 28 respective equivalent discrete generations. At the last generation, the mean values of each line were f = 0.41 and 0.22, F = 0.35 and 0.18, fX  = 0.46 and 0.22 and FX  = 0.47 and 0.19, respectively. Other calculated parameters were the effective number of founders (final values, 6.8 and 35.2) and non-founders (1.5 and 2.4), founder genome equivalents (1.2 and 2.3) and effective population size (16.0 and 57.7). Measures of Torbiscal effective size based on rates of coancestry (66.1), inbreeding (65.0) and linkage disequilibrium (71.0) were estimated from whole-genome SNP genotyping data. Values of new and old inbreeding and their respective rates by generation were computed to detect purging effects of natural selection. The analysis of 6854 Torbiscal litters showed significant negative impacts of new and fast inbreeding on litter size, as expected from the purging hypothesis: -0.20 born piglets per litter by a 10% of new inbreeding, and -0.03 and -0.02 piglets by 1% of total and new inbreeding rates, respectively. The analysis performed on 1274 litters of the Gamito line failed to show purging effects. The only significant results were reductions in -0.91 and -0.17 piglets by a 10% of old and X-linked genes inbreeding, respectively. These results may be useful for some practical issues in conservation programs of farm or captive wild animals.

Molecular Mechanisms and Evolutionary Processes Contributing to Accelerated Divergence of Gene Expression on the Drosophila X Chromosome

In species with a heterogametic sex, population genetics theory predicts that DNA sequences on the X chromosome can evolve faster than comparable sequences on autosomes. Both neutral and nonneutral evolutionary processes can generate this pattern. Complex traits like gene expression are not predicted to have accelerated evolution by these theories, yet a "faster-X" pattern of gene expression divergence has recently been reported for both Drosophila and mammals. Here, we test the hypothesis that accelerated adaptive evolution of cis-regulatory sequences on the X chromosome is responsible for this pattern by comparing the relative contributions of cis- and trans-regulatory changes to patterns of faster-X expression divergence observed between strains and species of Drosophila with a range of divergence times. We find support for this hypothesis, especially among male-biased genes, when comparing different species. However, we also find evidence that trans-regulatory differences contribute to a faster-X pattern of expression divergence both within and between species. This contribution is surprising because trans-acting regulators of X-linked genes are generally assumed to be randomly distributed throughout the genome. We found, however, that X-linked transcription factors appear to preferentially regulate expression of X-linked genes, providing a potential mechanistic explanation for this result. The contribution of trans-regulatory variation to faster-X expression divergence was larger within than between species, suggesting that it is more likely to result from neutral processes than positive selection. These data show how accelerated evolution of both coding and noncoding sequences on the X chromosome can lead to accelerated expression divergence on the X chromosome relative to autosomes.

Metapopulation inbreeding dynamics, effective size and subpopulation differentiation--A general analytical approach for diploid organisms

Motivated by problems in conservation biology we study genetic dynamics in structured populations of diploid organisms (monoecious or dioecious). Our analysis provides an analytical framework that unifies substantial parts of previous work in terms of exact identity by descent (IBD) and identity by state (IBS) recursions. We provide exact conditions under which two structured haploid and diploid populations are equivalent, and some sufficient conditions under which a dioecious diploid population can be treated as a monoecious diploid one. The IBD recursions are used for computing local and metapopulation inbreeding and coancestry effective population sizes and for predictions of several types of fixation indices over different time horizons.

A new general analytical approach for modeling patterns of genetic differentiation and effective size of subdivided populations over time

The main purpose of this paper is to develop a theoretical framework for assessing effective population size and genetic divergence in situations with structured populations that consist of various numbers of more or less interconnected subpopulations. We introduce a general infinite allele model for a diploid, monoecious and subdivided population, with subpopulation sizes varying over time, including local subpopulation extinction and recolonization, bottlenecks, cyclic census size changes or exponential growth. Exact matrix analytic formulas are derived for recursions of predicted (expected) gene identities and gene diversities, identity by descent and coalescence probabilities, and standardized variances of allele frequency change. This enables us to compute and put into a general framework a number of different types of genetically effective population sizes (Ne) including variance, inbreeding, nucleotide diversity, and eigenvalue effective size. General expressions for predictions (gST) of the coefficient of gene differentiation GST are also derived. We suggest that in order to adequately describe important properties of a subdivided population with respect to allele frequency change and maintenance of genetic variation over time, single values of gST and Ne are not enough. Rather, the temporal dynamic patterns of these properties are important to consider. We introduce several schemes for weighting subpopulations that enable effective size and expected genetic divergence to be calculated and described as functions of time, globally for the whole population and locally for any group of subpopulations. The traditional concept of effective size is generalized to situations where genetic drift is confounded by external sources, such as immigration and mutation. Finally, we introduce a general methodology for state space reduction, which greatly decreases the computational complexity of the matrix analytic formulas.

On the eigenvalue effective size of structured populations

A general theory is developed for the eigenvalue effective size (N(e)E) of structured populations in which a gene with two alleles segregates in discrete time. Generalizing results of Ewens (Theor Popul Biol 21:373-378, 1982), we characterize N(e)E in terms of the largest non-unit eigenvalue of the transition matrix of a Markov chain of allele frequencies. We use Perron-Frobenius Theorem to prove that the same eigenvalue appears in a linear recursion of predicted gene diversities between all pairs of subpopulations. Coalescence theory is employed in order to characterize this recursion, so that explicit novel expressions for N(e)E can be derived. We then study N(e)E asymptotically, when either the inverse size and/or the overall migration rate between subpopulations tend to zero. It is demonstrated that several previously known results can be deduced as special cases. In particular when the coalescence effective size N(e)C exists, it is an asymptotic version of N(e)E in the limit of large populations.

Natural selection reduced diversity on human y chromosomes

The human Y chromosome exhibits surprisingly low levels of genetic diversity. This could result from neutral processes if the effective population size of males is reduced relative to females due to a higher variance in the number of offspring from males than from females. Alternatively, selection acting on new mutations, and affecting linked neutral sites, could reduce variability on the Y chromosome. Here, using genome-wide analyses of X, Y, autosomal and mitochondrial DNA, in combination with extensive population genetic simulations, we show that low observed Y chromosome variability is not consistent with a purely neutral model. Instead, we show that models of purifying selection are consistent with observed Y diversity. Further, the number of sites estimated to be under purifying selection greatly exceeds the number of Y-linked coding sites, suggesting the importance of the highly repetitive ampliconic regions. While we show that purifying selection removing deleterious mutations can explain the low diversity on the Y chromosome, we cannot exclude the possibility that positive selection acting on beneficial mutations could have also reduced diversity in linked neutral regions, and may have contributed to lowering human Y chromosome diversity. Because the functional significance of the ampliconic regions is poorly understood, our findings should motivate future research in this area.

The human Y chromosome is found only in males, and exhibits surprisingly low levels of genetic diversity. This low diversity could result from neutral processes, for example, if there are fewer males successfully mating (and thus fewer Y chromosomes being inherited) relative to the number of females who successfully mate. Alternatively, natural selection may act on mutations on the Y chromosome to reduce genetic diversity. Because there is no recombination across most of the Y chromosome all sites on the Y are effectively linked together. Thus, selection acting on any one site will affect all sites on the Y indirectly. Here, studying the X, Y, autosomal and mitochondrial DNA, in combination with population genetic simulations, we show that low observed Y chromosome variability is consistent with models of purifying selection removing deleterious mutations and linked variation, although positive selection may also be acting. We further infer that the number of sites affected by selection likely includes some proportion of the highly repetitive ampliconic regions on the Y. Because the functional significance of the ampliconic regions is poorly understood, our findings should motivate future research in this area.

Comparative sex chromosome genomics in snakes: differentiation, evolutionary strata, and lack of global dosage compensation

Snakes exhibit genetic sex determination, with female heterogametic sex chromosomes (ZZ males, ZW females). Extensive cytogenetic work has suggested that the level of sex chromosome heteromorphism varies among species, with Boidae having entirely homomorphic sex chromosomes, Viperidae having completely heteromorphic sex chromosomes, and Colubridae showing partial differentiation. Here, we take a genomic approach to compare sex chromosome differentiation in these three snake families. We identify homomorphic sex chromosomes in boas (Boidae), but completely heteromorphic sex chromosomes in both garter snakes (Colubridae) and pygmy rattlesnake (Viperidae). Detection of W-linked gametologs enables us to establish the presence of evolutionary strata on garter and pygmy rattlesnake sex chromosomes where recombination was abolished at different time points. Sequence analysis shows that all strata are shared between pygmy rattlesnake and garter snake, i.e., recombination was abolished between the sex chromosomes before the two lineages diverged. The sex-biased transmission of the Z and its hemizygosity in females can impact patterns of molecular evolution, and we show that rates of evolution for Z-linked genes are increased relative to their pseudoautosomal homologs, both at synonymous and amino acid sites (even after controlling for mutational biases). This demonstrates that mutation rates are male-biased in snakes (male-driven evolution), but also supports faster-Z evolution due to differential selective effects on the Z. Finally, we perform a transcriptome analysis in boa and pygmy rattlesnake to establish baseline levels of sex-biased expression in homomorphic sex chromosomes, and show that heteromorphic ZW chromosomes in rattlesnakes lack chromosome-wide dosage compensation. Our study provides the first full scale overview of the evolution of snake sex chromosomes at the genomic level, thus greatly expanding our knowledge of reptilian and vertebrate sex chromosomes evolution.

The effect of nonindependent mate pairing on the effective population size

The effective population size (N(e)) quantifies the effectiveness of genetic drift in finite populations. When generations overlap, theoretical expectations for N(e) typically assume that the sampling of offspring genotypes from a given individual is independent among successive breeding events, even though this is not true in many species, including humans. To explore the effects on N(e) of nonindependent mate pairing across breeding events, we simulated the genetic drift of mitochondrial DNA, autosomal DNA, and sex chromosome DNA under three mating systems. Nonindependent mate pairing across breeding seasons has no effect when all adults mate pair for life, a small or moderate effect when females reuse stored sperm, and a large effect when there is intense male-male competition for reproduction in a harem social system. If adult females reproduce at a constant rate irrespective of the type of mate pairing, the general effect of nonindependent mate pairing is to decrease N(e) for paternally inherited components of the genome. These findings have significant implications for the relative N(e) values of different genomic regions, and hence for the expected levels of DNA sequence diversity in these regions.

Genomic conflict drives patterns of X-linked population structure in Drosophila neotestacea

Intragenomic conflict has the potential to cause widespread changes in patterns of genetic diversity and genome evolution. In this study, we investigate the consequences of sex-ratio (SR) drive on the population genetic patterns of the X-chromosome in Drosophila neotestacea. An SR X-chromosome prevents the maturation of Y-bearing sperm during male spermatogenesis and thus is transmitted to ~100% of the offspring, nearly all of which are daughters. Selection on the rest of the genome to suppress SR can be strong, and the resulting conflict over the offspring sex ratio can result in the accumulation of multiple loci on the X-chromosome that are necessary for the expression of drive. We surveyed variation at 12 random X-linked microsatellites across 16 populations of D. neotestacea that range in SR frequency from 0% to 30%. First, every locus was differentiated between SR and wild-type chromosomes, and this drives genetic structure at the X-chromosome. Once the association with SR is accounted for, the patterns of differentiation among populations are similar to the autosomes. Second, within wild-type chromosomes, the relative heterozygosity is reduced in populations with an increased prevalence of drive, and the heterozygosity of SR chromosomes is higher than expected based on its prevalence. The combination of the relatively high prevalence of SR drive and the structuring of polymorphism between the SR and wild-type chromosomes suggests that genetic conflict because of SR drive has had significant consequences on the patterns of X-linked polymorphism and thus also probably affects the tempo of X-chromosome evolution in D. neotestacea.

Surprising differences in the variability of Y chromosomes in African and cosmopolitan populations of Drosophila melanogaster

The nonrecombining Drosophila melanogaster Y chromosome is heterochromatic and has few genes. Despite these limitations, there remains ample opportunity for natural selection to act on the genes that are vital for male fertility and on Y factors that modulate gene expression elsewhere in the genome. Y chromosomes of many organisms have low levels of nucleotide variability, but a formal survey of D. melanogaster Y chromosome variation had yet to be performed. Here we surveyed Y-linked variation in six populations of D. melanogaster spread across the globe. We find surprisingly low levels of variability in African relative to Cosmopolitan (i.e., non-African) populations. While the low levels of Cosmopolitan Y chromosome polymorphism can be explained by the demographic histories of these populations, the staggeringly low polymorphism of African Y chromosomes cannot be explained by demographic history. An explanation that is entirely consistent with the data is that the Y chromosomes of Zimbabwe and Uganda populations have experienced recent selective sweeps. Interestingly, the Zimbabwe and Uganda Y chromosomes differ: in Zimbabwe, a European Y chromosome appears to have swept through the population.

The effects of selection and genetic drift on the genomic distribution of sexually antagonistic alleles

Sexual antagonism (SA) occurs when an allele that is beneficial to one sex, is detrimental to the other. This conflict can result in balancing, directional, or disruptive selection acting on SA alleles. A body of theory predicts the conditions under which sexually antagonistic mutants will invade and be maintained in stable polymorphism under balancing selection. There remains, however, considerable debate over the distribution of SA genetic variation across autosomes and sex chromosomes, with contradictory evidence coming from data and theory. In this article, we investigate how the interplay between selection and genetic drift will affect the genomic distribution of sexually antagonistic alleles. The effective population sizes can differ between the autosomes and the sex chromosomes due to a number of ecological factors and, consequently, the distribution of SA genetic variation in genomes. In general, we predict the interplay of SA selection and genetic drift should lead to the accumulation of SA alleles on the X in male heterogametic (XY) species and, on the autosomes in female heterogametic (ZW) species, especially when sexual competition is strong among males.

Pedigree- and marker-based methods in the estimation of genetic diversity in small groups of Holstein cattle

Genetic diversity is often evaluated using pedigree information. Currently, diversity can be evaluated in more detail over the genome based on large numbers of SNP markers. Pedigree- and SNP-based diversity were compared for two small related groups of Holstein animals genotyped with the 50 k SNP chip, genome-wide, per chromosome and for part of the genome examined. Diversity was estimated with coefficient of kinship (pedigree) and expected heterozygosity (SNP). SNP-based diversity at chromosome regions was determined using 5-Mb sliding windows, and significance of difference between groups was determined by bootstrapping. Both pedigree- and SNP-based diversity indicated more diversity in one of the groups; 26 of the 30 chromosomes showed significantly more diversity for the same group, as did 25.9% of the chromosome regions. Even in small populations that are genetically close, differences in diversity can be detected. Pedigree- and SNP-based diversity give comparable differences, but SNP-based diversity shows on which chromosome regions these differences are based. For maintaining diversity in a gene bank, SNP-based diversity gives a more detailed picture than pedigree-based diversity.

Accelerated evolution of sex chromosomes in aphids, an x0 system

Sex chromosomes play a role in many important biological processes, including sex determination, genomic conflicts, imprinting, and speciation. In particular, they exhibit several unusual properties such as inheritance pattern, hemizygosity, and reduced recombination, which influence their response to evolutionary factors (e.g., drift, selection, and demography). Here, we examine the evolutionary forces driving X chromosome evolution in aphids, an XO system where females are homozygous (XX) and males are hemizygous (X0) at sex chromosomes. We show by simulations that the unusual mode of transmission of the X chromosome in aphids, coupled with cyclical parthenogenesis, results in similar effective population sizes and predicted levels of genetic diversity for X chromosomes and autosomes under neutral evolution. These results contrast with expectations from standard XX/XY or XX/X0 systems (where the effective population size of the X is three-fourths that of autosomes) and have deep consequences for aphid X chromosome evolution. We then localized 52 microsatellite markers on the X and 351 on autosomes. We genotyped 167 individuals with 356 of these loci and found similar levels of allelic richness on the X and on the autosomes, as predicted by our simulations. In contrast, we detected higher dN and dN/dS ratio for X-linked genes compared with autosomal genes, a pattern compatible with either positive or relaxed selection. Given that both types of chromosomes have similar effective population sizes and that the single copy of the X chromosome of male aphids exposes its recessive genes to selection, some degree of positive selection seems to best explain the higher rates of evolution of X-linked genes. Overall, this study highlights the particular relevance of aphids to study the evolutionary factors driving sex chromosomes and genome evolution.

Genetic diversity, population structure and Wolbachia infection status in a worldwide sample of Drosophila melanogaster and D. simulans populations

Drosophila melanogaster and its close relatives have been extremely important model species in the development of population genetic models that serve to explain patterns of diversity in natural populations, a major goal of evolutionary biology. A detailed picture of the evolutionary history of these species is beginning to emerge, as the relative importance of forces including demographic changes and natural selection is established. A continuing aim is to characterise levels of genetic diversity in a large number of populations of these species, covering a wide geographic area. We have used collections from five previously un-sampled wild populations of D. melanogaster and two of D. simulans, across three continents. We estimated levels of genetic diversity within, and divergence between, these populations, and looked for evidence of genetic structure both between ancestral and derived populations, and amongst derived populations. We also investigated the prevalence of infection with the bacterial endosymbiont Wolbachia. We found that D. melanogaster populations from Sub-Saharan Africa are the most diverse, and that divergence is highest between these and non-Sub-Saharan populations. There is strong evidence for structuring of populations between Sub-Saharan Africa and the rest of the world, and some evidence for weak structure amongst derived populations. Populations from Sub-Saharan Africa also differ in the prevalence of Wolbachia infection, with very low levels of infection compared to populations from the rest of the world.

Coalescent theory for age-structured random mating populations with two sexes

Consider a population that develops over units of time labelled by zero and the negative integers. It is assumed that at any time r≤0 there are respectively N(m,a) males of age a and N(f,b) of females of age b, where a=1,…,A(m) and b=1,…,A(f). At time 0 a sample of n copies of a gene are assumed to be observed, where n≪min(ab)(N(m,a),N(f,b)). It is assumed that at any particular time r any possible mating is equally probable and that numbers of gametes contributed to offspring of age 1 and sex u by parents of sex v are exchangeable within age groups and independently distributed among age-groups. Coalescent theory is then derived, with time measured in multiples 2N(eC) of the effective population size N(eC), which depends on a measure T of the generation interval. Theory is developed for both autosomal and sex-linked loci in two special cases.

Contrasting population genetic patterns within the white-throated sparrow genome (Zonotrichia albicollis)

Background

The level of nucleotide diversity observed across the genome is positively correlated with the local rate of recombination. Avian karyotypes are typified by large variation in chromosome size and the rate of recombination in birds has been shown to be negatively correlated with chromosome size. It has thus been predicted that nucleotide diversity is negatively correlated with chromosome size in aves. However, there is limited empirical evidence to support this prediction.

Results

Here we sequenced 27 autosomal and 12 sex chromosome-linked loci in the white-throated sparrow (Zonotrichia albicollis) to quantify and compare patterns of recombination, linkage disequilibrium (LD), and genetic diversity across the genome of this North American songbird. Genetic diversity on the autosomes varied up to 8-fold, with the lowest diversity observed on the macrochromosomes and the highest diversity on the microchromosomes. Genetic diversity on the sex chromosomes was reduced compared to the autosomes, the most extreme difference being a ~300-fold difference between the W chromosome and the microchromosomes. LD and population structure associated with a common inversion polymorphism (ZAL2/2^m) in this species were found to be atypical compared to other macrochromosomes, and nucleotide diversity within this inversion on the two chromosome arrangements was more similar to that observed on the Z chromosome.

Conclusions

A negative correlation between nucleotide diversity and autosome size was observed in the white-throated sparrow genome, as well as low levels of diversity on the sex chromosomes comparable to those reported in other birds. The population structure and extended LD associated with the ZAL2/2^m chromosomal polymorphism are exceptional compared to the rest of the white-throated sparrow genome.

Can a sex-biased human demography account for the reduced effective population size of chromosome X in non-Africans?

Sex-biased demographic events can result in asymmetries in female and male effective population size that can lead to different patterns of genetic variation on chromosome X than are expected based on the patterns on the autosomes. Previous studies point to a period around the time of the dispersal of anatomically modern humans out of Africa when chromosome X experienced a significant reduction in effective population size relative to the autosomes. Here, we explore whether a sex-biased demographic history could explain these observations. We use coalescent simulations to show that a model of primarily male migration during the out-of-Africa dispersal can produce the striking patterns that are observed when comparing patterns of genetic variation on the autosomes and chromosome X. The model involves a history in which after the founder population of non-Africans lost much of its genetic diversity, subsequent mostly male gene flow from an African source brought new diversity into the population. We also explore two additional models, one of sex-biased generation time and one of a substructured population during the dispersal out of Africa with primarily female migration among demes. These latter models cannot account for the magnitude of the observed reduction in chromosome X effective population size, although it is plausible that they played a more minor role in producing the striking chromosome X/autosome patterns.

Chromosome-wide linkage disequilibrium caused by an inversion polymorphism in the white-throated sparrow (Zonotrichia albicollis)

Chromosomal inversions have been of long-standing interest to geneticists because they are capable of suppressing recombination and facilitating the formation of adaptive gene complexes. An exceptional inversion polymorphism (ZAL2(m)) in the white-throated sparrow (Zonotrichia albicollis) is linked to variation in plumage, social behavior and mate choice, and is maintained in the population by negative assortative mating. The ZAL2(m) polymorphism is a complex inversion spanning > 100 Mb and has been proposed to be a strong suppressor of recombination, as well as a potential model for studying neo-sex chromosome evolution. To quantify and evaluate these features of the ZAL2(m) polymorphism, we generated sequence from 8 ZAL2(m) and 16 ZAL2 chromosomes at 58 loci inside and 4 loci outside the inversion. Inside the inversion we found that recombination was completely suppressed between ZAL2 and ZAL2(m), resulting in uniformly high levels of genetic differentiation (F(ST)=0.94), the formation of two distinct haplotype groups representing the alternate chromosome arrangements and extensive linkage disequilibrium spanning ~104 Mb within the inversion, whereas gene flow was not suppressed outside the inversion. Finally, although ZAL2(m) homozygotes are exceedingly rare in the population, occurring at a frequency of < 1%, we detected evidence of historical recombination between ZAL2(m) chromosomes inside the inversion, refuting its potential status as a non-recombining autosome.