Chromosomal speciation in mice: a cytogenetic analysis of recombination

Chromosome-level genomics and historical museum collections reveal new insights into the population structure and chromosome evolution of waterbuck

Advances in the sequencing and assembly of chromosome-level genome assemblies has enabled the study of non-model animals, providing further insights into the evolution of genomes and chromosomes. Here, we present the waterbuck ( Kobus ellipsiprymnus ) as an emerging model antelope for studying population dynamics and chromosome evolution. Antelope evolutionary history has been shaped by Robertsonian (Rb) fusions, with waterbuck also showing variation in karyotype due to two polymorphic Rb fusions. These polymorphisms are variable between and within the two recognised subspecies, the common and defassa waterbuck. To provide new insights into waterbuck evolution, we firstly assembled a chromosome-level genome assembly for the defassa subspecies using PacBio HiFi and Hi-C sequencing. We then utilised museum collections to carry out whole genome sequencing (WGS) of 24 historical waterbuck skins from both subspecies. Combined with a previous WGS dataset (n = 119), this represents the largest study of waterbuck populations to date. We found novel population structure and gene flow between waterbuck populations and regions across the genome with high genomic differentiation between the two subspecies. Several of these regions were found around the centromeres of fixed and polymorphic Rb fusions, exhibiting signatures of low recombination and local population structure. Interestingly, these regions contain genes involved in development, fertility, and recombination. Our results highlight the importance of assembling genomes to the chromosome-level, the utility and value of historical collections in sampling a wide-ranging species to uncover fine-scale population structure, and the potential impacts of Rb fusions on genomic differentiation and the recombination landscape.

Making sense of chromosome polymorphisms in two leptysmine grasshoppers

The touchstone of the 'New Synthesis' was population cytogenetics -rather than genetics - due to the abundant polymorphic inversions in the genus Drosophila. Grasshoppers were not a material of choice because of their conservative karyotypes. However, nowadays seven species of Acrididae were described for polymorphic centric fusions, five of them in South-America. Leptysma argentina and the likely biocontrol of water-hyacinth Cornops aquaticum are semiaquatic Leptysminae (Acrididae: Orthoptera), polymorphic for centric fusions, supernumerary segments and a B-chromosome. We sought to demonstrate the operation of natural selection on them, by detecting: (I) latitudinal clines; (II) regression on environmental variables; (III) deviation from null models, such as linkage equilibrium; (IV) seasonal variation; (V) comparison between age classes and (VI) selection component analyses. All of them were confirmed in L. argentina, just (I) and (II) in C. aquaticum. Furthermore, the relationship between karyotype, phenotype and recombination was confirmed in both species. Karyotype-phenotype relationship may be due to the body enlargement the fusions are associated with, along with a latitudinal transition in voltinism. Karyotype-related recombination reduction in both species may help explain all fusion clines, although there is probably more than one factor at work. No effects were noticed for a supernumerary segment in L. argentina, but it is ubiquitous and certainly non-neutral. C. aquaticum is poised for introduction in South-Africa as a biocontrol of water-hyacinths; the recent discovery of four more segment polymorphisms may imply more chromosomal markers to make sense of its genetic system.

Comparative karyotype analysis of the red brocket deer (M. americana sensu lato and M. rufa) complex: evidence of drastic chromosomal evolution and implications on speciation process

Chromosomal rearrangements are often associated with playing a role in the speciation process. However, the underlying mechanism that favors the genetic isolation associated with chromosomal changes remains elusive. In this sense, the genus Mazama is recognized by its high level of karyotype diversity among species with similar morphology. A cryptic species complex has been identified within the genus, with the red brocket deer (Mazama americana and Mazama rufa) being the most impressive example. The chromosome variation was clustered in cytotypes with diploid numbers ranging from 42 to 53 and was correlated with geographical location. We conducted an analysis of chromosome evolution of the red brocket deer complex using comparative chromosome painting and Bacterial Artificial Chromosome (BAC) clones among different cytotypes. The aim was to deepen our understanding of the karyotypic relationships within the red brocket, thereby elucidating the significant chromosome variation among closely related species. This underscores the significance of chromosome changes as a key evolutionary process shaping their genomes. The results revealed the presence of three distinct cytogenetic lineages characterized by significant karyotypic divergence, suggesting the existence of efficient post-zygotic barriers. Tandem fusions constitute the main mechanism driving karyotype evolution, following a few centric fusions, inversion X-autosomal fusions. The BAC mapping has improved our comprehension of the karyotypic relationships within the red brocket deer complex, prompting questions regarding the role of these changes in the speciation process. We propose the red brocket as a model group to investigate how chromosomal changes contribute to isolation and explore the implications of these changes in taxonomy and conservation.

