Single Origin of Sex Chromosomes and Multiple Origins of B Chromosomes in Fish Genus Characidium

Abstract

Chromosome painting with DNA probes obtained from supernumerary (B) and sex chromosomes in three species of fish genus Characidium (C. gomesi, C. pterostictum and C. oiticicai) showed a close resemblance in repetitive DNA content between B and sex chromosomes in C. gomesi and C. pterostictum. This suggests an intraspecific origin for B chromosomes in these two species, probably deriving from sex chromosomes. In C. oiticicai, however, a DNA probe obtained from its B chromosome hybridized with the B but not with the A chromosomes, suggesting that the B chromosome in this species could have arisen interspecifically, although this hypothesis needs further investigation. A molecular phylogenetic analysis performed on nine Characidium species, with two mtDNA genes, showed that the presence of heteromorphic sex chromosomes in these species is a derived condition, and that their origin could have been unique, a conclusion also supported by interspecific chromosome painting with a CgW probe derived from the W chromosome in C. gomesi. Summing up, our results indicate that whereas heteromorphic sex chromosomes in the genus Characidium appear to have had a common and unique origin, B chromosomes may have had independent origins in different species. Our results also show that molecular phylogenetic analysis is an excellent complement for cytogenetic studies by unveiling the direction of evolutionary chromosome changes.

About the Authors:

José Carlos Pansonato-Alves

* E-mail: [email protected]

Affiliation: Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil

Érica Alves Serrano

Affiliation: Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil

Ricardo Utsunomia

Affiliation: Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil

Juan Pedro M. Camacho

Affiliation: Universidad de Granada, Departamento de Genetica, Granada, Spain

Guilherme José da Costa Silva

Affiliation: Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil

Marcelo Ricardo Vicari

Affiliation: Universidade Estadual de Ponta Grossa, Departamento de Biologia Estrutural, Molecular e Genética, Ponta Grossa, Paraná, Brazil

Roberto Ferreira Artoni

Affiliation: Universidade Estadual de Ponta Grossa, Departamento de Biologia Estrutural, Molecular e Genética, Ponta Grossa, Paraná, Brazil

Cláudio Oliveira

Affiliation: Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil

Fausto Foresti

Affiliation: Universidade Estadual Paulista (UNESP), Instituto de Biociências/IB, Departamento de Morfologia, Botucatu, São Paulo, Brazil

Introduction

Supernumerary (B) chromosomes are extra elements found in many eukaryotic genomes in addition to the standard (A) complement of chromosomes. They usually consist of highly repetitive DNA sequences, such as ribosomal DNA, satellite DNA, and mobile elements [1], [2], [3], while examples of B chromosomes carrying coding sequences [4], [5], [6], [7], functional ribosomal DNA [8], [9], and influencing sex determination [10] have recently been reported. In some cases, DNA sequences in the B chromosomes resemble those in the A chromosomes of the same genome (intraspecific origin) whereas, in others, they are more similar to DNA sequences in the genome of a relative species (interspecific origin) [11]. B chromosomes of intraspecific origin may be derived either from autosomes, as it is the case in Locusta migratoria [12], or from sex chromosomes, as in the frog Leiopelma hochstetteri [13], the fish Lithocromis rubripinnis [10] and the fly Drosophila albomicans [14]. In fact, sex and B chromosomes may share some features related to the distribution of chromatin, the accumulation of repetitive DNA and the loss of gene activity [11].

The presence of B chromosomes has been reported in several fish species [15], [16], [17], [10], [18]. Morphologically, these chromosomes are highly variable and range in size from very small, as in Moenkhausia sanctafilomenae [19] and Prochilodus scrofa [20], to very large, as in Astyanax scabripinnis [21], [22] and Alburnus alburnus [15]. In some cases, the repetitive sequences present on B chromosomes are shared with autosomal chromosomes, as in P. lineatus [23], [24] and A. scabripinnis [25], [26]. The B chromosomes of other species such as A. alburnus contain unique repetitive sequences [15].

Among Neotropical fish, the genus Characidium (of the family Crenuchidae) provides an interesting model for cytogenetic and evolutionary studies, particularly because of the presence of differentiated sex chromosome systems [27] and supernumerary chromosomes [28], [29], [30]. Despite the conserved diploid number of 50 chromosomes in all species in this group [29], interspecific and interpopulational differences have been reported in relation to either sex chromosomes of the ZZ-ZW type [31], [32], the location and number of rDNA sites [27], the occurrence of natural triploidy [33], [34] and the presence of B chromosomes [28], [29], [30].

In Characidium, it has been postulated that the sex chromosome systems found in several species have a common origin and that the observed differences between the Z and W chromosomes in different species and populations were caused by a variety of structural rearrangements, mainly involving rDNA sites and heterochromatic regions [31], [32]. However, though most species show a single rDNA-carrying chromosome pair, its chromosome location also shows variation [31]. Some species and populations are particularly interesting because they are highly differentiated: (i) the population of C. gomesi at Pardo River shows cells with up to four B chromosomes and lacks rDNA in the sex chromosomes [35]; (ii) the population of C. gomesi from the Tietê River has rDNA on the Z and W chromosomes and does not have any B chromosomes [30]; (iii) the population of C. oiticicai from the Paraitinga River has cells with up to 3 small extra chromosomes and shows NORs on the long arms of the W chromosome and on the short arms of the Z chromosome [29]; and (iv) a population of C. pterostictum from the Betari River exhibits cells with up to 3 small B chromosomes, NORs on the long arms of the W chromosome and on the short arms of the Z chromosome, and an acrocentric pair that is exclusive of this species [29]. Here we analyse the origin of B chromosomes in these three Characidium species by means of chromosome painting with DNA probes obtained from the B and W chromosomes, and infer the direction of chromosome changes by a mitochondrial DNA phylogeny built with these three Characidium species and six close relatives. The results have shown a unique origin for the sex chromosome systems in this genus, but multiple origins for B chromosomes.

