Fast-X on the Z: rapid evolution of sex-linked genes in birds

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.

Gene gain and loss from the Asian corn borer W chromosome

BACKGROUND

Sex-limited chromosomes Y and W share some characteristics, including the degeneration of protein-coding genes, enrichment of repetitive elements, and heterochromatin. However, although many studies have suggested that Y chromosomes retain genes related to male function, far less is known about W chromosomes and whether they retain genes related to female-specific function.

RESULTS

Here, we built a chromosome-level genome assembly of the Asian corn borer, Ostrinia furnacalis Guenée (Lepidoptera: Crambidae, Pyraloidea), an economically important pest in corn, from a female, including both the Z and W chromosome. Despite deep conservation of the Z chromosome across Lepidoptera, our chromosome-level W assembly reveals little conservation with available W chromosome sequence in related species or with the Z chromosome, consistent with a non-canonical origin of the W chromosome. The W chromosome has accumulated significant repetitive elements and experienced rapid gene gain from the remainder of the genome, with most genes exhibiting pseudogenization after duplication to the W. The genes that retain significant expression are largely enriched for functions in DNA recombination, the nucleosome, chromatin, and DNA binding, likely related to meiotic and mitotic processes within the female gonad.

CONCLUSIONS

Overall, our chromosome-level genome assembly supports the non-canonical origin of the W chromosome in O. furnacalis, which experienced rapid gene gain and loss, with the retention of genes related to female-specific function.

Evolution and expression patterns of the neo-sex chromosomes of the crested ibis

Bird sex chromosomes play a unique role in sex-determination, and affect the sexual morphology and behavior of bird species. Core waterbirds, a major clade of birds, share the common characteristics of being sexually monomorphic and having lower levels of inter-sexual conflict, yet their sex chromosome evolution remains poorly understood. Here, by we analyse of a chromosome-level assembly of a female crested ibis (Nipponia nippon), a typical core waterbird. We identify neo-sex chromosomes resulting from fusion of microchromosomes with ancient sex chromosomes. These fusion events likely occurred following the divergence of Threskiornithidae and Ardeidae. The neo-W chromosome of the crested ibis exhibits the characteristics of slow degradation, which is reflected in its retention of abundant gametologous genes. Neo-W chromosome genes display an apparent ovary-biased gene expression, which is largely driven by genes that are retained on the crested ibis W chromosome but lost in other bird species. These results provide new insights into the evolutionary history and expression patterns for the sex chromosomes of bird species.

Elasmobranch genome sequencing reveals evolutionary trends of vertebrate karyotype organization

Genomic studies of vertebrate chromosome evolution have long been hindered by the scarcity of chromosome-scale DNA sequences of some key taxa. One of those limiting taxa has been the elasmobranchs (sharks and rays), which harbor species often with numerous chromosomes and enlarged genomes. Here, we report the chromosome-scale genome assembly for the zebra shark Stegostoma tigrinum, an endangered species that has a relatively small genome among sharks (3.71 Gb), as well as for the whale shark Rhincodon typus Our analysis using a male-female comparison identified an X Chromosome, the first genomically characterized shark sex chromosome. The X Chromosome harbors the Hox C cluster whose intact linkage has not been shown for an elasmobranch fish. The sequenced shark genomes show a gradualism of chromosome length with remarkable length-dependent characteristics-shorter chromosomes tend to have higher GC content, gene density, synonymous substitution rate, and simple tandem repeat content as well as smaller gene length and lower interspersed repeat content. We challenge the traditional binary classification of karyotypes as with and without so-called microchromosomes. Even without microchromosomes, the length-dependent characteristics persist widely in nonmammalian vertebrates. Our investigation of elasmobranch karyotypes underpins their unique characteristics and provides clues for understanding how vertebrate karyotypes accommodate intragenomic heterogeneity to realize a complex readout. It also paves the way to dissecting more genomes with variable sizes to be sequenced at high quality.

Chromosome-Aware Phylogenomics of Assassin Bugs (Hemiptera: Reduvioidea) Elucidates Ancient Gene Conflict

Though the phylogenetic signal of loci on sex chromosomes can differ from those on autosomes, chromosomal-level genome assemblies for nonvertebrates are still relatively scarce and conservation of chromosomal gene content across deep phylogenetic scales has therefore remained largely unexplored. We here assemble a uniquely large and diverse set of samples (17 anchored hybrid enrichment, 24 RNA-seq, and 70 whole-genome sequencing samples of variable depth) for the medically important assassin bugs (Reduvioidea). We assess the performance of genes based on multiple features (e.g., nucleotide vs. amino acid, nuclear vs. mitochondrial, and autosomal vs. X chromosomal) and employ different methods (concatenation and coalescence analyses) to reconstruct the unresolved phylogeny of this diverse (∼7,000 spp.) and old (>180 Ma) group. Our results show that genes on the X chromosome are more likely to have discordant phylogenies than those on autosomes. We find that the X chromosome conflict is driven by high gene substitution rates that impact the accuracy of phylogenetic inference. However, gene tree clustering showed strong conflict even after discounting variable third codon positions. Alternative topologies were not particularly enriched for sex chromosome loci, but spread across the genome. We conclude that binning genes to autosomal or sex chromosomes may result in a more accurate picture of the complex evolutionary history of a clade.

Sex chromosome heteromorphism and the Fast-X effect in poeciliids

Fast-X evolution has been observed in a range of heteromorphic sex chromosomes. However, it remains unclear how early in the process of sex chromosome differentiation the Fast-X effect becomes detectible. Recently, we uncovered an extreme variation in sex chromosome heteromorphism across poeciliid fish species. The common guppy, Poecilia reticulata, Endler's guppy, P. wingei, swamp guppy, P. picta and para guppy, P. parae, appear to share the same XY system and exhibit a remarkable range of heteromorphism. Species outside this group lack this sex chromosome system. We combined analyses of sequence divergence and polymorphism data across poeciliids to investigate X chromosome evolution as a function of hemizygosity and reveal the causes for Fast-X effects. Consistent with the extent of Y degeneration in each species, we detect higher rates of divergence on the X relative to autosomes, a signal of Fast-X evolution, in P. picta and P. parae, species with high levels of X hemizygosity in males. In P. reticulata, which exhibits largely homomorphic sex chromosomes and little evidence of hemizygosity, we observe no change in the rate of evolution of X-linked relative to autosomal genes. In P. wingei, the species with intermediate sex chromosome differentiation, we see an increase in the rate of nonsynonymous substitutions on the older stratum of divergence only. We also use our comparative approach to test for the time of origin of the sex chromosomes in this clade. Taken together, our study reveals an important role of hemizygosity in Fast-X evolution.

