DNA methylation in reptiles

Comparative analysis of genome-scale, base-resolution DNA methylation profiles across 580 animal species

Methylation of cytosines is a prototypic epigenetic modification of the DNA. It has been implicated in various regulatory mechanisms across the animal kingdom and particularly in vertebrates. We mapped DNA methylation in 580 animal species (535 vertebrates, 45 invertebrates), resulting in 2443 genome-scale DNA methylation profiles of multiple organs. Bioinformatic analysis of this large dataset quantified the association of DNA methylation with the underlying genomic DNA sequence throughout vertebrate evolution. We observed a broadly conserved link with two major transitions-once in the first vertebrates and again with the emergence of reptiles. Cross-species comparisons focusing on individual organs supported a deeply conserved association of DNA methylation with tissue type, and cross-mapping analysis of DNA methylation at gene promoters revealed evolutionary changes for orthologous genes. In summary, this study establishes a large resource of vertebrate and invertebrate DNA methylomes, it showcases the power of reference-free epigenome analysis in species for which no reference genomes are available, and it contributes an epigenetic perspective to the study of vertebrate evolution.

A second type of N7-guanine RNA cap methyltransferase in an unusual locus of a large RNA virus genome

The order Nidovirales is a diverse group of (+)RNA viruses, with a common genome organization and conserved set of replicative and editing enzymes. In particular, RNA methyltransferases play a central role in mRNA stability and immune escape. However, their presence and distribution in different Nidovirales families is not homogeneous. In Coronaviridae, the best characterized family, two distinct methytransferases perform methylation of the N7-guanine and 2'-OH of the RNA-cap to generate a cap-1 structure (m7GpppNm). The genes of both of these enzymes are located in the ORF1b genomic region. While 2'-O-MTases can be identified for most other families based on conservation of both sequence motifs and genetic loci, identification of the N7-guanine methyltransferase has proved more challenging. Recently, we identified a putative N7-MTase domain in the ORF1a region (N7-MT-1a) of certain members of the large genome Tobaniviridae family. Here, we demonstrate that this domain indeed harbors N7-specific methyltransferase activity. We present its structure as the first N7-specific Rossmann-fold (RF) MTase identified for (+)RNA viruses, making it remarkably different from that of the known Coronaviridae ORF1b N7-MTase gene. We discuss the evolutionary implications of such an appearance in this unexpected location in the genome, which introduces a split-off in the classification of Tobaniviridae.

Hierarchical dinucleotide distribution in genome along evolution and its effect on chromatin packing

It describes how hierarchical CpG distribution on the genomes of species change in the evolution and the correlation with the chromatin structure.

Dinucleotide densities and their distribution patterns vary significantly among species. Previous studies revealed that CpG is susceptible to methylation, enriched at topologically associating domain boundaries and its distribution along the genome correlates with chromatin compartmentalization. However, the multi-scale organizations of CpG in the linear genome, their role in chromatin organization, and how they change along the evolution are only partially understood. By comparing the CpG distribution at different genomic length scales, we quantify the difference between the CpG distributions of different species and evaluate how the hierarchical uneven CpG distribution appears in evolution. The clustering of species based on the CpG distribution is consistent with the phylogenetic tree. Interestingly, we found the CpG distribution and chromatin structure to be correlated in many different length scales, especially for mammals and avians, consistent with the mosaic CpG distribution in the genomes of these species.

Embryonic Temperature Programs Phenotype in Reptiles