Materials and Methods

Ethics Statement

Samples were collected in accordance with Brazilian environmental protection legislation (collection permission MMA/IBAMA/SISBIO - number 3245), and the procedures for sampling, maintenance and analysis of the samples were performed in compliance with the Brazilian College of Animal Experimentation (COBEA) and was approved (protocol 595) by the BIOSCIENCE INSTITUTE/UNESP ETHICS COMMITTEE ON USE OF ANIMALS (CEUA).

Origin of samples and cytogenetic analysis

For the phylogenetic analysis, we assessed a population of Crenuchus spilurus as an outgroup and the ingroup was composed by C. zebra, C. gomesi, C. alipioi, C. lauroi, C. oiticicai, C. pterostictum, C. schubarti, C. lanei and Characidium sp. (fig. 1). The animals were collected in accordance with Brazilian environmental laws (collection permission from MMA/IBAMA/SISBIO, number 3245). The collection, maintenance and analysis of the animals were performed in accordance with the international regulations on animal experimentation followed by the University of the State of São Paulo (CEEAA/IBB/UNESP Protocol number 304). For the cytogenetic analysis, the animals were anesthetised and dissected, and the mitotic chromosomes were obtained from renal tissue and gills using the technique described by Foresti et al. [36]. C-banding was performed following the protocol described by Sumner [37]. The samples to be analysed were identified, fixed and deposited in the fish collection of the Laboratório de Biologia e Genética de Peixes, Botucatu, São Paulo, Brazil (table 1).

[Figure omitted. See PDF.]

Figure 1. Map showing collection localities for each Characidium species.

https://doi.org/10.1371/journal.pone.0107169.g001

[Figure omitted. See PDF.]

Table 1. Characidium species and populations analyzed.

https://doi.org/10.1371/journal.pone.0107169.t001

Microdissection and FISH

For the microdissection of B chromosomes, we utilised cell suspensions from C. gomesi (Pardo River, registration number 33637), C. pterostictum (Betari River - registration number 37903) and C. oiticicai (Tietê River, registration number 31212) whose metaphases possessed a single B chromosome. For the microdissection of the W chromosome, we used a cell suspension from C. gomesi (Tietê River - registration number 31099). The B chromosomes were easily identified because they are smaller than the other chromosomes of the standard complement (A chromosomes) and have a distinct morphology (fig. 2B–D). The W chromosome from the Tietê River C. gomesi population has a secondary constriction in its long arm (fig. 2A). The cell suspensions were dropped on 24 mm×60 mm glass coverslips and stained with 5% Giemsa for 5 minutes, and the microdissection was performed using a micromanipulator with a glass needle (Eppendorf-5171) coupled to an inverted microscope (Zeiss Axiovert 100). Ten B chromosomes from each population and 12 W chromosomes from C. gomesi (Tietê River) were microdissected. These chromosomes were placed in four different microtubes (0.2 mL) with 9 µL of ultrapure water, and the DNA was amplified using the GenomePlex Single Cell Whole Genome Amplification Kit (WGA4 - Sigma) [38]. The reaction products were visualised on agarose gels. We generated probes for the B chromosomes (CgB - C. gomesi, CpB - C. pterostictum, CoB - C. oiticicai) and for the W chromosome (CgW - C. gomesi) from a re-amplification reaction of the amplified DNA and the GenomePlex WGA Re-amplification Kit (WGA3 - Sigma) using the modified nucleotide digoxigenin-11-dUTP (Roche Applied Science).

[Figure omitted. See PDF.]

Figure 2. Metaphase chromosome spreads after C-banding.

a) C. gomesi (Alambari River), b) C. gomesi (Paranapanema River), c) C. pterostictum (Betari River), d) C. oiticicai (Paraitinguinha River). The scale bar equals 10 µm.

https://doi.org/10.1371/journal.pone.0107169.g002

Fluorescent in situ hybridisation (FISH) experiments using probes specific to the B chromosomes were performed on four individuals per species, including B-carrying and B-lacking males and females. The in situ hybridisation was performed as described by Pinkel et al. [39], under high stringency conditions (mix containing 200 ng of each probe, 50% formamide, 10% dextran sulphate, 2xSSC at 37°C overnight). The chromosomes were counterstained with 4′, 6-diamidino-2-phenylindole dihydrochloride (DAPI) and analysed using an optical photomicroscope (Olympus BX61). The images were captured using the Image Pro plus 6.0 software (Media Cybernetics) and processed using the Adobe Photoshop CS2 program to improve their brightness and contrast.

Analysis of mitochondrial DNA

To extract DNA, the PromegaWizardGenomic DNA Purification Kit was used according to the manufacturer's instructions. We used partial Cytochrome oxidase C subunit1 (COI) gene and Cytochrome B gene sequences (Cyt B) for molecular phylogenetic analyses. For the amplification and sequencing of the gene segments, we used the primers FishF1 5′TCA ACC AAC CAC AAA GAC ATT GGC AC3′ and FishR1 5′TAG ACT TCT GGG TGG CCA AAG AAT CA 3′ [40] to COI and Cyt B L 14841 5′ CCA TCC AAC ATC TCA GCA TGA TGA AA 3′ and Cyt B H 15915b 5′ AAC CTC CGA TCT TCG GAT TAC AAG AC 3′ [41] to Cyt B.

