cDNA-based gene mapping and GC3 profiling in the soft-shelled turtle suggest a chromosomal size-dependent GC bias shared by sauropsids

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.

The Cytogenetic Map of the Nile Crocodile (Crocodylus niloticus, Crocodylidae, Reptilia) with Fluorescence In Situ Localization of Major Repetitive DNAs

Tandemly arranged and dispersed repetitive DNA sequences are important structural and functional elements that make up a significant portion of vertebrate genomes. Using high throughput, low coverage whole genome sequencing followed by bioinformatics analysis, we have identified seven major tandem repetitive DNAs and two fragments of LTR retrotransposons in the genome of the Nile crocodile (Crocodylus niloticus, 2n = 32). The repeats showed great variability in structure, genomic organization, and chromosomal distribution as revealed by fluorescence in situ hybridization (FISH). We found that centromeric and pericentromeric heterochromatin of C. niloticus is composed of previously described in Crocodylus siamensis CSI-HindIII and CSI-DraI repetitive sequence families, a satellite revealed in Crocodylus porosus, and additionally contains at least three previously unannotated tandem repeats. Both LTR sequences identified here belong to the ERV1 family of endogenous retroviruses. Each pericentromeric region was characterized by a diverse set of repeats, with the exception of chromosome pair 4, in which we found only one type of satellite. Only a few repeats showed non-centromeric signals in addition to their centromeric localization. Mapping of 18S-28S ribosomal RNA genes and telomeric sequences (TTAGGG)n did not demonstrate any co-localization of these sequences with revealed centromeric and pericentromeric heterochromatic blocks.

A chromosome-level genome assembly and annotation of the desert horned lizard, Phrynosoma platyrhinos, provides insight into chromosomal rearrangements among reptiles

BACKGROUND

The increasing number of chromosome-level genome assemblies has advanced our knowledge and understanding of macroevolutionary processes. Here, we introduce the genome of the desert horned lizard, Phrynosoma platyrhinos, an iguanid lizard occupying extreme desert conditions of the American southwest. We conduct analysis of the chromosomal structure and composition of this species and compare these features across genomes of 12 other reptiles (5 species of lizards, 3 snakes, 3 turtles, and 1 bird).

FINDINGS

The desert horned lizard genome was sequenced using Illumina paired-end reads and assembled and scaffolded using Dovetail Genomics Hi-C and Chicago long-range contact data. The resulting genome assembly has a total length of 1,901.85 Mb, scaffold N50 length of 273.213 Mb, and includes 5,294 scaffolds. The chromosome-level assembly is composed of 6 macrochromosomes and 11 microchromosomes. A total of 20,764 genes were annotated in the assembly. GC content and gene density are higher for microchromosomes than macrochromosomes, while repeat element distributions show the opposite trend. Pathway analyses provide preliminary evidence that microchromosome and macrochromosome gene content are functionally distinct. Synteny analysis indicates that large microchromosome blocks are conserved among closely related species, whereas macrochromosomes show evidence of frequent fusion and fission events among reptiles, even between closely related species.

CONCLUSIONS

Our results demonstrate dynamic karyotypic evolution across Reptilia, with frequent inferred splits, fusions, and rearrangements that have resulted in shuffling of chromosomal blocks between macrochromosomes and microchromosomes. Our analyses also provide new evidence for distinct gene content and chromosomal structure between microchromosomes and macrochromosomes within reptiles.

Revisiting the Karyotypes of Alligators and Caimans (Crocodylia, Alligatoridae) after a Half-Century Delay: Bridging the Gap in the Chromosomal Evolution of Reptiles

Although crocodilians have attracted enormous attention in other research fields, from the cytogenetic point of view, this group remains understudied. Here, we analyzed the karyotypes of eight species formally described from the Alligatoridae family using differential staining, fluorescence in situ hybridization with rDNA and repetitive motifs as a probe, whole chromosome painting (WCP), and comparative genome hybridization. All Caimaninae species have a diploid chromosome number (2n) 42 and karyotypes dominated by acrocentric chromosomes, in contrast to both species of Alligatorinae, which have 2n = 32 and karyotypes that are predominantly metacentric, suggesting fusion/fission rearrangements. Our WCP results supported this scenario by revealing the homeology of the largest metacentric pair present in both Alligator spp. with two smaller pairs of acrocentrics in Caimaninae species. The clusters of 18S rDNA were found on one chromosome pair in all species, except for Paleosuchus spp., which possessed three chromosome pairs bearing these sites. Similarly, comparative genomic hybridization demonstrated an advanced stage of sequence divergence among the caiman genomes, with Paleosuchus standing out as the most divergent. Thus, although Alligatoridae exhibited rather low species diversity and some level of karyotype stasis, their genomic content indicates that they are not as conserved as previously thought. These new data deepen the discussion of cytotaxonomy in this family.