100 years of Haldane's rule

Haldane's rule is one of the 'two rules of speciation'. It states that if one sex is 'absent, rare or sterile' in a hybrid population, then that sex will be heterogametic. Since Haldane first made this observation, 100 years have passed and still questions arise over how many independent examples exist and what the underlying causes of Haldane's rule are. This review aims to examine research that has occurred over the last century. It seeks to do so by discussing possible causes of Haldane's rule, as well as gaps in the research of these causes that could be readily addressed today. After 100 years of research, it can be concluded that Haldane's rule is a complicated one, and much current knowledge has been accrued by studying the model organisms of speciation. This has led to the primacy of dominance theory and faster-male theory as explanations for Haldane's rule. However, some of the most interesting findings of the 21st century with regard to Haldane's rule have involved investigating a wider range of taxa emphasizing the need to continue using comparative methods, including ever more taxa as new cases are discovered.

Demography and linked selection interact to shape the genomic landscape of codistributed woodpeckers during the Ice Age

The influence of genetic drift on population dynamics during Pleistocene glacial cycles is well understood, but the role of selection in shaping patterns of genomic variation during these events is less explored. We resequenced whole genomes to investigate how demography and natural selection interact to generate the genomic landscapes of Downy and Hairy Woodpecker, species codistributed in previously glaciated North America. First, we explored the spatial and temporal patterns of genomic diversity produced by neutral evolution. Next, we tested (i) whether levels of nucleotide diversity along the genome are correlated with intrinsic genomic properties, such as recombination rate and gene density, and (ii) whether different demographic trajectories impacted the efficacy of selection. Our results revealed cycles of bottleneck and expansion, and genetic structure associated with glacial refugia. Nucleotide diversity varied widely along the genome, but this variation was highly correlated between the species, suggesting the presence of conserved genomic features. In both taxa, nucleotide diversity was positively correlated with recombination rate and negatively correlated with gene density, suggesting that linked selection played a role in reducing diversity. Despite strong fluctuations in effective population size, the maintenance of relatively large populations during glaciations may have facilitated selection. Under these conditions, we found evidence that the individual demographic trajectory of populations modulated linked selection, with purifying selection being more efficient in removing deleterious alleles in large populations. These results highlight that while genome-wide variation reflects the expected signature of demographic change during climatic perturbations, the interaction of multiple processes produces a predictable and highly heterogeneous genomic landscape.

Widespread genomic signatures of reproductive isolation and sex-specific selection in the Eastern Yellow Robin, Eopsaltria australis

How new species evolve is one of the most fundamental questions in biology. Population divergence, which may lead to speciation, may be occurring in the Eastern Yellow Robin, a common passerine that lives along the eastern coast of Australia. This species is composed of 2 parapatric lineages that have highly divergent mitochondrial DNA; however, similar levels of divergence have not been observed in the nuclear genome. Here we re-examine the nuclear genomes of these mitolineages to test potential mechanisms underlying the discordance between nuclear and mitochondrial divergence. We find that nuclear admixture occurs in a narrow hybrid zone, although the majority of markers across the genome show evidence of reproductive isolation between populations of opposing mitolineages. There is an 8 MB section of a previously identified putative neo-sex chromosome that is highly diverged between allopatric but not parapatric populations, which may be the result of a chromosomal inversion. The neo-sex chromosomal nature of this region, as well as the geographic patterns in which it exhibits divergence, suggest it is unlikely to be contributing to reproductive isolation through mitonuclear incompatibilities as reported in earlier studies. In addition, there are sex differences in the number of markers that are differentiated between populations of opposite mitolineages, with greater differentiation occurring in females, which are heterozygous, than males. These results suggest that, despite the absence of previously observed assortative mating, mitolineages of Eastern Yellow Robin experience at least some postzygotic isolation from each other, in a pattern consistent with Haldane’s Rule.

Avian introgression patterns are consistent with Haldane's Rule

According to Haldane’s Rule, the heterogametic sex will show the greatest fitness reduction in a hybrid cross. In birds, where sex is determined by a ZW system, female hybrids are expected to experience lower fitness compared to male hybrids. This pattern has indeed been observed in several bird groups, but it is unknown whether the generality of Haldane’s Rule also extends to the molecular level. First, given the lower fitness of female hybrids, we can expect maternally inherited loci (i.e., mitochondrial and W-linked loci) to show lower introgression rates than biparentally inherited loci (i.e., autosomal loci) in females. Second, the faster evolution of Z-linked loci compared to autosomal loci and the hemizygosity of the Z-chromosome in females might speed up the accumulation of incompatible alleles on this sex chromosome, resulting in lower introgression rates for Z-linked loci than for autosomal loci. I tested these expectations by conducting a literature review which focused on studies that directly quantified introgression rates for autosomal, sex-linked, and mitochondrial loci. Although most studies reported introgression rates in line with Haldane’s Rule, it remains important to validate these genetic patterns with estimates of hybrid fitness and supporting field observations to rule out alternative explanations. Genomic data provide exciting opportunities to obtain a more fine-grained picture of introgression rates across the genome, which can consequently be linked to ecological and behavioral observations, potentially leading to novel insights into the genetic mechanisms underpinning Haldane’s Rule.

Dynamic Patterns of Sex Chromosome Evolution in Neognath Birds: Many Independent Barriers to Recombination at the ATP5F1A Locus

