Different genomic evolutionary rates in the various reptile lineages

Evolution of ancient satellite DNAs in extant alligators and caimans (Crocodylia, Reptilia)

BACKGROUND

Crocodilians are one of the oldest extant vertebrate lineages, exhibiting a combination of evolutionary success and morphological resilience that has persisted throughout the history of life on Earth. This ability to endure over such a long geological time span is of great evolutionary importance. Here, we have utilized the combination of genomic and chromosomal data to identify and compare the full catalogs of satellite DNA families (satDNAs, i.e., the satellitomes) of 5 out of the 8 extant Alligatoridae species. As crocodilian genomes reveal ancestral patterns of evolution, by employing this multispecies data collection, we can investigate and assess how satDNA families evolve over time.

RESULTS

Alligators and caimans displayed a small number of satDNA families, ranging from 3 to 13 satDNAs in A. sinensis and C. latirostris, respectively. Together with little variation both within and between species it highlighted long-term conservation of satDNA elements throughout evolution. Furthermore, we traced the origin of the ancestral forms of all satDNAs belonging to the common ancestor of Caimaninae and Alligatorinae. Fluorescence in situ experiments showed distinct hybridization patterns for identical orthologous satDNAs, indicating their dynamic genomic placement.

CONCLUSIONS

Alligators and caimans possess one of the smallest satDNA libraries ever reported, comprising only four sets of satDNAs that are shared by all species. Besides, our findings indicated limited intraspecific variation in satellite DNA, suggesting that the majority of new satellite sequences likely evolved from pre-existing ones.

Cytogenetic Insights into the Evolution of Chromosomes and Sex Determination Reveal Striking Homology of Turtle Sex Chromosomes to Amphibian Autosomes

Turtle karyotypes are highly conserved compared to other vertebrates; yet, variation in diploid number (2n = 26-68) reflects profound genomic reorganization, which correlates with evolutionary turnovers in sex determination. We evaluate the published literature and newly collected comparative cytogenetic data (G- and C-banding, 18S-NOR, and telomere-FISH mapping) from 13 species spanning 2n = 28-68 to revisit turtle genome evolution and sex determination. Interstitial telomeric sites were detected in multiple lineages that underwent diploid number and sex determination turnovers, suggesting chromosomal rearrangements. C-banding revealed potential interspecific variation in centromere composition and interstitial heterochromatin at secondary constrictions. 18S-NORs were detected in secondary constrictions in a single chromosomal pair per species, refuting previous reports of multiple NORs in turtles. 18S-NORs are linked to ZW chromosomes in Apalone and Pelodiscus and to X (not Y) in Staurotypus. Notably, comparative genomics across amniotes revealed that the sex chromosomes of several turtles, as well as mammals and some lizards, are homologous to components of Xenopus tropicalis XTR1 (carrying Dmrt1). Other turtle sex chromosomes are homologous to XTR4 (carrying Wt1). Interestingly, all known turtle sex chromosomes, except in Trionychidae, evolved via inversions around Dmrt1 or Wt1. Thus, XTR1 appears to represent an amniote proto-sex chromosome (perhaps linked ancestrally to XTR4) that gave rise to turtle and other amniote sex chromosomes.

Physical Mapping and Refinement of the Painted Turtle Genome (Chrysemys picta) Inform Amniote Genome Evolution and Challenge Turtle-Bird Chromosomal Conservation

Comparative genomics continues illuminating amniote genome evolution, but for many lineages our understanding remains incomplete. Here, we refine the assembly (CPI 3.0.3 NCBI AHGY00000000.2) and develop a cytogenetic map of the painted turtle (Chrysemys picta-CPI) genome, the first in turtles and in vertebrates with temperature-dependent sex determination. A comparison of turtle genomes with those of chicken, selected nonavian reptiles, and human revealed shared and novel genomic features, such as numerous chromosomal rearrangements. The largest conserved syntenic blocks between birds and turtles exist in four macrochromosomes, whereas rearrangements were evident in these and other chromosomes, disproving that turtles and birds retain fully conserved macrochromosomes for greater than 300 Myr. C-banding revealed large heterochromatic blocks in the centromeric region of only few chromosomes. The nucleolar-organizing region (NOR) mapped to a single CPI microchromosome, whereas in some turtles and lizards the NOR maps to nonhomologous sex-chromosomes, thus revealing independent translocations of the NOR in various reptilian lineages. There was no evidence for recent chromosomal fusions as interstitial telomeric-DNA was absent. Some repeat elements (CR1-like, Gypsy) were enriched in the centromeres of five chromosomes, whereas others were widespread in the CPI genome. Bacterial artificial chromosome (BAC) clones were hybridized to 18 of the 25 CPI chromosomes and anchored to a G-banded ideogram. Several CPI sex-determining genes mapped to five chromosomes, and homology was detected between yet other CPI autosomes and the globally nonhomologous sex chromosomes of chicken, other turtles, and squamates, underscoring the independent evolution of vertebrate sex-determining mechanisms.

