Characterization of a satellite DNA in the genera Lacerta andTimon(Reptilia, Lacertidae) and its role in the differentiation of the W chromosome

Variants of a major DNA satellite discriminate parental subgenomes in a hybrid parthenogenetic lizard Darevskia unisexualis (Darevsky, 1966)

Hybrid parthenogenetic animals are an exceptionally interesting model for studying the mechanisms and evolution of sexual and asexual reproduction. A diploid parthenogenetic lizard Darevskia unisexualis is a result of an ancestral cross between a maternal species Darevskia raddei nairensis and a paternal species Darevskia valentini and presents a unique opportunity for a cytogenetic and computational analysis of a hybrid karyotype. Our previous results demonstrated a significant divergence between the pericentromeric DNA sequences of the parental Darevskia species; however, an in-depth comparative study of their pericentromeres is still lacking. Here, using target sequencing of microdissected pericentromeric regions, we reveal and compare the repertoires of the pericentromeric tandem repeats of the parental Darevskia lizards. We found species-specific sequences of the major pericentromeric tandem repeat CLsat, which allowed computational prediction and experimental validation of fluorescent DNA probes discriminating parental chromosomes within the hybrid karyotype of D. unisexualis. Moreover, we have implemented a generalizable computational method, based on the optimization of the Levenshtein distance between tandem repeat monomers, for finding species-specific fluorescent probes for pericentromere staining. In total, we anticipate that our comparative analysis of Darevskia pericentromeric repeats, the species-specific fluorescent probes that we found and the pipeline that we developed will form a basis for the future detailed cytogenomic studies of a wide range of natural and laboratory hybrids.

Cytogenetic Analysis of Satellitome of Madagascar Leaf-Tailed Geckos

Satellite DNA (satDNA) consists of sequences of DNA that form tandem repetitions across the genome, and it is notorious for its diversity and fast evolutionary rate. Despite its importance, satDNA has been only sporadically studied in reptile lineages. Here, we sequenced genomic DNA and PCR-amplified microdissected W chromosomes on the Illumina platform in order to characterize the monomers of satDNA from the Henkel's leaf-tailed gecko U. henkeli and to compare their topology by in situ hybridization in the karyotypes of the closely related Günther's flat-tail gecko U. guentheri and gold dust day gecko P. laticauda. We identified seventeen different satDNAs; twelve of them seem to accumulate in centromeres, telomeres and/or the W chromosome. Notably, centromeric and telomeric regions seem to share similar types of satDNAs, and we found two that seem to accumulate at both edges of all chromosomes in all three species. We speculate that the long-term stability of all-acrocentric karyotypes in geckos might be explained from the presence of specific satDNAs at the centromeric regions that are strong meiotic drivers, a hypothesis that should be further tested.

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.

Conservation of Major Satellite DNAs in Snake Heterochromatin

Repetitive DNA sequences constitute a sizeable portion of animal genomes, and tandemly organized satellite DNAs are a major part of them. They are usually located in constitutive heterochromatin clusters in or near the centromeres or telomeres, and less frequently in the interstitial parts of chromosome arms. They are also frequently accumulated in sex chromosomes. The function of these clusters is to sustain the architecture of the chromosomes and the nucleus, and to regulate chromosome behavior during mitosis and meiosis. The study of satellite DNA diversity is important for understanding sex chromosome evolution, interspecific hybridization, and speciation. In this work, we identified four satellite DNA families in the genomes of two snakes from different families: Daboia russelii (Viperidae) and Pantherophis guttatus (Colubridae) and determine their chromosomal localization. We found that one family is localized in the centromeres of both species, whereas the others form clusters in certain chromosomes or subsets of chromosomes. BLAST with snake genome assemblies showed the conservation of such clusters, as well as a subtle presence of the satellites in the interspersed manner outside the clusters. Overall, our results show high conservation of satellite DNA in snakes and confirm the "library" model of satellite DNA evolution.

Decoding sex: Elucidating sex determination and how high-quality genome assemblies are untangling the evolutionary dynamics of sex chromosomes

