Search

Abstract Background A detailed knowledge about spatial and temporal gene expression is important for understanding both the function of genes and their evolution. For the vast majority of species, transcriptomes are still largely uncharacterized and even in those where substantial information is available it is often in the form of partially sequenced transcriptomes. With the development of next generation sequencing, a single experiment can now simultaneously identify the transcribed part of a species genome and estimate levels of gene expression. Results mRNA from actively growing needles of Norway spruce (Picea abies) was sequenced using next generation sequencing technology. In total, close to 70 million fragments with a length of 76 bp were sequenced resulting in 5 Gbp of raw data. A de novo assembly of these reads, together with publicly available expressed sequence tag (EST) data from Norway spruce, was used to create a reference transcriptome. Of the 38,419 PUTs (putative unique transcripts) longer than 150 bp in this reference assembly, 83.5% show similarity to ESTs from other spruce species and of the remaining PUTs, 3,704 show similarity to protein sequences from other plant species, leaving 4,167 PUTs with limited similarity to currently available plant proteins. By predicting coding frames and comparing not only the Norway spruce PUTs, but also PUTs from the close relatives Picea glauca and Picea sitchensis to both Pinus taeda and Taxus mairei, we obtained estimates of synonymous and non-synonymous divergence among conifer species. In addition, we detected close to 15,000 SNPs of high quality and estimated gene expression differences between samples collected under dark and light conditions. Conclusions Our study yielded a large number of single nucleotide polymorphisms as well as estimates of gene expression on transcriptome scale. In agreement with a recent study we find that the synonymous substitution rate per year (0.6 × 10−09and 1.1 × 10−09) is an order of magnitude smaller than values reported for angiosperm herbs. However, if one takes generation time into account, most of this difference disappears. The estimates of the dN/dS ratio (non-synonymous over synonymous divergence) reported here are in general much lower than 1 and only a few genes showed a ratio larger than 1.

Transcriptional enhancers orchestrate cell type- and time point-specific gene expression programs. Genetic variation within enhancer sequences is an important contributor to phenotypic variation including evolutionary adaptations and human disease. Certain genes and pathways may be more prone to regulatory evolution than others, with different patterns across diverse organisms, but whether such patterns exist has not been investigated at a sufficient scale. To address this question, we identified signatures of accelerated sequence evolution in conserved enhancer elements throughout the mammalian phylogeny at an unprecedented scale. While different genes and pathways were enriched for regulatory evolution in different parts of the tree, we found a striking overall pattern of pleiotropic genes involved in gene regulatory and developmental processes being enriched for accelerated enhancer evolution. These genes were connected to more enhancers than other genes, which was the basis for having an increased amount of sequence acceleration over all their enhancers combined. We provide evidence that sequence acceleration is associated with turnover of regulatory function. Detailed study of one acceleration event in an enhancer of HES1 revealed that sequence evolution led to a new activity domain in the developing limb that emerged concurrently with the evolution of digit reduction in hoofed mammals. Our results provide evidence that enhancer evolution has been a frequent contributor to regulatory innovation at conserved developmental signaling genes in mammals., Competing Interests: Conflict of Interest The authors declare no conflicts of interest., (© The Author(s) 2024. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Human Accelerated Regions (HARs) are highly conserved across species but exhibit a significant excess of human-specific sequence changes, suggesting they may have gained novel functions in human evolution. HARs include transcriptional enhancers with human-specific activity and have been implicated in the evolution of the human brain. However, our understanding of how HARs contributed to uniquely human features of the brain is hindered by a lack of insight into the genes and pathways that HARs regulate. It is unclear whether HARs acted by altering the expression of gene targets conserved between HARs and their chimpanzee orthologs or by gaining new gene targets in human, a mechanism termed enhancer hijacking. We generated a high-resolution map of chromatin interactions for 1,590 HARs and their orthologs in human and chimpanzee neural stem cells (NSCs) to comprehensively identify gene targets in both species. HARs and their chimpanzee orthologs targeted a conserved set of 2,963 genes enriched for neurodevelopmental processes including neurogenesis and synaptic transmission. Changes in HAR enhancer activity were correlated with changes in conserved gene target expression. Conserved targets were enriched among genes differentially expressed between human and chimpanzee NSCs or between human and non-human primate developing and adult brain. Species-specific HAR gene targets did not converge on known biological functions and were not significantly enriched among differentially expressed genes, suggesting that HARs did not alter gene expression via enhancer hijacking. HAR gene targets, including differentially expressed targets, also showed cell type-specific expression patterns in the developing human brain, including outer radial glia, which are hypothesized to contribute to human cortical expansion. Our findings support that HARs influenced human brain evolution by altering the expression of conserved gene targets and provide the means to functionally link HARs with novel human brain features., Competing Interests: Competing Interests: The authors declare no competing financial interests.

Background: Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. CpG islands (CGIs) have recently been shown to influence enhancer activity, and here we test how their turnover across species contributes to enhancer evolution., Results: We integrate maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and find that CGI content in enhancers is strongly associated with increased histone modification levels. CGIs show widespread turnover across species and species-specific CGIs are strongly enriched for enhancers exhibiting species-specific activity across all tissues and species. Genes associated with enhancers with species-specific CGIs show concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon., Conclusions: Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals., (© 2024. The Author(s).)

Changes in gene regulation have been linked to the expansion of the human cerebral cortex and to neurodevelopmental disorders, potentially by altering neural progenitor proliferation. However, the effects of genetic variation within regulatory elements on neural progenitors remain obscure. We use sgRNA-Cas9 screens in human neural stem cells (hNSCs) to disrupt 10,674 genes and 26,385 conserved regions in 2,227 enhancers active in the developing human cortex and determine effects on proliferation. Genes with proliferation phenotypes are associated with neurodevelopmental disorders and show biased expression in specific fetal human brain neural progenitor populations. Although enhancer disruptions overall have weaker effects than gene disruptions, we identify enhancer disruptions that severely alter hNSC self-renewal. Disruptions in human accelerated regions, implicated in human brain evolution, also alter proliferation. Integrating proliferation phenotypes with chromatin interactions reveals regulatory relationships between enhancers and their target genes contributing to neurogenesis and potentially to human cortical evolution., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2024 The Authors. Published by Elsevier Inc. All rights reserved.)