Complete mitochondrial genome of Cuora trifasciata (Chinese three-striped box turtle), and a comparative analysis with other box turtles

Cuora trifasciata has become one of the most critically endangered species in the world. The complete mitochondrial genome of C. trifasciata (Chinese three-striped box turtle) was determined in this study. Its mitochondrial genome is a 16,575-bp-long circular molecule that consists of 37 genes that are typically found in other vertebrates. And the basic characteristics of the C. trifasciata mitochondrial genome were also determined. Moreover, a comparison of C. trifasciata with Cuora cyclornata, Cuora pani and Cuora aurocapitata indicated that the four mitogenomics differed in length, codons, overlaps, 13 protein-coding genes (PCGs), ND3, rRNA genes, control region, and other aspects. Phylogenetic analysis with Bayesian inference and maximum likelihood based on 12 protein-coding genes of the genus Cuora indicated the phylogenetic position of C. trifasciata within Cuora. The phylogenetic analysis also showed that C. trifasciata from Vietnam and China formed separate monophyletic clades with different Cuora species. The results of nucleotide base compositions, protein-coding genes and phylogenetic analysis showed that C. trifasciata from these two countries may represent different Cuora species.

Compositional properties and thermal adaptation of 18S rRNA in vertebrates

In order to investigate the influence of temperature on the GC level of the paired sequences of ribosomal 18S RNAs in vertebrates, we have studied their base composition in cold- and warm-blooded vertebrates using a stem-by-stem comparison. We observed that a number of stems of 18S ribosomal RNAs (rRNAs) are variable among species and that the majority of such stems are GC richer in warm-blooded than in cold-blooded vertebrates. We also constructed the secondary structures of the 18S rRNAs of a polar fish, a marsupial, and a monotreme to compare them with those of temperate/tropical fishes and of eutherians, respectively. In these cases, differences similar to those already mentioned were found. We conclude that there is a correlation between stem stability and body temperature even within the relatively limited temperature range of vertebrates.

Three tiers of genome evolution in reptiles

Characterization of reptilian genomes is essential for understanding the overall diversity and evolution of amniote genomes, because reptiles, which include birds, constitute a major fraction of the amniote evolutionary tree. To better understand the evolution and diversity of genomic characteristics in Reptilia, we conducted comparative analyses of online sequence data from Alligator mississippiensis (alligator) and Sphenodon punctatus (tuatara) as well as genome size and karyological data from a wide range of reptilian species. At the whole-genome and chromosomal tiers of organization, we find that reptilian genome size distribution is consistent with a model of continuous gradual evolution while genomic compartmentalization, as manifested in the number of microchromosomes and macrochromosomes, appears to have undergone early rapid change. At the sequence level, the third genomic tier, we find that exon size in Alligator is distributed in a pattern matching that of exons in Gallus (chicken), especially in the 101-200 bp size class. A small spike in the fraction of exons in the 301 bp-1 kb size class is also observed for Alligator, but more so for Sphenodon. For introns, we find that members of Reptilia have a larger fraction of introns within the 101 bp-2 kb size class and a lower fraction of introns within the 5-30 kb size class than do mammals. These findings suggest that the mode of reptilian genome evolution varies across three hierarchical levels of the genome, a pattern consistent with a mosaic model of genomic evolution.

Trends in the evolution of reptilian chromosomes

