HALDANE'S RULE IN AN AVIAN SYSTEM: USING CLINE THEORY AND DIVERGENCE POPULATION GENETICS TO TEST FOR DIFFERENTIAL INTROGRESSION OF MITOCHONDRIAL, AUTOSOMAL, AND SEX-LINKED LOCI ACROSS THEPASSERINABUNTING HYBRID ZONE

Ongoing introgression of a secondary sexual plumage trait in a stable avian hybrid zone

Hybrid zones are dynamic systems where natural selection, sexual selection, and other evolutionary forces can act on reshuffled combinations of distinct genomes. The movement of hybrid zones, individual traits, or both are of particular interest for understanding the interplay between selective processes. In a hybrid zone involving two lek-breeding birds, secondary sexual plumage traits of Manacus vitellinus, including bright yellow collar and olive belly color, have introgressed ~50 km asymmetrically across the genomic center of the zone into populations more genetically similar to Manacus candei. Males with yellow collars are preferred by females and are more aggressive than parental M. candei, suggesting that sexual selection was responsible for the introgression of male traits. We assessed the spatial and temporal dynamics of this hybrid zone using historical (1989-1994) and contemporary (2017-2020) transect samples to survey both morphological and genetic variation. Genome-wide single nucleotide polymorphism data and several male phenotypic traits show that the genomic center of the zone has remained spatially stable, whereas the olive belly color of male M. vitellinus has continued to introgress over this time period. Our data suggest that sexual selection can continue to shape phenotypes dynamically, independent of a stable genomic transition between species.

Quantitative Analyses of Coupling in Hybrid Zones

In hybrid zones, whether barrier loci experience selection mostly independently or as a unit depends on the ratio of selection to recombination as captured by the coupling coefficient. Theory predicts a sharper transition between an uncoupled and coupled system when more loci affect hybrid fitness. However, the extent of coupling in hybrid zones has rarely been quantified. Here, we use simulations to characterize the relationship between the coupling coefficient and variance in clines across genetic loci. We then reanalyze 25 hybrid zone data sets and find that cline variances and estimated coupling coefficients form a smooth continuum from high variance and weak coupling to low variance and strong coupling. Our results are consistent with low rates of hybridization and a strong genome-wide barrier to gene flow when the coupling coefficient is much greater than 1, but also suggest that this boundary might be approached gradually and at a near constant rate over time.

The population bottleneck of the Iberian wolf impacted genetic diversity but not admixture with domestic dogs: A temporal genomic approach

After decades of intense persecution, the Iberian wolf subspecies faced a severe bottleneck in the 1970s that considerably reduced its range and population size, nearly leading to its extinction in central and southern Iberian Peninsula. Such population decline could have impacted the genetic diversity of Iberian wolves through different processes, namely genetic drift and dynamics of hybridization with domestic dogs. By contrasting the genomes of 68 contemporary with 54 historical samples spanning the periods before and immediately after the 1970s bottleneck, we found evidence of its impact on genetic diversity and dynamics of wolf-dog hybridization. Our genome-wide assessment revealed that wolves and dogs form two well-differentiated genetic groups in Iberia and that hybridization rates did not increase during the bottleneck. However, an increased number of hybrid individuals was found over time during the population re-expansion, particularly at the edge of the wolf range. We estimated a low percentage of dog ancestry (~1.4%) in historical samples, suggesting that dog introgression was not a key driver for wolf extinction in central and southern Iberia. Our findings also unveil a significant decline in genetic diversity in contemporary samples, with the highest proportion of homozygous segments in the genome being recently inherited. Overall, our study provides unprecedented insight into the impact of a sharp decline on the Iberian wolf genome and refines our understanding of the ecological and evolutionary drivers of wolf-dog hybridization in the wild.

Hybridization in late stages of speciation: Strong but incomplete genome-wide reproductive isolation and 'large Z-effect' in a moving hybrid zone

In organisms reproducing sexually, speciation occurs when increasing divergence results in pre- or post-zygotic reproductive isolation between lineages. Studies focusing on reproductive isolation origin in early stages of speciation are common and many rely on genomic scans to infer introgression providing limited information on the genomic architecture of reproductive isolation long-term maintenance. This study analyses a natural hybrid zone between two species in a late stage of speciation. We used ddRADseq genotyping in the contact between Podarcis bocagei and P. carbonelli to examine admixture extent, analyse hybrid zone stability and assess genome-wide variation in selection against introgression. We confirmed strong but incomplete reproductive isolation in a bimodal hybrid zone. New findings revealed population genetic structure within P. carbonelli in the contact zone; geographical and genomic clines analysis suggested strong selection against gene flow, but a relatively small proportion of the loci can introgress, mostly within the narrow contact zone. However, geographical clines revealed that a few introgressed loci show signs of potential positive selection, particularly into P. bocagei. Geographical clines also detected a signal of hybrid zone movement towards P. bocagei distribution. Genomic cline analysis revealed heterogeneous patterns of introgression among loci within the syntopy zone, but the majority maintain a strong association with the genomic background of origin. However, incongruences between both cline approaches were found, potentially driven by confounding effects on genomic clines. Last, an important role of the Z chromosome in reproductive isolation is suggested. Importantly, overall patterns of restricted introgression seem to result from numerous strong intrinsic barriers across the genome.

Temporal stability of the hybrid zone between Calocitta magpie-jays revealed through comparison of museum specimens and iNaturalist photos