Reptiles are critically affected by temperature throughout their lifespan, but especially so during early development. Temperature-induced changes in phenotype are a specific example of a broader phenomenon called phenotypic plasticity in which a single individual is able to develop different phenotypes when exposed to different environments. With climate change occurring at an unprecedented rate, it is important to study temperature effects on reptiles. For example, the potential impact of global warming is especially pronounced in species with temperature-dependent sex determination (TSD) because temperature has a direct effect on a key phenotypic (sex) and demographic (population sex ratios) trait. Reptiles with TSD also serve as models for studying temperature effects on the development of other traits that display continuous variation. Temperature directly influences metabolic and developmental rate of embryos and can have permanent effects on phenotype that last beyond the embryonic period. For instance, incubation temperature programs post-hatching hormone production and growth physiology, which can profoundly influence fitness. Here, we review current knowledge of temperature effects on phenotypic and developmental plasticity in reptiles. First, we examine the direct effect of temperature on biophysical processes, the concept of thermal performance curves, and the process of thermal acclimation. After discussing these reversible temperature effects, we focus the bulk of the review on developmental programming of phenotype by temperature during embryogenesis (i.e., permanent developmental effects). We focus on oviparous species because eggs are especially susceptible to changes in ambient temperature. We then discuss recent work probing the role of epigenetic mechanisms in mediating temperature effects on phenotype. Based on phenotypic effects of temperature, we return to the potential impact of global warming on reptiles. Finally, we highlight key areas for future research, including the identification of temperature sensors and assessment of genetic variation for thermosensitivity.

The untapped potential of reptile biodiversity for understanding how and why animals age

1. The field of comparative aging biology has greatly expanded in the past 20 years. Longitudinal studies of populations of reptiles with a range of maximum lifespans have accumulated and been analyzed for evidence of mortality senescence and reproductive decline. While not as well represented in studies of amniote senescence, reptiles have been the subjects of many recent demographic and mechanistic studies of the biology of aging. 2. We review recent literature on reptile demographic senescence, mechanisms of senescence, and identify unanswered questions. Given the ecophysiological and demographic diversity of reptiles, what is the expected range of reptile senescence rates? Are known mechanisms of aging in reptiles consistent with canonical hallmarks of aging in model systems? What are the knowledge gaps in our understanding of reptile aging? 3. We find ample evidence of increasing mortality with advancing age in many reptiles. Testudines stand out as slower aging than other orders, but data on crocodilians and tuatara are sparse. Sex-specific analyses are generally not available. Studies of female reproduction suggest that reptiles are less likely to have reproductive decline with advancing age than mammals. 4. Reptiles share many physiological and molecular pathways of aging with mammals, birds, and laboratory model organisms. Adaptations related to stress physiology coupled with reptilian ectothermy suggest novel comparisons and contrasts that can be made with canonical aging phenotypes in mammals. These include stem cell and regeneration biology, homeostatic mechanisms, IIS/TOR signaling, and DNA repair. 5. To overcome challenges to the study of reptile aging, we recommend extending and expanding long-term monitoring of reptile populations, developing reptile cell lines to aid cellular biology, conducting more comparative studies of reptile morphology and physiology sampled along relevant life-history axes, and sequencing more reptile genomes for comparative genomics. Given the diversity of reptile life histories and adaptations, achieving these directives will likely greatly benefit all aging biology.

Differential hypomethylation of the repetitive Tol2/Alu-rich sequences in the genome of Bodianus species (Labriformes, Labridae)

Representatives of the order Labriformes show karyotypes of extreme conservatism together with others with high chromosomal diversification. However, the cytological characterization of epigenetic modifications remains unknown for the majority of the species. In the family Labridae, the most abundant fishes on tropical reefs, the genomes of the genus Bodianus Bloch, 1790 have been characterized by the occurrence of a peculiar chromosomal region, here denominated BOD. This region is exceptionally decondensed, heterochromatic, argentophilic, GC-neutral and, in contrast to classical secondary constrictions, shows no signals of hybridization with 18S rDNA probes. In order to characterize the BOD region, the methylation pattern, the distribution of Alu and Tol2 retrotransposons and of 18S and 5S rDNA sites, respectively, were analyzed by Fluorescence In Situ Hybridization (FISH) on metaphase chromosomes of two Bodianus species, B. insularis Gomon & Lubbock, 1980 and B. pulchellus (Poey, 1860). Immunolocalization of the 5-methylcytosine revealed hypermethylated chromosomal regions, dispersed along the entire length of the chromosomes of both species, while the BOD regions exhibited a hypomethylated pattern. Hypomethylation of the BOD region is associated with the precise co-location of Tol2 and Alu elements, suggesting their active participation in the regulatory epigenetic process. This evidence underscores a probable differential methylation action during the cell cycle, as well as the role of Tol2/Alu elements in functional processes of fish genomes.