Sexual reproduction is a diverse and widespread process. In gonochoristic species, the differentiation of sexes occurs through diverse mechanisms, influenced by environmental and genetic factors. In most vertebrates, a master-switch gene is responsible for triggering a sex determination network. However, only a few genes have acquired master-switch functions, and this process is associated with the evolution of sex-chromosomes, which have a significant influence in evolution. Additionally, their highly repetitive regions impose challenges for high-quality sequencing, even using high-throughput, state-of-the-art techniques. Here, we review the mechanisms involved in sex determination and their role in the evolution of species, particularly vertebrates, focusing on sex chromosomes and the challenges involved in sequencing these genomic elements. We also address the improvements provided by the growth of sequencing projects, by generating a massive number of near-gapless, telomere-to-telomere, chromosome-level, phased assemblies, increasing the number and quality of sex-chromosome sequences available for further studies.

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.

Amplified Fragments of an Autosome-Borne Gene Constitute a Significant Component of the W Sex Chromosome of Eremias velox (Reptilia, Lacertidae)

Heteromorphic W and Y sex chromosomes often experience gene loss and heterochromatinization, which is frequently viewed as their "degeneration". However, the evolutionary trajectories of the heterochromosomes are in fact more complex since they may not only lose but also acquire new sequences. Previously, we found that the heterochromatic W chromosome of a lizard Eremias velox (Lacertidae) is decondensed and thus transcriptionally active during the lampbrush stage. To determine possible sources of this transcription, we sequenced DNA from a microdissected W chromosome sample and a total female DNA sample and analyzed the results of reference-based and de novo assembly. We found a new repetitive sequence, consisting of fragments of an autosomal protein-coding gene ATF7IP2, several SINE elements, and sequences of unknown origin. This repetitive element is distributed across the whole length of the W chromosome, except the centromeric region. Since it retained only 3 out of 10 original ATF7IP2 exons, it remains unclear whether it is able to produce a protein product. Subsequent studies are required to test the presence of this element in other species of Lacertidae and possible functionality. Our results provide further evidence for the view of W and Y chromosomes as not just "degraded" copies of Z and X chromosomes but independent genomic segments in which novel genetic elements may arise.

The Role of Satellite DNAs in Genome Architecture and Sex Chromosome Evolution in Crambidae Moths

Tandem repeats are important parts of eukaryotic genomes being crucial e.g., for centromere and telomere function and chromatin modulation. In Lepidoptera, knowledge of tandem repeats is very limited despite the growing number of sequenced genomes. Here we introduce seven new satellite DNAs (satDNAs), which more than doubles the number of currently known lepidopteran satDNAs. The satDNAs were identified in genomes of three species of Crambidae moths, namely Ostrinia nubilalis, Cydalima perspectalis, and Diatraea postlineella, using graph-based computational pipeline RepeatExplorer. These repeats varied in their abundance and showed high variability within and between species, although some degree of conservation was noted. The satDNAs showed a scattered distribution, often on both autosomes and sex chromosomes, with the exception of both satellites in D. postlineella, in which the satDNAs were located at a single autosomal locus. Three satDNAs were abundant on the W chromosomes of O. nubilalis and C. perspectalis, thus contributing to their differentiation from the Z chromosomes. To provide background for the in situ localization of the satDNAs, we performed a detailed cytogenetic analysis of the karyotypes of all three species. This comparative analysis revealed differences in chromosome number, number and location of rDNA clusters, and molecular differentiation of sex chromosomes.

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

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

High dynamism for neo-sex chromosomes: satellite DNAs reveal complex evolution in a grasshopper

A common characteristic of sex chromosomes is the accumulation of repetitive DNA, which accounts for their diversification and degeneration. In grasshoppers, the X0 sex-determining system in males is considered ancestral. However, in some species, derived variants like neo-XY in males evolved several times independently by Robertsonian translocation. This is the case of Ronderosia bergii, in which further large pericentromeric inversion in the neo-Y also took place, making this species particularly interesting for investigating sex chromosome evolution. Here, we characterized the satellite DNAs (satDNAs) and transposable elements (TEs) of the species to investigate the quantitative differences in repeat composition between male and female genomes putatively associated with sex chromosomes. We found a total of 53 satDNA families and 56 families of TEs. The satDNAs were 13.5% more abundant in males than in females, while TEs were just 1.02% more abundant in females. These results imply differential amplification of satDNAs on neo-Y chromosome and a minor role of TEs in sex chromosome differentiation. We showed highly differentiated neo-XY sex chromosomes owing to major amplification of satDNAs in neo-Y. Furthermore, chromosomal mapping of satDNAs suggests high turnover of neo-sex chromosomes in R. bergii at the intrapopulation level, caused by multiple paracentric inversions, amplifications, and transpositions. Finally, the species is an example of the action of repetitive DNAs in the generation of variability for sex chromosomes after the suppression of recombination, and helps understand sex chromosome evolution at the intrapopulation level.

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.