Hybrid zones are natural experiments for the study of avian evolution. Hybrid zones can be dynamic, moving as species adjust to new climates and habitats, with unknown implications for species and speciation. There are relatively few studies that have comparable modern and historic sampling to assess change in hybrid zone location and width over time, and those studies have generally found mixed results, with many hybrid zones showing change over time, but others showing stability. The white-throated magpie-jay (Calocitta formosa) and black-throated magpie-jay (Calocitta colliei) occur along the western coast of Mexico and Central America. The two species differ markedly in throat color and tail length, and prior observation suggests a narrow hybrid zone in southern Jalisco where individuals have mixed throat color. This study aims to assess the existence and temporal stability of this putative hybrid zone by comparing throat color between georeferenced historical museum specimens and modern photos from iNaturalist with precise locality information. Our results confirm the existence of a narrow hybrid zone in Jalisco, with modern throat scores gradually increasing from the parental ends of the cline toward the cline center in a sigmoidal curve characteristic of hybrid zones. Our temporal comparison suggests that the hybrid zone has not shifted its position between historical (pre-1973) and modern (post-2005) time periods-a surprising result given the grand scale of habitat change to the western Mexican lowlands during this time. An anomalous pocket of white-throated individuals in the northern range of the black-throated magpie-jay hints at the possibility of prehistorical long-distance introduction. Future genomic data will help disentangle the evolutionary history of these lineages and better characterize how secondary contact is affecting both the DNA and the phenotype of these species.

Population genomic evidence of selection on structural variants in a natural hybrid zone

Structural variants (SVs) can promote speciation by directly causing reproductive isolation or by suppressing recombination across large genomic regions. Whereas examples of each mechanism have been documented, systematic tests of the role of SVs in speciation are lacking. Here, we take advantage of long-read (Oxford nanopore) whole-genome sequencing and a hybrid zone between two Lycaeides butterfly taxa (L. melissa and Jackson Hole Lycaeides) to comprehensively evaluate genome-wide patterns of introgression for SVs and relate these patterns to hypotheses about speciation. We found >100,000 SVs segregating within or between the two hybridizing species. SVs and SNPs exhibited similar levels of genetic differentiation between species, with the exception of inversions, which were more differentiated. We detected credible variation in patterns of introgression among SV loci in the hybrid zone, with 562 of 1419 ancestry-informative SVs exhibiting genomic clines that deviated from null expectations based on genome-average ancestry. Overall, hybrids exhibited a directional shift towards Jackson Hole Lycaeides ancestry at SV loci, consistent with the hypothesis that these loci experienced more selection on average than SNP loci. Surprisingly, we found that deletions, rather than inversions, showed the highest skew towards excess ancestry from Jackson Hole Lycaeides. Excess Jackson Hole Lycaeides ancestry in hybrids was also especially pronounced for Z-linked SVs and inversions containing many genes. In conclusion, our results show that SVs are ubiquitous and suggest that SVs in general, but especially deletions, might disproportionately affect hybrid fitness and thus contribute to reproductive isolation.

Limited movement of an avian hybrid zone in relation to regional variation in magnitude of climate change

Studies of natural hybrid zones can provide documentation of range shifts in response to climate change and identify loci important to reproductive isolation. Using a temporal (36-38 years) comparison of the black-capped (Poecile atricapillus) and Carolina (P. carolinensis) chickadee hybrid zone, we investigated movement of the western portion of the zone (western Missouri) and assessed whether loci and pathways underpinning reproductive isolation were similar to those in the eastern portion of the hybrid zone. Using 92 birds sampled along the hybrid zone transect in 2016 and 68 birds sampled between 1978 and 1980, we generated 11,669 SNPs via ddRADseq. These SNPs were used to assess movement of the hybrid zone through time and to evaluate variation in introgression among loci. We demonstrate that the interface has moved ~5 km to the northwest over the last 36-38 years, that is, at only one-fifth the rate at which the eastern portion (e.g., Pennsylvania, Ohio) of the hybrid zone has moved. Temperature trends over the last 38 years reveal that eastern areas have warmed 50% more than western areas in terms of annual mean temperature, possibly providing an explanation for the slower movement of the hybrid zone in Missouri. Our results suggest hybrid zone movement in broadly distributed species, such as chickadees, will vary between areas in response to local differences in the impacts of climate change.

Genotype-environment associations across spatial scales reveal the importance of putative adaptive genetic variation in divergence

Identifying areas of high evolutionary potential is a judicious strategy for developing conservation priorities in the face of environmental change. For wide-ranging species occupying heterogeneous environments, the evolutionary forces that shape distinct populations can vary spatially. Here, we investigate patterns of genomic variation and genotype-environment associations in the hermit thrush (Catharus guttatus), a North American songbird, at broad (across the breeding range) and narrow spatial scales (at a hybrid zone). We begin by building a genoscape or map of genetic variation across the breeding range and find five distinct genetic clusters within the species, with the greatest variation occurring in the western portion of the range. Genotype-environment association analyses indicate higher allelic turnover in the west than in the east, with measures of temperature surfacing as key predictors of putative adaptive genomic variation rangewide. Since broad patterns detected across a species' range represent the aggregate of many locally adapted populations, we investigate whether our broadscale analysis is consistent with a finer scale analysis. We find that top rangewide temperature-associated loci vary in their clinal patterns (e.g., steep clines vs. fixed allele frequencies) across a hybrid zone in British Columbia, suggesting that the environmental predictors and the associated candidate loci identified in the rangewide analysis are of variable importance in this particular region. However, two candidate loci exhibit strong concordance with the temperature gradient in British Columbia, suggesting a potential role for temperature-related barriers to gene flow and/or temperature-driven ecological selection in maintaining putative local adaptation. This study demonstrates how patterns identified at the broad (macrogeographic) scale can be validated by investigating genotype-environment correlations at the local (microgeographic) scale. Furthermore, our results highlight the importance of considering the spatial distribution of putative adaptive variation when assessing population-level sensitivity to climate change and other stressors.

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.

Sex chromosomes as supergenes of speciation: why amphibians defy the rules?