Genetic variants in the third intron of the PRDM6 gene have been associated with BP traits in multiple GWAS. By combining fine mapping, massively parallel reporter assays, and gene editing, we identified super enhancers that drive the expression of PRDM6 and are partly regulated by STAT1 as the causal variants for hypertension. The heterozygous disruption of Prdm6 in mice expressing Cre recombinase under the control of mouse smooth muscle cell protein 22-α promoter (Prdm6fl/+ SM22-Cre) exhibited a markedly higher number of renin-producing cells in the kidneys at E18.5 compared with WT littermates and developed salt-induced systemic hypertension that was completely responsive to the renin inhibitor aliskiren. Strikingly, RNA-Seq analysis of the mouse aortas identified a network of PRDM6-regulated genes that are located in GWAS-associated loci for blood pressure, most notably Sox6, which modulates renin expression in the kidney. Accordingly, the smooth muscle cell-specific disruption of Sox6 in Prdm6fl/+ SM22-Cre mice resulted in a dramatic reduction of renin. Fate mapping and histological studies also showed increased numbers of neural crest-derived cells accompanied by increased collagen deposition in the kidneys of Prdm6fl/+ Wnt1Cre-ZsGreen1Cre mice compared with WT mice. These findings establish the role of PRDM6 as a regulator of renin-producing cell differentiation into smooth muscle cells and as an attractive target for the development of antihypertensive drugs.

The evolution of uniquely human traits likely entailed changes in developmental gene regulation. Human Accelerated Regions (HARs), which include transcriptional enhancers harboring a significant excess of human-specific sequence changes, are leading candidates for driving gene regulatory modifications in human development. However, insight into whether HARs alter the level, distribution, and timing of endogenous gene expression remains limited. We examined the role of the HAR HACNS1 (HAR2) in human evolution by interrogating its molecular functions in a genetically humanized mouse model. We find that HACNS1 maintains its human-specific enhancer activity in the mouse embryo and modifies expression of Gbx2, which encodes a transcription factor, during limb development. Using single-cell RNA-sequencing, we demonstrate that Gbx2 is upregulated in the limb chondrogenic mesenchyme of HACNS1 homozygous embryos, supporting that HACNS1 alters gene expression in cell types involved in skeletal patterning. Our findings illustrate that humanized mouse models provide mechanistic insight into how HARs modified gene expression in human evolution., (© 2022. The Author(s).)