Avian sex chromosomes evolved after the divergence of birds and crocodilians from their common ancestor, so they are younger than the better-studied chromosomes of mammals. It has long been recognized that there may have been several stages to the evolution of avian sex chromosomes. For example, the CHD1 undergoes recombination in paleognaths but not neognaths. Genome assemblies have suggested that there may be variation in the timing of barriers to recombination among Neognathae, but there remains little understanding of the extent of this variability. Here, we look at partial sequences of ATP5F1A, which is on the avian Z and W chromosomes. It is known that recombination of this gene has independently ceased in Galliformes, Anseriformes, and at least five neoavian orders, but whether there are other independent cessations of recombination among Neoaves is not understood. We analyzed a combination of data extracted from published chromosomal-level genomes with data collected using PCR and cloning to identify Z and W copies in 22 orders. Our results suggest that there may be at least 19 independent cessations of recombination within Neognathae, and 3 clades that may still be undergoing recombination (or have only recently ceased recombination). Analyses of ATP5F1A protein sequences revealed an increased amino acid substitution rate for W chromosome gametologs, suggesting relaxed purifying selection on the W chromosome. Supporting this hypothesis, we found that the increased substitution rate was particularly pronounced for buried residues, which are expected to be more strongly constrained by purifying selection. This highlights the dynamic nature of avian sex chromosomes, and that this level of variation among clades means they should be a good system to understand sex chromosome evolution.

The functional genetic architecture of egg-laying and live-bearing reproduction in common lizards

All amniotes reproduce either by egg-laying (oviparity), which is ancestral to vertebrates or by live-bearing (viviparity), which has evolved many times independently. However, the genetic basis of these parity modes has never been resolved and, consequently, its convergence across evolutionary scales is currently unknown. Here, we leveraged natural hybridizations between oviparous and viviparous common lizards (Zootoca vivipara) to describe the functional genes and genetic architecture of parity mode and its key traits, eggshell and gestation length, and compared our findings across vertebrates. In these lizards, parity trait genes were associated with progesterone-binding functions and enriched for tissue remodelling and immune system pathways. Viviparity involved more genes and complex gene networks than did oviparity. Angiogenesis, vascular endothelial growth and adrenoreceptor pathways were enriched in the viviparous female reproductive tissue, while pathways for transforming growth factor were enriched in the oviparous. Natural selection on these parity mode genes was evident genome-wide. Our comparison to seven independent origins of viviparity in mammals, squamates and fish showed that genes active in pregnancy were related to immunity, tissue remodelling and blood vessel generation. Therefore, our results suggest that pre-established regulatory networks are repeatedly recruited for viviparity and that these are shared at deep evolutionary scales.

Autosomal, sex-linked and mitochondrial loci resolve evolutionary relationships among wrens in the genus Campylorhynchus

Although there is general consensus that sampling of multiple genetic loci is critical in accurate reconstruction of species trees, the exact numbers and the best types of molecular markers remain an open question. In particular, the phylogenetic utility of sex-linked loci is underexplored. Here, we sample all species and 70% of the named diversity of the New World wren genus Campylorhynchus using sequences from 23 loci, to evaluate the effects of linkage on efficiency in recovering a well-supported tree for the group. At a tree-wide level, we found that most loci supported fewer than half the possible clades and that sex-linked loci produced similar resolution to slower-coalescing autosomal markers, controlling for locus length. By contrast, we did find evidence that linkage affected the efficiency of recovery of individual relationships; as few as two sex-linked loci were necessary to resolve a selection of clades with long to medium subtending branches, whereas 4-6 autosomal loci were necessary to achieve comparable results. These results support an expanded role for sampling of the avian Z chromosome in phylogenetic studies, including target enrichment approaches. Our concatenated and species tree analyses represent significant improvements in our understanding of diversification in Campylorhynchus, and suggest a relatively complex scenario for its radiation across the Miocene/Pliocene boundary, with multiple invasions of South America.

Sex-Biased Gene Expression and Evolution in the Cerebrum and Syrinx of Chinese Hwamei (Garrulax canorus)

It is common that males and females display sexual dimorphisms, which usually result from sex-biased gene expression. Chinese hwamei (Garrulax canorus) is a good model for studying sex-biased gene expression because the song between the sexes is quite different. In this study, we analyze cerebrum and syrinx sex-biased gene expression and evolution using the de novo assembled Chinese hwamei transcriptome. In both the cerebrum and syrinx, our study revealed that most female-biased genes were actively expressed in females only, while most male-biased genes were actively expressed in both sexes. In addition, both male- and female-biased genes were enriched on the putative Z chromosome, suggesting the existence of sexually antagonistic genes and the insufficient dosage compensation of the Z-linked genes. We also identified a 9 Mb sex linkage region on the putative 4A chromosome which enriched more than 20% of female-biased genes. Resultantly, male-biased genes in both tissues had significantly higher Ka/Ks and effective number of codons (ENCs) than unbiased genes, and this suggested that male-biased genes which exhibit accelerated divergence may have resulted from positive selection. Taken together, our results initially revealed the reasons for the differences in singing behavior between males and females of Chinese hwamei.

Determinants of genome-wide distribution and evolution of uORFs in eukaryotes

Upstream open reading frames (uORFs) play widespread regulatory functions in modulating mRNA translation in eukaryotes, but the principles underlying the genomic distribution and evolution of uORFs remain poorly understood. Here, we analyze ~17 million putative canonical uORFs in 478 eukaryotic species that span most of the extant taxa of eukaryotes. We demonstrate how positive and purifying selection, coupled with differences in effective population size (Ne), has shaped the contents of uORFs in eukaryotes. Besides, gene expression level is important in influencing uORF occurrences across genes in a species. Our analyses suggest that most uORFs might play regulatory roles rather than encode functional peptides. We also show that the Kozak sequence context of uORFs has evolved across eukaryotic clades, and that noncanonical uORFs tend to have weaker suppressive effects than canonical uORFs in translation regulation. This study provides insights into the driving forces underlying uORF evolution in eukaryotes.

A Study of Faster-Z Evolution in the Great Tit (Parus major)

Sex chromosomes contribute substantially to key evolutionary processes such as speciation and adaptation. Several theories suggest that evolution could occur more rapidly on sex chromosomes, but currently our understanding of whether and how this occurs is limited. Here, we present an analysis of the great tit (Parus major) genome, aiming to detect signals of faster-Z evolution. We find mixed evidence of faster divergence on the Z chromosome than autosomes, with significantly higher divergence being found in ancestral repeats, but not at 4- or 0-fold degenerate sites. Interestingly, some 4-fold sites appear to be selectively constrained, which may mislead analyses that use these sites as the neutral reference (e.g., d_N/d_S). Consistent with other studies in birds, the mutation rate is significantly higher in males than females, and the long-term Z-to-autosome effective population size ratio is only 0.5, significantly lower than the expected value of 0.75. These are indicative of male-driven evolution and high variance in male reproductive success, respectively. We find no evidence for an increased efficacy of positive selection on the Z chromosome. In contrast, the Z chromosome in great tits appears to be affected by increased genetic drift, which has led to detectable signals of weakened intensity of purifying selection. These results provide further evidence that the Z chromosome often has a low effective population size, and that this has important consequences for its evolution. They also highlight the importance of considering multiple factors that can affect the rate of evolution and effective population sizes of sex chromosomes.

