Intra-genomic GC heterogeneity in sauropsids: evolutionary insights from cDNA mapping and GC3 profiling in snake
© Matsubara et al.; licensee BioMed Central Ltd. 2012
Received: 28 January 2012
Accepted: 24 October 2012
Published: 9 November 2012
Extant sauropsids (reptiles and birds) are divided into two major lineages, the lineage of Testudines (turtles) and Archosauria (crocodilians and birds) and the lineage of Lepidosauria (tuatara, lizards, worm lizards and snakes). Karyotypes of these sauropsidan groups generally consist of macrochromosomes and microchromosomes. In chicken, microchromosomes exhibit a higher GC-content than macrochromosomes. To examine the pattern of intra-genomic GC heterogeneity in lepidosaurian genomes, we constructed a cytogenetic map of the Japanese four-striped rat snake (Elaphe quadrivirgata) with 183 cDNA clones by fluorescence in situ hybridization, and examined the correlation between the GC-content of exonic third codon positions (GC3) of the genes and the size of chromosomes on which the genes were localized.
Although GC3 distribution of snake genes was relatively homogeneous compared with those of the other amniotes, microchromosomal genes showed significantly higher GC3 than macrochromosomal genes as in chicken. Our snake cytogenetic map also identified several conserved segments between the snake macrochromosomes and the chicken microchromosomes. Cross-species comparisons revealed that GC3 of most snake orthologs in such macrochromosomal segments were GC-poor (GC3 < 50%) whereas those of chicken orthologs in microchromosomes were relatively GC-rich (GC3 ≥ 50%).
Our results suggest that the chromosome size-dependent GC heterogeneity had already occurred before the lepidosaur-archosaur split, 275 million years ago. This character was probably present in the common ancestor of lepidosaurs and but lost in the lineage leading to Anolis during the diversification of lepidosaurs. We also identified several genes whose GC-content might have been influenced by the size of the chromosomes on which they were harbored over the course of sauropsid evolution.
KeywordsGC-content Lepidosauria Snake Macrochromosome Microchromosome
Molecular phylogenetic analyses have suggested that extant sauropsids (reptiles and birds) are divided into two major groups, the lineage of Testudines (turtles) and Archosauria (crocodilians and birds) and the lineage of Lepidosauria (tuatara, lizards, worm lizards and snakes) although phylogenetic position of Testudines is still debatable [1–7]. The divergence time between the two lineages has been estimated at around 275 million years [3, 7–9]. Most sauropsidan species have karyotypes consisting of macrochromosomes and microchromosomes, as for birds [10–17], except for crocodilian species, whose karyotypes contain no microchromosomes [18, 19].
Whole genome sequencing of chicken revealed that the overall GC-content of chromosomes increases as chromosomal size decreases, that is, microchromosomes exhibit a higher GC-content than macrochromosomes [20, 21]. In a compositional map of GC-content constructed by 100-kb window analysis for the chicken whole genome sequence, most microchromosomes were occupied by GC-rich DNA segments, whereas GC-poor segments were more common in macrochromosomes . The differences of other features such as gene density, distribution in interphase nuclei and rate of nucleotide divergence were also identified between the two chromosomal groups of birds [23–27].
Reptiles are crucial taxon for tracking genome evolution in amniotes [21, 28, 29]. Intra-genomic GC heterogeneity has been found in reptiles by calculating GC-content in exonic third positions (GC3) [21, 30–33]. Although the use of GC3 as a proxy for genomic GC-content has been controversial , it is known that GC3 generally reflects the local GC-content of the introns and flanking regions of a gene [21, 35–37]. Chojnowski et al.  analyzed the GC3 of more than 6,000 ESTs in the American alligator (Alligator mississippiensis) and suggested that the alligator genome has a certain level of GC heterogeneity. They also examined the isochore structure of the red-eared slider turtle ( Trachemys scripta ) and suggested that the isochore structure of the turtle is intermediate between that of a frog and the GC-rich isochore structures of archosaurs and mammals . However, the chromosomal distribution of the GC heterogeneity has not been fully investigated in reptiles.
We previously constructed a cytogenetic map with 90 cDNA clones for the Chinese soft-shelled turtle (Pelodiscus sinensis), which revealed that the chromosomes have been highly conserved between the turtle and chicken, with the six largest chromosomes being almost equivalent to each other . GC3 of the mapped genes showed a heterogeneous distribution, and orthologs exhibited similar GC3 levels between the turtle, chicken and human, suggesting that the intra-genome GC heterogeneity had already occurred in the last common ancestor of extant amniotes . Furthermore, our results suggested that the turtle microchromosomes tend to contain more GC-rich genes than GC-poor genes, as in chicken .
The green anole lizard (Anolis carolinensis) is the first reptilian species for which whole genomic sequence has been released . Anolis has a homogeneous genome composition compared with other amniotes [37, 39] and, unlike chicken, the GC-content is similar between macro- and microchromosomes . However, it remains unknown whether these genomic characteristics are common to other lepidosaurs or not. Snake karyotypes have been highly conserved within the group, and the usual diploid number is 2n = 36, consisting of eight pairs of macrochromosomes and 10 pairs of microchromosomes [10, 40, 41]. The chromosome number is largely different from the chicken karyotype (2n=78) because of the remarkable difference in the number of microchromosomes. The snake therefore provides an ideal system for exploring changes in GC-content between macro- and microchromosomes over the course of sauropsid evolution.
Previously we constructed a cytogenetic map with 109 cDNA clones for the Japanese four-striped rat snake, Elaphe quadrivirgata (Serpentes, Colubridae) [38, 42]. In this study, we have extended cDNA-based chromosome mapping of the snake genes and consequently constructed a cytogenetic map with a total of 183 genes. We compared GC3 of the mapped snake genes with GC3 of their orthologs of chicken, green anole lizard, Chinese soft-shelled turtle, human and Xenopus tropicalis. This highlighted the chromosome size-dependent GC heterogeneity in the snake genome and the shift of GC-content possibly caused by chromosome rearrangements during sauropsid evolution.
Selection of EST clones
A cytogenetic map with 109 cDNA clones was constructed in our previous study [33, 36]. In the present study, we searched the snake EST clones isolated from the cDNA library constructed from brain , selected clones with significant similarity (E-value < 2e-35) to human and/or chicken genes in BLASTX , and used them for chromosome mapping (Additional file 1 and Additional file 2).
We rigorously confirmed orthologies of the snake sequences to their homologs of other vertebrates by constructing molecular phylogenetic trees with the neighbor-joining method  using XCED in which the alignment algorithm MAFFT is implemented  and with the maximum-likelihood method using PhyML . Sequence IDs of orthologs in six species (Anolis carolinensis, Gallus gallus, Pelodiscus sinensis, Homo sapiens, Mus musculus and Xenopus tropicalis) are included in Additional file 2. When multiple more than one potential ortholog was detected for a snake gene, we used the sequence with the greatest similarity to snake for cross-species comparison of GC-content.
Chromosome preparation and fluorescence in situ hybridization
Cell culture, preparation of R-banded chromosomes and fluorescence in situ hybridization (FISH) were performed as described previously [38, 47]. Fibroblast cells derived from lung tissues of the Japanese four-striped rat snake were cultured and used for chromosome preparations. DNA probes were labeled by nick translation with biotin-16-dUTP (Roche Diagnostics). The hybridized cDNA probes were reacted with goat anti-biotin antibodies (Vector Laboratories), and then stained with Alexa488-labeled donkey anti-goat IgG (Molecular Probes).
Calculation of GC-content
GC3 and GC-content at fourfold degenerate sites (GC4) were calculated using an original Perl script with the Bioperl module . The calculation was automatically processed on the basis of the open reading frame identified by a pairwise alignment between a translated nucleotide sequence and amino acid sequences of orthologs using BLASTX . When multiple alternative splicing variants were found for one gene, we used only the one that had the longest stretch of sequence homology with its orthologs of other species. We arbitrarily classified genes into GC-rich (GC3 ≥ 50%) and GC-poor (GC3 < 50%) genes based on the GC3.
Identification of orthologous sequences in the Burmese python
We conducted nucleotide BLAST for whole genome shotgun sequence of Burmese python, Python molurus bivittatus using rat snake ESTs as queries. We selected python sequences (consisting of exons, introns and flanking regions) that exhibited high similarities for rat snake ESTs (Additional file 3). We deduced the non-coding regions and the protein coding regions within each python genomic region using Wise2 program . We then calculated GC-content of non-coding regions and GC3 of coding regions.
Gene location in sequenced genomes
Chromosome locations of chicken, human and mouse orthologs were retrieved from Ensembl (version 56)  and NCBI Entrez Gene . In this study, we categorized chicken chromosomes 1–8, Z and W as macrochromosomes, and the remaining chromosomes as microchromosomes. Recently the second version of the genome assembly and annotation of green anole lizard (Anolis carolinensis) was released in Ensembl (version 61) . However, locations of about half of all A. carolinensis genes are at present annotated only at the scaffold level. We therefore did not compare the chromosome locations of orthologs between the snake and green anole lizard.
Comparison of gene characteristics between chicken macrochromosomes and microchromosomes
We classified the chicken genes into two groups, macrochromosomal genes and microchromosomal genes, and examined the over- and under-representation of gene functions between the two groups by FatiGO . FatiGO detects over-represented functional categories of Gene Ontology (GO), KEGG pathway, InterPro motif and Swissprot in either group between two gene lists using Fisher’s exact test. The Ensembl IDs of all chicken genes whose chromosome locations are known were downloaded and used for comparison by FatiGO.
Cytogenetic map of the Japanese four-striped rat snake
Of the 183 genes, 144 genes were mapped to macrochromosomes (chromosomes 1–7, Z and W chromosomes), and the others were mapped to microchromosomes (Figure 1). Twenty-nine segments in the snake chromosomes 1–7 and the Z chromosome were conserved between the snake and chicken. Most of them had a one to one correspondence to a particular region of chicken chromosomes. However, chromosomal homology for each of chicken chromosomes 1, 2, 4, 6 and 7 was found on more than two snake macrochromosomes, indicating that some inter-chromosomal rearrangements occurred between the snake and chicken macrochromosomes (Figure 1).
Chicken orthologs of 36 snake microchromosomal genes were located on chicken microchromosomes (chromosomes 10, 11, 14, 15, 17, 21, 22, 24 and 25). PGK1, ATRX and STAG2 genes on the snake microchromosomes are localized to the short arm of chicken chromosome 4, which was derived from a microchromosome fused with acrocentric chromosome 4 of the avian ancestor [16, 17, 54]. Since all the snake microchromosomes corresponded to avian microchromosomes, they have likely been retained from the ancestral karyotype of extant sauropsids without dynamic chromosome rearrangements.
Linkage homologies with chicken microchromosomes 9, 10, 12, 13, 18–20 and 26–28 were found on the snake macrochromosomes (Figure 1). For example, chicken chromosomes 19 and 28 were homologous to the distal segments of the short and long arms of snake chromosome 1, respectively. These results confirmed our previous assumption  that the large differences of chromosome numbers between the snakes (2n = 36) and chicken (2n = 78) resulted from frequent chromosome rearrangements containing fusions between macro- and microchromosomes and also between microchromosomes in the lepidosaurian lineage. An alternative explanation is that fissions of macrochromosomes, which increase microchromosomes, may also have occurred in the lineage leading to birds.
Intra-genomic heterogeneity of GC3 in snake
Relationships between GC-content and chromosome location of the snake genes
GC-rich gene (GC3≥50)
GC-poor gene (GC3≤50)
GC-rich isochores are known to have a clear association with R-bands (and particularly T bands) in mammals and birds [55–58]. In R-banded metaphases of rat snake, most microchromosomes showed R-positive bands (Additional file 5). However, R-positive bands also observed on most macrochromosomal regions. Thus, there was no clear correlation between R-band and GC3 of the mapped genes in rat snake.
Correlation between GC3 and GC-content of non-coding regions in snakes
The GC-content of the python non-coding regions showed a narrow distribution with a low average (37.9 ± 4.3%, mean ± standard deviation) in contrast to GC3 (Figure 3A and C). We divided the python non-coding regions into macrochromosomal (n = 139) and microchromosomal regions (n = 37) on the postulate that the chromosome locations of all orthologs were conserved between the rat snake and the python, and compared the GC-content between the two chromosomal groups (Figure 3D). The average GC-content of microchromosomal regions (42.5 ± 4.0%, mean ± standard deviation) were significantly higher than those of macrochromosomal regions (36.8 ± 3.5%) (Mann-Whitney’s U-test, P < 0.01).
Cross-species comparison of GC3 between orthologs of amniotes
We examined the frequency distribution of GC3 of orthologs in green anole lizard, chicken, Chinese soft-shelled turtle, human, mouse and an amphibian species X. tropicalis (Additional file 4). In previous study, we analyzed the frequency distribution of GC3 with massive genes in these species except for anole lizard. The distribution patterns of GC3 of this study were similar to those in Figure 2 of our previous study  although the gene set of this study contains a somewhat higher proportion of GC-poor genes. Thus we thought that GC3 of the gene set of this study could be used as a representative of GC3 of the whole genes in each species.
Comparison of gene functions between macro- and microchromosomes
Chromosome size-dependent GC heterogeneity is a pan-Sauropsida characteristic
We constructed a cytogenetic map of the Japanese four-striped rat snake, which contained 183 genes. Our cytogenetic map covered most macrochromosomal regions and at least one gene was mapped to each microchromosome on the basis of the homologies with chicken chromosomes. Our map also showed linkage homologies with most chicken chromosomes, namely 1–15, 17–22, 24–28 and Z. These results make it possible to infer the global GC heterogeneity of the snake genome and the shift of GC-content caused by chromosome rearrangements during sauropsid evolution.
We calculated GC3 to investigate the intra-genomic GC heterogeneity in the rat snake. The GC3 of the snake genes exhibited a bimodal distribution (Figure 2A). Such a bimodal distribution of GC3 was also observed in the genomes of the Chinese soft-shelled turtle (P. sinensis), chicken and non-rodent mammals . This result suggests that the GC3 heterogeneity is a common feature of amniote genomes. However, the standard deviation of GC3 in snake genes was somewhat lower than that of other amniotes , and GC-content of python non-coding regions also showed a narrow distribution, as observed in the green anole lizard . Thus the heterogeneities of base composition have probably decreased in lepidosaurian lineages over evolutionary time.
Although the standard deviation of GC3 of snake genes was relatively small, our results suggest that snake microchromosomes contain a higher proportion of GC-rich genes than macrochromosomes, as observed in both the Chinese soft-shelled turtle and chicken . Recently, chromosome size-dependent GC heterogeneity was also identified in the red-eared slider turtle (Trachemys scripta elegans) and the Nile crocodile (Crocodylus niloticus) using chromosome flow sorting technique . Chromosome size-dependent GC heterogeneity therefore seems to be a widespread characteristic in sauropsids whose karyotypes consist of macrochromosomes and microchromosomes, and possibly originated in the common ancestor of sauropsids. Interestingly, the green anole lizard, whose karyotype consists of 6 pairs of macrochromosomes and 12 pairs of microchromosomes [39, 62], does not show such marked biases in GC-content between macro- and microchromosomes. This suggests that the chromosome size-dependent GC heterogeneity has disappeared in the lineage leading to the anole lizard. Lepidosauria is a species-rich group consisting of about 8,000 species, and the karyotypes are also diversified within the group. Further investigation in various lepidosaurian species may help clarify the relationship between GC-content and the karyotype.
Disparity between lepidosaurs and the turtle-archosaurs
The correlation coefficient of GC3 between the rat snake and chicken orthologs are lower than that between the Chinese soft-shelled turtle and chicken (Figure 4). One explanation for the lower correlation is that the phylogenetic distance between the snake and chicken is larger than between the turtle and chicken. However, the divergence time between turtles and birds is estimated to be more than 231 million years, which is not largely different from the time of the lepidosauria-archosauria split, 275 million years ago [3, 7–9]. Therefore we consider the effect of the large differences of karyotypes, especially the number of microchromosomes, between the snake and the other two species.
The chromosome numbers are largely different between the rat snake (2n = 36) and chicken (2n = 78). In contrast, the karyotype of the Chinese soft-shelled turtle, which consists of nine pairs of macrochromosomes and 24 pairs of microchromosomes (2n = 66), is very similar to the chicken karyotype . Our previous study demonstrated by comparative gene mapping that the chromosomes have been highly conserved between chicken and the turtle [21, 38]. Chicken microchromosomes were considered to extensively retain the ancestral linkage groups of genes , and Nakatani et al. also suggested that many chicken microchromosomes (i.e., chromosomes 11, 15, 19, 20, 21, 22, 23, 24, 27 and 28) have one-to-one correspondence to ancestral proto-chromosomes of the gnathostome ancestor . These results lead us to infer that chromosome rearrangements have occurred more frequently in the snake lineage than in the chicken lineage, that is, chromosome number has been reduced by frequent chromosome fusions between macro- and microchromosomes and also between microchromosomes in the snake lineage.
Impact of chromosome fissions/fusions on GC-content
What mechanisms were involved in the changes of the GC-content of the genes after the fusion of microchromosomes into macrochromosomal complement? It has been suggested that the GC-content is primarily influenced by local recombination rates via GC-biased gene conversion [65, 66]. Under this model, A or T is displaced by G or C through mismatch repair when an AT/GC heteroduplex is formed at recombining regions. Accordingly, AT/GC heterozygotes produce more GC than AT gametes, thus conferring predominance of GC alleles in frequently recombining regions.
Recombination rate is negatively correlated with the size of chromosome arms in human and chicken [20, 67]. In chicken in particular, recombination rate per unit physical length is much higher in microchromosomes than in macrochromosomes . Recombination rates per physical length are thus expected to be lower in the snake macrochromosomal segments derived from the ancestral microchromosomes than those in their homologous chicken microchromosomes. The chromosome size-dependent difference in the recombination rate thus seems to have caused the decrease of GC-content in the snake macrochromosomal genes derived from the microchromosomes of the common ancestor of sauropsids.
Empirical evidence of chromosome length-driven evolution of GC-content has been shown in marsupial and monotreme species . In contrast, chromosome size-dependent GC heterogeneity was not clearly demonstrated in eutherian species, likely because frequent chromosome rearrangements have obscured the history of GC-content changes . In many eutherian lineages, however, intra-chromosomal GC heterogeneity has been reported [21, 36, 68, 69]. In chicken, intra-chromosomal GC heterogeneity is not as prominent as in eutherians , and there is no intra-chromosomal GC heterogeneity known in the anole genome [37, 39]. Our approach could not demonstrate whether there was intra-chromosomal GC heterogeneity in snake genome because of the insufficient sequence data provided by cDNA sequencing. It is necessary to conduct whole genome sequencing and assemble the sequences into a chromosome scale in order to clarify the relationship between intra- and inter-chromosomal GC heterogeneity in the snake genomes.
When overall GC-content and recombination rates were compared among various vertebrate species, there was no clear correlation between GC-content and recombination rates . For example, the recombination rate of chicken is about two times higher than that of zebra finch, although GC-content is almost equal between the two species [70, 71]. Whereas rat snake microchromosomal genes show similar level of GC3 to chicken microchromosomal genes, GC3 of rat snake macrochromosomal genes is significantly lower than chicken macrochromosomal genes (Figure 5). Thus other factors might exert influences on overall genomic GC-content. In bacteria, genomic GC-content have been subject to natural selection but not to biased gene conversion . Modes other than biased gene conversion were also proposed for evolution of genome composition in vertebrates . Further consideration is therefore necessary for the evolution of overall genome composition in snakes.
Any biological partitioning between macrochromosomes and microchromosomes?
Many literatures reported existing correlations between gene function and base compositions of the genes, the genomes and the promoter regions [73–75]. The difference of global GC-content between the macro- and microchromosomes may potentially cause the biased distribution of gene functions between the chromosomes: some proteins containing more amino acids for GC-rich codons due to functional constraints may be more advantageous in being encoded in microchromosomes than in macrochromosomes. To test this possibility, we investigated functional difference between chicken macro- and microchromosomal genes on the basis of the frequencies of appearance of functional categories (Figure 6). Over-representations found in six categories could be caused by clustering of particular gene families in short genomic stretches. These gene families may be less informative to test the hypothesis that genes are selected based on the size of chromosomes on which they are encoded. The over-representations of the other eight categories were independent from an excess of genes on particular chromosomes, implying that there could be some functional differences between macro- and microchromosomal genes.
In order to examine the relationship between the apparent localization of functional gene categories and the chromosome size-dependent GC-content, we investigated the GC-content of the gene sets assorted in the eight over-represented categories: GO:0004871, GO:0051234, GO:0006810, GO:0005215, IPR007110, GO:0044430, GO:0015075 and GO:0003774 (Additional file 6: Table S16). As shown in Figure 6 and Additional file 6: Table S1, seven of the eight categories were over-represented in microchromosomes. If most genes assorted in these GO or InterPro categories have high GC-content, it is likely that functional compartmentalization of genes is influenced by the difference of global GC-content between macro- and microchromosomes, i.e., some functionally categorized gene groups with high GC-content have been selectively preserved on microchromosomes. The average GC-content of the genes of the seven categories ranged from 43.9 to 48.1% and most of them were lower than the average of all chicken genes (47.4%) (Additional file 6: Table S16). In the other six categories (GO:0042611, GO:0044427, GO:0000785, GO:0000786, IPR003461 and IPR003597), there was no clear correlation between abundance in macro- or microchromosomes and GC-content. For an example, GO:0044427, GO:0000785 and GO:0000786 were over-represented in macrochromosomes, but their GC-contents were relatively high (Additional file 6: Table S1 and S16). These results imply that functional differences of genes did not correlate with the global difference of GC-content between macro- and microchromosomes. Further characterization of functional difference between macro- and microchromosomal genes, as well as its correlation with the general GC trend in chromosomal environments should await more extensive analyses using multiple species in future.
In this study, we constructed a cytogenetic map with 183 genes in the Japanese four-striped rat snake, and calculated GC3 across all chromosomes. Our results revealed cytogenetic evidence that snake microchromosomal genes tend to have higher GC3 than macrochromosomal genes, as found in the chicken and the Chinese soft-shelled turtle, a feature apparently lost in the genome of anole lizard. By comparing GC3 of orthologs between snake and chicken, we show that the GC-content of genes is correlated with the size of chromosomes on which the genes are harbored. This chromosome size-dependent GC heterogeneity is particularly apparent in snake genes that have been translocated from microchromosomes to macrochromosomes since snakes and birds shared a common ancestor, some 275 million years. The addition of whole genome sequencing and karyotypes from wide variety of sauropsidan species will provide the fine-scale picture of timing and mode of GC shift accompanying karyotypic evolution in this important group of vertebrates.
GC-content of exonic third codon positions
Expressed sequence tag
Fluorescence in situ hybridization
We thank Yasuyuki Hashiguchi for advice on identification of coding and non-coding regions in the python contigs, Yukuto Sato for advice on comparison of the gene functions, Yosuke Kawai for discussion and Denis O’meally for English correction. This work was supported by Grants-in-Aid for Scientific Research (No. 15370001, No. 16086201 and No. 19770001) from the Ministry of Education, Culture, Sports, Science and Technology in Japan, and by JSPS Research Fellowships for Young Scientists to K.M. (No. 215132).
- Zardoya R, Meyer A: Complete mitochondrial genome suggests diapsid affinities of turtles. Proc Natl Acad Sci USA. 1998, 95: 14226-14231. 10.1073/pnas.95.24.14226.PubMed CentralView ArticlePubMedGoogle Scholar
- Hedges SB, Poling LL: A molecular phylogeny of reptiles. Science. 1999, 283: 998-1001. 10.1126/science.283.5404.998.View ArticlePubMedGoogle Scholar
- Kumazawa Y, Nishida M: Complete mitochondrial DNA sequences of the green turtle and blue-tailed mole skink: statistical evidence for archosaurian affinity of turtles. Mol Biol Evol. 1999, 16: 784-792. 10.1093/oxfordjournals.molbev.a026163.View ArticlePubMedGoogle Scholar
- Cao Y, Sorenson MD, Kumazawa Y, Mindell DP, Hasegawa M: Phylogenetic position of turtles among amniotes: evidence from mitochondrial and nuclear genes. Gene. 2000, 259: 139-148. 10.1016/S0378-1119(00)00425-X.View ArticlePubMedGoogle Scholar
- Rest JS, Ast JC, Austin CC, Waddell PJ, Tibbetts EA, Hay JM, Mindell DP: Molecular systematics of primary reptilian lineages and the tuatara mitochondrial genome. Mol Phylogenet Evol. 2003, 29: 289-297. 10.1016/S1055-7903(03)00108-8.View ArticlePubMedGoogle Scholar
- Iwabe N, Hara Y, Kumazawa Y, Shibamoto K, Saito Y, Miyata T, Katoh K: Sister group relationship of turtles to the bird-crocodilian clade revealed by nuclear DNA-coded proteins. Mol Biol Evol. 2005, 22: 810-813. 10.1093/molbev/msi075.View ArticlePubMedGoogle Scholar
- Kumazawa Y: Mitochondrial genomes from major lizard families suggest their phylogenetic relationships and ancient radiations. Gene. 2007, 388: 19-26. 10.1016/j.gene.2006.09.026.View ArticlePubMedGoogle Scholar
- Benton MJ, Donoghue PCJ: Paleontological evidence to date the tree of life. Mol Biol Evol. 2007, 24: 26-53.View ArticlePubMedGoogle Scholar
- Shedlock AM, Edwards SV: Amniotes (Amniota). The timetree of life. Edited by: Hedges SB, Kumar S. 2009, New York: Oxford University Press, 375-379.Google Scholar
- Beçak W, Beçak ML, Nazareth HRS, Ohno S: Close karyological kinship between the reptilian suborder Serpentes and the class Aves. Chromosoma. 1964, 15: 606-617. 10.1007/BF00319994.View ArticlePubMedGoogle Scholar
- Bickham JW, Baker RJ: Chromosome homology and evolution of emydid turtles. Chromosoma. 1976, 54: 201-219. 10.1007/BF00293451.View ArticlePubMedGoogle Scholar
- Bickham JW, Bull JJ, Legler JM: Karyotypes and evolutionary relationships of trinychoid turtles. Cytologia. 1983, 48: 177-183. 10.1508/cytologia.48.177.View ArticleGoogle Scholar
- Norris TB, Rickards GK, Daugherty CH: Chromosomes of tuatara, Sphenodon, a chromosome heteromorphism and an archaic reptilian karyotype. Cytogenet Genome Res. 2004, 105: 93-99. 10.1159/000078014.View ArticlePubMedGoogle Scholar
- Takagi N, Sasaki M: A phylogenetic study of bird karyotypes. Chromosoma. 1974, 46: 91-120. 10.1007/BF00332341.View ArticlePubMedGoogle Scholar
- Belterman RHR, de Boer LEM: A karyological study of 55 species of birds, including karyotypes of 39 species new to cytology. Genetica. 1984, 65: 39-82. 10.1007/BF00056765.View ArticleGoogle Scholar
- Shetty S, Griffin DK, Graves JAM: Comparative painting reveals strong chromosome homology over 80 million years of bird evolution. Chromosome Res. 1999, 7: 289-295. 10.1023/A:1009278914829.View ArticlePubMedGoogle Scholar
- Nishida-Umehara C, Tsuda Y, Ishijima J, Ando J, Fujiwara A, Matsuda Y, Griffin DK: The molecular basis of chromosome orthologies and sex chromosomal differentiation in palaeognathous birds. Chromosome Res. 2007, 15: 721-734. 10.1007/s10577-007-1157-7.View ArticlePubMedGoogle Scholar
- Cohen MM, Gans C: The chromosomes of the order Crocodylia. Cytogenetics. 1970, 9: 81-105. 10.1159/000130080.View ArticlePubMedGoogle Scholar
- Kawagoshi T, Nishida C, Ota H, Kumazawa Y, Endo H, Matsuda Y: Molecular structures of centromeric heterochromatin and karyotypic evolution in the Siamese crocodile (Crocodylus siamensis) (Crocodylidae, Crocodylia). Chromosome Res. 2008, 16: 1119-1132. 10.1007/s10577-008-1263-1.View ArticlePubMedGoogle Scholar
- International Chicken Genome Sequencing Consortium (ICGSC): Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432: 695-716. 10.1038/nature03154.View ArticleGoogle Scholar
- Kuraku S, Ishijima J, Nishida-Umehara C, Agata K, Kuratani S, Yoichi M: cDNA-based gene mapping and GC3 profiling in the soft-shelled turtle suggest a chromosomal size-dependent GC bias shared by sauropsids. Chromosome Res. 2006, 14: 187-202. 10.1007/s10577-006-1035-8.View ArticlePubMedGoogle Scholar
- Costantini M, Filippo MD, Auletta F, Bernardi G: Isochore pattern and gene distribution in the chicken genome. Gene. 2007, 400: 9-15. 10.1016/j.gene.2007.05.025.View ArticlePubMedGoogle Scholar
- McQueen HA, Siriaco G, Bird AP: Chicken microchromosomes are hyperacetylated, early replicating, and gene rich. Genome Res. 1998, 8: 621-630.PubMed CentralPubMedGoogle Scholar
- Smith J, Bruley CK, Paton IR, Dunn I, Jones CT, Windsor D, Morrice DR, Law AS, Masabanda J, Sazanov A, Waddington D, Fries R, Burt DW: Differences in gene density on chicken macrochromosomes and microchromosomes. Anim Genet. 2000, 31: 96-103. 10.1046/j.1365-2052.2000.00565.x.View ArticlePubMedGoogle Scholar
- Habermann FA, Cremer M, Walter J, Kreth G, von Hase J, Bauer K, Wienberg J, Cremer C, Cremer T, Solovei I: Arrangements of macro- and microchromosomes in chicken cells. Chromosome Res. 2001, 9: 569-584. 10.1023/A:1012447318535.View ArticlePubMedGoogle Scholar
- Axelsson E, Webster MT, Smith NG, Burt DW, Ellegren H: Comparison of the chicken and turkey genomes reveals a higher rate of nucleotide divergence on microchromosomes than macrochromosomes. Genome Res. 2005, 15: 120-125. 10.1101/gr.3021305.PubMed CentralView ArticlePubMedGoogle Scholar
- Federico C, Cantarella CD, Scavo C, Saccone S, Bed'Hom B, Bernardi G: Avian genomes: different karyotypes but a similar distribution of the GC-richest chromosome regions at interphase. Chromosome Res. 2005, 13: 785-793. 10.1007/s10577-005-1012-7.View ArticlePubMedGoogle Scholar
- Shedlock AM, Botka CW, Zhao S, Shetty J, Zhang T, Liu JS, Deschavanne PJ, Edwards SV: Phylogenomics of nonavian reptiles and the structure of the ancestral amniote genome. Proc Natl Acad Sci USA. 2007, 104: 2767-2772. 10.1073/pnas.0606204104.PubMed CentralView ArticlePubMedGoogle Scholar
- Janes DE, Organ CL, Fujita MK, Shedlock AM, Edwards SV: Genome evolution in Reptilia, the sister groups of mammals. Annu Rev Genomics Hum Genet. 2010, 11: 239-264. 10.1146/annurev-genom-082509-141646.View ArticlePubMedGoogle Scholar
- Belle EM, Smith N, Eyre-Walker A: Analysis of the phylogenetic distribution of isochores in vertebrates and a test of the thermal stability hypothesis. J Mol Evol. 2002, 55: 356-363. 10.1007/s00239-002-2333-1.View ArticlePubMedGoogle Scholar
- Hamada K, Horiike T, Ota H, Mizuno K, Shinozawa T: Presence of isochore structures in reptile genomes suggested by the relationship between GC content of intron regions and those of coding regions. Genes Genet Syst. 2003, 78: 195-198. 10.1266/ggs.78.195.View ArticlePubMedGoogle Scholar
- Chojnowski JL, Franklin J, Katsu Y, Iguchi T, Guillette LJ, Kimball RT, Braun EL: Patterns of vertebrate isochore evolution revealed by comparison of expressed mammalian, avian, and crocodilian genes. J Mol Evol. 2007, 65: 259-266. 10.1007/s00239-007-9003-2.View ArticlePubMedGoogle Scholar
- Chojnowski JL, Braun EL: Turtle isochore structure is intermediate between amphibians and other amniotes. Integr Comp Biol. 2008, 48: 454-462. 10.1093/icb/icn062.View ArticlePubMedGoogle Scholar
- Elhaik E, Landan G, Graur D: Can GC content at third-codon positions be used as a proxy for isochore composition?. Mol Biol Evol. 2009, 26: 1829-1833. 10.1093/molbev/msp100.View ArticlePubMedGoogle Scholar
- Clay OK, Bernardi G: GC3 of genes can be used as a proxy for isochore base composition: A reply to Elhaik et al. Mol Biol Evol. 2011, 28: 21-23. 10.1093/molbev/msq222.View ArticlePubMedGoogle Scholar
- Romiguier J, Ranwez V, Douzery EJP, Galtier N: Contrasting GC-content dynamics across 33 mammalian genomes: Relationship with life-history traits and chromosome sizes. Genome Res. 2010, 20: 1001-1009. 10.1101/gr.104372.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Fujita MK, Edwards SV, Ponting CP: The Anolis lizard genome: an amniote genome without isochores. Genome Biol Evol. in press
- Matsuda Y, Nishida-Umehara C, Tarui H, Kuroiwa A, Yamada K, Isobe T, Ando J, Fujiwara A, Hirao Y, Nishimura O, Ishijima J, Hayashi A, Saito T, Murakami T, Murakami Y, Kuratani S, Agata K: Highly conserved linkage homology between birds and turtles: Bird and turtle chromosomes are precise counterparts of each other. Chromosome Res. 2005, 13: 601-615. 10.1007/s10577-005-0986-5.View ArticlePubMedGoogle Scholar
- Alföldi J, Di Palma F, Grabherr M, Williams C, Kong L, Mauceli E, Russell P, Lowe CB, Glor RE, Jaffe JD, Ray DA, Boissinot S, Shedlock AM, Botka C, Castoe TA, Colbourne JK, Fujita MK, Moreno RG, Ten Hallers BF, Haussler D, Heger A, Heiman D, Janes DE, Johnson J, de Jong PJ, Koriabine MY, Lara M, Novick PA, Organ CL, Peach SE, Poe S, Pollock DD, de Queiroz K, Sanger T, Searle S, Smith JD, Smith Z, Swofford R, Turner-Maier J, Wade J, Young S, Zadissa A, Edwards SV, Glenn TC, Schneider CJ, Losos JB, Lander ES, Breen M, Ponting CP, Lindblad-Toh K: The genome of the green anole lizard and a comparative analysis with birds and mammals. Nature. 2011, 477: 587-591. 10.1038/nature10390.PubMed CentralView ArticlePubMedGoogle Scholar
- Beçak W, Beçak ML: Cytotaxonomy and chromosomal evolution in Serpentes. Cytogenetics. 1969, 8: 247-262. 10.1159/000130037.View ArticlePubMedGoogle Scholar
- Singh L: Evolution of karyotypes in snakes. Chromosoma. 1972, 38: 185-236. 10.1007/BF00326193.View ArticlePubMedGoogle Scholar
- Matsubara K, Tarui H, Toriba M, Yamada K, Nishida-Umehara C, Agata K, Matsuda Y: Evidence for different origin of sex chromosomes in snakes, birds, and mammals and step-wise differentiation of snake sex chromosomes. Proc Natl Acad Sci USA. 2006, 103: 18190-18195. 10.1073/pnas.0605274103.PubMed CentralView ArticlePubMedGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMed CentralView ArticlePubMedGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.PubMedGoogle Scholar
- Katoh K, Kuma K, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33: 511-518. 10.1093/nar/gki198.PubMed CentralView ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52: 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Matsuda Y, Chapman VM: Application of fluorescence in situ hybridization in genome analysis of the mouse. Electrophoresis. 1995, 16: 261-272. 10.1002/elps.1150160142.View ArticlePubMedGoogle Scholar
- Stajich JE, Block D, Boulez K, Brenner SE, Chervitz SA, Dagdigian C, Fuellen G, Gilbert JG, Korf I, Lapp H, Lehväslaiho H, Matsalla C, Mungall CJ, Osborne BI, Pocock MR, Schattner P, Senger M, Stein LD, Stupka E, Wilkinson MD, Birney E: The Bioperl toolkit: Perl modules for the life sciences. Genome Res. 2002, 12: 1611-1618. 10.1101/gr.361602.PubMed CentralView ArticlePubMedGoogle Scholar
- Castoe TA, de Koning AJ, Hall KT, Yokoyama KD, Gu W, Smith EN, Feschotte C, Uetz P, Ray DA, Dobry J, Bogden R, Mackessy SP, Bronikowski AM, Warren WC, Secor SM, Pollock DD: Sequencing the genome of the Burmese python (Python molurus bivittatus) as a model for studying extreme adaptations in snakes. Genome Biol. 2011, 12: 406-10.1186/gb-2011-12-7-406. URL: http://www.ncbi.nlm.nih.gov/Traces/wgs/?val=AEQUPubMed CentralView ArticlePubMedGoogle Scholar
- Wise2 - Intelligent algorithms for DNA searches. URL: http://www.ebi.ac.uk/Tools/Wise2/index.html
- Hubbard T, Andrews D, Caccamo M, Cameron G, Chen Y, Clamp M, Clarke L, Coates G, Cox T, Cunningham F, Curwen V, Cutts T, Down T, Durbin R, Fernandez-Suarez XM, Gilbert J, Hammond M, Herrero J, Hotz H, Howe K, Iyer V, Jekosch K, Kahari A, Kasprzyk A, Keefe D, Keenan S, Kokocinsci F, London D, Longden I, McVicker G, Melsopp C, Meidl P, Potter S, Proctor G, Rae M, Rios D, Schuster M, Searle S, Severin J, Slater G, Smedley D, Smith J, Spooner W, Stabenau A, Stalker J, Storey R, Trevanion S, Ureta-Vidal A, Vogel J, White S, Woodwark C, Birney E: Ensembl 2005. Nucleic Acids Res. 2005, 33: D447-453. 10.1093/nar/gki378. URL: http://www.ensembl.org/PubMed CentralView ArticlePubMedGoogle Scholar
- NCBI Entrez Gene. URL: http://www.ncbi.nlm.nih.gov/gene/
- Al-Shahrour F, Dîaz-Uriarte R, Dopazo J: FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics. 2004, 20: 578-580. 10.1093/bioinformatics/btg455. URL: http://babelomics.bioinfo.cipf.es/functional.htmlView ArticlePubMedGoogle Scholar
- Shibusawa M, Nishibori M, Nishida-Umehara C, Tsudzuki M, Masabanda J, Griffin DK, Matsuda Y: Karyotypic evolution in the Galliformes: an examination of the process of karyotypic evolution by comparison of the molecular cytogenetic findings with the molecular phylogeny. Cytogenet Genome Res. 2004, 106: 111-119. 10.1159/000078570.View ArticlePubMedGoogle Scholar
- Holmquist GP: Chromosome bands, their chromatin flavors, and their functional features. Am J Hum Genet. 1992, 51: 17-37.PubMed CentralPubMedGoogle Scholar
- Saccone S, De Sario A, Della Valle G, Bernardi G: The highest gene concentrations in the human genome are in telomeric bands of metaphase chromosomes. Proc Natl Acad Sci USA. 1992, 89: 4913-4917. 10.1073/pnas.89.11.4913.PubMed CentralView ArticlePubMedGoogle Scholar
- Saccone S, De Sario A, Wiegant J, Raap AK, Della Valle G, Bernardi G: Correlations between isochores and chromosomal bands in the human genome. Proc Natl Acad Sci USA. 1993, 90: 11929-11933. 10.1073/pnas.90.24.11929.PubMed CentralView ArticlePubMedGoogle Scholar
- Saccone S, Cacciò S, Kusuda J, Andreozzi L, Bernardi G: Identification of the gene-richest bands in human chromosomes. Gene. 1996, 174: 85-94. 10.1016/0378-1119(96)00392-7.View ArticlePubMedGoogle Scholar
- Nabholz B, Künstner A, Wang R, Jarvis ED, Ellegren H: Dynamic evolution of base composition: causes and consequences in avian phylogenomics. Mol Biol Evol. 2011, 28: 2197-2210. 10.1093/molbev/msr047.PubMed CentralView ArticlePubMedGoogle Scholar
- Bucciarelli G, Di Filippo M, Costagliola D, Alvarez-Valin F, Bernardi G, Bernardi G: Environmental genomics: a tale of two fishes. Mol Biol Evol. 2009, 26: 1235-1243. 10.1093/molbev/msp041.View ArticlePubMedGoogle Scholar
- Kasai F, O'Brien PC, Ferguson-Smith MA: Reassessment of genome size in turtle and crocodile based on chromosome measurement by flow karyotyping: close similarity to chicken. Biol Lett. 2012, 8: 631-635. 10.1098/rsbl.2012.0141.PubMed CentralView ArticlePubMedGoogle Scholar
- Benirschke K, Hsu TC: Chromosome atlas: fish, amphibians, reptiles and birds volume 1. 1971, New York: SpringerGoogle Scholar
- Burt DW: Origin and evolution of avian microchromosomes. Cytogenet Genome Res. 2002, 96: 97-112. 10.1159/000063018.View ArticlePubMedGoogle Scholar
- Nakatani Y, Takeda H, Kohara Y, Morishita S: Reconstruction of the vertebrate ancestral genome reveals dynamic genome reorganization in early vertebrates. Genome Res. 2007, 17: 1254-1265. 10.1101/gr.6316407.PubMed CentralView ArticlePubMedGoogle Scholar
- Duret L, Eyre-Walker A, Galtier N: A new perspective on isochore evolution. Gene. 2006, 385: 71-74.View ArticlePubMedGoogle Scholar
- Duret L, Galtier N: Biased gene conversion and the evolution of mammalian genomic landscapes. Annu Rev Genomics Hum Genet. 2009, 10: 285-311. 10.1146/annurev-genom-082908-150001.View ArticlePubMedGoogle Scholar
- Meunier J, Duret L: Recombination drives the evolution of GC-content in the human genome. Mol Biol Evol. 2004, 21: 984-990. 10.1093/molbev/msh070.View ArticlePubMedGoogle Scholar
- Costantini M, Cammarano R, Bernardi G: The evolution of isochore patterns in vertebrate genomes. BMC Genomics. 2009, 10: 146-10.1186/1471-2164-10-146.PubMed CentralView ArticlePubMedGoogle Scholar
- Costantini M, Clay O, Auletta F, Bernardi1 G: An isochore map of human chromosomes. Genome Res. 2006, 16: 536-541. 10.1101/gr.4910606.PubMed CentralView ArticlePubMedGoogle Scholar
- Kai W, Kikuchi K, Tohari S, Chew AK, Tay A, Fujiwara A, Hosoya S, Suetake H, Naruse K, Brenner S, Suzuki Y, Venkatesh B: Integration of the genetic map and genome assembly of fugu facilitates insights into distinct features of genome evolution in teleosts and mammals. Genome Biol Evol. 2011, 3: 424-442. 10.1093/gbe/evr041.View ArticlePubMedGoogle Scholar
- Backström N, Forstmeier W, Schielzeth H, Mellenius H, Nam K, Bolund E, Webster MT, Ost T, Schneider M, Kempenaers B, Ellegren H: The recombination landscape of the zebra finch Taeniopygia guttata genome. Genome Res. 2010, 20: 485-495. 10.1101/gr.101410.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Hildebrand F, Meyer A, Eyre-Walker A: Evidence of selection upon genomic GC-content in bacteria. PLoS Genet. 2010, 6: e1001107-10.1371/journal.pgen.1001107.PubMed CentralView ArticlePubMedGoogle Scholar
- Bajic VB, Tan SL, Christoffels A, Schönbach C, Lipovich L, Yang L, Hofmann O, Kruger A, Hide W, Kai C, Kawai J, Hume DA, Carninci P, Hayashizaki Y: Mice and men: their promoter properties. PLoS Genet. 2006, 2: e54-10.1371/journal.pgen.0020054.PubMed CentralView ArticlePubMedGoogle Scholar
- D'Onofrio G, Ghosh TC, Saccone S: Different functional classes of genes are characterized by different compositional properties. FEBS Lett. 2007, 581: 5819-5824. 10.1016/j.febslet.2007.11.052.View ArticlePubMedGoogle Scholar
- Ren L, Gao G, Zhao D, Ding M, Luo J, Deng H: Developmental stage related patterns of codon usage and genomic GC content: searching for evolutionary fingerprints with models of stem cell differentiation. Genome Biol. 2007, 8: R35-10.1186/gb-2007-8-3-r35.PubMed CentralView ArticlePubMedGoogle Scholar