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

Unleashing diversity through flexibility: The evolutionary journey of sex chromosomes in amphibians and reptiles

Sex determination systems have greatly diversified between amphibians and reptiles, with such as the different sex chromosome compositions within a single species and transition between temperature-dependent sex determination (TSD) and genetic sex determination (GSD). In most sex chromosome studies on amphibians and reptiles, the whole-genome sequence of Xenopous tropicalis and chicken have been used as references to compare the chromosome homology of sex chromosomes among each of these taxonomic groups, respectively. In the present study, we reviewed existing reports on sex chromosomes, including karyotypes, in amphibians and reptiles. Furthermore, we compared the identified genetic linkages of sex chromosomes in amphibians and reptiles with the chicken genome as a reference, which is believed to resemble the ancestral tetrapod karyotype. Our findings revealed that sex chromosomes in amphibians are derived from genetic linkages homologous to various chicken chromosomes, even among several frogs within single families, such as Ranidae and Pipidae. In contrast, sex chromosomes in reptiles exhibit conserved genetic linkages with chicken chromosomes, not only across most species within a single family, but also within closely related families. The diversity of sex chromosomes in amphibians and reptiles may be attributed to the flexibility of their sex determination systems, including the ease of sex reversal in these animals.

Comparative genomics reveals the dynamics of chromosome evolution in Lepidoptera

Chromosomes are a central unit of genome organization. One-tenth of all described species on Earth are butterflies and moths, the Lepidoptera, which generally possess 31 chromosomes. However, some species display dramatic variation in chromosome number. Here we analyse 210 chromosomally complete lepidopteran genomes and show that the chromosomes of extant lepidopterans are derived from 32 ancestral linkage groups, which we term Merian elements. Merian elements have remained largely intact through 250 million years of evolution and diversification. Against this stable background, eight lineages have undergone extensive reorganization either through numerous fissions or a combination of fusion and fission events. Outside these lineages, fusions are rare and fissions are rarer still. Fusions often involve small, repeat-rich Merian elements and the sex-linked element. Our results reveal the constraints on genome architecture in Lepidoptera and provide a deeper understanding of chromosomal rearrangements in eukaryotic genome evolution.

Karyotype Diversification and Chromosome Rearrangements in Squamate Reptiles

Karyotype diversification represents an important, yet poorly understood, driver of evolution. Squamate reptiles are characterized by a high taxonomic diversity which is reflected at the karyotype level in terms of general structure, chromosome number and morphology, and insurgence of differentiated simple or multiple-sex-chromosome systems with either male or female heterogamety. The potential of squamate reptiles as unique model organisms in evolutionary cytogenetics has been recognised in recent years in several studies, which have provided novel insights into the chromosome evolutionary dynamics of different taxonomic groups. Here, we review and summarize the resulting complex, but promising, general picture from a systematic perspective, mapping some of the main squamate karyological characteristics onto their phylogenetic relationships. We highlight how all the major categories of balanced chromosome rearrangements contributed to the karyotype evolution in different taxonomic groups. We show that distinct karyotype evolutionary trends may occur, and coexist, with different frequencies in different clades. Finally, in light of the known squamate chromosome diversity and recent research advances, we discuss traditional and novel hypotheses on karyotype evolution and propose a scenario of circular karyotype evolution.

Conserved chromatin and repetitive patterns reveal slow genome evolution in frogs

Frogs are an ecologically diverse and phylogenetically ancient group of anuran amphibians that include important vertebrate cell and developmental model systems, notably the genus Xenopus. Here we report a high-quality reference genome sequence for the western clawed frog, Xenopus tropicalis, along with draft chromosome-scale sequences of three distantly related emerging model frog species, Eleutherodactylus coqui, Engystomops pustulosus, and Hymenochirus boettgeri. Frog chromosomes have remained remarkably stable since the Mesozoic Era, with limited Robertsonian (i.e., arm-preserving) translocations and end-to-end fusions found among the smaller chromosomes. Conservation of synteny includes conservation of centromere locations, marked by centromeric tandem repeats associated with Cenp-a binding surrounded by pericentromeric LINE/L1 elements. This work explores the structure of chromosomes across frogs, using a dense meiotic linkage map for X. tropicalis and chromatin conformation capture (Hi-C) data for all species. Abundant satellite repeats occupy the unusually long (~20 megabase) terminal regions of each chromosome that coincide with high rates of recombination. Both embryonic and differentiated cells show reproducible associations of centromeric chromatin and of telomeres, reflecting a Rabl-like configuration. Our comparative analyses reveal 13 conserved ancestral anuran chromosomes from which contemporary frog genomes were constructed.

Cross-species chromosome painting and repetitive DNA mapping illuminate the karyotype evolution in true crocodiles (Crocodylidae)

Crocodilians have maintained very similar karyotype structures and diploid chromosome numbers for around 100 million years, with only minor variations in collinearity. Why this karyotype structure has largely stayed unaltered for so long is unclear. In this study, we analyzed the karyotypes of six species belonging to the genera Crocodylus and Osteolaemus (Crocodylidae, true crocodiles), among which the Congolian endemic O. osborni was included and investigated. We utilized various techniques (differential staining, fluorescence in situ hybridization with repetitive DNA and rDNA probes, whole chromosome painting, and comparative genomic hybridization) to better understand how crocodile chromosomes evolved. We studied representatives of three of the four main diploid chromosome numbers found in crocodiles (2n = 30/32/38). Our data provided new information about the species studied, including the identification of four major chromosomal rearrangements that occurred during the karyotype diversification process in crocodiles. These changes led to the current diploid chromosome numbers of 2n = 30 (fusion) and 2n = 38 (fissions), derived from the ancestral state of 2n = 32. The conserved cytogenetic tendency in crocodilians, where extant species keep near-ancestral state, contrasts with the more dynamic karyotype evolution seen in other major reptile groups.