Transcriptomic and Epigenetic Preservation of Genetic Sex Identity in Estrogen-feminized Male Chicken Embryonic Gonads

Whereas in ovo exposure of genetically male (ZZ) chicken embryos to exogenous estrogens temporarily feminizes gonads at the time of hatching, the morphologically ovarian ZZ-gonads (FemZZs for feminized ZZ gonads) are masculinized back to testes within 1 year. To identify the feminization-resistant "memory" of genetic male sex, FemZZs showing varying degrees of feminization were subjected to transcriptomic, DNA methylome, and immunofluorescence analyses. Protein-coding genes were classified based on their relative mRNA expression across normal ZZ-testes, genetically female (ZW) ovaries, and FemZZs. We identified a group of 25 genes that were strongly expressed in both ZZ-testes and FemZZs but dramatically suppressed in ZW-ovaries. Interestingly, 84% (21/25) of these feminization-resistant testicular marker genes, including the DMRT1 master masculinizing gene, were located in chromosome Z. Expression of representative marker genes of germline cells (eg, DAZL or DDX4/VASA) was stronger in FemZZs than normal ZZ-testes or ZW-ovaries. We also identified 231 repetitive sequences (RSs) that were strongly expressed in both ZZ-testes and FemZZs, but these RSs were not enriched in chromosome Z. Although 94% (165/176) of RSs exclusively expressed in ZW-ovaries were located in chromosome W, no feminization-inducible RS was detected in FemZZs. DNA methylome analysis distinguished FemZZs from normal ZZ- and ZW-gonads. Immunofluorescence analysis of FemZZ gonads revealed expression of DMRT1 protein in medullary SOX9+ somatic cells and apparent germline cell populations in both medulla and cortex. Taken together, our study provides evidence that both somatic and germline cell populations in morphologically feminized FemZZs maintain significant transcriptomic and epigenetic memories of genetic sex.

Dense sampling of bird diversity increases power of comparative genomics

Whole-genome sequencing projects are increasingly populating the tree of life and characterizing biodiversity^1–4. Sparse taxon sampling has previously been proposed to confound phylogenetic inference⁵, and captures only a fraction of the genomic diversity. Here we report a substantial step towards the dense representation of avian phylogenetic and molecular diversity, by analysing 363 genomes from 92.4% of bird families—including 267 newly sequenced genomes produced for phase II of the Bird 10,000 Genomes (B10K) Project. We use this comparative genome dataset in combination with a pipeline that leverages a reference-free whole-genome alignment to identify orthologous regions in greater numbers than has previously been possible and to recognize genomic novelties in particular bird lineages. The densely sampled alignment provides a single-base-pair map of selection, has more than doubled the fraction of bases that are confidently predicted to be under conservation and reveals extensive patterns of weak selection in predominantly non-coding DNA. Our results demonstrate that increasing the diversity of genomes used in comparative studies can reveal more shared and lineage-specific variation, and improve the investigation of genomic characteristics. We anticipate that this genomic resource will offer new perspectives on evolutionary processes in cross-species comparative analyses and assist in efforts to conserve species.

A dataset of the genomes of 363 species from the Bird 10,000 Genomes Project shows increased power to detect shared and lineage-specific variation, demonstrating the importance of phylogenetically diverse taxon sampling in whole-genome sequencing.

Molecular sex identification of Malaysian White-Nest Swiftlet (Aerodramus fuciphagus Thunberg, 1812)

The difficulty in differentiating the sex of monomorphic bird species has made molecular sexing an important tool in addressing this problem. This method uses noninvasively collected materials such as feathers and may be advantageous for sexing endangered as well as commercialized bird species. In this study, seven primer sets for sexing birds were screened in Aerodramus fuciphagus using a total of 13 feather samples that were randomly selected from the state of Perak, Malaysia. From the screening analysis, only one primer set (P8/WZ/W) successfully differentiated the sex of A. fuciphagus. PCR amplification produced a single 255‐bp DNA fragment for males which was derived from CHD‐Z (CHD gene region in the sex chromosome Z), while for the females it produced two fragments (144 and 255 bp). The 144‐bp fragment was from CHD‐W (CHD gene region in the sex chromosome W). Results from sequencing showed no variations in the base sequences of the CHD‐W and CHD‐Z amplified fragments within the same sexes, except for one male sample (A23) where at position 166, a base substitution occurred (G → A). Phylogenetic analysis of CHD‐W showed that four (Apodiformes; Gruiformes; Passeriformes; and Pelecaniformes) out of the five orders investigated had formed four clear clusters within their orders, including the studied order: Apodiformes. Whereas in CHD‐Z, four (Accipitriformes; Columbiformes; Galliformes; and Passeriformes) out of five orders investigated formed four clear clusters within their orders, excluding the studied order. In addition, A. fuciphagus and Apus apus (both Apodiformes) showed less divergence in CHD‐W than CHD‐Z (0% c.f. 9%). The result suggests that in A. fuciphagus, CHD gene evolution occurred at a higher rate in males (CHD‐Z) compared to females (CHD‐W). This finding may be useful for further studies on sex ratio and breeding management of A. fuciphagus.

Turtle Insights into the Evolution of the Reptilian Karyotype and the Genomic Architecture of Sex Determination

Sex chromosome evolution remains an evolutionary puzzle despite its importance in understanding sexual development and genome evolution. The seemingly random distribution of sex-determining systems in reptiles offers a unique opportunity to study sex chromosome evolution not afforded by mammals or birds. These reptilian systems derive from multiple transitions in sex determination, some independent, some convergent, that lead to the birth and death of sex chromosomes in various lineages. Here we focus on turtles, an emerging model group with growing genomic resources. We review karyotypic changes that accompanied the evolution of chromosomal systems of genotypic sex determination (GSD) in chelonians from systems under the control of environmental temperature (TSD). These transitions gave rise to 31 GSD species identified thus far (out of 101 turtles with known sex determination), 27 with a characterized sex chromosome system (13 of those karyotypically). These sex chromosomes are varied in terms of the ancestral autosome they co-opted and thus in their homology, as well as in their size (some are macro-, some are micro-chromosomes), heterogamety (some are XX/XY, some ZZ/ZW), dimorphism (some are virtually homomorphic, some heteromorphic with larger-X, larger W, or smaller-Y), age (the oldest system could be ~195 My old and the youngest < 25 My old). Combined, all data indicate that turtles follow some tenets of classic theoretical models of sex chromosome evolution while countering others. Finally, although the study of dosage compensation and molecular divergence of turtle sex chromosomes has lagged behind research on other aspects of their evolution, this gap is rapidly decreasing with the acceleration of ongoing research and growing genomic resources in this group.