The miRNA biogenesis is tightly regulated to avoid dysfunction and consequent disease development. Here, we describe modulation of miRNA processing as a novel noncanonical function of the 5-lipoxygenase (5-LO) enzyme in monocytic cells. In differentiated Mono Mac 6 (MM6) cells, we found an in situ interaction of 5-LO with Dicer, a key enzyme in miRNA biogenesis. RNA sequencing of small noncoding RNAs revealed a functional impact, knockout of 5-LO altered the expression profile of several miRNAs. Effects of 5-LO could be observed at two levels. qPCR analyses thus indicated that (a) 5-LO promotes the transcription of the evolutionarily conserved miR-99b/let-7e/miR-125a cluster and (b) the 5-LO-Dicer interaction downregulates the processing of pre-let-7e, resulting in an increase in miR-125a and miR-99b levels by 5-LO without concomitant changes in let-7e levels in differentiated MM6 cells. Our observations suggest that 5-LO regulates the miRNA profile by modulating the Dicer-mediated processing of distinct pre-miRNAs. 5-LO inhibits the formation of let-7e which is a well-known inducer of cell differentiation, but promotes the generation of miR-99b and miR-125a known to induce cell proliferation and the maintenance of leukemic stem cell functions., (© 2020 The Authors. The FASEB Journal published by Wiley Periodicals LLC on behalf of Federation of American Societies for Experimental Biology.)

Genetic changes that altered the function of gene regulatory elements have been implicated in the evolution of human traits such as the expansion of the cerebral cortex. However, identifying the particular changes that modified regulatory activity during human evolution remain challenging. Here we used massively parallel enhancer assays in neural stem cells to quantify the functional impact of >32,000 human-specific substitutions in >4,300 human accelerated regions (HARs) and human gain enhancers (HGEs), which include enhancers with novel activities in humans. We found that >30% of active HARs and HGEs exhibited differential activity between human and chimpanzee. We isolated the effects of human-specific substitutions from background genetic variation to identify the effects of genetic changes most relevant to human evolution. We found that substitutions interacted in both additive and nonadditive ways to modify enhancer function. Substitutions within HARs, which are highly constrained compared to HGEs, showed smaller effects on enhancer activity, suggesting that the impact of human-specific substitutions is buffered in enhancers with constrained ancestral functions. Our findings yield insight into how human-specific genetic changes altered enhancer function and provide a rich set of candidates for studies of regulatory evolution in humans., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

Dicer is a ribonuclease III enzyme in biosynthesis of micro-RNAs (miRNAs). Here we describe a regulation of Dicer expression in monocytic cells, based on proteolysis. In undifferentiated Mono Mac 6 (MM6) cells, full-length Dicer was undetectable; only an ∼50-kDa fragment appeared in Western blots. However, when MM6 cells were treated with zymosan or LPS during differentiation with TGF-β and 1,25diOHvitD3, full-length Dicer became abundant together with varying amounts of ∼170- and ∼50-kDa Dicer fragments. Mass spectrometry identified the Dicer fragments and showed cleavage about 450 residues upstream from the C terminus. Also, PGE 2 (prostaglandin E2) added to differentiating MM6 cells up-regulated full-length Dicer, through EP2/EP4 and cAMP. The TLR stimuli strongly induced miR-146a-5p, while PGE 2 increased miR-99a-5p and miR-125a-5p, both implicated in down-regulation of TNFα. The Ser protease inhibitor AEBSF (4-[2-aminoethyl] benzene sulfonyl fluoride) up-regulated full-length Dicer, both in MM6 cells and in primary human blood monocytes, indicating a specific proteolytic degradation. However, AEBSF alone did not lead to a general increase in miR expression, indicating that additional mechanisms are required to increase miRNA biosynthesis. Finally, differentiation of monocytes to macrophages with M-CSF or GM-CSF strongly up-regulated full-length Dicer. Our results suggest that differentiation regimens, both in the MM6 cell line and of peripheral blood monocytes, inhibit an apparently constitutive Dicer proteolysis, allowing for increased formation of miRNAs., Competing Interests: The authors declare no competing interest., (Copyright © 2020 the Author(s). Published by PNAS.)