Chromosome evolution in Iberolacerta, a genus that deviates from the standard karyotype formula of Lacertidae

This paper describes the preparation of flow-sorted chromosome paints from the Iberian Rock lizard Iberolacerta monticola, exemplifying their subsequent use in cross-species comparisons of chromosome painting. We carried out comparative analyses of chromosome evolution in the congeneric species I. galani and I. bonnali, as well as in two other species of Lacertini (Lacerta schreiberi and Timon lepidus) whose sex chromosomes were also studied through comparative genomic hybridization. Most species of Lacertini possess a diplod number of 2n = 38, with 36 acrocentric macrochromosomes and 2 microchromosomes. However, the nine species included in the genus Iberolacerta do not possess microchromosomes. Furthermore, very conspicuous differences from the standard Lacertini karyotype were observed in the three Pyrenean species of this genus, which included several biarmed metacentrics and a Z1Z2W multiple sex-chromosome system. With the possible exception of L. schreiberi, all the species of the family Lacertidae described to date appear to share homologous Z chromosomes, which date back to the last common ancestor of the whole group. We provide conclusive evidence that L. schreiberi should no longer be considered an exception to this rule, and demonstrate that the loss of microchromosomes in Iberolacerta was produced by their fusion to a middle-sized chromosome. Furthermore, we show that the multiple sex-chromosome system of the Pyrenean species of Iberolacerta originated from the fusion of the ancestral W chromosome with one of the shortest autosomes, and provide additional evidence of the fast evolution of DNA sequences linked to the W chromosome in Lacertini.

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.

Three amphioxus reference genomes reveal gene and chromosome evolution of chordates

The slow-evolving invertebrate amphioxus has an irreplaceable role in advancing our understanding of the vertebrate origin and innovations. Here we resolve the nearly complete chromosomal genomes of three amphioxus species, one of which best recapitulates the 17 chordate ancestor linkage groups. We reconstruct the fusions, retention, or rearrangements between descendants of whole-genome duplications, which gave rise to the extant microchromosomes likely existed in the vertebrate ancestor. Similar to vertebrates, the amphioxus genome gradually establishes its three-dimensional chromatin architecture at the onset of zygotic activation and forms two topologically associated domains at the Hox gene cluster. We find that all three amphioxus species have ZW sex chromosomes with little sequence differentiation, and their putative sex-determining regions are nonhomologous to each other. Our results illuminate the unappreciated interspecific diversity and developmental dynamics of amphioxus genomes and provide high-quality references for understanding the mechanisms of chordate functional genome evolution.

Evolutionary analysis of a complete chicken genome

Microchromosomes are prevalent in nonmammalian vertebrates [P. D. Waters et al., Proc. Natl. Acad. Sci. U.S.A. 118 (2021)], but a few of them are missing in bird genome assemblies. Here, we present a new chicken reference genome containing all autosomes, a Z and a W chromosome, with all gaps closed except for the W. We identified ten small microchromosomes (termed dot chromosomes) with distinct sequence and epigenetic features, among which six were newly assembled. Those dot chromosomes exhibit extremely high GC content and a high level of DNA methylation and are enriched for housekeeping genes. The pericentromeric heterochromatin of dot chromosomes is disproportionately large and continues to expand with the proliferation of satellite DNA and testis-expressed genes. Our analyses revealed that the 41-bp CNM repeat frequently forms higher-order repeats (HORs) at the centromeres of acrocentric chromosomes. The centromere core regions where the kinetochore attaches often encompass telomeric sequence (TTAGGG)n, and in a one of the dot chromosomes, the centromere core recruits an endogenous retrovirus (ERV). We further demonstrate that the W chromosome shares some common features with dot chromosomes, having large arrays of hypermethylated tandem repeats. Finally, using the complete chicken chromosome models, we reconstructed a fine picture of chordate karyotype evolution, revealing frequent chromosomal fusions before and after vertebrate whole-genome duplications. Our sequence and epigenetic characterization of chicken chromosomes shed insights into the understanding of vertebrate genome evolution and chromosome biology.

Dinosaurs: Comparative Cytogenomics of Their Reptile Cousins and Avian Descendants

Reptiles known as dinosaurs pervade scientific and popular culture, while interest in their genomics has increased since the 1990s. Birds (part of the crown group Reptilia) are living theropod dinosaurs. Chromosome-level genome assemblies cannot be made from long-extinct biological material, but dinosaur genome organization can be inferred through comparative genomics of related extant species. Most reptiles apart from crocodilians have both macro- and microchromosomes; comparative genomics involving molecular cytogenetics and bioinformatics has established chromosomal relationships between many species. The capacity of dinosaurs to survive multiple extinction events is now well established, and birds now have more species in comparison with any other terrestrial vertebrate. This may be due, in part, to their karyotypic features, including a distinctive karyotype of around n = 40 (~10 macro and 30 microchromosomes). Similarity in genome organization in distantly related species suggests that the common avian ancestor had a similar karyotype to e.g., the chicken/emu/zebra finch. The close karyotypic similarity to the soft-shelled turtle (n = 33) suggests that this basic pattern was mostly established before the Testudine-Archosaur divergence, ~255 MYA. That is, dinosaurs most likely had similar karyotypes and their extensive phenotypic variation may have been mediated by increased random chromosome segregation and genetic recombination, which is inherently higher in karyotypes with more and smaller chromosomes.

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.