Absence of a Faster-X Effect in Beetles (Tribolium, Coleoptera)

The faster-X effect, namely the rapid evolution of protein-coding genes on the X chromosome, has been widely reported in metazoans. However, the prevalence of this phenomenon across diverse systems and its potential causes remain largely unresolved. Analysis of sex-biased genes may elucidate its possible mechanisms: for example, in systems with X/Y males a more pronounced faster-X effect in male-biased genes than in female-biased or unbiased genes may suggest fixation of recessive beneficial mutations rather than genetic drift. Further, theory predicts that the faster-X effect should be promoted by X chromosome dosage compensation. Here, we asked whether we could detect a faster-X effect in genes of the beetle Tribolium castaneum (and T. freemani orthologs), which has X/Y sex-determination and heterogametic males. Our comparison of protein sequence divergence (dN/dS) on the X chromosome vs. autosomes indicated a rarely observed absence of a faster-X effect in this organism. Further, analyses of sex-biased gene expression revealed that the X chromosome was particularly highly enriched for ovary-biased genes, which evolved slowly. In addition, an evaluation of male X chromosome dosage compensation in the gonads and in non-gonadal somatic tissues indicated a striking lack of compensation in the testis. This under-expression in testis may limit fixation of recessive beneficial X-linked mutations in genes transcribed in these male sex organs. Taken together, these beetles provide an example of the absence of a faster-X effect on protein evolution in a metazoan, that may result from two plausible factors, strong constraint on abundant X-linked ovary-biased genes and a lack of gonadal dosage compensation.

Evolutionary Dynamics of Sex Chromosomes of Paleognathous Birds

Standard models of sex chromosome evolution propose that recombination suppression leads to the degeneration of the heterogametic chromosome, as is seen for the Y chromosome in mammals and the W chromosome in most birds. Unlike other birds, paleognaths (ratites and tinamous) possess large nondegenerate regions on their sex chromosomes (PARs or pseudoautosomal regions). It remains unclear why these large PARs are retained over >100 Myr, and how this retention impacts the evolution of sex chromosomes within this system. To address this puzzle, we analyzed Z chromosome evolution and gene expression across 12 paleognaths, several of whose genomes have recently been sequenced. We confirm at the genomic level that most paleognaths retain large PARs. As in other birds, we find that all paleognaths have incomplete dosage compensation on the regions of the Z chromosome homologous to degenerated portions of the W (differentiated regions), but we find no evidence for enrichments of male-biased genes in PARs. We find limited evidence for increased evolutionary rates (faster-Z) either across the chromosome or in differentiated regions for most paleognaths with large PARs, but do recover signals of faster-Z evolution in tinamou species with mostly degenerated W chromosomes, similar to the pattern seen in neognaths. Unexpectedly, in some species, PAR-linked genes evolve faster on average than genes on autosomes, suggested by diverse genomic features to be due to reduced efficacy of selection in paleognath PARs. Our analysis shows that paleognath Z chromosomes are atypical at the genomic level, but the evolutionary forces maintaining largely homomorphic sex chromosomes in these species remain elusive.

Sex linkage of the skeletal muscle sodium channel gene (SCN4A) explains apparent deviations from Hardy-Weinberg equilibrium of tetrodotoxin-resistance alleles in garter snakes (Thamnophis sirtalis)

The arms race between tetrodotoxin-bearing Pacific newts (Taricha) and their garter snake predators (Thamnophis) in western North America has become a classic example of coevolution, shedding light on predator-prey dynamics, the molecular basis of adaptation, and patterns of convergent evolution. Newts are defended by tetrodotoxin (TTX), a neurotoxin that binds to voltage-gated sodium channels (Na_v proteins), arresting electrical activity in nerves and muscles and paralyzing would-be predators. However, populations of the common garter snake (T. sirtalis) have overcome this defense, largely through polymorphism at the locus SCN4A, which renders the encoded protein (Na_v1.4) less vulnerable to TTX. Previous work suggests that SCN4A commonly shows extreme deviations from Hardy-Weinberg equilibrium (HWE) in these populations, which has been interpreted as the result of intense selection imposed by newts. Here we show that much of this apparent deviation can be attributed to sex linkage of SCN4A. Using genomic data and quantitative PCR, we show that SCN4A is on the Z chromosome in Thamnophis and other advanced snakes. Taking Z-linkage into account, we find that most apparent deviations from HWE can be explained by female hemizygosity rather than low heterozygosity. Sex linkage can affect mutation rates, selection, and drift, and our results suggest that Z-linkage of SCN4A may make significant contributions to the overall dynamics of the coevolutionary arms race between newts and snakes.

Population differentiation and historical demography of the threatened snowy plover Charadrius nivosus (Cassin, 1858)

Delineating conservation units is a complex and often controversial process that is particularly challenging for highly vagile species. Here, we reassess population genetic structure and identify those populations of highest conservation value in the threatened snowy plover (Charadrius nivosus, Cassin, 1858), a partial migrant shorebird endemic to the Americas. We use four categories of genetic data—mitochondrial DNA (mtDNA), microsatellites, Z-linked and autosomal single nucleotide polymorphisms (SNPs)—to: (1) assess subspecies delineation and examine population structure (2) compare the sensitivity of the different types of genetic data to detect spatial genetic patterns, and (3) reconstruct demographic history of the populations analysed. Delineation of two traditionally recognised subspecies was broadly supported by all data. In addition, microsatellite and SNPs but not mtDNA supported the recognition of Caribbean snowy plovers (C. n. tenuirostris) and Floridian populations (eastern C. n. nivosus) as distinct genetic lineage and deme, respectively. Low migration rates estimated from autosomal SNPs (m < 0.03) reflect a general paucity of exchange between genetic lineages. In contrast, we detected strong unidirectional migration (m = 0.26) from the western into the eastern nivosus deme. Within western nivosus, we found no genetic differentiation between coastal Pacific and inland populations. The correlation between geographic and genetic distances was weak but significant for all genetic data sets. All demes showed signatures of bottlenecks occurring during the past 1000 years. We conclude that at least four snowy plover conservation units are warranted: in addition to subspecies nivosus and occidentalis, a third unit comprises the Caribbean tenuirostris lineage and a fourth unit the distinct eastern nivosus deme.

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.

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.

