Intra-genomic GC heterogeneity in sauropsids: evolutionary insights from cDNA mapping and GC(3) profiling in snake

Insights into avian molecular cytogenetics-with reptilian comparisons

In last 100 years or so, much information has been accumulated on avian karyology, genetics, physiology, biochemistry and evolution. The chicken genome project generated genomic resources used in comparative studies, elucidating fundamental evolutionary processes, much of it funded by the economic importance of domestic fowl (which are also excellent model species in many areas). Studying karyotypes and whole genome sequences revealed population processes, evolutionary biology, and genome function, uncovering the role of repetitive sequences, transposable elements and gene family expansion. Knowledge of the function of many genes and non-expressed or identified regulatory components is however still lacking. Birds (Aves) are diverse, have striking adaptations for flight, migration and survival and inhabit all continents most islands. They also have a unique karyotype with ~ 10 macrochromosomes and ~ 30 microchromosomes that are smaller than other reptiles. Classified into Palaeognathae and Neognathae they are evolutionarily close, and a subset of reptiles. Here we overview avian molecular cytogenetics with reptilian comparisons, shedding light on their karyotypes and genome structure features. We consider avian evolution, then avian (followed by reptilian) karyotypes and genomic features. We consider synteny disruptions, centromere repositioning, and repetitive elements before turning to comparative avian and reptilian genomics. In this context, we review comparative cytogenetics and genome mapping in birds as well as Z- and W-chromosomes and sex determination. Finally, we give examples of pivotal research areas in avian and reptilian cytogenomics, particularly physical mapping and map integration of sex chromosomal genes, comparative genomics of chicken, turkey and zebra finch, California condor cytogenomics as well as some peculiar cytogenetic and evolutionary examples. We conclude that comparative molecular studies and improving resources continually contribute to new approaches in population biology, developmental biology, physiology, disease ecology, systematics, evolution and phylogenetic systematics orientation. This also produces genetic mapping information for chromosomes active in rearrangements during the course of evolution. Further insights into mutation, selection and adaptation of vertebrate genomes will benefit from these studies including physical and online resources for the further elaboration of comparative genomics approaches for many fundamental biological questions.

Rapid Targeted Assembly of the Proteome Reveals Evolutionary Variation of GC Content in Avian Lice

Nucleotide base composition plays an influential role in the molecular mechanisms involved in gene function, phenotype, and amino acid composition. GC content (proportion of guanine and cytosine in DNA sequences) shows a high level of variation within and among species. Many studies measure GC content in a small number of genes, which may not be representative of genome-wide GC variation. One challenge when assembling extensive genomic data sets for these studies is the significant amount of resources (monetary and computational) associated with data processing, and many bioinformatic tools have not been optimized for resource efficiency. Using a high-performance computing (HPC) cluster, we manipulated resources provided to the targeted gene assembly program, automated target restricted assembly method (aTRAM), to determine an optimum way to run the program to maximize resource use. Using our optimum assembly approach, we assembled and measured GC content of all of the protein-coding genes of a diverse group of parasitic feather lice. Of the 499 426 genes assembled across 57 species, feather lice were GC-poor (mean GC = 42.96%) with a significant amount of variation within and between species (GC range = 19.57%-73.33%). We found a significant correlation between GC content and standard deviation per taxon for overall GC and GC3, which could indicate selection for G and C nucleotides in some species. Phylogenetic signal of GC content was detected in both GC and GC3. This research provides a large-scale investigation of GC content in parasitic lice laying the foundation for understanding the basis of variation in base composition across species.

GC Content Across Insect Genomes: Phylogenetic Patterns, Causes and Consequences

The proportions of A:T and G:C nucleotide pairs are often unequal and can vary greatly between animal species and along chromosomes. The causes and consequences of this variation are incompletely understood. The recent release of high-quality genome sequences from the Darwin Tree of Life and other large-scale genome projects provides an opportunity for GC heterogeneity to be compared across a large number of insect species. Here we analyse GC content along chromosomes, and within protein-coding genes and codons, of 150 insect species from four holometabolous orders: Coleoptera, Diptera, Hymenoptera, and Lepidoptera. We find that protein-coding sequences have higher GC content than the genome average, and that Lepidoptera generally have higher GC content than the other three insect orders examined. GC content is higher in small chromosomes in most Lepidoptera species, but this pattern is less consistent in other orders. GC content also increases towards subtelomeric regions within protein-coding genes in Diptera, Coleoptera and Lepidoptera. Two species of Diptera, Bombylius major and B. discolor, have very atypical genomes with ubiquitous increase in AT content, especially at third codon positions. Despite dramatic AT-biased codon usage, we find no evidence that this has driven divergent protein evolution. We argue that the GC landscape of Lepidoptera, Diptera and Coleoptera genomes is influenced by GC-biased gene conversion, strongest in Lepidoptera, with some outlier taxa affected drastically by counteracting processes.

Elasmobranch genome sequencing reveals evolutionary trends of vertebrate karyotype organization

Genomic studies of vertebrate chromosome evolution have long been hindered by the scarcity of chromosome-scale DNA sequences of some key taxa. One of those limiting taxa has been the elasmobranchs (sharks and rays), which harbor species often with numerous chromosomes and enlarged genomes. Here, we report the chromosome-scale genome assembly for the zebra shark Stegostoma tigrinum, an endangered species that has a relatively small genome among sharks (3.71 Gb), as well as for the whale shark Rhincodon typus Our analysis using a male-female comparison identified an X Chromosome, the first genomically characterized shark sex chromosome. The X Chromosome harbors the Hox C cluster whose intact linkage has not been shown for an elasmobranch fish. The sequenced shark genomes show a gradualism of chromosome length with remarkable length-dependent characteristics-shorter chromosomes tend to have higher GC content, gene density, synonymous substitution rate, and simple tandem repeat content as well as smaller gene length and lower interspersed repeat content. We challenge the traditional binary classification of karyotypes as with and without so-called microchromosomes. Even without microchromosomes, the length-dependent characteristics persist widely in nonmammalian vertebrates. Our investigation of elasmobranch karyotypes underpins their unique characteristics and provides clues for understanding how vertebrate karyotypes accommodate intragenomic heterogeneity to realize a complex readout. It also paves the way to dissecting more genomes with variable sizes to be sequenced at high quality.

