- Research article
Are ribosomal DNA clusters rearrangement hotspots? A case study in the genus Mus (Rodentia, Muridae)
BMC Evolutionary Biologyvolume 11, Article number: 124 (2011)
Recent advances in comparative genomics have considerably improved our knowledge of the evolution of mammalian karyotype architecture. One of the breakthroughs was the preferential localization of evolutionary breakpoints in regions enriched in repetitive sequences (segmental duplications, telomeres and centromeres). In this context, we investigated the contribution of ribosomal genes to genome reshuffling since they are generally located in pericentromeric or subtelomeric regions, and form repeat clusters on different chromosomes. The target model was the genus Mus which exhibits a high rate of karyotypic change, a large fraction of which involves centromeres.
The chromosomal distribution of rDNA clusters was determined by in situ hybridization of mouse probes in 19 species. Using a molecular-based reference tree, the phylogenetic distribution of clusters within the genus was reconstructed, and the temporal association between rDNA clusters, breakpoints and centromeres was tested by maximum likelihood analyses. Our results highlighted the following features of rDNA cluster dynamics in the genus Mus: i) rDNA clusters showed extensive diversity in number between species and an almost exclusive pericentromeric location, ii) a strong association between rDNA sites and centromeres was retrieved which may be related to their shared constraint of concerted evolution, iii) 24% of the observed breakpoints mapped near an rDNA cluster, and iv) a substantial rate of rDNA cluster change (insertion, deletion) also occurred in the absence of chromosomal rearrangements.
This study on the dynamics of rDNA clusters within the genus Mus has revealed a strong evolutionary relationship between rDNA clusters and centromeres. Both of these genomic structures coincide with breakpoints in the genus Mus, suggesting that the accumulation of a large number of repeats in the centromeric region may contribute to the high level of chromosome repatterning observed in this group. However, the elevated rate of rDNA change observed in the chromosomally invariant clade indicates that the presence of these sequences is insufficient to lead to genome instability. In agreement with recent studies, these results suggest that additional factors such as modifications of the epigenetic state of DNA may be required to trigger evolutionary plasticity.
The extensive advances in comparative cytogenomics in recent years have considerably enhanced our knowledge of the evolution of mammalian genomic architecture . In particular, the identification of syntenic associations of homologous chromosomal segments has led to the delimitation of breakpoint sites, i.e., regions where genome synteny has been disrupted by chromosomal rearrangements. Complementary to the cytogenetic approach, whole genome comparisons between different mammalian groups have highlighted the extensive reuse of breakpoints during chromosome evolution supporting the notion of evolutionary breakpoints as hotspots of genome repatterning . This predisposition of certain genomic regions to instability has been formalized as the 'fragile breakage model' . Moreover, comparative genomics, particularly in primates, has provided greater resolution of breakpoints. These studies have shown that breakpoint regions are enriched with different types of repetitive sequences such as SINEs, LINEs, LTRs and, in particular, segmental duplications [2, 4–6]. The presence of repetitive sequences at evolutionary breakpoints is thought to be related to the role tandem repeats play as a substrate for non-homologous recombination (i.e., exchanges between two different chromosomes) thereby promoting chromosomal rearrangements . However, several analyses suggest that rearrangement breakpoints and repeat sequences (flanking regions of segmental duplications) share similar physicochemical properties (high DNA flexibility and low stability) characteristic of fragile sites [1, 8]. Centromeric and subtelomeric domains appear as the key regions in chromosome evolution, since both accumulate repeat sequences, harbour many breakpoints and engage in non-homologous recombination [9, 10]. Several studies have shown a significant association between breakpoint reuse and centromere repositioning in different mammalian groups such as Marsupialia, Muridae, Equidae, and Primates [2, 11–13], underscoring the pivotal role of these genomic regions in chromosome restructuring.
The rRNA genes represent another family of tandem repeat sequences. These genes code for the ribosomal subunits that are essential for the cellular translation machinery. In mammals, each unit is composed of three genes coding for 18S, 5.8S and 28S ribosomal RNA; these genes are separated by two intergenic spacers and an external transcribed spacer . Whereas the sequences of spacers are generally highly divergent, the ribosomal coding elements show in some regions a remarkable sequence conservation within species and even among distantly related organisms . The tandem repeats of units are further organized into clusters present on one to several chromosome pairs where they are most often located in pericentromeric or subtelomeric regions [4, 16–18]. Several studies have highlighted the species-specific chromosomal distribution of rDNA clusters even between closely related species, suggesting that the location of rDNA clusters can rapidly change through transposition. The rDNA clusters therefore show several features in common with breakpoint regions: they are tandemly repeated; they are generally located in pericentromeric and subtelomeric regions; they transpose; they are subject to high rates of intra- and inter-chromosomal recombination. Indeed, the coincidence of rDNA clusters with chromosomal breakpoints has been highlighted in vitro in plants . Here, we investigate the evolutionary dynamics of rDNA clusters in relation to genome reshuffling in an emblematic mammalian model, the genus Mus which includes the house mouse.
Mus is an ideal biological model to investigate processes of chromosomal evolution since it exhibits one of the highest rates of karyotypic repatterning documented in mammals . The genus is species-rich with more than 40 species distributed among four subgenera . These include the subgenus Coelomys (shrew mice) with four species and a South-East Asian distribution, the subgenus Pyromys (spiny mice) with five species restricted to the Indian subcontinent, and the African subgenus Nannomys (pygmy mice) the most species-rich (with 18 species) and karyotypically diverse of the four subgenera [21, 22]. Finally, the Eurasian subgenus Mus comprises 16 species, one of which is the house mouse [21, 23]. Comparative cytogenomics has revealed a 10- to 30-fold acceleration in chromosomal change that was coincidental with the subgeneric cladogenesis, followed by a remarkable stasis in the subgenus Mus in which all 16 species share the same 2n = 40 karyotype. Moreover, studies have shown that more than half of the observed rearrangements involved centromeres (reactivation of latent sites or neocentromerization events) [13, 24].
The aim of the present study was to investigate if the evolutionary dynamics of rDNA clusters played a role in the extensive genomic reshuffling of the genus Mus. To do so, we first established the variation in location and number of the rDNA clusters between species belonging to different subgenera. Second, we used the available molecular phylogenies to reconstruct and track the evolution of these clusters in the genus. Finally, building upon the existing chromosomal phylogeny, we tested the temporal association between rDNA clusters, breakpoint regions, and centromeres.
The 19 species and subspecies as well as the number of specimens studied are listed in Table 1. Animals were either obtained from the Conversatoire Génétique de la Souris Sauvage (Institut des Sciences de l'Evolution, Montpellier, France) or collected in the wild. Four additional species (Mus booduga, M. terricolor, Apodemus sylvaticus, Rattus rattus) for which rDNA cluster data were available, were included in the analyses [25–27].
For all taxa, mitotic metaphases were obtained by the air-drying method from bone marrow cells after yeast stimulation . Identification of chromosomes was performed by DAPI-banding following the nomenclature of Cowell  for the subgenus Mus, Veyrunes et al. [22, 24] for the subgenera Nannomys and Coelomys, and Matsubara et al.  for the subgenus Pyromys. At least five metaphases per specimen were analysed. All observations were made with a Zeiss Axiophot fluorescence microscope equipped with an image analyser (Cytovision 3.93.2, Genetix).
Fluorescence in-situ hybridization
The chromosomal location and number of rDNA clusters was determined by fluorescence in-situ hybridization (FISH). Cloned gene fragments of the house mouse 28S rDNA [BE-2-pSP64, 1.5kb; ] and 18S rDNA [SalC-pSP64, 2kb; ] were labelled separately with Digoxigenin-11-dUTP by nick translation according to the Roche Protocol and added to the same hybridization solution. The chromosome slides were treated with RNase for 1 hour at 37°C to remove cellular rRNA, dehydrated in a series of ethanol washes and then air-dried. They were denatured for 2 min at 72°C in 70% formamide, 2XSSC, dehydrated in a series of ice-cold ethanol washes and air-dried. The probes were denatured for 10 min at 72°C. The chromosome slides were hybridized overnight with the two probes (150 ng/slide). The slides were washed at 37°C for 2 min in 2XSSC and 4XT (4x SSC, 0.05% Tween 20, pH = 7) before they were incubated with FITC conjugated with anti-digoxigenin antibody (Roche). The slides were mounted in a Vectashield antifade solution containing DAPI (4', 6'-diamidino-2-phenylindole; Vector Laboratories). A cluster was considered as present on a chromosome pair when a signal was observed on at least one of the homologs. Three positions were possible: i) pericentromeric when the cluster was observed on the proximal region of the chromosome, i.e., adjacent to the centromere, ii) subtelomeric when on the distal end of the chromosome, and iii) interstitial when neither proximal nor distal.
Phylogenetic inference of rDNA cluster evolution
The determination of the phylogenetic distribution of the rDNA clusters throughout the genus Mus first required that the chromosomal orthology among all taxa be established. For this, we relied on the comparative chromosome maps for the subgenera established by Veyrunes et al.  based on Zoo-FISH and the chromosome nomenclature of the house mouse as reference. The Veyrunes et al.  study concluded that the ancestral Mus karyotype comprised 30 syntenic segments (Figure 1). A matrix was then constructed using these 30 segments as characters and the following states: the presence or absence of rDNA clusters on the proximal region of the segment was coded 1 or 0 respectively, and a distal localisation as 2. Apodemus sylvaticus and Rattus rattus were used as outgroups in all analyses. The chromosomal homology between these species and the house mouse was previously determined by Matsubara et al.  and Cavagna et al.  respectively. However, as Matsubara et al.  did not discriminate a distal and proximal segment corresponding to mouse chromosome 1 in A. sylvaticus, the presence or absence of clusters on these segments was coded as "?".
Next, a reference tree was needed providing a topology and branch lengths and including all of the species studied. As none of the published phylogenies met this requirement, we reconstructed a molecular phylogeny of the genus Mus using the nuclear interstitial retinol-binding protein 3 (RBP3 = IRBP) and mitochondrial cytochrome b gene sequences. These data were available in GenBank for all species (Table 2) except Mus cypriacus: the full-length sequence (1140 bp) of this gene was sequenced for M. cypriacus following Montgelard et al.  [EMBL: FR751074]. The sequences were aligned using Bioedit (v.7.0.5). The best fitting model was determined by Modeltest  to be the GTR+G+I model (-LnL = 11707.61; A = 26%, C = 30%, G = 21%, T = 23%, proportion of invariable sites = 62%, gamma shape parameter = 0.95). The phylogenetic analyses were performed on a constrained topology and the combined dataset (the missing data were coded by "?") using maximum likelihood (ML) as implemented in PAUP* v.4b.0. Two nodes of the topology were constrained according to published data: (i) the relationships between the subgenera were unambiguously resolved by Veyrunes et al.  with the subgenus Coelomys diverging first, followed by Nannomys as the sister group of the Mus-Pyromys clade.; (ii) M. cypriacus and M. macedonicus were considered as sister-species as proposed by Cucchi et al. . According to the highest-likelihood topology, the subgenus Mus was divided into three clades: an Southeast Asian group with three species (M. caroli, M. cervicolor and M. cooki), an Indian group (M. terricolor, M. booduga, M. fragilicauda and M. famulus), and a Palearctic group (M. spretus, M. spicilegus, M. macedonicus, M. cypriacus, and the three subspecies of M. musculus [35, 36]). The sister-species relationships within some of the Paleartic taxa were not resolved. However, slight changes in the topology of the reference tree did not affect the results, since the chromosomal distribution of clusters within these groups was quite homogeneous.
We used an ML approach to reconstruct the ancestral states of the rDNA clusters since we had no a priori knowledge on their mode of evolution. Using the previously determined reference tree (topology, branch length), the probability of each state (absent, pericentromeric or distal location) of the rDNA clusters was calculated at all nodes for each orthologous segment. This analysis was performed using R  and the function ACE (ancestral character estimation) from the package APE . We used the ARD model (All-Rates-Different) where all change rates were different among the three rDNA cluster states.
Association between rDNA clusters, breakpoints and centromeres
The number and positions of breakpoints were inferred from the chromosomal phylogeny derived by Veyrunes et al. , and comprised the sites involved in fissions, fusions, translocations and inversions. The events (breakpoints, loss or emergence of a centromere) were then mapped onto the reconstructed rDNA trees for each segment. The localization and co-occurrence of breakpoints and rDNA clusters were then determined.
The co-evolution between centromeres and rDNA clusters was next investigated by BayesTraits  using the previous phylogenetic framework. The correlated evolution between pairs of discrete binary traits (centromere, rDNA cluster) was analysed by ML. A matrix was constructed by coding 1/0 the presence/absence of a centromere and of an rDNA cluster for each orthologous segment and species. This analysis provided the likelihoods of an independent (4 parameters) and a dependent model (8 parameters) of evolution of centromeres and clusters. Likelihood ratio tests were performed to test the significance of the correlated evolution between centromeres and rDNA clusters with the statistics distributed as a chi square (df = 4).
The chromosomal location of the 18S-28S rRNA genes for all species and subspecies is shown in Table 1. The distribution of rDNA clusters in this genus was highly variable ranging from 1 to 21 chromosome pairs. The Asian species harboured the highest number of clusters: 26-42 in Mus caroli, M. cervicolor, M. cooki, M. fragilicauda, M. famulus and M. pahari. The exception was M. platythrix for which only three chromosome clusters were observed. This pattern contrasted markedly with the African subgenus Nannomys in which the number of clusters varied between two and 14.
The rDNA clusters were always located in the pericentromeric region (see insert Figure 1). This organization differed in two species (i) in M. spretus, the rDNA clusters were subtelomeric and (ii) in M. terricolor, the clusters were located between the telomere and the centromere . The clusters were exclusively present on autosomal chromosomes with the exception of M. caroli which showed a signal on the X chromosome.
Eight of the studied taxa showed a distribution of rDNA clusters that was identical to published data [see Table 1; [26, 41, 42]]. This included the rare subtelomeric position noted on chromosome 4 in M. m. domesticus . New rDNA data are provided for six species (Figure 2). rRNA genes were present on all M. caroli chromosomes except the Y chromosome (i.e., all autosomes and the X). Among the Asian species, the rDNA clusters were located in the pericentromeric region of chromosomes 7-19 in M. cooki, of chromosomes 2-19 in M. famulus, and of chromosomes 4-19 in M. cervicolor. In M. fragilicauda, all autosomes harboured clusters. The 18S-28S rRNA genes were located on chromosomes 1, 4, 11, 15-19 in the endemic M. cypriacus. We detected variation in the number of clusters among the three specimens studied from different locations in Cyprus that showed five, seven and eight rDNA clusters respectively.
The 18S-28S rRNA genes were located in the pericentromeric regions of chromosomes: 1, 4, 8, 11, 13, 14 and 17 in M. mattheyi and chromosomes 14, 15, 17 in M. indutus. Only one signal was detected in M. haussa and this was on chromosome 15 which has a diagnostic pericentric inversion . All autosomes in M. musculoides and the South African sample of M. minutoides result from centric fusions of ancestral acrocentric chromosomes . The poor resolution of the centromeric areas in these metacentrics made the precise arm-localization of the hybridization signal problematic. Thus, the clusters were simply assigned to metacentrics: Rb(4.13) and Rb(8.15) in M. musculoides, and Rb(4.7), Rb(12.17), Rb(13.16) and Rb(14.15) in M. minutoides. The clusters were considered as present on both chromosomes of the metacentrics in the phylogenetic analyses.
M. pahari showed 34 signals on the centromeric regions of all but three chromosome pairs (6, 9 and X ).
Our results were similar to those of Matsubara et al. : rDNA clusters were present on three chromosomes 5, 8 and 12.
Inference of the ancestral state of clusters
The results of the rDNA study allowed us to identify a total of 30 clusters in the genus, 27 of which were located in the pericentromeric area and three were telomeric. These clusters were mapped onto the 30 orthologous segments of the ancestral Mus karyotype  (Figure 1). The phylogenetic distribution of the rDNA clusters was established for each of the orthologous segments except one, since an rDNA cluster was never observed on the segment corresponding to chromosome 10d (see additional files 1, 2, 3 and 4). At each node of the tree, the ACE analysis assigned a probability (0-100%) for the absence/presence of a cluster at a pericentromeric/subtelomeric position. From the probability distribution of all segments and nodes, we determined a threshold probability value of 75%, i.e., each state associated with a probability ≥75% was considered as known. Using this threshold value, the state of the rDNA clusters for a segment could be ascertained for 65% of the nodes over all segments i.e., a cluster was absent (or present) in 42% (or 23%) of the nodes. Among the 29 segments, the state of the cluster could not be determined for four of them since all nodal character states had a probability of 50%. In seven additional segments, clusters were never present, whereas for one segment (Chro 18, see Figure 3), a cluster was present at all nodes. Overall, the results indicated that the rDNA clusters showed a high lability: three insertions and 10 deletions could be validated in the deeper nodes and 11 insertions/25 deletions in the terminal branches. The deletions occurred in all subgenera, whereas 86% of the insertions concerned the subgenera Coelomys and Mus. The recurrence of events on segments 6, 7, 12 and 14 (41% of insertions and deletions) is noteworthy. For example, the cluster on the orthologous segment to chromosome 14 was present at the ancestral node of the genus Mus. This cluster disappeared at the nodes leading to the subgenus Nannomys and to the European species of the subgenus Mus (M. spretus, M. musculus, M. spicilegus, M. cypriacus, M. macedonicus), and then reappeared at the node leading to M. cypriacus + M. macedonicus (Figure 3).
Association between breakpoints, rDNA clusters and centromeres
The 42 breakpoints identified by Veyrunes et al.  were mapped onto the ancestral karyotype as well as the location of the rDNA clusters and centromeres observed in the genus (Figure 1). The comparative distribution of these three genomic structures highlighted that 26/30 of the rDNA clusters occurred in the vicinity of breakpoint regions and that 96% of these clusters were also associated with a centromere. Thus, more than half (25/42) of the rearrangement breakpoints in the genus are localized at sites where both an rDNA cluster and a centromere were present. The phylogenetic analyses performed herein together with the published chromosomal phylogeny  provided the opportunity to investigate the temporal dynamics of these genomic structures. In other words, was an rDNA cluster present on the orthologous segment when the rearrangement occurred? When the breakpoints and centromeres were mapped onto the rDNA trees of each orthologous segment, the results showed that the state of clusters for nine breakpoints could not be determined, and that 13 breakpoints did not involve a cluster or a centromere. In all, eight breakpoints occurred in a region where an rDNA cluster as well as a centromere were present. In two cases, the rearrangement resulted in the loss of both the centromere and the cluster suggesting that the break occurred at the distal end of the cluster (see Figure 4). In the six remaining cases, the break was situated between the centromere and the rDNA cluster. These rearrangements led to the loss of the centromere, four of which were subsequently reacquired. Discounting the undetermined clusters, these data suggest that 24% of breakpoints occurred on either side of a cluster.
Given the predominantly pericentromeric location of rDNA clusters and the co-occurrence of centromeres and clusters at breakpoints, the evolutionary association between centromeres and rDNA clusters was tested. The results of this analysis are summarized in Table 3. Of the 30 orthologous segments, one showed no rDNA cluster and one exhibited a cluster at a telomeric location only. These segments were thus eliminated from the analysis. Of the 28 remaining segments, four showed significant evidence of co-evolution between centromeres and rDNA (p-value = 0.03-0.04). For example, the centromere of the segment corresponding to chromosome 19 appeared at the ancestral node of the subgenus Mus. At the same node, an rDNA cluster appeared on this segment and persisted in all species of the subgenus (see Figure 3). In five other cases, co-evolution approached significant values (p = 0.06-0.11). Finally, in the remaining 19 segments, the independent evolution of centromeres and clusters could not be rejected. This is illustrated by the evolution of the rDNA cluster on the segment orthologous to chromosome 7 (see Figure 3). A centromere is always present on this segment but the cluster which was present on the ancestral node to the genus, is lost on the node leading to Nannomys and on the node leading to the European species (Mus spretus, M. musculus, M. spicilegus, M. cypriacus, M. macedonicus).
Extensive variation in number and location of rDNA clusters
This study is the first to report such an impressive variation in number and localization of rDNA clusters within a mammalian genus. The rDNA clusters in the Mus species were present on 1 to 21 pairs of chromosomes with important numerical differences between the taxa of the different subgenera. Overall, the subgenera Coelomys and Mus possessed a much larger number of clusters than species within the subgenera Nannomys and Pyromys. This difference is reflected by that of the diploid numbers since the karyotypes in the former subgenera have more chromosomes (2n = 48, 40) than those in the latter (2n = 18 to 36). As there is generally only one cluster per chromosome pair, a lower diploid number entails fewer chromosomes and thus a reduction in the number of available sites for rRNA genes. Notwithstanding the relation between rDNA cluster and chromosome number, however, variation in the number of clusters was also considerable within the subgenus Mus, although all species share the same 2n = 40 karyotype. In this group, the Palaearctic species have less than half the number of clusters present in the Asian and Indian species. Given that the origin of Mus lies in southern Asia , the decrease in cluster number may be compatible with the loss of rDNA sites by stochastic processes during the westward colonization of the ancestral taxa.
In the genus Mus, the high number of clusters observed raises questions on the processes involved in maintaining such a large number of repeats. Whereas cellular life requires a minimum of one cluster of rRNA genes to construct the ribosome, our results suggest that there may be no upper constraint on the maximal number of clusters in a genome. In a study on rDNA cluster number variation in 40 species of rodents, a mean of 4.2 chromosome pairs carried rDNA clusters (range 1 to 5 ). This is far below the value observed in the present study (mean: 10.1; range 1 to 21). However, the number of clusters provides no information on the number of copies within a cluster, or on their transcriptional activity. It is possible that when many clusters are present, the number of repeat units at each chromosomal site may in fact be small (not studied). With respect to transcriptional activity of repeats within and between clusters, cytogenetic methods involving Ag-staining (silver nitrate stain) may be useful. This approach identifies rRNA genes that were transcribed during the previous interphase. Comparison between published Ag-staining data vs our own FISH analyses for nine Mus species revealed discrepancies in some taxa [41–43]. This was the case for example in M. caroli (15 Ag- staining/20 FISH), M. cervicolor (9/19) and M. fragilicauda (16/19). Differences such as these have similarly been observed in other vertebrate groups such as teleosts, bats, and horses [44–47]. Determining the frequency of these silent rRNA genes on a larger sample of mice will provide clues to their regulation patterns .
Another notable characteristic of the rDNA clusters in the genus Mus is their almost exclusive pericentromeric location since they were subtelomeric in only one species (M. spretus). In rodents, both positions are commonly observed (58% pericentromeric/31% distal vs 11% interstitial) [25, 27, 49–56]. Few mammal species (8/126) exhibit clusters in an interstitial position and in several cases, these result from chromosomal rearrangements [50, 57]. It should be noted that none of the extant Mus species harboured rDNA clusters in an interstitial position, although the phylogenetic reconstruction inferred four instances in which they occupied a transient interstitial position following a rearrangement (Figure 4). The low frequency of these clusters suggests that an interstitial location may not be evolutionary stable and this may be due to two possibilities: (i) the location may be deleterious, or, (ii) as suggested by our results, it may contribute to genomic instability and predispose the chromatin to centromere formation .
The location of clusters between Mus species differed, each taxon showing its own chromosomal distribution. Thus, the chromosomal distribution of rDNA clusters appears as a useful cytogenetic marker (number and position on the chromosome) to discriminate species . However, as our results have shown that rDNA clusters may be labile, their use in inferring orthologous chromosomal sites between species must be treated with caution [e.g. Anura; ].
Do rDNA clusters contribute to rearrangements?
Our study is the first to reconstruct the phylogeny of rDNA clusters. By combining this original approach and a published chromosomal phylogeny, the evolutionary relationship between rDNA clusters and rearrangements could be evaluated and the association with the accelerated rate of genome repatterning of the genus Mus explored. Previous studies have underscored the pivotal role of centromere change in the high rate of chromosomal evolution in this group of rodents [13, 24]. The phylogenetic analyses performed in this present study provided a unique opportunity to reconstruct the timing and infer the position of breakpoints relative to that of the centromere and the rDNA site. Two patterns were evidenced. First, we were able to determine that 24% of the breakpoints occurred at the proximal or distal end of the rDNA clusters. The corresponding rearrangements all involved the coincidental loss of linked centromeres, some of which subsequently re-emerged in the same position. This association is further supported by the significant co-evolutionary behaviour observed between some of the rDNA clusters and centromeres suggesting that the joint presence of these two genomic structures may lead to genome instability and predispose to chromosomal rearrangements. Second, we identified an impressive rate of rDNA cluster change (deletions and insertions), most of which somewhat surprisingly occurred in the subgenus Mus, an evolutionary clade in which all species share the same chromosomal complement. These results suggest that rDNA clusters can move from one chromosome pair to another with no other modifications of the karyotype.
A role of rDNA genes as promoters of genome reorganization, particularly when located in the pericentromeric region, has previously been the focus of studies on the mechanisms of centric or Robertsonian (Rb) fusions [61, 62]. This type of rearrangement involves the joining by the centromere of two non-homologous chromosomes. The rationale of a direct or indirect role in the Rb fusion was based on the fact that rDNA clusters group together during interphase to form one or several nucleoli. The physical proximity between rDNA-bearing chromosomes in interphase was predicted to increase the probability of their being involved in Rb fusions. Strikingly, this rearrangement is prominent within two of the present subgenera studied herein. In two of the samples studied of the subgenus Nannomys, all chromosomes were the product of Rb fusions . The hybridization signal of the rDNA probe seemed to co-localize with the junction between the two chromosomal arms, suggesting that at least one of the breakpoints may have occurred in or close to the rDNA cluster. In the subgenus Mus, populations within two taxa carry Rb fusions: M. m. domesticus and M. terricolor. In the latter, only two rare Rb fusions are documented . In contrast, M. m. domesticus shows a very high diversity of populations carrying different numbers and combinations of fusions, although only five pairs of chromosomes harbour pericentromeric rDNA clusters [, this study]. The fusion mechanism is well documented in M. m. domesticus and has been shown to involve breaks in the centromeric satellite sequences  and not in the rDNA clusters which are conserved intact. In addition, the prediction of a higher involvement of the rDNA-bearing chromosomes in the Rb fusions was not confirmed in the house mouse, as a statistical analysis indicated that they had the lowest frequency of fusion . These findings suggest that centric fusions may involve different sequences in different genomes, and determining the precise mechanism involved requires high resolution sequence analysis.
Is there a functional association between centromeres and rDNA clusters?
Our results highlight two characteristic features of rDNA cluster distribution in the genus Mus: (i) rDNA clusters are preferentially located in pericentromeric regions, and (ii) a change in chromosomal distribution of rDNA clusters always occurs between pericentromeric regions of different chromosomes. This physical linkage between rDNA clusters and centromeres may be driven by their genomic structures, since both are subject to sequence homogenization. In the case of rRNA genes, several studies have underscored the high degree in sequence similarity within and between genomes . Two models of evolution have been proposed to account for this observation. These are the Birth and Death model and the concerted evolution model [66, 67]. By comparing the level of intragenomic variation of rDNA sequences, Ganley and Kobayashi  were able to confirm that rRNA genes evolved via concerted evolution. This model involves two mechanisms of sequence homogenization: gene conversion and non-homologous recombination. Gene conversion corresponds to a non-reciprocal transfer of an allelic difference from one chromosome to its homologue, whereas non-homologous recombination consists in a recombination event between non-sister chromatids. In humans, studies have shown that rDNA clusters were involved in meiotic rearrangements at a frequency >10% per cluster and per meiosis , and that the regulatory and coding sequences were highly homogenized . Similar observations of high sequence homogeneity exist for both subtelomeric and particularly centromeric regions in several mammalian species [71–73].
How and when would concerted evolution take place? During early prophase, all the chromosomes migrate into one area of the nucleus and adopt a particular orientation known as the bouquet in which all telomeres attach to the nuclear membrane (Figure 5) . Double-strand breaks also appear at this stage that are programmed to be repaired by recombination be it by reciprocal exchange (cross-over) or gene conversion . Thus, the formation of the bouquet provides the physical opportunity for concerted evolution to occur between similar sequences on non-homologous chromosomes (such as centromeres, and the associated rDNA genes [71, 72, 76]). Reciprocally, the presence of clusters could increase the rate of non-homologous recombination in the pericentromeric or subtelomeric regions. Several consequences of rDNA evolution are expected. First, rDNA clusters would be predicted to localize preferentially in centromeric or subtelomeric regions. This is in agreement with available data in mammalian species. Second, whereas the telomere orientation in the bouquet obviously leads to the close proximity between centromeres on acrocentric chromosomes during meiosis, this is less evident where metacentric chromosomes are involved (Figure 5). This difference is likely to influence the rate of homogenization between centromeres (and linked rRNA genes) of acrocentric vs metacentric chromosomes. Such a pattern has in fact been confirmed in humans and recently in the pig in which the percentage similarity of centromeric sequences is higher in acrocentric than metacentric chromosomes [71, 77]. Finally, the mechanism of concerted evolution in a bouquet context paves the way for the occurrence of exchanges of rDNA clusters between non-homologous chromosomes.
In conclusion, our results agree with a main outcome of recent studies indicating that breakpoints overlap with a diversity of repetitive families among different groups , which thus may include rDNA clusters. The present analysis on the dynamics of rDNA clusters within the genus Mus, has revealed an impressive variation in the number and location of clusters between taxa, and highlighted the strong evolutionary relationship between rDNA clusters and centromeres. Both of these genomic structures coincide with breakpoints in the genus Mus, suggesting that the accumulation in the centromeric region of a large number of repeats subject to concerted evolution may contribute to the high rates of chromosome repatterning observed in this group. However, the high rate of rDNA change observed in the chromosomally invariant subgenus Mus indicates that the presence of these sequences is insufficient to lead to genome instability. Thus, an additional factor would be required to trigger evolutionary plasticity. Emerging clues now point to modifications of the epigenetic state of DNA, particularly of interspersed repeats, as a prime source of instability [1, 6, 78]. In plants and insects, it has been shown that transposable elements may be associated with rDNA clusters, and thus, may be responsible for changes in their chromosomal distribution [15, 79]. However, no direct association between rRNA genes and transposable elements has so far been described in mammals. Although this study deals with large-scale rearrangements and low-resolution breakpoints, our results may serve to focus genomic investigations on the factors triggering genome plasticity and evolution.
short interspersed nuclear element
long interspersed nuclear element
long terminal repeat
fluorescence in situ hybridization
Brown JD, O'Neill RJ: Chromosomes, conflict, and epigenetics: chromosomes speciation revisited. Annual Review Genomics Human Genetics. 2010, 11: 13.11-13.26.
Murphy WJ, Larkin DM, der Wind A, Bourque G, Tesler G, Auvil L, Beever JE, Chowdhary BP, Galibert F, Gatzke L, et al: Dynamics of mammalian chromosome evolution inferred from multispecies comparative maps. Science. 2005, 309: 613-617. 10.1126/science.1111387.
Pevzner P, Tesler G: Genome rearrangements in mammalian evolution: lessons from human and mouse genomes. Genome Research. 2003, 13: 37-45. 10.1101/gr.757503.
Eichler EE, Sankoff D: Structural dynamics of eukaryotic chromosome evolution. Science. 2003, 301: 793-797. 10.1126/science.1086132.
Kehrer-Sawatzki H, Cooper DN: Molecular mechanisms of chromosomal rearrangement during primate evolution. Chromosome Research. 2008, 16: 41-56. 10.1007/s10577-007-1207-1.
Carbone L, Harris RA, Vessere GM, Mootnick AR, Humphray S, Rogers J, Kim SK, Wall JD, Martin D, Jurka J, et al: Evolutionary breakpoints in the gibbon suggest association between cytosine methylation and karyotype evolution. PLoS Genet. 2009, 5: e1000538-10.1371/journal.pgen.1000538.
Froenicke L, Lyons AL: Hotspots of mammalian chromosome evolution. 2008, Chichester
Glover TW: Common fragile sites. Cancer Letters. 2006, 232: 4-12. 10.1016/j.canlet.2005.08.032.
Bailey JA, Yavor AM, Massa HF, Trask BJ, Eichler EE: Segmental duplications: organization and impact within the current human genome project assembly. Genome Research. 2001, 11: 1005-1017. 10.1101/gr.GR-1871R.
Stankiewicz P, Lupski JR: Genome architecture, rearrangements and genomic disorders. Trends in Genetics. 2002, 18: 74-81. 10.1016/S0168-9525(02)02592-1.
O'Neill RJ, Eldridge BMD, Metcalfe CJ: Centromere dynamics and chromosome evolution in marsupials. Journal of Heredity. 2004, 95: 375-381. 10.1093/jhered/esh063.
Robinson TJ, Ruiz-Herrera A, Froenicke L: Dissecting the mammalian genome - new insights into chromosomal evolution. Trends in Genetics. 2006, 22: 297-301. 10.1016/j.tig.2006.04.002.
Mlynarski EE, Obergfell CJ, O'Neill MJ, O'Neill RJ: Divergent patterns of breakpoint reuse in Muroid rodents. Mammalian Genome. 2010, 21: 77-87. 10.1007/s00335-009-9242-1.
Kominami R, Mishima Y, Urano Y, Sakai M, Muramatsu M: Cloning and determination of the transcription site of ribosomal RNA gene of the mouse. Nucleic Acids Research. 1982, 10: 1963-1979. 10.1093/nar/10.6.1963.
Eickbush TH, Eickbush DG: Finely orchestrated movements: evolution of the ribosomal RNA genes. Genetics. 2007, 175: 477-485. 10.1534/genetics.107.071399.
Baker RJ, Maltbie M, Owen JG, Hamilton MJ, Bradley RD: Reduced number of ribosomal sites in bats: evidence for a mechanism to contain genome size. Journal of Mammalogy. 1992, 73: 847-858. 10.2307/1382206.
Averbeck KT, Eickbush TH: Monitoring the mode and tempo of concerted evolution in the Drosophila melanogaster rDNA locus. Genetics. 2005, 171: 1837-1846. 10.1534/genetics.105.047670.
Datson PM, Murray BG: Ribsomal DNA locus evolution in Nemesia: transposition rather than structural rearrangement as the key mechanism?. Chromosome Research. 2006, 14: 845-857. 10.1007/s10577-006-1092-z.
Huang J, Ma L, Fei S, Li L: 45S rDNA regions are chromosome fragile sites expressed as gaps in vitro on metaphase chromosomes of root-tip meristematic celles in Lolium spp. Plos One. 2008, 3: e2167-10.1371/journal.pone.0002167.
Trifonov V, Stanyon R, Nesterenko A, Fu B, Perelman P, O'Brien P, Stone G, Rubtsova N, Houck M, Robinson T, et al: Multidirectional cross-species painting illuminates the history of karyotypic evolution in Perissodactyla. Chromosome Research. 2008, 16: 89-107. 10.1007/s10577-007-1201-7.
Musser GG, Carleton MD: Family Muridae. 2005, Washington: Smithsonian Institution
Veyrunes F, Catalan J, Sicard B, Robinson TJ, Duplantier JM, Granjon L, Dobigny G, Britton-Davidian J: Autosome and sex chromosome diversity among the African pygmy mice, subgenus Nannomys (Murinae; Mus). Chromosome Research. 2004, 12: 369-382.
Shimada T, Aplin KP, Suzuki H: Mus lepidoides (Muridae, Rodentia) of central burma is a distinct species of potentially great evolutionary and biogeographic significance. Zoological Science. 2010, 27: 449-459. 10.2108/zsj.27.449.
Veyrunes F, Dobigny G, Yang F, O'Brien PCM, Catalan J, Robinson TJ, Britton-Davidian J: Phylogenomics of the genus Mus (Rodentia; Muridae): extensive genome repatterning is not restricted to the house mouse. Proceeding of the Royal Society London B. 2006, 273: 2925-2934. 10.1098/rspb.2006.3670.
Cavagna P, Stone G, Stanyon R: Black rat (Rattus rattus) genomic variability characterized by chromosome painting. Mammalian Genome. 2002, 13: 157-163.
Matsubara K, Nishida-Umehara C, Kuroiwa A, Tsuchiya K, Matsuda Y: Identification of chromosome rearrangements between the laboratory mouse (Mus musculus) and the Indian Spiny mouse (Mus platythrix) by comparative FISH analysis. Chromosome Research. 2003, 11: 57-64. 10.1023/A:1022010116287.
Matsubara K, Nishida-Umehara C, Tsuchiya K, Nukaya D, Matsuda Y: Karyotypic evolution of Apodemus (Muridae, Rodentia) inferred from comparative FISH analyses. Chromosome Research. 2004, 12: 383-395.
Lee MR, Elder FFB: Yeast simulation of bone marrow mitosis for cytogenetic investigations. Cytogenetic Cellular Genetics. 1980, 26: 36-40. 10.1159/000131419.
Cowell JK: A photographic representation of the variability of G-banded structure of the chromosomes of the mouse karyotype. Chromosoma. 1984, 89: 294-320. 10.1007/BF00292478.
Hassouna N, Michot B, Bachellerie JP: The complete nucleotide sequence of mouse 28S rRNA gene. Implications for the process of size increase of the large sub-unit rRNA in higher eukaryotes. Nucleic Acids Research. 1984, 12: 3563-3585. 10.1093/nar/12.8.3563.
Raynal F, Michot B, Bachellerie JP: Complete nucleotide sequence of mouse 18S rRNA gene: comparison with other available analogs. FEBS Letters. 1984, 167: 263-268. 10.1016/0014-5793(84)80139-8.
Montgelard C, Bentz S, Tirard C, Verneau O, Catzeflis FM: Molecular systematics of sciurognathi (Rodentia):the mitochondrial cytochrome b and 12S rRNA genes support the Anomaluroidea (Pedetidae and Anomaluridae). Molecular Phylogenetics and Evolution. 2002, 22: 220-233. 10.1006/mpev.2001.1056.
Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-817. 10.1093/bioinformatics/14.9.817.
Cucchi T, Orth A, Auffray JC, Renaud S, Fabre L, Catalan J, Hadjisterkotis E, Bonhomme F, Vigne JD: A new endemic species of the subgenus Mus (Rodentia, Mammalia) on the island of Cyprus. Zootaxa. 2006, 1241: 1-36.
Suzuki H, Shimada T, Terashima M, Tsuchiya K, Aplin K: Temporal, spatial, and ecological modes of evolution of Eurasian Mus based on mitochondrial and nuclear gene sequences. Molecular Phylogenetics and Evolution. 2004, 33: 626-646. 10.1016/j.ympev.2004.08.003.
Chevret P, Veyrunes F, Britton-Davidian J: Molecular phylogeny of the genus Mus (Rodentia: Murinae) based on mitochondrial and nuclear data. Biological Journal of the Linnean Society. 2005, 84: 417-427. 10.1111/j.1095-8312.2005.00444.x.
The R Development Core Team: R: A Language and Environment for Statistical Computing. 2008, Version 2.8.0
Paradis E, Claude J, Strimmer K: APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics. 2004, 20: 289-290. 10.1093/bioinformatics/btg412.
Pagel M, Meade A: Bayesian analysis of correlated evolution of discrete charcaters by reversible-jump Markov chain monte Carlo. The American Naturalist. 2006, 167: 808-825. 10.1086/503444.
Bardhan A, Sharma T: Dynamic of nucleolus organizer regions and karyotype evolution in Indian pygmy field mice. Cytogenetics and Cell Genetics. 2000, 91: 47-51. 10.1159/000056817.
Winking H, Nielsen K, Gropp A: Variable positions of NORs in Mus musculus. Cytogenetics Cell Genetics. 1980, 26: 158-164. 10.1159/000131437.
Suzuki H, Kurihara Y, Kanehisa T, Moriwaki T: Variation in the distribution of silver-staining Nucleolar Organizer Regions on the chromosomes of the wild mouse, Mus musculus. Molecular Biology and Evolution. 1990, 7: 271-282.
Elsevier SM, Ruddle FH: Location of genes coding for 18s and 28s ribosomal RNA within the genome of Mus musculus. Chromosoma. 1975, 52: 219-228. 10.1007/BF00332112.
Derjusheva SE, Loginova J, Parada R, Chiryaeva O, Smirnov A, Jaszczak K: The comparative analysis of NOR polymorphism detected by FISH and Ag-staining on horse chromosomes. Caryologia. 1998, 51: 1-11.
Ferro DAD, Neo DM, Moreira O, Bertollo LAC: Nucleolar organizing regions, 18S and 5S rDNA in Astyanax scabripinnis (Pisces, Characidae): populations distribution and functional diversity. Genetica. 2000, 110: 55-62. 10.1023/A:1017963217795.
Santos N, Fagundes V, Yonenaga-Yassuda Y, De Souza M: Localization of rDNA genes in Phyllostomidae bats reveals silent NORs in Artibeus cinereus. Hereditas. 2002, 136: 137-143. 10.1034/j.1601-5223.2002.1360208.x.
Singh M, Kumar R, Nagpure NS, Kushwaha B, Gond I, Lakra WS: Chromosomal localization of the 18S and 5S rDNA using FISH in the genus Tor (Pisces, Cyprinidae). Genetica. 2009, 137: 245-252. 10.1007/s10709-009-9367-x.
Ide S, Miyazaki T, Maki H, Kobayashi T: Abundance of ribosomal RNA gene copies maintains genome integrity. Science. 2010, 327: 693-696. 10.1126/science.1179044.
Djalali M, Hameister H, Vogel W: Further chromosomal studies on Ellobius lutescens: Heteromorphism of chromosome No. 1 is not associated with sex determination. Experientia. 1986, 42: 1281-1282. 10.1007/BF01946423.
Britton-Davidian J, Catalan J, Granjon L, Duplantier JM: Chromosomal phylogeny and evolution in the genus Mastomys (Mammalia, Rodentia). Journal of Mammalogy. 1995, 76: 248-262. 10.2307/1382333.
Stitou S, Burgos M, Zurita F, Jimenez R, Sanchez A, Diaz de la Guardia R: Recent evolution of NOR-bearing and sex chromosomes of the North African rodent Lemniscomys barbarus. Chromosome Research. 1997, 5: 481-485. 10.1023/A:1018473131446.
Fagundes V, Christoff AU, Amaro-Ghilard RC, Scheibler DR, Yonenaga-Yassuda Y: Multiple interstitial ribosomal sites (NORs) in the brazilian squirrel Sciurus aestuans ingrami (Rodentia, Sciuridae) with 2n = 40. Genetics and Molecular Biology. 2003, 26: 253-257. 10.1590/S1415-47572003000300007.
Zambelli A, Catanesi CI, Vidal-Rioja L: Autosomal rearrangements in Graomys griseoflavus (Rodentia): a model of non-random divergence. Hereditas. 2003, 139: 167-173.
Dobigny G, Ozouf-Costaz C, Bonillo C, Volobouev V: Evolution of rRNA gene clusters and telomeric repaets during explosive genome repatterning in Taterillus X (Rodentia, Gerbillinae). Cytogenetic and Genome Research. 2003, 103: 94-103. 10.1159/000076296.
Arslan A, Albayrak I, Pamukoglu N, Yorulmaz T: Nucleolar organizer regions (NORs) of the spiny mouse Acomys cilicicus (Mammalia: Rodentia) in Turkey. Turk Journal Zoology. 2008, 32: 75-78.
Gornung E, Castiglia R, Rovatsos MT, Marchal JA, diaz De La Guardia R, Sanchez A: Comparative cytogenetic study of two sister species of Iberian ground voles, Microtus (Terricola) duodecimcostatus and M. (T). lusitanicus (Rodentia, Cricetidae). Cytogenetic and Genome Research. 2011, 132: 144-150. 10.1159/000321572.
Ropiquet A, Hassanin A, Pagacova E, Gerbault-Seureau M, Cernohorska H, Kubickova S, Bonillo C, Rubes J, Robinson T: A paradox revealed: karyotype evolution in the four-horned antelope occurs by tandem fusion (Mammalia, Bovidae, Tetracerus quadricornis). Chromosome Research. 2010, 18: 277-286. 10.1007/s10577-010-9115-1.
Marshall OJ, Chueh AC, Wong LH, Choo AKH: Neocentromeres: new insights into centromere structure, disease development, and karyotype evolution. The American Journal of Human Genetics. 2008, 82: 261-282. 10.1016/j.ajhg.2007.11.009.
Nguyen TT, Aniskin VM, Gerbault-Seureau M, Planton H, Renard JP, Nguyen BX, Hassanin A, V.T V: Phylogenetic position of the saola (Pseudoryx nghetinhensis) inferred from cytogenetic analysis of elenen species of Bovidae. Cytogenetic and Genome Research. 2008, 122: 41-54. 10.1159/000151315.
Carvalho KA, Garcia PCA, Recco-Pimentel SM: NOR dispersion, telomeric sequence detection in centromeric region and meiotic multivalent configurations in species of the Aplastodiscus albofrenatus group (Anura, Hyldae). Cytogenetic and Genome Research. 2009, 126: 359-367. 10.1159/000264179.
Winking H, Dulic B, Bulfield G: Robertsonian karyotype variation in the European house mouse, Mus musculus; survey of present knowledge and new observations. Zeitschrift für Säugetierkunde. 1988, 53: 148-161.
Sullivan BA, Jenkins LS, Karson EM, Leana-Cox J, Schwartz S: Evidence for structural heterogeneity from molecular cytogenetic analysis of dicentric robertsonian translocations. American Journal Human Genetics. 1996, 59: 167-175.
Pialek J, Hauffe HC, Searle JB: Chromosomal variation in the house mouse. Biological Journal of the Linnean Society. 2005, 84: 535-563. 10.1111/j.1095-8312.2005.00454.x.
Garagna S, Marziliano N, Zuccotti M, Searle JB, Capanna E, Redi CA: Pericentromeric organization at the fusion point of mouse Robertsonian translocation chromosomes. Proceedings of the National Academy of Sciences of the United States of America. 2001, 98: 171-175. 10.1073/pnas.98.1.171.
Gazave E, Catalan J, Da Graça Ramalhinho M, Da Luz Mathias M, Nunes AC, Dumas D, Britton-Davidian J, Auffray JC: The non-random occurrence of Robertsonian fusion in the house mouse. Genetical Research. 2003, 81: 33-42. 10.1017/S001667230200602X.
Nei M, Rooney AP: Concerted and birth-and-death evolution of multigene families. Annual Review Genetics. 2005, 39: 121-152. 10.1146/annurev.genet.39.073003.112240.
Arnheim N, Krystal M, Schmickel R, Wilson G, Ryder O, Zimmer E: Molecular evidence for genetic exchanges among ribosomal genes on nonhomologous chromosomes in man and apes. Proceedings of the National Academy of Sciences. 1980, 77: 7323-7327. 10.1073/pnas.77.12.7323.
Ganley ARD, Kobayashi T: Highly efficient concerted evolution in the ribosomal DNA repeats: Total rDNA repeat variation revealed by whole-genome shotgun sequence data. Genome Research. 2007, 17: 184-191. 10.1101/gr.5457707.
Stults DM, Killen MW, Pierce HH, Pierce AJ: Genomic architecture and inheritance of human ribosomal DNA. Genome Research. 2008, 18: 13-18.
Gonzalez IL, Sylvester JE: Human rDNA: evolutionary patterns within the genes and tandem arrays derived from multiple chromosomes. Genomics. 2001, 73: 255-263. 10.1006/geno.2001.6540.
Hirai H, Taguchi T, Godwin AK: Genomic differentiation of 18S ribosomal DNA and b-satellite DNA in the hominoid and its evolutionary aspects. Chromosome Research. 1999, 7: 531-540. 10.1023/A:1009237412155.
Linardopoulou EV, Williams EM, Fan Y, Friedman C, Young JM, Trask BJ: Human subtelomeres are hot spots of interchromosomal recombination and segmental duplication. Nature. 2005, 437: 94-100. 10.1038/nature04029.
Gauthier P, Hima K, Dobigny G: Robertsonian fusions, pericentromeric repeat organization and evolution: a case study within a highly polymorphic rodent species, Gerbillus nigeriae. Chromosome Research. 2010, 18: 473-486. 10.1007/s10577-010-9128-9.
Scherthan H, Weich S, Schwegler H, Heyting C, Härle M, Cremer T: Centromere and telomere movements during early meiotic prophase of mouse and man are associated with the onset of chromosome pairing. The Journal of Cell Biology. 1996, 134: 1109-1125. 10.1083/jcb.134.5.1109.
Roeder GS: Meiotic chromosomes: it takes two to tango. Genes & Development. 1997, 11: 2600-2621. 10.1101/gad.11.20.2600.
Brannan CI, Disteche CM, Park LS, Copeland NG, Jenkins NA: Autosomal telomere exchange results in the the rapid amplification and dispersion and Csf2ra genes in wild-derived mice. Mammalian Genome. 2001, 12: 882-886. 10.1007/s00335-001-2084-0.
Adega F, Chaves R, Guedes-Pinto H: Suiformes orthologous orthogous satellite DNAs as a hallmark of Pecari tajacu and Tayassu pecari (Tayassuidae) evolutionary rearrangements. Micron. 2008, 39: 1281-1287. 10.1016/j.micron.2008.03.004.
O'Neill RJW, O'Neill MJ, Graves JAM: Undermethylation associated with retroelement activation and chromosome remodeling in an interspecific mammalian hybrid. Nature. 1998, 393: 68-72. 10.1038/29985.
Raskina Q, Barber JC, Nevo E, Belyayev A: Repetitive DNA and chromosomal rearrangements: speciation-related events in plant genomes. Cytogenetic and Genome Research. 2008, 120: 351-10.1159/000121084.
Robins JH, McLenachan PA, Phillips MJ, Craig L, Ross HA, Matisoo-Smith E: Dating of divergences within the Rattus genus phylogeny using whole mitochondrial genomes. Molecular Phylogenetic Evolution. 2008, 49: 460-466. 10.1016/j.ympev.2008.08.001.
Suzuki H, Tsuchiya K, Takezaki N: A molecular phylogenetic framework for the ryukyu endemic rodents Tokudaia osimensis and Diplothrix legata. Molecular Phylogenetics and Evolution. 2000, 15: 15-24. 10.1006/mpev.1999.0732.
Serizawa K, Tsuchiya K: A phylogenetic view on species radiation in Apodemus inferred from variation of nuclear and mitochondrial genes. Biochemical Genetics. 2000, 38: 27-39. 10.1023/A:1001828203201.
Lundrigan BL, Jansa SA, Tucker PK: Phylogenetic relationship in the genus Mus, based on paternally maternally and biparentally inherited characters. Systematic Biology. 2002, 51: 410-431. 10.1080/10635150290069878.
Veyrunes F, Britton-Davidian J, Robinson TJ, Calvet E, Denys C, Chevret P: Molecular phylogeny of the African pygmy mice, subgenus Nannomys (Rodentia, Murinae, Mus): implications for chromosomal evolution. Molecular Phylogenetics and Evolution. 2005, 36: 358-369. 10.1016/j.ympev.2005.02.011.
Martin Y, Gerlach G, Schlotterer C, Meyer A: Molecular phylogeny of European muroid rodents based on complete cytochrome b sequences. Molecular Phylogenetic Evolution. 2000, 16: 37-47. 10.1006/mpev.1999.0760.
Bibb MJ, Van Etten RA, Wright CT, Walberg MW, Clayton DA: Sequence and gene organization of mouse mitochondrial DNA. Cell. 1981, 26: 167-180. 10.1016/0092-8674(81)90300-7.
Acknowledgements and funding
We are grateful to A. Orth for animal samples and to J.-P. Bachellerie for providing the rDNA probes. We would like to thank P. Boursot and T. Robinson for their critical comments of earlier versions of this manuscript. This study was supported by recurrent funding from the CNRS and UM2. This is publication ISEM n° 2011-21.
BC carried out the laboratory work, and participated in all analyses. JC participated in the acquisition of data, and drew the figures. FV participated in the design of the study. EJPD participated in the analyses. JBD conceived the study. BC and JBD wrote the paper. All authors read and approved the final manuscript.