Evolution of an X-Linked miRNA Family Predominantly Expressed in Mammalian Male Germ Cells

MicroRNAs (miRNAs) are important posttranscriptional regulators of gene expression. However, comprehensive expression profiles of miRNAs during mammalian spermatogenesis are lacking. Herein, we sequenced small RNAs in highly purified mouse spermatogenic cells at different stages. We found that a family of X-linked miRNAs named spermatogenesis-related miRNAs (spermiRs) is predominantly expressed in the early meiotic phases and has a conserved testis-specific high expression pattern in different mammals. We identified one spermiR homolog in opossum; this homolog might originate from THER1, a retrotransposon that is active in marsupials but extinct in current placental mammals. SpermiRs have expanded rapidly with mammalian evolution and are diverged into two clades, spermiR-L and spermiR-R, which are likely to have been generated at least in part by tandem duplication mediated by flanking retrotransposable elements. Notably, despite having undergone highly frequent lineage-specific duplication events, the sequences encoding all spermiR family members are strictly located between two protein-coding genes, Slitrk2 and Fmr1. Moreover, spermiR-Ls and spermiR-Rs have evolved different expression patterns during spermatogenesis in different mammals. Intriguingly, the seed sequences of spermiRs, which are critical for the recognition of target genes, are highly divergent within and among mammals, whereas spermiR target genes largely overlap. When miR-741, the most highly expressed spermiR, is knocked out in cultured mouse spermatogonial stem cells (SSCs), another spermiR, miR-465a-5p, is dramatically upregulated and becomes the most abundant miRNA. Notably, miR-741^−/− SSCs grow normally, and the genome-wide expression levels of mRNAs remain unchanged. All these observations indicate functional compensation between spermiR family members and strong coevolution between spermiRs and their targets.

Support for the Dominance Theory in Drosophila Transcriptomes

Interactions among divergent elements of transcriptional networks from different species can lead to misexpression in hybrids through regulatory incompatibilities, some with the potential to generate sterility. While the possible contribution of faster-male evolution to this misexpression has been explored, the role of the hemizygous X chromosome (i.e., the dominance theory for transcriptomes) remains yet to be determined. Here, we study genome-wide patterns of gene expression in females and males of Drosophila yakuba, Drosophila santomea and their hybrids. We used attached-X stocks to specifically test the dominance theory, and we uncovered a significant contribution of recessive alleles on the X chromosome to hybrid misexpression. Our analyses also suggest a contribution of weakly deleterious regulatory mutations to gene expression divergence in genes with sex-biased expression, but only in the sex toward which the expression is biased (e.g., genes with female-biased expression when analyzed in females). In the opposite sex, we found stronger selective constraints on gene expression divergence. Although genes with a high degree of male-biased expression show a clear signal of faster-X evolution of gene expression, we also detected slower-X evolution in other gene classes (e.g., female-biased genes). This slower-X effect is mediated by significant decreases in cis- and trans-regulatory divergence. The distinct behavior of X-linked genes with a high degree of male-biased expression is consistent with these genes experiencing a higher incidence of positively selected regulatory mutations than their autosomal counterparts.

Chromosomal Context Affects the Molecular Evolution of Sex-linked Genes and Their Autosomal Counterparts in Turtles and Other Vertebrates

Sex chromosomes evolve differently from autosomes because natural selection acts distinctly on them given their reduced recombination and smaller population size. Various studies of sex-linked genes compared with different autosomal genes within species support these predictions. Here, we take a novel alternative approach by comparing the rate of evolution between subsets of genes that are sex-linked in selected reptiles/vertebrates and the same genes located in autosomes in other amniotes. We report for the first time the faster evolution of Z-linked genes in a turtle (the Chinese softshell turtle Pelodiscus sinensis) relative to autosomal orthologs in other taxa, including turtles with temperature-dependent sex determination (TSD). This faster rate was absent in its close relative, the spiny softshell turtle (Apalone spinifera), thus revealing important lineage effects, and was only surpassed by mammalian-X linked genes. In contrast, we found slower evolution of X-linked genes in the musk turtle Staurotypus triporcatus (XX/XY) and homologous Z-linked chicken genes. TSD lineages displayed overall faster sequence evolution than taxa with genotypic sex determination (GSD), ruling out global effects of GSD on molecular evolution beyond those by sex-linkage. Notably, results revealed a putative selective sweep around two turtle genes involved in vertebrate gonadogenesis (Pelodiscus-Z-linked Nf2 and Chrysemys-autosomal Tspan7). Our observations reveal important evolutionary changes at the gene level mediated by chromosomal context in turtles despite their low overall evolutionary rate and illuminate sex chromosome evolution by empirically testing expectations from theoretical models. Genome-wide analyses are warranted to test the generality and prevalence of the observed patterns.

Evolution of sex-specific pace-of-life syndromes: genetic architecture and physiological mechanisms

Sex differences in life history, physiology, and behavior are nearly ubiquitous across taxa, owing to sex-specific selection that arises from different reproductive strategies of the sexes. The pace-of-life syndrome (POLS) hypothesis predicts that most variation in such traits among individuals, populations, and species falls along a slow-fast pace-of-life continuum. As a result of their different reproductive roles and environment, the sexes also commonly differ in pace-of-life, with important consequences for the evolution of POLS. Here, we outline mechanisms for how males and females can evolve differences in POLS traits and in how such traits can covary differently despite constraints resulting from a shared genome. We review the current knowledge of the genetic basis of POLS traits and suggest candidate genes and pathways for future studies. Pleiotropic effects may govern many of the genetic correlations, but little is still known about the mechanisms involved in trade-offs between current and future reproduction and their integration with behavioral variation. We highlight the importance of metabolic and hormonal pathways in mediating sex differences in POLS traits; however, there is still a shortage of studies that test for sex specificity in molecular effects and their evolutionary causes. Considering whether and how sexual dimorphism evolves in POLS traits provides a more holistic framework to understand how behavioral variation is integrated with life histories and physiology, and we call for studies that focus on examining the sex-specific genetic architecture of this integration.