Irregularities in Meiotic Prophase I as Prerequisites for Reproductive Isolation in Experimental Hybrids Carrying Robertsonian Translocations

The basic causes of postzygotic isolation can be elucidated if gametogenesis is studied, which is a drastically different process in males and females. As a step toward clarifying this problem, we obtained an experimental inbred lineage of the eastern mole vole Ellobius tancrei, whose founder animals were animals with identical diploid numbers 2n = 50 but with different Robertsonian translocations (Rb), namely 2Rb4.12 and 2Rb9.13 in the female and 2Rb.2.18 and 2Rb5.9 in the male. Here, we analyzed strictly inbred hybrids (F1, fertile and F10, sterile) using immunocytochemical methods in order to study spermatocytes during the meiotic prophase I. Previously, the presence of trivalents was assumed to have no significant effect on spermatogenesis and fertility in hybrids, but we demonstrated that spermatogenesis might be disturbed due to the cumulative effects of the retarded synapses of Rb bivalents as well as trivalents and their associations with XX sex bivalents. Alterations in the number of gametes due to the described processes led to a decrease in reproductive capacity up to sterility and can be examined as a mechanism for reproductive isolation, thus starting speciation.

Sex chromosome quadrivalents in oocytes of the African pygmy mouse Mus minutoides that harbors non-conventional sex chromosomes

Eutherian mammals have an extremely conserved sex-determining system controlled by highly differentiated sex chromosomes. Females are XX and males XY, and any deviation generally leads to infertility, mainly due to meiosis disruption. The African pygmy mouse (Mus minutoides) presents an atypical sex determination system with three sex chromosomes: the classical X and Y chromosomes and a feminizing X chromosome variant, called X*. Thus, three types of females coexist (XX, XX*, and X*Y) that all show normal fertility. Moreover, the three chromosomes (X and Y on one side and X* on the other side) are fused to different autosomes, which results in the inclusion of the sex chromosomes in a quadrivalent in XX* and X*Y females at meiotic prophase. Here, we characterized the configurations adopted by these sex chromosome quadrivalents during meiotic prophase. The XX* quadrivalent displayed a closed structure in which all homologous chromosome arms were fully synapsed and with sufficient crossovers to ensure the reductional segregation of all chromosomes at the first meiotic division. Conversely, the X*Y quadrivalents adopted either a closed configuration with non-homologous synapsis of the X* and Y chromosomes or an open chain configuration in which X* and Y remained asynapsed and possibly transcriptionally silenced. Moreover, the number of crossovers was insufficient to ensure chromosome segregation in a significant fraction of nuclei. Together, these findings raise questions about the mechanisms allowing X*Y females to have a level of fertility as good as that of XX and XX* females, if not higher.

Meiotic behavior of a complex hexavalent in heterozygous mice for Robertsonian translocations: insights for synapsis dynamics

Natural populations of the house mouse Mus musculus domesticus show great diversity in chromosomal number due to the presence of chromosomal rearrangements, mainly Robertsonian translocations. Breeding between two populations with different chromosomal configurations generates subfertile or sterile hybrid individuals due to impaired meiotic development. In this study, we have analyzed prophase-I spermatocytes of hybrids formed by crossing mice from Vulcano and Lipari island populations. Both populations have a 2n = 26 karyotype but different combinations of Robertsonian translocations. We studied the progress of synapsis, recombination, and meiotic silencing of unsynapsed chromosomes during prophase-I through the immunolocalization of the proteins SYCP3, SYCP1, γH2AX, RAD51, and MLH1. In these hybrids, a hexavalent is formed that, depending on the degree of synapsis between chromosomes, can adopt an open chain, a ring, or a closed configuration. The frequency of these configurations varies throughout meiosis, with the maximum degree of synapsis occurring at mid pachytene. In addition, we observed the appearance of heterologous synapsis between telocentric and metacentric chromosomes; however, this synapsis seems to be transient and unstable and unsynapsed regions are frequently observed in mid-late pachytene. Interestingly, we found that chiasmata are frequently located at the boundaries of unsynapsed chromosomal regions in the hexavalent during late pachytene. These results provide new clues about synapsis dynamics during meiosis. We propose that mechanical forces generated along chromosomes may induce premature desynapsis, which, in turn, might be counteracted by the location of chiasmata. Despite these and additional meiotic features, such as the accumulation of γH2AX on unsynapsed chromosome regions, we observed a large number of cells that progressed to late stages of prophase-I, indicating that synapsis defects may not trigger a meiotic crisis in these hybrids.