Interstitial Telomeric Repeats Are Rare in Turtles

Telomeres are nucleoprotein complexes protecting chromosome ends in most eukaryotic organisms. In addition to chromosome ends, telomeric-like motifs can be accumulated in centromeric, pericentromeric and intermediate (i.e., between centromeres and telomeres) positions as so-called interstitial telomeric repeats (ITRs). We mapped the distribution of (TTAGGG)_n repeats in the karyotypes of 30 species from nine families of turtles using fluorescence in situ hybridization. All examined species showed the expected terminal topology of telomeric motifs at the edges of chromosomes. We detected ITRs in only five species from three families. Combining our and literature data, we inferred seven independent origins of ITRs among turtles. ITRs occurred in turtles in centromeric positions, often in several chromosomal pairs, in a given species. Their distribution does not correspond directly to interchromosomal rearrangements. Our findings support that centromeres and non-recombining parts of sex chromosomes are very dynamic genomic regions, even in turtles, a group generally thought to be slowly evolving. However, in contrast to squamate reptiles (lizards and snakes), where ITRs were found in more than half of the examined species, and birds, the presence of ITRs is generally rare in turtles, which agrees with the expected low rates of chromosomal rearrangements and rather slow karyotype evolution in this group.

Cytogenetic mechanisms of unisexuality in rock lizards

Darevskia rock lizards is a unique complex taxa, including more than thirty species, seven of which are parthenogenetic. In mixed populations of Darevskia lizards, tri- and tetraploid forms can be found. The most important issues in the theory of reticulate evolution of Darevskia lizards are the origin of parthenogenetic species and their taxonomic position. However, there is little data on how meiosis proceeds in these species. The present work reports the complex results of cytogenetics in a diploid parthenogenetic species – D. unisexualis. Here we detail the meiotic prophase I progression and the specific features оf mitotic chromosomes organization. The stages of meiosis prophase I were investigated by immunocytochemical analysis of preparations obtained from isolated primary oocytes of D. unisexualis in comparison with maternal species D. raddei nairensis. It has been shown that in D. unisexualis at the leptotene-zygotene stages the axial elements and the synaptonemal complex (SC) form typical “bouquets”. At the pachytene-diplotene stage, 18 autosomal SC-bivalents and thickened asynapted sex Z and w univalents were observed. The presence of SYCP1 protein between the lateral elements of autosomal chromosomes proved the formation of assembled SCs. Comparative genomic hybridization (CGH) on the mitotic metaphase chromosomes of D. unisexualis was carried out using the genomic DNA isolated from the parental species D. raddei nairensis and D. valentini. In the pericentromeric regions of half of the mitotic chromosomes of D. unisexualis, specific regions inherited from maternal species have been found. Following our results, we suggest a model for diploid germ cells formation from diploid oocytes without premeiotic duplication of chromosomes in the oogenesis of diploid parthenogenetic lizards D. unisexualis. Taken as a whole, our findings confirm the hybrid nature of D. unisexualis and shed light on heterozygosity and automixis in diploid parthenogenetic forms.

Evolutionary Variability of W-Linked Repetitive Content in Lacertid Lizards

Lacertid lizards are a widely radiated group of squamate reptiles with long-term stable ZZ/ZW sex chromosomes. Despite their family-wide homology of Z-specific gene content, previous cytogenetic studies revealed significant variability in the size, morphology, and heterochromatin distribution of their W chromosome. However, there is little evidence about the accumulation and distribution of repetitive content on lacertid chromosomes, especially on their W chromosome. In order to expand our knowledge of the evolution of sex chromosome repetitive content, we examined the topology of telomeric and microsatellite motifs that tend to often accumulate on the sex chromosomes of reptiles in the karyotypes of 15 species of lacertids by fluorescence in situ hybridization (FISH). The topology of the above-mentioned motifs was compared to the pattern of heterochromatin distribution, as revealed by C-banding. Our results show that the topologies of the examined motifs on the W chromosome do not seem to follow a strong phylogenetic signal, indicating independent and species-specific accumulations. In addition, the degeneration of the W chromosome can also affect the Z chromosome and potentially also other parts of the genome. Our study provides solid evidence that the repetitive content of the degenerated sex chromosomes is one of the most evolutionary dynamic parts of the genome.

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.