Recent Sex Chromosome Divergence despite Ancient Dioecy in the Willow Salix viminalis

Sex chromosomes can evolve when recombination is halted between a pair of chromosomes, and this can lead to degeneration of the sex-limited chromosome. In the early stages of differentiation sex chromosomes are homomorphic, and even though homomorphic sex chromosomes are very common throughout animals and plants, we know little about the evolutionary forces shaping these types of sex chromosomes. We used DNA- and RNA-Seq data from females and males to explore the sex chromosomes in the female heterogametic willow, Salix viminalis, a species with ancient dioecy but with homomorphic sex chromosomes. We detected no major sex differences in read coverage in the sex determination (SD) region, indicating that the W region has not significantly degenerated. However, single nucleotide polymorphism densities in the SD region are higher in females compared with males, indicating very recent recombination suppression, followed by the accumulation of sex-specific single nucleotide polymorphisms. Interestingly, we identified two female-specific scaffolds that likely represent W-chromosome-specific sequence. We show that genes located in the SD region display a mild excess of male-biased expression in sex-specific tissue, and we use allele-specific gene expression analysis to show that this is the result of masculinization of expression on the Z chromosome rather than degeneration of female-expression on the W chromosome. Together, our results demonstrate that insertion of small DNA fragments and accumulation of sex-biased gene expression can occur before the detectable decay of the sex-limited chromosome.

Evolution of Dosage Compensation in Anolis carolinensis, a Reptile with XX/XY Chromosomal Sex Determination

In species with highly heteromorphic sex chromosomes, the degradation of one of the sex chromosomes will result in unequal gene expression between the sexes (e.g. between XX females and XY males) and between the sex chromosomes and the autosomes. Dosage compensation is a process whereby genes on the sex chromosomes achieve equal gene expression. We compared genome-wide levels of transcription between males and females, and between the X chromosome and the autosomes in the green anole, Anolis carolinensis. We present evidence for dosage compensation between the sexes, and between the sex chromosomes and the autosomes. When dividing the X chromosome into regions based on linkage groups, we discovered that genes in the first reported X-linked region, anole linkage group b (LGb), exhibit complete dosage compensation, although the rest of the X-linked genes exhibit incomplete dosage compensation. Our data further suggest that the mechanism of this dosage compensation is upregulation of the X chromosome in males. We report that approximately 10% of coding genes, most of which are on the autosomes, are differentially expressed between males and females. In addition, genes on the X chromosome exhibited higher ratios of nonsynonymous to synonymous substitution than autosomal genes, consistent with the fast-X effect. Our results from the green anole add an additional observation of dosage compensation in a species with XX/XY sex determination.

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.

Heterogeneous Patterns of Genetic Diversity and Differentiation in European and Siberian Chiffchaff (Phylloscopus collybita abietinus/P. tristis)

Identification of candidate genes for trait variation in diverging lineages and characterization of mechanistic underpinnings of genome differentiation are key steps toward understanding the processes underlying the formation of new species. Hybrid zones provide a valuable resource for such investigations, since they allow us to study how genomes evolve as species exchange genetic material and to associate particular genetic regions with phenotypic traits of interest. Here, we use whole-genome resequencing of both allopatric and hybridizing populations of the European (Phylloscopus collybita abietinus) and the Siberian chiffchaff (P. tristis)—two recently diverged species which differ in morphology, plumage, song, habitat, and migration—to quantify the regional variation in genome-wide genetic diversity and differentiation, and to identify candidate regions for trait variation. We find that the levels of diversity, differentiation, and divergence are highly heterogeneous, with significantly reduced global differentiation, and more pronounced differentiation peaks in sympatry than in allopatry. This pattern is consistent with regional differences in effective population size and recurrent background selection or selective sweeps reducing the genetic diversity in specific regions prior to lineage divergence, but the data also suggest that postdivergence selection has resulted in increased differentiation and fixed differences in specific regions. We find that hybridization and backcrossing is common in sympatry, and that phenotype is a poor predictor of the genomic composition of sympatric birds. The combination of a differentiation scan approach with identification of fixed differences pinpoint a handful of candidate regions that might be important for trait variation between the two species.

Global Dosage Compensation Is Ubiquitous in Lepidoptera, but Counteracted by the Masculinization of the Z Chromosome

While chromosome-wide dosage compensation of the X chromosome has been found in many species, studies in ZW clades have indicated that compensation of the Z is more localized and/or incomplete. In the ZW Lepidoptera, some species show complete compensation of the Z chromosome, while others lack full equalization, but what drives these inconsistencies is unclear. Here, we compare patterns of male and female gene expression on the Z chromosome of two closely related butterfly species, Papilio xuthus and Papilio machaon, and in multiple tissues of two moths species, Plodia interpunctella and Bombyx mori, which were previously found to differ in the extent to which they equalize Z-linked gene expression between the sexes. We find that, while some species and tissues seem to have incomplete dosage compensation, this is in fact due to the accumulation of male-biased genes and the depletion of female-biased genes on the Z chromosome. Once this is accounted for, the Z chromosome is fully compensated in all four species, through the up-regulation of Z expression in females and in some cases additional down-regulation in males. We further find that both sex-biased genes and Z-linked genes have increased rates of expression divergence in this clade, and that this can lead to fast shifts in patterns of gene expression even between closely related species. Taken together, these results show that the uneven distribution of sex-biased genes on sex chromosomes can confound conclusions about dosage compensation and that Z chromosome-wide dosage compensation is not only possible but ubiquitous among Lepidoptera.

The genomic mosaicism of hybrid speciation

Hybridization is widespread in nature and, in some instances, can result in the formation of a new hybrid species. We investigate the genetic foundation of this poorly understood process through whole-genome analysis of the hybrid Italian sparrow and its progenitors. We find overall balanced yet heterogeneous levels of contribution from each parent species throughout the hybrid genome and identify areas of novel divergence in the hybrid species exhibiting signals consistent with balancing selection. High-divergence areas are disproportionately located on the Z chromosome and overrepresented in gene networks relating to key traits separating the focal species, which are likely involved in reproductive barriers and/or species-specific adaptations. Of special interest are genes and functional groups known to affect body patterning, beak morphology, and the immune system, which are important features of diversification and fitness. We show that a combination of mosaic parental inheritance and novel divergence within the hybrid lineage has facilitated the origin and maintenance of an avian hybrid species.