Reptiles are a karyologically heterogeneous group, where some orders and suborders exhibit characteristics similar to those of anamniotes and others share similarities with homeotherms. The class also shows different evolutionary trends, for instance in genome and chromosome size and composition. The turtle DNA base composition is similar to that of mammals, whereas that of lizards and snakes is more similar to that of anamniotes. The major karyological differences between turtles and squamates are the size and composition of the genome and the rate at which chromosomes change. Turtles have larger and more variable genome sizes, and a greater amount of middle repetitive DNA that differs even among related species. In lizards and snakes size of the genome are smaller, single-copy DNA is constant within each suborder, and differences in repetitive DNA involve fractions that become increasingly heterogeneous with widening phylogenetic distance. With regard to variation in karyotype morphology, turtles and crocodiles show low variability in chromosome number, morphology, and G-banding pattern. Greater variability is found among squamates, which have a similar degree of karyotypic change-as do some mammals, such as carnivores and bats-and in which there are also differences among congeneric species. An interesting relationship has been highlighted in the entire class Reptilia between rates of change in chromosomes, number of living species, and rate of extinction. However, different situations obtain in turtles and crocodiles on the one hand, and squamates on the other. In the former, the rate of change in chromosomes is lower and the various evolutionary steps do not seem to have entailed marked chromosomal variation, whereas squamates have a higher rate of change in chromosomes clearly related to the number of living species, and chromosomal variation seems to have played an important role in the evolution of several taxa. The different evolutionary trends in chromosomes observed between turtles and crocodiles on the one hand and squamates on the other might depend on their different patterns of G-banding.

Molecular cloning of natriuretic peptides from the heart of reptiles: loss of ANP in diapsid reptiles and birds

Atrial natriuretic peptide (ANP) and B-type NP (BNP) are hormones involved in homeostatic control of body fluid and cardiovascular regulation. Both ANP and BNP have been cloned from the heart of mammals, amphibians, and teleost fishes, while an additional cardiac peptide, ventricular NP, has been found in selected species of teleost fish. However, in chicken, BNP is the primary cardiac peptide identified thus far. In contrast, the types of NP/s present in the reptilian heart are unknown, representing a considerable gap in our understanding of NP evolution. In the present study, we cloned and sequenced a BNP cDNA from the atria of representative species of reptile, including crocodile, lizard, snake, and tortoise. In addition, we cloned BNP from the pigeon atria. The reptilian and pigeon BNP cDNAs had ATTTA repeats in the 3' untranslated region, as observed in all vertebrate BNP mRNAs. A high sequence homology was evident when comparing reptile and pigeon preproBNP with the previously identified chicken preproBNP. In particular, the predicted mature BNP-29 was identical between crocodile, tortoise, and chicken, with pigeon having a single amino acid substitution; lizard and snake BNP had seven and nine substitutions, respectively. Furthermore, an ANP cDNA could only be cloned from the tortoise atria. Since ANP was not isolated from the heart of any non-chelonian reptile and appears to be absent in birds, we propose that the ANP gene has been lost after branching of the turtles in the amniote line. This data provides new avenues for research on NP function in reptiles.

Characteristics of neural and humoral systems involved in the regulation of blood pressure in snakes

Cardiovascular function is affected by many mechanisms, including the autonomic system, the kallikrein-kinin system (KKS), the renin-angiotensin system (RAS) and the endothelin system. The function of these systems seems to be fairly well preserved throughout the vertebrate scale, but evolution required several adaptations. Snakes are particularly interesting for studies related to the cardiovascular function because of their elongated shape, their wide variation in size and length, and because they had to adapt to extremely different habitats and gravitational influences. To keep the normal cardiovascular control the snakes developed anatomical and functional adaptations and interesting structural peculiarities are found in their autonomic, KKS, RAS and endothelin systems. Our laboratory has characterized some biochemical, pharmacological and physiological properties of these systems in South American snakes. This review compares the components and function of these systems in snakes and other vertebrates, and focuses on differences found in snakes, related with receptor or ligand structure and/or function in autonomic system, RAS and KKS, absence of components in KKS and the intriguing identity between a venom and a plasma component in the endothelin system.

Karyotypic diversity and speciation in Agrodiaetus butterflies

That chromosomal rearrangements may play an important role in maintaining postzygotic isolation between well-established species is part of the standard theory of speciation. However, little evidence exists on the role of karyotypic change in speciation itself--in the establishment of reproductive barriers between previously interbreeding populations. The large genus Agrodiaetus (Lepidoptera: Lycaenidae) provides a model system to study this question. Agrodiaetus butterflies exhibit unusual interspecific diversity in chromosome number, from n= 10 to n= 134; in contrast, the majority of lycaenid butterflies have n= 23/24. We analyzed the evolution of karyotypic diversity by mapping chromosome numbers on a thoroughly sampled mitochondrial phylogeny of the genus. Karyotypic differences accumulate gradually between allopatric sister taxa, but more rapidly between sympatric sister taxa. Overall, sympatric sister taxa have a higher average karyotypic diversity than allopatric sister taxa. Differential fusion of diverged populations may account for this pattern because the degree of karyotypic difference acquired between allopatric populations may determine whether they will persist as nascent biological species in secondary sympatry. This study therefore finds evidence of a direct role for chromosomal rearrangements in the final stages of animal speciation. Rapid karyotypic diversification is likely to have contributed to the explosive speciation rate observed in Agrodiaetus, 1.6 species per million years.