Evolution of the ancestral mammalian karyotype and syntenic regions

Computational reconstruction of ancestral mammalian karyotypes revealed a comprehensive picture of the chromosome rearrangements that occurred over the evolutionary history of mammals. Ancient gene order, in some cases extending to full chromosomes, was found conserved for more than 300 My, demonstrating strong evolutionary constraint against rearrangements in some regions. Conserved segments of chromosomes are enriched for genes that control developmental processes. Therefore, Darwinian selection likely maintains ancient gene combinations while allowing for genomic innovations within or near chromosomal sites that break and rearrange over evolutionary time. The revealed relationship between the three-dimensional structure of chromosomes and the evolutionary stability of chromosome segments provides additional insights into the mechanisms of chromosome evolution and diseases associated with genome rearrangements.

Decrypting the rearrangements that drive mammalian chromosome evolution is critical to understanding the molecular bases of speciation, adaptation, and disease susceptibility. Using 8 scaffolded and 26 chromosome-scale genome assemblies representing 23/26 mammal orders, we computationally reconstructed ancestral karyotypes and syntenic relationships at 16 nodes along the mammalian phylogeny. Three different reference genomes (human, sloth, and cattle) representing phylogenetically distinct mammalian superorders were used to assess reference bias in the reconstructed ancestral karyotypes and to expand the number of clades with reconstructed genomes. The mammalian ancestor likely had 19 pairs of autosomes, with nine of the smallest chromosomes shared with the common ancestor of all amniotes (three still conserved in extant mammals), demonstrating a striking conservation of synteny for ∼320 My of vertebrate evolution. The numbers and types of chromosome rearrangements were classified for transitions between the ancestral mammalian karyotype, descendent ancestors, and extant species. For example, 94 inversions, 16 fissions, and 14 fusions that occurred over 53 My differentiated the therian from the descendent eutherian ancestor. The highest breakpoint rate was observed between the mammalian and therian ancestors (3.9 breakpoints/My). Reconstructed mammalian ancestor chromosomes were found to have distinct evolutionary histories reflected in their rates and types of rearrangements. The distributions of genes, repetitive elements, topologically associating domains, and actively transcribed regions in multispecies homologous synteny blocks and evolutionary breakpoint regions indicate that purifying selection acted over millions of years of vertebrate evolution to maintain syntenic relationships of developmentally important genes and regulatory landscapes of gene-dense chromosomes.

Expression and Characterization of the Spats1 Gene and Its Response to E2/MT Treatment in the Chinese Soft-Shelled Turtle (Pelodiscus sinensis)

Spats1 (spermatogenesis-associated, serinerich 1) has been characterized as a male-biased gene which acts an important role in the germ cell differentiation of mammals. Nevertheless, the function of Spats1 in the Chinese soft-shelled turtle (P. sinensis) has not yet been reported. To initially explore the expression of Spats1 in P. sinensis and its response to sex steroid treatment, we cloned the CDS of Spats1 for the first time and analyzed its expression profile in different tissues, including the testes in different seasons. The Spats1 cDNA fragment is 1201 base pairs (bp) in length and contains an open reading frame (ORF) of 849 bp, which codes for 283 amino acids. Spats1 mRNA was highly expressed in the testes (p < 0.01) and barely detectable in other tissues. In P. sinensis, the relative expression of Spats1 also responsive to seasonal changes in testis development. In summer (July) and autumn (October), Spats1 gene expression was significantly higher in the testes than in other seasons (p < 0.05). Spats1 mRNA was found to be specifically expressed in germ cells by chemical in situ hybridization (CISH), and it was mainly located in primary spermatocytes (Sc1), secondary spermatocytes (Sc2) and spermatozoa (St). Spats1 expression in embryos was not significantly changed after 17α-methyltestosterone (MT)and 17β-estradiol (E2) treatment. In adults, MT significantly induced Spats1 expression in male P. sinensis. However, the expression of Spats1 in testes was not responsive to E2 treatment. In addition, the expression of Spats1 in females was not affected by either MT or E2 treatment. These results imply that Spats1 is a male-specific expressed gene that is mainly regulated by MT and is closely linked to spermatogenesis and release in P. sinensis.

Sequencing of laser captured Z and W chromosomes of the tocantins paradoxical frog (Pseudis tocantins) provides insights on repeatome and chromosomal homology

Pseudis tocantins is the only frog species of the hylid genus Pseudis that possesses highly heteromorphic sex chromosomes. Z and W chromosomes of Ps. tocantins differ in size, morphology, position of the nucleolar organizer region (NOR) and the amount and distribution of heterochromatin. A chromosomal inversion and heterochromatin amplification on the W chromosome were previously inferred to be involved in the evolution of this sex chromosome pair. Despite these findings, knowledge related to the molecular composition of the large heterochromatic band of this W chromosome is restricted to the PcP190 satellite DNA, and no data are available regarding the gene content of either the W or the Z chromosome of Ps. tocantins. Here, we sequenced microdissected Z and W chromosomes of this species to further resolve their molecular composition. Comparative genomic analysis suggests that Ps. tocantins sex chromosomes are likely homologous to chromosomes 4 and 10 of Xenopus tropicalis. Analyses of the repetitive DNA landscape in the Z and W assemblies allowed for the identification of several transposable elements and putative satellite DNA sequences. Finally, some transposable elements from the W assembly were found to be highly diverse and divergent from elements found elsewhere in the genome, suggesting a rapid amplification of these elements on the W chromosome.

