Evolutionary Origins and Patterns of Diversification in Animal Brood Parasitism.

AbstractBrood parasitism involves the exploitation of host parental care rather than the extraction of resources directly from hosts. We identify defining characteristics of this strategy and consider its position along continua with adjacent behaviors but focus on canonical brood parasites, where parasitism is obligate and hosts are noneusocial (thereby distinguishing from social parasitism). A systematic literature survey revealed 59 independently derived brood parasitic lineages with most origins (49) in insects, particularly among bees and wasps, and other origins in birds (seven) and fish (three). Insects account for more than 98% of brood parasitic species, with much of that diversity reflecting ancient (≥100-million-year-old) brood parasitic lineages. Brood parasites usually, but not always, evolve from forms that show parental care. In insects, brood parasitism often first evolves through exploitation of a closely related species, following Emery's rule, but this is less typical in birds, which we discuss. We conducted lineage-level comparisons between brood parasitic clades and their sister groups, finding mixed results but an overall neutral to negative effect of brood parasitism on species richness and diversification. Our review of brood parasites reveals many unanswered questions requiring new research, including further modeling of the coevolutionary dynamics of brood parasites and their hosts.

Resource competition drives an invasion-replacement event among shrew species on an island.

Invasive mammals are responsible for the majority of native species extinctions on islands. While most of these extinction events will be due to novel interactions between species (e.g. exotic predators and naive prey), it is more unusual to find incidences where a newly invasive species causes the decline/extinction of a native species on an island when they normally coexist elsewhere in their overlapping mainland ranges. We investigated if resource competition between two insectivorous small mammals was playing a significant role in the rapid replacement of the native pygmy shrew Sorex minutus in the presence of the recently invading greater white-toothed shrew Crocidura russula on the island of Ireland. We used DNA metabarcoding of gut contents from >300 individuals of both species to determine each species' diet and measured the body size (weight and length) during different stages of the invasion in Ireland (before, during and after the species come into contact with one another) and on a French island where both species have long coexisted (acting as a natural 'control' site). Dietary composition, niche width and overlap and body size were compared in these different stages. The body size of the invasive C. russula and composition of its diet changes between when it first invades an area and after it becomes established. During the initial stages of the invasion, individual shrews are larger and consume larger sized invertebrate prey species. During later stages of the invasion, C. russula switches to consuming smaller prey taxa that are more essential for the native species. As a result, the level of interspecific dietary overlap increases from between 11% and 14% when they first come into contact with each other to between 39% and 46% after the invasion. Here we show that an invasive species can quickly alter its dietary niche in a new environment, ultimately causing the replacement of a native species. In addition, the invasive shrew could also be potentially exhausting local resources of larger invertebrate species. These subsequent changes in terrestrial invertebrate communities could have severe impacts further downstream on ecosystem functioning and services.

Tracking Chromosomal Origins in the Northern Italy System of Metacentric Races of the House Mouse.

The Western European house mouse is chromosomally diverse, with diploid karyotypes ranging from the standard 40 telocentric chromosomes down to 22 chromosomes. Karyotypes are modified through Robertsonian (Rb) fusion of 2 telocentrics into a single metacentric, occurring repeatedly with fixation, and whole-arm reciprocal translocations (WARTs) generating additional novel karyotypes. Over 100 metacentric populations (chromosomal races) have been identified, geographically clustered into "systems." Chromosomal races within systems often hybridise, and new races may emerge through this hybridisation ("zonal raciation"). We wished to determine the degree to which chromosomal races in a system have evolved independently or share common ancestry. Recombination between chromosomes from hybridising chromosomal races can erase the signals associated with a particular metacentric of interest, making inferences challenging. However, reduced recombination near the centromeres of chromosomal race-specific metacentrics makes centromere-adjacent markers ideal for solving this problem. For the Northern Italy System (NIS), we used microsatellite markers near the centromere to test previous hypotheses about evolutionary relationships of 5 chromosomal races. We chose markers from chromosomes 1, 3, 4, and 6, all of which comprise one arm of a metacentric in at least 2 of these NIS metacentric populations. We used estimates of FST and RST, as well as principal components analyses and neighbour-joining phylogenetic analyses, to infer evolutionary relationships between these 5 chromosomal races and neighbouring mice with the standard karyotype. We showed that the metacentric populations form a single grouping distinct from the standard populations, consistent with their common origin and consistent with a parsimonious sequence of chromosomal rearrangements to explain the relationship of the chromosomal races. That origin and evolution of the chromosomal races in the system would have involved Rb fusions, explaining the occurrence of chromosomal races with diploid numbers as low as 22. However, WARTs and zonal raciation have also been inferred, and the rare occurrence of chromosome 1 in different metacentrics in closely related chromosomal races is almost certainly explained by a WART. Our results with centromeric microsatellites are consistent with the above scenarios, illustrating, once again, the value of markers in the centromeric region to test evolutionary hypotheses in house mouse chromosomal systems.