Why Do Some Vertebrates Have Microchromosomes?

With more than 70,000 living species, vertebrates have a huge impact on the field of biology and research, including karyotype evolution. One prominent aspect of many vertebrate karyotypes is the enigmatic occurrence of tiny and often cytogenetically indistinguishable microchromosomes, which possess distinctive features compared to macrochromosomes. Why certain vertebrate species carry these microchromosomes in some lineages while others do not, and how they evolve remain open questions. New studies have shown that microchromosomes exhibit certain unique characteristics of genome structure and organization, such as high gene densities, low heterochromatin levels, and high rates of recombination. Our review focuses on recent concepts to expand current knowledge on the dynamic nature of karyotype evolution in vertebrates, raising important questions regarding the evolutionary origins and ramifications of microchromosomes. We introduce the basic karyotypic features to clarify the size, shape, and morphology of macro- and microchromosomes and report their distribution across different lineages. Finally, we characterize the mechanisms of different evolutionary forces underlying the origin and evolution of microchromosomes.

Quantitative Approach to Fish Cytogenetics in the Context of Vertebrate Genome Evolution

Our novel Python-based tool EVANGELIST allows the visualization of GC and repeats percentages along chromosomes in sequenced genomes and has enabled us to perform quantitative large-scale analyses on the chromosome level in fish and other vertebrates. This is a different approach from the prevailing analyses, i.e., analyses of GC% in the coding sequences that make up not more than 2% in human. We identified GC content (GC%) elevations in microchromosomes in ancient fish lineages similar to avian microchromosomes and a large variability in the relationship between the chromosome size and their GC% across fish lineages. This raises the question as to what extent does the chromosome size drive GC% as posited by the currently accepted explanation based on the recombination rate. We ascribe the differences found across fishes to varying GC% of repetitive sequences. Generally, our results suggest that the GC% of repeats and proportion of repeats are independent of the chromosome size. This leaves an open space for another mechanism driving the GC evolution in vertebrates.

Snake W Sex Chromosome: The Shadow of Ancestral Amniote Super-Sex Chromosome

: Heteromorphic sex chromosomes, particularly the ZZ/ZW sex chromosome system of birds and some reptiles, undergo evolutionary dynamics distinct from those of autosomes. The W sex chromosome is a unique karyological member of this heteromorphic pair, which has been extensively studied in snakes to explore the origin, evolution, and genetic diversity of amniote sex chromosomes. The snake W sex chromosome offers a fascinating model system to elucidate ancestral trajectories that have resulted in genetic divergence of amniote sex chromosomes. Although the principal mechanism driving evolution of the amniote sex chromosome remains obscure, an emerging hypothesis, supported by studies of W sex chromosomes of squamate reptiles and snakes, suggests that sex chromosomes share varied genomic blocks across several amniote lineages. This implies the possible split of an ancestral super-sex chromosome via chromosomal rearrangements. We review the major findings pertaining to sex chromosomal profiles in amniotes and discuss the evolution of an ancestral super-sex chromosome by collating recent evidence sourced mainly from the snake W sex chromosome analysis. We highlight the role of repeat-mediated sex chromosome conformation and present a genomic landscape of snake Z and W chromosomes, which reveals the relative abundance of major repeats, and identifies the expansion of certain transposable elements. The latest revolution in chromosomics, i.e., complete telomere-to-telomere assembly, offers mechanistic insights into the evolutionary origin of sex chromosomes.

Genome Complexity Reduction High-Throughput Genome Sequencing of Green Iguana (Iguana iguana) Reveal a Paradigm Shift in Understanding Sex-Chromosomal Linkages on Homomorphic X and Y Sex Chromosomes

The majority of lizards classified in the superfamily Iguanoidea have an XX/XY sex-determination system in which sex-chromosomal linkage shows homology with chicken (Gallus gallus) chromosome 15 (GGA15). However, the genomics of sex chromosomes remain largely unexplored owing to the presence of homomorphic sex chromosomes in majority of the species. Recent advances in high-throughput genome complexity reduction sequencing provide an effective approach to the identification of sex-specific loci with both single-nucleotide polymorphisms (SNPs) and restriction fragment presence/absence (PA), and a better understanding of sex chromosome dynamics in Iguanoidea. In this study, we applied Diversity Arrays Technology (DArTseq^TM) in 29 phenotypic sex assignments (14 males and 15 females) of green iguana (Iguana iguana). We confirmed a male heterogametic (XX/XY) sex determination mode in this species, identifying 29 perfectly sex-linked SNP/PA loci and 164 moderately sex-linked SNP/PA loci, providing evidence probably indicative of XY recombination. Three loci from among the perfectly sex-linked SNP/PA loci showed partial homology with several amniote sex chromosomal linkages. The results support the hypothesis of an ancestral super-sex chromosome with overlaps of partial sex-chromosomal linkages. However, only one locus among the moderately sex-linked loci showed homology with GGA15, which suggests that the specific region homologous to GGA15 was located outside the non-recombination region but in close proximity to this region of the sex chromosome in green iguana. Therefore, the location of GGA15 might be further from the putative sex-determination locus in green iguana. This is a paradigm shift in understanding linkages on homomorphic X and Y sex chromosomes. The DArTseq platform provides an easy-to-use strategy for future research on the evolution of sex chromosomes in Iguanoidea, particularly for non-model species with homomorphic or highly cryptic sex chromosomes.