As reflected by the two rules of speciation (Haldane's rule and the large X-/Z-effect), sex chromosomes are expected to behave like supergenes of speciation: they recombine only in one sex (XX females or ZZ males), supposedly recruit sexually antagonistic genes and evolve faster than autosomes, which can all contribute to pre-zygotic and post-zygotic isolation. While this has been mainly studied in organisms with conserved sex-determining systems and highly differentiated (heteromorphic) sex chromosomes like mammals, birds and some insects, these expectations are less clear in organismal groups where sex chromosomes repeatedly change and remain mostly homomorphic, like amphibians. In this article, we review the proposed roles of sex-linked genes in isolating nascent lineages throughout the speciation continuum and discuss their support in amphibians given current knowledge of sex chromosome evolution and speciation modes. Given their frequent recombination and lack of differentiation, we argue that amphibian sex chromosomes are not expected to become supergenes of speciation, which is reflected by the rarity of empirical studies consistent with a 'large sex chromosome effect' in frogs and toads. The diversity of sex chromosome systems in amphibians has a high potential to disentangle the evolutionary mechanisms responsible for the emergence of sex-linked speciation genes in other organisms. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.

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.

Displaced clines in an avian hybrid zone (Thamnophilidae: Rhegmatorhina) within an Amazonian interfluve

Secondary contact between species often results in the formation of a hybrid zone, with the eventual fates of the hybridizing species dependent on evolutionary and ecological forces. We examine this process in the Amazon Basin by conducting the first genomic and phenotypic characterization of the hybrid zone formed after secondary contact between two obligate army-ant-followers: the White-breasted Antbird (Rhegmatorhina hoffmannsi) and the Harlequin Antbird (Rhegmatorhina berlepschi). We found a major geographic displacement (∼120 km) between the mitochondrial and nuclear clines, and we explore potential hypotheses for the displacement, including sampling error, genetic drift, and asymmetric cytonuclear incompatibilities. We cannot exclude roles for sampling error and genetic drift in contributing to the discordance; however, the data suggest expansion and unidirectional introgression of hoffmannsi into the distribution of berlepschi.

Sex chromosomes in meiotic, hemiclonal, clonal and polyploid hybrid vertebrates: along the 'extended speciation continuum'

We review knowledge about the roles of sex chromosomes in vertebrate hybridization and speciation, exploring a gradient of divergences with increasing reproductive isolation (speciation continuum). Under early divergence, well-differentiated sex chromosomes in meiotic hybrids may cause Haldane-effects and introgress less easily than autosomes. Undifferentiated sex chromosomes are more susceptible to introgression and form multiple (or new) sex chromosome systems with hardly predictable dominance hierarchies. Under increased divergence, most vertebrates reach complete intrinsic reproductive isolation. Slightly earlier, some hybrids (linked in 'the extended speciation continuum') exhibit aberrant gametogenesis, leading towards female clonality. This facilitates the evolution of various allodiploid and allopolyploid clonal ('asexual') hybrid vertebrates, where 'asexuality' might be a form of intrinsic reproductive isolation. A comprehensive list of 'asexual' hybrid vertebrates shows that they all evolved from parents with divergences that were greater than at the intraspecific level (K2P-distances of greater than 5-22% based on mtDNA). These 'asexual' taxa inherited genetic sex determination by mostly undifferentiated sex chromosomes. Among the few known sex-determining systems in hybrid 'asexuals', female heterogamety (ZW) occurred about twice as often as male heterogamety (XY). We hypothesize that pre-/meiotic aberrations in all-female ZW-hybrids present Haldane-effects promoting their evolution. Understanding the preconditions to produce various clonal or meiotic allopolyploids appears crucial for insights into the evolution of sex, 'asexuality' and polyploidy. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Mass of genes rather than master genes underlie the genomic architecture of amphibian speciation

The genetic architecture of speciation, i.e., how intrinsic genomic incompatibilities promote reproductive isolation (RI) between diverging lineages, is one of the best-kept secrets of evolution. To directly assess whether incompatibilities arise in a limited set of large-effect speciation genes, or in a multitude of loci, we examined the geographic and genomic landscapes of introgression across the hybrid zones of 41 pairs of frog and toad lineages in the Western Palearctic region. As the divergence between lineages increases, phylogeographic transitions progressively become narrower, and larger parts of the genome resist introgression. This suggests that anuran speciation proceeds through a gradual accumulation of multiple barrier loci scattered across the genome, which ultimately deplete hybrid fitness by intrinsic postzygotic isolation, with behavioral isolation being achieved only at later stages. Moreover, these loci were disproportionately sex linked in one group (Hyla) but not in others (Rana and Bufotes), implying that large X-effects are not necessarily a rule of speciation with undifferentiated sex chromosomes. The highly polygenic nature of RI and the lack of hemizygous X/Z chromosomes could explain why the speciation clock ticks slower in amphibians compared to other vertebrates. The clock-like dynamics of speciation combined with the analytical focus on hybrid zones offer perspectives for more standardized practices of species delimitation.

X-Linked Signature of Reproductive Isolation in Humans is Mirrored in a Howler Monkey Hybrid Zone

Reproductive isolation is a fundamental step in speciation. While sex chromosomes have been linked to reproductive isolation in many model systems, including hominids, genetic studies of the contribution of sex chromosome loci to speciation for natural populations are relatively sparse. Natural hybrid zones can help identify genomic regions contributing to reproductive isolation, like hybrid incompatibility loci, since these regions exhibit reduced introgression between parental species. Here, we use a primate hybrid zone (Alouatta palliata × Alouatta pigra) to test for reduced introgression of X-linked SNPs compared to autosomal SNPs. To identify X-linked sequence in A. palliata, we used a sex-biased mapping approach with whole-genome re-sequencing data. We then used genomic cline analysis with reduced-representation sequence data for parental A. palliata and A. pigra individuals and hybrids (n = 88) to identify regions with non-neutral introgression. We identified ~26 Mb of non-repetitive, putatively X-linked genomic sequence in A. palliata, most of which mapped collinearly to the marmoset and human X chromosomes. We found that X-linked SNPs had reduced introgression and an excess of ancestry from A. palliata as compared to autosomal SNPs. One outlier region with reduced introgression overlaps a previously described "desert" of archaic hominin ancestry on the human X chromosome. These results are consistent with a large role for the X chromosome in speciation across animal taxa and further, suggest shared features in the genomic basis of the evolution of reproductive isolation in primates.

Evaluating evidence of mitonuclear incompatibilities with the sex chromosomes in an avian hybrid zone

The exploration of hybrid zones and the intergenomic conflicts exposed through hybridization provide windows into the processes of divergence and speciation. Sex chromosomes and mitonuclear incompatibilities have strong associations with the genetics of hybrid dysfunction. In ZW sex-determining systems, maternal co-inheritance of the mitochondrial and W chromosomes immediately exposes incompatibilities between these maternal contributions of one species and the Z chromosome of another. We analyze mitochondrial and Z chromosome admixture in the long-tailed finch (Poephila acuticauda) of Australia, where hybridizing subspecies differ prominently in Z chromosome genotype and in bill color, yet the respective centers of geographic admixture for these two traits are offset by 350 km. We report two well-defined mitochondrial clades that diverged ∼0.5 million years ago. Mitochondrial contact is geographically co-located within a hybrid zone of Z chromosome admixture and is displaced from bill color admixture by nearly 400 km. Consistent with Haldane's rule expectations, hybrid zone females are significantly less likely than males to carry an admixed Z chromosome or have mismatched Z-mitochondrial genotypes. Furthermore, there are significantly fewer than expected mitonuclear mismatches in hybrid zone females and paternal backcross males. Results suggest a potential for mitonuclear/sex chromosome incompatibilities in the emergence of reproductive isolation in this system.

New view on the organization and evolution of Palaeognathae mitogenomes poses the question on the ancestral gene rearrangement in Aves

BACKGROUND

Bird mitogenomes differ from other vertebrates in gene rearrangement. The most common avian gene order, identified first in Gallus gallus, is considered ancestral for all Aves. However, other rearrangements including a duplicated control region and neighboring genes have been reported in many representatives of avian orders. The repeated regions can be easily overlooked due to inappropriate DNA amplification or genome sequencing. This raises a question about the actual prevalence of mitogenomic duplications and the validity of the current view on the avian mitogenome evolution. In this context, Palaeognathae is especially interesting because is sister to all other living birds, i.e. Neognathae. So far, a unique duplicated region has been found in one palaeognath mitogenome, that of Eudromia elegans.

RESULTS

Therefore, we applied an appropriate PCR strategy to look for omitted duplications in other palaeognaths. The analyses revealed the duplicated control regions with adjacent genes in Crypturellus, Rhea and Struthio as well as ND6 pseudogene in three moas. The copies are very similar and were subjected to concerted evolution. Mapping the presence and absence of duplication onto the Palaeognathae phylogeny indicates that the duplication was an ancestral state for this avian group. This feature was inherited by early diverged lineages and lost two times in others. Comparison of incongruent phylogenetic trees based on mitochondrial and nuclear sequences showed that two variants of mitogenomes could exist in the evolution of palaeognaths. Data collected for other avian mitogenomes revealed that the last common ancestor of all birds and early diverging lineages of Neoaves could also possess the mitogenomic duplication.

CONCLUSIONS

The duplicated control regions with adjacent genes are more common in avian mitochondrial genomes than it was previously thought. These two regions could increase effectiveness of replication and transcription as well as the number of replicating mitogenomes per organelle. In consequence, energy production by mitochondria may be also more efficient. However, further physiological and molecular analyses are necessary to assess the potential selective advantages of the mitogenome duplications.

Avian Diversity: Speciation, Macroevolution, and Ecological Function

The origin, distribution, and function of biological diversity are fundamental themes of ecology and evolutionary biology. Research on birds has played a major role in the history and development of these ideas, yet progress was for many decades limited by a focus on patterns of current diversity, often restricted to particular clades or regions. Deeper insight is now emerging from a recent wave of integrative studies combining comprehensive phylogenetic, environmental, and functional trait data at unprecedented scales. We review these empirical advances and describe how they are reshaping our understanding of global patterns of bird diversity and the processes by which it arises, with implications for avian biogeography and functional ecology. Further expansion and integration of data sets may help to resolve longstanding debates about the evolutionary origins of biodiversity and offer a framework for understanding and predicting the response of ecosystems to environmental change.

Dispersal Predicts Hybrid Zone Widths across Animal Diversity: Implications for Species Borders under Incomplete Reproductive Isolation

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

Differential divergence in autosomes and sex chromosomes is associated with intra-island diversification at a very small spatial scale in a songbird lineage

Recently diverged taxa showing marked phenotypic and ecological diversity provide optimal systems to understand the genetic processes underlying speciation. We used genome-wide markers to investigate the diversification of the Reunion grey white-eye (Zosterops borbonicus) on the small volcanic island of Reunion (Mascarene archipelago), where this species complex exhibits four geographical forms that are parapatrically distributed across the island and differ strikingly in plumage colour. One form restricted to the highlands is separated by a steep ecological gradient from three distinct lowland forms which meet at narrow hybrid zones that are not associated with environmental variables. Analyses of genomic variation based on single nucleotide polymorphism data from genotyping-by-sequencing and pooled RAD-seq approaches show that signatures of selection associated with elevation can be found at multiple regions across the genome, whereas most loci associated with the lowland forms are located on the Z sex chromosome. We identified TYRP1, a Z-linked colour gene, as a likely candidate locus underlying colour variation among lowland forms. Tests of demographic models revealed that highland and lowland forms diverged in the presence of gene flow, and divergence has progressed as gene flow was restricted by selection at loci across the genome. This system holds promise for investigating how adaptation and reproductive isolation shape the genomic landscape of divergence at multiple stages of the speciation process.

Cryptic and extensive hybridization between ancient lineages of American crows

Most species and therefore most hybrid zones have historically been defined using phenotypic characters. However, both speciation and hybridization can occur with negligible morphological differentiation. Recently developed genomic tools provide the means to better understand cryptic speciation and hybridization. The Northwestern Crow (Corvus caurinus) and American Crow (Corvus brachyrhynchos) are continuously distributed sister taxa that lack reliable traditional characters for identification. In this first population genomic study of Northwestern and American crows, we use genomic SNPs (nuDNA) and mtDNA to investigate the degree of genetic differentiation between these crows and the extent to which they may hybridize. Our results indicate that American and Northwestern crows have distinct evolutionary histories, supported by two nuDNA ancestry clusters and two 1.1%-divergent mtDNA clades dating to the late Pleistocene, when glacial advances may have isolated crow populations in separate refugia. We document extensive hybridization, with geographic overlap of mtDNA clades and admixture of nuDNA across >900 km of western Washington and western British Columbia. This broad hybrid zone consists of late-generation hybrids and backcrosses, but not recent (e.g., F1) hybrids. Nuclear DNA and mtDNA clines had concordant widths and were both centred in southwestern British Columbia, farther north than previously postulated. Overall, our results suggest a history of reticulate evolution in American and Northwestern crows, perhaps due to recurring neutral expansion(s) from Pleistocene glacial refugia followed by lineage fusion(s). However, we do not rule out a contributing role for more recent potential drivers of hybridization, such as expansion into human-modified habitats.

The Origin of a New Sex Chromosome by Introgression between Two Stickleback Fishes

Introgression is increasingly recognized as a source of genetic diversity that fuels adaptation. Its role in the evolution of sex chromosomes, however, is not well known. Here, we confirm the hypothesis that the Y chromosome in the ninespine stickleback, Pungitius pungitius, was established by introgression from the Amur stickleback, P. sinensis. Using whole genome resequencing, we identified a large region of Chr 12 in P. pungitius that is diverged between males and females. Within but not outside of this region, several lines of evidence show that the Y chromosome of P. pungitius shares a most recent common ancestor not with the X chromosome, but with the homologous chromosome in P. sinensis. Accumulation of repetitive elements and gene expression changes on the new Y are consistent with a young sex chromosome in early stages of degeneration, but other hallmarks of Y chromosomes have not yet appeared. Our findings indicate that porous species boundaries can trigger rapid sex chromosome evolution.

Historical climate changes and hybridization shaped the evolution of Atlantic Forest spinetails (Aves: Furnariidae)

Combining phylogeographic approaches and hybrid zone inference in a single framework is a robust way to depict respectively the biogeographic history of lineages and the evolutionary processes responsible for speciation. Here, we studied the spatiotemporal patterns of diversification and characterize the hybrid zone between two Atlantic Forest spinetails (Synallaxis ruficapilla and Synallaxis cinerea) using mitochondrial DNA and nuclear (autosomal and Z-linked) genes. We consistently recovered divergence between and within the two species during the late Pliocene and Pleistocene using an isolation with migration model. Also, our results indicate distinct levels of introgression among lineages. Ecological niche models and demographic inferences, used to infer range distributions throughout the late Quaternary, were not consistent with the hypothesis of a large river as a primary barrier responsible for the divergence of the two species. Instead, a scenario of isolation and divergence followed by geographic expansion and admixture as a consequence of Quaternary climatic oscillations was supported. Paleomodels also were not consistent with the idea that the hybrid zone originated in primary differentiation and favor a secondary contact scenario. Model fitting indicated that clines of different loci spanning the hybrid zone are coincident and concordant. The narrow cline for one Z-linked locus could be indicative of some form of post-zygotic selection hindering genetic homogenization between the two species.

Current geography masks dynamic history of gene flow during speciation in northern Australian birds

Genome divergence is greatly influenced by gene flow during early stages of speciation. As populations differentiate, geographic barriers can constrain gene flow and so affect the dynamics of divergence and speciation. Current geography, specifically disjunction and continuity of ranges, is often used to predict the historical gene flow during the divergence process. We test this prediction in eight meliphagoid bird species complexes codistributed in four regions. These regions are separated by known biogeographical barriers across northern Australia and Papua New Guinea. We find that bird populations currently separated by terrestrial habitat barriers within Australia and marine barriers between Australia and Papua New Guinea have a range of divergence levels and probability of gene flow not associated with current range connectivity. Instead, geographic distance and historical range connectivity better predict divergence and probability of gene flow. In this dynamic environmental context, we also find support for a nonlinear decrease of the probability of gene flow during the divergence process. The probability of gene flow initially decreases gradually after a certain level of divergence is reached. Its decrease then accelerates until the probability is close to zero. This implies that although geographic connectivity may have more of an effect early in speciation, other factors associated with higher divergence may play a more important role in influencing gene flow midway through and later in speciation. Current geographic connectivity may then mislead inferences regarding potential for gene flow during speciation under a complex and dynamic history of geographic and reproductive isolation.

Reduced Introgression of Sex Chromosome Markers in the Mexican Howler Monkey (Alouatta palliata × A. pigra) Hybrid Zone

Interspecific hybridization allows the introgression or movement of alleles from one genome to another. While some genomic regions freely exchange alleles during hybridization, loci associated with reproductive isolation do not intermix. In many model organisms, the X chromosome displays limited introgression compared to autosomes owing to the presence of multiple loci associated with hybrid sterility or inviability (the "large X-effect"). Similarly, if hybrids are produced, the heterogametic sex is usually inviable or sterile, a pattern known as Haldane's rule. We analyzed the patterns of introgression of genetic markers located in the mitochondrial (control region) and nuclear (autosomal microsatellites and sex chromosome genes) genomes of two howler monkey species (Alouatta palliata and A. pigra) that form a natural hybrid zone in southern Mexico, to evaluate whether the large X-effect and Haldane's rule affect the outcomes of hybridization between these sister species. To identify the level of admixture of each individual in the hybrid zone (N = 254) we analyzed individuals sampled outside the hybrid zone (109 A. pigra and 39 A. palliata) to determine allele frequencies of parental species and estimated a hybrid index based on nuclear markers. We then performed a cline analysis using individuals in the hybrid zone to determine patterns of introgression for each locus. Our analyses show that although the hybrid zone is bimodal (with no known F1 s and few recent generation hybrids) and quite narrow, there has been extensive introgression in both directions, and there is a large array of admixed individuals in the hybrid zone. Mitochondrial and most autosomal markers showed bidirectional introgression, but some had biased introgression toward one species or the other. All markers on the sex chromosomes and a few autosomal markers showed highly restricted introgression. This pattern is consistent with the hypothesis that the sex chromosomes make a disproportionate contribution to reproductive isolation, and our results broaden the taxonomic representation of these patterns across animal taxa.

Sex chromosome inversions enforce reproductive isolation across an avian hybrid zone

Across hybrid zones, the sex chromosomes are often more strongly differentiated than the autosomes. This is regularly attributed to the greater frequency of reproductive incompatibilities accumulating on sex chromosomes and their exposure in the heterogametic sex. Working within an avian hybrid zone, we explore the possibility that chromosome inversions differentially accumulate on the Z chromosome compared to the autosomes and thereby contribute to Z chromosome differentiation. We analyse the northern Australian hybrid zone between two subspecies of the long-tailed finch (Poephila acuticauda), first described based on differences in bill colour, using reduced-representation genomic sequencing for 293 individuals over a 1,530-km transect. Autosomal differentiation between subspecies is minimal. In contrast, 75% of the Z chromosome is highly differentiated and shows a steep genomic cline, which is displaced 350 km to the west of the cline in bill colour. Differentiation is associated with two or more putative chromosomal inversions, each predominating in one subspecies. If inversions reduce recombination between hybrid incompatibilities, they are selectively favoured and should therefore accumulate in hybrid zones. We argue that this predisposes inversions to differentially accumulate on the Z chromosome. One genomic region affecting bill colour is on the Z, but the main candidates are on chromosome 8. This and the displacement of the bill colour and Z chromosome cline centres suggest that bill colour has not strongly contributed to inversion accumulation. Based on cline width, however, the Z chromosome and bill colour both contribute to reproductive isolation established between this pair of subspecies.

Dobzhansky-Muller incompatibilities, dominance drive, and sex-chromosome introgression at secondary contact zones: A simulation study

Dobzhansky-Muller (DM) incompatibilities involving sex chromosomes have been proposed to account for Haldane's rule (lowered fitness among hybrid offspring of the heterogametic sex) as well as Darwin's corollary (asymmetric fitness costs with respect to the direction of the cross). We performed simulation studies of a hybrid zone to investigate the effects of different types of DM incompatibilities on cline widths and positions of sex-linked markers. From our simulations, X-Y incompatibilities generate steep clines for both X-linked and Y-linked markers; random effects may produce strong noise in cline center positions when migration is high relative to fitness costs, but X- and Y-centers always coincide strictly. X-autosome and Y-autosome incompatibilities also generate steep clines, but systematic shifts in cline centers occur when migration is high relative to selection, as a result of a dominance drive linked to Darwin's corollary. Interestingly, sex-linked genes always show farther introgression than the associated autosomal genes. We discuss ways of disentangling the potentially confounding effects of sex biases in migration, we compare our results to those of a few documented contact zones, and we stress the need to study independent replicates of the same contact zone.

Sex chromosomes and speciation in birds and other ZW systems

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

Interspecific competition, hybridization, and reproductive isolation in secondary contact: missing perspectives on males and females

Research on sexual selection and hybridization has focused on female mate choice and male-male competition. While the evolutionary outcomes of interspecific female preference have been well explored, we are now gaining a better understanding of the processes by which male-male competition between species in secondary contact promotes reproductive isolation versus hybridization. What is relatively unexplored is the interaction between female choice and male competition, as they can oppose one another or align with similar outcomes for reproductive isolation. The role of female-female competition in hybridization is also not well understood, but could operate similarly to male-male competition in polyandrous and other systems where costs to heterospecific mating are low for females. Reproductive competition between either sex of sympatric species can cause the divergence and/or convergence of sexual signals and recognition, which in turn influences the likelihood for interspecific mating. Future work on species interactions in secondary contact should test the relative influences of both mate choice and competition for mates on hybridization outcomes, and should not ignore the possibilities that females can compete over mating resources, and males can exercise mate choice.

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.

Reproductive isolation and environmental adaptation shape the phylogeography of mountain pine beetle (Dendroctonus ponderosae)

Chromosomal rearrangement can be an important mechanism driving population differentiation and incipient speciation. In the mountain pine beetle (MPB, Dendroctonus ponderosae), deletions on the Y chromosome that are polymorphic among populations are associated with reproductive incompatibility. Here, we used RAD sequencing across the entire MPB range in western North America to reveal the extent of the phylogeographic differences between Y haplotypes compared to autosomal and X-linked loci. Clustering and geneflow analyses revealed three distinct Y haplogroups geographically positioned within and on either side of the Great Basin Desert. Despite close geographic proximity between populations on the boundaries of each Y haplogroup, there was extremely low Y haplogroup mixing among populations, and gene flow on the autosomes was reduced across Y haplogroup boundaries. These results are consistent with a previous study suggesting that independent degradation of a recently evolved neo-Y chromosome in previously isolated populations causes male sterility or inviability among Y haplotype lineages. Phylogeographic results supported historic contraction of MPB into three separate Pleistocene glacial refugia followed by postglacial range expansion and secondary contact. Distinct sets of SNPs were statistically associated with environmental data among the most genetically distinct sets of geographic populations. This finding suggests that the process of adaptation to local climatic conditions is influenced by population genetic structure, with evidence for largely independent evolution in the most genetically isolated Y haplogroup.

The influence of molecular markers and methods on inferring the phylogenetic relationships between the representatives of the Arini (parrots, Psittaciformes), determined on the basis of their complete mitochondrial genomes

BACKGROUND

Conures are a morphologically diverse group of Neotropical parrots classified as members of the tribe Arini, which has recently been subjected to a taxonomic revision. The previously broadly defined Aratinga genus of this tribe has been split into the 'true' Aratinga and three additional genera, Eupsittula, Psittacara and Thectocercus. Popular markers used in the reconstruction of the parrots' phylogenies derive from mitochondrial DNA. However, current phylogenetic analyses seem to indicate conflicting relationships between Aratinga and other conures, and also among other Arini members. Therefore, it is not clear if the mtDNA phylogenies can reliably define the species tree. The inconsistencies may result from the variable evolution rate of the markers used or their weak phylogenetic signal. To resolve these controversies and to assess to what extent the phylogenetic relationships in the tribe Arini can be inferred from mitochondrial genomes, we compared representative Arini mitogenomes as well as examined the usefulness of the individual mitochondrial markers and the efficiency of various phylogenetic methods.

RESULTS

Single molecular markers produced inconsistent tree topologies, while different methods offered various topologies even for the same marker. A significant disagreement in these tree topologies occurred for cytb, nd2 and nd6 genes, which are commonly used in parrot phylogenies. The strongest phylogenetic signal was found in the control region and RNA genes. However, these markers cannot be used alone in inferring Arini phylogenies because they do not provide fully resolved trees. The most reliable phylogeny of the parrots under study is obtained only on the concatenated set of all mitochondrial markers. The analyses established significantly resolved relationships within the former Aratinga representatives and the main genera of the tribe Arini. Such mtDNA phylogeny can be in agreement with the species tree, owing to its match with synapomorphic features in plumage colouration.

CONCLUSIONS

Phylogenetic relationships inferred from single mitochondrial markers can be incorrect and contradictory. Therefore, such phylogenies should be considered with caution. Reliable results can be produced by concatenated sets of all or at least the majority of mitochondrial genes and the control region. The results advance a new view on the relationships among the main genera of Arini and resolve the inconsistencies between the taxa that were previously classified as the broadly defined genus Aratinga. Although gene and species trees do not always have to be consistent, the mtDNA phylogenies for Arini can reflect the species tree.

A new parrot taxon from the Yucatán Peninsula, Mexico-its position within genus Amazona based on morphology and molecular phylogeny

Parrots (Psittaciformes) are a diverse group of birds which need urgent protection. However, many taxa from this order have an unresolved status, which makes their conservation difficult. One species-rich parrot genus is Amazona, which is widely distributed in the New World. Here we describe a new Amazona form, which is endemic to the Yucatán Peninsula. This parrot is clearly separable from other Amazona species in eleven morphometric characters as well as call and behavior. The clear differences in these features imply that the parrot most likely represents a new species. In contrast to this, the phylogenetic tree based on mitochondrial markers shows that this parrot groups with strong support within A. albifrons from Central America, which would suggest that it is a subspecies of A. albifrons. However, taken together tree topology tests and morphometric analyses, we can conclude that the new parrot represents a recently evolving species, whose taxonomic status should be further confirmed. This lineage diverged from its closest relative about 120,000 years ago and was subjected to accelerated morphological and behavioral changes like some other representatives of the genus Amazona. Our phylogenies, which are so far the most comprehensive for Amazona taxa enabled us to consider the most feasible scenarios about parrot colonization of the Greater and Lesser Antilles and Central America from South America mainland. The molecular dating of these migrations and diversification rate were correlated with climatic and geological events in the last five million years, giving an interesting insight into Amazon parrot phylogeography and their evolution in general.

Nuclear introgression without mitochondrial introgression in two turtle species exhibiting sex-specific trophic differentiation

Despite the presence of reproductive barriers between species, interspecific gene introgression has been documented in a range of natural systems. Comparing patterns of genetic introgression in biparental versus matrilineal markers can potentially reveal sex-specific barriers to interspecific gene flow. Hybridization has been documented in the freshwater turtles Graptemys geographica and G. pseudogeographica, whose ranges are largely sympatric. Morphological differentiation between the species is restricted to females, with female G. geographica possessing large heads and jaws compared to the narrow heads of G. pseudogeographica females. If hybrid females are morphologically intermediate, they may be less successful at exploiting parental feeding niches, thereby limiting the introgression of maternally inherited, but not biparental, molecular markers. We paired sequence data with stable isotope analysis and examined sex-specific genetic introgression and trophic differentiation in sympatric populations of G. geographica and G. pseudogeographica. We observed introgression from G. pseudogeographica into G. geographica at three nuclear loci, but not at the mitochondrial locus. Analysis of ∂(15)N and ∂(13)C was consistent with species differences in trophic positioning in females, but not males. These results suggest that ecological divergence in females may reduce the opportunity for gene flow in this system.

Differential introgression and the maintenance of species boundaries in an advanced generation avian hybrid zone

Background

Evolutionary processes, including selection and differential fitness, shape the introgression of genetic material across a hybrid zone, resulting in the exchange of some genes but not others. Differential introgression of molecular or phenotypic markers can thus provide insight into factors contributing to reproductive isolation. We characterized patterns of genetic variation across a hybrid zone between two tidal marsh birds, Saltmarsh (Ammodramus caudacutus) and Nelson’s (A. nelsoni) sparrows (n = 286), and compared patterns of introgression among multiple genetic markers and phenotypic traits.

Results

Geographic and genomic cline analyses revealed variable patterns of introgression among marker types. Most markers exhibited gradual clines and indicated that introgression exceeds the spatial extent of the previously documented hybrid zone. We found steeper clines, indicating strong selection for loci associated with traits related to tidal marsh adaptations, including for a marker linked to a gene region associated with metabolic functions, including an osmotic regulatory pathway, as well as for a marker related to melanin-based pigmentation, supporting an adaptive role of darker plumage (salt marsh melanism) in tidal marshes. Narrow clines at mitochondrial and sex-linked markers also offer support for Haldane’s rule. We detected patterns of asymmetrical introgression toward A. caudacutus, which may be driven by differences in mating strategy or differences in population density between the two species.

Conclusions

Our findings offer insight into the dynamics of a hybrid zone traversing a unique environmental gradient and provide evidence for a role of ecological divergence in the maintenance of pure species boundaries despite ongoing gene flow.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0635-y) contains supplementary material, which is available to authorized users.