Multifaceted Hi-C benchmarking: what makes a difference in chromosome-scale genome scaffolding?

Background

Hi-C is derived from chromosome conformation capture (3C) and targets chromatin contacts on a genomic scale. This method has also been used frequently in scaffolding nucleotide sequences obtained by de novo genome sequencing and assembly, in which the number of resultant sequences rarely converges to the chromosome number. Despite its prevalent use, the sample preparation methods for Hi-C have not been intensively discussed, especially from the standpoint of genome scaffolding.

Results

To gain insight into the best practice of Hi-C scaffolding, we performed a multifaceted methodological comparison using vertebrate samples and optimized various factors during sample preparation, sequencing, and computation. As a result, we identified several key factors that helped improve Hi-C scaffolding, including the choice and preparation of tissues, library preparation conditions, the choice of restriction enzyme(s), and the choice of scaffolding program and its usage.

Conclusions

This study provides the first comparison of multiple sample preparation kits/protocols and computational programs for Hi-C scaffolding by an academic third party. We introduce a customized protocol designated “inexpensive and controllable Hi-C (iconHi-C) protocol,” which incorporates the optimal conditions identified in this study, and demonstrate this technique on chromosome-scale genome sequences of the Chinese softshell turtle Pelodiscus sinensis.

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.

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.

Madagascar ground gecko genome analysis characterizes asymmetric fates of duplicated genes

Background

Conventionally, comparison among amniotes – birds, mammals, and reptiles – has often been approached through analyses of mammals and, for comparison, birds. However, birds are morphologically and physiologically derived and, moreover, some parts of their genomes are recognized as difficult to sequence and/or assemble and are thus missing in genome assemblies. Therefore, sequencing the genomes of reptiles would aid comparative studies on amniotes by providing more comprehensive coverage to help understand the molecular mechanisms underpinning evolutionary changes.

Results

Herein, we present the whole genome sequences of the Madagascar ground gecko (Paroedura picta), a promising study system especially in developmental biology, and used it to identify changes in gene repertoire across amniotes. The genome-wide analysis of the Madagascar ground gecko allowed us to reconstruct a comprehensive set of gene phylogenies comprising 13,043 ortholog groups from diverse amniotes. Our study revealed 469 genes retained by some reptiles but absent from available genome-wide sequence data of both mammals and birds. Importantly, these genes, herein collectively designated as ‘elusive’ genes, exhibited high nucleotide substitution rates and uneven intra-genomic distribution. Furthermore, the genomic regions flanking these elusive genes exhibited distinct characteristics that tended to be associated with increased gene density, repeat element density, and GC content.

Conclusion

This highly continuous and nearly complete genome assembly of the Madagascar ground gecko will facilitate the use of this species as an experimental animal in diverse fields of biology. Gene repertoire comparisons across amniotes further demonstrated that the fate of a duplicated gene can be affected by the intrinsic properties of its genomic location, which can persist for hundreds of millions of years.

Electronic supplementary material

The online version of this article (10.1186/s12915-018-0509-4) contains supplementary material, which is available to authorized users.

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.

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.

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.

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.

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.

The bat genome: GC-biased small chromosomes associated with reduction in genome size

Bats are distinct from other mammals in their small genome size as well as their high metabolic rate, possibly related to flight ability. Although the genome sequence has been published in two species, the data lack cytogenetic information. In this study, the size and GC content of each chromosome are measured from the flow karyotype of the mouse-eared bat, Myotis myotis (MMY). The smaller chromosomes are GC-rich compared to the larger chromosomes, and the relative proportions of homologous segments between MMY and human differ among the MMY chromosomes. The MMY genome size calculated from the sum of the chromosome sizes is 2.25 Gb, and the total GC content is 42.3%, compared to human and dog with 41.0 and 41.2%, respectively. The GC-rich small MMY genome is characterised by GC-biased smaller chromosomes resulting from preferential loss of AT-rich sequences. Although the association between GC-rich small chromosomes and small genome size has been reported only in birds so far, we show in this paper, for the first time, that the same phenomenon is observed in at least one group of mammals, implying that this may be a mechanism common to genome evolution in general.

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.