Genome-wide SNP analysis of Siamese cobra (Naja kaouthia) reveals the molecular basis of transitions between Z and W sex chromosomes and supports the presence of an ancestral super-sex chromosome in amniotes

Elucidation of the process of sex chromosome differentiation is necessary to understand the dynamics of evolutionary mechanisms in organisms. The W sex chromosome of the Siamese cobra (Naja kaouthia) contains a large number of repeats and shares amniote sex chromosomal linkages. Diversity Arrays Technology provides an effective approach to identify sex-specific loci that are epoch-making, to understand the dynamics of molecular transitions between the Z and W sex chromosomes in a snake lineage. From a total of 543 sex-specific loci, 90 showed partial homology with sex chromosomes of several amniotes and 89 loci were homologous to transposable elements. Two loci were confirmed as W-specific nucleotides after PCR amplification. These loci might result from a sex chromosome differentiation process and involve putative sex-determination regions in the Siamese cobra. Sex-specific loci shared linkage homologies among amniote sex chromosomes, supporting an ancestral super-sex chromosome.

Partial Amniote Sex Chromosomal Linkage Homologies Shared on Snake W Sex Chromosomes Support the Ancestral Super-Sex Chromosome Evolution in Amniotes

Squamate reptile chromosome 2 (SR2) is thought to be an important remnant of an ancestral amniote super-sex chromosome, but a recent study showed that the Siamese cobra W sex chromosome is also a part of this larger ancestral chromosome. To confirm the existence of an ancestral amniote super-sex chromosome and understand the mechanisms of amniote sex chromosome evolution, chromosome maps of two snake species [Russell's viper: Daboia russelii (DRU) and the common tiger snake: Notechis scutatus (NSC)] were constructed using bacterial artificial chromosomes (BACs) derived from chicken and zebra finch libraries containing amniote sex chromosomal linkages. Sixteen BACs were mapped on the W sex chromosome of DRU and/or NSC, suggesting that these BACs contained a common genomic region shared with the W sex chromosome of these snakes. Two of the sixteen BACs were co-localized to DRU2 and NSC2, corresponding to SR2. Prediction of genomic content from all BACs mapped on snake W sex chromosomes revealed a large proportion of long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) retrotransposons. These results led us to predict that amplification of LINE and SINE may have occurred on snake W chromosomes during evolution. Genome compartmentalization, such as transposon amplification, might be the key factor influencing chromosome structure and differentiation. Multiple sequence alignments of all BACs mapped on snake W sex chromosomes did not reveal common sequences. Our findings indicate that the SR2 and snake W sex chromosomes may have been part of a larger ancestral amniote super-sex chromosome, and support the view of sex chromosome evolution as a colorful myriad of situations and trajectories in which many diverse processes are in action.

Consequence of Paradigm Shift with Repeat Landscapes in Reptiles: Powerful Facilitators of Chromosomal Rearrangements for Diversity and Evolution

Reptiles are notable for the extensive genomic diversity and species richness among amniote classes, but there is nevertheless a need for detailed genome-scale studies. Although the monophyletic amniotes have recently been a focus of attention through an increasing number of genome sequencing projects, the abundant repetitive portion of the genome, termed the "repeatome", remains poorly understood across different lineages. Consisting predominantly of transposable elements or mobile and satellite sequences, these repeat elements are considered crucial in causing chromosomal rearrangements that lead to genomic diversity and evolution. Here, we propose major repeat landscapes in representative reptilian species, highlighting their evolutionary dynamics and role in mediating chromosomal rearrangements. Distinct karyotype variability, which is typically a conspicuous feature of reptile genomes, is discussed, with a particular focus on rearrangements correlated with evolutionary reorganization of micro- and macrochromosomes and sex chromosomes. The exceptional karyotype variation and extreme genomic diversity of reptiles are used to test several hypotheses concerning genomic structure, function, and evolution.

Molecular cytogenetic characterization of repetitive sequences comprising centromeric heterochromatin in three Anseriformes species

The highly repetitive DNA sequence of centromeric heterochromatin is an effective molecular cytogenetic marker for investigating genomic compartmentalization between macrochromosomes and microchromosomes in birds. We isolated four repetitive sequence families of centromeric heterochromatin from three Anseriformes species, viz., domestic duck (Anas platyrhynchos, APL), bean goose (Anser fabalis, AFA), and whooper swan (Cygnus cygnus, CCY), and characterized the sequences by molecular cytogenetic approach. The 190-bp APL-HaeIII and 101-bp AFA-HinfI-S sequences were localized in almost all chromosomes of A. platyrhynchos and A. fabalis, respectively. However, the 192-bp AFA-HinfI-L and 290-bp CCY-ApaI sequences were distributed in almost all microchromosomes of A. fabalis and in approximately 10 microchromosomes of C. cygnus, respectively. APL-HaeIII, AFA-HinfI-L, and CCY-ApaI showed partial sequence homology with the chicken nuclear-membrane-associated (CNM) repeat families, which were localized primarily to the centromeric regions of microchromosomes in Galliformes, suggesting that ancestral sequences of the CNM repeat families are observed in the common ancestors of Anseriformes and Galliformes. These results collectively provide the possibility that homogenization of centromeric heterochromatin occurred between microchromosomes in Anseriformes and Galliformes; however, homogenization between macrochromosomes and microchromosomes also occurred in some centromeric repetitive sequences.

Distribution and amplification of interstitial telomeric sequences (ITSs) in Australian dragon lizards support frequent chromosome fusions in Iguania

Telomeric sequences are generally located at the ends of chromosomes; however, they can also be found in non-terminal chromosomal regions when they are known as interstitial telomeric sequences (ITSs). Distribution of ITSs across closely related and divergent species elucidates karyotype evolution and speciation as ITSs provide evolutionary evidence for chromosome fusion. In this study, we performed physical mapping of telomeric repeats by fluorescence in situ hybridisation (FISH) in seven Australian dragon lizards thought to represent derived karyotypes of squamate reptiles and a gecko lizard with considerably different karyotypic feature. Telomeric repeats were present at both ends of all chromosomes in all species, while varying numbers of ITSs were also found on microchromosomes and in pericentromeric or centromeric regions on macrochromosomes in five lizard species examined. This suggests that chromosomal rearrangements from ancestral squamate reptiles to Iguania occurred mainly by fusion between ancestral types of acrocentric chromosomes and/or between microchromosomes, leading to appearance of bi-armed macrochromosomes, and in the reduction of microchromosome numbers. These results support the previously proposed hypothesis of karyotype evolution in squamate reptiles. In addition, we observed the presence of telomeric sequences in the similar regions to heterochromatin of the W microchromosome in Pogona barbata and Doporiphora nobbi, while sex chromosomes for the two species contained part of the nucleolar organiser regions (NORs). This likely implies that these ITSs are a part of the satellite DNA and not relics of chromosome fusions. Amplification of telomeric repeats may have involved heterochromatinisation of sex-specific W chromosomes and play a role in the organisation of the nucleolus.

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?

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.

The habu genome reveals accelerated evolution of venom protein genes

Evolution of novel traits is a challenging subject in biological research. Several snake lineages developed elaborate venom systems to deliver complex protein mixtures for prey capture. To understand mechanisms involved in snake venom evolution, we decoded here the ~1.4-Gb genome of a habu, Protobothrops flavoviridis. We identified 60 snake venom protein genes (SV) and 224 non-venom paralogs (NV), belonging to 18 gene families. Molecular phylogeny reveals early divergence of SV and NV genes, suggesting that one of the four copies generated through two rounds of whole-genome duplication was modified for use as a toxin. Among them, both SV and NV genes in four major components were extensively duplicated after their diversification, but accelerated evolution is evident exclusively in the SV genes. Both venom-related SV and NV genes are significantly enriched in microchromosomes. The present study thus provides a genetic background for evolution of snake venom composition.

Physical mapping of repetitive DNA suggests 2n reduction in Amazon turtles Podocnemis (Testudines: Podocnemididae)

Cytogenetic studies show that there is great karyotypic diversity in order Testudines (2n = 26-68), and that this may be mainly attributed to the presence/absence of microchromosomes. Members of the Podocnemididae family have the smallest diploid numbers of this order (2n = 26-28), which may be a derived condition of the group. Diverse studies suggest that repetitive-DNA-rich sites generally act as hotspots for double-strand breaks and chromosomal reorganization. In this context, we used fluorescent in situ hybridization (FISH) to map telomeric sequences (TTAGGG)n, 45S rDNA, and the genes encoding histones H1 and H3 in two species of genus Podocnemis. We also observed conservation of the 45S rDNA and H1 histone sequences (probable case of conserved synteny), but multiple conserved and non-conserved clusters of H3 genes, which colocalized with the interstitial telomeric sequences in the Podocnemis genome. Our results suggest that fusions have occurred between macro and microchromosomes or between microchromosomes, leading to the observed reduction in diploid number in the family Podocnemididae.

Evolutionary Dynamics of the Gametologous CTNNB1 Gene on the Z and W Chromosomes of Snakes

Snakes exhibit genotypic sex determination with female heterogamety (ZZ males and ZW females), and the state of sex chromosome differentiation also varies among lineages. To investigate the evolutionary history of homologous genes located in the nonrecombining region of differentiated sex chromosomes in snakes, partial sequences of the gametologous CTNNB1 gene were analyzed for 12 species belonging to henophid (Cylindrophiidae, Xenopeltidae, and Pythonidae) and caenophid snakes (Viperidae, Elapidae, and Colubridae). Nonsynonymous/synonymous substitution ratios (Ka/Ks) in coding sequences were low (Ka/Ks < 1) between CTNNB1Z and CTNNB1W, suggesting that these 2 genes may have similar functional properties. However, frequencies of intron sequence substitutions and insertion–deletions were higher in CTNNB1Z than CTNNB1W, suggesting that Z-linked sequences evolved faster than W-linked sequences. Molecular phylogeny based on both intron and exon sequences showed the presence of 2 major clades: 1) Z-linked sequences of Caenophidia and 2) W-linked sequences of Caenophidia clustered with Z-linked sequences of Henophidia, which suggests that the sequence divergence between CTNNB1Z and CTNNB1W in Caenophidia may have occurred by the cessation of recombination after the split from Henophidia.