Phylogenetic relationships and the evolution of host preferences in the largest clade of brood parasitic bees (Apidae: Nomadinae).

Brood parasites (also known as cleptoparasites) represent a substantial fraction of global bee diversity. Rather than constructing their own nests, these species instead invade those of host bees to lay their eggs. Larvae then hatch and consume the food provisions intended for the host's offspring. While this life history strategy has evolved numerous times across the phylogeny of bees, the oldest and most speciose parasitic clade is the subfamily Nomadinae (Apidae). However, the phylogenetic relationships among brood parasitic apids both within and outside the Nomadinae have not been fully resolved. Here, we present new findings on the phylogeny of this diverse group of brood parasites based on ultraconserved element (UCE) sequence data and extensive taxon sampling with 114 nomadine species representing all tribes. We suggest a broader definition of the subfamily Nomadinae to describe a clade that includes almost all parasitic members of the family Apidae. The tribe Melectini forms the sister group to all other Nomadinae, while the remainder of the subfamily is composed of two sister clades: a "nomadine line" representing the former Nomadinae sensu stricto, and an "ericrocidine line" that unites several mostly Neotropical lineages. We find the tribe Osirini Handlirsch to be polyphyletic, and divide it into three lineages, including the newly described Parepeolini trib. nov. In addition to our taxonomic findings, we use our phylogeny to explore the evolution of different modes of parasitism, detecting two independent transitions from closed-cell to open-cell parasitism. Finally, we examine how nomadine host-parasite associations have evolved over time. In support of Emery's rule, which suggests close relationships between hosts and parasites, we confirm that the earliest nomadines were parasites of their close free-living relatives within the family Apidae, but that over time their host range broadened to include more distantly related hosts spanning the diversity of bees. This expanded breadth of host taxa may also be associated with the transition to open-cell parasitism.

The complete mitochondrial genome of Talpa martinorum (Mammalia: Talpidae), a mole species endemic to Thrace: genome content and phylogenetic considerations.

The complete mitogenome sequence of Talpa martinorum, a recently described Balkan endemic mole, was assembled from next generation sequence data. The mitogenome is similar to that of the three other Talpa species sequenced to date, being 16,835 bp in length, and containing 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, an origin of L-strand replication, and a control region or D-loop. Compared to other Talpa mitogenomes sequenced to date, that of T. martinorum differs in the length of D-loop and stop codon usage. TAG and T-- are the stop codons for the ND1 and ATP8 genes, respectively, in T. martinorum, whilst TAA acts as a stop codon for both ND1 and ATP8 in the other three Talpa species sequenced. Phylogeny reconstructions based on Maximum Likelihood and Bayesian inference analyses yielded phylogenies with similar topologies, demonstrating that T. martinorum nests within the western lineage of the genus, being closely related to T. aquitania and T. occidentalis.

Population structure and inbreeding in wild house mice (Mus musculus) at different geographic scales.

House mice (Mus musculus) have spread globally as a result of their commensal relationship with humans. In the form of laboratory strains, both inbred and outbred, they are also among the most widely used model organisms in biomedical research. Although the general outlines of house mouse dispersal and population structure are well known, details have been obscured by either limited sample size or small numbers of markers. Here we examine ancestry, population structure, and inbreeding using SNP microarray genotypes in a cohort of 814 wild mice spanning five continents and all major subspecies of Mus, with a focus on M. m. domesticus. We find that the major axis of genetic variation in M. m. domesticus is a south-to-north gradient within Europe and the Mediterranean. The dominant ancestry component in North America, Australia, New Zealand, and various small offshore islands are of northern European origin. Next we show that inbreeding is surprisingly pervasive and highly variable, even between nearby populations. By inspecting the length distribution of homozygous segments in individual genomes, we find that inbreeding in commensal populations is mostly due to consanguinity. Our results offer new insight into the natural history of an important model organism for medicine and evolutionary biology.