Measuring Animal Age with DNA Methylation: From Humans to Wild Animals

DNA methylation (DNAm) is a key mechanism for regulating gene expression in animals and levels are known to change with age. Recent studies have used DNAm changes as a biomarker to estimate chronological age in humans and these techniques are now also being applied to domestic and wild animals. Animal age is widely used to track ongoing changes in ecosystems, however chronological age information is often unavailable for wild animals. An ability to estimate age would lead to improved monitoring of (i) population trends and status and (ii) demographic properties such as age structure and reproductive performance. Recent studies have revealed new examples of DNAm age association in several new species increasing the potential for developing DNAm age biomarkers for a broad range of wild animals. Emerging technologies for measuring DNAm will also enhance our ability to study age-related DNAm changes and to develop new molecular age biomarkers.

MeDIP-seq and nCpG analyses illuminate sexually dimorphic methylation of gonadal development genes with high historic methylation in turtle hatchlings with temperature-dependent sex determination

Background

DNA methylation alters gene expression but not DNA sequence and mediates some cases of phenotypic plasticity. Temperature-dependent sex determination (TSD) epitomizes phenotypic plasticity where environmental temperature drives embryonic sexual fate, as occurs commonly in turtles. Importantly, the temperature-specific transcription of two genes underlying gonadal differentiation is known to be induced by differential methylation in TSD fish, turtle and alligator. Yet, how extensive is the link between DNA methylation and TSD remains unclear. Here we test for broad differences in genome-wide DNA methylation between male and female hatchling gonads of the TSD painted turtle Chrysemys picta using methyl DNA immunoprecipitation sequencing, to identify differentially methylated candidates for future study. We also examine the genome-wide nCpG distribution (which affects DNA methylation) in painted turtles and test for historic methylation in genes regulating vertebrate gonadogenesis.

Results

Turtle global methylation was consistent with other vertebrates (57% of the genome, 78% of all CpG dinucleotides). Numerous genes predicted to regulate turtle gonadogenesis exhibited sex-specific methylation and were proximal to methylated repeats. nCpG distribution predicted actual turtle DNA methylation and was bimodal in gene promoters (as other vertebrates) and introns (unlike other vertebrates). Differentially methylated genes, including regulators of sexual development, had lower nCpG content indicative of higher historic methylation.

Conclusions

Ours is the first evidence suggesting that sexually dimorphic DNA methylation is pervasive in turtle gonads (perhaps mediated by repeat methylation) and that it targets numerous regulators of gonadal development, consistent with the hypothesis that it may regulate thermosensitive transcription in TSD vertebrates. However, further research during embryogenesis will help test this hypothesis and the alternative that instead, most differential methylation observed in hatchlings is the by-product of sexual differentiation and not its cause.

Electronic supplementary material

The online version of this article (doi:10.1186/s13072-017-0136-2) contains supplementary material, which is available to authorized users.

ALS and FTD: an epigenetic perspective

Amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD) are two fatal neurodegenerative diseases seen in comorbidity in up to 50 % of cases. Despite tremendous efforts over the last two decades, no biomarkers or effective therapeutics have been identified to prevent, decelerate, or stop neuronal death in patients. While the identification of multiple mutations in more than two dozen genes elucidated the involvement of several mechanisms in the pathogenesis of both diseases, identifying the hexanucleotide repeat expansion in C9orf72, the most common genetic abnormality in ALS and FTD, opened the door to the discovery of several novel pathogenic biological routes, including chromatin remodeling and transcriptome alteration. Epigenetic processes regulate DNA replication and repair, RNA transcription, and chromatin conformation, which in turn further dictate transcriptional regulation and protein translation. Transcriptional and post-transcriptional epigenetic regulation is mediated by enzymes and chromatin-modifying complexes that control DNA methylation, histone modifications, and RNA editing. While the alteration of DNA methylation and histone modification has recently been reported in ALS and FTD, the assessment of epigenetic involvement in both diseases is still at an early stage, and the involvement of multiple epigenetic players still needs to be evaluated. As the epigenome serves as a way to alter genetic information not only during aging, but also following environmental signals, epigenetic mechanisms might play a central role in initiating ALS and FTD, especially for sporadic cases. Here, we provide a review of what is currently known about altered epigenetic processes in both ALS and FTD and discuss potential therapeutic strategies targeting epigenetic mechanisms. As approximately 85 % of ALS and FTD cases are still genetically unexplained, epigenetic therapeutics explored for other diseases might represent a profitable direction for the field.