A novel method of caenophidian snake sex identification using molecular markers based on two gametologous genes

Sex identification provides important information for ecological and evolutionary studies, as well as benefiting snake conservation management. Traditional methods such as cloacal probing or cloacal popping are counterproductive for sex identification concerning very small species, resulting in difficulties in the management of their breeding programs. In this study, the nucleotide sequences of gametologous genes (CTNNB1 and WAC genes) were used for the development of molecular sexing markers in caenophidian snakes. Two candidate markers were developed with the two primer sets, and successfully amplified by a single band on the agarose gel in male (ZZ) and two bands, differing in fragment sizes, in female (ZW) of 16 caenophidian snakes for CTNNB1 and 12 caenophidian snakes for WAC. Another candidate marker was developed with the primer set to amplify the specific sequence for CTNNB1W homolog, and the PCR products were successfully obtained in a female‐specific 250‐bp DNA bands. The three candidate PCR sexing markers provide a simple sex identification method based on the amplification of gametologous genes, and they can be used to facilitate effective caenophidian snake conservation and management programs.

Sex chromosome evolution in snakes inferred from divergence patterns of two gametologous genes and chromosome distribution of sex chromosome-linked repetitive sequences

BACKGROUND

The discovery of differentially organized sex chromosome systems suggests that heteromorphic sex chromosomes evolved from a pair of homologous chromosomes. Whereas karyotypes are highly conserved in alethinophidian snakes, the degeneration status of the W chromosomes varies among species. The Z and W chromosomes are morphologically homomorphic in henophidian species, whereas in snakes belonging to caenophidian families the W chromosomes are highly degenerated. Snakes therefore are excellent animal models in which to study sex chromosome evolution. Herein, we investigated the differentiation processes for snake sex chromosomes using both coding and repetitive sequences. We analyzed phylogenetic relationships of CTNNB1 and WAC genes, localized to the centromeric and telomeric regions, respectively, of the long arms on snake sex chromosomes, and chromosome distribution of sex chromosome-linked repetitive sequences in several henophidian and caenophidian species.