Convergent recombination suppression suggests role of sexual selection in guppy sex chromosome formation

Sex chromosomes evolve once recombination is halted between a homologous pair of chromosomes. The dominant model of sex chromosome evolution posits that recombination is suppressed between emerging X and Y chromosomes in order to resolve sexual conflict. Here we test this model using whole genome and transcriptome resequencing data in the guppy, a model for sexual selection with many Y-linked colour traits. We show that although the nascent Y chromosome encompasses nearly half of the linkage group, there has been no perceptible degradation of Y chromosome gene content or activity. Using replicate wild populations with differing levels of sexually antagonistic selection for colour, we also show that sexual selection leads to greater expansion of the non-recombining region and increased Y chromosome divergence. These results provide empirical support for longstanding models of sex chromosome catalysis, and suggest an important role for sexual selection and sexual conflict in genome evolution.

Interspecific hybrids of dwarf hamsters and Phasianidae birds as animal models for studying the genetic and developmental basis of hybrid incompatibility

Hybrid incompatibility is important in speciation as it prevents gene flow between closely related populations. Reduced fitness from hybrid incompatibility may also reinforce prezygotic reproductive isolation between sympatric populations. However, the genetic and developmental basis of hybrid incompatibility in higher vertebrates remains poorly understood. Mammals and birds, both amniotes, have similar developmental processes, but marked differences in development such as the XY/ZW sex determination systems and the presence or absence of genomic imprinting. Here, we review the sterile phenotype of hybrids between the Phodopus dwarf hamsters P. campbelli and P. sungorus, and the inviable phenotype of hybrids between two birds of the family Phasianidae, chicken (Gallus gallus domesticus) and Japanese quail (Coturnix japonica). We propose hypotheses for developmental defects that are associated with these hybrid incompatibilities. In addition, we discuss the genetic and developmental basis for these defects in conjunction with recent findings from mouse and avian models of genetics, reproductive biology and genomics. We suggest that these hybrids are ideal animal models for studying the genetic and developmental basis of hybrid incompatibility in amniotes.

Contrasting Levels of Molecular Evolution on the Mouse X Chromosome

The mammalian X chromosome has unusual evolutionary dynamics compared to autosomes. Faster-X evolution of spermatogenic protein-coding genes is known to be most pronounced for genes expressed late in spermatogenesis, but it is unclear if these patterns extend to other forms of molecular divergence. We tested for faster-X evolution in mice spanning three different forms of molecular evolution-divergence in protein sequence, gene expression, and DNA methylation-across different developmental stages of spermatogenesis. We used FACS to isolate individual cell populations and then generated cell-specific transcriptome profiles across different stages of spermatogenesis in two subspecies of house mice (Mus musculus), thereby overcoming a fundamental limitation of previous studies on whole tissues. We found faster-X protein evolution at all stages of spermatogenesis and faster-late protein evolution for both X-linked and autosomal genes. In contrast, there was less expression divergence late in spermatogenesis (slower late) on the X chromosome and for autosomal genes expressed primarily in testis (testis-biased). We argue that slower-late expression divergence reflects strong regulatory constraints imposed during this critical stage of sperm development and that these constraints are particularly acute on the tightly regulated sex chromosomes. We also found slower-X DNA methylation divergence based on genome-wide bisulfite sequencing of sperm from two species of mice (M. musculus and M. spretus), although it is unclear whether slower-X DNA methylation reflects development constraints in sperm or other X-linked phenomena. Our study clarifies key differences in patterns of regulatory and protein evolution across spermatogenesis that are likely to have important consequences for mammalian sex chromosome evolution, male fertility, and speciation.

Variation in promiscuity and sexual selection drives avian rate of Faster-Z evolution

Higher rates of coding sequence evolution have been observed on the Z chromosome relative to the autosomes across a wide range of species. However, despite a considerable body of theory, we lack empirical evidence explaining variation in the strength of the Faster-Z Effect. To assess the magnitude and drivers of Faster-Z Evolution, we assembled six de novo transcriptomes, spanning 90 million years of avian evolution. Our analysis combines expression, sequence and polymorphism data with measures of sperm competition and promiscuity. In doing so, we present the first empirical evidence demonstrating the positive relationship between Faster-Z Effect and measures of promiscuity, and therefore variance in male mating success. Our results from multiple lines of evidence indicate that selection is less effective on the Z chromosome, particularly in promiscuous species, and that Faster-Z Evolution in birds is due primarily to genetic drift. Our results reveal the power of mating system and sexual selection in shaping broad patterns in genome evolution.

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.

Purifying Selection Maintains Dosage-Sensitive Genes during Degeneration of the Threespine Stickleback Y Chromosome

Sex chromosomes are subject to unique evolutionary forces that cause suppression of recombination, leading to sequence degeneration and the formation of heteromorphic chromosome pairs (i.e., XY or ZW). Although progress has been made in characterizing the outcomes of these evolutionary processes on vertebrate sex chromosomes, it is still unclear how recombination suppression and sequence divergence typically occur and how gene dosage imbalances are resolved in the heterogametic sex. The threespine stickleback fish (Gasterosteus aculeatus) is a powerful model system to explore vertebrate sex chromosome evolution, as it possesses an XY sex chromosome pair at relatively early stages of differentiation. Using a combination of whole-genome and transcriptome sequencing, we characterized sequence evolution and gene expression across the sex chromosomes. We uncovered two distinct evolutionary strata that correspond with known structural rearrangements on the Y chromosome. In the oldest stratum, only a handful of genes remain, and these genes are under strong purifying selection. By comparing sex-linked gene expression with expression of autosomal orthologs in an outgroup, we show that dosage compensation has not evolved in threespine sticklebacks through upregulation of the X chromosome in males. Instead, in the oldest stratum, the genes that still possess a Y chromosome allele are enriched for genes predicted to be dosage sensitive in mammals and yeast. Our results suggest that dosage imbalances may have been avoided at haploinsufficient genes by retaining function of the Y chromosome allele through strong purifying selection.