Sympatric speciation remains a central focus of evolutionary biology. Although some evidence shows speciation occurring in this way, little is known about the gene expression evolution and the characteristics of population genetics as species diverge. Two closely related Gymnocypris fish (Gymnocypris chui and Gymnocypris scleracanthus), which come from a small glacier lake in the Tibetan Plateau, Lake Langcuo, exist a possible incipient sympatric adaptive ecological speciation. We generated large amounts of RNA-Seq data from multiple individuals and tissues from each of the two species and compared gene expression patterns and genetic polymorphisms between them. Ordination analysis separated samples by organ rather than by species. The degree of expression difference between organs within and between species was different. Phylogenetic analyses indicated that the two closely related taxa formed a monophyletic complex. Population structure analysis displayed two distinctly divergent clusters of G. chui and G. scleracanthus populations. By contrast, G. scleracanthus population genetic diversity is higher than that of G. chui. Considerable sites of the two populations were differentiated with a coefficient of FST = 0.25-0.50, implying that a small proportion of loci nevertheless exhibited deep divergence in two comparisons. Concomitantly, putatively selected genes during speciation revealed functional categories are enriched in bone morphogenesis, cell growth, neurogenetics, enzyme activity, and binding activity in G. chui population. In contrast, nutrition and localization were highlighted in G. scleracanthus. Collectively, morphological traits and dietary preference combine with genetic variation and expression variation, probably contributed to the incipient speciation of two sympatric populations., (© The Author(s) 2020. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Variation in gene expression is believed to make a significant contribution to phenotypic diversity and divergence. The analysis of allele-specific expression (ASE) can reveal important insights into gene expression regulation. We developed a novel method called RPASE (Read-backed Phasing-based ASE detection) to test for genes that show ASE. With mapped RNA-seq data from a single individual and a list of SNPs from the same individual as the only input, RPASE is capable of aggregating information across multiple dependent SNPs and producing individual-based gene-level tests for ASE. RPASE performs well in simulations and comparisons. We applied RPASE to multiple bird species and found a potentially rich landscape of ASE., (© 2018 John Wiley & Sons Ltd.)

Developmental gene expression patterns are orchestrated by thousands of distant-acting transcriptional enhancers. However, identifying enhancers essential for the expression of their target genes has proven challenging. Maps of long-range regulatory interactions may provide the means to identify enhancers crucial for developmental gene expression. To investigate this hypothesis, we used circular chromosome conformation capture coupled with interaction maps in the mouse limb to characterize the regulatory topology of Pitx1 , which is essential for hindlimb development. We identified a robust hindlimb-specific interaction between Pitx1 and a putative hindlimb-specific enhancer. To interrogate the role of this interaction in Pitx1 regulation, we used genome editing to delete this enhancer in mouse. Although deletion of the enhancer completely disrupts the interaction, Pitx1 expression in the hindlimb is only mildly affected, without any detectable compensatory interactions between the Pitx1 promoter and potentially redundant enhancers. Pitx1 enhancer null mice did not exhibit any of the characteristic morphological defects of the Pitx1 -/- mutant. Our results suggest that robust, tissue-specific physical interactions at essential developmental genes have limited predictive value for identifying enhancer mutations with strong loss-of-function phenotypes., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2018. Published by The Company of Biologists Ltd.)