RESULTS

Partial or full-length coding sequences of CTNNB1 and WAC were identified for Z homologs of henophidian species from Tropidophiidae, Boidae, Cylindrophiidae, Xenopeltidae, and Pythonidae, and for Z and W homologs of caenophidian species from Acrochordidae, Viperidae, Elapidae, and Colubridae. Female-specific sequences for the two genes were not found in the henophidian (boid and pythonid) species examined. Phylogenetic trees constructed using each gene showed that the Z and W homologs of the caenophidian species cluster separately. The repetitive sequence isolated from the W chromosome heterochromatin of the colubrid Elaphe quadrivirgata and a microsatellite motif (AGAT)8 were strongly hybridized with W chromosomes of the viperid and colubrid species examined.

CONCLUSION

Our phylogenetic analyses suggest that the cessation of recombination between the Z and W homologs of CTNNB1 and WAC predated the diversification of the caenophidian families. As the repetitive sequences on the W chromosomes were shared among viperid and colubrid species, heterochromatinization of the proto-W chromosome appears to have occurred before the splitting of these two groups. These results collectively suggest that differentiation of the proto-Z and proto-W chromosomes extended to wide regions on the sex chromosomes in the common ancestor of caenophidian families during a relatively short period.

Anchoring genome sequence to chromosomes of the central bearded dragon (Pogona vitticeps) enables reconstruction of ancestral squamate macrochromosomes and identifies sequence content of the Z chromosome

Background

Squamates (lizards and snakes) are a speciose lineage of reptiles displaying considerable karyotypic diversity, particularly among lizards. Understanding the evolution of this diversity requires comparison of genome organisation between species. Although the genomes of several squamate species have now been sequenced, only the green anole lizard has any sequence anchored to chromosomes. There is only limited gene mapping data available for five other squamates. This makes it difficult to reconstruct the events that have led to extant squamate karyotypic diversity. The purpose of this study was to anchor the recently sequenced central bearded dragon (Pogona vitticeps) genome to chromosomes to trace the evolution of squamate chromosomes. Assigning sequence to sex chromosomes was of particular interest for identifying candidate sex determining genes.

Results

By using two different approaches to map conserved blocks of genes, we were able to anchor approximately 42 % of the dragon genome sequence to chromosomes. We constructed detailed comparative maps between dragon, anole and chicken genomes, and where possible, made broader comparisons across Squamata using cytogenetic mapping information for five other species. We show that squamate macrochromosomes are relatively well conserved between species, supporting findings from previous molecular cytogenetic studies. Macrochromosome diversity between members of the Toxicofera clade has been generated by intrachromosomal, and a small number of interchromosomal, rearrangements. We reconstructed the ancestral squamate macrochromosomes by drawing upon comparative cytogenetic mapping data from seven squamate species and propose the events leading to the arrangements observed in representative species. In addition, we assigned over 8 Mbp of sequence containing 219 genes to the Z chromosome, providing a list of genes to begin testing as candidate sex determining genes.

Conclusions

Anchoring of the dragon genome has provided substantial insight into the evolution of squamate genomes, enabling us to reconstruct ancestral macrochromosome arrangements at key positions in the squamate phylogeny, demonstrating that fusions between macrochromosomes or fusions of macrochromosomes and microchromosomes, have played an important role during the evolution of squamate genomes. Assigning sequence to the sex chromosomes has identified NR5A1 as a promising candidate sex determining gene in the dragon.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2774-3) contains supplementary material, which is available to authorized users.

Role of Chromosome Changes in Crocodylus Evolution and Diversity

The karyotypes of most species of crocodilians were studied using conventional and molecular cytogenetics. These provided an important contribution of chromosomal rearrangements for the evolutionary processes of Crocodylia and Sauropsida (birds and reptiles). The karyotypic features of crocodilians contain small diploid chromosome numbers (30~42), with little interspecific variation of the chromosome arm number (fundamental number) among crocodiles (56~60). This suggested that centric fusion and/or fission events occurred in the lineage, leading to crocodilian evolution and diversity. The chromosome numbers of Alligator, Caiman, Melanosuchus, Paleosuchus, Gavialis, Tomistoma, Mecistops, and Osteolaemus were stable within each genus, whereas those of Crocodylus (crocodylians) varied within the taxa. This agreed with molecular phylogeny that suggested a highly recent radiation of Crocodylus species. Karyotype analysis also suggests the direction of molecular phylogenetic placement among Crocodylus species and their migration from the Indo-Pacific to Africa and The New World. Crocodylus species originated from an ancestor in the Indo-Pacific around 9~16 million years ago (MYA) in the mid-Miocene, with a rapid radiation and dispersion into Africa 8~12 MYA. This was followed by a trans-Atlantic dispersion to the New World between 4~8 MYA in the Pliocene. The chromosomes provided a better understanding of crocodilian evolution and diversity, which will be useful for further study of the genome evolution in Crocodylia.