Nonmodal scutes patterns in the Loggerhead Sea Turtle (Caretta caretta): a possible epigenetic effect?

Eleven specimens of the threatened Loggerhead Sea Turtle (Caretta caretta (L., 1758)) were caught accidentally by fishermen in different parts of Sicily (Mediterranean Sea). Five of them showed an atypical number of carapacial and plastron scutes, making the immediate identification of the specimens as C. caretta difficult. Both genetic and epigenetic analysis were carried out on these specimens. Sequencing of a 649 bp sequence of the mitochondrial cytochrome c oxidase I (COI) gene allowed us to classify all the individuals as C. caretta. Epigenetic analysis, performed by evaluating the total level of DNA cytosine methylation, showed a reduced and significant (F = 72.65, p < 0.01) global level of methylated cytosines in the turtles with nonmodal scutes compared with the normal turtles. Our results suggest that the variability in the number of scutes could be dependent on the environmental conditions during embryonic incubation, via an epigenetic mechanism. This finding could have implications in our understanding of the pathways of morphological evolution and diversification in the chelonians.

DNA methylation and temperature stress in an Antarctic polychaete, Spiophanes tcherniai

Epigenetic modifications of DNA and histones are a primary mechanism by which gene expression activities may be modified in response to environmental stimuli. Here we characterize patterns of methyl-cytosine composition in the marine polychaete Spiophanes tcherniai from McMurdo Sound, Antarctica. We cultured adult worms at two temperatures, -1.5°C (ambient control) and +4°C (warm treatment), for 4 weeks. We observed a rapid capacity for S. tcherniai organismal respiration rates and underlying catalytic rates of citrate synthase at +4°C to return to control levels in less than 4 weeks. We profiled changes in the methylation states of CpG sites in these treatments using an NGS strategy to computationally reconstruct and quantify methylation status across the genome. In our analysis we recovered 120,000 CpG sites in assembled contigs from both treatments. Of those, we were able to align 28,000 CpG sites in common between the two sample groups. In comparing these aligned sites between treatments, only 3000 (11%) evidenced a change in methylation state, but over 85% of changes involved a gain of a 5-methyl group on a CpG site (net increase in methyation). The ability to score CpG sites as partially methylated among gDNA copies in a sample opens up a new avenue for assessing DNA methylation responses to changing environments. By quantitatively distinguishing a "mixed" population of copies of one CpG site, we can begin to identify dynamic, non-binary, continuous-response reactions in DNA methylation intensity or density that previously may have been overlooked as noise.

Patterns of DNA methylation in animals: an ecotoxicological perspective