Empirical evidence for large X-effects in animals with undifferentiated sex chromosomes

Reproductive isolation is crucial for the process of speciation to progress. Sex chromosomes have been assigned a key role in driving reproductive isolation but empirical evidence from natural population processes has been restricted to organisms with degenerated sex chromosomes such as mammals and birds. Here we report restricted introgression at sex-linked compared to autosomal markers in a hybrid zone between two incipient species of European tree frog, Hyla arborea and H. orientalis, whose homologous X and Y sex chromosomes are undifferentiated. This large X-effect cannot result from the dominance or faster-X aspects of Haldane's rule, which are specific to degenerated sex chromosomes, but rather supports a role for faster-heterogametic-sex or faster-male evolutionary processes. Our data suggest a prominent contribution of undifferentiated sex chromosomes to speciation.

Biased introgression of mitochondrial and nuclear genes: a comparison of diploid and haplodiploid systems

Hybridization between recently diverged species, even if infrequent, can lead to the introgression of genes from one species into another. The rates of mitochondrial and nuclear introgression often differ, with some taxa showing biases for mitochondrial introgression and others for nuclear introgression. Several hypotheses exist to explain such biases, including adaptive introgression, sex differences in dispersal rates, sex-specific prezygotic isolation and sex-specific fitness of hybrids (e.g. Haldane's rule). We derive a simple population genetic model that permits an analysis of sex-specific demographic and fitness parameters and measures the relative rates of mitochondrial and nuclear introgression between hybridizing pairs. We do this separately for diploid and haplodiploid species. For diploid taxa, we recover results consistent with previous hypotheses: an excess of one sex among the hybridizing migrants or sex-specific prezygotic isolation causes a bias for one type of marker or the other; when Haldane's rule is obeyed, we find a mitochondrial bias in XY systems and a nuclear bias in ZW systems. For haplodiploid taxa, the model reveals that owing to their unique transmission genetics, they are seemingly assured of strong mitochondrial biases in introgression rates, unlike diploid taxa, where the relative fitness of male and female hybrids can tip the bias in either direction. This heretofore overlooked aspect of hybridization in haplodiploids provides what is perhaps the most likely explanation for differential introgression of mitochondrial and nuclear markers and raises concerns about the use of mitochondrial DNA barcodes for species delimitation in these taxa.

Dynamic karyotype evolution and unique sex determination systems in Leptidea wood white butterflies

Background

Chromosomal rearrangements have the potential to limit the rate and pattern of gene flow within and between species and thus play a direct role in promoting and maintaining speciation. Wood white butterflies of the genus Leptidea are excellent models to study the role of chromosome rearrangements in speciation because they show karyotype variability not only among but also within species. In this work, we investigated genome architecture of three cryptic Leptidea species (L. juvernica, L. sinapis and L. reali) by standard and molecular cytogenetic techniques in order to reveal causes of the karyotype variability.

Results

Chromosome numbers ranged from 2n = 85 to 91 in L. juvernica and 2n = 69 to 73 in L. sinapis (both from Czech populations) to 2n = 51 to 55 in L. reali (Spanish population). We observed significant differences in chromosome numbers and localization of cytogenetic markers (rDNA and H3 histone genes) within the offspring of individual females. Using FISH with the (TTAGG)_n telomeric probe we also documented the presence of multiple chromosome fusions and/or fissions and other complex rearrangements. Thus, the intraspecific karyotype variability is likely due to irregular chromosome segregation of multivalent meiotic configurations. The analysis of female meiotic chromosomes by GISH and CGH revealed multiple sex chromosomes: W₁W₂W₃Z₁Z₂Z₃Z₄ in L. juvernica, W₁W₂W₃Z₁Z₂Z₃ in L. sinapis and W₁W₂W₃W₄Z₁Z₂Z₃Z₄ in L. reali.

Conclusions

Our results suggest a dynamic karyotype evolution and point to the role of chromosomal rearrangements in the speciation of Leptidea butterflies. Moreover, our study revealed a curious sex determination system with 3–4 W and 3–4 Z chromosomes, which is unique in the Lepidoptera and which could also have played a role in the speciation process of the three Leptidea species.

Electronic supplementary material

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