Immunofluorescent staining reveals hypermethylation of microchromosomes in the central bearded dragon, Pogona vitticeps

Background

Studies of model organisms have demonstrated that DNA cytosine methylation and histone modifications are key regulators of gene expression in biological processes. Comparatively little is known about the presence and distribution of epigenetic marks in non-model amniotes such as non-avian reptiles whose genomes are typically packaged into chromosomes of distinct size classes. Studies of chicken karyotypes have associated the gene-richness and high GC content of microchromosomes with a distinct epigenetic landscape. To determine whether this is likely to be a common feature of amniote microchromosomes, we have analysed the distribution of epigenetic marks using immunofluorescence on metaphase chromosomes of the central bearded dragon (Pogona vitticeps). This study is the first to study the distribution of epigenetic marks on non-avian reptile chromosomes.

Results

We observed an enrichment of DNA cytosine methylation, active modifications H3K4me2 and H3K4me3, as well as the repressive mark H3K27me3 in telomeric regions on macro and microchromosomes. Microchromosomes were hypermethylated compared to macrochromosomes, as they are in chicken. However, differences between macro- and microchromosomes for histone modifications associated with actively transcribed or repressed DNA were either less distinct or not detectable.

Conclusions

Hypermethylation of microchromosomes compared to macrochromosomes is a shared feature between P. vitticeps and avian species. The lack of the clear distinction between macro- and microchromosome staining patterns for active and repressive histone modifications makes it difficult to determine at this stage whether microchrosome hypermethylation is correlated with greater gene density as it is in aves, or associated with the greater GC content of P. vitticeps microchromosomes compared to macrochromosomes.

Karyotype Reorganization in the Hokou Gecko (Gekko hokouensis, Gekkonidae): The Process of Microchromosome Disappearance in Gekkota

The Hokou gecko (Gekko hokouensis: Gekkonidae, Gekkota, Squamata) has the chromosome number 2n = 38, with no microchromosomes. For molecular cytogenetic characterization of the gekkotan karyotype, we constructed a cytogenetic map for G. hokouensis, which retains the ancestral karyotype of Gekkota, with 86 functional genes, and compared it with cytogenetic maps for four Toxicofera species that have many microchromosomes (Elaphe quadrivirgata, Varanus salvator macromaculatus, Leiolepis reevesii rubritaeniata, and Anolis carolinensis) and that for a lacertid species (Lacerta agilis) with only one pair of autosomal microchromosomes. Ten pairs of G. hokouensis chromosomes [GHO1, 2, 3, Z(4), 6, 7, 8, 13, 14, and 15] showed highly conserved linkage homology with macrochromosomes and/or macrochromosome arms of the four Toxicofera species and corresponded to eight L. agilis macrochromosomes (LAG). However, GHO5, GHO9, GHO10, GHO11, and LAG6 were composed of chromosome segments that have a homology with Toxicofera microchromosomes, and no homology was found in the chromosomes between G. hokouensis and L. agilis. These results suggest that repeated fusions of microchromosomes may have occurred independently in each lineage of Gekkota and Lacertidae, leading to the disappearance of microchromosomes and appearance of small-sized macrochromosomes.

Biased gene conversion and GC-content evolution in the coding sequences of reptiles and vertebrates

Mammalian and avian genomes are characterized by a substantial spatial heterogeneity of GC-content, which is often interpreted as reflecting the effect of local GC-biased gene conversion (gBGC), a meiotic repair bias that favors G and C over A and T alleles in high-recombining genomic regions. Surprisingly, the first fully sequenced nonavian sauropsid (i.e., reptile), the green anole Anolis carolinensis, revealed a highly homogeneous genomic GC-content landscape, suggesting the possibility that gBGC might not be at work in this lineage. Here, we analyze GC-content evolution at third-codon positions (GC3) in 44 vertebrates species, including eight newly sequenced transcriptomes, with a specific focus on nonavian sauropsids. We report that reptiles, including the green anole, have a genome-wide distribution of GC3 similar to that of mammals and birds, and we infer a strong GC3-heterogeneity to be already present in the tetrapod ancestor. We further show that the dynamic of coding sequence GC-content is largely governed by karyotypic features in vertebrates, notably in the green anole, in agreement with the gBGC hypothesis. The discrepancy between third-codon positions and noncoding DNA regarding GC-content dynamics in the green anole could not be explained by the activity of transposable elements or selection on codon usage. This analysis highlights the unique value of third-codon positions as an insertion/deletion-free marker of nucleotide substitution biases that ultimately affect the evolution of proteins.

Tracing the evolution of amniote chromosomes

A great deal of diversity in chromosome number and arrangement is observed across the amniote phylogeny. Understanding how this diversity is generated is important for determining the role of chromosomal rearrangements in generating phenotypic variation and speciation. Gaining this understanding is achieved by reconstructing the ancestral genome arrangement based on comparisons of genome organization of extant species. Ancestral karyotypes for several amniote lineages have been reconstructed, mainly from cross-species chromosome painting data. The availability of anchored whole genome sequences for amniote species has increased the evolutionary depth and confidence of ancestral reconstructions from those made solely from chromosome painting data. Nonetheless, there are still several key lineages where the appropriate data required for ancestral reconstructions is lacking. This review highlights the progress that has been made towards understanding the chromosomal changes that have occurred during amniote evolution and the reconstruction of ancestral karyotypes.