Microchromosomes are building blocks of bird, reptile, and mammal chromosomes

Microchromosomes, once considered unimportant shreds of the chicken genome, are gene-rich elements with a high GC content and few transposable elements. Their origin has been debated for decades. We used cytological and whole-genome sequence comparisons, and chromosome conformation capture, to trace their origin and fate in genomes of reptiles, birds, and mammals. We find that microchromosomes as well as macrochromosomes are highly conserved across birds and share synteny with single small chromosomes of the chordate amphioxus, attesting to their origin as elements of an ancient animal genome. Turtles and squamates (snakes and lizards) share different subsets of ancestral microchromosomes, having independently lost microchromosomes by fusion with other microchromosomes or macrochromosomes. Patterns of fusions were quite different in different lineages. Cytological observations show that microchromosomes in all lineages are spatially separated into a central compartment at interphase and during mitosis and meiosis. This reflects higher interaction between microchromosomes than with macrochromosomes, as observed by chromosome conformation capture, and suggests some functional coherence. In highly rearranged genomes fused microchromosomes retain most ancestral characteristics, but these may erode over evolutionary time; surprisingly, de novo microchromosomes have rapidly adopted high interaction. Some chromosomes of early-branching monotreme mammals align to several bird microchromosomes, suggesting multiple microchromosome fusions in a mammalian ancestor. Subsequently, multiple rearrangements fueled the extraordinary karyotypic diversity of therian mammals. Thus, microchromosomes, far from being aberrant genetic elements, represent fundamental building blocks of amniote chromosomes, and it is mammals, rather than reptiles and birds, that are atypical.

Remnant of Unrelated Amniote Sex Chromosomal Linkage Sharing on the Same Chromosome in House Gecko Lizards, Providing a Better Understanding of the Ancestral Super-Sex Chromosome

Comparative chromosome maps investigating sex chromosomal linkage groups in amniotes and microsatellite repeat motifs of a male house gecko lizard (Hemidactylus frenatus, HFR) and a flat-tailed house gecko lizard (H. platyurus, HPL) of unknown sex were examined using 75 bacterial artificial chromosomes (BACs) from chicken and zebra finch genomes. No massive accumulations of microsatellite repeat motifs were found in either of the gecko lizards, but 10 out of 13 BACs mapped on HPL chromosomes were associated with other amniote sex chromosomes. Hybridization of the same BACs onto multiple different chromosome pairs suggested transitions to sex chromosomes across amniotes. No BAC hybridization signals were found on HFR chromosomes. However, HFR diverged from HPL about 30 million years ago, possibly due to intrachromosomal rearrangements occurring in the HFR lineage. By contrast, heterochromatin likely reshuffled patterns between HPL and HFR, as observed from C-positive heterochromatin distribution. Six out of ten BACs showed partial homology with squamate reptile chromosome 2 (SR2) and snake Z and/or W sex chromosomes. The gecko lizard showed shared unrelated sex chromosomal linkages-the remnants of a super-sex chromosome. A large ancestral super-sex chromosome showed a correlation between SR2 and snake W sex chromosomes.

Lizards as Model Organisms of Sex Chromosome Evolution: What We Really Know from a Systematic Distribution of Available Data?

Lizards represent unique model organisms in the study of sex determination and sex chromosome evolution. Among tetrapods, they are characterized by an unparalleled diversity of sex determination systems, including temperature-dependent sex determination (TSD) and genetic sex determination (GSD) under either male or female heterogamety. Sex chromosome systems are also extremely variable in lizards. They include simple (XY and ZW) and multiple (X1X2Y and Z1Z2W) sex chromosome systems and encompass all the different hypothesized stages of diversification of heterogametic chromosomes, from homomorphic to heteromorphic and completely heterochromatic sex chromosomes. The co-occurrence of TSD, GSD and different sex chromosome systems also characterizes different lizard taxa, which represent ideal models to study the emergence and the evolutionary drivers of sex reversal and sex chromosome turnover. In this review, we present a synthesis of general genome and karyotype features of non-snakes squamates and discuss the main theories and evidences on the evolution and diversification of their different sex determination and sex chromosome systems. We here provide a systematic assessment of the available data on lizard sex chromosome systems and an overview of the main cytogenetic and molecular methods used for their identification, using a qualitative and quantitative approach.

Why Do Some Vertebrates Have Microchromosomes?

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

Whole-chromosome fusions in the karyotype evolution of Sceloporus (Iguania, Reptilia) are more frequent in sex chromosomes than autosomes