The corresponding sequences were amplified by PCR in a total reaction volume of 12.5 µl with 35 cycles (240 s at 95°C, 45 s at 48–54°C, and 60 s at 72°C). For the PCR, we employed the Gotaq kit (Promega) according to the manufacturer's recommendations. The PCR products were analysed on 1% agarose gel and purified using ExoSAP-IT reagent (USB Corporation) according to the manufacturer's instructions. The purified products were then used to assemble sequencing reactions with the BigDye TM Terminator v3.1 Cycle Sequencing Ready Reaction Kit (Applied Biosystems) in accordance with the instructions provided in the user's manual. The reactions were purified, and the sequences were obtained using an ABI3130 Automated Capillary DNA Sequencer. The sequences were initially analysed using the ATGC program (Genetix Inc.) to obtain consensus sequences. The consensus sequences were aligned using the MUSCLE algorithm [42]. Each gene block was aligned separately and then concatenated into a single matrix. The matrix was partitioned into six partitions, one partition for each codon position of each gene. A matrix saturation test was conducted with the DAMBE program version 5.1.1 using the method described by Xia et al. [43] and in the same program by the transition/transversion rate graphic. The relationship between the sampled taxa using Maximum Likelihood (ML) analysis was performed using the ‘RAxML-HPC2onXSEDE’ tool [44]. ML analyses were conducted under GTR+G since RAxML only applies this model [44]. A bootstrap test with 1,000 pseudoreplicates [45] was used as a statistical test of the phylogeny, nodes with bootstrap values below 70% were collapsed. An independent Bayesian analysis was conducted. Two independent MCMC chains were run with 5,000,000 replicates each, sampling one tree every 1000 steps. The search for the best model of nucleotide evolution for each partition was performed using the model test 3.6 [46] using the Akaike information criterion [47]. The distribution of log likelihood scores was examined to determine stationarity for each search and to decide if extra runs were required to achieve convergence, using the program Tracer 1.4 [48]. Initial trees estimated prior to convergence were discarded as part of a burn-in procedure, and the remaining trees were used to construct a 70% majority rule consensus tree in Paup* [49].

Statistical Analysis

We compared the prevalence of B chromosomes between males and females, in each species, by means of contingency chi square tests, with Yates' correction, using the GraphPad PRISM 5.1 software.

Results

The standard karyotype of all Characidium species hitherto analysed consists of 50 standard (A) chromosomes. In addition, the three species analyzed here (C. gomesi - Pardo River, C. oiticicai- Paraitinguinha River and C. Pterostictum- Betari River) harbour mitotically unstable B chromosomes, i.e. varying in number among the cells within a single individual. In all three species, the B chromosomes were the smallest members of the karyotype. C. gomesi and C. oiticicai showed very similar size groups for the A chromosomes, with 32 metacentric (m) and 18 submetacentric (sm) chromosomes, whereas C. pterostictum showed 32 m, 16 sm and 2 acrocentric (a) chromosomes. All three species showed a ZZ/ZW sex chromosome system (Table 1), similar to most other Characidium species hitherto analysed.

Only one of the four C. gomesi populations analyzed, collected in the Pardo River (Table 1) showed B chromosomes, with 72% of individuals carrying them (table 2). Likewise, only one of the five C. pterostictum populations analyzed, from the Betari River, carried B chromosomes (62% of B-carriers) (table 2). Finally, 50% of the individuals analyzed from the single C. oiticicai population sampled (Paraitinga River) carried B chromosomes (table 2). Contingency chi square tests failed to show significant differences in B chromosome prevalence between females and males, in any of the three species (table 2).

[Figure omitted. See PDF.]

Table 2. Proportion of individuals carrying B chromosomes (prevalence).

https://doi.org/10.1371/journal.pone.0107169.t002

The C-banding technique revealed very similar patterns in the three species, especially for the W and B chromosomes, which were completely darkly C-banded in all three species (fig. 2). The Z chromosome showed a large pericentromeric C-band in all analysed populations of C. gomesi and C. oiticicai, but it showed an additional distal C-band in C. pterostictum. The autosomes showed pericentromeric C-bands in all three species. In C. pterostictum, there were also distal C-bands in the short arm of the submetacentric pair no. 22 and the long arm of the acrocentric pair no. 25.

Intraspecific chromosome painting with the probes obtained from B chromosome microdissection (table 3) showed remarkable differences between the species. In C. gomesi, the CgB probe hybridised with the B chromosomes, the pericentromeric region of the Z chromosome, and the entire W chromosome (fig. 3). Likewise, in C. pterostictum, the CpB probe hybridised with the B chromosomes, the distal region of the Z chromosome long arm, the complete W chromosome and the distal region of the acrocentric pair no. 25 (fig. 3). In C. oiticicai, however, the CoB probe only hybridised with the B chromosomes (fig. 3). Hybridization signals near the telomere of some A chromosomes without a defined pattern have been observed in experiments with B chromosome probes for C. gomesi and C. oiticicai, presumably due to the presence, in the B-probes, of some repetitive DNAs being also present in some A chromosomes.