Highly species-specific centromeric repetitive DNA sequences in lizards: molecular cytogenetic characterization of a novel family of satellite DNA sequences isolated from the water monitor lizard (Varanus salvator macromaculatus, Platynota)

Two novel repetitive DNA sequences, VSAREP1 and VSAREP2, were isolated from the water monitor lizard (Varanus salvator macromaculatus, Platynota) and characterized using molecular cytogenetics. The respective lengths and guanine-cytosine (GC) contents of the sequences were 190 bp and 57.5% for VSAREP1 and 185 bp and 59.7% for VSAREP2, and both elements were tandemly arrayed as satellite DNA in the genome. VSAREP1 and VSAREP2 were each located at the C-positive heterochromatin in the pericentromeric region of chromosome 2q, the centromeric region of chromosome 5, and 3 pairs of microchromosomes. This suggests that genomic compartmentalization between macro- and microchromosomes might not have occurred in the centromeric repetitive sequences of V. salvator macromaculatus. These 2 sequences did only hybridize to genomic DNA of V. salvator macromaculatus, but no signal was observed even for other squamate reptiles, including Varanus exanthematicus, which is a closely related species of V. salvator macromaculatus. These results suggest that these sequences were differentiated rapidly or were specifically amplified in the V. salvator macromaculatus genome.

Genomic structures of the kW1 loci on the Z and W chromosomes in ratite birds: structural changes at an early stage of W chromosome differentiation

The W chromosome of ratite birds shows minimal morphological differentiation and retains homology of genetic linkage and gene order with a substantial stretch of the Z chromosome; however, the molecular structure in the differentiated region is still not well known. The kW1 sequence was isolated from the kiwi as a W-specific DNA marker for PCR-based molecular sexing of ratite birds. In ratite W chromosomes, this sequence commonly contains a ∼200-bp deletion. To characterize the very early event of avian sex chromosome differentiation, we performed molecular cytogenetic analyses of kW1 and its flanking sequences in paleognathous and neognathous birds and reptiles. Female-specific repeats were found in the kW1-flanking sequence of the cassowary (Casuarius casuarius), and the repeats have been amplified in the pericentromeric region of the W chromosomes of ratites, which may have resulted from the cessation of meiotic recombination between the Z and W chromosomes at an early stage of sex chromosome differentiation. The presence of the kW1 sequence in neognathous birds and a crocodilian species suggests that the kW1 sequence was present in the ancestral genome of Archosauria; however, it disappeared in other reptilian taxa and several lineages of neognathous birds.

Dynamic chromosome reorganization in the osprey ( Pandion haliaetus , Pandionidae, Falconiformes): relationship between chromosome size and the chromosomal distribution of centromeric repetitive DNA sequences

The osprey (Pandion haliaetus) has a diploid number of 74 chromosomes, consisting of a large number of medium-sized macrochromosomes and relatively few microchromosomes; this differs greatly from the typical avian karyotype. Chromosome painting with chicken DNA probes revealed that the karyotype of P. haliaetus differs from the chicken karyotype by at least 14 fission events involving macrochromosomes (chicken chromosomes 1-9 and Z) and at most 15 fusions of microchromosomes, suggesting that considerable karyotype reorganization occurred in P. haliaetus in a similar manner previously reported for Accipitridae. A distinct difference was observed, however, between Accipitridae and Pandionidae with respect to the pattern of chromosome rearrangements that occurred after fissions of macrochromosomes. Metacentric or submetacentric chromosomes 1-5 in P. haliaetus appear to have been formed by centric fusion of chromosome segments derived from macrochromosomal fissions. By contrast, many pairs of bi-armed chromosomes in Accipitridae species seem to result from pericentric inversions that occurred in the fission-derived chromosomes. Two families of repetitive sequences were isolated; the 173-bp PHA-HaeIII sequence occurred on all chromosomes, whereas intense signals from the 742-bp PHA-NsiI sequence were localized to all acrocentric chromosomes, with weak signals on most of the bi-armed chromosomes. Two repetitive sequences cohybridized in the centromeric heterochromatin; however, the sequences differed in unit size, nucleotide sequence and GC content. The results suggest that the 2 sequence families originated from different ancestral sequences and were homogenized independently in centromeres, and that a chromosome size-dependent compartmentalization may have been lost in P. haliaetus.

Karyotype evolution in monitor lizards: cross-species chromosome mapping of cDNA reveals highly conserved synteny and gene order in the Toxicofera clade

The water monitor lizard (Varanus salvator macromaculatus (VSA), Platynota) has a chromosome number of 2n = 40: its karyotype consists of 16 macrochromosomes and 24 microchromosomes. To delineate the process of karyotype evolution in V. salvator macromaculatus, we constructed a cytogenetic map with 86 functional genes and compared it with those of the butterfly lizard (Leiolepis reevesii rubritaeniata (LRE); 2n = 36) and Japanese four-striped rat snake (Elaphe quadrivirgata (EQU); 2n = 36), members of the Toxicofera clade. The syntenies and gene orders of macrochromosomes were highly conserved between these species except for several chromosomal rearrangements: eight pairs of VSA macrochromosomes and/or chromosome arms exhibited homology with six pairs of LRE macrochromosomes and eight pairs of EQU macrochromosomes. Furthermore, the genes mapped to microchromosomes of three species were all located on chicken microchromosomes or chromosome 4p. No reciprocal translocations were found in the species, and their karyotypic differences were caused by: low frequencies of interchromosomal rearrangements, such as tandem fusions, or centric fissions/fusions between macrochromosomes and between macro- and microchromosomes; and intrachromosomal rearrangements, such as paracentric inversions or centromere repositioning. The chromosomal rearrangements that occurred in macrochromosomes of the Varanus lineage were also identified through comparative cytogenetic mapping of V. salvator macromaculatus and V. exanthematicus. Morphologic differences in chromosomes 6-8 between the two species could have resulted from pericentric inversion or centromere repositioning.