Whole-chromosome fusions play a major role in the karyotypic evolution of reptiles. It has been suggested that certain chromosomes tend to fuse with sex chromosomes more frequently than others. However, the comparative genomic synteny data are too scarce to draw strong conclusions. We obtained and sequenced chromosome-specific DNA pools of Sceloporus malachiticus, an iguanian species which has experienced many chromosome fusions. We found that four of seven lineage-specific fusions involved sex chromosomes, and that certain syntenic blocks which constitute the sex chromosomes, such as the homologues of the Anolis carolinensis chromosomes 11 and 16, are repeatedly involved in sex chromosome formation in different squamate species. To test the hypothesis that the karyotypic shift could be associated with changes in recombination patterns, we performed a synaptonemal complex analysis in this species and in Sceloporus variabilis (2n = 34). It revealed that the sex chromosomes in S. malachiticus had two distal pseudoautosomal regions and a medial differentiated region. We found that multiple fusions little affected the recombination rate in S. malachiticus. Our data confirm more frequent involvement of certain chromosomes in sex chromosome formation, but do not reveal a connection between the gonosome-autosome fusions and the evolution of recombination rate. This article is part of the theme issue 'Challenging the paradigm in sex chromosome evolution: empirical and theoretical insights with a focus on vertebrates (Part II)'.

Reconstruction of proto-vertebrate, proto-cyclostome and proto-gnathostome genomes provides new insights into early vertebrate evolution

Ancient polyploidization events have had a lasting impact on vertebrate genome structure, organization and function. Some key questions regarding the number of ancient polyploidization events and their timing in relation to the cyclostome-gnathostome divergence have remained contentious. Here we generate de novo long-read-based chromosome-scale genome assemblies for the Japanese lamprey and elephant shark. Using these and other representative genomes and developing algorithms for the probabilistic macrosynteny model, we reconstruct high-resolution proto-vertebrate, proto-cyclostome and proto-gnathostome genomes. Our reconstructions resolve key questions regarding the early evolutionary history of vertebrates. First, cyclostomes diverged from the lineage leading to gnathostomes after a shared tetraploidization (1R) but before a gnathostome-specific tetraploidization (2R). Second, the cyclostome lineage experienced an additional hexaploidization. Third, 2R in the gnathostome lineage was an allotetraploidization event, and biased gene loss from one of the subgenomes shaped the gnathostome genome by giving rise to remarkably conserved microchromosomes. Thus, our reconstructions reveal the major evolutionary events and offer new insights into the origin and evolution of vertebrate genomes.

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.

Universal nomenclature for oxytocin-vasotocin ligand and receptor families

Oxytocin (OXT; hereafter OT) and arginine vasopressin or vasotocin (AVP or VT; hereafter VT) are neurotransmitter ligands that function through specific receptors to control diverse functions^1,2. Here we performed genomic analyses on 35 species that span all major vertebrate lineages, including newly generated high-contiguity assemblies from the Vertebrate Genomes Project^3,4. Our findings support the claim⁵ that OT (also known as OXT) and VT (also known as AVP) are adjacent paralogous genes that have resulted from a local duplication, which we infer was through DNA transposable elements near the origin of vertebrates and in which VT retained more of the parental sequence. We identified six major oxytocin-vasotocin receptors among vertebrates. We propose that all six of these receptors arose from a single receptor that was shared with the common ancestor of invertebrates, through a combination of whole-genome and large segmental duplications. We propose a universal nomenclature based on evolutionary relationships for the genes that encode these receptors, in which the genes are given the same orthologous names across vertebrates and paralogous names relative to each other. This nomenclature avoids confusion due to differential naming in the pre-genomic era and incomplete genome assemblies, furthers our understanding of the evolution of these genes, aids in the translation of findings across species and serves as a model for other gene families.

Analysis of Paralogons, Origin of the Vertebrate Karyotype, and Ancient Chromosomes Retained in Extant Species

A manually curated set of ohnolog families has been assembled, for seven species of bony vertebrates, that includes 255 four-member families and 631 three-member families, encompassing over 2,900 ohnologs. Across species, the patterns of chromosomes upon which the ohnologs reside fall into 17 distinct categories. These 17 paralogons reflect the 17 ancestral chromosomes that existed in our chordate ancestor immediately prior to the two rounds of whole-genome duplication (2R-WGD) that occurred around 600 Ma. Within each paralogon, it has now been possible to assign those pairs of ohnologs that diverged from each other at the first round of duplication, through analysis of the molecular phylogeny of four-member families. Comparison with another recent analysis has identified four apparently incorrect assignments of pairings following 2R, along with several omissions, in that study. By comparison of the patterns between paralogons, it has also been possible to identify nine chromosomal fusions that occurred between 1R and 2R, and three chromosomal fusions that occurred after 2R, that generated an ancestral bony-vertebrate karyotype comprising 47 chromosomes. At least 27 of those ancestral bony-vertebrate chromosomes can, in some extant species, be shown not to have undergone any fusion or fission events. Such chromosomes are here termed “archeochromosomes,” and have each survived essentially unchanged in their content of genes for some 400 Myr. Their utility lies in their potential for tracking the various fusion and fission events that have occurred in different lineages throughout the expansion of bony vertebrates.

A new emu genome illuminates the evolution of genome configuration and nuclear architecture of avian chromosomes