DNA methylation refers to the addition of a methyl group to nucleotides within DNA. As with other epigenetic endpoints, patterns of DNA methylation are susceptible to alterations due to exposure to environmental stressors, including contaminants. These alterations can persist in the absence of the initial stressor as cells divide, and can even be inherited between generations if they occur in the germ line. Although our knowledge concerning patterns of DNA methylation in animals is increasing, there remains a gap in the literature when it comes to species outside of those typically used for biomedical research. Here, I review the literature relating to DNA methylation in an array of taxa (mammals, fish, birds, amphibians, reptiles, and invertebrates) and discuss these data from an ecotoxicological perspective. The pattern and extent of DNA methylation is well conserved across species of vertebrates; methylation appears mainly on cytosine residues within a CpG context, and much of the genome is methylated, with the notable exception of cytosines within CpG islands in the promoters of genes. Highly methylated genes in vertebrates tend to be transcriptionally repressed. However, large differences occur between classes of vertebrates in terms of the timing and nature of reprogramming and genomic imprinting: epigenetic processes that establish patterns of DNA methylation in the early embryo and which are sensitive to environmental stress. In invertebrates, patterns of DNA methylation are extremely variable and differ significantly from the condition observed in vertebrates. Some invertebrate genomes exhibit no DNA methylation while others are methylated to a level that is comparable to vertebrates. Additionally, DNA methylation may have different functions in invertebrates, e.g., alternative splicing. This variability in basic patterns of DNA methylation among species during sensitive periods of development suggests that responses to epigenetically active environmental contaminants may be similarly variable. For example, the timing of exposure to a contaminant may be a critical factor when considered in the light of variable reprogramming schedules among species. With this in mind, I review data relating to the effects of contaminants on DNA methylation in animals, focusing on non-model organisms and on exposures in natural environments, when possible. An ecotoxicological perspective on patterns of DNA methylation in animals may improve our understanding of the range and diversity of epigenetic phenomena in the natural world.

Application of the LUminometric Methylation Assay to ecological species: tissue quality requirements and a survey of DNA methylation levels in animals

The LUminometric Methylation Assay (LUMA) measures global DNA methylation. LUMA depends on digestion of DNA with methyl-sensitive and methyl-insensitive restriction enzymes, followed by pyrosequencing. Until recently, LUMA has been principally used for biomedical research. Here, we use chickens as a model to investigate sample quality issues relating to LUMA and then apply the method to ecological species. First, we assessed the effect of tissue storage conditions on DNA methylation values. This is an important consideration for ecological species because samples are not always ideally preserved and LUMA is sensitive to poor DNA quality. We found that good quality LUMA data could be obtained from chicken liver and brain tissues stored at 21 °C for at least 2 and 12 h, respectively. Longer storage times introduced nonspecific peaks to pyrograms which were associated with reduced DNA methylation. Repeatedly, freezing and thawing the tissues did not affect LUMA data. Second, we measured DNA methylation in 12 species representing five animal classes: amphibians (African and Western clawed frog), reptiles (green anole lizard), fish (yellow perch, goldfish, lake trout), mammals (American mink, polar bear, short-beaked common dolphin, Atlantic white-sided dolphin) and birds (chicken, Japanese quail). We saw a pattern of high DNA methylation in fish (84-87%), and intermediate levels in mammals (68-72%) and birds (52-71%). This pattern corresponds well with previous measures of DNA methylation generated by HPLC. Our data represent the first CpG methylation values to be reported in several species and provide a basis for studying patterns of epigenetic inheritance in an ecological context.

DNA methylation, epigenetics, and evolution in vertebrates: facts and challenges

DNA methylation is a key epigenetic modification in the vertebrate genomes known to be involved in biological processes such as regulation of gene expression, DNA structure and control of transposable elements. Despite increasing knowledge about DNA methylation, we still lack a complete understanding of its specific functions and correlation with environment and gene expression in diverse organisms. To understand how global DNA methylation levels changed under environmental influence during vertebrate evolution, we analyzed its distribution pattern along the whole genome in mammals, reptiles and fishes showing that it is correlated with temperature, independently on phylogenetic inheritance. Other studies in mammals and plants have evidenced that environmental stimuli can promote epigenetic changes that, in turn, might generate localized changes in DNA sequence resulting in phenotypic effects. All these observations suggest that environment can affect the epigenome of vertebrates by generating hugely different methylation patterns that could, possibly, reflect in phenotypic differences. We are at the first steps towards the understanding of mechanisms that underlie the role of environment in molding the entire genome over evolutionary times. The next challenge will be to map similarities and differences of DNA methylation in vertebrates and to associate them with environmental adaptation and evolution.

The correlation of genome size and DNA methylation rate in metazoans