Genome size influences adaptive plasticity of water loss, but not metabolic rates, in lungless salamanders.

Many expressions of phenotype, such as physiological performance, integrate multiple underlying traits to function. Linking component traits to adaptive physiology thus gives insight into mechanisms of selection acting on performance. Genome size (C-value) is a trait that influences physiology in multiple taxa by exerting a nucleotypic effect, constraining cell size and cellular physiology such that whole-organism mass-specific metabolism is reduced with increasing C-value. We tested for this mechanism of C-value function acting in lungless salamanders, plus an unexplored potential mechanism of C-value effects constraining water transport across the body surface to influence cutaneous water loss rates. We found no evidence for a nucleotypic effect on metabolic rates, but we demonstrate a relationship between C-value and water loss physiology. Under warmer experimental conditions, C-value was inversely correlated with water loss and positively correlated with resistance to water loss, which demonstrated adaptive plasticity at higher temperatures. We hypothesize that this pattern results from differences in cell size constraining diffusion and evaporation of water from the skin under warm conditions when cutaneous perfusion is reduced. Testing this hypothesis may confirm a previously unappreciated adaptive role for C-value variation in this group, and reveals the possibility that genome size influences physiological exchange across transport barriers more broadly.

PRDM9 Diversity at Fine Geographical Scale Reveals Contrasting Evolutionary Patterns and Functional Constraints in Natural Populations of House Mice.

One of the major challenges in evolutionary biology is the identification of the genetic basis of postzygotic reproductive isolation. Given its pivotal role in this process, here we explore the drivers that may account for the evolutionary dynamics of the PRDM9 gene between continental and island systems of chromosomal variation in house mice. Using a data set of nearly 400 wild-caught mice of Robertsonian systems, we identify the extent of PRDM9 diversity in natural house mouse populations, determine the phylogeography of PRDM9 at a local and global scale based on a new measure of pairwise genetic divergence, and analyze selective constraints. We find 57 newly described PRDM9 variants, this diversity being especially high on Madeira Island, a result that is contrary to the expectations of reduced variation for island populations. Our analysis suggest that the PRDM9 allelic variability observed in Madeira mice might be influenced by the presence of distinct chromosomal fusions resulting from a complex pattern of introgression or multiple colonization events onto the island. Importantly, we detect a significant reduction in the proportion of PRDM9 heterozygotes in Robertsonian mice, which showed a high degree of similarity in the amino acids responsible for protein-DNA binding. Our results suggest that despite the rapid evolution of PRDM9 and the variability detected in natural populations, functional constraints could facilitate the accumulation of allelic combinations that maintain recombination hotspot symmetry. We anticipate that our study will provide the basis for examining the role of different PRDM9 genetic backgrounds in reproductive isolation in natural populations.

High genomic diversity in the bank vole at the northern apex of a range expansion: The role of multiple colonizations and end-glacial refugia.

The history of repeated northern glacial cycling and southern climatic stability has long dominated explanations for how genetic diversity is distributed within temperate species in Eurasia and North America. However, growing evidence indicates the importance of cryptic refugia for northern colonization dynamics. An important geographic region to assess this is Fennoscandia, where recolonization at the end of the last glaciation was restricted to specific routes and temporal windows. We used genomic data to analyse genetic diversity and colonization history of the bank vole (Myodes glareolus) throughout Europe (>800 samples) with Fennoscandia as the northern apex. We inferred that bank voles colonized Fennoscandia multiple times by two different routes; with three separate colonizations via a southern land-bridge route deriving from a "Carpathian" glacial refugium and one via a north-eastern route from an "Eastern" glacial refugium near the Ural Mountains. Clustering of genome-wide SNPs revealed high diversity in Fennoscandia, with eight genomic clusters: three of Carpathian origin and five Eastern. Time estimates revealed that the first of the Carpathian colonizations occurred before the Younger Dryas (YD), meaning that the first colonists survived the YD in Fennoscandia. Results also indicated that introgression between bank and northern red-backed voles (Myodes rutilus) took place in Fennoscandia just after end-glacial colonization. Therefore, multiple colonizations from the same and different cryptic refugia, temporal and spatial separations and interspecific introgression have shaped bank vole genetic variability in Fennoscandia. Together, these processes drive high genetic diversity at the apex of the northern expansion in this emerging model species.