[Figure omitted. See PDF.]

Figure 3. Karyotypes of Characidium species after chromosome painting with the CgW, CpB, CgB and CoB probes and counterstained with DAPI.

m =  Metacentric; sm =  Submetacentric; a =  Acrocentric.

https://doi.org/10.1371/journal.pone.0107169.g003

[Figure omitted. See PDF.]

Table 3. Summary of in situ hybridisation experiments with the B and W chromosome probes.

https://doi.org/10.1371/journal.pone.0107169.t003

Interspecific chromosome painting (table 3) with the CgB probe showed hybridisation with the B, Z, and W chromosomes and the acrocentric autosome pair no. 25 in C. pterostictum, but only with the sex chromosomes in C. oiticicai (fig. 3). Likewise, the CpB probe showed hybridisation signals on the sex and B chromosomes in C. gomesi, but only on the sex chromosomes in C. oiticicai (fig. 3). In high contrast, however, the CoB probe did not show any hybridisation signal on the C. gomesi and C. pterostictum chromosomes (table 3). Finally, the experiments using the CgW probe yielded the same results as those obtained with the CgB and CpB probes (fig. 3).

To infer the evolutionary direction of the observed chromosome changes, we build a molecular phylogeny for mitochondrial DNA using the Cyt B (GenBank acession numbers from KF914671 to KF914692) and COI (GenBank accession numbers from KF914693 to KF914710) genes on nine Characidium species and an outgroup. We analysed a total of 1,694 nucleotides, 610 of which were variable, and 458 were parsimoniously informative. There was no saturation since the Iss.c value was greater than the Iss value, and the information found in this matrix was used in the subsequent phylogenetic analysis [43], [50]. The graphic analysis of transition and transversion vs. genetic distance implemented by DAMBE 5.2.31 also indicated that the data were not saturated (R² = 0.79 for transition; R² = 0.90 for transversion). The best nucleotide substitution models selected for each partition after model test analysis are shown in Table S1. The topologies yielded by ML and Bayesian analyses were very similar, with high statistical support and posterior probability. The final topology, resulting from the strict consensus between the topologies of both analyses (fig. 4), showed the existence of two main clades, a basal one including all populations of C. zebra (a species lacking sex chromosomes), and the other clade including the remaining species (all of them carrying sex chromosomes). This finding strongly suggests that the origin of sex chromosomes in these species was unique, and that their presence is a derived state. In addition, the tree in fig. 4 shows that C. gomesi belongs to a clade that is rather distant from the one including C. pterostictum and C. oiticicai. This result indicates that the B chromosomes in C. gomesi and C. pterostictum presumably showed independent origins despite the high resemblance in their DNA content suggested by chromosome painting.

[Figure omitted. See PDF.]

Figure 4. Consensus topology obtained by maximum likelihood and Bayesian analysis of the concatenated mtDNA dataset.

The two numbers over each node represent the percentage of bootstrap obtained by ML and the posterior probability for that split obtained in the Bayesian analysis, respectively. When one of these indices is equal to 100, it is replaced by an asterisk.

https://doi.org/10.1371/journal.pone.0107169.g004

Discussion

Sex chromosome evolution

The molecular phylogenetic analysis performed on the mtDNA of nine Characidium species revealed the existence of two main clades, with C. zebra (a species lacking sex chromosomes) constituting the basal clade, and a derived clade including eight species showing wide geographical distribution [27], [51], [29], [31], [32] and carrying ZW♀/ZZ♂ sex chromosomes (see fig. 4). This suggests a unique origin for sex chromosomes in this genus, and this conclusion is reinforced by the fact that the CgW probe hybridized with the sex chromosomes of the three species analyzed here. The existence of the same repetitive DNAs in the Z and W chromosomes of, at least, the three Characidium species analyzed, thus suggests the possibility that these sex chromosomes arose in an ancestor species. A single origin has also been suggested in fish genus Triportheus [52], but other genera of Neotropical fish, such as Eigenmannia, exhibit the opposite situation, with different sex chromosome systems being observed in closely related species, with presumable different and recent origins [53]. The joint cytogenetic and phylogenetic analysis of more representatives of Characidium, and related genera, also including DNA probes obtained from the Z chromosome, would provide a more general picture on sex chromosome evolution in this group of species.

B chromosome evolution

The presence of B chromosomes in the genus Characidium appears to be limited to only a few populations in several species. Out of nine species with known cytogenetic data, B chromosomes have been found in a single population of C. gomesi (out of nine populations analysed), C. pterostictum (out of five populations analysed), C. zebra (out of seven populations analysed) and C. oiticicai (only one population analysed) [28], [33], [54], [55], [29], [30], [31]. The B chromosomes observed in these species were all mitotically unstable, but exhibited some differences regarding shape and heterochromatin distribution [28], [30], as expected for genomic elements evolving in separate species.