Total DNA methylation rates are well known to vary widely between different metazoans. The phylogenetic distribution of this variation, however, has not been investigated systematically. We combine here publicly available data on methylcytosine content with the analysis of nucleotide compositions of genomes and transcriptomes of 78 metazoan species to trace the evolution of abundance and distribution of DNA methylation. The depletion of CpG and the associated enrichment of TpG and CpA dinucleotides are used to infer the intensity and localization of germline CpG methylation and to estimate its evolutionary dynamics. We observe a positive correlation of the relative methylation of CpG motifs with genome size. We tested this trend successfully by measuring total DNA methylation with LC/MS in orthopteran insects with very different genome sizes: house crickets, migratory locusts and meadow grasshoppers. We hypothesize that the observed correlation between methylation rate and genome size is due to a dependence of both variables from long-term effective population size and is driven by the accumulation of repetitive sequences that are typically methylated during periods of small population sizes. This process may result in generally methylated, large genomes such as those of jawed vertebrates. In this case, the emergence of a novel demethylation pathway and of novel reader proteins for methylcytosine may have enabled the usage of cytosine methylation for promoter-based gene regulation. On the other hand, persistently large populations may lead to a compression of the genome and to the loss of the DNA methylation machinery, as observed, e.g., in nematodes.

Single amino acid and trinucleotide repeats: function and evolution

The most well known effect of single amino acid repeat expansion, beyond a certain threshold, is the development of a specific disease, depending on the protein in which the expansion has occurred. For example, the expansion of the glutamine repeat in huntingtin leads to the debilitating neurodegenerative disease, Huntington's disease. Similarly, there are a range of other disorders caused by trinucleotide repeat expansions encoding polyglutamine or polyalanine tracts. The age of onset of the polyglutamine-induced neurodegenerative diseases is usually negatively correlated with the length of expanded CAG/glutamine repeat. However, recent studies have given evidence that single amino acid repeats may also play critical roles in normal protein function and that changes in the length of single amino acid repeats is likely to play a beneficial role in evolution. This chapter will look at the prevalence, function and possible role single amino acid repeats have in evolution and other biological processes.

Ecogeographical associations between climate and human body composition: analyses based on anthropometry and skinfolds

In the 19th century, two "ecogeographical rules" were proposed hypothesizing associations of climate with mammalian body size and proportions. Data on human body weight and relative leg length support these rules; however, it is unknown whether such associations are attributable to lean tissue (the heat-producing component) or fat (energy stores). Data on weight, height, and two skinfold thickness were obtained from the literature for 137 nonindustrialized populations, providing 145 male and 115 female individual samples. A variety of indices of adiposity and lean mass were analyzed. Preliminary analyses indicated secular increases in skinfolds in men but not women, and associations of age and height with lean mass in both sexes. Decreasing annual temperature was associated with increasing body mass index (BMI), and increasing triceps but not subscapular skinfold. After adjusting for skinfolds, decreasing temperature remained associated with increasing BMI. These results indicate that colder environments favor both greater peripheral energy stores, and greater lean mass. Contrasting results for triceps and subscapular skinfolds might be due to adaptive strategies either constraining central adiposity in cold environments to reduce cardiovascular risk, or favoring central adiposity in warmer environments to maintain energetic support of the immune system. Polynesian populations were analyzed separately and contradicted all of the climate trends, indicating support for the hypothesis that they are cold-adapted despite occupying a tropical region. It is unclear whether such associations emerge through natural selection or through trans-generational and life-course plasticity. These findings nevertheless aid understanding of the wide variability in human physique and adiposity.

Evolutionary genomics implies a specific function of Ant4 in mammalian and anole lizard male germ cells