Measurements of hybrid fertility and a test of mate preference for two house mouse races with massive chromosomal divergence.

BACKGROUND: Western house mice Mus musculus domesticus are among the most important mammalian model species for chromosomal speciation. Hybrids between chromosomal races of M. m. domesticus suffer various degrees of fertility reduction between full fertility and complete sterility, depending on the complexity of the chromosomal differences between the races. This complexity presents itself in hybrids as meiotic configurations of chromosome chains and rings, with longer configurations having a stronger impact on fertility. While hybrids with short configurations have been intensively studied, less work has been done on hybrids with very long configurations. In this study, we investigated laboratory-reared wild mice from two chromosomally very different races in Switzerland found in close proximity. Hybrids between these races form a meiotic chain of fifteen chromosomes. We performed a detailed analysis of male and female hybrid fertility, including three generations of female backcrosses to one of the parental races. We also tested for possible divergence of mate preference in females. RESULTS: While all male F1 hybrids were sterile with sperm counts of zero, 48% of female F1 hybrids produced offspring. Their litter sizes ranged from one to three which is significantly lower than the litter size of parental race females. When hybrid females were backcrossed to a parental race, half of the offspring resembled the parental race in karyotype and fertility, while the other half resembled the F1 hybrids. In the preference test, females of both races indicated a lack of a preference for males of their own karyotype. CONCLUSIONS: Although the fertility of the F1 hybrids was extremely low because of the complexity of the chromosomal differences between the races, reproductive isolation was not complete. As we did not find assortative female preferences, we expect that contact between these races would lead to the production of hybrids and that gene flow would occur eventually, as fertility can be restored fully after one backcross generation.

Isotope values of California vole (Microtus californicus) hair relate to historical drought and land use patterns in California, USA.

Increased drought frequency and intensity and agricultural intensification have been key stressors to ecological systems over the past century. Biological proxies (e.g., pollen, tree rings) have been used to track this environmental change; however, linking these changes to the ecology of organisms remains challenging. Here, we link historical drought records to conditions of high water-stress in grassland habitats through the stable isotope analysis of California vole museum specimens (Microtus californicus). Using museum collections spanning 118-years (1891-2009), isotope values of dated hair tissues were associated with statewide drought metrics on the Palmer Drought Severity Index. We observed a positive correlation between δ15N and δ18O values and drought severity. The range in δ15N values (~ 18‰) is greater than what would be expected as a result of dietary shifts across the landscape (~ 3‰), and is likely attributed to the combined effects of physiological responses of M. californicus and isotopic shifts in plant resources with increased water-stress. Geospatial patterns in δ34S values of hair tissues reflect higher baseline isotope values in coastal habitats. However, comparably high δ34S values in the southern-most inland localities suggest sulfur fertilization of croplands and subsequent transfer to surrounding grassland habitats in 34S enriched forms. A broad δ13C range (- 28.7 to - 14.3‰) further suggests the consumption of C3 and C4 plant-based dietary proteins. As shown here, stable isotope analysis of museum collections can provide a climate and land use record based on the physiological performance and ecology of a study species in a region affected intensely by anthropogenic activities.

Genomics of end-Pleistocene population replacement in a small mammal.

Current species distributions at high latitudes are the product of expansion from glacial refugia into previously uninhabitable areas at the end of the last glaciation. The traditional view of postglacial colonization is that southern populations expanded their ranges into unoccupied northern territories. Recent findings on mitochondrial DNA (mtDNA) of British small mammals have challenged this simple colonization scenario by demonstrating a more complex genetic turnover in Britain during the Pleistocene-Holocene transition where one mtDNA clade of each species was replaced by another mtDNA clade of the same species. Here, we provide evidence from one of those small mammals, the bank vole (Clethrionomys glareolus), that the replacement was genome-wide. Using more than 10 000 autosomal SNPs we found that similar to mtDNA, bank vole genomes in Britain form two (north and south) clusters which admix. Therefore, the genome of the original postglacial colonists (the northern cluster) was probably replaced by another wave of migration from a different continental European population (the southern cluster), and we gained support for this by modelling with approximate Bayesian computation. This finding emphasizes the importance of analysis of genome-wide diversity within species under changing climate in creating opportunities for sophisticated testing of population history scenarios.

Synchronous diversification of Sulawesi's iconic artiodactyls driven by recent geological events.