The DNA probes obtained from the B chromosomes of C. gomesi and C. pterostictum, and the W chromosome of C. gomesi, showed the same hybridisation pattern on the B, W and Z chromosomes of these two species, indicating that similar types of repetitive DNA make up their B and sex chromosomes. This probably reflects a descent relationship between these chromosomes and thus an intraspecific origin for these B chromosomes. As our molecular phylogenetic analysis has clearly shown, sex chromosomes were present in eight species constituting a well-defined clade, whereas B chromosomes were present in a single population of only three of these species. Therefore, sex chromosome origin preceded B chromosome origin, for which reason we can infer that the B chromosomes sharing DNA sequences with sex chromosomes in C. gomesi and C. pterostictum probably derived from sex chromosomes in each species. A possible caveat to this conclusion could be due to the fact that the DNA probes obtained by the GenomePlex method appear to contain a predominance of A+T rich DNA [56], and most painting signals observed by us were located on DAPI⁺ (i.e. A+T rich) regions. In fact, the probe obtained from the W chromosome (which contains rDNA, a kind of G+C rich DNA sequences) in C. gomesi (CgW) showed no hybridisation signals on the distal region of the autosome pair no. 17 (also carrying rDNA) of C. gomesi from the Pardo River population (fig. 3), suggesting that the probe did not contain rDNA sequences. But, in spite of this bias, it is clear that the repetitive DNAs included in the probes suggested a high resemblance in DNA content between B and sex chromosomes in C. gomesi and C. pterostictum, and thus these Bs could have conceivably arisen from sex chromosomes. Similar derivation from sex chromosomes has been shown in the frog Leiopelma hochstetteri [13]. Similarity of DNA sequences between B and A chromosomes, thus suggesting the intraspecific origin of B chromosomes, has also been reported in the fish species Astyanax scabripinnis [25], [26], Prochilodus lineatus [23] and Lithocromis rubripinnis [10], the grasshoppers Eyprepocnemis plorans [57], [58], [59] Podisma sapporensis [60] and Locusta migratoria [12] and, the mouse Apodemus peninsulae [61].

Alternatively, B chromosomes in these species could have derived from autosomes and later acquired the repetitive elements they share with sex chromosomes, but this event should have occurred in the two species separately, given their phylogenetic distance (see Fig. 4) and is thus unlikely.

In C. pterostictum, the CpB, CgB and CgW probes showed also hybridization on the heterochromatic blocks of the autosomal acrocentric pair no. 25. Therefore, in this species, B chromosomes could have also derived from this A chromosome. Nevertheless, the similarity of DNA sequences between the acrocentric pair and sex chromosomes in populations lacking B chromosomes, such as Miracatu, Morretes, Paranaguá and Jaraguá do Sul (results not shown) indicates that these shared sequences had settled in the acrocentric pair prior to B origin in the Betari River population. Therefore, given the highly dynamic nature of repetitive DNA sequences and their ability to spread across chromosomes within a same genome, the former conclusions should be taken with caution since chromosome painting results, based only in sequence similarity, could in fact be supporting a false idea of common descent between sex and B chromosomes. Comparison of nucleotide sequence will provide more precise conclusions at this respect.

Our molecular phylogenetic analysis suggests that the B chromosomes in C. gomesi and C. pterostictum most likely had independent origins, as these two species belong to very different clades. The fact that B chromosomes in both species are found in only one population could point to a recent origin which has not yet been followed by spread to other populations. An estimation of gene flow among populations in these two species would throw much light on this hypothesis. Alternatively, the B chromosomes in these two species could have derived from a B chromosome already present in a common ancestor species. But this would require the loss of the ancestor B chromosome in some species and populations, which is hardly parsimonious.

In contrast to the B probes obtained in the two other species, that obtained in C. oiticicai failed to hybridize with any chromosome other than the B chromosome itself, The lack of hybridization with any A chromosome has already been observed in chromosome painting experiments for a human supernumerary microchromosome, indicating the absence of detectable homologous sequences in the normal chromosomes, thus being the only supernumerary marker chromosome whose ancestry could not be determined [62]. These authors suggested possible mechanisms to explain the origin of this extra chromosome, including neocentromeric activity and either a complex amplicon of different genomic regions or the amplification of a very small region. Of course, these kinds of events cannot be ruled out to explain the origin of the B chromosome in C. oiticicai but, in fish, an alternative possibility is the interspecific origin through hybrydization. This has been suggested in several cases, such as the fish Poecilia formosa [63] or the wasp Nasonia vitripennis [64], and experimentally shown by Perfectti and Werren [65] inthe wasp Nasonia. The B chromosomes found in C. oiticicai showed some cytological similarities with those found in C. pterostictum, such as very similar size and heterochromatin content. But our chromosome painting and molecular phylogenetic results show that they can be B chromosomes independently originated in each species through different evolutionary pathways: intraspecific origin in C. pterostictum but interspecific origin in C. oiticicai. In other organisms, e.g. wasps, functionally equivalent B chromosomes showed independent origins, as shown for the paternal sex ratio (PSR) chromosomes in Nasonia vitripennis and Trichogramma kaykai [64], [66]. Ascertaining the precise introgressive hybridisation event that gave rise to the B chromosome in C. oiticicai is not an easy task, but a comparative analysis of DNA sequences contained in this B with those contained in the A chromosomes of this and other relative species, could be helpful [12], [67]. Two of the most interesting species for this research are C. zebra (which lives in sympatry with C. oiticicai) and C. lauroi, which is its closest relative species (see fig. 4). However, we cannot rule out that the amplification method used here may have privileged specific repetitive DNA sequences of B chromosomes being absent in the autosomes, and thus other DNA sequences presumably shared between As and Bs would have gone unnoticed with the employed methods.

Conclusions

Our present results indicate that whereas sex chromosomes in the genus Characidium appear to have had a common and unique origin, B chromosomes in C. gomesi, C. pterostictum and C. oiticicai may have had an independent origin in every species, in the two former most likely derived intraspecifically from sex chromosomes, whereas in C. oiticicai they most likely arose through interspecific hybridization.

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Nucleotide substitution models for each partition used in the phylogenetic analyses of each program.

