Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Chromosome map of the Siamese cobra: did partial synteny of sex chromosomes in the amniote represent “a hypothetical ancestral super-sex chromosome” or random distribution?

BMC Genomics volume 19, Article number: 939 (2018) Cite this article

Abstract

Background

Unlike the chromosome constitution of most snakes (2n=36), the cobra karyotype shows a diploid chromosome number of 38 with a highly heterochromatic W chromosome and a large morphologically different chromosome 2. To investigate the process of sex chromosome differentiation and evolution between cobras, most snakes, and other amniotes, we constructed a chromosome map of the Siamese cobra (Naja kaouthia) with 43 bacterial artificial chromosomes (BACs) derived from the chicken and zebra finch libraries using the fluorescence in situ hybridization (FISH) technique, and compared it with those of the chicken, the zebra finch, and other amniotes.

Results

We produced a detailed chromosome map of the Siamese cobra genome, focusing on chromosome 2 and sex chromosomes. Synteny of the Siamese cobra chromosome 2 (NKA2) and NKAZ were highly conserved among snakes and other squamate reptiles, except for intrachromosomal rearrangements occurring in NKA2. Interestingly, twelve BACs that had partial homology with sex chromosomes of several amniotes were mapped on the heterochromatic NKAW as hybridization signals such as repeat sequences. Sequence analysis showed that most of these BACs contained high proportions of transposable elements. In addition, hybridization signals of telomeric repeat (TTAGGG)n and six microsatellite repeat motifs ((AAGG)8, (AGAT)8, (AAAC)8, (ACAG)8, (AATC)8, and (AAAAT)6) were observed on NKAW, and most of these were also found on other amniote sex chromosomes.

Conclusions

The frequent amplification of repeats might involve heterochromatinization and promote sex chromosome differentiation in the Siamese cobra W sex chromosome. Repeat sequences are also shared among amniote sex chromosomes, which supports the hypothesis of an ancestral super-sex chromosome with overlaps of partial syntenies. Alternatively, amplification of microsatellite repeat motifs could have occurred independently in each lineage, representing convergent sex chromosomal differentiation among amniote sex chromosomes.

Background

Sex chromosomes evolved independently from a homologous autosomal pair in both plants and animals when one chromosomal partner acquired a sex-determining allele [1, 2]. Differentiated sex chromosomes then developed from the accumulation of sexually antagonistic alleles at loci linked to the sex-specific region. To maintain linkage of these genes, meiotic recombination between proto-sex chromosomes was suppressed around the heterologous region in one sex, which in turn promoted the loss of active genes or segmental insertions and deletions, and caused this region to extend along part or all of the chromosome. This produced chromosomal morphological differences between the X and Y (male heterogamety) or Z and W (female heterogamety) sex chromosomes. It has also been proposed that structural changes in sex chromosomes might expedite the suppression of recombination, favoring the further amplification of repeat sequences and leading to heterochromatinization [3]. It is thus of great importance to understand the mechanism and role of sex chromosome differentiation with regard to sex determination.

Sex chromosomes of amniotes are highly diverse, ranging from cryptic to highly heteromorphic with XY or ZW systems. Using comparative gene mapping (chromosome mapping via a cytogenetic technique) and whole-genome sequencing, genomic convergence was discovered, in which unrelated sex chromosomes share syntenies across distantly related taxa [4]. Highly conserved synteny of the chicken Z chromosome has been found in the X chromosomes of the Mexican giant musk turtle (Staurotypus triporcatus) and the giant musk turtle (S. salvinii) [5]. However, the X chromosomes of the marsh turtle (Siebenrockiella crassicollis) and wood turtle (Glyptemys insculpta) correspond to chicken chromosome 5, while the Z chromosome of the Chinese softshell turtle (Pelodiscus sinensis) and the spiny softshell turtle (Apalone spinifera) corresponds to chicken chromosome 15 [6,7,8,9]. The Hokou gecko (Gekko hokouenesis, Gekkota) Z chromosome is homologous with the chicken Z chromosome [10, 11], whereas the micro-X chromosome of the Anolis lizard (Anolis carolinensis, Iguania), the micro-Z chromosome of the dragon lizard (Pogona vitticeps, Iguania), and the Z chromosome of the sand lizard (Lacerta agilis, Lacertoidea) have homology with chicken chromosomes 15, 17 and 23, and 6 and 9, respectively [4, 12,13,14]. Notably, monotreme sex chromosomes also share partial synteny with the chicken Z chromosome [15], which corresponds partially to squamate chromosome 2. Comparative genomics, based on chromosome painting and gene mapping for several squamate reptiles, have revealed that the chicken Z chromosome is homologous with the short arm of the bi-armed chromosome 2 in most squamate reptiles [11, 13, 14, 16,17,18,19]. This suggests the possibility that squamate reptile chromosome 2 is part of a larger ancestral amniote super-sex chromosome with overlaps of partial syntenies among amniotes. However, the hypothesis of an ancestral super-sex chromosome does not just relate to the homology of chicken Z chromosome, but states that unrelated sex chromosomes share syntenies across distantly related taxa. Ancestral super-sex chromosomes probably exist in amniotes. Multiple translocations have led to the appearance of sauropsid and diapsid sex chromosomes [4].

Snakes (Serpentes), a species-rich lineage of extant reptiles, exhibit conserved ZZ/ZW-type sex chromosomes in most species, except for Boa imperator, Python bivittatus, P. regius [20,21,22,23,24]. The Z chromosome, which is similar in size across species, is the fourth or fifth largest metacentric chromosome, whereas the W chromosome varies from being homomorphic in most henophidians (primitive snakes) to highly differentiated in caenophidians (advanced snakes), as a consequence of different centromere positions and/or amounts of heterochromatin [25,26,27,28,29,30,31]. Large amplification of microsatellite repeat motifs or telomeric repeats have been found in the W chromosome of several colubroid snakes [25, 26], and genomic comparisons have revealed that syntenies with the snake sex chromosome are also present in a large variety of chromosomes in squamate reptiles [4, 11,12,13,14, 16,17,18, 32,33,34]. However, the chicken Z chromosome is homologous with the short arm of snake chromosome 2, as also found in other squamate reptiles, while the snake Z chromosome is homologous with chicken chromosomes 2 and 27 [21,22,23,24, 32,33,34]. Interestingly, the non-homologous W chromosomes of chickens and the common tiger snake (Notechis scutatus) share common repeat sequences that are not present elsewhere in the genomes analyzed to date [30]. This suggests that repeat sequences are partially shared between the sex chromosomes of the chicken and snake, and supports the hypothesis that squamate chromosome 2 is part of an ancestral super-sex chromosome in amniotes [4].