Reptiles: a new model for brain evo-devo research

Vertebrate brains exhibit vast amounts of anatomical diversity. In particular, the elaborate and complex nervous system of amniotes is correlated with the size of their behavioral repertoire. However, the evolutionary mechanisms underlying species-specific brain morphogenesis remain elusive. In this review we introduce reptiles as a new model organism for understanding brain evolution. These animal groups inherited ancestral traits of brain architectures. We will describe several unique aspects of the reptilian nervous system with a special focus on the telencephalon, and discuss the genetic mechanisms underlying reptile-specific brain morphology. The establishment of experimental evo-devo approaches to studying reptiles will help to shed light on the origin of the amniote brains.

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.

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

BACKGROUND

Extant sauropsids (reptiles and birds) are divided into two major lineages, the lineage of Testudines (turtles) and Archosauria (crocodilians and birds) and the lineage of Lepidosauria (tuatara, lizards, worm lizards and snakes). Karyotypes of these sauropsidan groups generally consist of macrochromosomes and microchromosomes. In chicken, microchromosomes exhibit a higher GC-content than macrochromosomes. To examine the pattern of intra-genomic GC heterogeneity in lepidosaurian genomes, we constructed a cytogenetic map of the Japanese four-striped rat snake (Elaphe quadrivirgata) with 183 cDNA clones by fluorescence in situ hybridization, and examined the correlation between the GC-content of exonic third codon positions (GC3) of the genes and the size of chromosomes on which the genes were localized.

RESULTS

Although GC3 distribution of snake genes was relatively homogeneous compared with those of the other amniotes, microchromosomal genes showed significantly higher GC3 than macrochromosomal genes as in chicken. Our snake cytogenetic map also identified several conserved segments between the snake macrochromosomes and the chicken microchromosomes. Cross-species comparisons revealed that GC3 of most snake orthologs in such macrochromosomal segments were GC-poor (GC3 < 50%) whereas those of chicken orthologs in microchromosomes were relatively GC-rich (GC3 ≥ 50%).

CONCLUSION

Our results suggest that the chromosome size-dependent GC heterogeneity had already occurred before the lepidosaur-archosaur split, 275 million years ago. This character was probably present in the common ancestor of lepidosaurs and but lost in the lineage leading to Anolis during the diversification of lepidosaurs. We also identified several genes whose GC-content might have been influenced by the size of the chromosomes on which they were harbored over the course of sauropsid evolution.

Different patterns of Robertsonian fusion pairing in Bovidae and the house mouse: the relationship between chromosome size and nuclear territories

Using a dataset of karyotypic changes reported for bovids and the house mouse (Mus musculus domesticus) together with information from the cattle (Bos taurus) and mouse genomes, we examined two principal variables that have been proposed to predict chromosomal positioning in the nucleus, chromosome size and GC content. These were expected to influence the distribution of Robertsonian (Rb) fusions, the predominant mode of chromosomal change in both taxa. We found the largest chromosomes to be most frequently involved in fusions in bovids, and confirm earlier reports that chromosomes of intermediate size were the most frequent fusers in mice. We then tested whether chromosomal positioning can explain Rb fusion frequencies. We classified chromosomes into groups by size and considered the frequency of interactions between specific groups. Among the interactions, mouse chromosomes showed a slight tendency to fuse with neighbouring chromosomes, in line with expectations of chromosomal positioning, but also resembling predictions from meiotic spindle-induced bias. Bovids, on the other hand, showed no trend in interactions, with small chromosomes being the least frequent partner for all size classes. We discuss the results in terms of nuclear organization at various cell cycle stages and the proposed mechanisms of Rb fusion formation, and note that the difference can be explained by (i) considering bovid species generally to be characterized by a greater intermingling of chromosomal size classes than the house mouse, or (ii) by the vastly different timescales underpinning their evolutionary histories.

Reassessment of genome size in turtle and crocodile based on chromosome measurement by flow karyotyping: close similarity to chicken