https://doi.org/10.1371/journal.pone.0107169.s001

(XLSX)

Acknowledgments

We thank Renato Devidé for support during sample collection and the anonymous reviewers for their substantial contributions.

Author Contributions

Conceived and designed the experiments: JCPA EAS RU JPMC GJCS MRV RFA CO FF. Performed the experiments: JCPA EAS RU GJCS. Analyzed the data: JCPA ES RU GJCS JPMC MRV RFA. Contributed reagents/materials/analysis tools: JPMC CO FF. Contributed to the writing of the manuscript: JCPA EAS RU GJCS JPMC MRV RFA CO FF.

Citation: Pansonato-Alves JC, Serrano ÉA, Utsunomia R, Camacho JPM, Costa Silva GJd, Vicari MR, et al. (2014) Single Origin of Sex Chromosomes and Multiple Origins of B Chromosomes in Fish Genus Characidium. PLoS ONE 9(9): e107169. https://doi.org/10.1371/journal.pone.0107169

References

1. Green DM, Zeyl CW, Sharbel TF (1993) The evolution of hypervariable sex and supernumerary (B) chromosomes in the relict New Zealand frog, Leiopelma hochstetteri. J Evol Biol 6: 417–441.

2. Bougourd SM, Jones RNB (1997) Chromosomes: A Physiological Enigma. New Phytologist 137: 43–54.

3. Camacho JPM (2005) B chromosomes. In: Gregory TR editor. The evolution of the genome. Elsevier, San Diego pp. 223–286.

4. Graphodatsky AS, Kukekova AV, Yudkin DV, Trifonov VA, Vorobieva NV, et al. (2005) The proto-oncogene C-KIT maps to canid B chromosomes. Chromosome Res 13: 113–122.

5. Yudkin DV, Trifonov VA, Kukekova AV, Vorobieva NV, Rubtsova NV, et al. (2007) Mapping of KIT adjacent sequence on canid autosomes and B chromosomes. Cytogenet Genome Res 116: 100–103.

6. Banaei-Moghaddam AM, Meier K, Karimi-Ashtiyani R, Houben A (2013) Formation and expression of pseudogenes on the B chromosome of rye. Plant Cell 25(7): 2536–2544.

7. Trifonov VA, Dementyea PV, Larkin DM, O'Brien PCM, Perelma PL, et al. (2013) Transcription of a protein-coding gene on B chromosomes of the Siberian roe deer (Capreolus pygargus). BMC Biology 11: 90.

8. Ruíz-Estévez M, López-Léon MD, Cabrero J, Camacho JPM (2012) B-Chromosome Ribosomal DNA is functional in the grasshopper Eyprepocnemis plorans. PLoS ONE 7(5): e36600.

9. Ruíz-Estévez M, López-Léon MD, Cabrero J, Camacho JPM (2013) Ribosomal DNA is active in different B chromosome variants of the ghrasshopper Eyprepocnemis plorans. Genetica 141: 337–345.

10. Yoshida K, Terai Y, Mizoiri S, Aibara M, Nishihara H, et al. (2011) B Chromosomes have a functional effect on female sex determination in Lake Victoria cichlid fishes. PLoS Gen 7: e1002203.

11. Camacho JPM, Sharbel TF, Beukboom LW (2000) B-chromosome evolution. Phil. Trans. R. Soc. Lond 355: 163–178.

12. Teruel M, Cabrero J, Perfectti F, Camacho JPM (2010) B chromosome ancestry revealed by histone genes in the migratory locust. Chromosoma 119: 217–225.

13. Sharbel TF, Green DM, Houben A (1998) B-chromosome origin in the endemic New Zealand frog Leiopelma hochstetteri through sex chromosome devolution. Genome 41: 14–22.

14. Zhou Q, Zhu HM, Huang QF, Zhao L, Zhang GJ, et al. (2012) Deciphering neo-sex and B chromosome evolution by the draft genome of Drosophila albomicans. BMC Genomics 13: 109.

15. Ziegler CG, Lamatsch DK, Steinlein C, Engel W, Schartl M, et al. (2003) The giant B chromosome of the cyprinid fish Alburnus alburnus harbours a retrotransposon-derived repetitive DNA sequence. Chromosome Res 11: 23–35.

16. Oliveira C, Foresti F, Almeida-Toledo LF (2007) Karyotypic evolution in Neotropical fishes. In: Pisano E, Ozouf-Costaz C, Foresti F, Kapoor BG, editors. Fish Cytogenetics. Enfield NH, Science Publisher, pp. 111–164.

17. Poletto AB, Ferreira IA, Martins C (2010) The B chromosome of the African ciclid fish Haplochromis obliquidens harbor 18S rRNA gene copies. BMC Gen 11: 1–8.

18. Fantinatti BEA, Mazzuchelli J, Valente GT, Cabral-de-Mello DC, Martins C (2012) Genomic content and new insights on the origin of the B chromosome of the cichlid fish Astatotilapia latifasciata. Genetica 139: 1273.

19. Foresti F, Almeida-Toledo LF, Toledo AS (1989) Supranumerary chromosome system, C-banding pattern characterization and multiple organizer regions in Moenkhausia sanctaefilomenae (Pisces, Characidae). Genetica 79: 107–114.