Most vertebrates have three paralogous genes with identical intron-exon structures and a high degree of sequence identity that encode mitochondrial adenine nucleotide translocase (Ant) proteins, Ant1 (Slc25a4), Ant2 (Slc25a5) and Ant3 (Slc25a6). Recently, we and others identified a fourth mammalian Ant paralog, Ant4 (Slc25a31), with a distinct intron-exon structure and a lower degree of sequence identity. Ant4 was expressed selectively in testis and sperm in adult mammals and was indeed essential for mouse spermatogenesis, but it was absent in birds, fish and frogs. Since Ant2 is X-linked in mammalian genomes, we hypothesized that the autosomal Ant4 gene may compensate for the loss of Ant2 gene expression during male meiosis in mammals. Here we report that the Ant4 ortholog is conserved in green anole lizard (Anolis carolinensis) and demonstrate that it is expressed in the anole testis. Further, a degenerate DNA fragment of putative Ant4 gene was identified in syntenic regions of avian genomes, indicating that Ant4 was present in the common amniote ancestor. Phylogenetic analyses suggest an even more ancient origin of the Ant4 gene. Although anole lizards are presumed male (XY) heterogametic, like mammals, copy numbers of the Ant2 as well as its neighboring gene were similar between male and female anole genomes, indicating that the anole Ant2 gene is either autosomal or located in the pseudoautosomal region of the sex chromosomes, in contrast to the case to mammals. These results imply the conservation of Ant4 is not likely simply driven by the sex chromosomal localization of the Ant2 gene and its subsequent inactivation during male meiosis. Taken together with the fact that Ant4 protein has a uniquely conserved structure when compared to other somatic Ant1, 2 and 3, there may be a specific advantage for mammals and lizards to express Ant4 in their male germ cells.

Distribution of DNA methylation, CpGs, and CpG islands in human isochores

DNA methylation is a major epigenetic modification of the genome that affects basic biological functions, such as gene expression and cell development. We used the human genome sequences and the DNA methylation data that are available in order to establish a map of the levels of GC and methylation in isochores. We also looked for the correlations that hold between GC levels and the distribution of the (1) dinucleotide CpG, (2) ratio 5mC/CpG, and (3) CpG islands. Our results show that methylation levels, CpG frequencies, and the density of CpG islands are positively correlated with the GC level of isochores. In contrast, the correlation between the 5mC/CpG ratio and GC is a negative one because the increase in methylation lags behind that of CpG, to reach a plateau in the GC-richest, gene-richest isochore families H2 and H3. In conclusion, there are more CpG targets that remain unmethylated in the GC-richest, gene-richest isochores in comparison with the other isochores. This conclusion supports the idea that the widespread methylation under consideration here has a general inhibitory effect on gene expression.

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.

Turtle isochore structure is intermediate between amphibians and other amniotes

Vertebrate genomes are comprised of isochores that are relatively long (>100 kb) regions with a relatively homogenous (either GC-rich or AT-rich) base composition and with rather sharp boundaries with neighboring isochores. Mammals and living archosaurs (birds and crocodilians) have heterogeneous genomes that include very GC-rich isochores. In sharp contrast, the genomes of amphibians and fishes are more homogeneous and they have a lower overall GC content. Because DNA with higher GC content is more thermostable, the elevated GC content of mammalian and archosaurian DNA has been hypothesized to be an adaptation to higher body temperatures. This hypothesis can be tested by examining structure of isochores across the reptilian clade, which includes the archosaurs, testudines (turtles), and lepidosaurs (lizards and snakes), because reptiles exhibit diverse body sizes, metabolic rates, and patterns of thermoregulation. This study focuses on a comparative analysis of a new set of expressed genes of the red-eared slider turtle and orthologs of the turtle genes in mammalian (human, mouse, dog, and opossum), archosaurian (chicken and alligator), and amphibian (western clawed frog) genomes. EST (expressed sequence tag) data from a turtle cDNA library enriched for genes that have specialized functions (developmental genes) revealed using the GC content of the third-codon-position to examine isochore structure requires careful consideration of the types of genes examined. The more highly expressed genes (e.g., housekeeping genes) are more likely to be GC-rich than are genes with specialized functions. However, the set of highly expressed turtle genes demonstrated that the turtle genome has a GC content that is intermediate between the GC-poor amphibians and the GC-rich mammals and archosaurs. There was a strong correlation between the GC content of all turtle genes and the GC content of other vertebrate genes, with the slope of the line describing this relationship also indicating that the isochore structure of turtles is intermediate between that of amphibians and other amniotes. These data are consistent with some thermal hypotheses of isochore evolution, but we believe that the credible set of models for isochore evolution still includes a variety of models. These data expand the amount of genomic data available from reptiles upon which future studies of reptilian genomics can build.

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.