Unlike most snakes, which have a chromosome number of 2n=36, consisting of 16 macrochromosomes and 20 microchromosomes [29, 35], cobra (Naja spp.) karyotypes show 2n=38, comprising five pairs of macrochromosomes and 14 pairs of microchromosomes [35]. In the Siamese cobra (Naja kaouthia), the W chromosome exhibits highly remarkable amplified telomeric repeats, indicating the role of repeat sequences for sex chromosome differentiation. Whereas most snakes have a large submetacentric chromosome 2, which is also highly conserved in other squamate reptiles, cobras have a large subtelocentric chromosome 2 [35]. Examination of cobra karyotypes points to several fusions and fissions in the Naja lineage, suggesting that the subtelocentric chromosome 2 in the Naja lineage appeared by centric fission of a large submetacentric chromosome 2 in ancestral snake karyotypes. By contrast, examples of chromosome 2 resembling those of Naja are not found in the majority of squamate reptiles; however, several examples of lizard subtelocentric/acrocentric chromosome 2 show homology with chromosome 2 of most squamate reptiles [11, 14, 16,17,18,19]. An alternative explanation should be considered, whereby the subtelocentric chromosome 2 appeared by intrachromosomal rearrangements in the Naja lineage, although no reports have confirmed this. With these two possible scenarios involving the W chromosome and structural change in snake chromosome 2, cobras are very good models with which to study the process of sex chromosome differentiation and evolution, including the hypothesis of an ancestral super-sex chromosome in amniotes. In this study, we constructed a comparative chromosome map with 43 bacterial artificial chromosomes (BACs) derived from the chicken and zebra finch genomes using fluorescence in situ hybridization (FISH) [36,37,38], and compared the chromosome homology of the Siamese cobra (N. kaouthia) with those of the Japanese striped snake (Elaphe quadrivirgata) and other amniotes. We also screened the Siamese cobra chromosomes using FISH with 19 microsatellite repeat motifs and determined in silico the copy number of microsatellite repeats on the chicken and zebra finch BAC sequences for comparison. This allowed us to delineate the synteny between species, to investigate its significance, and to discuss the hypothesis of an ancestral super-sex chromosome in amniotes.

Results

Karyotype and Chromosomal Locations of Telomeric Repeat and Microsatellite Repeat Motifs

We examined more than 20 DAPI-stained metaphase spreads from two males and two females. Chromosome numbers were 2n=38, comprising 10 macrochromosomes and 28 microchromosomes in both sexes. The macrochromosomes comprised one pair of large metacentrics (first), one pair of large subtelocentrics (second), two pairs of medium-sized metacentrics (third and fourth), and one pair of submetacentric (fifth). The fourth largest chromosomes were sex chromosomes, in which the female karyotype contained the metacentric Z chromosome and the submetacentric W chromosome (Figs. 1a and b). The DAPI-stained chromosome banding pattern showed large DAPI-negative bands were on the long arm of the W chromosome, indicating GC richness in this region. C-banding revealed that the W chromosome exhibited C-positive heterochromatin in the entire region (Fig. 1c). Hybridization signals of hexamer repeat sequences (TTAGGG)n were found at the ends of all chromosomes. High intensity signals were observed on two microchromosome pairs. Interstitial telomeric sites were found on three macrochromosome pairs in either sex, while two large-sized (TTAGGG)n signals on the q arm of the W chromosome were also observed, but not on the Z chromosome (Fig. 1d). Eleven of nineteen microsatellite repeat motifs were successfully mapped onto most N. kaouthia chromosomes (NKA) (Fig. 2, Additional file 1: Figure S1); however, only six microsatellite repeat motifs showed high amplification on the W chromosome. Strong hybridization signals of tetranucleotide (AAGG)8 were localized to NKAWq in addition to the predominant telomeric sites of most chromosomes. By contrast, (AGAT)8, (AAAC)8, (ACAG)8, (AATC)8, and (AAAAT)6 repeat motifs were specifically mapped to NKAWq (Fig. 2).

Fig. 1
figure 1

Karyotype and metaphase spread of the Siamese cobra (Naja kaouthia). Giemsa- (a), grey images of DAPI-stained karyotypes (b), and C-banded metaphase spread (c) of female Siamese cobra (N. kaouthia). FISH patterns of the telomeric (TTAGGG)n sequence on DAPI-stained metaphase chromosome spreads of female Siamese cobra (N. kaouthia) (d). Arrows indicate the centromeric region of chromosome 2 and arrowheads indicate signals of interstitial telomeric sites. “Z” indicates Z sex chromosome, and “W” indicates W sex chromosome. Scale bar represents 10 μm

Full size image
Fig. 2
figure 2

Chromosomal locations of microsatellite repeat motifs in the Siamese cobra (Naja kaouthia). Hybridization patterns of FITC-labeled (AGAT)8 (a), (AAGG)8 (b), (ACAG)8 (c), (AATC)8 (d), (AAAC)8 (e), and (AAAAT)6 (f) on DAPI-stained chromosomes. Arrows indicate the hybridization signals. “Z” indicates Z sex chromosome, and “W” indicates W sex chromosome. Scale bar represents 10 μm

Full size image

Chromosome Homology among the Siamese cobra, Chicken, and Zebra finch