Molecular cytogenetics in artificial hybrid and highly polyploid sturgeons: An evolutionary story narrated by repetitive sequences

We applied comparative genomic hybridization (CGH) and genomic in situ hybridization (GISH) to examine genomes of artificially produced sturgeon hybrids between sterlet, Acipenser ruthenus female (∼120 chromosomes) or Russian sturgeon, A. gueldenstaedtii female (∼240 chromosomes) and a spontaneous triploid Siberian sturgeon A. baerii male (∼360 chromosomes), respectively. The ploidy levels of progenies were analyzed by karyotyping and flow cytometry. We found that the species-specific regions were surprisingly identifiable only on some micro- and small(er) macrochromosomes in hybrid metaphases. We hypothesize that these distinguishable regions are represented by species-specific repetitive sequences driven by more dynamic molecular evolutionary mechanisms. On larger chromosomes, GISH faintly visualized only blocks of pericentromeric and telomeric repetitive sequences, remaining regions were equally shared by both parental species. We concluded that the interspecies hybridization producing viable and even fertile progeny is enabled by the fact that genomes of the species involved are likely divergent at the level of the repetitive sequences only and probably highly conserved in the coding sequences. These small differences of coding sequences are in concordance with previous estimations of relatedness of examined species producing artificial as well as natural hybrids. CGH and GISH represent a challenge in sturgeon cytogenetics as a valuable though technically not simple tool to discriminate chromosomes of parental species in hybrids. The potentials and drawbacks of CGH and GISH application in sturgeons are discussed and further experimental possibilities are proposed.

Molecular cytogenetic map of the central bearded dragon, Pogona vitticeps (Squamata: Agamidae)

Reptiles, as the sister group to birds and mammals, are particularly valuable for comparative genomic studies among amniotes. The Australian central bearded dragon (Pogona vitticeps) is being developed as a reptilian model for such comparisons, with whole-genome sequencing near completion. The karyotype consists of 6 pairs of macrochromosomes and 10 pairs microchromosomes (2n = 32), including a female heterogametic ZW sex microchromosome pair. Here, we present a molecular cytogenetic map for P. vitticeps comprising 87 anchor bacterial artificial chromosome clones that together span each macro- and microchromosome. It is the first comprehensive cytogenetic map for any non-avian reptile. We identified an active nucleolus organizer region (NOR) on the sub-telomeric region of 2q by mapping 18S rDNA and Ag-NOR staining. We identified interstitial telomeric sequences in two microchromosome pairs and the W chromosome, indicating that microchromosome fusion has been a mechanism of karyotypic evolution in Australian agamids within the last 21 to 19 million years. Orthology searches against the chicken genome revealed an intrachromosomal rearrangement of P. vitticeps 1q, identified regions orthologous to chicken Z on P. vitticeps 2q, snake Z on P. vitticeps 6q and the autosomal microchromosome pair in P. vitticeps orthologous to turtle Pelodiscus sinensis ZW and lizard Anolis carolinensis XY. This cytogenetic map will be a valuable reference tool for future gene mapping studies and will provide the framework for the work currently underway to physically anchor genome sequences to chromosomes for this model Australian squamate.

Inference of the protokaryotypes of amniotes and tetrapods and the evolutionary processes of microchromosomes from comparative gene mapping

Comparative genome analysis of non-avian reptiles and amphibians provides important clues about the process of genome evolution in tetrapods. However, there is still only limited information available on the genome structures of these organisms. Consequently, the protokaryotypes of amniotes and tetrapods and the evolutionary processes of microchromosomes in tetrapods remain poorly understood. We constructed chromosome maps of functional genes for the Chinese soft-shelled turtle (Pelodiscus sinensis), the Siamese crocodile (Crocodylus siamensis), and the Western clawed frog (Xenopus tropicalis) and compared them with genome and/or chromosome maps of other tetrapod species (salamander, lizard, snake, chicken, and human). This is the first report on the protokaryotypes of amniotes and tetrapods and the evolutionary processes of microchromosomes inferred from comparative genomic analysis of vertebrates, which cover all major non-avian reptilian taxa (Squamata, Crocodilia, Testudines). The eight largest macrochromosomes of the turtle and chicken were equivalent, and 11 linkage groups had also remained intact in the crocodile. Linkage groups of the chicken macrochromosomes were also highly conserved in X. tropicalis, two squamates, and the salamander, but not in human. Chicken microchromosomal linkages were conserved in the squamates, which have fewer microchromosomes than chicken, and also in Xenopus and the salamander, which both lack microchromosomes; in the latter, the chicken microchromosomal segments have been integrated into macrochromosomes. Our present findings open up the possibility that the ancestral amniotes and tetrapods had at least 10 large genetic linkage groups and many microchromosomes, which corresponded to the chicken macro- and microchromosomes, respectively. The turtle and chicken might retain the microchromosomes of the amniote protokaryotype almost intact. The decrease in number and/or disappearance of microchromosomes by repeated chromosomal fusions probably occurred independently in the amphibian, squamate, crocodilian, and mammalian lineages.