The Karyotype of Salvator merianae (Squamata, Teiidae): Analyses by Classical and Molecular Cytogenetic Techniques

In this study, we analyzed the karyotype of Salvator merianae (Teiidae) from the Brazilian semiarid region using different cytogenetic markers. Chromosomes were examined by classical (Giemsa and AgNOR staining) and molecular (FISH with ribosomal, telomeric, and microsatellite probes) cytogenetic approaches. S. merianae showed a diploid chromosome number of 2n = 38 (10 biarmed macrochromosomes + 28 microchromosomes). No sex-linked chromosome heteromorphisms were observed. Clusters of 18S/28S rDNA were localized in the terminal region of the long arm of pair 2. In addition to the typical telomeric signals, (TTAGGG)n repeats were detected in the pericentromeric region of some macrochromosome pairs, which might indicate the occurrence of chromosomal rearrangements via chromosome fusions. Hybridization signals of the microsatellite probes (GA)n, (GAA)n, and (GAG)n were uniformly distributed across all chromosomes, while (CA)n, (CAA)n, and (CAC)n produced brighter signals in the telomeric and pericentromeric regions of specific chromosome pairs. The comparison with previous studies demonstrates that, despite the wide distribution of S. merianae, the macrostructure organization of the karyotype remained unchanged, showing stability in diploid number and chromosome morphology.

Isolation and characterization of two satellite DNAs in Atlantolacerta andreanskyi (Werner, 1929) (Reptilia, Lacertidae)

Two satellite DNAs (satDNAs) have been isolated and characterized from three populations of Atlantolacerta andreanskyi. One satDNA (AAN-TaqI) has been isolated here from the first time. It is characterized by a tendency to AT enrichment (AT = 54.2%) and monomer length ranging from 187 to 199 bp. FISH experiments showed that this element occurs in subterminal position on the short arms of all chromosomes of the complement. The analyses of genetic variability of AAN-TaqI showed that the concerted evolution is acting effectively on these repeats that form separate clusters consistent with the geographic origin in the phylogenetic tree, thus supporting the hypothesis that A. andreanskyi would be a species complex. In addition, in the population from Jbel Aoulime this satDNA is already differentiated into two subfamilies. The other satDNA belongs to the family of IMO-TaqI already isolated in other lacertids. Differently from AAN-TaqI, concerted evolution does not seem to act effectively on this element that is not differentiated between populations. These results confirm that IMO-TaqI (AT = 53.4%) is conserved in both chromosomal position and most of its sequence in the lacertids from which it has been characterized so far. Its remarkable evolutionary conservation for about 45 million years could indicate that this satDNA may have a functional role that future investigations could unveil. Once again, this study shows how satDNAs coexisting in the same genome may differ in their evolutionary pattern, even though the reasons underlying this phenomenon in the species here studied have still to be fully understood.

Diversity of PBI-DdeI satellite DNA in snakes correlates with rapid independent evolution and different functional roles

To better understand PBI-DdeI satellite DNA located in the centromeric region of python, molecular evolution analysis was conducted on 40 snake species. A ladder-like pattern of DNA bands with repetition of the 194–210 bp monomer was observed in 15 species using PCR. Molecular cloning was performed to obtain 97 AT-rich monomer sequences. Phylogenetic and network analyses showed three PBI-DdeI subfamilies with sequences grouped in species-specific clusters, suggesting rapid evolution. Slow evolution was found in eight species with shared PBI-DdeI sequences, suggesting recent species diversification, allowing PBI-DdeI no time to diverge, with limited homogenization and fixation processes. Quantitative real-time PCR showed large differences in copy number between Python bivittatus and other snakes, consistent with repeat scanning of whole genome sequences. Copy numbers were significantly higher in female Naja kaouthia than in males, concurring with chromosomal distribution of PBI-DdeI specifically localized to female W chromosomes. PBI-DdeI might act as an evolutionary driver with several repeats to promote W chromosome differentiation and heterochromatinization in N. kaouthia. Analysis revealed PBI-DdeI with a reduced copy number, compared to P. bivittatus, in most snakes studied, and it is possible that it subsequently dispersed and amplified on W chromosomes with different functional roles in N. kaouthia.