The genome size in turtles and crocodiles is thought to be much larger than the 1.2 Gb of the chicken (Gallus gallus domesticus, GGA), according to the animal genome size database. However, GGA macrochromosomes show extensive homology in the karyotypes of the red eared slider (Trachemys scripta elegans, TSC) and the Nile crocodile (Crocodylus niloticus, CNI), and bird and reptile genomes have been highly conserved during evolution. In this study, size and GC content of all chromosomes are measured from the flow karyotypes of GGA, TSC and CNI. Genome sizes estimated from the total chromosome size demonstrate that TSC and CNI are 1.21 Gb and 1.29 Gb, respectively. This refines previous overestimations and reveals similar genome sizes in chicken, turtle and crocodile. Analysis of chromosome GC content in each of these three species shows a higher GC content in smaller chromosomes than in larger chromosomes. This contrasts with mammals and squamates in which GC content does not correlate with chromosome size. These data suggest that a common ancestor of birds, turtles and crocodiles had a small genome size and a chromosomal size-dependent GC bias, distinct from the squamate lineage.

Unresolved orthology and peculiar coding sequence properties of lamprey genes: the KCNA gene family as test case

BACKGROUND

In understanding the evolutionary process of vertebrates, cyclostomes (hagfishes and lamprey) occupy crucial positions. Resolving molecular phylogenetic relationships of cyclostome genes with gnathostomes (jawed vertebrates) genes is indispensable in deciphering both the species tree and gene trees. However, molecular phylogenetic analyses, especially those including lamprey genes, have produced highly discordant results between gene families. To efficiently scrutinize this problem using partial genome assemblies of early vertebrates, we focused on the potassium voltage-gated channel, shaker-related (KCNA) family, whose members are mostly single-exon.

RESULTS

Seven sea lamprey KCNA genes as well as six elephant shark genes were identified, and their orthologies to bony vertebrate subgroups were assessed. In contrast to robustly supported orthology of the elephant shark genes to gnathostome subgroups, clear orthology of any sea lamprey gene could not be established. Notably, sea lamprey KCNA sequences displayed unique codon usage pattern and amino acid composition, probably associated with exceptionally high GC-content in their coding regions. This lamprey-specific property of coding sequences was also observed generally for genes outside this gene family.

CONCLUSIONS

Our results suggest that secondary modifications of sequence properties unique to the lamprey lineage may be one of the factors preventing robust orthology assessments of lamprey genes, which deserves further genome-wide validation. The lamprey lineage-specific alteration of protein-coding sequence properties needs to be taken into consideration in tackling the key questions about early vertebrate evolution.

Substitution patterns are GC-biased in divergent sequences across the metazoans

The fastest-evolving regions in the human and chimpanzee genomes show a remarkable excess of weak (A,T) to strong (G,C) nucleotide substitutions since divergence from their common ancestor. We investigated the phylogenetic extent and possible causes of this weak to strong (W→S) bias in divergent sequences (BDS) using recently sequenced genomes and recombination maps from eight trios of eukaryotic species. To quantify evidence for BDS, we inferred substitution histories using an efficient maximum likelihood approach with a context-dependent evolutionary model. We then annotated all lineage-specific substitutions in terms of W→S bias and density on the chromosomes. Finally, we used the inferred substitutions to calculate a BDS score—a log odds ratio between substitution type and density—and assessed its statistical significance with Fisher's exact test. Applying this approach, we found significant BDS in the coding and noncoding sequence of human, mouse, dog, stickleback, fruit fly, and worm. We also observed a significant lack of W→S BDS in chicken and yeast. The BDS score varies between species and across the chromosomes within each species. It is most strongly correlated with different genomic features in different species, but a strong correlation with recombination rates is found in several species. Our results demonstrate that a W→S substitution bias in fast-evolving sequences is a widespread phenomenon. The patterns of BDS observed suggest that a recombination-associated process, such as GC-biased gene conversion, is involved in the production of the bias in many species, but the strength of the BDS likely depends on many factors, including genome stability, variability in recombination rate over time and across the genome, the frequency of meiosis, and the amount of outcrossing in each species.

GC-biased gene conversion and selection affect GC content in the Oryza genus (rice)

