- Research article
- Open Access
Integrating cytogenetics and genomics in comparative evolutionary studies of cichlid fish
BMC Genomics volume 13, Article number: 463 (2012)
The availability of a large number of recently sequenced vertebrate genomes opens new avenues to integrate cytogenetics and genomics in comparative and evolutionary studies. Cytogenetic mapping can offer alternative means to identify conserved synteny shared by distinct genomes and also to define genome regions that are still not fine characterized even after wide-ranging nucleotide sequence efforts. An efficient way to perform comparative cytogenetic mapping is based on BAC clones mapping by fluorescence in situ hybridization. In this report, to address the knowledge gap on the genome evolution in cichlid fishes, BAC clones of an Oreochromis niloticus library covering the linkage groups (LG) 1, 3, 5, and 7 were mapped onto the chromosomes of 9 African cichlid species. The cytogenetic mapping data were also integrated with BAC-end sequences information of O. niloticus and comparatively analyzed against the genome of other fish species and vertebrates.
The location of BACs from LG1, 3, 5, and 7 revealed a strong chromosomal conservation among the analyzed cichlid species genomes, which evidenced a synteny of the markers of each LG. Comparative in silico analysis also identified large genomic blocks that were conserved in distantly related fish groups and also in other vertebrates.
Although it has been suggested that fishes contain plastic genomes with high rates of chromosomal rearrangements and probably low rates of synteny conservation, our results evidence that large syntenic chromosome segments have been maintained conserved during evolution, at least for the considered markers. Additionally, our current cytogenetic mapping efforts integrated with genomic approaches conduct to a new perspective to address important questions involving chromosome evolution in fishes.
Integrated genome maps became a powerful tool to fill the gaps generated by the low resolution of linkage mapping and the problems of genome sequencing and assembly, providing a more accurate scenario of the genome structure. The application of fluorescent in situ hybridization (FISH) based in the use of bacterial artificial chromosome (BAC) clones as probes represents an efficient approach to anchor genomic and linkage data on physical chromosomes. BAC libraries have been explored for many aspects of molecular and genomic studies, such as the positional cloning of genes , comparative studies of synteny and gene organization among different species , as well as for local or whole genome physical and genetic mapping and sequencing . In cytogenetic research and chromosome mapping, the potential of BACs for animal genome analyses has improved, since several good quality genomes are already available in public databases  bringing up the possibility to refine the chromosome analysis under the focus of BAC-FISH mapping [5–12].
Cichlid fishes represent a monophyletic group classified in 4 subfamilies: Etroplinae (India and Madagascar), Ptychochrominae (Madagascar), Cichlinae (Neotropical region) and Pseudocrenilabrinae (Africa) [13, 14]. Although African Pseudocrenilabrinae cichlids can be separated in three major groups (hemichromines, haplochromines, and tilapiines), these groups are not recognized as valid taxonomic units . The African groups of lakes Victoria, Malawi and Tanganika represent a classical example of extensive and rapid radiation and therefore are highly interesting for evolutionary biology analyses . The karyotypes of more than 135 cichlid species have been determined, evidencing a diploid chromosome number ranging from 2n = 32 to 2n = 60. African cichlids have a modal diploid chromosome number of 44 chromosomes whereas the Neotropical cichlids contain 2n = 48 chromosomes reviewed by .
Genetic maps containing a few hundred polymorphic markers provide a starting point in order to resolve the chromosomal location of cloned genes or markers. Genetic linkage maps have been developed for a number of fish species, including zebrafish , medaka , catfish , rainbow trout , Atlantic salmon [22, 23], and also cichlids as the Nile tilapia Oreochromis niloticus (tilapiine) , Lake Malawi haplochromines , and Astatotilapia burtoni. There are also other genetic/genomic resources for cichlids including extensive collections of expressed sequence tags (ESTs) for Lake Victoria haplochromines [27, 28], A. burtoni[29, 30], and O. niloticus, BAC libraries for Nile tilapia , and haplochromines from lakes Malawi , Victoria , and Tanganika , and a high-resolution map for Nile tilapia . Altogether, these genomic resources have driven the investigation of several aspects of cichlid’s biology, including sex determination [24, 37–39]. Tilapiines have an XY sex chromosome system on linkage group (LG) 1 or a ZW system on LG3 [7, 40]. However, at least two distinct genetic sex determination systems in the Lake Malawi cichlids were found: a XY sex-determination system on LG7 and a ZW system on LG5 . Furthermore, some genes observed in LG1, CYP19a and WT1, are involved in the sexual differentiation of mammals [42, 43]; the CLCN5 gene (associated with renal disorder in humans), located on LG3, turns to be interesting since it was detected in the human X chromosome ; the opsin genes responsible for the color-spectrum vision of cichlids and thus involved with sexual selection and adaptation to new environments are located in LG5 . Additionally, the genome of O. niloticus has been recently sequenced , leading to the opportunity to integrate nucleotide sequence information and other genetic data. In this context, we address the comparative analysis of LG1, 3, 5, and 7 of cichlids based on the integration of genomic and molecular cytogenetic. Furthermore, the combined genomic/cytogenetic information obtained for cichlids was also comparatively analyzed with other vertebrates. Our results provide evidence for extensive synteny conservation of segments among Pseudocrenilabrinae cichlids and also between cichlids and other vertebrates. Such information is promising in the establishment of a framework for additional genome-wide studies.
The chromosome diploid number and morphology were analyzed in Giemsa-stained metaphases (Table 1). Such analysis confirmed the previous karyotype data available for cichlids composed of meta and submetacentric (m/sm), and telo/acrocentric (t/a) chromosomes [17, 47]. Metaphases of males and females of O. niloticus (the species source of the BAC clones used for FISH) were used to set up the BAC-FISH mapping experiments (Figure 1). After that, a comparative mapping, including eight other Pseudocrenilabrinae species, was carried out to elucidate the evolutionary history of chromosomes carrying LG1, 3, 5, and 7 (Figure 2, Additional file 1).
Four BACs from LG1 mapped on the long arm of a small t/a chromosome pair of O. niloticus (Figure 1A). The largest chromosome of O. niloticus evidenced labeling signals of eight BACs of LG3 (Figure 1B) and the LG3 mapping confirmed the presence of a lot of repetitive DNA in the end of the largest chromosome of all analyzed species (Figure 2).
Three BACs from LG5 mapped on a medium t/a chromosome different from that containing the LG1 (Figure 1C), and three BACs from LG7 mapped on the second largest pair of O. niloticus (Figures 1D). No differences were found between males and females for any of the hybridized BAC-probes.
At least two individuals of each species were used in BAC-FISH experiments and the signal patterns observed for LG1, 3, 5, and 7 were conserved in number and position when compared to the reference species O. niloticus (Figure 2; See Additional file 1). No results were obtained for BACs of LG7 in Labeotropheus trewavasae due to the low quality chromosome preparations of this species for FISH procedures (Figure 2). In South American cichlids (Cichlinae), none of the BAC probes produced identifiable chromosomal signals. Furthermore, no hybridization signal was observed in the Asian Etroplinae species, Etroplus maculatus.
Besides the conservation of the studied LGs, some differences were observed in the chromosomes’ morphology of the species (Figure 2): Astatotilapia burtoni showed a different pattern for LG5 which is located on the short arm of a large m/sm chromosome instead of a smal t/a chromosome as observed in the other cichlids. The single BAC signal observed in the short arms of the chromosomes was related to repeated sequences as observed for the UNH115 marker in LG3 of Metriaclima lombardoi (Figure 2; Additional file 1). Additionally, small variations in the chromosomal position of the markers were also observed (See Additional file 1), but these could be just a consequence of variation in the chromosome condensation of the analyzed individuals.
The chromosomal positions of all selected BACs were correlated with their corresponding order on their linkage maps. This enabled the association of LGs with chromosomes and to orientate their position according to the short/long arm of the corresponding chromosomes. Using the convention that the short arm is north oriented to and the long arm is south oriented, the north/south orientation of the current Oreochromis LG1 matches with the short/long arm (Figure 3); however, the orientation of LG3, 5, and 7 should be inverted in relation to the order of previously published molecular markers (Figures 4, 5 and 6).
In silico comparative analyses
Comparative genomic analyzes of LG1 at Bouillabase genome browser suggested a conserved pattern of these markers in medaka chromosome 3, stickleback chromosome 2, and pufferfish chromosome 5 (Figure 3) (Table 2). Additionally, using Ensembl and Genomicus databases, the chromosomal position of the genes WT1b and CYP19A1, mapped in LG1, was comparatively analyzed among fishes and also other vertebrates (Figure 3). The present analyses showed that WT1b and CYP19A1 genes are located at chromosome 18 of Danio rerio (zebrafish) (Figure 3) and separated by 6Mb. In medaka and pufferfish, these genes are separated by 10Mb and 5Mb, respectively (Table 2). A divergent scenario occur in other vertebrates, in which these genes appear in different chromosomes as observed in chicken and primates, where both markers are in chromosomes 10 and 5, and chromosomes 15 and 11, respectively (Figure 3).
LG3 proved to be conserved with at least two correspondent markers between the fish species medaka (chromosome 18) and stickleback (chromosome 7). In the pufferfish genome, markers of LG3 mapped at a non identified region (Table 2). Despite the synteny of LG3 markers identified in medaka and stickleback, an inversion of the markers UNH115 and UNH180, separeted by 8 Mb, was detected in medaka in relation to Oreochromis (Figure 4).
LG3 contains few known genes that could be used for comparative analyses in other vertebrates. CLCN5 was the only gene identified in other vertebrates. Using CLCN5 in a comparative analysis against the zebrafish genome, this gene appears duplicated and located at chromosomes 7 and 21. On the other hand, in all mammal species, CLCN5 has been conserved in X chromosome (Figure 4).
The majority of LG5 markers corresponds to chromosome 5 in medaka, chromosome 17 in stickleback, and chromosome 11 in pufferfish. Although other LG5 markers (Green OPSIN and c-Ski) were found in stickleback and medaka genome, it was not possible to categorize their chromosome position (Table 2). The position of the MME gene differs, as it can appear in chromosome 13, 1, or 16, respectively, in the three mentioned species (Table 2) (Figure 5). With the addition of zebrafish in Ensembl datasheet, an extremely conserved pattern of LG5 in chromosome 11 was evidenced. Once more, a duplicated marker was found in zebrafish (see c-Ski, for instance, in Figure 5). Genes observed in LG5 (AKR, ATP and c-Ski) have their positions conserved in chromosome 1 of primates. Additionally, this conservation is also related to horse chromosome 2, mouse chromosome 4, chicken chromosome 21, and turkey chromosome 23 (Figure 5). Despite located in the same chromosome, the position of these genes are inconsistent - in O. niloticus these genes appear tight together at 17 cM; however, in other species they appear separated, sometimes around 70Mb apart (data not show). The OPN gene (Blue opsin), present in the end of O. niloticus LG5, matched in the X chromosome for all analyzed mammals (Figure 5).
Analysis of LG7 markers in BouillaBase showed a block of synteny between Oreochromis, stickleback chromosome 19 and pufferfish chromosome 13, even taking into account that the distances between these markers are similar among the species. In medaka genome, only the UNH179 marker was identified in chromosome 6 (Figure 6). It was impossible to identify the position of the IGFBP2 marker in medaka and stickleback genomes, but in pufferfish it is located together with DUSP6 (Table 2), as also observed for O. niloticus (Table 2). In Ensembl comparative studies, the analysis of a small region (about 3cM) of LG7 that contains the genes RERG, KCNE1L, DUSP6, and IGFBP2, revealed a conserved position of RERG and DUSP6 at the same chromosome in vertebrates: Homo sapiens chromosome 12, Macaca mulatta chromosome 11, Bos taurus chromosome 5, Gallus gallus and Meleagris gallopavo chromosome 1, and indeed in the reptile Anolis carolinensis chromosome 5. However, these genes are separated by 10 to 75 Mb (data not shown) in non-fish species (Figure 6). Moreover, the gene KCNE1L was conserved in the X chromosome of all mammals (Figure 6).
Cytogenetic mapping and chromosome stability in Pseudocrenilabrinae
The most important find in our study is the extreme chromosome conservation observed within African cichlids. It is already known that African cichlids (Pseudocrenilabrinae) comprise about 1.400 species  and that their modal chromosome number is 2n = 44 . Extensive comparative mapping has demonstrated that the genetic maps for tilapia and Malawi cichlids are almost perfectly collinear . This conservation could be evidenced through comparative physical chromosomal mapping of linkage groups as shown in the present study. On the other hand, variability at the cytogenetic level among cichlids was mostly observed based in repetitive DNA chromosome mapping [49, 50] and seems to reflect the evolutionary dynamics of the repetitive genomic fraction and not broad processes acting in the whole karyotype.
Within Pseudocrenilabrinae, the major karyotype feature in the tilapiines is the presence of one large t/a chromosome pair (LG3), which is significantly larger than all other elements of the karyotype [17, 51]. On the other hand, the haplochromine and hemichromine cichlids have shown two outstanding chromosome pairs [17, 52] - the first and the second largest elements that correspond to LG3 and LG7, respectively (Figure 2). The chromosome mapping of LG3 and LG7 markers (present report) in Pseudocrenilabrinae species confirms the proposed homology between the two largest chromosome pairs of tilapiines and haplochromines/hemichromines . However, the hypothesis concerning an independent chromosome fusion that originated the largest chromosomes in tilapiines and haplochromines  is not in agreement with the results observed for the cytogenetic mapping of LG3 and LG7 markers.
It is known that the largest chromosome pair of O. niloticus originated by a centric fusion event of three other pairs of the ancestral cichlid karyotype composed of 48 acrocentric chromosomes . A previous study  proposed that a first chromosome fusion took place before the divergence of the main East African cichlid groups. The second chromosome fusion occurred independently in the tilapiines and non-tilapiines. In the tilapiines, a new chromosome was fused to the largest pair, and in the non-tilapiines, the second fusion did not involve the largest chromosome, but two other chromosomes, which gave rise to the chromosome pair 2. However, when we compared the distribution of BAC signals (LG3 and LG7) through the long arm of the two largest chromosomes in tilapiines and non-tilapiines species, we detected that the second chromosome fusion seems to be identical in these two groups, producing a chromosome with the same genomic content that remained conserved in all Pseudocrenilabrinae species so far analyzed (Figure 2).
The differences observed in the chromosome size of the first pair between tilapiines and haplochromines could be due to the intense dynamics of repeated DNAs located in the entire long arm instead of differences related to possible chromosome fusions. A range of studies has shown that the largest chromosome in O. niloticus presents a great amount of heterochromatin, microsatellites, transposable elements, LINES and SINES, and non-LTR retrotransposons [52–59]. The positive results using LG3-BACs (enriched of repetitive DNAs, see Table 3) evidenced the accumulation of repeated DNAs in the LG3 chromosome of all Pseudocrenilabrinae, specially in tilapiine species.
Despite displaying a highly conserved karyotype structure, many events of duplication, inversion, and fusion occurred during the diversification of African cichlids and resulted in diploid chromosome numbers divergent from the pattern of 44 chromosomes, as observed in Tilapia mariae with 2n = 40 [52, 60], T. sparrmanii with 2n = 42 , Oreochromis alcalicus with 2n = 48 , O. karongae with 2n = 38 , and Astatotilapia burtoni with 2n = 40 , among others.
A. burtoni has two metacentric chromosome pairs, which are probably the result of centric fusions of two small t/a chromosome pairs. According to our results, the ancestral t/a chromosome that contains the LG5 (observed in other haplochromines) might be involved in one of these events resulting in a m/sm chromosome. Such assumption could not be inferred for T. mariae (tilapiine), in which the diploid chromosome number reduction to 2n = 40 must be the result of rearrangements that did not involve the linkage groups investigated here (Figure 2).
Additionally, the reduction of chromosome number observed in O. karongae (2n = 38) and the presence of three pairs of medium- sized chromosomes, that are not found in the typical Oreochromis species, were originated by chromosome fusion events involving LG1 [10, 61] and could represent recent chromosomal rearrangements that have occurred independently in tilapiine and haplochromine groups. The additional signal observed in the 1p arm (LG3 – UNH115) of the largest chromosome of males and females of M. lombardoi could be a consequence of specific rearrangements or, even more, the presence of chimeric BAC inserts.
The absence of BAC-FISH signals in South American cichlids (Cichlinae) belonging to different tribes (Table 1) could be associated to genomic rearrangements that have disrupted in a small-scale level the genomic blocks carried by the BAC clones of O. niloticus in relation to these cichlines. Instead of being a monophyletic group , Neotropical cichlids harbor significantly higher levels of genetic variation compared to the African Pseudocrenilabrinae group . Although genomic rearrangements seem to have occurred differentiating South American and African cichlids, it is plausible that large genomic blocks are still conserved between them as well as in relation to other fish groups. Unfortunately, deep analyses integrating cytogenetics and genomic data were not possible for Cichlinae species since there is no available large-scale genome information for this subfamily so far.
Comparative analyses of vertebrates
Using BAC-end sequence data available at BouillaBase, it was possible to perform comparative analyses and detect regions of synteny among African cichlids and model fish species, and also vertebrates. In this work, conserved (syntenic) chromosome segments have been successfully identified by means of comparative cytogenetics between cichlid markers (mainly LG5 and 7) and mammals and birds chromosomes. The marker distances are very similar among the analyzed teleosts. However, considering non-fish groups, the markers’ distances are, in most cases, very divergent, evidencing that genomes are suffering rearrangements, despite the maintenance of conserved large genomic blocks as part of the same LG. Although large genomic blocks are conserved among vertebrates, variability in the genes/DNA sequences harbored in these regions is expected.
Comparative BAC mapping has already shown that some stickleback chromosomes have a nearly complete synteny with those of cichlids. Although stickleback does not belong to the order Perciformes, it is considered to be the closest related genome to cichlids among the currently available sequenced teleost genomes  and is currently the best reference sequence for assembling comparative maps of tilapia .
A range of studies using comparative cytogenetics has demonstrated that karyotypes are conserved in a macro scale throughout vertebrates. In Canidae, for example, despite the extensive variation in chromosome numbers and morphology, the majority of conserved chromosome segments appear to have remained largely intact in the karyotypes of extant canid species, although the relative orientation and distances are not always conserved . Conserved karyotypes are also observed in birds, which show a slow rate of interchromosomal rearrangements , and within the reptile family Scincidae, a character that was shown by cross-species chromosome painting . Similarly, at different taxonomic levels, comparative gene mapping has revealed a highly conserved linkage homology between an agamid lizard (Leiolepis reevesii) and a snake (Elaphe quadrivirgata) .
Conservation of large genomic blocks were also identified among different species of the Salmonidae family based on BAC-FISH blocks , and on SOX genes regions in cichlids that revealed a large genomic block that was conserved through vertebrates . These conservations of gene orders at scales of several Mb in diverse vertebrate groups permit the use of relatively complete sequences of model fish species to accelerate gene discovery and positional cloning of non-model species [63, 68]. The synteny observed in different vertebrates could be due to intrinsic chromosomal properties that confer selective pressure on large genomic blocks that should be preserved from major changes .
There is a classical idea that fish genomes have high rates of chromosomal rearrangements compared to other vertebrates and then, probably, evidence low rates of synteny [70–73]. However, our data associated to previous genomic studies are revealing that when we compare the genomic synteny among teleosts and other vertebrates in a macro scale level, large syntenic blocks can be clearly identified.
Comparative analysis of sex chromosome regions
Three genes from different LGs of cichlids (CLCN5-LG3, OPSIN-LG5, and KCNE1L-LG7) are located on the X chromosome of human and other mammals, highlighting the conservation of this sex chromosome through mammals . However, in bird clade these genes do not correspond to the Z or W chromosomes and, instead of this, are located in autosomes. This is explained by the fact that the ZW and XY chromosomes have no homology and both sexual chromosome pair systems were derived from different autosomes from their common ancestors .
In contrast to mammals and birds, the pufferfish, like most fish, does not possess heteromorphic sex chromosomes, and the genetic mechanism(s) of sex determination is still unclear. However, the human X is a mosaic of orthologs from three chromosomes of the pufferfish: most human Xp and Xq genes are syntenic on pufferfish chromosomes 1, 2, and 7 . In addition, zebrafish LGs 9 and 23 are also related with both human Xp and Xq orthologs [18, 76].
The largest chromosome of Nile tilapia that contains the LG3 was previously thought to be the sex chromosome of this species [77–80]. However, the major sex-determining region in the Nile tilapia was mapped on LG1  located in a small t/a chromosome [present work, 7]. LG1 also contains the WT1 and CYP19a genes involved in mammalian sex differentiation , although they are not considered anymore to be candidate genes for sex determination in African cichlids . Even though our results are still limited, the integrated comparative analysis approach seems to be promising in the clarification of the complex evolutionary dynamics of sex chromosomes among fishes.
Although African cichlids present karyotype variations related to number and morphology of the chromosomes, the linkage groups that were currently investigated (LG1, 3, 5, and 7) were preserved from major changes during their chromosomal diversification. The linkage of large chromosome blocks among cichlids was also preserved in other fish and vertebrates.
The use of BACs containing genes/markers represents a promising alternative for a better physical mapping for cichlids, and the integration of BAC-FISH maps to genomic data stands for a powerful tool to support a better assembling of genomes, contributing in the establishment of a framework for comparative genome-wide studies.
Animals and sampling
Cichlids from Lake Malawi were collected from the wild from 2005–2008 and maintained in the Tropical Aquaculture Facility (TAF) of the University of Maryland (UMD), College Park, MD, USA. Additional African and Asian species of uncertain origin were obtained from commercial sources in Botucatu, SP, Brazil, and South American species were collected from the wild in Brazilian rivers (Table 1) and maintained in the fish room of the Laboratory of Genômica Integrativa (FR-LGI) at Sao Paulo State University (UNESP), Botucatu, SP, Brazil. All the examined specimens were fixed in formaldehyde and then stored in alcohol in the fish collections of TAF-UMD and FR-LGI.
BAC clones and probes labeling
BAC clones containing specific markers of LG1, 3, 5, and 7 (Table 3) were obtained from a BAC library of the Nile tilapia, O. niloticus, which was previously developed , and were used as probes for FISH. BAC extraction was conducted using the PhasePrep®TM BAC DNA Kit (Sigma-Aldrich, St Louis, MO, USA) according to supplier’s protocol. The BAC clones were labeled with biotin, digoxigenin coupled nucleotides (Roche Applied Sciences, Indianapolis, IN, USA), CY3- and CY5-avidin (GE-Healthcare, UK) using whole genome amplification (WGA2 & 3 kits) (Sigma-Aldrich), according to the supplier’s protocol. After that, DNase-I (Sigma-Aldrich) concentration was titrated to yield labeled DNA fragments ranging from 100 to 500 base pair products. For multicolor FISH, we used 16 μl of probe mixture containing: 10 μl of hybridization mixture (62.5% deionized formamide, 12.5% of 20XSSC, 12.5% of dextran sulfate 50%), 4 μg of blocking DNA (salmon/herring sperm DNA or Cot-1 DNA resuspended in the hybridization mixture) and 100 nanograms of each probe. The probe mixture was denatured for 10 min at 65°C and immediately cooled on ice.
Chromosome preparation and FISH procedure
Mitotic chromosomes of cichlid species belonging to the Pseudocrenilabrinae (including representatives of tilapiine, haplochromine, and hemichromine groups), Cichlinae and Etroplinae subfamilies (Table 1) were prepared from anterior kidney cells with in vivo colchicine treatment . The slides with chromosomes were air-dried, treated with pepsin (0.01% in 10 mM HCl) and dehydrated in an ethanol series one day before use. The slides were denatured in 70% formamide/2xSSC, pH 7.0 for 40 s, and dehydrated in an ice-cold ethanol series. The probe mixture was hybridized under a 24 x 50 mm cover slip in a 37°C moist chamber for 48 h. Slides were washed two times for 5 min each in 50% formamide/2xSSC, pH 7.0 at 43°C under agitation, then 10 min in 2xSSC, pH 7.0 at 42°C under continuous agitation. For undirected labeled probes, the hybridization signals were detected with avidin-FITC and rhodamine-anti-DIG (Roche Applied Sciences, Indianapolis, IN, USA), according to the supplier’s protocol. After three washes of 2 min in phosphate buffer detergent (4xSSC/1% Tween-20), slides were mounted with antifade solution containing DAPI (Vectashield mounting medium). Results were recorded with an Olympus BX61 microscope equipped with an Olympus digital camera DP71 and the software Image-Pro MC 6.0.
The selected BACs were checked in order to identify if the FISH mapping results were in agreement to the anchoring marker positions on LGs. Each group of BACs from a particular linkage group was checked to ensure if it hybridizes to the same chromosome using at least two markers from each LG.
Comparative genomic database analyses
The comparative analyses among O. niloticus BAC-end sequences and Oryzias latipes (medaka), Tetraodon nigroviridis (pufferfish) and Gasterosteus aculeatus (stickleback) were done based on the BAC ID (Table 3) of the clones used for FISH against the updated version of the comparative genome browsers of the BouillaBase database . The BAC ID was the starting point to detect a landmark or region of similarity in the databank source. The identification of genomic positions of genes located in LG1, 3, 5, and 7 of O. niloticus in other vertebrates was determined using the currently available public genomic databases Sanger Institute Ensembl Database , Genomicus genome browser , and NCBI Map Viewer . Once a query gene name has been entered in the search box, the position of the gene in all species that have their genome available and annotated was retrieved. The accession numbers of all genes used in the comparative analyze are available as Additional file 2.
JM carried out chromosome preparations, molecular cytogenetics experiments with BAC probes, in silico analyses and data interpretation , and drafted the manuscript. FY conceived the initial chromosome hybridization experiments with BAC probes, analyses and data interpretation and revised the manuscript. TDK helped in the analyses and data interpretation and revised the manuscript. CM conceived the study, participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.
ATPase, Na+/K + transporting
Bacterial artificial chromosomes
Chloride channel 5
v-ski sarcoma viral oncogene homolog
Cytochrome P450, family 19, subfamily A, polypeptide 1
Dual specificity phosphatase 6
Expressed sequence tags
Fluorescence in situ hybridization
Insulin-like growth factor binding protein 2
Potassium voltage-gated channel, Isk-related family, member 1-like
Opsin – cone pigments
RAS-like, estrogen-regulated, growth inhibitor
Sex determining region Y
Wilms tumor 1.
Stiglec R, Ezaz T, Graves JAM: Reassignment of chicken W chromosome sequences to the Z chromosome by fluorescence in situ hybridization (FISH). Cytogenet Genome Res. 2007, 116: 132-134. 10.1159/000097432.
Goldammer T, Brunner RM, Rebl A, Wu CH, Nomura K, Hadfield T, Maddox JF, Cockett NE: Cytogenetic anchoring of radiation hybrid and virtual maps of sheep chromosome X and comparison of X chromosomes in sheep, cattle, and human. Chromosome Res. 2009, 17: 497-506. 10.1007/s10577-009-9047-9.
Katagiri T, Kidd C, Tomasino E, Davis JT, Wishon C, Stern JE, Carleton KL, Howe AE, Kocher TD: A BAC-based physical map of the Nile tilapia genome. BMC Genomics. 2005, 89: 1-6.
NCBI: National Center for Biotechnology Information databases.Http://www.ncbi.nlm.nih.gov/genbank/,
Sahara K, Yoshido A, Marec F, Fukova I, Zhang HB, Wu CC, Goldsmith MR, Yasukochi Y: Conserved synteny of genes between chromosome 15 of Bombyx mori and a chromosome of Manduca sexta shown by five-color BAC-FISH. Genome. 2007, 50: 1061-1065. 10.1139/G07-082.
Freeman JL, Adeniyi A, Banerjee R, Dallaire S, Maguire SF, Chi J, Ng BL, Zepeda C, Scott CE, Humphray S, Rogers J, Zhou Y, Zon LI, Carter NP, Yang F, Lee C: Definition of the zebrafish genome using flow cytometry and cytogenetic mapping. BMC Genomics. 2007, 8: 195-10.1186/1471-2164-8-195.
Cnaani A, Lee B-Y, Zilberman N, Ozouf-Costaz C, Hulata G, Ron M, D’Hont A, Baroiller J-F, D’Cotta H, Penman DJ, Tomasino E, Coutanceau J-P, Pepey E, Shirak A, Kocher TD: Genetics of sex determination in Tilapiine species. Sex Dev. 2008, 1: 43-54.
Phillips RB, Keatley KA, Morasch MR, Ventura AB, Lubieniecki KP, Koop BF, Danzmann RG, Davidson WS: Assignment of Atlantic salmon (Salmo salar) linkage groups to specific chromosomes: Conservation of large syntenic blocks corresponding to whole chromosome arms in rainbow trout (Oncorhynchus mykiss). BMC Genet. 2009, 10: 46-
Ellegren H: Evolutionary stasis: the stable chromosomes of birds. Trends Ecol Evol. 2010, 25: 283-291. 10.1016/j.tree.2009.12.004.
Mota-Velasco JC, Ferreira IA, Cioffi MB, Ocalewicz K, Campos-Ramos R, Shirak A, Lee BY, Martins C, Penman DJ: Characterization of the chromosome fusions in Oreochromis karongae. Chromosome Res. 2010, 18: 575-586. 10.1007/s10577-010-9141-z.
Di Meo GP, Goldammer T, Perucatti A, Genualdo V, Iannuzzi A, Incarnato D, Rebl A, Di Berardino D, Iannuzzi L: Extended cytogenetic maps of sheep chromosome 1 and their cattle and river buffalo homoeologues: comparison with the oar1 RH map and human chromosomes 2, 3, 21 and 1q. Cytogenet Genome Res. 2011, 133: 16-24. 10.1159/000323796.
Mazzuchelli J, Yang F, Kocher TD, Martins C: Comparative cytogenetic mapping of Sox2 and Sox14 in cichlid fishes and inferences on the genomic organization of both genes in vertebrates. Chromosome Res. 2011, 19: 657-667. 10.1007/s10577-011-9225-4.
Sparks JS, Smith WL: Phylogeny and biogeography of cichlid fishes (Teleostei, Perciformes, Cichlidae). Cladistics. 2004, 20: 501-517. 10.1111/j.1096-0031.2004.00038.x.
Smith WL, Chakrabarty P, Sparks JS: Phylogeny, taxonomy, and evolution of Neotropical cichlids (Teleostei: Cichlidae: Cichlinae). Cladistics. 2008, 24: 625-641. 10.1111/j.1096-0031.2008.00210.x.
Lowe-McConnell RH: Ecology of cichlids in South American and African waters, excluding the African Great Lakes. Cichlid fishes: behaviour, ecology and evolution. Edited by: Keenleyside MHA. 1991, London: Chapman and Hall, 60-85.
Kocher TD: Adaptive evolution and explosive speciation: the cichlid fish model. Nat Rev Genet. 2004, 5: 288-298. 10.1038/nrg1316.
Poletto AB, Ferreira IA, Cabral-de-Mello DC, Nakajima RT, Mazzuchelli J, Ribeiro HB, Venere PC, Nirchio M, Kocher TD, Martins C: Chromosome differentiation patterns during cichlid fish evolution. BMC Genet. 2010, 11: 50-
Woods IG, Kelly PD, Chu F, Ngo-Hazelett P, Yan YL, Huang H, Postlethwait JH, Talbot WS: A comparative map of the zebrafish genome. Genome Res. 2000, 10: 1903-1914. 10.1101/gr.10.12.1903.
Naruse K, Fukamachi S, Mitani H, Kondo M, Matsuoka T, Kondo S, Hanamura N, Morita Y, Hasegawa K, Nishigaki R, Shimada A, Wada H, Kusakabe T, Suzuki N, Kinoshita M, Kanamori A, Terado T, Kimura H, Nonaka M, Shima A: A detailed linkage map of medaka, Oryzias latipes: comparative genomics and genome evolution. Genetics. 2000, 154: 1773-1784.
Waldbieser GC, Bosworth BG, Nonneman DJ, Wolters WR: A microsatellite based genetic linkage map for channel catfish, Ictalurus punctatus. Genetics. 2000, 158: 727-734.
Young WP, Wheeler PA, Coryell VH, Keim P, Thorgaard GH: A detailed linkage map of rainbow trout produced using doubled haploids. Genetics. 1998, 148: 839-850.
Moen T, Hoyheim B, Munck H, Gomez-Raya L: A linkage map of Atlantic salmon (Salmo salar) reveals an uncommonly large difference in recombination rate between the sexes. Anim Genet. 2004, 35: 81-92. 10.1111/j.1365-2052.2004.01097.x.
Gilbey J, Verspoor E, McLay A, Houlihan D: A microsatellite linkage map for Atlantic salmon (Salmo salar). Anim Genet. 2004, 35: 98-105. 10.1111/j.1365-2052.2004.01091.x.
Lee BY, Lee WJ, Streelman JT, Carleton KL, Howe AE, Hulata G, Slettan A, Stern JE, Terai Y, Kocher TD: A second-generation genetic linkage map of tilapia (Oreochromis spp.). Genetics. 2005, 170: 237-244. 10.1534/genetics.104.035022.
Albertson RC, Streelman JT, Kocher TD: Directional selection has shaped the oral jaws of Lake Malawi cichlid fishes. Proc Natl Acad Sci USA. 2003, 100: 5252-5257. 10.1073/pnas.0930235100.
Sanetra M, Henning F, Fukamachi S, Meyer A: A microsatellite-based genetic linkage map of the cichlid fish, Astatotilapia burtoni (Teleostei): a comparison of genomic architectures among rapidly speciating cichlids. Genetics. 2009, 182: 387-97. 10.1534/genetics.108.089367.
Watanabe M, Kobayashi N, Shin-i T, Horiike T, Tateno Y, Kohara Y, Okada N: Extensive analysis of ORF sequences from two different cichlid species in Lake Victoria provides molecular evidence for a recent radiation event of the Victoria species flock: identity of EST sequences between Haplochromis chilotes and Haplochromis sp. “Redtailsheller”. Gene. 2004, 343: 263-269. 10.1016/j.gene.2004.09.013.
Kobayashi N, Watanabe M, Horiike T, Kohara Y, Okada N: Extensive analysis of EST sequences reveals that all cichlid species in Lake Victoria share almost identical transcript sets. Gene. 2009, 441: 187-191. 10.1016/j.gene.2008.11.023.
Renn SC, Aubin-Horth N, Hofmann HA: Biologically meaningful expression profiling across species using heterologous hybridization to a cDNA microarray. BMC Genomics. 2004, 5: 42-10.1186/1471-2164-5-42.
Salzburger W, Renn SC, Steinke D, Braasch I, Hofmann HA, Meyer A: Annotation of expressed sequence tags for the East African cichlid fish Astatotilapia burtoni and evolutionary analyses of cichlid ORFs. BMC Genomics. 2008, 9: 96-10.1186/1471-2164-9-96.
Lee BY, Howe AE, Conte MA, D’Cotta H, Pepey E, Baroiller JF, di Palma F, Carleton KL, Kocher TD: An EST resource for tilapia based on 17 normalized libraries and assembly of 116,899 sequence tags. BMC Genomics. 2010, 11: 278-10.1186/1471-2164-11-278.
Katagiri T, Asakawa S, Minagawa S, Shimizu N, Hirono I, Aoki T: Construction and characterization of BAC libraries for three fish species; rainbow trout, carp and tilapia. Anim Genet. 2001, 32: 200-204. 10.1046/j.1365-2052.2001.00764.x.
Di Palma F, Kidd C, Borowsky R, Kocher TD: Construction of bacterial artificial chromosome libraries for the Lake Malawi cichlid (Metriaclima zebra), and the blind cavefish (Astyanax mexicanus). Zebrafish. 2007, 4: 41-47. 10.1089/zeb.2006.9996.
Watanabe M, Kobayashi N, Fujiyama A, Okada N: Construction of a BAC library for Haplochromis chilotes, a cichlid fish from Lake Victoria. Genes Genet Syst. 2003, 78: 103-105. 10.1266/ggs.78.103.
Lang M, Miyake T, Braasch I, Tinnemore D, Siegel N, Salzburger W, Amemiya CT, Meyer A: A BAC library of the East African haplochromine cichlid fish Astatotilapia burtoni. J Exp Zool B Mol Dev Evol. 2006, 306: 35-44.
Guyon R, Rakotomanga M, Azzouzi N, Coutanceau J-P, Bonillo C, D'Cotta H, Pepey E, Soler L, Rodier-Goud M, D'Hont A, Conte MA, van Bers NEM, Penman DJ, Hitte C, Crooijmans RPMA, Kocher TDK, Ozouf-Costaz C, Baroiller J-F, Galibert F: A high-resolution map of the Nile tilapia genome: a resource for studying cichlids and other percomorphs. BMC Genomics. 2012, 13: 222-10.1186/1471-2164-13-222.
Ezaz MT, Harvey SC, Boonphakdee C, Teale AJ, McAndrew BJ, Penman DJ: Isolation and physical mapping of sex-linked AFLP markers in Nile tilapia (Oreochromis niloticus L.). Mar Biotechnol. 2004, 6: 435-445. 10.1007/s10126-004-3004-6.
Lee BY, Kocher TD: Exclusion of Wilms tumor (WT1b) and ovarian cytochrome P450 aromatase (CYP19A1) as candidates for sex determination genes in Nile tilapia (Oreochromis niloticus). Anim Genet. 2007, 38: 85-86. 10.1111/j.1365-2052.2007.01563.x.
Lee BY, Coutanceau JP, Ozouf-Costaz C, D’Cotta H, Baroiller JF, Kocher T: Genetic and physical mapping of sex linked AFLP markers in Nile tilapia (Oreochromis niloticus). Mar Biotechnol. 2011, 13: 557-562. 10.1007/s10126-010-9326-7.
Cnaani A, Kocher TD: Sex-linked markers and microsatellite locus duplication in the cichlid species Oreochromis tanganicae. Biol Lett. 2008, 4: 700-703. 10.1098/rsbl.2008.0286.
Ser JR, Roberts RB, Kocher TD: Multiple interacting loci control sex determination in Lake Malawi Cichlid fish. Evolution. 2009, 64: 486-501.
Hossain A, Saunders GF: The human sex-determining gene SRY is a direct target of WT1*. J Biol Chem. 2001, 276: 16817-16823. 10.1074/jbc.M009056200.
Chang X, Kobayashi T, Senthilkumaran B, Kobayashi-Kajura H, Sudhakumari CC, Nagahama Y: Two types of aromatase with different encoding genes, tissue distribution and developmental expression in Nile tilapia (Oreochromis niloticus). Gen Comp Endocrinol. 2005, 141: 101-15. 10.1016/j.ygcen.2004.11.020.
Levtchenko EN, Monnens LA, Bökenkamp A, Knoers NV: From gene to disease; Dent's disease caused by abnormalities in the CLCN5 and OCRL1 genes. Ned Tijdschr Geneeskd. 2007, 151: 2377-2380.
Carleton KL, Hofmann CM, Klisz C, Patel Z, Chircus LM, Simenauer LH, Soodoo N, Albertson RC, Ser JR: Genetic basis of differential opsin gene expression in cichlid fishes. J Evol Biol. 2010, 23: 840-853. 10.1111/j.1420-9101.2010.01954.x.
Comparative Genome Browsers.Http://www.BouillaBase.org,
Feldberg E, Porto JIR, Bertollo LAC: Chromosomal changes and adaptation of cichlid fishes during evolution. Fish adaptations. Edited by: Val AL, Kapoor BG. 2003, New Dehli and New York: Science Publishers, 285-308.
Turner GF: Adaptive radiation of cichlid fish. Curr Biol. 2007, 17: 827-831. 10.1016/j.cub.2007.07.026.
Gross MC, Schneider CH, Valente GT, Porto JI, Martins C, Feldberg E: Comparative cytogenetic analysis of the genus Symphysodon (discus fishes, cichlidae): chromosomal characteristics of retrotransposons and minor ribosomal DNA. Cytogenet Genome Res. 2009, 127: 43-53. 10.1159/000279443.
Gross MC, Schneider CH, Valente GT, Martins C, Feldberg E: Variability of 18S rDNA locus among Symphysodon fishes: chromosomal rearrangements. J Fish Biol. 2010, 76: 1117-1127. 10.1111/j.1095-8649.2010.02550.x.
Majumdar KC, McAndrew BJ: Relative DNA content of somatic nuclei and chromosomal studies in three genera, Tilapia, Sarotherodon and Oreochromis of the tribe Tilapiini (Pisces, Cichlidae). Genetica. 1986, 68: 175-188. 10.1007/BF02424441.
Ferreira IA, Poletto AB, Mota-Velasco J, Penman DJ, Kocher TD, Martins C: Chromosome evolution in African cichlid fish: contributions from the physical mapping of repeated DNAs. Cytogenet Genome Res. 2010, 129: 314-22. 10.1159/000315895.
Chew JSK, Oliveira C, Wright JM, Dobson MJ: Molecular and cytogenetic analysis of the telomeric (TTAGGG)n repetitive sequences in the Nile tilapia, Oreochromis niloticus (Teleostei: Cichlidae). Chromosoma. 2002, 111: 45-52. 10.1007/s00412-002-0187-3.
Oliveira C, Chew JSK, Foresti F, Dobson M, Wright JM: A LINE2 repetitive DNA sequence from the cichlid fish, Oreochromis niloticus: sequence analysis and chromosomal distribution. Chromosoma. 1999, 108: 457-468. 10.1007/s004120050397.
Oliveira C, Wang Y, Bryden J, Wright JM: Short interspersed repetitive elements (SINEs) from the cichlid fish, Oreochromis niloticus, and their chromosomal localization by fluorescent in situ hybridization. Caryologia. 2003, 56: 181-189.
Martins C, Oliveira C, Wasko AP, Wright JM: Physical mapping of the Nile tilapia (Oreochromis niloticus) genome by fluorescent in situ hybridization of repetitive DNAs to metaphase chromosomes - A review. Aquaculture. 2004, 231: 37-49. 10.1016/j.aquaculture.2003.08.017.
Ferreira IA, Martins C: Physical chromosome mapping of the Nile tilapia Oreochromis niloticus using repetitive DNA sequences. Micron. 2008, 39: 411-418. 10.1016/j.micron.2007.02.010.
Valente GT, Schneider CH, Gross MC, Feldberg E, Martins C: Comparative cytogenetics of cichlid fishes through genomic in-situ hybridization (GISH) with emphasis on Oreochromis niloticus. Chromosome Res. 2009, 17: 791-799. 10.1007/s10577-009-9067-5.
Valente GT, Mazzuchelli J, Ferreira IA, Poletto AB, Fantinatti BEA, Martins C: Cytogenetic mapping of the retroelements Rex1, Rex3 and Rex6 among cichlid fish: new insights on the chromosomal distribution of transposable elements. Cytogenet Genome Res. 2011, 133: 34-42. 10.1159/000322888.
Thompson KW: Karyotypes of six species of African Cichlidae (Pisces: Perciformes). Experientia. 1981, 37: 351-352. 10.1007/BF01959857.
Harvey SC, Campos-Ramos R, Kennedy DD, Ezaz MT, Bromage NR, Griffin DK, Penman DJ: Karyotype evolution in Tilapia: mitotic and meiotic chromosome analysis of Oreochromis karongae and O. niloticus × O. karongae hybrids. Genetica. 2002, 115: 169-177. 10.1023/A:1020190918431.
Farias IP, Ortí G, Sampaio I, Schneider H, Meyer A: Mitochondrial DNA phylogeny of the family Cichlidae: monophyly and fast molecular evolution of the Neotropical assemblage. J Mol Evol. 1999, 48: 703-711. 10.1007/PL00006514.
Soler L, Conte MA, Katagiri T, Howe AE, Lee BY, Amemiya C, Stuart A, Dossat C, Poulain J, Johnson J, Di Palma F, Lindblad-Toh K, Baroiller JF, D'Cotta H, Ozouf-Costaz C, Kocher TD: Comparative physical maps derived from BAC end sequences of tilapia (Oreochromis niloticus). BMC Genomics. 2010, 11: 636-10.1186/1471-2164-11-636.
Sarropoulou E, Nousdili A, Magoulas A, Kotoulas G: Linking the genomes of nonmodel teleosts through comparative genomics. Mar Biotechnol. 2008, 10: 227-233. 10.1007/s10126-007-9066-5.
Becker SED, Thomas R, Trifonov VA, Wayne RK, Graphodatsky AS, Breen M: Anchoring the dog to its relatives reveals new evolutionary breakpoints across 11 species of the Canidae and provides new clues for the role of B chromosomes. Chromosome Res. 2011, 19: 685-708. 10.1007/s10577-011-9233-4.
Giovannotti M, Caputo V, O’Brien PCM, Lovell FL, Trifonov V, Nisi Cerioni P, Olmo E, Ferguson-Smith MA, Rens W: Skinks (Reptilia: Scincidae) have highly conserved karyotypes as revealed by chromosome painting. Cytogenet Genome Res. 2009, 127: 224-231. 10.1159/000295002.
Srikulnath K, Nishida C, Matsubara K, Uno Y, Thongpan A, Suputtitada S, Apisitwanich S, Matsuda Y: Karyotypic evolution in squamate reptiles: comparative gene mapping revealed highly conserved linkage homology between the butterfly lizard (Leiolepis reevesii rubritaeniata, Agamidae, Lacertilia) and the Japanese four-striped rat snake (Elaphe quadrivirgata, Colubridae, Serpentes). Chromosome Res. 2009, 17: 975-986. 10.1007/s10577-009-9101-7.
Lee BY, Kocher TD: Comparative genomics and positional cloning. Aquaculture Genome Technologies. Edited by: Liu J. 2007, Ames, IA: Blackwell, 325-337.
Grützner F, Crollius HR, Lütjens G, Jaillon O, Weissenbach J, Ropers HH, Haaf T: Four-hundred million years of conserved synteny of human Xp and Xq genes on three Tetraodon chromosomes. Genome Res. 2002, 12: 1316-1322. 10.1101/gr.222402.
Robinson-Rechavi M, Marchand O, Escriva H, Bardet P, Zelus D, Hughes S, Laudet V: Euteleost fish genomes are characterized by expansion of gene families. Genome Res. 2001, 11: 781-788. 10.1101/gr.165601.
Robinson-Rechavi M, Marchand O, Escriva H, Laudet V: An ancient whole-genome duplication may not have been responsible for the abundance of duplicated fish genes. Curr Biol. 2001, 11: 458-459. 10.1016/S0960-9822(01)00280-9.
Venkatesh B: Evolution and diversity of fish genomes. Curr Opin Genet Dev. 2003, 13: 588-592. 10.1016/j.gde.2003.09.001.
Ravi V, Venkatesh B: Rapidly evolving fish genomes and teleost diversity. Curr Opin Genet Dev. 2008, 18: 544-50. 10.1016/j.gde.2008.11.001.
Murphy WJ, Sun S, Chen ZQ, Pecon-Slattery J, O'Brien SJ: Extensive conservation of sex chromosome organization between cat and human revealed by parallel radiation hybrid mapping. Genome Res. 1999, 9: 1223-1230. 10.1101/gr.9.12.1223.
Kawai A, Nishida-Umehara C, Ishijima J, Tsuda Y, Ota H, Matsuda Y: Different origins of bird and reptile sex chromosomes inferred from comparative mapping of chicken Z linked genes. Cytogenet Genome Res. 2007, 117: 92-102. 10.1159/000103169.
Barbazuk WB, Korf I, Kadavi C, Heyen J, Tate S, Wun E, Bedell JA, McPherson JD, Johnson SL: Syntenic relationship of the zebrafish and human genomes. Genome Res. 2000, 10: 1351-1358. 10.1101/gr.144700.
Foresti F, Oliveira C, Galetti PM, de Almeida-Toledo LF: Synaptonemal complex analysis in spermatocytes of tilapia, Oreochromis niloticus (Pisces, Cichlidae). Genome. 1993, 36: 1124-1128. 10.1139/g93-150.
Carrasco LAP, Penman DJ, Bromage N: Evidence for the presence of sex chromosomes in the Nile tilapia (Oreochromis niloticus) from synaptonemal complex analysis of XX, XY and YY genotypes. Aquaculture. 1999, 173: 207-218. 10.1016/S0044-8486(98)00488-8.
Harvey SC, Boonphakdee C, Campos-Ramos R, Ezaz MT, Griffin DK, Bromage NR, Penman P: Analysis of repetitive DNA sequences in the sex chromosomes of Oreochromis niloticus. Cytogenet Genome Res. 2003, 101: 314-319. 10.1159/000074355.
Ocalewicz K, Mota-Velasco JC, Campos-Ramos R, Penman DJ: FISH and DAPI staining of the synaptonemal complex of the Nile tilapia (Oreochromis niloticus) allow orientation of the unpaired region of bivalent 1 observed during early pachytene. Chromosome Res. 2009, 17: 773-782. 10.1007/s10577-009-9071-9.
Lee BY, Penman DJ, Kocher TD: Identificationof a sex-determining region in Nile tilapia (Oreochromis niloticus) using bulked segregant analysis. Anim Genet. 2003, 34: 379-383. 10.1046/j.1365-2052.2003.01035.x.
Bertollo LAC, Takahashi CS, Moreira-Filho O: Cytotaxonomic consideration on Hoplias lacerdae (Pisces, Erythrinidae). Braz J Genet. 1978, 1: 103-120.
Sanger Institute Ensembl Database.http://www.ensembl.org,
Genomicus genome browser.http://www.dyogen.ens.fr/genomicus,
NCBI Map Viewer.http://www.ncbi.nlm.nih.gov/genomes,
This study was supported by Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), Coordenadoria de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) from Brazil. FY was supported by a Wellcome Trust grant (WT098051). The authors are grateful to Beiyuan Fu of Wellcome Trust Sanger Institute for technical assistant in BAC-FISH, and AP Wasko of Sao Paulo State University for the careful revision of the manuscript.
The authors declare that they have no competing interests.