DNA methylation in reptiles

Very recent investigations have provided evidence for a higher DNA methylation level in polar and sub-antarctic fishes compared to temperate/tropical fishes, the latter being in turn higher than the DNA methylation level of warm-blooded vertebrates. These results confirm and extend the finding [Jabbari, K., Cacciò, S., Pais de Barros, J.P., Desgres, J., Bernardi G., 1997. Evolutionary changes in CpG and methylation levels in the genome of vertebrates. Gene 205, 109-118] that DNA methylation level of vertebrates is inversely related to body temperature. Here we studied the methylation level of reptilian genomes. The species previously analyzed exhibited methylation levels closer to those of mammals and birds rather than to those of fishes and amphibians. The sample was, however, too small to reach a final conclusion. Here we used Reversed-Phase-High-Performance Liquid Chromatography (RP-HPLC) to analyze the DNA methylation levels of 43 reptiles representing three out of four orders and 20 families. Such analysis has shown that snakes and lizards exhibit methylation levels covering the whole range comprised between those of temperate/tropical fish and mammals, while turtles, and, more so, crocodiles are close to mammals. We discuss some ecological and physiological data that explain these results.

Rate of Chromosome changes and Speciation in Reptiles

The chromosome changing rate (i.e. the number of chromosome rearrangements per million years) was studied in 1,329 reptile species in order to evaluate the karyological evolutionary trend and the existence of possible correlations between chromosome mutations and some aspects of the evolution of this class. The results obtained highlight the existence of a general direct correlation between chromosome changing rate and number of living species, although different trends can be observed in the different orders and suborders. In turtles, the separation of pleurodires from cryptodires was accompanied by a considerable karyological diversification. Among pleurodires, the evolution of the Chelidae and Pelomedusidae was also characterised by chromosome variation, while in cryptodires a marked karyological homogeneity is observed between and within infra-orders. Similarly there is no correlation between changing rate and species number in crocodiles, where the evolution of the families and genera has entailed few chromosome mutations. Chromosome variability was greater in lizards and snakes. In the formers variations in chromosome changing rate accompanied the separation of the infra-orders and the evolution of most of the families and of some genera. The origin of snakes has also been accompanied by a marked karyological diversification, while the subsequent evolution of the infra-orders and families has entailed a high level of chromosome variability only in colubroids. The karyological evolution in reptiles generally entailed a progressive reduction in chromosome changing rate, albeit with differences in the diverse orders and suborders. This trend seems to be consistent with the "canalization model" as originally proposed by Bickham and Baker in [Bickham, J.W. & R J. Baker, 1979. Bull. Carnegie Mus. Nat. Hist. 13: 70-84.] However, several inconsistencies have been found excluding that in this class the ultimate goal of chromosome variations was the achievement of a so-called "optimum karyotype'' as suggested by the above-mentioned theory. Other mechanisms could underpin chromosome variability in Reptiles. Among them a genomic composition more or less favourable to promoting chromosome rearrangements and factors favouring the fixation of a mutant karyotype in condition of homozygosis. Turtles and crocodiles would have a genome characterised by large chromosomes and a low level of chromosome compartmentalisation limiting the recombination and the frequency of rearrangements. A low rate of chromosome variability modifying little if at all the gene linkage groups would have favoured a conservative evolutionary strategy. In the course of evolution, lizards and snakes could have achieved a genome characterised by smaller chromosomes and a higher level of compartmentalisation. This would have raised the frequency of recombination and consequently an evolutionary strategy promoting a higher degree of variability and a greater level of speciation.

Chromosomes of tuatara, Sphenodon, a chromosome heteromorphism and an archaic reptilian karyotype