Base composition varies among and within eukaryote genomes. Although mutational bias and selection have initially been invoked, more recently GC-biased gene conversion (gBGC) has been proposed to play a central role in shaping nucleotide landscapes, especially in yeast, mammals, and birds. gBGC is a kind of meiotic drive in favor of G and C alleles, associated with recombination. Previous studies have also suggested that gBGC could be at work in grass genomes. However, these studies were carried on third codon positions that can undergo selection on codon usage. As most preferred codons end in G or C in grasses, gBGC and selection can be confounded. Here we investigated further the forces that might drive GC content evolution in the rice genus using both coding and noncoding sequences. We found that recombination rates correlate positively with equilibrium GC content and that selfing species (Oryza sativa and O. glaberrima) have significantly lower equilibrium GC content compared with more outcrossing species. As recombination is less efficient in selfing species, these results suggest that recombination drives GC content. We also detected a positive relationship between expression levels and GC content in third codon positions, suggesting that selection favors codons ending with G or C bases. However, the correlation between GC content and recombination cannot be explained by selection on codon usage alone as it was also observed in noncoding positions. Finally, analyses of polymorphism data ruled out the hypothesis that genomic variation in GC content is due to mutational processes. Our results suggest that both gBGC and selection on codon usage affect GC content in the Oryza genus and likely in other grass species.

Contrasting GC-content dynamics across 33 mammalian genomes: relationship with life-history traits and chromosome sizes

The origin, evolution, and functional relevance of genomic variations in GC content are a long-debated topic, especially in mammals. Most of the existing literature, however, has focused on a small number of model species and/or limited sequence data sets. We analyzed more than 1000 orthologous genes in 33 fully sequenced mammalian genomes, reconstructed their ancestral isochore organization in the maximum likelihood framework, and explored the evolution of third-codon position GC content in representatives of 16 orders and 27 families. We showed that the previously reported erosion of GC-rich isochores is not a general trend. Several species (e.g., shrew, microbat, tenrec, rabbit) have independently undergone a marked increase in GC content, with a widening gap between the GC-poorest and GC-richest classes of genes. The intensively studied apes and (especially) murids do not reflect the general placental pattern. We correlated GC-content evolution with species life-history traits and cytology. Significant effects of body mass and genome size were detected, with each being consistent with the GC-biased gene conversion model.

Surprising fitness consequences of GC-biased gene conversion: I. Mutation load and inbreeding depression

GC-biased gene conversion (gBGC) is a recombination-associated process mimicking selection in favor of G and C alleles. It is increasingly recognized as a widespread force in shaping the genomic nucleotide landscape. In recombination hotspots, gBGC can lead to bursts of fixation of GC nucleotides and to accelerated nucleotide substitution rates. It was recently shown that these episodes of strong gBGC could give spurious signatures of adaptation and/or relaxed selection. There is also evidence that gBGC could drive the fixation of deleterious amino acid mutations in some primate genes. This raises the question of the potential fitness effects of gBGC. While gBGC has been metaphorically termed the "Achilles' heel" of our genome, we do not know whether interference between gBGC and selection merely has practical consequences for the analysis of sequence data or whether it has broader fundamental implications for individuals and populations. I developed a population genetics model to predict the consequences of gBGC on the mutation load and inbreeding depression. I also used estimates available for humans to quantitatively evaluate the fitness impact of gBGC. Surprising features emerged from this model: (i) Contrary to classical mutation load models, gBGC generates a fixation load independent of population size and could contribute to a significant part of the load; (ii) gBGC can maintain recessive deleterious mutations for a long time at intermediate frequency, in a similar way to overdominance, and these mutations generate high inbreeding depression, even if they are slightly deleterious; (iii) since mating systems affect both the selection efficacy and gBGC intensity, gBGC challenges classical predictions concerning the interaction between mating systems and deleterious mutations, and gBGC could constitute an additional cost of outcrossing; and (iv) if mutations are biased toward A and T alleles, very low gBGC levels can reduce the load. A robust prediction is that the gBGC level minimizing the load depends only on the mutational bias and population size. These surprising results suggest that gBGC may have nonnegligible fitness consequences and could play a significant role in the evolution of genetic systems. They also shed light on the evolution of gBGC itself.

Biased gene conversion and the evolution of mammalian genomic landscapes