Chromosome Synapsis and Recombination in Male Hybrids between Two Chromosome Races of the Common Shrew (Sorex araneus L., Soricidae, Eulipotyphla)

Hybrid zones between chromosome races of the common shrew (Sorex araneus) provide exceptional models to study the potential role of chromosome rearrangements in the initial steps of speciation. The Novosibirsk and Tomsk races differ by a series of Robertsonian fusions with monobrachial homology. They form a narrow hybrid zone and generate hybrids with both simple (chain of three chromosomes) and complex (chain of eight or nine) synaptic configurations. Using immunolocalisation of the meiotic proteins, we examined chromosome pairing and recombination in males from the hybrid zone. Homozygotes and simple heterozygotes for Robertsonian fusions showed a low frequency of synaptic aberrations (<10%). The carriers of complex synaptic configurations showed multiple pairing abnormalities, which might lead to reduced fertility. The recombination frequency in the proximal regions of most chromosomes of all karyotypes was much lower than in the other regions. The strong suppression of recombination in the pericentromeric regions and co-segregation of race specific chromosomes involved in the long chains would be expected to lead to linkage disequilibrium between genes located there. Genic differentiation, together with the high frequency of pairing aberrations in male carriers of the long chains, might contribute to maintenance of the narrow hybrid zone.

A Half-Century of Studies on a Chromosomal Hybrid Zone of the House Mouse

The first natural chromosomal variation in the house mouse was described nearly 50 years ago in Val Poschiavo on the Swiss side of the Swiss-Italian border in the Central Eastern Alps. Studies have extended into neighboring Valtellina, and the house mice of the Poschiavo-Valtellina area have been subject to detailed analysis, reviewed here. The maximum extent of this area is 70 km, yet it has 4 metacentric races and the standard 40-chromosome telocentric race distributed in a patchwork fashion. The metacentric races are characterized by highly reduced diploid numbers (2n = 22-26) resulting from Robertsonian fusions, perhaps modified by whole-arm reciprocal translocations. The races hybridize and the whole Poschiavo-Valtellina area can be considered a "hybrid zone." The studies of this area have provided insights into origin of races within hybrid zones, gene flow within hybrid zones and the possibility of speciation in hybrid zones. This provides a case study of how chromosomal rearrangements may impact the genetic structure of populations and their diversification.

Conservation of genetic uniqueness of populations may increase extinction likelihood of endangered species: the case of Australian mammals

BACKGROUND

As increasingly fragmented and isolated populations of threatened species become subjected to climate change, invasive species and other stressors, there is an urgent need to consider adaptive potential when making conservation decisions rather than focussing on past processes. In many cases, populations identified as unique and currently managed separately suffer increased risk of extinction through demographic and genetic processes. Other populations currently not at risk are likely to be on a trajectory where declines in population size and fitness soon appear inevitable.

RESULTS

Using datasets from natural Australian mammal populations, we show that drift processes are likely to be driving uniqueness in populations of many threatened species as a result of small population size and fragmentation. Conserving and managing such remnant populations separately will therefore often decrease their adaptive potential and increase species extinction risk.

CONCLUSIONS

These results highlight the need for a paradigm shift in conservation biology practise; strategies need to focus on the preservation of genetic diversity at the species level, rather than population, subspecies or evolutionary significant unit. The introduction of new genetic variants into populations through in situ translocation needs to be considered more broadly in conservation programs as a way of decreasing extinction risk by increasing neutral genetic diversity which may increase the adaptive potential of populations if adaptive variation is also increased.

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.

Variation of Meiotic Recombination Rates and MLH1 Foci Distribution in Spermatocytes of Cattle, Sheep and Goats

Despite similar genome sizes, a great variability in recombination rates is observed in mammals. We used antibodies against SYCP3, MLH1 and centromeres to compare crossover frequency, position along chromosome arms and the effect of crossover interference in spermatocytes of 4 species from the family Bovidae (Bos taurus, 2n = 60, tribe Bovini; Ovis aries, 2n = 54, Capra hircus, 2n = 60 and Ammotragus lervia, 2n = 58, tribe Caprini). Despite significant individual variability, our results also show significant differences in both recombination rates and the total length of autosomal synaptonemal complexes (SC) between cattle (47.53 MLH1 foci/cell, 244.59 µm) and members of the tribe Caprini (61.83 MLH1 foci, 296.19 µm) which can be explained by the length of time that has passed since their evolutionary divergence. Sheep displayed the highest number of MLH1 foci per cell and recombination density, although they have a lower diploid chromosome number caused by centric fusions corresponding to cattle chromosomes 1;3, 2;8 and 5;11. However, the proportion of MLH1 foci observed on the fused chromosomes in sheep (26.14%) was significantly lower than on the orthologous acrocentrics in cattle (27.6%) and goats (28.2%), and their distribution along the SC arms differed significantly. The reduced recombination rate in metacentrics is probably caused by interference acting across the centromere.