Patterns of Vertebrate Isochore Evolution Revealed by Comparison of Expressed Mammalian, Avian, and Crocodilian Genes

Vertebrate genomes are mosaics of isochores, defined as long (>100 kb) regions with relatively homogeneous within-region base composition. Birds and mammals have more GC-rich isochores than amphibians and fish, and the GC-rich isochores of birds and mammals have been suggested to be an adaptation to homeothermy. If this hypothesis is correct, all poikilothermic (cold-blooded) vertebrates, including the nonavian reptiles, are expected to lack a GC-rich isochore structure. Previous studies using various methods to examine isochore structure in crocodilians, turtles, and squamates have led to different conclusions. We collected more than 6000 expressed sequence tags (ESTs) from the American alligator to overcome sample size limitations suggested to be the fundamental problem in the previous reptilian studies. The alligator ESTs were assembled and aligned with their human, mouse, chicken, and western clawed frog orthologs, resulting in 366 alignments. Analyses of third-codon-position GC content provided conclusive evidence that the poikilothermic alligator has GC-rich isochores, like homeothermic birds and mammals. We placed these results in a theoretical framework able to unify available models of isochore evolution. The data collected for this study allowed us to reject the models that explain the evolution of GC content using changes in body temperature associated with the transition from poikilothermy to homeothermy. Falsification of these models places fundamental constraints upon the plausible pathways for the evolution of isochores.

The neoselectionist theory of genome evolution

The vertebrate genome is a mosaic of GC-poor and GC-rich isochores, megabase-sized DNA regions of fairly homogeneous base composition that differ in relative amount, gene density, gene expression, replication timing, and recombination frequency. At the emergence of warm-blooded vertebrates, the gene-rich, moderately GC-rich isochores of the cold-blooded ancestors underwent a GC increase. This increase was similar in mammals and birds and was maintained during the evolution of mammalian and avian orders. Neither the GC increase nor its conservation can be accounted for by the random fixation of neutral or nearly neutral single-nucleotide changes (i.e., the vast majority of nucleotide substitutions) or by a biased gene conversion process occurring at random genome locations. Both phenomena can be explained, however, by the neoselectionist theory of genome evolution that is presented here. This theory fully accepts Ohta's nearly neutral view of point mutations but proposes in addition (i) that the AT-biased mutational input present in vertebrates pushes some DNA regions below a certain GC threshold; (ii) that these lower GC levels cause regional changes in chromatin structure that lead to deleterious effects on replication and transcription; and (iii) that the carriers of these changes undergo negative (purifying) selection, the final result being a compositional conservation of the original isochore pattern in the surviving population. Negative selection may also largely explain the GC increase accompanying the emergence of warm-blooded vertebrates. In conclusion, the neoselectionist theory not only provides a solution to the neutralist/selectionist debate but also introduces an epigenomic component in genome evolution.

DNA methylation and body temperature in fishes

Previous investigations from our laboratory [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.] led to the discovery of two different methylation levels in the genomes of vertebrates, a higher one exhibited by fishes and amphibians and a lower one shown by mammals and birds. It was also noted that data from the literature indicated a higher CpG level in fishes and amphibians compared to mammals and birds. Such observations led to suggesting the existence of two equilibria and to speculate that the transitions between the two equilibria in DNA methylation and CpG levels were due to a higher deamination rate in warm-blooded vertebrates related to their higher body temperature. Here we used Reverse-Phase High-Performance Liquid Chromatography (RP-HPLC) analysis to study methylation levels in a number of fish genomes living at different temperatures. We found that polar fishes exhibit DNA methylation levels that are higher than those of tropical and temperate fishes, the latter being in turn higher than the methylation levels of warm-blooded vertebrates, as expected from previous work. A closer analysis of the data revealed that, among Antarctic fishes, the Channichthyidae (the icefishes, deprived of haemoglobin) had the highest methylation level, and that, among temperate and tropical fishes the latter showed the lowest methylation level. These results confirm the existence of an inverse relationship between DNA methylation and body temperature, when the latter is maintained over evolutionary times.