Recombination is typically thought of as a symmetrical process resulting in large-scale reciprocal genetic exchanges between homologous chromosomes. Recombination events, however, are also accompanied by short-scale, unidirectional exchanges known as gene conversion in the neighborhood of the initiating double-strand break. A large body of evidence suggests that gene conversion is GC-biased in many eukaryotes, including mammals and human. AT/GC heterozygotes produce more GC- than AT-gametes, thus conferring a population advantage to GC-alleles in high-recombining regions. This apparently unimportant feature of our molecular machinery has major evolutionary consequences. Structurally, GC-biased gene conversion explains the spatial distribution of GC-content in mammalian genomes-the so-called isochore structure. Functionally, GC-biased gene conversion promotes the segregation and fixation of deleterious AT --> GC mutations, thus increasing our genomic mutation load. Here we review the recent evidence for a GC-biased gene conversion process in mammals, and its consequences for genomic landscapes, molecular evolution, and human functional genomics.

The GC-heterogeneity of teleost fishes

BACKGROUND

One of the most striking features of mammalian and birds chromosomes is the variation in the guanine-cytosine (GC) content that occurs over scales of hundreds of kilobases to megabases; this is known as the "isochore" structure. Among other vertebrates the presence of isochores depends upon the taxon; isochore are clearly present in Crocodiles and turtles but fish genome seems very homogeneous on GC content. This has suggested a unique isochore origin after the divergence between Sarcopterygii and Actinopterygii, but before that between Sauropsida and mammals. However during more than 30 years of analysis, isochore characteristics have been studied and many important biological properties have been associated with the isochore structure of human genomes. For instance, the genes are more compact and their density is highest in GC rich isochores.

RESULTS

This paper shows in teleost fish genomes the existence of "GC segmentation" sharing some of the characteristics of isochores although teleost fish genomes presenting a particular homogeneity in CG content. The entire genomes of T nigroviridis and D rerio are now available, and this has made it possible to check whether a mosaic structure associated with isochore properties can be found in these fishes. In this study, hidden Markov models were trained on fish genes (T nigroviridis and D rerio) which were classified by using the isochore class of their human orthologous. A clear segmentation of these genomes was detected.

CONCLUSION

The GC content is an excellent indicator of isochores in heterogeneous genomes as mammals. The segmentation we obtained were well correlated with GC content and other properties associated to GC content such as gene density, the number of exons per gene and the length of introns. Therefore, the GC content is the main property that allows the detection of isochore but more biological properties have to be taken into account. This method allows detecting isochores in homogeneous genomes.

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.

The impact of recombination on nucleotide substitutions in the human genome

Unraveling the evolutionary forces responsible for variations of neutral substitution patterns among taxa or along genomes is a major issue for detecting selection within sequences. Mammalian genomes show large-scale regional variations of GC-content (the isochores), but the substitution processes at the origin of this structure are poorly understood. We analyzed the pattern of neutral substitutions in 1 Gb of primate non-coding regions. We show that the GC-content toward which sequences are evolving is strongly negatively correlated to the distance to telomeres and positively correlated to the rate of crossovers (R² = 47%). This demonstrates that recombination has a major impact on substitution patterns in human, driving the evolution of GC-content. The evolution of GC-content correlates much more strongly with male than with female crossover rate, which rules out selectionist models for the evolution of isochores. This effect of recombination is most probably a consequence of the neutral process of biased gene conversion (BGC) occurring within recombination hotspots. We show that the predictions of this model fit very well with the observed substitution patterns in the human genome. This model notably explains the positive correlation between substitution rate and recombination rate. Theoretical calculations indicate that variations in population size or density in recombination hotspots can have a very strong impact on the evolution of base composition. Furthermore, recombination hotspots can create strong substitution hotspots. This molecular drive affects both coding and non-coding regions. We therefore conclude that along with mutation, selection and drift, BGC is one of the major factors driving genome evolution. Our results also shed light on variations in the rate of crossover relative to non-crossover events, along chromosomes and according to sex, and also on the conservation of hotspot density between human and chimp.

Mammalian genomes show a very strong heterogeneity of base composition along chromosomes (the so-called isochores). The functional significance of these peculiar genomic landscapes is highly debated: do isochores confer some selective advantage, or are they simply the by-product of neutral evolutionary processes? To resolve this issue, we analyzed the pattern of substitution in the human genome by comparison with chimpanzee and macaque. We show that the evolution of base composition (GC-content) is essentially determined by the rate of recombination. This effect appears to be much stronger in male than in female germline, which rules out selective explanations for the evolution of isochores. We show that this impact of recombination is most probably a consequence of the process of biased gene conversion (BGC). This neutral process mimics the action of selection and can induce strong substitution hotspots within recombination hotspots, sometimes leading to the fixation of deleterious mutations. BGC appears to be one of the major factors driving genome evolution. It is therefore essential to take this process into account if we want to be able to interpret genome sequences.