We examined karyotypes of the endemic New Zealand reptile genus Sphenodon (tuatara) from five populations, finding a karyotype unchanged for at least one million years. Animals karyotyped were from five geographically distinct populations, representing three groups, namely S. guntheri, S. punctatus (Cook Strait group), and S. punctatus (northeastern North Island group). All five populations have a diploid chromosome number of 2n = 36, consisting of 14 pairs of macrochromosomes and four pairs of microchromosomes. Chromosomal differences were not found between the five populations nor between female and male animals, except for one animal with a structural heteromorphism. Similarity between Sphenodon and Testudine karyotypes suggests an ancestral karyotype with a macrochromosome complement of 14 pairs and the ability to accumulate variable numbers of microchromosome pairs. Our research supports molecular phylogenies of the Reptilia.

Molecular cloning, sequence analysis and expression of the snake follicle-stimulating hormone receptor

Follicle-stimulating hormone (FSH) and luteinizing hormone (LH) control gonadal function in mammalian and many non-mammalian vertebrates through the interaction with their receptors, FSHR and LHR. Although the same is true for some reptilian species, in Squamata (lizards and snakes) there is no definitive evidence for the presence of either two distinct gonadotropins or two distinct gonadotropin receptors. Our aim was to characterize the gonadotropin receptor(s) of the Bothrops jararaca snake. Using a cDNA library from snake testis and amplification of the 5'-cDNA ending, we cloned a cDNA related to FSHR. Attempts to clone a cDNA more closely related to LHR were unsuccessful. Expression of FSHR mRNA was restricted to gonadal tissues. The snake FSHR is a G protein-coupled receptor with 673 amino acids, and the aminoterminal domain with 346 amino acids consists of a nine leucine-rich repeat-containing subdomain (LRR) flanked by two cysteine-rich subdomains. The beta-strands in the LRR are conserved with exception of the third, a region that may be important for FSH binding. In contrast with mammalian, avian and amphibian FSHRs, the snake FSHR presents amino acid deletions in the carboxyterminal region of the extracellular domain which are also seen in fish and lizard FSHRs. cAMP assays with the recombinant protein transiently expressed in HEK-293 cells showed that the snake FSHR is more sensitive to human FSH (hFSH) than to human chorionic gonadotropin. Phylogenetic analysis indicated that the squamate FSHRs group separately from mammalian FSHRs. Our data are consistent with the apparently unique gonadotropin-receptor system in Squamata reptilian subgroup. Knowledge about the snake FSHR structure may help identify structural determinants for receptor function.

CALIBRATION AGE AND QUARTET DIVERGENCE DATE ESTIMATION

The date of a single divergence point--between living alligators and crocodiles--was estimated with quartet dating using calibrations of widely divergent ages. For five mitochondrial sequence datasets, there is a clear relationship between calibration age and quartet estimate--quartets based on two relatively recent calibrations support younger divergence estimates than do quartets based on two older calibrations. Some of the estimates supported by young quartets are impossibly young and exclude the first appearance of the group in the fossil record as too old. The older estimates--those based on two relatively old calibrations--may be overestimates, and those based on one old and one recent calibration support divergence estimates very close to fossil data. This suggests that quartet dating methods may be most effective when calibrations are applied from different parts of a clade's history.

Reptiles: a group of transition in the evolution of genome size and of the nucleotypic effect

A comparison between genome size and some phenotypic parameters, such as developmental length and metabolic rate, showed in reptiles a nucleotypic correlation similar to the one observed in birds and mammals. Indeed, like homeotherms, reptiles exhibit a highly significant, inverse correlation of genome size with metabolic rate but unlike amphibians, no relationship with developmental length. Several lines of evidence suggest that these nucleotypic correlations are influenced by body temperature, which also affects the guanine + cytosine nuclear percentage, and that they play an important role in the adaptation of these amniotes. However, the reptilian suborders exhibit differences in the quantitative and compositional characters of the genome that do not completely correspond to differences in the phenotypic parameters commonly involved in the nucleotypic effect. Thus, additional factors could have influenced genome size in this class. These data could be explained with the model of Hartl and Petrov, who observed an inverse correlation between genome size, non-coding portion of the genome and rate of DNA loss and hypothesized a strong role for different spectra of spontaneous insertions and deletions (indels) in the variations of genome size. It is thus reasonable to surmise that variations in the reptilian genome were initially influenced by different indels spectra typical of the diverse lineages, possibly related to different chromosome compartmentalizations. The consequent size increases or decreases would have influenced various morphological and functional cell parameters, and through these some phenotypic characteristics of the whole organism, especially the metabolic rate, very important for environmental adaptation and thus subject to natural selection. Through this "nucleotypic" bond, natural selection would also have controlled genome size variations.