Chromosomal polymorphism in mammals: an evolutionary perspective

Although chromosome rearrangements (CRs) are central to studies of genome evolution, our understanding of the evolutionary consequences of the early stages of karyotypic differentiation (i.e. polymorphism), especially the non-meiotic impacts, is surprisingly limited. We review the available data on chromosomal polymorphisms in mammals so as to identify taxa that hold promise for developing a more comprehensive understanding of chromosomal change. In doing so, we address several key questions: (i) to what extent are mammalian karyotypes polymorphic, and what types of rearrangements are principally involved? (ii) Are some mammalian lineages more prone to chromosomal polymorphism than others? More specifically, do (karyotypically) polymorphic mammalian species belong to lineages that are also characterized by past, extensive karyotype repatterning? (iii) How long can chromosomal polymorphisms persist in mammals? We discuss the evolutionary implications of these questions and propose several research avenues that may shed light on the role of chromosome change in the diversification of mammalian populations and species.

Effect of species-specific differences in chromosome morphology on chromatin compaction and the frequency and distribution of RAD51 and MLH1 foci in two bovid species: cattle (Bos taurus) and the common eland (Taurotragus oryx)

Meiotic recombination between homologous chromosomes is crucial for their correct segregation into gametes and for generating diversity. We compared the frequency and distribution of MLH1 foci and RAD51 foci, synaptonemal complex (SC) length and DNA loop size in two related Bovidae species that share chromosome arm homology but show an extreme difference in their diploid chromosome number: cattle (Bos taurus, 2n = 60) and the common eland (Taurotragus oryx, 2nmale = 31). Compared to cattle, significantly fewer MLH1 foci per cell were observed in the common eland, which can be attributed to the lower number of initial double-strand breaks (DSBs) detected as RAD51 foci in leptonema. Despite the significantly shorter total autosomal SC length and longer DNA loop size of the common eland bi-armed chromosomes compared to those of bovine acrocentrics, the overall crossover density in the common eland was still lower than in cattle, probably due to the reduction in the number of MLH1 foci in the proximal regions of the bi-armed chromosomes. The formation of centric fusions during karyotype evolution of the common eland accompanied by meiotic chromatin compaction has greater implications in the reduction in the number of DSBs in leptonema than in the decrease of MLH1 foci number in pachynema.

Impact of the number of Robertsonian chromosomes on germ cell death in wild male house mice

Previous studies in the house mouse have shown that the presence of Robertsonian (Rb) metacentric chromosomes in heterozygous condition affects the process of spermatogenesis. This detrimental effect mainly depends on the number of metacentrics involved and the complexity of the resulting meiotic figures. In this study, we aimed at elucidating the relationship between the chromosomal composition and spermatogenesis impairment in mice present in an area of chromosomal polymorphism (the so-called Barcelona system BRbS) in which Rb mice are surrounded by all acrocentric animals, no established metacentric races are present and the level of structural heterozygosity is relatively low. Using the terminal deoxinucleotidyl transferase-mediated dUTP nick end labelling (TUNEL) assay, we report higher frequency of apoptotic spermatogenetic cells in mice carrying six pairs of metacentrics at the homozygous state than in those carrying two or three fusions at the heterozygous state. Specifically, we detected a higher frequency of TUNEL-positive (T+) tubules and of T+ cells per tubule cross section and also a lower spermatid/spermatocyte ratio. These results indicate that the number of metacentrics at the homozygous state is more influential in determining apoptotic germ cell death than that of moderate chromosome heterozygosity. The percentage of germ cell death lower than 50 % found in our samples and the geographic distribution of the set of metacentrics within the BRbS indicate that although the spermatogenic alterations detected in this area could act as a partial barrier to gene flow, they are not sufficient to prevent Rb chromosomes from spreading in nature.

Genetic recombination variation in wild Robertsonian mice: on the role of chromosomal fusions and Prdm9 allelic background