Twenty-six chicken BACs and seventeen zebra finch BACs were mapped to the Siamese cobra chromosomes. We constructed a chromosome map for the Siamese cobra with 43 BACs (markers) (Table 1, Figs. 3 and 4, Additional file 2: Figure S2). More than 20 metaphase spreads were observed for each BAC, with hybridization efficiencies ranging from approximately 70% to 90%. Chromosome homology between N. kaouthia, chicken and zebra finch was analyzed using the chicken genome database (http://www.ncbi.nlm.nih.gov/genome/guide/chicken/) and zebra finch genome database (http://www.ncbi.nlm.nih.gov/genome/367). Seven BACs mapped to NKA2 were localized to GGA1, GGA18 and GGAZ and TGU12 and TGUZ. Three BACs mapped to NKAZ were localized to GGA2p, GGA27, and TGU27. However, twelve BACs on NKAW were located on GGA1, GGA4, GGA9, GGAZ, TGU1B, TGU9, TGU13, TGU17, TGU23, TGU27, and TGUZ. The remaining 22 BACs on NKA1 and microchromosomes were located on GGA1, GGA4, GGA5, GGA9, GGA15, GGA17, GGA23, TGU4A, TGU5, TGU15 and TGU23 (Table 1).

Table 1 List of chicken and zebra finch BACs mapped to the Siamese cobra chromosomes and their chromosomal location in the chicken (Gallus gallus)
Full size table
Fig. 3
figure 3

Chromosomal locations of chicken and zebra finch BACs in female Siamese cobra (Naja kaouthia). GGAZ BAC (CH261-133M4) (a), GGA18 BAC (CH261-67N15) (b), GGA18 BAC (CH261-72B18) (c), and GGA12 BAC (TGMCBA-305E19) (d) were localized to chromosome 2 (NKA2). Arrows indicate the hybridization signals. “Z” indicates Z sex chromosome, and “W” indicates W sex chromosome. Scale bar represents 10 μm

Full size image
Fig. 4
figure 4

Chromosomal locations of chicken and zebra finch BAC clones in female Siamese cobra (Naja kaouthia). GGA2p BAC (CH261-177K1) (a), TGU27 BAC (TGMCBA-23C5) (b), GGA9 BAC (CH261-95N3) (c), and GGAZ BAC (CH261-129A16) (d) were localized to W chromosome (NKAW). Arrows indicate the hybridization signals. “Z” indicates Z sex chromosome, and “W” indicates W sex chromosome. Scale bar represents 10 μm

Full size image

Sequence Analysis of Chicken and Zebra finch BACs Located on the Siamese cobra W sex chromosomes

The content of repeat sequences in chicken and zebra finch BACs mapped to NKAW was examined using Repbase-GIRI and MSATCOMMANDER (Tables 2 and 3). All BACs contained 3.04–24.85% repeat sequences, of which 0.07–1.26% were transposable elements. The BACs mapped to NKAW tended to have twice as many repeat sequences as BACs mapped to NKAZ, NKA2, or other chromosomes (Additional files 3, 4, 5, 6, 7 and 8: Tables S1, S2, S3, S4, S5, and S6), showing no preferential accumulation of any specific repeat sequences. Our analysis identified several microsatellite repeat motifs within BACs mapped to NKAW, most of which were not the same as those found using microsatellite repeat motif mapping (Table 3, Fig. 2). Additionally, BLAST analysis of BAC sequences mapped to NKAW revealed the predictions of five functional genes, including DDX26B, DHX36, FHL1, GPR149, and MAP7D3, which were related to sex determination and sex development in vertebrates (Table 4).

Table 2 Comparison of major classes of repeat sequences in chicken and zebra finch BACs mapped on the Siamese cobra W chromosome
Full size table
Table 3 Comparison of frequencies of microsatellite repeat motifs in chicken and zebra finch BACs mapped on the Siamese cobra W chromosome
Full size table
Table 4 Gene list relating to sex determination and sex developmental pathways for chicken and zebra finch BACs mapped on the Siamese cobra W chromosome
Full size table

Discussion

Comparison of the draft genome assembly and chromosome maps among amniotes revealed a high level of conservation for synteny and allows us to deduce the process of chromosomal rearrangement over millions of years [11,12,13,14, 16,17,18, 21, 33, 34]. However, sex chromosomes are a more dynamic entity in the genome than autosomes [3, 39]. The W chromosome of snakes is the most variable element in the snake genome, with differences in the number of genes and amplification of several repeat motifs, even between closely related species [21, 22, 25, 26]. This suggests that molecular differentiation has often occurred on the W chromosome in snakes, after divergence from a common ancestor. Amplification of repeat sequences might also have resulted in the large size of the W chromosome found in the tiger snake (N. scutatus) and the tiger keelback snake (Rhabdophis tigrinus), arising from a recent addition of repeats [26, 30, 40]. Conversely, the small size of the W chromosome in the Japanese striped snake (E. quadrivirgata) may be a consequence of degradation [21].

Amplification of Microsatellite Repeat Motifs on the Siamese cobra W Chromosome

The W chromosome was smaller than the Z chromosome in the Siamese cobra, while C-banding revealed almost the entire heterochromatic nature of the W chromosome. This suggests that size reduction and heterochromatinization are correlated with W chromosome differentiation in the Siamese cobra. Our FISH analyses revealed that six microsatellite repeat motifs ((AAGG)8, (AGAT)8, (AAAC)8, (ACAG)8, (AATC)8, and (AAAAT)6) were amplified on NKAWq. Hybridization signals of six microsatellite repeat motifs also overlapped with telomeric (TTAGGG)n repeats identified on the NKAWq. These results collectively suggest that the NKAW are structural complexes containing various repeat sequences on the female-specific region. A different hybridization pattern among repeat motifs indicates that each repeat type is distinct. Such amplification of microsatellite repeats and telomeric repeats in the heterochromatic region on W chromosomes has been reported in several caenophidian snakes and other squamate reptiles [25,26,27, 41] (Fig. 5). One microsatellite repeat motif amplified on the W chromosome in several caenophidian snakes though not for basal lineage is the banded krait minor satellite (Bkm), consisting of microsatellite repeat motif (AGAT)n or (GACA)n sequences, and related to the degree of ZW differentiation [25, 26, 42, 43]. This also suggests that frequent amplification of the repeats has a structural role in heterochromatinization and promotes further sex chromosome differentiation.

Fig. 5
figure 5

Schematic representation of microsatellite repeat motif amplification on the Siamese cobra W chromosome and other amniote sex chromosomes sharing the same motifs. Phylogeny was partially derived from Pyron et al. [43]. Regions where the microsatellite repeat motifs were hybridized are shown to the right of the chromosomes. Chromosome numbers indicate the pink-tailed worm-lizard (Aprasia parapulchella, APA), the Japanese striped snake (Elaphe quadrivirgata, EQU), the tiger keelback snake (Rhabdophis tigrinus, RTI), the habu Protobothrops flavoviridis, (PFL), the Japanese moccasin Gloydius blomhoffii, (GBL), the eastern three-lined skink (Bassiana duperreyi, BDU), the central bearded dragon (Pogona vitticeps, PVI), the eastern long-necked turtle (Chelodina longicollis, CLO), the tiger snake (Notechis scutatus, NSC), the tiger keelback snake (Rhabdophis tigrinus, RTI), the dragon snake (Xenodermus javanicus, XJA), the South American rattlesnake (Crotalus durissus, CDU), the masked water snake (Homalopsis buccata, HBU), the African house snake (Boaedon fuliginosus, BFU and Boaedon olivaceus, BOL), the milk snake (Lampropeltis triangulum, LTR), the Ruthven's kingsnake (Lampropeltis ruthveni, LRU), the corn snake (Pantherophis guttatus, PGU), the twin-spotted ratsnake (Elaphe bimaculata), the leopard snake (Zamenis situla, ZSI), the grass snake (Natrix natrix, NNA), and the garter snake (Thamnophis sirtalis, TSI) showing shared microsatellite repeat motifs with the Siamese cobra W chromosome. Chromosomal locations of repeat motifs in the amniotes were obtained from the following sources: A. parapulchella from Matsubara et al. [53], E. quadrivirgata, R. tigrinus, P. flavoviridis and G. blomhoffii from Matsubara et al. [27, 28], B. duperreyi, P. vitticeps, C. longicollis, N. scutatus and R. tigrinus from Matsubara et al. [26, 27], X. javanicus from Rovatsos et al. [54], C. durissus, H. buccata, B. fuliginosus, B. olivaceus, L. triangulum, L. ruthveni, P. guttatus, E. bimaculata, Z. situla and N. natrix from Augstenová et al. [25], and T. sirtalis from Perry et al. [42]

Full size image

Amplification of (ACAG)n and (AATC)n was found on the W chromosome of tiger snakes (N. scutatus), while (AGAT)n was also observed in the Japanese striped snake (E. quadrivirgata), and the tiger keelback snake (R. tigrinus) [26] (Fig. 5). However, the sets of microsatellite repeat motifs used in this study generally differed from those in previous reports. Therefore, it cannot be ruled out that other microsatellite repeat motifs might be amplified in the Siamese cobra. Interestingly, both (AAGG)n and telomeric repeats identified on the central bearded dragon (P. vitticeps) W chromosome were also identified on NKAW, even though the species are only distantly related ([26, 41, 44]; Srikulnath et al. unpublished data). Similar results were also found in other amniotes (Fig. 5). This suggests that amplification of microsatellite repeat motifs occurred independently in each lineage and might represent convergent sex chromosomal differentiation among amniote sex chromosomes. It should be considered that although it is known that sex chromosomes share no homology among amniotes, evidence of synteny from several amniotes showed some overlap of partial synteny likely to have been part of an ancestral super-sex chromosome [4]. An example of such non-homologous W chromosomes of the chicken and snake sharing several repeat groups has been found in the common tiger snake (N. scutatus) [30]. Is it possible that amplified microsatellite repeat motifs were retained in the sex chromosomes of a common ancestor and subsequent reshuffling provided the appearance of sex chromosomes in each lineage?

Is squamate reptile chromosome 2 part of an ancestral super-sex chromosome in amniotes?

Except for gecko lizards and lacertid lizards, synteny on the snake Z chromosome and chromosome 2 are remarkably well conserved in most squamate reptiles [11,12,13,14, 16,17,18, 30]. Comparison of the chromosome maps of the Siamese cobra and other squamate reptiles revealed that subtelocentric NKA2 was homologous with submetacentric chromosome 2 of most squamate reptiles [11, 14, 16,17,18,19, 21]. This suggests that NKA2 resulted from intrachromosomal rearrangements (such as pericentric inversion or centromere repositioning) in an ancestral bi-armed macrochromosome. However, large hybridization signals like repeat sequences of one GGAZ BAC (CH261-129A16) and one TGUZ BAC (TGMCBA-200J22) observed on the Siamese cobra W chromosome were not present in other Siamese cobra chromosomes, even though these two BACs were expected to be located on NKA2. Similar results were also found in an additional eight BACs derived from GGA1, GGA4, TGU1B, TGU13, TGU17, TGU23, TGU27, which showed partial homology with various amniote sex chromosomes (Fig. 6) [4]. This suggests that these BACs contained repeat sequences shared with NKAW. Sequence analysis using RepeatMasker and MSATCOMMANDER showed 1.95–23.26% of repeat sequences for retrotransposon and 0.58–1.66 % for microsatellites repeats, also found in BACs mapped on NKAW. Transposable elements form the majority of shared elements in NKAW, supporting the co-location of heterochromatin and putative transposable elements as a result of the unchecked activity of retrotransposons in cells [45]. These repeat sequences probably drive these regions as derived sex chromosomes [3]. In agreement with the hypothesis of an ancestral amniote super-sex chromosome, observations of hybridization patterns from chicken and zebra finch BACs mapped to NKA2 and NKAW suggest that a large ancestral super-sex chromosome has strong correlation between squamate chromosome 2 (NKA2) and snake sex chromosomes. Alternatively, the shared partial synteny of amniote sex chromosomes could follow a stochastic pattern within a restricted set of species. It remains possible that the partial synteny of squamate chromosome 2 to several amniote sex chromosomes represents random homologies, given that only small sets of genes in a restricted set of species are involved. However, we also performed chromosome mapping using chicken and zebra finch BACs derived from other chromosomes and did not observe hybridization signals on NKA2, NKAZ, or NKAW. The same BACs were also used to perform chromosome mapping in other amniotes, which did not show hybridization signals such as repeat sequences ([36,37,38]; Singchat et al. unpublished data). Hence, mapping data on NKAW with several chicken and zebra finch BACs might not be random. However, limitations on the resolution of the FISH technique may indicate repeat sequences on NKAW, with repeats of other locations also escaping detection. At this stage, theories of independent rearrangements involving chromosomes 2, Z, and W in the snake lineage, and independent amplification of the same repeat sequences on snake and amniote sex chromosomes require further analysis of data from other squamate reptiles using FISH mapping and whole genome sequencing in order to investigate this hypothesis.

Fig. 6
figure 6

Chromosome maps of the Siamese cobra chromosome 2, and Z and W chromosomes showing homologies with the chicken, zebra finch, and several other amniotes. This map was constructed using 43 BACs derived from the chicken and zebra finch. Locations of the BACs on Siamese cobra chromosomes are shown to the right of the chromosomes. Chromosome numbers indicate the chicken (Gallus gallus, GGA), the green anole (Anolis carolinensis, ACA), the bearded dragon lizard (Pogona vitticeps, PVI), the Japanese four-striped rat snake (Elaphe quadrivirgata, EQU), humans (Homo sapiens, HSA), the tammar wallaby (Macropus eugnii, MEU), the duck-billed platypus (Ornithorhynchus anatinus, OAN), the giant musk turtle (Staurotypus salvinii, SSA), the mexican musk turtle (Staurotypus triporcatus, STR), the Hokou gecko (Gekko hokouenesis, GHO), and the sand lizard (Lacerta agilis, LAG) showing partial homologies with the Siamese cobra chromosomes. Partial syntenies are shown using the same color. Chromosomal locations of genes in the amniotes were obtained from the following sources: G. gallus from Matsuda et al. [33], A. carolinensis from Alföldi et al. [12], E. quadrivirgata from Matsubara et al. [21], G. hokouenesis from Kawai et al. [10], P. vitticeps from Ezaz et al. [13], H. sapiens and M. eugnii from Grützner et al. [55], O. anatinus from Veyrunes et al. [15], S. salvinii from Kawagoshi et al. [5] S. triporcatus from Montiel et al. [56], and L. agilis from Srikulnath et al. [14]

Full size image

Sex-determining genes are well established in several amniotes, and orthologs or paralogs of the functional genes appear repeatedly in distantly related amniotes. This suggests the possibility of cooperative favorable linked sex chromosomes with hitchhiking effects (genes physically linked to one of the chosen sex-determining genes undergoing selective sweep) or interaction between genes [4]. The sequence of chicken and zebra finch BACs mapped to NKAW contained at least five known protein-coding genes involving sex determination and sex development pathways (Table 4). Two genes (GPR149 and DHX36) located on TGU9 BAC (TGMCBA-217A3) were mapped to NKAW as a twin dot hybridization signal. This hybridization pattern is unlikely to be caused by repeat sequences. However, we have no direct evidence that these genes are involved in snake sex determination or sex development. It is likely that convergent evolution of sex chromosomes across distantly related taxa led to genomic elements that are particularly adept in a sex-determining role. Are these genes coincidental, or are there sequences in these regions that have a selectable function in sex determination? More information regarding genomic analysis and transcriptomic activity from squamate reptiles is required to test this hypothesis.

Conclusions

This study suggests the fundamental basis on which the Siamese cobra chromosome 2 and the W sex chromosome share partial synteny and some characteristics of sex chromosomes observed in other amniotes, for example the amplification of repeat sequences. This in turn supports the hypothesis of ancestral super-sex chromosomes with overlaps of partial syntenies among amniote sex chromosomes. This might result from random distribution of microsatellite repeat motifs in each lineage as an alternative hypothesis, representing convergent sex chromosomal differentiation among amniote sex chromosomes. Further research is required to compare additional maps and sequences to explore syntenies in several squamate reptiles, using recent advances in comparative gene mapping and whole-genome sequencing. This may provide a more comprehensive understanding of the complexity of evolutionary sex determination and sex chromosomes in amniotes.

Methods

Specimen Collection

Two male and two female Siamese cobras (N. kaouthia) were collected from Snake Farm, Queen Saovabha Memorial Institute (QSMI, Bangkok). The sex of each individual was identified morphologically and confirmed using a molecular sexing approach [22,23,24]. Blood samples were collected from the ventral tail vein using a 23-gauge needle attached to 2-ml disposable syringes. These contained 10 mM ethylenediaminetetraacetic acid (EDTA) for DNA extraction or 75 USP unit/ml heparin for cell culture. Whole genomic DNA was extracted following the standard salting-out protocol as described previously by Supikamolseni et al. [46] and used as templates for molecular sexing. DNA quality and concentration were determined using 1% agarose gel electrophoresis and spectrophotometric analysis.

Lymphocyte Cell culture and Chromosome Preparation

Lymphocytes from two male and two female Siamese cobras were isolated from peripheral blood, and then cultured for 5 days in RPMI 1640 medium supplemented with 15% fetal bovine serum (FBS), 3 μg/ml concanavalin A (type IV-S) (Sigma-Aldrich, St. Louis, MO, USA), 10 μg/ml lipopolysaccharide (Sigma-Aldrich), 1% phytohaemagglutinin (HA15) (Remel, Lenexa, KS, USA) and 1% Antibiotic-Antimycotic (Life Technologies-Gibco, Carlsbad, CA, USA). After 5 days lymphocytes were subject to colcemid treatment (100 ng/ml) for 60 min and fixed (3:1 methanol/acetic acid) after hypotonic treatment in 0.075 M KCl, before being harvested. The cell suspension was dropped onto clean glass slides and air-dried. The slides were kept at -80 °C until required for use.

Karyotyping and C-banding

Morphology and size of macrochromosomes were characterized according to Turpin and Lejeune [47] and Levan et al. [48]. To examine the chromosomal distribution of constitutive heterochromatin, C-banding was performed using the standard barium hydroxide/saline/Giemsa method [49] with slight modification as follows: chromosome slides were treated with 0.2 N HCl at room temperature for 60 min and then with 5% Ba (OH)2 at 50 °C for 30 s, followed by 2× SSC at 65 °C for 60 min.

Fluorescence in situ Hybridization (FISH) Mapping of Telomeric Repeat and Microsatellite Repeat Motifs

The chromosomal locations of telomeric (TTAGGG)n sequences and 19 microsatellite repeat motifs: (CA)15, (GC)15, (GA)15, (AT)15, (CAA)10, (CAG)10, (CAT)10, (CGG)10, (GAG)10, (AAT)10, (AAGG)8, (AATC)8, (AGAT)8, (ACGC)8, (AAAT)8, (AAAC)8, (AATG)8, (AAATC)6, and (AAAAT)6 were determined using FISH as described previously [16, 50]. We used commercially biotin-labeled 42-bp oligonucleotide complementary to (TTAGGG)n sequences and 19 commercially biotin-labeled oligonucleotide microsatellite repeat probes (Macrogen Co., Seoul, Korea), ethanol-precipitated with salmon sperm DNA, and Escherichia coli tRNA. After the hybridization of biotin-labeled probes to the Siamese cobra chromosomes, the probes were detected by incubating the chromosome slides with avidin labeled with fluorescein isothiocyanate (avidin-FITC; Invitrogen, CA, USA). Slides were subsequently stained with 1 μg/ml DAPI (4′, 6′-diamidino-2-phenylindole). Fluorescence hybridization signals were captured using a cooled Charge-Coupled Device (CCD) camera mounted on a ZEISS Axioplan2 microscope and processed using MetaSystems ISIS v.5.2.8 software (MetaSystems, Altlussheim, Germany).

Isolation, Amplification and Labeling of Chicken and Zebra finch Bacterial Artificial Chromosomes (BACs)

Chicken (Gallus gallus) and zebra finch (Taeniopygia guttata) BACs were applied for cross-species FISH mapping according to a range of the proportion of conserved elements shared across multiple species. Due to the high degree of apparent genome conservation between avian and reptilian species ([36,37,38]; Singchat et al. unpublished data), these sets of BACs were applied to the Siamese cobras. All BACs were anchored to the chicken and zebra finch genome assembly by linkage and sequencing. Twenty-six chicken and seventeen zebra finch BACs, comprising chicken chromosome 1 (GGA1), GGA2p, GGA4, GGA5, GGA9, GGA15, GGA17, GGA18, GGA23, GGA27 and GGAZ, and zebra finch chromosome 1 (TGU1B), TGU4A, TGU5, TGU9, TGU12, TGU13, TGU15, TGU17, TGU23, TGU27 and TGUZ, were isolated using the Qiagen Miniprep Kit (Qiagen, Manchester, UK) prior to amplification and direct labeling by nick translation (Roche, Welwyn Garden City, UK). Probes were labeled with Texas Red 12-dUTP (Invitrogen Corporation and Applied Biosystems Inc., Carlsbad, California, US) and FITC-12-UTP (Roche) prior to purification using the Qiagen Nucleotide Removal Kit (Qiagen).

Cross-Species Bacterial Artificial Chromosome (BAC) Fluorescence in situ Hybridization (FISH) Mapping

Chromosome slides were dehydrated through an ethanol series (2 min each in 2×SSC, 70%, 85%, and 100% ethanol at room temperature). Probes were diluted in a Hybridization Solution I (Cytocell Ltd., Cambridge, UK) with Chicken Hybloc (Insight Biotechnology limited, Wembley, London, UK) and applied to the Siamese cobra chromosomes on a 37 °C hotplate before being sealed with rubber cement. Probe and target DNA were simultaneously denatured on a 75 °C hotplate for 2 min prior to hybridization in a humidified chamber at 37 °C for 72 h. Slides were washed post-hybridization for 30 s in 2 × SSC/0.05% Tween 20 at room temperature and then counterstained using VECTASHIELD antifade mounting medium with DAPI (Vector Laboratories, Inc., Burlingame, CA, USA). Images were captured using an Olympus BX61 epifluorescence microscope with a cooled CCD camera and SmartCapture system (Digital Scientific UK Ltd., Cambridge, UK). Confirmation of BAC signal order was achieved by dual color experiments where Texas Red 12-dUTP- and FITC-12-UTP- labeled probes were hybridized simultaneously.

Candidate Bacterial Artificial Chromosome (BAC) Sequence Analysis

Gene position and gene structure in each candidate BAC were derived from ENSEMBL (ensembl.org/Gallus_gallus/Info/Annotation). The CpG content and repeat information were retrieved from the UCSC Genome Browser http://genome.ucsc.edu/ for the corresponding gene build. Simultaneously, each BAC sequence was subjected to homology search using Blastx and Blastn (RefSeq G. gallus and T. guttata). Repeat data were identified using Arian Smit's Repeat Masker program [51], which screens DNA sequences for interspersed repeats and low complexity DNA sequences. Query species were assumed to be Vertebrata, Metazoa, Repeat Masker Combined Database; Dfam_Consensus-20170127 and RepBase-20170127, run with rmblastn version 2.2.27+ [51]. Microsatellite repeat motifs were scanned from candidate BAC sequences using MSATCOMMANDER [52]. Parameters of the minimum number of repetitions were set at 10, 6, 4, 4, 4, and 4 for mono, di, tri, tetra, penta, and hexanucleotide motifs, respectively.

Abbreviations

ACA:

Anolis carolinensis

APA:

Aprasia parapulchella

BACs:

bacterial artificial chromosomes

BDU:

Bassiana duperreyi

BFU:

Boaedon fuliginosus

Bkm:

banded krait minor satellite

BOL:

Boaedon olivaceus

CCD:

Charge-Coupled Device

CDU:

Crotalus durissus

CLO:

Chelodina longicollis

DAPI:

4′, 6-diamidino-2-phenylindole

EBI:

Elaphe bimaculata

EQU:

Elaphe quadrivirgata

FBS:

fetal bovine serum

FISH:

Fluorescent in situ hybridization

FITC:

Fluorescein isothiocyanate

GBL:

Gloydius blomhoffii

GGA:

Gallus gallus

GHO:

Gekko hokouenesis

HBU:

Homalopsis buccata

HSA:

Homo sapiens

LAG:

Lacerta agilis

LRU:

Lampropeltis ruthveni

LTR:

Lampropeltis triangulum

MEU:

Macropus eugnii

NKA:

Naja kaouthia

NNA:

Natrix natrix

NSC:

Notechis scutatus

OAN:

Ornithorhynchus anatinus

PFL:

Protobothrops flavoviridis

PGU:

Pantherophis guttatus

PVI:

Pogona vitticeps

QSMI:

Queen Saovabha Memorial Institute

RTI:

Rhabdophis tigrinus

SSA:

Staurotypus salvinii

STR:

Staurotypus triporcatus

TGU:

Taeniopygia guttata

TSI:

Thamnophis sirtalis

XJA:

Xenodermus javanicus

ZSI:

Zamenis situla

References

  1. Charlesworth B. The evolution of sex chromosomes. Science. 1991;251:1030–3.

    Article  CAS  Google Scholar 

  2. Bachtrog D, Mank JE, Peichel CL, et al. Sex Determination: why so many ways of doing it? PLoS Biol. 2014;12:e1001899.

  3. Charlesworth D, Charlesworth B, Marais G. Steps in the evolution of heteromorphic sex chromosomes. Heredity (Edinb). 2005;95:118–28.

    Article  CAS  Google Scholar 

  4. Ezaz T, Srikulnath K, Graves JA. Origin of amniote sex chromosomes: an ancestral super-sex chromosome, or common requirements? J Hered. 2017;108:94–105.

    Article  Google Scholar 

  5. Kawagoshi T, Uno Y, Nishida C, Matsuda Y. The Staurotypus Turtles and Aves share the same origin of sex chromosomes but evolved different types of heterogametic sex determination. PLoS One. 2014;9:e105315.

    Article  Google Scholar 

  6. Kawagoshi T, Uno Y, Matsubara K, Matsuda Y, Nishida C. The ZW micro-sex chromosomes of the Chinese soft-shelled turtle (Pelodiscus sinensis, Trionychidae, Testudines) have the same origin as chicken chromosome 15. Cytogenet Genome Res. 2009;125:125–31.

    Article  CAS  Google Scholar 

  7. Kawagoshi T, Nishida C, Matsuda Y. The origin and differentiation process of X and Y chromosomes of the black marsh turtle (Siebenrockiella crassicollis, Geoemydidae, Testudines). Chromosome Res. 2012;20:95–110.

    Article  CAS  Google Scholar 

  8. Badenhorst D, Stanyon R, Engstrom T, Valenzuela N. A ZZ/ZW microchromosome system in the spiny softshell turtle, Apalone spinifera, reveals an intriguing sex chromosome conservation in Trionychidae. Chromosome Res. 2013;21:137–47.

    Article  CAS  Google Scholar 

  9. Montiel EE, Badenhorst D, Tamplin J, Burke RL, Valenzuela N. Discovery of the youngest sex chromosomes reveals first case of convergent co-option of ancestral autosomes in turtles. Chromosoma. 2016;126:105–13.

    Article  Google Scholar 

  10. Kawai A, Ishijima J, Nishida C, et al. The ZW sex chromosomes of Gekko hokouenesis (Gekkonidae, Squamata) represent highly conserved homology with those of avian species. Chromosoma. 2009;118:43–51.

    Article  Google Scholar 

  11. Srikulnath K, Uno Y, Nishida C, Ota H, Matsuda Y. Karyotype reorganization in the Hokou gecko (Gekko hokouenesis, Gekkonidae): the process of microchromosome disappearance in Gekkota. PLoS One. 2015;10:e0134829.

  12. Alföldi J, Palma FD, Grabherr M, et al. The genome of the green anole lizard and a comparative analysis with birds and mammals. Nature. 2011;477:587.

    Article  Google Scholar 

  13. Ezaz T, Azad B, O’Meally D, et al. Sequence and gene content of a large fragment of a lizard sex chromosome and evaluation of candidate sex differentiating gene R-spondin 1. BMC Genomics. 2013;14:899.

    Article  Google Scholar 

  14. Srikulnath K, Matsubara K, Uno Y, Nishida C, Olsson M, Matsuda Y. Identification of the linkage group of the Z sex chromosomes of the sand lizard (Lacerta agilis, Lacertidae) and elucidation of karyotype evolution in lacertid lizards. Chromosoma. 2014;123:563–75.

    Article  Google Scholar 

  15. Veyrunes F, Waters PD, Miethke P, et al. Bird-like sex chromosomes of platypus imply recent origin of mammal sex chromosomes. Genome Res. 2008;18:965–73.

    Article  CAS  Google Scholar 

  16. Srikulnath K, Matsubara K, Uno Y, et al. Karyological characterization of the butterfly lizard (Leiolepis reevesii rubritaeniata, Agamidae, Squamata) by molecular cytogenetic approach. Cytogenet Genome Res. 2009;125:213–23.

    Article  CAS  Google Scholar 

  17. Srikulnath K, Nishida C, Matsubara K, et al. 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–86.

    Article  CAS  Google Scholar 

  18. Srikulnath K, Uno Y, Nishida C, Matsuda Y. Karyotype evolution in monitor lizards: cross-species chromosome mapping of cDNA reveals highly conserved synteny and gene order in the Toxicofera clade. Chromosome Res. 2013;21:805–19.

    Article  CAS  Google Scholar 

  19. Pokorná M, Giovannotti M, Kratochvíl L, et al. Strong conservation of the bird Z chromosome in reptilian genomes is revealed by comparative painting despite 275 million years divergence. Chromosoma. 2011;120:455–68.

    Article  Google Scholar 

  20. Gamble T, Castoe TA, Nielsen SV, et al. The discovery of XY sex chromosomes in a Boa and Python. Curr Biol. 2017;27:2148–2153.e4.

    Article  CAS  Google Scholar 

  21. Matsubara K, Tarui H, Toriba M, et al. Evidence for different origin of sex chromosomes in snakes, birds, and mammals and step-wise differentiation of snake sex chromosomes. Proc Natl Acad Sci U S A. 2006;103:18190–5.

    Article  CAS  Google Scholar 

  22. Laopichienpong N, Muangmai N, Chanhome L, et al. Evolutionary dynamics of the gametologous CTNNB1 gene on the Z and W chromosomes of snakes. J Hered. 2017;108:142–51.

    PubMed  Google Scholar 

  23. Laopichienpong N, Tawichasri P, Chanhome L, et al. A novel method of caenophidian snake sex identification using molecular markers based on two gametologous genes. Ecol Evol. 2017;7:4661–9.

    Article  Google Scholar 

  24. Tawichasri P, Laopichienpong N, Chanhome L, et al. Using blood and non-invasive shed skin samples to identify sex of caenophidian snakes based on multiplex PCR assay. Zool Anz. 2017;271:6–14.

    Article  Google Scholar 

  25. Augstenová B, Mazzoleni S, Kratochvíl L, Rovatsos M. Evolutionary dynamics of the W chromosome in caenophidian snakes. Genes. 2018;9:5.

    Article  Google Scholar 

  26. Matsubara K, O'Meally D, Azad B, et al. Amplification of microsatellite repeat motifs is associated with the evolutionary differentiation and heterochromatinization of sex chromosomes in Sauropsida. Chromosoma. 2015a;125:111–23.

    Article  Google Scholar 

  27. Matsubara K, Uno Y, Srikulnath K, Matsuda Y, Miller E, Olsson M. No interstitial telomeres on autosomes but remarkable amplification of telomeric repeats on the W sex chromosome in the sand lizard (Lacerta agilis). J Hered. 2015b;106:753–7.

    Article  CAS  Google Scholar 

  28. Matsubara K, Nishida C, Matsuda Y, Kumazawa Y. Sex chromosome evolution in snakes inferred from divergence patterns of two gametologous genes and chromosome distribution of sex chromosome-linked repetitive sequences. Zool Lett. 2016;2:19.

    Article  Google Scholar 

  29. Oguiura N, Ferrarezzi H, Batistic RF. Cytogenetics and molecular data in snakes: a phylogenetic approach. Cytogenet Genome Res. 2009;127:128–42.

    Article  CAS  Google Scholar 

  30. O’Meally D, Patel HR, Stiglec R, et al. Non-homologous sex chromosomes of birds and snakes share repetitive sequences. Chromosome Res. 2010;18:787–800.

    Article  Google Scholar 

  31. Epplen JT, McCarrey JR, Sutou S, Ohno S. Base sequence of a cloned snake W-chromosome DNA fragment and identification of a male-specific putative mRNA in the mouse. Proc Natl Acad Sci U S A. 1982;79:3798–802.

    Article  CAS  Google Scholar 

  32. Vicoso B, Emerson JJ, Zektser Y, Mahajan S, Bachtrog D. Comparative sex chromosome genomics in snakes: differentiation, evolutionary strata, and lack of global dosage compensation. PLoS Biol. 2013;11:e1001643.

    Article  CAS  Google Scholar 

  33. Matsuda Y, Nishida-Umehara C, Tarui H, et al. Highly conserved linkage homology between birds and turtles: bird and turtle chromosomes are precise counterparts of each other. Chromosome Res. 2005;13:601–15.

    Article  CAS  Google Scholar 

  34. Matsubara K, Kuraku S, Tarui H, et al. Intra-genomic GC heterogeneity in sauropsids: evolutionary insights from cDNA mapping and GC (3) profiling in snake. BMC Genomics. 2012;13:604.

    Article  CAS  Google Scholar 

  35. Olmo E, Signorino G. Chromorep: a reptile chromosomes. database. 2005; http://chromorep.univpm.it. Accessed 26 July 2018.

  36. Damas J, O'Connor R, Farré M, et al. Upgrading short-read animal genome assemblies to chromosome level using comparative genomics and a universal probe set. Genome Res. 2017;27:875–84.

    Article  CAS  Google Scholar 

  37. O’Connor RE, Romanov MN, Kiazim LG, et al. Reconstruction of the diapsid ancestral genome permits chromosome evolution tracing in avian and non-avian dinosaurs. Nat Commun. 2018;9:1883.

    Article  Google Scholar 

  38. O'Connor RE, Farré M, Joseph S, et al. Chromosome-level assembly reveals extensive rearrangement in saker falcon and budgerigar, but not ostrich, genomes. Genome Biol. 2018;19:171.

    Article  Google Scholar 

  39. Graves JAM. Sex chromosome specialization and degeneration in mammals. Cell. 2006;124:901–14.

    Article  Google Scholar 

  40. Schartl M, Schmid M, Nanda I. Dynamics of vertebrate sex chromosome evolution: from equal size to giants and dwarfs. Chromosoma. 2016;125:553–71.

    Article  Google Scholar 

  41. Young MJ, O’Meally D, Sarre SD, Georges A, Ezaz T. Molecular cytogenetic map of the central bearded dragon, Pogona vitticeps (Squamata: Agamidae). Chromosome Res. 2013;21:361–74.

    Article  CAS  Google Scholar 

  42. Perry BW, Card DC, McGlothlin JW, et al. Molecular adaptations for sensing and securing prey and insight into amniote genome diversity from the garter snake genome. Genome Biol Evol. 2018;10:2110–29.

    Article  Google Scholar 

  43. Rovatsos M, Altmanová M, Pokorná JM, Augstenová B, Kratochvíl L. Cytogenetics of the Javan file snake (Acrochordus javanicus) and the evolution of snake sex chromosomes. J Zool Syst Evol Res. 2018;56:117–25.

    Article  Google Scholar 

  44. Pyron RA, Burbrink FT, Wiens JJ. A phylogeny and revised classification of Squamata, including 4161 species of lizards and snakes. BMC Evol Biol. 2013;13:93.

    Article  Google Scholar 

  45. Steinemann S, Steinemann M. Retroelements: tools for sex chromosome evolution. Cytogenetic Genom Res. 2005;110:134–43.

    Article  CAS  Google Scholar 

  46. Supikamolseni A, Ngaoburanawit N, Sumontha M, et al. Molecular barcoding of venomous snakes and species-specific multiplex PCR assay to identify snake groups for which antivenom is available in Thailand. Genet Mol Res. 2015;14:13981–97.

    Article  CAS  Google Scholar 

  47. Turpin R, Lejeune J. Les Chromosomes Humains. Paris: Gauthier-Villars; 1965.

    Google Scholar 

  48. Levan A, Fredga K, Sandberg A. Nomenclature for centromeric position on chromosomes. Hereditas. 1964;52:201–20.

    Article  Google Scholar 

  49. Sumner AT. A simple technique for demonstrating centromeric heterochromatin. Exp Cell Res. 1972;75:304–6.

    Article  CAS  Google Scholar 

  50. Matsuda Y, Chapman VM. Application of fluorescence in situ hybridization in genome analysis of the mouse. Electrophoresis. 1995;16:261–72.

    Article  CAS  Google Scholar 

  51. Smit AFA, Hubley R, Green P. RepeatMasker Open-4.0. 2015. http://www.repeatmasker.org. Accessed 26 July 2018.

  52. Faircloth BC. MSATCOMMANDER: detection of microsatellite repeat arrays and automated, locus-specific primer design. Mol Ecol Resour. 2008;8:92–4.

    Article  CAS  Google Scholar 

  53. Matsubara K, Knopp T, Sarre SD, Georges A, Ezaz T. Karyotypic analysis and FISH mapping of microsatellite motifs reveal highly differentiated XX/XY sex chromosomes in the pink-tailed worm-lizard (Aprasia parapulchella, Pygopodidae, Squamata). Mol Cytogenet. 2013;6:60.

    Article  Google Scholar 

  54. Rovatsos M, Johnson Pokorná M, Kratochvíl L. Differentiation of sex chromosomes and karyotype characterisation in the dragonsnake Xenodermus javanicus (Squamata: Xenodermatidae). Cytogenet. Genome Res. 2015;147:48–54.

    Article  Google Scholar 

  55. Grützner F, Rens W, Tsend-Ayush E, et al. In the platypus a meiotic chain of ten sex chromosomes shares genes with the bird Z and mammal X chromosomes. Nature. 2004;432:913–7.

    Article  Google Scholar 

  56. Montiel EE, Badenhorst D, Lee LS, Literman R, Trifonov V, Valenzuela N. Cytogenetic insights into the evolution of chromosomes and sex determination reveal striking homology of turtle sex chromosomes to amphibian autosomes. Cytogenet Genome Res. 2016;148:292–304.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

We would like to thank Tanapong Tawan, Taksa Vasaruchapong and Panithi Laoungbua (Snake Farm, Queen Saovabha Memorial Institute, Thai Red Cross Society, Thailand), and Lucas Kiazim (University of Kent, UK), and the Conservation Research and Education Division, Zoological Park Organization, Dusit, Bangkok, Thailand for help with sample collection.

Funding

This study was financially supported by grants from the Thailand Research Fund (TRF; no. RSA6180075, PHD60I0014, PHD60I0082, MSD60I0035) awarded to KS, WS, AS and PT, the Thailand Research Fund - Newton Fund Placement, Travel Grant for PhD Supervisors (GA/PhD/Sup/Year4/003) awarded to KS, Kasetsart University Research and Development Institute (KURDI; grant no. 80.60) awarded to KS, the Fellowship of Capacity Building for Kasetsart University on Internationalization at Kasetsart University (No.0513.10109/8384) awarded to WA and KS, the Center for Advanced Studies in Tropical Natural Resources, National Research University-Kasetsart University (CASTNAR, NRU-KU, Thailand) awarded to KS, the National Research Council of Thailand (NRCT; No. 2560096003012) awarded to SB and KS, the Center of Excellence on Agricultural Biotechnology, Science and Technology Postgraduate Education and Research Development Office, Office of Higher Education Commission, Ministry of Education (AG-BIO/PERDO-CHE) awarded to KS. ROC was supported by Biotechnology and Biological Sciences Research Council (BBSRC) grant (BB/K008161/1) awarded to DK.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Author information

Authors and Affiliations

  1. Laboratory of Animal Cytogenetics and Comparative Genomics (ACCG), Department of Genetics, Faculty of Science, Kasetsart University, Bangkok, 10900, Thailand

    Worapong Singchat, Panupong Tawichasri, Aorarat Suntronpong, Siwapech Sillapaprayoon & Kornsorn Srikulnath

  2. Animal Breeding and Genetics Consortium of Kasetsart University (ABG-KU), Bangkok, 10900, Thailand

    Worapong Singchat, Panupong Tawichasri, Aorarat Suntronpong, Siwapech Sillapaprayoon, Narongrit Muangmai & Kornsorn Srikulnath

  3. School of Biosciences, University of Kent, Canterbury, CT2 7NY, United Kingdom

    Rebecca E. O’Connor & Darren Griffin

  4. Department of Research and Development, Queen Saovabha Memorial Institute, The Thai Red Cross Society, Bangkok, 10330, Thailand

    Sunutcha Suntrarachun

  5. Department of Fishery Biology, Faculty of Fisheries, Kasetsart University, Bangkok, 10900, Thailand

    Narongrit Muangmai

  6. Bureau of Conservation and Research, Zoological Park Organization under the Royal Patronage of His Majesty the King, Bangkok, Thailand

    Sudarath Baicharoen

  7. Department of Biology, Faculty of Science, Naresuan University, Phitsanulok, 65000, Thailand

    Surin Peyachoknagul

  8. Snake Farm, Queen Saovabha Memorial Institute, The Thai Red Cross Society, Bangkok, 10330, Thailand

    Lawan Chanhome

  9. Center for Advanced Studies in Tropical Natural Resources, National Research University-Kasetsart University, Kasetsart University, Bangkok, 10900, Thailand

    Kornsorn Srikulnath

  10. Center of Excellence on Agricultural Biotechnology (AG-BIO/PERDO-CHE), Bangkok, 10900, Thailand

    Kornsorn Srikulnath

Authors
  1. Worapong Singchat
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Rebecca E. O’Connor
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Panupong Tawichasri
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Aorarat Suntronpong
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Siwapech Sillapaprayoon
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Sunutcha Suntrarachun
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Narongrit Muangmai
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Sudarath Baicharoen
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Surin Peyachoknagul
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Lawan Chanhome
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Darren Griffin
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Kornsorn Srikulnath
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

W.S. and K.S. drafted the manuscript; W.S., R.O., D.G., and K.S. conceived the ideas and designed methodology; W.S., R.O., P.T., A.S., and K.S. carried out the lab work; W.S. and K.S. participated in data analysis; W.S., R.O., P.T., A.S., S.S., S.S., N.M., S.B., S.P., L.C., D.G., and K.S. reviewed the data and the manuscript. All authors gave final approval for publication.

Corresponding author

Correspondence to Kornsorn Srikulnath.

Ethics declarations

Ethics approval and consent to participate

Animal care and all experimental procedures were approved by the Animal Experiment Committee, Kasetsart University, Thailand (approval no. ACKU59-SCI-034 and ACKU61-SCI-024) and conducted in accordance with the Regulations on Animal Experiments at Kasetsart University.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional information

Worapong Singchat and Rebecca E. O’Connor are joint first authors.

Darren Griffin and Kornsorn Srikulnath are joint last authors.

Additional files

Additional file 1:

Figure S1. Chromosomal locations of microsatellite repeat motifs in the Siamese cobra (Naja kaouthia). Hybridization patterns of FITC-labeled (CA)15 (a), (AT)15 (b), (GC)15 (c), (CAT)10 (d), and (AAT)10 (e) on DAPI-stained chromosomes. Scale bar represents 10 μm. (JPG 175 kb)

Additional file 2:

Figure S2. Cytogenetic map of the Siamese cobra (Naja kaouthia), which shows chromosome homologies with chicken and zebra finch. This map was constructed with 25 chicken and zebra finch BACs mapped on the Siamese cobra macrochromosomes. Locations of BACs are shown to the right of the Siamese cobra chromosomes. The chromosome numbers show the chromosomes of the chicken (Gallus gallus, GGA) and zebra finch (Taeniopygia guttata, TGU), which show homologies with the Siamese cobra chromosomes. (JPG 118 kb)

Additional file 3:

Table S1. Comparison of major classes of repeat sequences in chicken and zebra finch BACs mapped on the Siamese cobra Z chromosome. (DOCX 14 kb)

Additional file 4:

Table S2. Comparison of frequencies of microsatellite repeat motifs in chicken and zebra finch BACs mapped on the Siamese cobra Z chromosome. (DOCX 13 kb)

Additional file 5:

Table S3. Comparison of major classes of repeat sequences in chicken and zebra finch BACs mapped on the Siamese cobra chromosome 2. (DOCX 15 kb)

Additional file 6:

Table S4. Comparison of frequencies of microsatellite repeat motifs in chicken and zebra finch BACs mapped on the Siamese cobra chromosome 2. (DOCX 14 kb)

Additional file 7:

Table S5. Comparison of major classes of repeat sequences in 22 chicken and zebra finch BACs mapped on the Siamese cobra chromosome 1 and microchromosomes. (DOCX 23 kb)

Additional file 8:

Table S6. Comparison of frequencies of microsatellite repeat motifs in 22 chicken and zebra finch BACs mapped on the Siamese cobra chromosome 1 and microchromosomes. (DOCX 19 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Singchat, W., O’Connor, R.E., Tawichasri, P. et al. Chromosome map of the Siamese cobra: did partial synteny of sex chromosomes in the amniote represent “a hypothetical ancestral super-sex chromosome” or random distribution?. BMC Genomics 19, 939 (2018). https://doi.org/10.1186/s12864-018-5293-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-018-5293-6

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us