Emu and other ratites are more informative than any other birds in reconstructing the evolution of the ancestral avian or vertebrate karyotype because of their much slower rate of genome evolution. Here, we generated a new chromosome-level genome assembly of a female emu, and estimated the tempo of chromosome evolution across major avian phylogenetic branches, by comparing it to chromosome-level genome assemblies of 11 other bird and one turtle species. We found ratites exhibited the lowest numbers of intra- and inter-chromosomal changes among birds since their divergence with turtles. The small-sized and gene-rich emu microchromosomes have frequent inter-chromosomal contacts that are associated with housekeeping genes, which appears to be driven by clustering their centromeres in the nuclear interior, away from the macrochromosomes in the nuclear periphery. Unlike nonratite birds, only less than one-third of the emu W Chromosome regions have lost homologous recombination and diverged between the sexes. The emu W is demarcated into a highly heterochromatic region (WS0) and another recently evolved region (WS1) with only moderate sequence divergence with the Z Chromosome. WS1 has expanded its inactive chromatin compartment, increased chromatin contacts within the region, and decreased contacts with the nearby regions, possibly influenced by the spreading of heterochromatin from WS0. These patterns suggest that alteration of chromatin conformation comprises an important early step of sex chromosome evolution. Overall, our results provide novel insights into the evolution of avian genome structure and sex chromosomes in three-dimensional space.

Vertebrate Chromosome Evolution

The study of chromosome evolution is undergoing a resurgence of interest owing to advances in DNA sequencing technology that facilitate the production of chromosome-scale whole-genome assemblies de novo. This review focuses on the history, methods, discoveries, and current challenges facing the field, with an emphasis on vertebrate genomes. A detailed examination of the literature on the biology of chromosome rearrangements is presented, specifically the relationship between chromosome rearrangements and phenotypic evolution, adaptation, and speciation. A critical review of the methods for identifying, characterizing, and visualizing chromosome rearrangements and computationally reconstructing ancestral karyotypes is presented. We conclude by looking to the future, identifying the enormous technical and scientific challenges presented by the accumulation of hundreds and eventually thousands of chromosome-scale assemblies.

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

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

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.

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

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

Karyotypic Evolution of Sauropsid Vertebrates Illuminated by Optical and Physical Mapping of the Painted Turtle and Slider Turtle Genomes

Recent sequencing and software enhancements have advanced our understanding of the evolution of genomic structure and function, especially addressing novel evolutionary biology questions. Yet fragmentary turtle genome assemblies remain a challenge to fully decipher the genetic architecture of adaptive evolution. Here, we use optical mapping to improve the contiguity of the painted turtle (Chrysemys picta) genome assembly and use de novo fluorescent in situ hybridization (FISH) of bacterial artificial chromosome (BAC) clones, BAC-FISH, to physically map the genomes of the painted and slider turtles (Trachemys scripta elegans). Optical mapping increased C. picta's N50 by ~242% compared to the previous assembly. Physical mapping permitted anchoring ~45% of the genome assembly, spanning 5544 genes (including 20 genes related to the sex determination network of turtles and vertebrates). BAC-FISH data revealed assembly errors in C. picta and T. s. elegans assemblies, highlighting the importance of molecular cytogenetic data to complement bioinformatic approaches. We also compared C. picta's anchored scaffolds to the genomes of other chelonians, chicken, lizards, and snake. Results revealed a mostly one-to-one correspondence between chromosomes of painted and slider turtles, and high homology among large syntenic blocks shared with other turtles and sauropsids. Yet, numerous chromosomal rearrangements were also evident across chelonians, between turtles and squamates, and between avian and non-avian reptiles.

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

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

Chromosome Painting Does Not Support a Sex Chromosome Turnover in Lacerta agilis Linnaeus, 1758

Reptiles show a remarkable diversity of sex determination mechanisms and sex chromosome systems, derived from different autosomal pairs. The origin of the ZW sex chromosomes of Lacerta agilis, a widespread Eurasian lizard species, is a matter of discussion: is it a small macrochromosome from the 11-18 group common to all lacertids, or does this species have a unique ZW pair derived from the large chromosome 5? Using independent molecular cytogenetic methods, we investigated the karyotype of L. agilis exigua from Siberia, Russia, to identify the sex chromosomes. FISH with a flow-sorted chromosome painting probe derived from L. strigata and specific to chromosomes 13, 14, and Z confirmed that the Z chromosome of L. agilis is a small macrochromosome, the same as in L. strigata. FISH with the telomeric probe showed an extensive accumulation of the telomere-like repeat in the W chromosome in agreement with previous studies, excluding the possibility that the lineages of L. agilis studied in different works could have different sex chromosome systems due to a putative intra-species polymorphism. Our results reinforce the idea of the stability of the sex chromosomes and lack of evidence for sex-chromosome turnovers in known species of Lacertidae.

Evolution of Marsupial Genomes

Marsupial genomes, which are packaged into large chromosomes, provide a powerful resource for studying the mechanisms of genome evolution. The extensive and valuable body of work on marsupial cytogenetics, combined more recently with genome sequence data, has enabled prediction of the 2n = 14 karyotype ancestral to all marsupial families. The application of both chromosome biology and genome sequencing, or chromosomics, has been a necessary approach for various aspects of mammalian genome evolution, such as understanding sex chromosome evolution and the origin and evolution of transmissible tumors in Tasmanian devils. The next phase of marsupial genome evolution research will employ chromosomics approaches to begin addressing fundamental questions in marsupial genome evolution and chromosome evolution more generally. The answers to these complex questions will impact our understanding across a broad range of fields, including the genetics of speciation, genome adaptation to environmental stressors, and species management.

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.

Time lapse: A glimpse into prehistoric genomics