The high degree of endemism on Sulawesi has previously been suggested to have vicariant origins, dating back to 40 Ma. Recent studies, however, suggest that much of Sulawesi's fauna assembled over the last 15 Myr. Here, we test the hypothesis that more recent uplift of previously submerged portions of land on Sulawesi promoted diversification and that much of its faunal assemblage is much younger than the island itself. To do so, we combined palaeogeographical reconstructions with genetic and morphometric datasets derived from Sulawesi's three largest mammals: the babirusa, anoa and Sulawesi warty pig. Our results indicate that although these species most likely colonized the area that is now Sulawesi at different times (14 Ma to 2-3 Ma), they experienced an almost synchronous expansion from the central part of the island. Geological reconstructions indicate that this area was above sea level for most of the last 4 Myr, unlike most parts of the island. We conclude that emergence of land on Sulawesi (approx. 1-2 Myr) may have allowed species to expand synchronously. Altogether, our results indicate that the establishment of the highly endemic faunal assemblage on Sulawesi was driven by geological events over the last few million years.

Integrative approaches to guide conservation decisions: Using genomics to define conservation units and functional corridors.

Climate change and increasing habitat loss greatly impact species survival, requiring range shifts, phenotypic plasticity and/or evolutionary change for long-term persistence, which may not readily occur unaided in threatened species. Therefore, defining conservation actions requires a detailed assessment of evolutionary factors. Existing genetic diversity needs to be thoroughly evaluated and spatially mapped to define conservation units (CUs) in an evolutionary context, and we address that here. We also propose a multidisciplinary approach to determine corridors and functional connectivity between CUs by including genetic diversity in the modelling while controlling for isolation by distance and phylogeographic history. We evaluate our approach on a Near Threatened Iberian endemic rodent by analysing genotyping-by-sequencing (GBS) genomic data from 107 Cabrera voles (Microtus cabrerae), screening the entire species distribution to define categories of CUs and their connectivity: We defined six management units (MUs) which can be grouped into four evolutionarily significant units (ESUs) and three (putatively) adaptive units (AUs). We demonstrate that the three different categories of CU can be objectively defined using genomic data, and their characteristics and connectivity can inform conservation decision-making. In particular, we show that connectivity of the Cabrera vole is very limited in eastern Iberia and that the pre-Pyrenean and part of the Betic geographic nuclei contribute the most to the species genetic diversity. We argue that a multidisciplinary framework for CU definition is essential and that this framework needs a strong evolutionary basis.

The Microtus voles: Resolving the phylogeny of one of the most speciose mammalian genera using genomics.

Sequential rapid radiations pose some of the greatest difficulties in phylogenetics, especially when analysing only a small number of genetic markers. Given that most of the speciation events occur in quick succession at various points in time, this creates particular challenges in determining phylogenetic relationships, i.e. branching order and divergence times. With the development of high throughput sequencing, thousands of markers can now readily be used to tackle these issues. Microtus is a speciose genus currently composed of 65 species that evolved over the last 2 million years. Although it is a well-studied group, there is still phylogenetic uncertainty at various divergence levels. Building upon previous studies that generally used small numbers of mitochondrial and/or nuclear loci, in this genomic-scale study we used both mitochondrial and nuclear data to study the rapid radiation within Microtus, using partial mitogenomes and genotyping-by-sequencing (GBS) on seven species representing five Microtus subgenera and the main biogeographic ranges where this group occurs. Both types of genome (mitochondrial and nuclear) generated similar tree topologies, with a basal split of the Nearctic (M. ochrogaster) and Holarctic (M. oeconomus) species, and then a subdivision of the five Palearctic species into two subgroups. These data support the occurrence of two European radiations, one North American radiation, and a later expansion of M. oeconomus from Asia to both Europe and North America. We further resolved the positioning of M. cabrerae as sister group of M. agrestis and refute the claim that M. cabrerae should be elevated to its own genus (Iberomys). Finally, the data support ongoing speciation events, especially within M. agrestis, with high levels of genetic divergence between the three Evolutionarily Significant Units (ESUs) previously identified. Similar high levels of divergence were also found among ESUs within M. oeconomus and M. arvalis.

Sperm competition as an under-appreciated factor in domestication.