20. Pauls E, Bertollo LAC (1981) Evidence for a system of supernumerary chromosomes in Prochilodus scrofa Steindachner (Pisces, Prochilodontidae). Caryologia 36: 307–314.

21. Maistro EL, Foresti F, Oliveira C, Almeida-Toledo LF (1992) Occurrence of macro B chromosomes in Astyanax scabripinnis paranae (Pisces, Characiformes, Characidae). Genetica 87: 101–106.

22. Salvador LB, Moreira-Filho O (1992) B chromosome in Astyanax scabripinnis (Pisces, Characidae). Heredity 69: 50–56.

23. Jesus CM, Galetti Jr PM, Valentini SR, Moreira-Filho O (2003) Molecular characterization and chromosomal localization of two families of satellite DNA in Prochilodus lineatus (Pisces, Prochilodontidae), A species with B chromosomes. Genetica 118: 25–3.

24. Artoni RF, Vicari MR, Endler AL, Cavallaro ZI, Jesus CM, et al. (2006) Banding pattern of A and B chromosomes of Prochilodus lineatus (Characiformes, Prochilodontidae), with comments on B chromosomes evolution. Genetica 127: 277–284.

25. Mestriner CA, Galetti Jr PM, Valentini SR, Ruiz IRG, Abel LDS, et al. (2000) Structural and functional evidence that a B chromosome in the characid fish Astyanax scabripinnis is an isochromosome. Heredity 85: 1–9.

26. Vicari MR, Pistune HFM, Castro JP, Almeida MC, Bertollo LAC, et al. (2011) New insights on the origin of B chromosomes in Astyanax scabripinnis obtained by chromosome painting and FISH. Genetica 139: 1073–1081.

27. Vicari MR, Artoni RF, Moreira-Filho O, Bertollo LAC (2008) Diversification of a ZZ/ZW sex chromosome system in Characidium fish (Crenuchidae, Characiformes). Genetica 134: 311–317.

28. Miyazawa CS, Galetti Jr PM (1994) First cytogenetical studies in Characidium species (Pisces: Characiformes, Characidiinae). Cytologia 59: 73–79.

29. Pansonato-Alves JC, Paiva LRS, Oliveira C, Foresti F (2010) Interspecific chromosomal divergences in the genus Characidium (Teleostei: Characiformes: Crenuchidae). Neotrop Ichthyol 8: 77–86.

30. Pansonato-Alves JC, Vicari MR, Oliveira C, Foresti F (2011a) Chromosomal diversification in populations of Characidium cf. gomesi (Teleostei, Crenuchidae). J Fish Biol 78: 183–194.

31. Machado TC, Pansonato-Alves JC, Pucci MB, Nogaroto V, Almeida MC, et al. (2011) Chromosomal painting and ZW sex chromosomes differentiation in Characidium (Characiformes, Crenuchidae). BMC Genet 12: 65.

32. Pucci MB, Barbosa P, Nogaroto V, Almeida MC, Artoni RF, et al.. (2014) Population differentiation and speciation in the genus Characidium (Characiformes: Crenuchidae): effects of reproductive and chromosomal barriers. Biol J Linn Soc (in press). DOI: 10.1111/bij.12218.

33. Centofante L, Bertollo LAC, Moreira-Filho O (2001) Comparative cytogenetics among sympatric species of Characidium (Pisces, Characiformes): Diversity analysis with the description of a ZW sex chromosome system and natural triploidy. Caryologia 5: 253–260.

34. Pansonato-Alves JC, Oliveira C, Foresti F (2011b) Karyotypic conservatism in samples of Characidium cf. zebra (Teleostei, Characiformes, Crenuchidae): Physical mapping of ribosomal genes and natural triploidy. Genet Mol Biol 34: 208–213.

35. Maistro EL, Jesus CM, Oliveira C, Moreira-Filho O, Foresti F (2004) Cytogenetic analysis of A-, B- chromosomes and ZZ/ZW sex chromosomes of Characidium gomesi (Teleostei, Characiformes, Crenuchidae). Cytologia 69: 181–186.

36. Foresti F, Almeida-Toledo LF, Toledo-Filho SA (1981) Polymorphic nature of nucleolus organizer regions in fishes. Cytogenet Cell Gen 31: 37–144.

37. Sumner AT (1972) A simple technique for demonstrating centromeric heterochromatin. Exp Cell Res75: 304–306.

38. Gribble S, Ng BL, Prigmore E, Burford DC, Carter NP (2004) Chromosome paints from single copies of chromosomes. Chromosome Res 12: 143–151.

39. Pinkel D, Straume T, Gray JW (1986) Cytogenetic analysis using quantitative, high-sensitivity, fluorescence hybridization. Proc Natl Acad Sci 83: 2934–2938.

40. Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PDN (2005) DNA barcoding Australia's fish species. Philos Trans R Soc Lond B Biol Sci 360: 1847–1857.

41. Kocher TD, Thomas WK, Meyer A, Edwards SV, Paabo S, et al. (1989) Dynamics of mitochondrial DNA evolution in animals: Amplification and sequencing with conserved primers. Proc Nat Acad Sci 86: 6196–6200.

42. Edgar RC (2004) MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res 32: 1792–1797.

43. Xia X, Xie Z, Salemi M, Chen L, Wang Y (2003) An index of substitution saturation and its application. Mol Phylogenet Evol 26: 1–7.

44. Stamatakis A, Hoover P, Rougemont J (2008) A Rapid Bootstrap Algorithm for the RAxML Web Servers. Syst Biol 57: 758–771.