For the purpose of this review, ‘time-lapse’ refers to the reconstruction of ancestral (in this case dinosaur) karyotypes using genome assemblies of extant species. Such reconstructions are only usually possible when genomes are assembled to ‘chromosome level’ i.e. a complete representation of all the sequences, correctly ordered contiguously on each of the chromosomes. Recent paleontological evidence is very clear that birds are living dinosaurs, the latest example of dinosaurs emerging from a catastrophic extinction event. Non-avian dinosaurs (ever present in the public imagination through art, and broadcast media) emerged some 240 million years ago and have displayed incredible phenotypic diversity. Here we report on our recent studies to infer the overall karyotype of the Theropod dinosaur lineage from extant avian chromosome level genome assemblies. Our work first focused on determining the likely karyotype of the avian ancestor (most likely a chicken-sized, two-legged, feathered, land dinosaur from the Jurassic period) finding karyotypic similarity to the chicken. We then took the work further to determine the likely karyotype of the bird-lizard ancestor and the chromosomal changes (chiefly translocations and inversions) that occurred between then and modern birds. A combination of bioinformatics and cross-species fluorescence in situ hybridization (zoo-FISH) uncovered a considerable number of translocations and fissions from a ‘lizard-like’ genome structure of 2n = 36–46 to one similar to that of soft-shelled turtles (2n = 66) from 275 to 255 million years ago (mya). Remarkable karyotypic similarities between some soft-shelled turtles and chicken suggests that there were few translocations from the bird-turtle ancestor (plus ∼7 fissions) through the dawn of the dinosaurs and pterosaurs, through the theropod linage and on to most to modern birds. In other words, an avian-like karyotype was in place about 240mya when the dinosaurs and pterosaurs first emerged. We mapped 49 chromosome inversions from then to the present day, uncovering some gene ontology enrichment in evolutionary breakpoint regions. This avian-like karyotype with its many (micro)chromosomes provides the basis for variation (the driver of natural selection) through increased random segregation and recombination. It may therefore contribute to the ability of dinosaurs to survive multiple extinction events, emerging each time as speciose and diverse.

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

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

Reconstruction of avian ancestral karyotypes reveals differences in the evolutionary history of macro- and microchromosomes

Background

Reconstruction of ancestral karyotypes is critical for our understanding of genome evolution, allowing for the identification of the gross changes that shaped extant genomes. The identification of such changes and their time of occurrence can shed light on the biology of each species, clade and their evolutionary history. However, this is impeded by both the fragmented nature of the majority of genome assemblies and the limitations of the available software to work with them. These limitations are particularly apparent in birds, with only 10 chromosome-level assemblies reported thus far. Algorithmic approaches applied to fragmented genome assemblies can nonetheless help define patterns of chromosomal change in defined taxonomic groups.

Results

Here, we make use of the DESCHRAMBLER algorithm to perform the first large-scale study of ancestral chromosome structure and evolution in birds. This algorithm allows us to reconstruct the overall genome structure of 14 key nodes of avian evolution from the Avian ancestor to the ancestor of the Estrildidae, Thraupidae and Fringillidae families.

Conclusions

Analysis of these reconstructions provides important insights into the variability of rearrangement rates during avian evolution and allows the detection of patterns related to the chromosome distribution of evolutionary breakpoint regions. Moreover, the inclusion of microchromosomes in our reconstructions allows us to provide novel insights into the evolution of these avian chromosomes, specifically.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1544-8) contains supplementary material, which is available to authorized users.

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

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

Reconstruction of the diapsid ancestral genome permits chromosome evolution tracing in avian and non-avian dinosaurs

Genomic organisation of extinct lineages can be inferred from extant chromosome-level genome assemblies. Here, we apply bioinformatic and molecular cytogenetic approaches to determine the genomic structure of the diapsid common ancestor. We then infer the events that likely occurred along this lineage from theropod dinosaurs through to modern birds. Our results suggest that most elements of a typical ‘avian-like’ karyotype (40 chromosome pairs, including 30 microchromosomes) were in place before the divergence of turtles from birds ~255 mya. This genome organisation therefore predates the emergence of early dinosaurs and pterosaurs and the evolution of flight. Remaining largely unchanged interchromosomally through the dinosaur–theropod route that led to modern birds, intrachromosomal changes nonetheless reveal evolutionary breakpoint regions enriched for genes with ontology terms related to chromatin organisation and transcription. This genomic structure therefore appears highly stable yet contributes to a large degree of phenotypic diversity, as well as underpinning adaptive responses to major environmental disruptions via intrachromosomal repatterning.

Ancient diapsids diverged into the lineages leading to turtles and birds over 250 million years ago. Here, the authors use genomic and molecular cytogenetic analyses of modern species to infer the genome structure of the diapsid common ancestor (DCA) and the changes occurring along the lineage to birds through theropod dinosaurs.

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

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

Amniogenic somatopleure: a novel origin of multiple cell lineages contributing to the cardiovascular system

The somatopleure is the amniotic primordium in amniote development, but its boundary to the embryonic body at early embryonic stages and the fate of cells constituting this structure are not well characterized. It also remains unclear how cells behave during the demarcation between intra- and extra-embryonic tissues. Here we identify cellular alignments, which indicate two streams towards the sites of dorsal amniotic closure and ventral thoracic wall formation. A subpopulation of mesodermal cells moving ventrally from the somatopleural region adjacent to the base of the head fold enter the body of the embryo and distribute to the thoracic wall, pharyngeal arches and heart. These cells are induced to differentiate into vascular endothelial cells and cardiomyocytes possibly by FGF and BMP signaling, respectively. These results indicate that the somatopleure acting as the amniotic primordium also serves as a source of embryonic cells, which may contribute to cardiovascular development.