Humans created an environment that increased selective pressures on subgroups of those species that became domestic. We propose that the domestication process may in some cases have been facilitated by changes in mating behaviour and resultant sperm competition. By adapting to sperm competition, proto-domestic animals could potentially have outcompeted their wild counterparts in human-constructed niches. This could have contributed to the restriction of gene flow between the proto-domesticates and their wild counterparts, thereby promoting the fixation of other domestication characteristics. Further to this novel perspective for domestication, we emphasize the general potential of postcopulatory sexual selection in the restriction of gene flow between populations, and urge more studies.

Exploring Rattus praetor (Rodentia, Muridae) as a possible species complex using geometric morphometrics on dental morphology.

Taxonomic uncertainties in the Rattus genus persist due to among-species morphological conservatism coupled with within-species environmental variation in morphology. As a result, this genus contains a number of possible cryptic species. One important example can be found in R. praetor, where morphological studies indicate it is a possible species complex. Genetic studies of R. praetor (limited to analysis of mitochondrial DNA) have been inconclusive, but do indicate such subdivision. Here we use geometric morphometrics to explore this possible species complex by analysing the dental traits of 48 specimens from New Guinea and neighbouring regions. We find separate molar morphologies for Bougainsville Island, central New Guinea and west New Guinea which cannot be easily explained by different environmental factors (climate, precipitation and altitude), strongly suggesting the existence of a number of evolutionarily distinct taxa within what is currently called R. praetor thus supporting previous suggestions that R. praetor is a species complex. Our findings demonstrate the potential of advanced morphological analyses in identifying separate species, contrary to the claims of morphological conservatism. Future analyses should combine geometric morphometrics with genetic analyses over the species range and include sub-fossil specimens from the Bismarck archipelago and Solomon Islands to resolve the evolutionary history of R. praetor.

R2d2 Drives Selfish Sweeps in the House Mouse.

A selective sweep is the result of strong positive selection driving newly occurring or standing genetic variants to fixation, and can dramatically alter the pattern and distribution of allelic diversity in a population. Population-level sequencing data have enabled discoveries of selective sweeps associated with genes involved in recent adaptations in many species. In contrast, much debate but little evidence addresses whether "selfish" genes are capable of fixation-thereby leaving signatures identical to classical selective sweeps-despite being neutral or deleterious to organismal fitness. We previously described R2d2, a large copy-number variant that causes nonrandom segregation of mouse Chromosome 2 in females due to meiotic drive. Here we show population-genetic data consistent with a selfish sweep driven by alleles of R2d2 with high copy number (R2d2(HC)) in natural populations. We replicate this finding in multiple closed breeding populations from six outbred backgrounds segregating for R2d2 alleles. We find that R2d2(HC) rapidly increases in frequency, and in most cases becomes fixed in significantly fewer generations than can be explained by genetic drift. R2d2(HC) is also associated with significantly reduced litter sizes in heterozygous mothers, making it a true selfish allele. Our data provide direct evidence of populations actively undergoing selfish sweeps, and demonstrate that meiotic drive can rapidly alter the genomic landscape in favor of mutations with neutral or even negative effects on overall Darwinian fitness. Further study will reveal the incidence of selfish sweeps, and will elucidate the relative contributions of selfish genes, adaptation and genetic drift to evolution.

Genetic differentiation within and away from the chromosomal rearrangements characterising hybridising chromosomal races of the western house mouse (Mus musculus domesticus).

The importance of chromosomal rearrangements for speciation can be inferred from studies of genetic exchange between hybridising chromosomal races within species. Reduced fertility or recombination suppression in karyotypic hybrids has the potential to maintain or promote genetic differentiation in genomic regions near rearrangement breakpoints. We studied genetic exchange between two hybridising groups of chromosomal races of house mouse in Upper Valtellina (Lombardy, Italy), using microsatellites. These groups differ by Robertsonian fusions and/or whole-arm reciprocal translocations such that F1 hybrids have a chain-of-five meiotic configuration. Previous studies showed genetic differentiation in two chromosomes in the chain-of-five (10 and 12) close to their centromeres (i.e. the rearrangement breakpoints); we have shown here that the centromeric regions of the other two chromosomes in the chain (2 and 8) are similarly differentiated. The internal chromosomes of the chain (8 and 12) show the greatest differentiation, which may reflect pairing and recombination properties of internal and external elements in a meiotic chain. Importantly, we found that centromeric regions of some non-rearranged chromosomes also showed genetic differentiation between the hybridising groups, indicating a complex interplay between chromosomal rearrangements and other parts of the genome in maintaining or promoting differentiation and potentially driving speciation between chromosomal races.