45. Felsenstein J (1985) ‘Confidence limits on phylogenies: An approach using the bootstrap’. Evolution Vol. 39: 783–791.

46. Posada D, Crandall KA (1998) MODELTEST: testing the model of DNA substitution. Bioinformatics 14: 817–818.

47. Posada D, Buckley TR (2004) Model selection and model averaging in phylogenetics: Advantages of Akaike Information Criterion and Bayesian approaches over likelihood ratio tests. Syst Biol 53: 793–808.

48. Rambaut A, Drummond AJ (2004) Tracer v1 4 [http://beast.bio.ed.ac.uk/Tracer].

49. Swofford DL (2003) PAUP*: Phylogenetic analysis using parsimony (*and other methods). Version 4. Massachusetts: Sinauer Associates.

50. Xia X, Lemey P (2009) Assessing substitution saturation with DAMBE. In The Phylogenetic Handbook: A Practical Approach to DNA and Protein Phylogeny. 2nd edition. Edited by Lemey P, Salemi M, Vandamme A-M. Cambridge: University Press: 615–630.

51. Noleto RB, Amorin, AP Vicari MR, Artoni RF, Cestari MM (2009) An unusual ZZ/ZW sex chromosome system in Characidium fishes (Crenuchidae, Characiformes) with the presence of rDNA sites. J Fish Biol 75: 448–453.

52. Artoni RF, Bertollo LAC (2002) Evolutionary aspects of the ZZ/ZW sex chromosome system in the Characidae fish, genus Triportheus. A monophyletic state and NOR. Heredity 89: 15–19.

53. Henning F, Moyses CB, Calcagnotto D, Meyer D, Almeida-Toledo LF (2011) Independent fusions and recent origins of sex chromosomes in the evolution and diversification of glass knife fishes (Eigenmannia). Heredity 106: 391–400.

54. Centofante L, Bertollo LAC, Buckup PA, Moreira-Filho O (2003) Chromosomal divergence and maintenance of sympatric Characidium fish species (Crenuchidae, Characidiinae). Hereditas 138: 213–218.

55. Silva AR, Maistro EL (2006) Cytogenetic divergence between two sympatric species of Characidium (Teleostei, Characiformes, Crenuchidae) from the Machado river, Minas Gerais, Brazil. Genet Mol Biol 29: 459–463.

56. Teruel M, Cabrero J, Montiel EE, Acosta MJ, Sánchez A, et al. (2009) Microdissection and chromosome painting of X and B chromosomes in Locusta migratoria. Chromosome Res 17: 11–18.

57. López-León MD, Neves N, Schwarzacher T, Heslop-Harrison TS, Hewitt GM, et al. (1994) Possible origin of a B chromosome deduced from its DNA composition using double FISH technique. Chromosome Res 2: 87–92.

58. Teruel M, Cabrero J, Perfectti F, Acosta MJ, Sánchez A, et al. (2009) Microdissection and Chromosome Painting of X and B Chromosomes in the Grasshopper Eyprepocnemis plorans. Cytogenet Genome Res 125: 286–291.

59. Teruel M, Ruiz-Ruano FJ, Marchal JA, Sanchez A, Cabrero J, et al. (2014) Disparate molecular evolution of two types of repetitive DNAs in the genome of the grasshopper Eyprepocnemis plorans. Heredity 112: 531–542.

60. Bugrov AG, Karamysheva TV, Pyatkova MS, Rubtsov DN, Andreenkova OV, et al. (2003) B chromosomes of the Podisma sapporensis Shir (Orthoptera, Acrididae) analyzed by chromosome microdissection and FISH. Folia Biol Krakow 51: 1–11.

61. Matsubara K, Yamada K, Umemoto S, Tsuchiya K, Ikeda N, et al. (2008) Molecular cloning and characterization of the repetitive DNA sequences that comprise the constitutive heterochromatin of the A and B chromosomes of the Korean field mouse (Apodemus peninsulae, Muridae, Rodentia). Chromosome Res 16: 1013–1026.

62. Mackie Ogilvie C, Harrison RH, Horsley SW, Hodgson SV, Kearney L (2001) A mitotically stable marker chromosome negative for whole chromosome libraries, centromere probes and chromosome specific telomere regions: a novel class of supernumerary marker chromosome? Cytogenet Cell Genet 92: 69–73.

63. Schartl M, Nanda I, Schlupp I, Wilde B, Epplen JT, et al. (1995) Incorporation of subgenomic amounts of DNA as compensation for mutational load in a gynogenetic fish. Nature 373: 68–71.

64. McAllister BF, Werren JH (1997) Hybrid origin of a B chromosome (PSR) in the parasitic wasp Nasonia vitripennis. Chromosoma 106: 243–253.

65. Perfectti F, Werren JH (2001) The interspecific origin of B-chromosomes: experimental evidence. Evolution 55: 1069–1073.

66. Van Vugt JFA, Jong JH, Stouthamer R (2009) The origin of a selfish B chromosome triggering paternal sex ratio in the parasitoid wasp Trichogramma kaykai. Proc R Soc B 276: 4149–4154.

67. Silva DMZA, Pansonato-Alves JC, Utsunomia R, Arraya-Jaime C, Ruiz-Ruano FJ, et al. (2014) Delimiting the origin of a B chromosome by FISH mapping, chromosome painting and DNA sequence analysis in Astyanax paranae (Teleostei, Characiformes). PLoS ONE 9(4): e94896.