Noncanonical role of Hox14 revealed by its expression patterns in lamprey and shark

Hox genes are arranged in uninterrupted clusters in vertebrate genomes, and the nested patterns of their expression define spatial identities in multiple embryonic tissues. The ancestral Hox cluster of vertebrates has long been thought to consist of, maximally, 13 Hox genes. However, recently, Hox14 genes were discovered in three chordate lineages, the coelacanth, cartilaginous fishes, and amphioxus, but their expression patterns have not yet been analyzed. We isolated Hox14 cDNAs from the Japanese lamprey and cloudy catshark. These genes were not expressed in the central nervous systems, somites, or fin buds/folds but were expressed in a restricted cell population surrounding the hindgut. The lack of Hox14 expression in most of the embryonic axial elements, where nested Hox expressions define spatial identities, suggests a decoupling of Hox14 genes' regulation from the ancestral regulatory mechanism. The relaxation of preexisting constraint for collinear expression may have permitted the secondary losses of this Hox member in the tetrapod and teleost lineages.

Different origins of bird and reptile sex chromosomes inferred from comparative mapping of chicken Z-linked genes

Recent progress of chicken genome projects has revealed that bird ZW and mammalian XY sex chromosomes were derived from different autosomal pairs of the common ancestor; however, the evolutionary relationship between bird and reptilian sex chromosomes is still unclear. The Chinese soft-shelled turtle (Pelodiscus sinensis) exhibits genetic sex determination, but no distinguishable (heteromorphic) sex chromosomes have been identified. In order to investigate this further, we performed molecular cytogenetic analyses of this species, and thereby identified ZZ/ZW-type micro-sex chromosomes. In addition, we cloned reptile homologues of chicken Z-linked genes from three reptilian species, the Chinese soft-shelled turtle and the Japanese four-striped rat snake (Elaphe quadrivirgata), which have heteromorphic sex chromosomes, and the Siam crocodile (Crocodylus siamensis), which exhibits temperature-dependent sex determination and lacks sex chromosomes. We then mapped them to chromosomes of each species using FISH. The linkage of the genes has been highly conserved in all species: the chicken Z chromosome corresponded to the turtle chromosome 6q, snake chromosome 2p and crocodile chromosome 3. The order of the genes was identical among the three species. The absence of homology between the bird Z chromosome and the snake and turtle Z sex chromosomes suggests that the origin of the sex chromosomes and the causative genes of sex determination are different between birds and reptiles.

Time scale for cyclostome evolution inferred with a phylogenetic diagnosis of hagfish and lamprey cDNA sequences

The Cyclostomata consists of the two orders Myxiniformes (hagfishes) and Petromyzoniformes (lampreys), and its monophyly has been unequivocally supported by recent molecular phylogenetic studies. Under this updated vertebrate phylogeny, we performed in silico evolutionary analyses using currently available cDNA sequences of cyclostomes. We first calculated the GC-content at four-fold degenerate sites (GC(4)), which revealed that an extremely high GC-content is shared by all the lamprey species we surveyed, whereas no striking pattern in GC-content was observed in any of the hagfish species surveyed. We then estimated the timing of diversification in cyclostome evolution using nucleotide and amino acid sequences. We obtained divergence times of 470-390 million years ago (Mya) in the Ordovician-Silurian-Devonian Periods for the interordinal split between Myxiniformes and Petromyzoniformes; 90-60 Mya in the Cretaceous-Tertiary Periods for the split between the two hagfish subfamilies, Myxininae and Eptatretinae; 280-220 Mya in the Permian-Triassic Periods for the split between the two lamprey subfamilies, Geotriinae and Petromyzoninae; and 30-10 Mya in the Tertiary Period for the split between the two lamprey genera, Petromyzon and Lethenteron. This evolutionary configuration indicates that Myxiniformes and Petromyzoniformes diverged shortly after the common ancestor of cyclostomes split from the future gnathostome lineage. Our results also suggest that intra-subfamilial diversification in hagfish and lamprey lineages (especially those distributed in the northern hemisphere) occurred in the Cretaceous or Tertiary Periods.

