Integration of molecular cytogenetics, dated molecular phylogeny, and model-based predictions to understand the extreme chromosome reorganization in the Neotropical genus Tonatia (Chiroptera: Phyllostomidae)
BMC Evolutionary Biology volume 15, Article number: 220 (2015)
Defining factors that contributed to the fixation of a high number of underdominant chromosomal rearrangements is a complex task because not only molecular mechanisms must be considered, but also the uniqueness of natural history attributes of each taxon. Ideally, detailed investigation of the chromosome architecture of an organism and related groups, placed within a phylogenetic context, is required. We used multiple approaches to investigate the dynamics of chromosomal evolution in lineages of bats with considerable karyotypic variation, focusing on the different facets contributing to fixation of the exceptional chromosomal changes in Tonatia saurophila. Integration of empirical data with proposed models of chromosome evolution was performed to understand the probable conditions for Tonatia’s karyotypic evolution.
The trajectory of reorganization of chromosome blocks since the common ancestor of Glossophaginae and Phyllostominae subfamilies suggests that multiple tandem fusions, as well as disruption and fusions of conserved phyllostomid chromosomes were major drivers of karyotypic reshuffling in Tonatia. Considerable variation in the rates of chromosomal evolution between phyllostomid lineages was observed. Thirty–nine unique fusions and fission events reached fixation in Tonatia over a short period of time, followed by ~12 million years of chromosomal stasis. Physical mapping of repetitive DNA revealed an unusual accumulation of LINE-1 sequences on centromeric regions, probably associated with the chromosomal dynamics of this genus.
Multiple rearrangements have reached fixation in a wave-like fashion in phyllostomid bats. Different biological features of Tonatia support distinct models of rearrangement fixation, and it is unlikely that the fixations were a result of solely stochastic processes in small ancient populations. Increased recombination rates were probably facilitated by expansion of repetitive DNA, reinforced by aspects of taxon reproduction and ecology.
How chromosomal rearrangements become established in natural populations has been a recurrent question in evolutionary cytogenetics for decades [1–3]. The absence of “the one model” of chromosomal fixation can be attributed to the amount of variables associated with cellular, molecular, ecological, and population distinctiveness of each taxonomic group where chromosome rearrangements have become established . Because individual models usually focus on specific aspects of chromosomal change, each containing its own set of variables and assumptions, none is widely accountable for all patterns observed in nature. The widespread model predicts that chromosome rearrangements resulting in underdominance require small population sizes to achieve fixation [1, 5]. However, some observed patterns are difficult to reconcile with this single model, such as cases of karyotypic megaevolution, where there is a major re-patterning of the chromosomal segments due to a large number of multiple types of chromosomal rearrangements . In order to understand karyotypic evolution in a given taxon, it is necessary to contrast the observed data with all available models, and then define what combination of circumstances most likely explains the observations. It is noteworthy that some models of chromosomal speciation intrinsically bear assumptions applicable to models of rearrangement fixation within a population (see  for a review).
Because of their high levels of karyotypic variation, phyllostomid bats are an ideal group to investigate the mechanisms responsible for chromosome evolution and stasis. This group arose in the Eocene, and has diversified into 11 distinct subfamilies, with over 57 genera and 200 species [8, 9]. The varying rates of karyotypic change predicted for morphologically divergent lineages, coupled with their ample set of biological differences, depict a scenario in which highly rearranged karyotypes might have been established by different means in species with distinct life histories. The subfamily Phyllostominae is an example of karyotypic variation among closely related species, with diploid numbers (2n) ranging from 16 to 34 . Most of the species within this group present a karyotype of 2n = 32; however, the two species in the genus Tonatia (T. bidens and T. saurophila) have a 2n = 16 karyotype that is heavily rearranged relative to other members of the subfamily, including the closely related genus, Lophostoma [11, 12]. Because no homology could be detected, based on G-band comparison, to the karyotypes of closely related genera, Tonatia was proposed as an example of a megaevolved karyotype .
Molecular cytogenetic techniques, such as cross-species chromosome painting, allow the detection of chromosome homologies based on conservation at DNA sequence level, and thus are robust tools to estimate the magnitude of chromosomal change, especially in groups for which classical G-banding comparison have failed to detect interspecific chromosomal homologies. Moreover, by integrating chromosome painting data, molecular dating, and an explicit phylogenetic framework, one can infer evolutionary trajectories and estimate rates of fixation of chromosome rearrangements within particular lineages . The resulting information, together with the mapping of other DNA sequences on the karyotype of closely related taxa, is valuable to test the predictions of different proposed models for fixation of chromosome rearrangements and provides a starting point to unravel the underlying forces shaping karyotypes of extant species.
The goal of this study was to investigate the radical reshuffling of the Tonatia karyotype. To do so, we 1) used cross-species chromosome painting to map the chromosomal homologies between T. saurophila and six phyllostomid species from four subfamilies (Macrotinae, Phyllostominae, Glossophaginae, and Lonchophyllinae); 2) reconstructed phylogenetic relationships using molecular data to define ancestral and derived syntenic associations on the karyotypes of these subfamilies; and, 3) used molecular time estimates to provide a temporal framework for the observed chromosomal changes. Our data allowed us to estimate the minimum number of rearrangements required to derive the extant karyotype of T. saurophila from the proposed ancestral Phyllostominae condition, as well as the number of unique fusion and fission events in each lineage since the divergence of two distinct subfamilies (Glossophaginae/Phyllostominae). In addition, probes of repetitive sequences were mapped onto the genomes of Phyllostominae bats to investigate their role as potential drivers of chromosome reshuffling within the group. The probability of fixation of Tonatia chromosome rearrangements is discussed in light of different models, integrating data on fluorescence in situ hybridization (FISH), long-term effective population sizes (Ne), selection coefficients for specific rearrangements, and natural history aspects of the genus.
Chromosome preparations and G-banding
Chromosome preparations were obtained from bone marrow or tissue culture following the methods described by Baker et al. , and in accordance with animal welfare guidelines established by the Texas Tech University Animal Care and Use Committee. Voucher specimens, tissues, and cell suspensions are deposited at the Natural Science Research Laboratory (NSRL) of the Museum of Texas Tech University, under the following identification numbers: Tonatia saurophila (TSA) from Honduras, Ecuador, and Costa Rica: TK101463♀, TK104519♂, TK104616♀, TK104655♀, TK167809♀; Lophostoma occidentalis (LOC) from Ecuador: TK104520♂ and TK104505♀; Mimon crenulatum (MCR) from Ecuador: TK104437♀, TK104615♀, TK104516♂, TK104620♀, TK135711♂, TK135714♀, TK135715♀. G-banding technique was carried out following Seabright , with trypsin (0.25 %), incubation times varying from 15 to 30 min at 37 °C.
Whole chromosome paint probes from Macrotus californicus (MCA), generated by DOP-PCR of flow-sorted chromosomes , were used for the chromosome painting experiments on metaphase chromosomes of T. saurophila, L. occidentalis, and M. crenulatum. In situ hybridizations were carried out according to Yang et al.  and Volleth et al. . Cy3-conjugated streptavidin (Amersham Biosciences, 1:1000) was used for detection after 72 h of hybridization. The chromosomes were counterstained with DAPI (4’,6-diamidino-2-phenylindole). To distinguish MCA chromosomes 8, 10, and 13, we have used Myotis myotis chromosome 12 specific paint, which corresponds to MCA 8. No hybridization with MCA Y chromosome paint was performed. Chromosome painting results for Glossophaginae and Lonchophyllinae species from Sotero-Caio et al.  were integrated in the analysis.
FISH with repetitive DNA
Our repetitive DNA FISH (fluorescence in situ hybridizations) included the use of 45S ribosomal DNA (rDNA), long interspersed element 1 (LINE-1) and telomeric (TTAGGG)n, sequences as probes. The rDNA sequences were amplified from genomic DNA of Noctilio albiventris, using the following18S and 28S specific primers: 28S F-1: 5’ GCC GAA ACG ATC TCA ACC TAT T 3’; 28S R-1: 5’ GAG CCA ATC CTT ATC CCG AA 3’; 18S F-1: 5’TCA ACT TTC GAT GGT AGT 3’; and 18S R-1: 5’GCA AGC TTA TGA CCC GCA CTT A 3’. The LINE-1 probe used (Tbid8b) was isolated from T. saurophila genome by degenerate PCR amplification of a 575 bp portion of ORF2 straddling the reverse transcriptase domain, followed by cloning and enrichment for LINE-1 fragments retaining a single open reading frame, i.e., the youngest elements in the T. saurophila genome [19, 20]. The rDNA PCR products and the isolated LINE-1 clones were labeled with biotin using the Bionick™ DNA Labeling System (Molecular Probes®) or with digoxigenin using the DIG-Nick Translation Mix (Roche). Hybridizations were performed following the FISH protocol described by Raudsepp and Chowdhary . Commercially available ready-to-use human chromosome Pan-telomeric probes (Star*FISH©, Cambio) were used according to manufacturer’s instructions.
We performed Bayesian and Maximum Likelihood (ML) phylogenetic analyses using available sequences from the mitochondrial cytochrome-b and 12S–16S rRNA fragment, as well as a portion of the nuclear RAG2 gene, using the Hoffmann et al.  dataset, adding representatives of Glossophaga to the ingroup, and representatives of Noctilio and Pteronotus, as outgroup sequences. Due to missing data for L. occidentalis, four other species of Lophostoma were used to estimate the relationships between this genus and other phyllostomine (Additional file 1: Table S1) Sequences for each fragment were aligned using MUSCLE , and concatenated prior to phylogenetic estimation. Bayesian estimation was done in MrBayes version 3.1.2 , running four simultaneous chains for 5 × 107 generations, sampling trees every 2,500 generations, and using default priors. Convergence was accessed by measuring the standard deviation of the split frequency among parallel chains. Chains were considered converged once the average split frequency was lower than 0.01. We recovered a majority-rule consensus of the last 2,500 trees collected after convergence was reached, and discarded trees collected before convergence. We implemented a 7-partition analysis where each codon position in the nuclear RAG2, the mitochondrial cytochrome-b, and the rDNA fragment had an independent GTR + Γ model of nucleotide substitution. The best-fitting models of nucleotide substitution for partition were independently selected using the corrected Akaike Information Criterion [25, 26]. Maximum-likelihood searches were done in Treefinder version June 2007 , following similar methods used in Bayesian analyses. Node support was estimated with 1,000 bootstrap replicates .
Molecular time estimation was performed in BEAST v2.1.3  according to Khan et al.  to estimate divergences within Noctilio. A lognormal molecular clock with an independent model of nucleotide substitution for each separate partition was implemented. Because the node corresponding to the Mormoopidae/Phyllostomidae split is dated to the late Oligocene, it was constrained to have a minimum age of 28.3 million years (myr) with an exponential mean of 3.2, allowing the maximum age of this divergence to occur ~40.6 myr ago [31, 32]. Node dates were estimated using a birth–death process prior as proposed for reconstructing phylogenies without fossil lineages . Analyses consisted of multiple independent runs for a total of 150 × 106 iterations, with every 1,000th iteration logged for all analyses. Independent runs were combined for each data set using LogCombiner version 2.1.3 . TRACER version 1.6  was used to determine appropriate burn-in (10 %), and examine convergence, effective sample sizes (ESSs), and 95 % highest probability density intervals (HPD) of constrained priors. Due to missing data, a shorter alignment from the RAG2 dataset was used to estimate divergence time since the split of the two Anoura from the Glossophaginae common ancestor.
Determination of ancestral karyotypes
We integrated the chromosome painting data of the present work and Sotero-Caio et al.  to determine the syntenic associations present in the karyotype of the Phyllostominae and Glossophaginae common ancestor (PGA), as well as to determine syntenic associations or chromosomal blocks that were present at the basal node of each subfamily. Using M. californicus as outgroup, syntenic associations present in species from the two subfamilies were parsimoniously placed in the PGA karyotype. Chromosome morphology and banding patterns were evaluated to uncover synapomorphic inversions.
Rates of chromosomal evolution
The rates of fixation of unique chromosome rearrangements for Phyllostominae and Glossophaginae species were calculated using the estimated number of fixed unique chromosomal changes divided by the amount of time required to derive the karyotypes in each particular lineage. The magnitude of fixed unique chromosomal changes per species was obtained using two major parameters: 1) number of unique fusions, calculated as the number of fusion events from the ancestral subfamilial state required to form unique syntenic associations in the species analyzed; and, 2) number of required fissions to derive the number of chromosome blocks on extant karyotypes, calculated using the syntenic associations of M. californicus and the number of blocks on the proposed ancestral subfamilial karyotype as the outgroups. We decided not to include other karyotypic change parameters in this analysis, such as inversions or centromeric shifts, due to the subjectivity of banding pattern comparison resulting from the differential G-banding quality of different preparations. The time component was calculated using internodal divergence times obtained from BEAST analysis. Standard errors for the time estimates were calculated from the upper and lower confidence intervals of each node, allowing the inference of the upper and lower rates of fixation of chromosome rearrangements.
Models of fixation of chromosome rearrangements and long-term effective population size estimations
The model of fixation of chromosome rearrangements from Lande  was used to estimate average long-term effective population size (Ne) from chromosomal data for TSA and also Anoura cultrata (ACU) for comparisons. We have used three selection coefficient values (s = 0.1, 0.3, and 0.5), for differential selective pressures/disadvantage on heterozygous conditions to account for a combination of different types of rearrangements presented by these species. A spontaneous chromosomal mutation rate (u) range of 10−3 to 10−4 was selected to accommodate for the variation proposed for mammals ( and references therein). Rates of fixation of chromosomal change (R) to be incorporated in this method were estimated using the size-corrected average generation length data estimated by Pacifici et al.  to calculate R in terms of fixation of rearrangements per species per generation. Since estimations of divergence times from molecular data are broad, we included upper and lower values of our date estimates to account for the obtained standard errors and find corresponding Ne variations. We contrasted the obtained results with the predictions of Chesser and Baker , which test the feasibility of Lande’s model, taking bat biological characteristics into consideration. Additional theoretical models were examined in light of the current knowledge of TSA biological features.
Hybridizations and chromosome characterization of Tonatia saurophila (TSA)
The karyotype of Tonatia saurophila from Ecuador and Costa Rica had diploid and fundamental numbers of 2n = 16, FN = 22, respectively (Fig. 1). This chromosomal formula differs from the reported karyotype of 2n = 16, FN = 20 of specimens from Trinidad and Brazil by a pericentric inversion on the fourth largest autosome [10, 12, 38]. The hybridizations with MCA whole chromosome paints detected a total of 36 pairs of homologous segments on the autosomes of TSA. With the exception of the allosomes, no chromosome of MCA was conserved as an individual chromosome in the TSA genome. Additionally, the syntenies of most MCA chromosomes appeared disrupted on the TSA karyotype, except for seven chromosome pairs, which seemed to have fused to other chromosomes as entire syntenic blocks (MCA 9, 11, 14, 15, 17, 18, and 19). The assignment of MCA homologous segments on the G-banded karyotype of TSA is shown in Fig. 1a.
The hybridizations with repetitive DNA probes (Fig. 1b–d) revealed different patterns: 1) the major cluster of 45S rDNA is located at the distal portion of the short arm of TSA chromosome 4; 2) telomeric DNA was detected distally on all chromosome arms; and, 3) the probe of young Tonatia LINE-1 sequences hybridized to the centromeres of all TSA chromosomes, to the negative G-band regions of some autosomes, apparently corresponding to breakpoint regions between different MCA syntenic blocks, as well as to the entire length of the sex chromosomes. Additionally, we found LINE-1 signal in a region on TSA 3, between MCA 5 and 16, whereby no MCA or PHA  homology have been detected.
Hybridizations and chromosome characterization of Lophostoma occidentalis (LOC)
The karyotype of Lophostoma occidentalis is comprised of 17 chromosome pairs (2n = 34) with FN = 62, as described in Baker et al.,  see Velazco and Cadenillas  for taxonomic status and species diagnosis. This karyotype is identical to the one reported for the species L. silvicolum . A total of 23 pairs of homologous segments were detected by FISH with the MCA autosome probes (Fig. 2a). The rDNA FISH has revealed nucleolar organizing regions (NORs) on two LOC chromosome pairs: on the distal portion of LOC 13q and LOC 16p. Telomeric sequences were located at all chromosome termini, whereas hybridizations with Tonatia LINE-1 have not revealed any differential accumulation of these sequences in any particular chromosome or chromosomal region (Fig. 2b–d).
Hybridizations and chromosome characterization of Mimon crenulatum (MCR)
The specimens of Mimon crenulatum analyzed herein had 2n = 32, with either FN = 59 or 60. The individual with a FN = 59 was a heterozygote for two chromosomal inversions, encompassing one of the chromosomes of pair five (acrocentric vs. submetacentric) and one of pair 6 (subtelocentric vs. submetacentric) (Fig. 3). The individuals with FN = 60 had a submetacentric morphology for both MCR 5 and 6 (data not shown). A single pair of NORs was detected distally on the long arm of MCR 14 (Fig. 3b) whereas the telomeric FISH detected the telomeres of all chromosomes, plus one non-telomeric site on the centromeric region of MCR 7 (Fig. 3c). Hybridizations with the young LINE Tbid8b did not reveal consistent accumulation patterns for these sequences (Fig. 3d). Representative images of the chromosome painting on the three species are presented on (Additional file 2: Figure S1).
Phylogenetic relationships and divergence dates of Phyllostominae
The relationships of the studied taxa revealed by the Bayesian analysis are similar to those reported by Hoffmann et al. . Within Phyllostominae, the nodes for tribes Phyllostomini, Macrophyllini and Vampyrini (sensu Baker et al. ) were recovered. Within Phyllostomini, Tonatia diverged first, followed by the radiation of lineages that gave rise to four genera, including Mimon and Lophostoma. The tree topology with the estimated divergence dates and lineage name abbreviations are presented in the (Additional file 3: Figure S2).
Phyllostominae, Glossophaginae, and PGA ancestral karyotypes
Several conserved chromosomes and chromosomal blocks were identified between the analyzed species (Fig. 4, Additional file 4: Table S2). Seven chromosomes (including a submetacentric X chromosome) were shared by all phyllostomids analyzed and were present at the most basal family node, as well as the PGA karyotype. The syntenic blocks corresponding to MCA 16 through 19 were prone to homoplasy, but were able to define their ancestral states using MCA as the outgroup. The different states of MCA 16, 17, 18, and 19 on the karyotypes of the analyzed species are presented in a phylogenetic context in Fig. 4. Seven additional syntenic segments were synapomorphic for Phyllostominae, Glossophaginae and Lonchophyllinae. The above information allows an estimation of a PGA karyotype comprised of 36 chromosomes.
The comparative analysis among the Phyllostominae (TSA, LOC, and MCR) has shown considerable karyotypic conservation between LOC and MCR. Their karyotypes share most of the syntenic blocks identified by chromosome painting, except for the differential linkage of MCA 16 and 19. Interestingly, the karyotype of Glossophaga soricina (GSO, Glossophaginae) also presents most of the syntenies retained by LOC and MCR. Because of the observed conservatism, we propose that the Phyllostominae and Glossophaginae ancestral karyotypes were essentially those of LOC and GSO, respectively (synapomorphies for each subfamily are shown in Fig. 4). Three Phyllostominae synapomorphies, inversion (i) of MCA 4, 12 + 2qi, and 17 + 18) were confirmed by G-band comparison with the karyotype of Trachops cirrhosus (Phyllostominae) as a Phyllostomini outgroup .
Rates of chromosomal evolution and long-term effective population size estimations: chromosome evolution of TSA (Phyllostominae) and ACU (Glossophaginae)
TSA and ACU had the highest values of unique fusions and fissions (39 and 14, respectively) required to derive their karyotypes from their respective ancestral subfamilial karyotypes. We estimate that at least 13 autapomorphic tandem fusions have contributed to the reorganization of TSA karyotype. The magnitude of unique fixed changes for MCR consisted of a single fission event, not including the two polymorphic inversions, whereas no changes from the ancestral subfamilial chromosome complement were observed for GSO and LOC (Table 1, Additional file 4: Tables S3–S5). The rates of fixation of chromosome rearrangements varied for the different branches across the phylogenetic tree, with periods of chromosomal stasis versus evolution (Table 1). The standard errors for the ages of distinct branches ranged from 1.2 to 3.4 myr and contributed to wide variations in R (fixation events over time) values, especially for lineages where chromosome stasis occurred. The fixation of a large number of chromosome rearrangements in the lineage that gave rise to the genus Tonatia occurred in a relatively short period of time (6 ± 2.7 myr), followed by chromosome stasis for the remaining 12 ± 1.7 myr after divergence of the two species within the genus. Similarly, a short period of time of 4.6 ± 3.2 myr was estimated for the fixation of 14 unique changes in the lineage that gave rise to the genus Anoura from the glossophagine ancestor. The two cases described above differ from the general trend observed for most of the branches analyzed, in which less than two fixation events occurred for periods of time spanning up to 12 myr.
The results of the Ne analyses for TSA and ACU are summarized in Fig. 5, and the detailed calculations are presented in the Additional file 5: Table S6. As expected from Lande’s model , long-term effective population sizes required for the fixation of the number of rearrangements presented by the two species are low (less than 70 and 80 individuals for TSA and ACU, respectively), even when the lower and upper values of fixation times are considered. For selection coefficients greater than 0.3, Ne values are significantly reduced to lower than 24 for TSA and 27 for ACU, considering the highest rates of mutation for the range used.
Extensive and conservative chromosomal evolution in Phyllostominae and Glossophaginae
Our results show that the rates of chromosomal evolution since the divergence of the lineages that gave rise to the analyzed species were not constant over time, and that the set of rearrangements leading to the chromosomal constitution of genera such as Tonatia (Phyllostominae) and Anoura (Glossophaginae), have occurred as the result of waves of rearrangements specific to those lineages. On the other hand, little chromosomal evolution was observed for the lineages leading to the genera Mimon and Lophostoma (Phyllostominae), and Glossophaga (Glossophaginae), since these groups have retained most of the syntenic blocks that were present in the PGA karyotype. G-banding and FISH analyses suggest that the karyotype of other phyllostomine species have not changed considerably from the ancestral Phyllostominae condition proposed herein [6, 12, 43, 44], and that the ancestral Phyllostomidae karyotype is not very different from that presented by the phyllostomine Phyllostomus hastatus in terms of syntenic associations [12, 45], implying an overall conservative trend in terms of chromosomal rearrangements within the subfamily. Similarly, Glossophaga is nested within a clade comprised of four other genera (Erophylla, Brachyphylla, Monophyllus, and Leptonycteris), all of which with identical karyotypes . The chromosomal variation within Glossophaginae seems to be restricted to a monophyletic clade comprising Anoura (older origin) and four other relatively recent genera: Hylonycteris, Choeroniscus, Choeronycteris, and Musonycteris [6, 8, 47]. These observations suggest that the processes responsible for the extensive chromosome reshuffling in Tonatia and Anoura might be the result of ephemeral molecular processes, occurring independently in a temporally punctuated fashion.
The unusual high rates of gross chromosomal change of Tonatia are comparable to the extreme rates of chromosomal evolution within gibbon genera (22 rearrangements/5myr in Hoolock) and can be placed among one of the highest among mammals [48, 49]. For Tonatia, a total of 39 unique fission and fusion events have been estimated in the present analysis, however the definition of specific rearrangements as well as determination of the order of fixation events for this lineage proved to be a difficult task. Thus, the number of rearrangements of Tonatia and Anoura calculated herein might be an overestimation if translocations were an important component of their karyotypic mutation framework (each accounting for a simultaneous fission + fusion). Nevertheless, we estimated that 13 tandem fusions have occurred to derive the karyotype of Tonatia from the Phyllostominae ancestor, which would still result in at least 26 rearrangements in that lineage if the remainder rearrangements were all translocation-derived. In addition, refined molecular cytogenetic techniques, such as BAC-FISH, microarrays, and selective sequencing, have shown that results from cross-species chromosome painting might underestimate the total number of rearrangements identified for a given species [48, 50]. Thus, since our analysis does not account for inversions, centromeric shifts, and smaller chromosome translocations, we conclude that the Tonatia lineage had one of the highest rates of chromosomal evolution among mammals.
An intriguing aspect of Tonatia karyotypic evolution is that regardless of the mechanism responsible for generating these high rates of change, it seems to have ceased after the establishment of the 2n = 16 karyotype and before speciation of the genus. That raises the question of how long natural populations can cope with high rates of change without extinction due to the negative effects on meiosis and before selection drives more intensive control of the underlying mechanisms. Perhaps, an even more relevant question would be how a large number of rearrangements like that observed for Tonatia could have achieved fixation within a relatively short period of time. Other bats of the family Phyllostomidae have also undergone multiple rearrangements, resulting in karyotypes with little identifiable G-band arm homologies . These include several non-related lineages from the subfamilies Carolliinae, Glossophaginae, and Stenodermatinae. Family-wide cross-species chromosome painting studies will help unravel lineage-specific rates of rearrangements and their implications for speciation and diversification among phyllostomids. Questions as to why closely related groups can present differential levels of rearrangements will then be able to be addressed.
Finally, we observed that some syntenic blocks are prone to convergence within phyllostomid bats. Examples are the multiple independent fusion/fission events, resulting in slightly different associations in MCR and LOC (Phyllostominae) versus LCO (Lonchophyllinae) and ACU (Glossophaginae). Conversely, the convergent states of MCA 18 might be derived from intense chromosome reorganization on ACU and TSA and less likely are an intrinsic feature of this chromosome.
Repetitive DNA and evolution of genome architecture in Tonatia and other phyllostomids
Repetitive DNA (and duplicated gene families) has been hypothesized to be one of the major drivers of chromosomal change due to their potential to promote non-homologous chromosome exchange through illegitimate recombination between homologous sequences [51–55]. Our chromosome painting and rates of chromosomal change analyses indicated that the patterns of rapid reshuffling of TSA karyotype might be the outcome of non-allelic recombination of such sequences. To better understand the contribution of repetitive sequences in the genomic architecture of TSA, we analyzed the patterns of distribution and accumulation of three classes of repetitive DNA: telomeric, rDNA, and young LINE-1 sequences.
Interstitial telomeric sites (ITS) have been regarded as remnants of chromosomal fusion promoted by repetitive subtelomeric elements [56, 57] as well as inversions [58, 59]. Our telomeric FISH revealed a telomere-only pattern, despite the 13 tandem fusions that contributed to the generation of the highly rearranged karyotype of TSA. This suggests that these sequences might either be absent or in low copy numbers at the identified tandem fusion regions. If subtelomeric DNA was responsible for the rapid reshuffling of Tonatia genome, the significant number reduction of these and of telomeric sequences after fusion events might be the critical factor resulting in karyotypic stasis for the last 12 myr. Interestingly, MCR presented telomeric sequences enrichment on the centromere of a chromosome pair that is conserved among several phyllostomine species (MCR 7/ LOC 7). This is a submetacentric element, derived from a metacentric ancestral phyllostomid chromosome (MCA 4). ITS amplification in this particular MCR chromosome is most likely correlated with other types of species-specific repetitive DNA, rolling circle replication, or mechanisms that allow tandem amplification of the canonical telomere sequence within satellite DNA, rather than be derived from gross chromosome rearrangements, such as inversions . Alternatively, the relocation of telomeres by at least two inversions of the ancestral metacentric element might have provided raw material for ITS amplification in MCR that was lost in other species of Phyllostominae.
LINE elements are ubiquitous transposable elements (TEs) in mammalian genomes, and our hybridization results have shown a substantial centromeric buildup of a recently transposed sequence corresponding to partial ORF2 of a Tonatia-specific LINE-1 , which is not shared by the other two phyllostomine bats analyzed. This LINE-1 centromeric enrichment is uncommon among mammals, for which there is a preferential LINE-1 accumulation on positive G-band regions throughout autosomal length, and especially in sex chromosomes. The centromeric regions, however, are usually devoid of these sequences. The only in situ hybridization study of LINE-1 distribution on bat chromosomes has shown that at least for one phyllostomid species, Carollia brevicauda, the distribution of LINE-1 sequences follows this widespread pattern [52, 61–65]. The pattern seen in Tonatia, however, resembles the LINE-1-Alu-SVA (LAVA) gibbon-specific TE massive centromeric accumulation detailed by Carbone et al. . Carbone et al.  have demonstrated for the first time that high rates of gross chromosome rearrangements can be a direct result of premature termination of transcripts of chromosome segregation genes due to insertion of these composite elements.
Accumulation of at least a portion of LINE-1 at TSA centromeric and some breakpoint regions provide a starting point to investigate their potential role as karyotype modifiers in this species. On one hand, it is possible that Tonatia LINE-1 sequences have played a role in promoting chromosome translocations through non-homologous recombination. Alternatively, because the probes used correspond solely to a portion of LINE-1, we cannot reject the hypothesis that these sequences are part of composite elements, which could increase levels of rearrangements by affecting expression patterns of chromosome segregation genes . Regardless of their role as karyotypic modifiers, an important aspect of centromeric LINE-1 accumulation observed in TSA chromosomes is their contribution as structural centromeric components and possibly epigenetic regulators of the centromere dynamics in TSA. If part of the centromeric satellite DNA, these sequences are expected to be stable structural components of heterochromatin, serving as binding sites for chromatin remodeling complexes and centromere-related binding factors and having specific amplification dynamics . Centromeric amplification of retroelement-derived repeats has been reported for other mammalian groups with varying degrees of rates of chromosomal evolution, which might display distinct levels of epigenetic repression of TEs [51, 68, 69]. There is no evidence that the probes used in our study differ substantially in structure and composition from sequences isolated from other bats [19, 20] or that phyllostomid bats present an unusual repetitive DNA landscape compared to other mammals . Thus, our results provide an insight into the complexity of TE roles for the chromosomal dynamics within TSA.
Ribosomal genes (rDNA) are unusual repetitive sequences when it comes to patterns of accumulation and dispersal (reviewed by McStay and Grummt ). Besides particularities of rDNA evolution, such as concerted evolution and paralog recombination , rDNA sites seem to be hotspots for double strand breaks, which are often in physical proximity on non-homologous chromosomes . Because of the abovementioned characteristics, rDNA-bearing chromosomes are proposed to be highly prone to translocations. Our results and the integration of previously published work on FISH with ribosomal DNA and on chromosome painting shows that the distributional pattern of rDNA in phyllostomid bats follows the trends described by Milhomem et al. , where lack of association of the major rDNA cluster to specific homologous chromosomes has been reported. We observed the association of NORs to MCA 19 in distantly-related species (LOC and GSO), to MCA 5 in the closely related LOC and MCR, to MCA 8 in TSA, and to MCA 9 in ACU. Similarly, published work using Phyllostomus hastatus (PHA) whole chromosome probes have shown rDNA association to different PHA chromosomes, namely PHA 9 and 12 in Desmodus rotundus, PHA 15 in PHA and possibly Diphylla ecaudata [44, 75, 76], PHA 2 and 13 in Micronycteris hirsuta , and PHA 5 and X in Carollia brevicauda [43, 78]. This trend indicates an independent amplification of rDNA in non-homologous chromosomes in different taxa, and does not support the hypothesis that the rearrangements observed for phyllostomid species were derived from translocations induced by rDNA recombination. Thus, these sequences might not be informative phylogenetic markers on the study of Phyllostomidae evolution.
Models for the fixation of chromosome rearrangements and their implications for chromosomal evolution in TSA
The integration of our data with a set of published models and the biological features of the studied species allows further investigation of the mechanisms responsible for the fixation of large numbers of chromosomal rearrangements in phyllostomid bats. We first used the classic model of Lande  of stochastic fixation of rearrangements to estimate the long-term effective population sizes required to fix the chromosomal rearrangements of TSA and ACU. To our knowledge, this is the first attempt to integrate chromosome painting data and molecular dating to this model. The most important assumption of Lande’s model is the lower fitness conferred to heterozygotes for chromosomal rearrangements, with tandem fusions having a higher associated negative selective pressure than Robertsonian rearrangements [1, 2]. One prediction from this model is that small population sizes would facilitate fixation of chromosomal rearrangements that have a significant negative heterotic effect (drift). Accordingly, TSA and ACU analysis output consisted of small Ne values regardless of the selective pressures for the set of rearrangements, underlying chromosomal mutation rates, or estimated fixation rates. The larger values of long-term population sizes, approximately 70 and 80 individuals for TSA and ACU, respectively, would have to be accompanied by large mutation rates and less than 10 % production of unviable gametes by heterozygotes to allow the fixation of the detected rearrangements. The Ne values obtained from this model are unrealistic to ensure long-term persistence of a species. Therefore, either rates of chromosome fixation are not appropriate predictors of Ne, or these species are exceptions with successful population recovery after extended periods of bottlenecks for over 6 myrs. We propose that because our results have shown that the rates of chromosomal rearrangements are not constant over time, they violate one of the major assumptions of Lande’s model (homogeneous rate of fixation of chromosome rearrangements), and thus, this model is a poor predictor of historical population sizes for phyllostomid bats.
According to the model of Chesser and Baker , the fixation of chromosomal rearrangements in small populations through drift would only be achieved under the conditions of extremely small population size (less than 10 founders), little underdominance of the rearrangements on heterozygotes, and a high number of offspring per mating. Our results present conflictive situations regarding the two first circumstances because the smaller population size values obtained here are associated with the strongest selection coefficient against the rearrangements (s > 0.5). Thus, if Chesser and Baker’s predictions are accurate, no fixation is possible for TSA for most of our Ne values, since population sizes smaller than 20 founders were associated with fecundity reductions greater than 0.25. Additionally, with a single offspring per mating , population recovery would not be possible for small demes of less than 10 TSA individuals.
Interestingly, different biological features of Tonatia are compatible with different models. For instance, TSA is considered locally rare in terms of abundance, and has only been found in groups of less than 11 individuals with high fidelity to their roosting sites [80, 81], which agrees with predictions of rearrangement fixation through drift in small demes. The potential for daily dispersal of these bats, however, is relatively high, which might indicate that local populations are actually larger due to area coverage . Additionally, some Tonatia extant features are compatible to the predictions of the adaptive rearrangement models, in which suppressed recombination might lead to adapted chromosomal forms [83, 84]. These include the omnivorous feeding habit, and current patterns of distribution for the two species which might be widely distributed and adapted to different types of climatic and ecological conditions as a result of phenotypic plasticity provided by ancestral rearrangements.
The simulations of fixation of chromosome rearrangements by Chesser and Baker  were performed considering a single initial chromosome rearrangement, and its trajectory of fixation/population size. In the cases of rapid and extensive chromosome reshuffling, new computational models allowing for an input of several concomitant and subsequent rearrangements using fixation rates established with the methods used in this work would be useful to understand for how long polymorphisms for different rearrangements can coexist.
In summary, it is unlikely that the great number and complexity of TSA chromosomal rearrangements would be fixed solely through drift in natural populations. Thus, we propose alternative scenarios might have played a role to circumvent the negative heterosis in these species. 1) Based on the observations of LINE-1 accumulation for this species, increased recombination [52, 85] among repetitive sequences in this lineage might have led to unprecedented mutation rates. The larger amplification in gene-poor regions, such as centromeres and perhaps telomeres might have led to recurrent breaks involving the same chromosomes in different individuals. This would allow an increased occurrence of homozygote individuals for specific rearrangements, facilitating fixation in large populations. 2) The current adaptive success of Tonatia species might be indicative of adaptive karyotypes. We hypothesize that the large number of disrupted euchromatic regions, coupled with suppressed recombination between locally adapted loci in different rearrangement regions would have promoted adaptation and speciation in Tonatia. 3) We cannot reject the hypothesis that other factors such as meiotic drive, a phenomenon poorly studied in bats, might have acted synergistically to allow the fixation of Tonatia rearrangements. Therefore, we conclude that outcomes of the different models discussed above bring features that act synergistically with ancestral biological features, facilitating the fixation of multiple rearrangements in TSA.
Degenerate oligonucleotide primed polymerase chain reaction
Fluorescence in situ hybridization
Fundamental number of autosomal arms
Interstitial telomeric site
LINE-1-Alu-SVA transposable element
Long interspersed element 1
Lophostoma + Mimon ancestor
- Ne :
Effective population size
Nucleolar organizing region
Phyllostominae + Glossophaginae ancestor
Rate of fixation of chromosome rearrangements
Spontaneous chromosome mutation rate
Wright S. On the probability of fixation of reciprocal translocations. Amer Nat. 1941;75(761):513–22.
White MJD. Animal cytology and evolution. 3rd ed. London: Cambridge University Press; 1973.
Rieseberg LH. Chromosomal rearrangements and speciation. Trends Ecol Evol. 2001;16(7):351–8.
King M. Species evolution: the role of chromosome change. New York: Cambridge University Press; 1993.
Lande R. Effective deme sizes during long-term evolution estimated from rates of chromosomal rearrangement. Evolution. 1979;33(1):234–51.
Baker RJ, Bickham JW. Karyotypic evolution in bats: evidence of extensive and conservative chromosomal evolution in closely related taxa. Syst Zool. 1980;29(3):239–53.
Faria R, Navarro A. Chromosomal speciation revisited: rearranging theory with pieces of evidence. Trends Ecol Evol. 2010;25(11):660–9.
Baker RJ, Bininda-Emonds ORP, Mantilla-Meluk H, Porter CA, Van Den Bussche RA. Molecular timescale of diversification of feeding strategy and morphology in new world leaf-nosed bats (Phyllostomidae): a phylogenetic perspective. In: Gunnell GF, Simmons NB, editors. Evolutionary history of bats: fossils, molecules and morphology. Cambridge University Press: Cambridge Studies in Molecules and Morphology–New Evolutionary Paradigms; 2012. p. 385–409.
Solari S, Martínez-Arias V. Cambios recientes en la sistemática y taxonomía de murciélagos Neotropicales (Mammalia: Chiroptera). Therya. 2014;5(1):167–96.
Baker RJ. Karyology. In: Baker RJ, Jones JK, Carter DC, editors. Biology of bats of the New World family Phyllostomatidae Part III. Special Publications, The Museum, Texas Tech University. Lubbock; 1979. p. 107–55.
Lee Jr TE, Hoofer SR, Van Den Bussche RA. Molecular phylogenetics and taxonomic revision of the genus Tonatia (Chiroptera: Phyllostomidae). J Mammal. 2002;83(1):49–57.
Ribas T, Rodrigues L, Nagamachi C, Gomes A, Rissino J, O’Brien P, et al. Phylogenetic reconstruction by cross-species chromosome painting and G-banding in four species of Phyllostomini tribe (Chiroptera, Phyllostomidae) in the Brazilian Amazon: an independent evidence for monophyly. PLoS ONE. 2015;10(3), e0122845.
Burt DW, Bruley C, Dunn IC, Jones CT, Ramage A, Law AS, et al. The dynamics of chromosome evolution in birds and mammals. Nature. 1999;402(6760):411–3.
Baker RJ, Hamilton MJ, Parish DA. Preparations of mammalian karyotypes under field conditions. Occas Papers Mus Texas Tech Univ. 2003;228:1–8.
Seabright M. A rapid banding technique for human chromosomes. Lancet. 1971;298(7731):971–2.
Sotero-Caio CG, Volleth M, Gollahon LS, Fu B, Cheng W, Ng BL, et al. Chromosomal evolution among leaf-nosed nectarivorous bats-evidence from cross-species chromosome painting (Phyllostomidae, Chiroptera). BMC Evol Biol. 2013;13(1):276.
Yang F, Carter NP, Shi L, Ferguson-Smith MA. A comparative study of karyotypes of muntjacs by chromosome painting. Chromosoma. 1995;103(9):642–52.
Volleth M, Heller K-G, Pfeiffer RA, Hameister H. A comparative ZOO-FISH analysis in bats elucidates the phylogenetic relationships between Megachiroptera and five microchiropteran families. Chromosome Res. 2002;10(6):477–97.
Cantrell MA, Grahn RA, Scott L, Wichman HA. Isolation of markers from recently transposed LINE-1 retrotransposons. Biotechniques. 2000;29(6):1310–7.
Cantrell MA, Scott L, Brown CJ, Martinez AR, Wichman HA. Loss of LINE-1 activity in the megabats. Genetics. 2008;178(1):393–404.
Raudsepp T, Chowdhary BP. FISH for mapping single copy genes. In: Murphy WJ, editor. Phylogenomics. New York: Humana Press; 2008. p. 31–49.
Hoffmann FG, Hoofer SR, Baker RJ. Molecular dating of the diversification of Phyllostominae bats based on nuclear and mitochondrial DNA sequences. Mol Phylogenet Evol. 2008;49(2):653–8.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32(5):1792–7.
Ronquist F, Huelsenbeck J. MRBAYES: Bayesian inference of phylogeny. Bioinformatics. 2003;19(12):1572–4.
Akaike H. A new look at the statistical model identification. IEEE Trans Autom Control. 1974;19(6):716–23.
Sugiura N. Further analysts of the data by Akaike’s Information Criterion and the finite corrections: Further analysts of the data by Akaike’s. Communications in Statistics-Theory and Methods. 1978;7(1):13–26.
Jobb G, Von Haeseler A, Strimmer K. TREEFINDER: a powerful graphical analysis environment for molecular phylogenetics. BMC Evol Biol. 2004;4(1):18.
Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39:783–91.
Bouckaert R, Heled J, Kühnert D, Vaughan T, Wu C-H, Xie D, et al. BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput Biol. 2014;10(4), e1003537.
Khan FAA, Phillips CD, Baker RJ. Timeframes of speciation, reticulation, and hybridization in the bulldog bat explained through phylogenetic analyses of all genetic transmission elements. Syst Biol. 2013;63(1):96–110.
Morgan G, Czaplewski N. New bats in the Neotropical families Emballonuridae and Mormoopidae from the Oligocene and Miocene of Florida, and the biochronology of Florida Whitneyan, Arikareean, and Hemingfordian faunas. J Vertebr Paleontol. 2002;22:90A.
Meredith RW, Janečka JE, Gatesy J, Ryder OA, Fisher CA, Teeling EC, et al. Impacts of the Cretaceous Terrestrial Revolution and KPg extinction on mammal diversification. Science. 2011;334(6055):521–4.
Gernhard T. The conditioned reconstructed process. J Theor Biol. 2008;253(4):769–78.
Drummond A, Rambaut A: LogCombiner v1. 4.8. In. http://beast.bio.ed.ac.uk/LogCombiner; 2007.
Rambaut A, Drummond A: Tracer v1. 4. In. http://tree.bio.ed.ac.uk/software/tracer/. 2007.
Pacifici M, Santini L, Di Marco M, Baisero D, Francucci L, Marasini GG, et al. Generation length for mammals. Nat Conserv. 2013;5:89–94.
Chesser RK, Baker RJ. On factors affecting the fixation of chromosomal rearrangements and neutral genes: computer simulations. Evolution. 1986;40(3):625–32.
Santos N, Fagundes V, Yonenaga‐Yassuda Y, Souza MJ. Localization of rRNA genes in Phyllostomidae bats reveals silent NORs in Artibeus cinereus. Hereditas. 2002;136(2):137–43.
Baker RJ, Fonseca RM, Parish DA, Phillips CA, Hoffmann FG. New bat of the genus Lophostoma (Phyllostomidae: Phyllostominae) from northwestern Ecuador. Occas Papers Mus Texas Tech Univ. 2004;232:1–16.
Velazco PM, Cadenillas R. On the identity of Lophostoma silvicolum occidentalis (Davis & Carter, 1978)(Chiroptera: Phyllostomidae). Zootaxa. 2011;2962:1–20.
Baker RJ, Hoofer SR, Porter CA, Van Den Bussche RA. Diversification among New World leaf-nosed bats: An evolutionary hypothesis and classification inferred from digenomic congruence of DNA sequence. Occas Papers Mus Texas Tech Univ. 2003;230:1–32.
Barros HMDR, Sotero-Caio CG, Santos N, Souza MJ. Comparative cytogenetic analysis between Lonchorhina aurita and Trachops cirrhosus (Chiroptera, Phyllostomidae). Genet Mol Biol. 2009;32(4):748–52.
Pieczarka JC, Nagamachi CY, O’Brien PC, Yang F, Rens W, Barros RMS, et al. Reciprocal chromosome painting between two South American bats: Carollia brevicauda and Phyllostomus hastatus (Phyllostomidae, Chiroptera). Chromosome Res. 2005;13(4):339–47.
Gomes AJB, Nagamachi CY, Rodrigues LRR, Farias SG, Rissino JD, Pieczarka JC. Karyotypic variation in Rhinophylla pumilio Peters, 1865 and comparative analysis with representatives of two subfamilies of Phyllostomidae (Chiroptera). Comparative Cytogenetics. 2012;6(2):213.
Pieczarka JC, Gomes AJ, Nagamachi CY, Rocha DC, Rissino JD, O’Brien PC, et al. A phylogenetic analysis using multidirectional chromosome painting of three species (Uroderma magnirostrum, U. bilobatum and Artibeus obscurus) of subfamily Stenodermatinae (Chiroptera-Phyllostomidae). Chromosome Res. 2013;21(4):383–92.
Baker RJ, Bass RA. Evolutionary relationship of the Brachyphyllinae to the glossophagine genera Glossophaga and Monophyllus. J Mammal. 1979;60(2):364–72.
Ribeiro NAB, Nagamachi CY, Pieczarka JC, Rissino JD, Neves ACB, Gonçalves ACO, et al. Cytogenetic analysis in species of the Subfamily Glossophaginae (Phyllostomidae, Chiroptera) supports a polyphyletic origin. Caryologia. 2003;56(1):85–96.
Capozzi O, Carbone L, Stanyon RR, Marra A, Yang F, Whelan CW, et al. A comprehensive molecular cytogenetic analysis of chromosome rearrangements in gibbons. Genome Res. 2012;22(12):2520–8.
Carbone L, Harris RA, Gnerre S, Veeramah KR, Lorente-Galdos B, Huddleston J, et al. Gibbon genome and the fast karyotype evolution of small apes. Nature. 2014;513(7517):195–201.
Roberto R, Capozzi O, Wilson RK, Mardis ER, Lomiento M, Tuzun E, et al. Molecular refinement of gibbon genome rearrangements. Genome Res. 2007;17(2):249–57.
Wichman H, Payne C, Ryder O, Hamilton M, Maltbie M, Baker R. Genomic distribution of heterochromatic sequences in equids: implications to rapid chromosomal evolution. J Hered. 1991;82(5):369–77.
Wichman HA, Van Den Bussche RA, Hamilton MJ, Baker RJ. Transposable elements and the evolution of genome organization in mammals. Genetica. 1992;86:287–93.
Eichler EE, Sankoff D. Structural dynamics of eukaryotic chromosome evolution. Science. 2003;301(5634):793–7.
Argueso JL, Westmoreland J, Mieczkowski PA, Gawel M, Petes TD, Resnick MA. Double-strand breaks associated with repetitive DNA can reshape the genome. Proc Natl Acad Sci. 2008;105(33):11845–50.
Anand RP, Tsaponina O, Greenwell PW, Lee C-S, Du W, Petes TD, et al. Chromosome rearrangements via template switching between diverged repeated sequences. Genes Dev. 2014;28(21):2394–406.
Meyne J, Baker RJ, Hobart HH, Hsu TC, Ryder OA, Ward OG, et al. Distribution of non-telomeric sites of the (TTAGGG) n telomeric sequence in vertebrate chromosomes. Chromosoma. 1990;99(1):3–10.
Lee C, Sasi R, Lin C. Interstitial localization of telomeric DNA sequences in the Indian muntjac chromosomes: further evidence for tandem chromosome fusions in the karyotypic evolution of the Asian muntjacs. Cytogenet Genome Res. 1993;63(3):156–9.
Paço A, Chaves R, Vieira-da-Silva A, Adega F. The involvement of repetitive sequences in the remodelling of karyotypes: The Phodopus genomes (Rodentia, Cricetidae). Micron. 2013;46:27–34.
Rovatsos MT, Marchal J, Romero-Fernández I, Fernández F, Giagia-Athanosopoulou E, Sánchez A. Rapid, independent, and extensive amplification of telomeric repeats in pericentromeric regions in karyotypes of arvicoline rodents. Chromosome Res. 2011;19(7):869–82.
Aksenova AY, Greenwell PW, Dominska M, Shishkin AA, Kim JC, Petes TD, et al. Genome rearrangements caused by interstitial telomeric sequences in yeast. Proc Natl Acad Sci. 2013;110(49):19866–71.
Chen TL, Manuelidis L. SINEs and LINEs cluster in distinct DNA fragments of Giemsa band size. Chromosoma. 1989;98(5):309–16.
Boyle AL, Ballard SG, Ward DC. Differential distribution of long and short interspersed element sequences in the mouse genome: chromosome karyotyping by fluorescence in situ hybridization. Proc Natl Acad Sci. 1990;87(19):7757–61.
Thomsen P, Miller J. Pig genome analysis: differential distribution of SINE and LINE sequences is less pronounced than in the human and mouse genomes. Mamm Genome. 1996;7(1):42–6.
Parish D, Vise P, Wichman H, Bull J, Baker R. Distribution of LINEs and other repetitive elements in the karyotype of the bat Carollia: implications for X-chromosome inactivation. Cytogenet Genome Res. 2002;96:191–7.
Scott L, Kuroiwa A, Matsuda Y, Wichman H. X accumulation of LINE-1 retrotransposons in Tokudaia osimensis, a spiny rat with the karyotype XO. Cytogenet Genome Res. 2006;112(3–4):261–9.
Carbone L, Harris RA, Mootnick AR, Milosavljevic A, Martin DI, Rocchi M, et al. Centromere remodeling in Hoolock leuconedys (Hylobatidae) by a new transposable element unique to the gibbons. Genome Biol Evol. 2012;4(7):760–70.
Laurent A, Puechberty J, Prades C, Gimenez S, Roizes G. Site-specific retrotransposition of L1 elements within human alphoid satellite sequences. Genomics. 1997;46(1):127–32.
Kapitonov VV, Holmquist GP, Jurka J. L1 repeat is a basic unit of heterochromatin satellites in cetaceans. Mol Biol Evol. 1998;15(5):611–2.
Acosta MJ, Marchal JA, Fernández-Espartero CH, Bullejos M, Sánchez A. Retroelements (LINEs and SINEs) in vole genomes: differential distribution in the constitutive heterochromatin. Chromosome Res. 2008;16(7):949–59.
Pagán HJ, Macas J, Novák P, McCulloch ES, Stevens RD, Ray DA. Survey sequencing reveals elevated DNA transposon activity, novel elements, and variation in repetitive landscapes among vesper bats. Genome Biol Evol. 2012;4(4):575–85.
McStay B, Grummt I. The epigenetics of rRNA genes: from molecular to chromosome biology. Annu Rev Cell Dev Biol. 2008;24:131–57.
Buckler ES, Ippolito A, Holtsford TP. The evolution of ribosomal DNA divergent paralogues and phylogenetic implications. Genetics. 1997;145(3):821–32.
Tchurikov NA, Fedoseeva DM, Sosin DV, Snezhkina AV, Melnikova NV, Kudryavtseva AV, et al. Hot spots of DNA double-strand breaks and genomic contacts of human rDNA units are involved in epigenetic regulation. J Mol Cell Biol. 2014;1–17.
Milhomem SS, Scacchetti PC, Pieczarka JC, Ferguson-Smith MA, Pansonato-Alves JC, O’Brien PC, et al. Are NORs always located on homeologous chromosomes? A FISH investigation with rDNA and whole chromosome probes in Gymnotus fishes (Gymnotiformes). PLoS ONE. 2013;8(2), e55608.
Santos N, Fagundes V, Yonenaga‐Yassuda Y, Souza MJ. Comparative karyology of Brazilian vampire bats Desmodus rotundus and Diphylla ecaudata (Phyllostomidae, Chiroptera): banding patterns, base‐specific fluorochromes and FISH of ribosomal genes. Hereditas. 2001;134(3):189–94.
Sotero-Caio CG, Pieczarka JC, Nagamachi CY, Gomes AJB, Lira TC, O’Brien PCM, et al. Chromosomal homologies among vampire bats revealed by chromosome painting (Phyllostomidae, Chiroptera). Cytogenet Genome Res. 2011;132(3):156–64.
Ribas TF, Rodrigues LR, Nagamachi CY, Gomes AJ, Benathar TC, Yang F, et al. Two new cytotypes reinforce that Micronycteris hirsuta Peters, 1869 does not represent a monotypic taxon. BMC Genet. 2013;14(1):119.
Noronha RCR, Nagamachi CY, Pieczarka JC, Marques-Aguiar S, Assis MDFLD, Barros RMDS. Meiotic analyses of the sex chromosomes in Carolliinae-Phyllostomidae (Chiroptera): NOR separates the XY1Y2 into two independent parts. Caryologia. 2004;57(1):1–9.
Gardner AL. The distributional status of some Peruvian mammals. Occasional papers of the museum of zoology. 1976;48:1–18.
Esbérard C, Bergallo H. Aspectos da biologia de Tonatia bidens (Spix) no estado do Rio de Janeiro, sudeste do Brasil (Mammalia, Chiroptera, Phyllostomidae). Rev Bras Biol. 2004;2:253–9.
Smith P, Pheasey H, Atkinson K, Miller J. Records of the Phyllostomine bats Tonatia bidens (Spix, 1823) and Lophostoma silvicolum d’ Orbigny, 1836 (Chiroptera, Phyllostomidae) associated with human dwellings in Paraguay. Chiroptera Neotropical. 2013;18(2):1139–43.
Bernard E, Fenton MB. Bat mobility and roosts in a fragmented landscape in central Amazonia. Brazil Biotropica. 2003;35(2):262–77.
Kirkpatrick M, Barton N. Chromosome inversions, local adaptation and speciation. Genetics. 2006;173(1):419–34.
Guerrero RF, Kirkpatrick M. Local adaptation and the evolution of chromosome fusions. Evolution. 2014;68(10):2747–56.
Livingstone K, Rieseberg L. Chromosomal evolution and speciation: a recombination‐based approach. New Phytol. 2004;161(1):107–12.
The authors are thankful to Lauren Gollahon for the cell culture facilities; Ron Chesser, Robert Bradley, Sean Rice, and Richard Strauss for insightful discussions and help with the statistical analysis; the Costa Rica field team: Caleb Phillips, Maria Sagot, Karina Monteiro, Julie Senawi, and Shannon Crawford; James Sowell and Santiago Burneo for funding and logistics of the Ecuador field expeditions, respectively, as well as the personnel associated with the field work; Texas Tech University and the Dept. of Biological Sciences for CGSC Doctoral Dissertation Completion Fellowship; TTUAB for the GIA fund. FY is supported by Wellcome Trust (WT098051).
The authors declare that they have no competing interests.
CGSC conducted the field work to Costa Rica to obtain TSA samples, carried out the FISH experiments, performed the analysis and interpretation of data, and drafted the manuscript. MV helped with the TSA FISH experiments, analysis and interpretation of data, and drafting of manuscript. FGH performed the BEAST analysis and drafted the manuscript. LAS generated the LINE-1 clone and helped drafting the manuscript, HAW helped with analysis and drafting of the manuscript; FY participated in the acquisition and interpretation of data and drafting of manuscript. RJB conceived the study, participated in its design and coordination, conducted the field work to collect the Ecuador samples, helped in the interpretation of data and with drafting the manuscript. All authors read and approved the manuscript.
Specimens used in the molecular analysis and respective GenBank accession numbers. (DOCX 24 kb)
Examples of in situ hybridizations using Macrotus californicus (MCA) chromosome paints on Phyllostominae metaphases. The MCA chromosome-specific probes used in TSA (a–c), LOC (d–f), and MCR (g–i) metaphases are indicated in pink in the upper-right corner of each picture. (JPEG 315 kb)
Molecular phylogram depicting the relationships among genera of bats in the subfamily Phyllostominae (highlighted in bold) and the divergence time estimates (node ages) from BEAST analysis of all genes. Nodes with posterior probability values higher than 0.98 in (a) are marked with an *. Abbreviations for node names are as follows: Phyllostominae + Glossophaginae ancestor (PGA), Phyllostominae ancestor (PA), Glossophaginae ancestor (GA); Phyllostomini ancestor (PiA); Lophostoma + Mimon ancestor (L/M); genus Tonatia (T); genus Anoura (A). (PDF 1529 kb)
Detailed calculations of unique fusions and fissions for ACU and TSA and chromosomal associations shared between the analyzed species. (DOCX 29 kb)
Long term effective population size (Ne) estimations using Lande (1979) model and the calculated rates of chromosomal rearrangement for TSA and ACU. (XLSX 140 kb)
About this article
Cite this article
Sotero-Caio, C.G., Volleth, M., Hoffmann, F.G. et al. Integration of molecular cytogenetics, dated molecular phylogeny, and model-based predictions to understand the extreme chromosome reorganization in the Neotropical genus Tonatia (Chiroptera: Phyllostomidae). BMC Evol Biol 15, 220 (2015). https://doi.org/10.1186/s12862-015-0494-y