Despite the existence of formal models to explain how chromosomal rearrangements can be fixed in a population in the presence of gene flow, few empirical data are available regarding the mechanisms by which genome shuffling contributes to speciation, especially in mammals. In order to shed light on this intriguing evolutionary process, here we present a detailed empirical study that shows how Robertsonian (Rb) fusions alter the chromosomal distribution of recombination events during the formation of the germline in a Rb system of the western house mouse (Mus musculus domesticus). Our results indicate that both the total number of meiotic crossovers and the chromosomal distribution of recombination events are reduced in mice with Rb fusions and that this can be related to alterations in epigenetic signatures for heterochromatinization. Furthermore, we detected novel house mouse Prdm9 allelic variants in the Rb system. Remarkably, mean recombination rates were positively correlated with a decrease in the number of ZnF domains in the Prdm9 gene. The suggestion that recombination can be modulated by both chromosomal reorganizations and genetic determinants that control the formation of double-stranded breaks during meiosis opens new avenues for understanding the role of recombination in chromosomal speciation.

Reduced recombination patterns in Robertsonian hybrids between chromosomal races of the house mouse: chiasma analyses

The recombination suppression models of chromosomal speciation posit that chromosomal rearrangements act as partial barriers to gene flow allowing these regions to accumulate genetic incompatibilities, thus contributing to the divergence of populations. Empirical and theoretical studies exploring the requirements of these models have mostly focused on the role of inversions. Here, the recombination landscape of heterozygosity for Robertsonian (Rb) fusions is investigated in the house mouse. Laboratory-bred F1 males and females between highly differentiated races from Tunisia (Rb: 2n=22, Standard, St: 2n=40) were produced in which all Rb fusions are present as trivalents in meiosis. Recombination patterns were determined by the analysis of chiasmata and compared with previous data on the Tunisian parental mice. A comparative analysis was performed on wild-caught male mice spanning the hybrid zone between two Italian races (2n=40, 2n=22). The results showed that the chiasma characteristics of both male and female Tunisian F1 and Italian hybrids clearly differed from those of Rb and St mice. Not only was the mean chiasma number (CN) intermediate between those of the parental mice in both geographic samples, but the distribution of chiasmata along the chromosomal arms of the F1 showed a distinct mosaic pattern. In short, the proximal region in the F1 exhibited a reduced CN similar to that observed in homozygous Rb, whereas distal regions more closely matched those in St mice. These results suggest that Rb rearrangements (homozygous or heterozygous) reduce recombination in the proximal regions of the chromosomes supporting their potential role in recombination-mediated speciation models.

The Robertsonian phenomenon in the house mouse: mutation, meiosis and speciation

Many different chromosomal races with reduced chromosome number due to the presence of Robertsonian fusion metacentrics have been described in western Europe and northern Africa, within the distribution area of the western house mouse Mus musculus domesticus. This subspecies of house mouse has become the ideal model for studies to elucidate the processes of chromosome mutation and fixation that lead to the formation of chromosomal races and for studies on the impact of chromosome heterozygosities on reproductive isolation and speciation. In this review, we briefly describe the history of the discovery of the first and subsequent metacentric races in house mice; then, we focus on the molecular composition of the centromeric regions involved in chromosome fusion to examine the molecular characteristics that may explain the great variability of the karyotype that house mice show. The influence that metacentrics exert on the nuclear architecture of the male meiocytes and the consequences on meiotic progression are described to illustrate the impact that chromosomal heterozygosities exert on fertility of house mice-of relevance to reproductive isolation and speciation. The evolutionary significance of the Robertsonian phenomenon in the house mouse is discussed in the final section of this review.

The role of chromosome variation in the speciation of the red brocket deer complex: the study of reproductive isolation in females

BACKGROUND

The red brocket deer, Mazama americana, has at least six distinct karyotypes in different regions of South America that suggest the existence of various species that are today all referred to as M. americana. From an evolutionary perspective, the red brockets are a relatively recent clade that has gone through intense diversification. This study sought to prove the existence of post-zygotic reproductive isolation in deer offspring between distinct chromosome lineages. To achieve this, inter-cytotype and intra-cytotype crosses were performed, which resulted in both F1 hybrid (n = 5) and pure offspring (n = 3) in captivity.

RESULTS

F1 females were analyzed in terms of their karyotypes, ovarian histology, estrous cycles and in vitro embryo production. Pure females presented parameters that were similar to those previously reported for M. Americana; however, the parameters for hybrid females were different. Two hybrids were determined to be sterile, while the remaining hybrids presented characteristics of subfertility.

CONCLUSIONS

The results support the existence of well-established reproductive isolation among the most distant karyotype lineages and elucidates the need to define all karyotype variants and their geographical ranges in order to define the number of species of red brocket.