Immunocytological analysis of meiotic recombination in two anole lizards (Squamata, Dactyloidae)

Although the evolutionary importance of meiotic recombination is not disputed, the significance of interspecies differences in the recombination rates and recombination landscapes remains under-appreciated. Recombination rates and distribution of chiasmata have been examined cytologically in many mammalian species, whereas data on other vertebrates are scarce. Immunolocalization of the protein of the synaptonemal complex (SYCP3), centromere proteins and the mismatch-repair protein MLH1 was used, which is associated with the most common type of recombination nodules, to analyze the pattern of meiotic recombination in the male of two species of iguanian lizards, Anolis carolinensis Voigt, 1832 and Deiroptyx coelestinus (Cope, 1862). These species are separated by a relatively long evolutionary history although they retain the ancestral iguanian karyotype. In both species similar and extremely uneven distributions of MLH1 foci along the macrochromosome bivalents were detected: approximately 90% of crossovers were located at the distal 20% of the chromosome arm length. Almost total suppression of recombination in the intermediate and proximal regions of the chromosome arms contradicts the hypothesis that "homogenous recombination" is responsible for the low variation in GC content across the anole genome. It also leads to strong linkage disequilibrium between the genes located in these regions, which may benefit conservation of co-adaptive gene arrays responsible for the ecological adaptations of the anoles.

Genome evolution in the allotetraploid frog Xenopus laevis

To explore the origins and consequences of tetraploidy in the African clawed frog, we sequenced the Xenopus laevis genome and compared it to the related diploid X. tropicalis genome. We demonstrate the allotetraploid origin of X. laevis by partitioning its genome into two homeologous subgenomes, marked by distinct families of “fossil” transposable elements. Based on the activity of these elements and the age of hundreds of unitary pseudogenes, we estimate that the two diploid progenitor species diverged ~34 million years ago (Mya) and combined to form an allotetraploid ~17–18 Mya. 56% of all genes are retained in two homeologous copies. Protein function, gene expression, and the amount of flanking conserved sequence all correlate with retention rates. The subgenomes have evolved asymmetrically, with one chromosome set more often preserving the ancestral state and the other experiencing more gene loss, deletion, rearrangement, and reduced gene expression.

Microchromosome polymorphism in the sand lizard, Lacerta agilis Linnaeus, 1758 (Reptilia, Squamata)

Most true lizards (Lacertidae) share a conservative karyotype, consisting of 18 pairs of macrochromosomes and one microchromosome pair. Homeologues of the microchromosome are present in other squamates and even in chickens. No structural autosomal microchromosome polymorphisms have been described previously in lizards. We found homozygous and heterozygous carriers of a microchromosome variant in a Siberian population of the sand lizard, Lacerta agilis Linnaeus, 1758. The variant microchromosome was almost twice as long as the standard one. In heterozygotes at pachytene, the microchromosomes firstly pair in proximal regions and the central part of the longer axial element undergoes foldback synapsis, then its distal region pairs with the distal region of the standard partner. At metaphase-I, the heteromorphic microchromosome bivalents have a proximal chiasma. The content of the additional segment was Ag-NOR, C-like DAPI, CMA3 negative. FISH with telomere PNA probe did not detect interstitial (TTAGGG)_n sequences in the heteromorphic and any other bivalents. Both homo- and heterozygous carriers were phenotypically normal. The presence of homozygotes shows that heterozygotes are fertile. Reduction in the number of microchromosomes is a clear trend in squamate evolution, as a result of microchromosomes fusing together or with macrochromosomes. Our findings indicate that gaining additional DNA may lead to a transformation of microchromosomes into small macrochromosomes without fusion.

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

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

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

Genomic Landscape of Long Terminal Repeat Retrotransposons (LTR-RTs) and Solo LTRs as Shaped by Ectopic Recombination in Chicken and Zebra Finch

Transposable elements (TEs) are nearly ubiquitous among eukaryotic genomes, but TE contents vary dramatically among phylogenetic lineages. Several mechanisms have been proposed as drivers of TE dynamics in genomes, including the fixation/loss of a particular TE insertion by selection or drift as well as structural changes in the genome due to mutation (e.g., recombination). In particular, recombination can have a significant and directional effect on the genomic TE landscape. For example, ectopic recombination removes internal regions of long terminal repeat retrotransposons (LTR-RTs) as well as one long terminal repeat (LTR), resulting in a solo LTR. In this study, we focus on the intra-species dynamics of LTR-RTs and solo LTRs in bird genomes. The distribution of LTR-RTs and solo LTRs in birds is intriguing because avian recombination rates vary widely within a given genome. We used published linkage maps and whole genome assemblies to study the relationship between recombination rates and LTR-removal events in the chicken and zebra finch. We hypothesized that regions with low recombination rates would harbor more full-length LTR-RTs (and fewer solo LTRs) than regions with high recombination rates. We tested this hypothesis by comparing the ratio of full-length LTR-RTs and solo LTRs across chromosomes, across non-overlapping megabase windows, and across physical features (i.e., centromeres and telomeres). The chicken data statistically supported the hypothesis that recombination rates are inversely correlated with the ratio of full-length to solo LTRs at both the chromosome level and in 1-Mb non-overlapping windows. We also found that the ratio of full-length to solo LTRs near chicken telomeres was significantly lower than those ratios near centromeres. Our results suggest a potential role of ectopic recombination in shaping the chicken LTR-RT genomic landscape.