Polymorphism in cis-regulatory sequences can lead to different levels of expression for the two alleles of a gene, providing a starting point for the evolution of gene expression. Little is known about the genome-wide abundance of genetic variation in gene regulation in natural populations but analysis of allele-specific expression (ASE) provides a means for investigating such variation. We performed RNA-seq of multiple tissues from population samples of two closely related flycatcher species and developed a Bayesian algorithm that maximizes data usage by borrowing information from the whole data set and combines several SNPs per transcript to detect ASE. Of 2,576 transcripts analyzed in collared flycatcher, ASE was detected in 185 (7.2%) and a similar frequency was seen in the pied flycatcher. Transcripts with statistically significant ASE commonly showed the major allele in >90% of the reads, reflecting that power was highest when expression was heavily biased toward one of the alleles. This would suggest that the observed frequencies of ASE likely are underestimates. The proportion of ASE transcripts varied among tissues, being lowest in testis and highest in muscle. Individuals often showed ASE of particular transcripts in more than one tissue (73.4%), consistent with a genetic basis for regulation of gene expression. The results suggest that genetic variation in regulatory sequences commonly affects gene expression in natural populations and that it provides a seedbed for phenotypic evolution via divergence in gene expression., (© The Author 2017. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Relatively little is known about the character of gene expression evolution as species diverge. It is for instance unclear if gene expression generally evolves in a clock-like manner (by stabilizing selection or neutral evolution) or if there are frequent episodes of directional selection. To gain insights into the evolutionary divergence of gene expression, we sequenced and compared the transcriptomes of multiple organs from population samples of collared (Ficedula albicollis) and pied flycatchers (F. hypoleuca), two species which diverged less than one million years ago. Ordination analysis separated samples by organ rather than by species. Organs differed in their degrees of expression variance within species and expression divergence between species. Variance was negatively correlated with expression breadth and protein interactivity, suggesting that pleiotropic constraints reduce gene expression variance within species. Variance was correlated with between-species divergence, consistent with a pattern expected from stabilizing selection and neutral evolution. Using an expression PST approach, we identified genes differentially expressed between species and found 16 genes uniquely expressed in one of the species. For one of these, DPP7, uniquely expressed in collared flycatcher, the absence of expression in pied flycatcher could be associated with a ≈20-kb deletion including 11 of 13 exons. This study of a young vertebrate speciation model system expands our knowledge of how gene expression evolves as natural populations become reproductively isolated., (© 2016 The Authors. Molecular Ecology Published by John Wiley & Sons Ltd.)

There is increasing evidence that dosage compensation is not a ubiquitous feature following sex chromosome evolution, especially not in organisms where females are the heterogametic sex, like in birds. Even when it occurs, compensation can be incomplete and limited to dosage-sensitive genes. However, previous work has mainly studied transcriptional regulation of sex-linked genes, which may not reflect expression at the protein level. Here, we used liquid chromatography-tandem mass spectrometry to detect and quantify expressed levels of more than 2,400 proteins in ten different tissues of male and female chicken embryos. For comparison, transcriptome sequencing was performed in the same individuals, five of each sex. The proteomic analysis revealed that dosage compensation was incomplete, with a mean male-to-female (M:F) expression ratio of Z-linked genes of 1.32 across tissues, similar to that at the RNA level (1.29). The mean Z chromosome-to-autosome expression ratio was close to 1 in males and lower than 1 in females, consistent with partly reduced Z chromosome expression in females. Although our results exclude a general mechanism for chromosome-wide dosage compensation at translation, 30% of all proteins encoded from Z-linked genes showed a significant change in the M:F ratio compared with the corresponding ratio at the RNA level. This resulted in a pattern where some genes showed balanced expression between sexes and some close to 2-fold higher expression in males. This suggests that proteomic analyses will be necessary to reveal a more complete picture of gene regulation and sex chromosome evolution., (© The Author 2015. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The typically repetitive nature of the sex-limited chromosome means that it is often excluded from or poorly covered in genome assemblies, hindering studies of evolutionary and population genomic processes in non-recombining chromosomes. Here, we present a draft assembly of the non-recombining region of the collared flycatcher W chromosome, containing 46 genes without evidence of female-specific functional differentiation. Survival of genes during W chromosome degeneration has been highly non-random and expression data suggest that this can be attributed to selection for maintaining gene dose and ancestral expression levels of essential genes. Re-sequencing of large population samples revealed dramatically reduced levels of within-species diversity and elevated rates of between-species differentiation (lineage sorting), consistent with low effective population size. Concordance between W chromosome and mitochondrial DNA phylogenetic trees demonstrates evolutionary stable matrilineal inheritance of this nuclear-cytonuclear pair of chromosomes. Our results show both commonalities and differences between W chromosome and Y chromosome evolution.

The molecular characteristics of the pseudoautosomal region (PAR) of sex chromosomes remain elusive. Despite significant genome-sequencing efforts, the PAR of highly differentiated avian sex chromosomes remains to be identified. Here we use linkage analysis together with whole-genome re-sequencing to uncover the 630-kb PAR of an ecological model species, the collared flycatcher. The PAR contains 22 protein-coding genes and is GC rich. The genetic length is 64 cM in female meiosis, consistent with an obligate crossing-over event. Recombination is concentrated to a hotspot region, with an extreme rate of >700 cM/Mb in a 67-kb segment. We find no signatures of sexual antagonism and propose that sexual antagonism may have limited influence on PAR sequences when sex chromosomes are nearly fully differentiated and when a recombination hotspot region is located close to the PAR boundary. Our results demonstrate that a very small PAR suffices to ensure homologous recombination and proper segregation of sex chromosomes during meiosis.

Sex chromosome divergence, which follows the cessation of recombination and degeneration of the sex-limited chromosome, can cause a reduction in expression level for sex-linked genes in the heterozygous sex, unless some mechanisms of dosage compensation develops to counter the reduction in gene dose. Because large-scale perturbations in expression levels arising from changes in gene dose might have strong deleterious effects, the evolutionary response should be strong. However, in birds and in at least some other female heterogametic organisms, wholesale sex chromosome dosage compensation does not seem to occur. Using RNA-seq of multiple tissues and individuals, we investigated male and female expression levels of Z-linked and autosomal genes in the collared flycatcher, a bird for which a draft genome sequence recently has been reported. We found that male expression of Z-linked genes was on average 50% higher than female expression, although there was considerable variation in the male-to-female ratio among genes. The ratio for individual genes was well correlated among tissues and there was also a correlation in the extent of compensation between flycatcher and chicken orthologs. The relative excess of male expression was positively correlated with expression breadth, expression level, and number of interacting proteins (protein connectivity), and negatively correlated with variance in expression. These observations lead to a model of compensation occurring on a gene-by-gene basis, supported by an absence of clustering of genes on the Z chromosome with respect to the extent of compensation. Equal mean expression level of autosomal and Z-linked genes in males, and 50% higher expression of autosomal than Z-linked genes in females, is compatible with that partial compensation is achieved by hypertranscription from females' single Z chromosome. A comparison with male-to-female expression ratios in orthologous Z-linked genes of ostriches, where Z-W recombination still occurs, suggests that male-biased expression of Z-linked genes is a derived trait after avian sex chromosome divergence.

Unravelling the genomic landscape of divergence between lineages is key to understanding speciation. The naturally hybridizing collared flycatcher and pied flycatcher are important avian speciation models that show pre- as well as postzygotic isolation. We sequenced and assembled the 1.1-Gb flycatcher genome, physically mapped the assembly to chromosomes using a low-density linkage map and re-sequenced population samples of each species. Here we show that the genomic landscape of species differentiation is highly heterogeneous with approximately 50 'divergence islands' showing up to 50-fold higher sequence divergence than the genomic background. These non-randomly distributed islands, with between one and three regions of elevated divergence per chromosome irrespective of chromosome size, are characterized by reduced levels of nucleotide diversity, skewed allele-frequency spectra, elevated levels of linkage disequilibrium and reduced proportions of shared polymorphisms in both species, indicative of parallel episodes of selection. Proximity of divergence peaks to genomic regions resistant to sequence assembly, potentially including centromeres and telomeres, indicate that complex repeat structures may drive species divergence. A much higher background level of species divergence of the Z chromosome, and a lower proportion of shared polymorphisms, indicate that sex chromosomes and autosomes are at different stages of speciation. This study provides a roadmap to the emerging field of speciation genomics.