Phylogenomics of nonavian reptiles and the structure of the ancestral amniote genome

We report results of a megabase-scale phylogenomic analysis of the Reptilia, the sister group of mammals. Large-scale end-sequence scanning of genomic clones of a turtle, alligator, and lizard reveals diverse, mammal-like landscapes of retroelements and simple sequence repeats (SSRs) not found in the chicken. Several global genomic traits, including distinctive phylogenetic lineages of CR1-like long interspersed elements (LINEs) and a paucity of A-T rich SSRs, characterize turtles and archosaur genomes, whereas higher frequencies of tandem repeats and a lower global GC content reveal mammal-like features in Anolis. Nonavian reptile genomes also possess a high frequency of diverse and novel 50-bp unit tandem duplications not found in chicken or mammals. The frequency distributions of approximately 65,000 8-mer oligonucleotides suggest that rates of DNA-word frequency change are an order of magnitude slower in reptiles than in mammals. These results suggest a diverse array of interspersed and SSRs in the common ancestor of amniotes and a genomic conservatism and gradual loss of retroelements in reptiles that culminated in the minimalist chicken genome. The sequences reported in this paper have been deposited in the GenBank database (accession nos. CZ 250707-CZ 257443 and DX 390731-DX 389174).

Tuatara (Sphenodon) Genomics: BAC Library Construction, Sequence Survey, and Application to the DMRT Gene Family

The tuatara (Sphenodon punctatus) is of "extraordinary biological interest" as the most distinctive surviving reptilian lineage (Rhyncocephalia) in the world. To provide a genomic resource for an understanding of genome evolution in reptiles, and as part of a larger project to produce genomic resources for various reptiles (evogen.jgi.doe.gov/second_levels/BACs/our_libraries.html), a large-insert bacterial artificial chromosome (BAC) library from a male tuatara was constructed. The library consists of 215 424 individual clones whose average insert size was empirically determined to be 145 kb, yielding a genomic coverage of approximately 6.3x. A BAC-end sequencing analysis of 121 420 bp of sequence revealed a genomic GC content of 46.8%, among the highest observed thus far for vertebrates, and identified several short interspersed repetitive elements (mammalian interspersed repeat-type repeats) and long interspersed repetitive elements, including chicken repeat 1 element. Finally, as a quality control measure the arrayed library was screened with probes corresponding to 2 conserved noncoding regions of the candidate sex-determining gene DMRT1 and the DM domain of the related DMRT2 gene. A deep coverage contig spanning nearly 300 kb was generated, supporting the deep coverage and utility of the library for exploring tuatara genomics.

Evidence for different origin of sex chromosomes in snakes, birds, and mammals and step-wise differentiation of snake sex chromosomes

All snake species exhibit genetic sex determination with the ZZ/ZW type of sex chromosomes. To investigate the origin and evolution of snake sex chromosomes, we constructed, by FISH, a cytogenetic map of the Japanese four-striped rat snake (Elaphe quadrivirgata) with 109 cDNA clones. Eleven of the 109 clones were localized to the Z chromosome. All human and chicken homologues of the snake Z-linked genes were located on autosomes, suggesting that the sex chromosomes of snakes, mammals, and birds were all derived from different autosomal pairs of the common ancestor. We mapped the 11 Z-linked genes of E. quadrivirgata to chromosomes of two other species, the Burmese python (Python molurus bivittatus) and the habu (Trimeresurus flavoviridis), to investigate the process of W chromosome differentiation. All and 3 of the 11 clones were localized to both the Z and W chromosomes in P. molurus and E. quadrivirgata, respectively, whereas no cDNA clones were mapped to the W chromosome in T. flavoviridis. Comparative mapping revealed that the sex chromosomes are only slightly differentiated in P. molurus, whereas they are fully differentiated in T. flavoviridis, and E. quadrivirgata is at a transitional stage of sex-chromosome differentiation. The differentiation of sex chromosomes was probably initiated from the distal region on the short arm of the protosex chromosome of the common ancestor, and then deletion and heterochromatization progressed on the sex-specific chromosome from the phylogenetically primitive boids to the more advanced viperids.