Next Article in Journal
Imprecise Medicine: BRCA2 Variants of Uncertain Significance (VUS), the Challenges and Benefits to Integrate a Functional Assay Workflow with Clinical Decision Rules
Next Article in Special Issue
Multi-Scale Organization of the Drosophila melanogaster Genome
Previous Article in Journal
Evolutionary Dynamics of Wild Populations
Previous Article in Special Issue
Bridging the Gap between Vertebrate Cytogenetics and Genomics with Single-Chromosome Sequencing (ChromSeq)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 5
10.3390/genes12050779
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Artem Lisachov
1,2,*,
Daria Andreyushkova
3,
Guzel Davletshina
2,3,
Dmitry Prokopov
3,
Svetlana Romanenko
3,
Svetlana Galkina
4,
Alsu Saifitdinova
5,
Evgeniy Simonov
1,
Pavel Borodin
2,6 and
Vladimir Trifonov
3,6
1
Institute of Environmental and Agricultural Biology (X-BIO), University of Tyumen, Lenina str. 23, 625003 Tyumen, Russia
2
Institute of Cytology and Genetics SB RAS, Acad. Lavrentiev Ave. 10, 630090 Novosibirsk, Russia
3
Institute of Molecular and Cellular Biology SB RAS, Acad. Lavrentiev Ave. 8/2, 630090 Novosibirsk, Russia
4
Department of Genetics and Biotechnology, Saint Petersburg State University, Universitetskaya Emb. 7–9, 199034 Saint Petersburg, Russia
5
Department of Human and Animal Anatomy and Physiology, Herzen State Pedagogical University of Russia, Moyka Emb. 48, 191186 Saint Petersburg, Russia
6
Novosibirsk State University, Pirogova str. 3, 630090 Novosibirsk, Russia
*
Author to whom correspondence should be addressed.
Genes 2021, 12(5), 779; https://doi.org/10.3390/genes12050779
Submission received: 2 April 2021 / Revised: 11 May 2021 / Accepted: 17 May 2021 / Published: 20 May 2021
(This article belongs to the Special Issue Chromosome-Centric View of the Genome Organization and Evolution)
Browse Figures
Versions Notes

Abstract

:
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.

1. Introduction

Sex chromosomes originate from autosomal pairs when one chromosome acquires a dominant sex-determining locus. In some cases, the product of this locus interacts directly with the molecular physiological pathways of gonadal development, such as in mammals and in the African clawed frog (Xenopus laevis) [1,2]. In other cases, this locus represents a null allele and acts via haploinsufficiency mechanisms in heterozygotes, such as in birds [3]. The emergence of sex chromosomes is frequently followed by a suppression of recombination between the homologs in the heterogametic sex. The non-recombining segments gradually spread from a narrow zone around the sex determination locus to almost the whole length of the chromosomes, except the small pseudoautosomal region required for proper segregation of the sex chromosomes during meiosis [4]. This pathway was followed by a majority of sufficiently old vertebrate sex chromosome systems, although there are some systems that do not show significant spreading of the non-recombining regions over millions of years [5].
There are different explanations for this suppression of recombination. According to the most popular concept, nascent sex chromosomes tend to accumulate sexually antagonistic genes (i.e., useful for one sex but deleterious for the other), and inhibition of recombination is favored because it prevents disruption of the linkage between the sex-determining locus and the sexually antagonistic alleles [6,7]. Indeed, it is known that the Y chromosomes of mammals carry genes that are involved in testis development (for example, ZFY) [8], and sex chromosomes of guppies (Poecilia reticulata) carry genes responsible for male coloration [9,10]. In the subgenus Micropoecilia, which is closely related to guppies, male coloration is almost completely determined by the Y chromosome haplotype, and, thus, there is a limited number of distinct male color morphs linked to specific Y haplotypes [11].
The suppression of recombination between the sex chromosomes leads to a divergence between the X/Z and Y/W chromosome sequences over evolutionary time. In particular, the non-recombining Y and W chromosomes show loss of protein-coding genes, heterochromatinization, and repetitive sequence accumulation, which is generally understood as their “degradation” [12,13]. However, this term oversimplifies the evolutionary trajectories of Y and W chromosomes, which in many taxa show complex genetic content. In addition to the sexually antagonistic genes, they may contain preserved or neofunctionalized ancestral genes, homologous to the genes on X and Z chromosomes (gametologs), rDNA gene clusters, and other functional sequences [14,15,16,17,18]. Studies of genetic contents of Y and W chromosomes are important to understand their underestimated functional role in sex-specific developmental and physiological processes.
Currently, the best studied Y and W chromosomes are those of mammals and birds, and, to a lesser extent, snakes [15,17,19,20,21,22,23,24]. There are also some important studies of sex chromosomes of fish and amphibians: for example, sex-determining genes are known for medaka (XY) [25] and the African clawed frog (ZW) [2]. Sex chromosome structures and patterns of evolutionary transitions are also extensively studied in other species of frogs [26,27,28]. An extensive body of literature exists on the sex chromosomes of guppy fish [29,30,31,32]. However, for many other vertebrate Y and W chromosomes, especially in reptiles, it is only known that they do not have certain genes that their X and Z counterparts possess [33,34], and that they accumulate certain simple repeat sequences [35].
Among the vertebrate groups with insufficiently studied “degenerated” sex chromosomes, the true, or lacertid, lizards (Lacertidae) are particularly interesting. First, their sex chromosomes are partially homologous to the sex chromosomes of therian mammals [36], although the lacertids are female heterogametic (ZZ/ZW) [37]. The sex chromosomes of lacertids are homologous across the whole family, but in outgroup taxa, such as teiids and amphisbaenians, the genomic segments homologous to the lacertid ZZ/ZW chromosomes remain autosomal [38].
Second, their W chromosomes are extremely diverse in size and genetic content, in contrast with Z chromosomes. W chromosomes of various lacertid species contain telomeric repeats, microsatellites, and IMO-TaqI satellites, which are usually located in pericentromeric heterochromatic blocks on lacertid autosomes [39,40,41].
Recently, we studied sex chromosome behavior in meiotic prophase at the lampbrush stage in a lacertid Eremias velox (rapid racerunner), and obtained a microdissected W-chromosome painting probe to identify the sex bivalent at the lampbrush chromosome spread [42]. We found that the E. velox W chromosome at the lampbrush stage is more decondensed than autosomes and possesses numerous lateral DNA loops. The lateral loops of lampbrush chromosomes are proven to be sites of transcriptional activity [43,44,45]. This is in sharp contrast to the bird W chromosome, which is highly condensed and silenced at the lampbrush stage [46]. Considering the fact that the bulk of lampbrush lateral loops is represented by sequences that do not encode proteins [45], the W chromosome of E.velox should be saturated with these elements (including tandem and interspersed repeats).
The genetic content of the W chromosome of E. velox was investigated previously using FISH with probes to diverse microsatellite sequences, and the accumulation of several satellite types was detected [47]. However, this approach is not suitable for uncovering unique features of W chromosomes, since it is limited to checking the presence or absence of a priori known sequences.
To investigate the genetic content of the E. velox W chromosome in more detail and identify possible sources of the transcription, in this study we sequenced the previously obtained W-specific probe, and analyzed it both by aligning it to the reference genome of a related species, the European wall lizard (Podarcis muralis), and by de novo assembly. Since microdissected chromosome probes usually contain only partial chromosomes (due to material loss during microdissection and to amplification artifacts), we also sequenced the total genomic DNA of a female E. velox. Since the W chromosome in this species is relatively large [42], W-derived reads should constitute a significant part of the total genomic sequences.

2. Materials and Methods

The samples, the methods of the preparation of mitotic and lampbrush chromosome samples, the extraction of genomic DNA, and the production of the W-specific microdissected probe were described in detail previously [42]. Briefly, mitotic chromosome spreads were prepared from fibroblast cell cultures obtained from muscle and connective tissues of female E. velox as described previously [48,49]. The W chromosome sample was microdissected from a Giemsa-stained metaphase plate, using an Olympus IX-51 microscope equipped with an Eppendorf Transferman NK2 micromanipulator, and amplified using the GenomePlex Whole Genome Amplification (WGA-1) kit (Sigma, Darmstadt, Germany). The lampbrush chromosome spreads were obtained from previtellogenic and early vitellogenic oocytes, using the standard avian lampbrush technique described by Saifitdinova et al. [50]. The total genomic DNA was extracted from ethanol-preserved blood, using the phenol–chloroform technique [51]. All manipulations with live animals and euthanasia were approved by the Saint Petersburg State University Ethics Committee (statement #131-03-2) and the Institute of Molecular and Cellular Biology Ethics Committee (statement #01/18 from 05.03.2018).
The libraries of the W-specific probe and the total genomic DNA of a female E. velox were prepared using the TruSeq Nano DNA Low Throughput Library Prep Kit (Illumina, San Diego, USA) following the manufacturer’s protocol and sequenced on the Illumina MiSeq platform with the MiSeq Reagent Kit v3 (600 cycles). Before the sequencing, the libraries were examined for insert size and adapter content by electrophoresis in 2% agarose gel, and the molar concentrations were measured via qPCR. The sequence pre-processing steps included the trimming of Illumina adapters, WGA primers and low-quality bases, and discarding low-quality and too short (below 50 bp) reads with CLC Genomic Workbench (QIAGEN, Venlo, Netherlands) software. The sequencing reads were deposited in the NCBI SRA database under the accession number PRJNA716491.
Mapping the reads from the W-specific sample to the reference genome of P. muralis (GCF_004329235.1) was performed via Unipro UGENE v. 36 (https://www.ugene.net/, accessed on 1 September 2019) using the BWA-MEM algorithm [52] with default parameters and independently via the dopseq pipeline [53]. The regions of the reference genome with high coverage were searched via the UGENE graphical interface and using the dopseq output. The consensus sequences of the region of interest were extracted using both UGENE and the CLC Genomic Workbench. The total genomic DNA reads were aligned to the P. muralis genome as described above, and the previously identified sequences, which were present in the W-specific sample, were analyzed for coverage, and consensus sequences were extracted.
De novo assemblies of the W chromosome sample and the total DNA sample were generated via CLC Genomic Workbench using the ‘de novo assembly’ tool with default parameters. The W chromosome contigs were searched and aligned via the CLC Genomic Workbench using BLASTn search with the reference-based consensus sequence as query and the de novo assemblies as databases. Then, the identified W chromosome contigs were characterized via NCBI BLASTn through the web interface using the “nr/nt” database for all taxa and the “refseq_genomes” database with taxon specified as “Lacertidae” (i.e., the assemblies GCF_009819535.1 for Lacerta agilis, GCF_004329235.1 for P. muralis, and GCF_011800845.1 for Zootoca vivipara). The obtained W chromosome contigs were deposited in GenBank (MW846328-MW846339). To reconstruct the whole repeat unit, the identified contigs, which contain partial overlaps, were aligned to each other using CLC and MEGA X (https://www.megasoftware.net/, accessed on 1 October 2020) programs.
Since the data on lacertid repeated sequences is limited, RepeatModeler 2.0.1 [54] with the “-LTRStruct” option was used for the creation of a de novo library of repeated elements in the genomes of L. agilis (GCF_009819535.1), P. muralis (GCF_004329235.1), and Z. vivipara (GCF_011800845.1). The resulting fasta files with classified repeats were clustered in a single file “lacertalib.fa” using cd-hit-est 4.8.1 [55] with default settings. Annotation of repeated DNA elements in the assembled contigs of E. velox using the de novo library obtained in the previous step was carried out using RepeatMasker 4.1.1 (http://www.repeatmasker.org/, accessed on 1 February 2021) with the parameters “-e rmblast -s -no_is -lib lacertalib.fa”.
To confirm that the identified sequences indeed belong to the W chromosome, we performed FISH with a PCR-amplified fragment of the reference-based consensus sequence (homologous to positions 46,390,515—46,391,408 of chromosome 14 in P. muralis assembly GCF_004329235.1). For this purpose, we designed the following primers: EVWF 5′-CCTGTAGGGTTCACCGTTCC-3′; EVWR 5′-GGGTTGAGTGACCTTCTCGG-3′. The primers were designed using the NCBI online tool Primer-BLAST (https://www.ncbi.nlm.nih.gov/tools/primer-blast/, accessed on 27 September 2019) based on the reference-based consensus sequence. We used the total genomic DNA of a female E. velox as a template. The PCR mix contained 10 μL of a PCR master mix (Biolabmix, Novosibirsk, Russia), 7 μL of MilliQ water, 1μL of 10 μM solutions of each primer, and 1 μL of the genomic DNA solution (50 ng/µL). The PCR conditions were as follows: initial denaturation at 95 °C for 5 min, 26 cycles of amplification (95 °C for 30 s, 55 °C for 30 s, 72 °C for 1 min), 72 °C for 5 min, and hold at 4 °C. The PCR results were checked by electrophoresis in 1% agarose gel. To prepare a fluorescently labeled amplicon, we carried out a secondary PCR round, using the amplicon as a template, with an addition of 1 μL of 1 mM TAMRA-dUTP (Biosan, Novosibirsk, Russia) to 20 µL of PCR solution. FISH with the amplified probe was carried out using standard techniques [56]. The preparations were visualized with an Axioplan 2 imaging microscope (Carl Zeiss, Oberkochen, Germany ) equipped with a CCD camera (CV M300, JAI), CHROMA filter sets, and ISIS4 image processing package (MetaSystems GmbH). The brightness and contrast of all images were enhanced using Corel PaintShop Photo Pro X6 (Corel Corp, Ottawa, ON, Canada).

3. Results

3.1. Reference-Based Assembly

In total, 1,575,398 paired-end (PE) reads were obtained for the W-specific DNA sample and 2,797,722 PE reads for the total DNA sample. After the pre-processing, 859,598 PE reads and 1,828,412 PE reads were retained for the respective samples. The dopseq analysis of the W-specific sample reads revealed exceptionally high coverage (up to 24,878×, mean 3376×) of the fragments with a total length of 6170 bp inside a larger fragment of 17,213 bp length, situated inside the ATF7IP2 gene on chromosome 14 (Supplementary Table S1, Figure 1a). The alignment of total DNA reads to P. muralis genome yielded high coverage in the same region, and although the maximum coverage was lower (up to 848×), the coverage graph was more even (Figure 1b). No significant homology between the W chromosome of E. velox and the Z chromosome was detected by the dopseq analysis.

3.2. De Novo Assembly

The de novo assembly statistics for the sequenced DNA samples were as follows: for the W-specific sample, N50 = 321, contig lengths 200–5386 bp, and 1438 contigs in total; for the total DNA sample, N50 = 360, contig lengths 200–9509 bp, and 250,177 contigs in total. From each assembly, we selected the contigs with the best BLAST alignment scores to the reference-based consensus sequence, which resulted in 2 sets of 6 contigs (12 contigs in total). The contigs from the W-specific sample assembly had a length distribution from 501 to 4126 bp and a total length of 11,035 bp. The contigs from the total DNA sample assembly had a length distribution from 485 to 3309 bp and a total length of 10,658 bp.
Mutual alignment of partially overlapping contigs from the two sets resulted in a final sequence 10,628 bp in length. Thus, the fragments derived from the ATF7IP2 gene constitute a little more than half of the repeat unit, and other parts have different origins. The NCBI BLAST search did not retrieve significant hits for most of these non-ATF7IP2 fragments, with few notable exceptions. The search in the lacertid genomic databases retrieved multiple hits for three small fragments of 200–400 bp. One of them is originally present in the ATF7IP2 gene, whereas two are not (Figure 2a). The search in the “nr/nt” databases identified one of these non-ATF7IP2 fragments as a SINE mobile element of the POM/Squam1 family, whereas two other repetitive fragments showed similarity to various unrelated lacertid microsatellite and mRNA sequences. The RepeatMasker search identified all three segments as SINE repeats. Interestingly, both assemblies had contigs that overlap both the beginning and the end of the reference-based sequence; however, there was a gap in this sequence that was not covered by any contigs (Figure 2b).
The ATF7IP2-derived fragments on the W chromosome of E. velox include both intron and exon sequences, although only 3 of the 10 original exons are present (Figure 1). The nucleotide analysis of the exon sequences showed one 3 bp insertion in one of the exons, in comparison with the reference genomes, but no premature stop codons.

3.3. FISH with the Probe to the Fragment of ATF7IP2

PCR amplification of the fragment of ATF7IP2 resulted in a single band in the agarose gel, which corresponded to the predicted length. FISH with the probe derived from this amplicon resulted in an intense signal across almost the whole length of the W chromosome, except the centromeric region (Figure 3).

4. Discussion

Translocations of autosomal sequences onto sex chromosomes are widespread in vertebrates. In humans and other placental mammals, the sex chromosomes are constituted of the ancestral X/Y chromosome region and the added autosome-derived region [57]. In reptiles, such cases are known, for example, in geckos, anoles, and fence lizards [58,59,60,61]. In lacertids, the most well-known example is the European common lizard (Z. vivipara), which acquired a multiple sex chromosome system (Z1Z2Z1Z2/Z1Z2W) via a sex chromosome–autosome fusion [62]. However, all of these cases represent chromosomal rearrangements, involving whole autosomes or large autosomal fragments. The copying or repositioning of individual genes and their fragments, although predicted by the classical theory of sex chromosome evolution [7], is rarely documented. In the American crow, an intronless fragment of the gene NARF was copied to the W chromosome but was not amplified there [63]. In humans, the DAZ gene was copied to the Y chromosome and amplified in it, similarly to ATF7IP2 in E. velox [64]. Three autosome-derived Y-linked genes were found in dogs [65]. Previously, we detected similar small autosome-derived segments in the sex chromosomes of anole lizards (Ctenonotus sabanus and C. pogus) by the dopseq analysis, but they were not studied in detail [66].
Here, for the first time, we minutely describe an autosome-derived sequence in a reptile sex chromosome, including its content and chromosomal localization. According to FISH, and based on high sequencing coverage, it is highly repetitive and dispersed throughout the whole W chromosome of E. velox, except its centromeric region. In lacertids, the centromeric heterochromatic blocks usually consist of centromere-specific satellites, which are common both to autosomes and sex chromosomes [40,67,68]. Thus, the centromeric region of the E. velox W chromosome is probably constituted of such a satellite, which is not related to the ATF7IP2-related repeat. In the main part of the W chromosome, different microsatellite sequences were previously localized [47]. Therefore, its sequence composition is quite complex, and the ATF7IP2-derived repeat is not the only repetitive sequence hosted in it. It was probably copied from the chromosome 14 and inserted to the W chromosome by a mobile element, and then amplified by subsequent mobile element activity or by intrachromosomal recombination. The absence of detectable homology between the W chromosome of E. velox and the Z chromosomes of reference lacertid genomes may be either because of a complete loss of any gametologous sequences or by the limited sensitivity of our method: since the W chromosome probe was prepared from a single microdissected chromosome, some parts of it may have been lost during microdissection and subsequent PCR amplification.
The reconstructed repeat unit is 10–11 kb long. The existence of contigs that span through the end and the beginning of the reference-based sequence suggests that it could be a tandem repeat, with “tail-to-head” organization. However, the gap inside the reference-based sequence, not covered by any contigs, suggests an alternative explanation: the ATF7IP2-derived sequences may be partially inverted or reshuffled in the repeat unit. Moreover, the different coverage in different parts of ATF7IP2 suggests a divergence between the repeat units in their structure and composition. To reconstruct the exact structure and context of the ATF7IP2-derived repeat in the W chromosome of E. velox, its genomic DNA should be sequenced with a long-read technology, such as PacBio or Nanopore [69].
It is unknown whether the W chromosome copies of ATF7IP2 fragments remain functional. The three preserved exons may potentially be translated, since they do not have stop codons, but it is unclear whether 3 of the 10 exons may still yield a functional product. Given the fact that the W chromosome of E. velox is transcribed, at least at the lampbrush stage, the identified AFT7IP2-derived repeat is the best candidate for the object of this transcriptional activity among the known W chromosome sequences. The proteins ATF7IP and ATF7IP2, also known as MCAF1 and MCAF2, are essential for trimethylation of lysine 9 in histone H3 and thus gene silencing and heterochromatin formation [70,71]. Thereby, if the W chromosome ATF7IP2-derived repeat is translated, its products may participate in chromatin regulation during oogenesis. Recently, non-coding RNAs derived from intronic sequences (sisRNAs, stable intronic sequence RNAa) were identified in the lampbrush-stage nuclei of a frog Xenopus tropicalis, and it was hypothesized that they could play important regulatory roles [72]. This suggests that the ATF7IP2-derived repeat could produce functional transcripts even if it does not encode a protein. Transcriptome and proteome analyses of lampbrush-stage oocytes of E. velox are required to clarify this question.
Insights into the possible functional significance of the ATF7IP2-derived repeat may also be gained from phylogenetic and evolutionary studies. The rapid racerunner is a widely distributed polymorphic species, which began to diversify around 6 MYA [73]. It would be interesting to study the W chromosome composition in different lineages of this species, as well as in other species of Eremias. Among the cytogenetically characterized Eremias species, E. multiocellata has a large heterochromatic W chromosome, similar to that of E. velox [74]. Another species, E. arguta, has a small W chromosome [41]. This suggests that the ATF7IP2-derived repeat may be shared by several species of Eremias. We demonstrated that it is easily found in low-coverage female genomic reads due to its high repetitiveness, and its fragments may be specifically amplified by PCR. This allows screening of the Eremias species for its presence even without cytogenetic analysis, via NGS or qPCR.
In conclusion, our results underline the complex structure and function of highly diverged Y and W chromosomes, which are often erroneously viewed as mere “degenerated” copies of X and Z chromosomes. In addition, our findings are a warning against making long-reaching general hypotheses on sex chromosome origins and evolution based on W and Y chromosome contents: we show that these chromosomes may accumulate sequences from other parts of the genome, unrelated to their origins.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12050779/s1, Table S1: The DOP-SEQ analysis of the sequences of the microdissected W chromosome probe of Eremias velox, aligned to the reference genome of Podarcis muralis.

Author Contributions

Conceptualization, A.L., V.T., and P.B.; methodology, A.L., V.T., and E.S.; software, D.P. and E.S.; validation, A.L. and D.P.; formal analysis, A.L., D.P., and E.S.; investigation, A.L., S.G., A.S., S.R., D.A., and G.D.; resources, V.T. and S.G.; data curation, A.L. and D.P.; writing—original draft preparation, A.L.; writing—review and editing, A.L., P.B., V.T., S.G., S.R., A.S., and E.S.; visualization, A.L.; supervision, A.L.; project administration, A.L., V.T., and P.B.; funding acquisition, V.T. and P.B. All authors have read and agreed to the published version of the manuscript.

Funding

The work was supported by the research grant #19-14-00050 from the Russian Science Foundation, the research grant #19-54-26017 from the Russian Foundation for Basic Research, the research grants #2019-0546 (FSUS-2020-0040) and #0324-2019-0042 from the Ministry of Science and Higher Education (Russia) via the Novosibirsk State University and the Institute of Cytology and Genetics.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Saint Petersburg State University Ethics Committee (statement #131-03-2) and the Institute of Molecular and Cellular Biology Ethics Committee (statement #01/18 from 05.03.2018).

Informed Consent Statement

Not applicable.

Data Availability Statement

The sequencing reads were deposited in the NCBI SRA database under the accession number PRJNA716491. The obtained W chromosome contigs were deposited in GenBank (MW846328-MW846339).

Acknowledgments

We thank the Microscopic Center of the Siberian Branch of the Russian Academy of Sciences and the Chromas Research Resource Center (Saint Petersburg State University) for granting access to microscopic equipment.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Berta, P.; Hawkins, J.R.; Sinclair, A.H.; Taylor, A.; Griffiths, B.L.; Goodfellow, P.N.; Fellous, M. Genetic evidence equating SRY and the testis-determining factor. Nature 1990, 348, 448–450. [Google Scholar] [CrossRef]
  2. Yoshimoto, S.; Okada, E.; Umemoto, H.; Tamura, K.; Uno, Y.; Nishida-Umehara, C.; Matsuda, Y.; Takamatsu, N.; Shiba, T.; Ito, M. A W-linked DM-domain gene, DM-W, participates in primary ovary development in Xenopus laevis. Proc. Natl. Acad. Sci. USA 2008, 105, 2469–2474. [Google Scholar] [CrossRef] [Green Version]
  3. Hirst, C.E.; Major, A.T.; Ayers, K.L.; Brown, R.J.; Mariette, M.; Sackton, T.B.; Smith, C.A. Sex Reversal and Comparative Data Undermine the W Chromosome and Support Z-linked DMRT1 as the Regulator of Gonadal Sex Differentiation in Birds. Endocrinology 2017, 158, 2970–2987. [Google Scholar] [CrossRef]
  4. Bachtrog, D.; Mank, J.E.; Peichel, C.L.; Kirkpatrick, M.; Otto, S.P.; Ashman, T.-L.; Hahn, M.W.; Kitano, J.; Mayrose, I.; Ming, R.; et al. Sex determination: Why so many ways of doing it? PLoS Biol. 2014, 12, e1001899. [Google Scholar] [CrossRef] [Green Version]
  5. Kostmann, A.; Kratochvíl, L.; Rovatsos, M. Poorly differentiated XX/XY sex chromosomes are widely shared across skink radiation. Proc. Biol. Sci. 2021, 288, 20202139. [Google Scholar] [CrossRef]
  6. Rice, W.R. The Accumulation of Sexually Antagonistic Genes as a Selective Agent Promoting the Evolution of Reduced Recombination between Primitive Sex Chromosomes. Evolution 1987, 41, 911. [Google Scholar] [CrossRef]
  7. Charlesworth, D. Evolution of recombination rates between sex chromosomes. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2017, 372. [Google Scholar] [CrossRef] [Green Version]
  8. Waters, P.D.; Ruiz-Herrera, A. Meiotic Executioner Genes Protect the Y from Extinction. Trends Genet. 2020, 36, 728–738. [Google Scholar] [CrossRef]
  9. Wright, A.E.; Darolti, I.; Bloch, N.I.; Oostra, V.; Sandkam, B.; Buechel, S.D.; Kolm, N.; Breden, F.; Vicoso, B.; Mank, J.E. Convergent recombination suppression suggests role of sexual selection in guppy sex chromosome formation. Nat. Commun. 2017, 8, 14251. [Google Scholar] [CrossRef]
  10. Charlesworth, D. The Guppy Sex Chromosome System and the Sexually Antagonistic Polymorphism Hypothesis for Y Chromosome Recombination Suppression. Genes 2018, 9, 264. [Google Scholar] [CrossRef] [Green Version]
  11. Lindholm, A.K.; Brooks, R.; Breden, F. Extreme polymorphism in a Y-linked sexually selected trait. Heredity 2004, 92, 156–162. [Google Scholar] [CrossRef] [Green Version]
  12. Charlesworth, B.; Charlesworth, D. The degeneration of Y chromosomes. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2000, 355, 1563–1572. [Google Scholar] [CrossRef] [Green Version]
  13. Schartl, M.; Schmid, M.; Nanda, I. Dynamics of vertebrate sex chromosome evolution: From equal size to giants and dwarfs. Chromosoma 2016, 125, 553–571. [Google Scholar] [CrossRef]
  14. Matsubara, K.; Nishida, C.; Matsuda, Y.; Kumazawa, Y. Sex chromosome evolution in snakes inferred from divergence patterns of two gametologous genes and chromosome distribution of sex chromosome-linked repetitive sequences. Zool. Lett. 2016, 2, 19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Cortez, D.; Marin, R.; Toledo-Flores, D.; Froidevaux, L.; Liechti, A.; Waters, P.D.; Grützner, F.; Kaessmann, H. Origins and functional evolution of Y chromosomes across mammals. Nature 2014, 508, 488–493. [Google Scholar] [CrossRef] [PubMed]
  16. Bellott, D.W.; Skaletsky, H.; Cho, T.-J.; Brown, L.; Locke, D.; Chen, N.; Galkina, S.; Pyntikova, T.; Koutseva, N.; Graves, T.; et al. Avian W and mammalian Y chromosomes convergently retained dosage-sensitive regulators. Nat. Genet. 2017, 49, 387–394. [Google Scholar] [CrossRef] [Green Version]
  17. Zhou, Q.; Zhang, J.; Bachtrog, D.; An, N.; Huang, Q.; Jarvis, E.D.; Gilbert, M.T.P.; Zhang, G. Complex evolutionary trajectories of sex chromosomes across bird taxa. Science 2014, 346, 1246338. [Google Scholar] [CrossRef] [Green Version]
  18. Gunski, R.J.; Kretschmer, R.; Santos de Souza, M.; de Oliveira Furo, I.; Barcellos, S.A.; Costa, A.L.; Cioffi, M.B.; de Oliveira, E.H.C.; Del Valle Garnero, A. Evolution of Bird Sex Chromosomes Narrated by Repetitive Sequences: Unusual W Chromosome Enlargement in Gallinula melanops (Aves: Gruiformes: Rallidae). Cytogenet. Genome Res. 2019, 158, 152–159. [Google Scholar] [CrossRef] [PubMed]
  19. Augstenová, B.; Mazzoleni, S.; Kratochvíl, L.; Rovatsos, M. Evolutionary Dynamics of the W Chromosome in Caenophidian Snakes. Genes 2017, 9, 5. [Google Scholar] [CrossRef] [Green Version]
  20. Laopichienpong, N.; Muangmai, N.; Chanhome, L.; Suntrarachun, S.; Twilprawat, P.; Peyachoknagul, S.; Srikulnath, K. Evolutionary Dynamics of the Gametologous CTNNB1 Gene on the Z and W Chromosomes of Snakes. J. Hered. 2017, 108, 142–151. [Google Scholar] [CrossRef]
  21. Cechova, M.; Vegesna, R.; Tomaszkiewicz, M.; Harris, R.S.; Chen, D.; Rangavittal, S.; Medvedev, P.; Makova, K.D. Dynamic evolution of great ape Y chromosomes. Proc. Natl. Acad. Sci. USA 2020, 117, 26273–26280. [Google Scholar] [CrossRef]
  22. Xu, L.; Zhou, Q. The Female-Specific W Chromosomes of Birds Have Conserved Gene Contents but Are Not Feminized. Genes 2020, 11, 1126. [Google Scholar] [CrossRef]
  23. Vicoso, B.; Emerson, J.J.; Zektser, Y.; Mahajan, S.; Bachtrog, D. Comparative sex chromosome genomics in snakes: Differentiation, evolutionary strata, and lack of global dosage compensation. PLoS Biol. 2013, 11, e1001643. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Schield, D.R.; Card, D.C.; Hales, N.R.; Perry, B.W.; Pasquesi, G.M.; Blackmon, H.; Adams, R.H.; Corbin, A.B.; Smith, C.F.; Ramesh, B.; et al. The origins and evolution of chromosomes, dosage compensation, and mechanisms underlying venom regulation in snakes. Genome Res. 2019, 29, 590–601. [Google Scholar] [CrossRef]
  25. Kondo, M.; Hornung, U.; Nanda, I.; Imai, S.; Sasaki, T.; Shimizu, A.; Asakawa, S.; Hori, H.; Schmid, M.; Shimizu, N.; et al. Genomic organization of the sex-determining and adjacent regions of the sex chromosomes of medaka. Genome Res. 2006, 16, 815–826. [Google Scholar] [CrossRef] [Green Version]
  26. Jeffries, D.L.; Lavanchy, G.; Sermier, R.; Sredl, M.J.; Miura, I.; Borzée, A.; Barrow, L.N.; Canestrelli, D.; Crochet, P.-A.; Dufresnes, C.; et al. A rapid rate of sex-chromosome turnover and non-random transitions in true frogs. Nat. Commun. 2018, 9, 4088. [Google Scholar] [CrossRef]
  27. Bewick, A.J.; Chain, F.J.J.; Zimmerman, L.B.; Sesay, A.; Gilchrist, M.J.; Owens, N.D.L.; Seifertova, E.; Krylov, V.; Macha, J.; Tlapakova, T.; et al. A large pseudoautosomal region on the sex chromosomes of the frog Silurana tropicalis. Genome Biol. Evol. 2013, 5, 1087–1098. [Google Scholar] [CrossRef] [Green Version]
  28. Stöck, M.; Horn, A.; Grossen, C.; Lindtke, D.; Sermier, R.; Betto-Colliard, C.; Dufresnes, C.; Bonjour, E.; Dumas, Z.; Luquet, E.; et al. Ever-young sex chromosomes in European tree frogs. PLoS Biol. 2011, 9, e1001062. [Google Scholar] [CrossRef] [Green Version]
  29. Tripathi, N.; Hoffmann, M.; Weigel, D.; Dreyer, C. Linkage analysis reveals the independent origin of Poeciliid sex chromosomes and a case of atypical sex inheritance in the guppy (Poecilia reticulata). Genetics 2009, 182, 365–374. [Google Scholar] [CrossRef] [Green Version]
  30. Nanda, I.; Schories, S.; Tripathi, N.; Dreyer, C.; Haaf, T.; Schmid, M.; Schartl, M. Sex chromosome polymorphism in guppies. Chromosoma 2014, 123, 373–383. [Google Scholar] [CrossRef]
  31. Charlesworth, D.; Zhang, Y.; Bergero, R.; Graham, C.; Gardner, J.; Yong, L. Using GC Content to Compare Recombination Patterns on the Sex Chromosomes and Autosomes of the Guppy, Poecilia reticulata, and Its Close Outgroup Species. Mol. Biol. Evol. 2020, 37, 3550–3562. [Google Scholar] [CrossRef] [PubMed]
  32. Lisachov, A.P.; Zadesenets, K.S.; Rubtsov, N.B.; Borodin, P.M. Sex chromosome synapsis and recombination in male guppies. Zebrafish 2015, 12, 174–180. [Google Scholar] [CrossRef] [PubMed]
  33. Marin, R.; Cortez, D.; Lamanna, F.; Pradeepa, M.M.; Leushkin, E.; Julien, P.; Liechti, A.; Halbert, J.; Brüning, T.; Mössinger, K.; et al. Convergent origination of a Drosophila-like dosage compensation mechanism in a reptile lineage. Genome Res. 2017, 27, 1974–1987. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Rovatsos, M.; Rehák, I.; Velenský, P.; Kratochvíl, L. Shared Ancient Sex Chromosomes in Varanids, Beaded Lizards, and Alligator Lizards. Mol. Biol. Evol. 2019, 36, 1113–1120. [Google Scholar] [CrossRef] [PubMed]
  35. Matsubara, K.; Sarre, S.D.; Georges, A.; Matsuda, Y.; Marshall Graves, J.A.; Ezaz, T. Highly differentiated ZW sex microchromosomes in the Australian Varanus species evolved through rapid amplification of repetitive sequences. PLoS ONE 2014, 9, e95226. [Google Scholar] [CrossRef]
  36. Rovatsos, M.; Vukić, J.; Kratochvíl, L. Mammalian X homolog acts as sex chromosome in lacertid lizards. Heredity 2016, 117, 8–13. [Google Scholar] [CrossRef]
  37. Olmo, E.; Odierna, G.; Capriglione, T. Evolution of sex-chromosomes in lacertid lizards. Chromosoma 1987, 96, 33–38. [Google Scholar] [CrossRef] [Green Version]
  38. Rovatsos, M.; Vukić, J.; Mrugała, A.; Suwala, G.; Lymberakis, P.; Kratochvíl, L. Little evidence for switches to environmental sex determination and turnover of sex chromosomes in lacertid lizards. Sci. Rep. 2019, 9, 7832. [Google Scholar] [CrossRef] [Green Version]
  39. Matsubara, K.; Uno, Y.; Srikulnath, K.; Matsuda, Y.; Miller, E.; Olsson, M. No Interstitial Telomeres on Autosomes but Remarkable Amplification of Telomeric Repeats on the W Sex Chromosome in the Sand Lizard (Lacerta agilis). J. Hered. 2015, 106, 753–757. [Google Scholar] [CrossRef] [Green Version]
  40. Giovannotti, M.; Nisi Cerioni, P.; Rojo, V.; Olmo, E.; Slimani, T.; Splendiani, A.; Caputo Barucchi, V. Characterization of a satellite DNA in the genera Lacerta and Timon (Reptilia, Lacertidae) and its role in the differentiation of the W chromosome. J. Exp. Zool. B Mol. Dev. Evol. 2018, 330, 83–95. [Google Scholar] [CrossRef]
  41. Suwala, G.; Altmanová, M.; Mazzoleni, S.; Karameta, E.; Pafilis, P.; Kratochvíl, L.; Rovatsos, M. Evolutionary Variability of W-Linked Repetitive Content in Lacertid Lizards. Genes 2020, 11, 531. [Google Scholar] [CrossRef] [PubMed]
  42. Lisachov, A.P.; Galkina, S.A.; Saifitdinova, A.F.; Svetlana, A.R.; Andreyushkova, D.A.; Trifonov, V.A.; Borodin, P.M. Identification of sex chromosomes in Eremiasvelox (Lacertidae, Reptilia) using lampbrush chromosome analysis. Comp. Cytogenet. 2019, 13, 121–132. [Google Scholar] [CrossRef]
  43. Morgan, G.T. Imaging the dynamics of transcription loops in living chromosomes. Chromosoma 2018, 127, 361–374. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Gall, J.G. Are lampbrush chromosomes unique to meiotic cells? Chromosome Res. 2012, 20, 905–909. [Google Scholar] [CrossRef] [PubMed]
  45. Gaginskaya, E.; Kulikova, T.; Krasikova, A. Avian lampbrush chromosomes: A powerful tool for exploration of genome expression. Cytogenet. Genome Res. 2009, 124, 251–267. [Google Scholar] [CrossRef]
  46. Solovei, I.; Gaginskaya, E.; Hutchison, N.; Macgregor, H. Avian sex chromosomes in the lampbrush form: The ZW lampbrush bivalents from six species of bird. Chromosome Res. 1993, 1, 153–166. [Google Scholar] [CrossRef] [PubMed]
  47. Pokorná, M.; Kratochvíl, L.; Kejnovský, E. Microsatellite distribution on sex chromosomes at different stages of heteromorphism and heterochromatinization in two lizard species (Squamata: Eublepharidae: Coleonyx elegans and lacertidae: Eremias velox). BMC Genet. 2011, 12, 90. [Google Scholar] [CrossRef] [Green Version]
  48. Stanyon, R.; Galleni, L. A rapid fibroblast culture technique for high resolution karyotypes. Bolletino Zool. 1991, 58, 81–83. [Google Scholar] [CrossRef] [Green Version]
  49. Yang, F.; O’Brien, P.C.; Milne, B.S.; Graphodatsky, A.S.; Solanky, N.; Trifonov, V.; Rens, W.; Sargan, D.; Ferguson-Smith, M.A. A complete comparative chromosome map for the dog, red fox, and human and its integration with canine genetic maps. Genomics 1999, 62, 189–202. [Google Scholar] [CrossRef]
  50. Saifitdinova, A.; Galkina, S.; Volodkina, V.; Gaginskaya, E. Preparation of lampbrush chromosomes dissected from avian and reptilian growing oocytes. Biol. Commun. 2017, 62, 165–168. [Google Scholar] [CrossRef]
  51. Sambrook, J.; Fritsch, E.F.; Maniatis, T. Molecular Cloning: A Laboratory Manual, 2nd ed.; Cold Spring Harbor Laboratory Press: New York, NY, USA, 1989; ISBN 0879693096. [Google Scholar]
  52. Li, H.; Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 2009, 25, 1754–1760. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Makunin, A.I.; Kichigin, I.G.; Larkin, D.M.; O’Brien, P.C.M.; Ferguson-Smith, M.A.; Yang, F.; Proskuryakova, A.A.; Vorobieva, N.V.; Chernyaeva, E.N.; O’Brien, S.J.; et al. Contrasting origin of B chromosomes in two cervids (Siberian roe deer and grey brocket deer) unravelled by chromosome-specific DNA sequencing. BMC Genom. 2016, 17, 618. [Google Scholar] [CrossRef] [Green Version]
  54. Flynn, J.M.; Hubley, R.; Goubert, C.; Rosen, J.; Clark, A.G.; Feschotte, C.; Smit, A.F. RepeatModeler2 for automated genomic discovery of transposable element families. Proc. Natl. Acad. Sci. USA 2020, 117, 9451–9457. [Google Scholar] [CrossRef] [PubMed]
  55. Li, W.; Godzik, A. Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 2006, 22, 1658–1659. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Liehr, T. Fluorescence In situ Hybridization (FISH): Application Guide, 2nd ed.; Springer: New York, NY, USA, 2017; ISBN 9783662529591. [Google Scholar]
  57. Waters, P.D.; Duffy, B.; Frost, C.J.; Delbridge, M.L.; Graves, J.A. The human Y chromosome derives largely from a single autosomal region added to the sex chromosomes 80–130 million years ago. Cytogenet. Cell Genet. 2001, 92, 74–79. [Google Scholar] [CrossRef] [PubMed]
  58. Pokorná, M.; Rábová, M.; Ráb, P.; Ferguson-Smith, M.A.; Rens, W.; Kratochvíl, L. Differentiation of sex chromosomes and karyotypic evolution in the eye-lid geckos (Squamata: Gekkota: Eublepharidae), a group with different modes of sex determination. Chromosome Res. 2010, 18, 809–820. [Google Scholar] [CrossRef]
  59. Rovatsos, M.; Johnson Pokorná, M.; Altmanová, M.; Kratochvíl, L. Mixed-Up Sex Chromosomes: Identification of Sex Chromosomes in the X1X1X2X2/X1X2Y System of the Legless Lizards of the Genus Lialis (Squamata: Gekkota: Pygopodidae). Cytogenet. Genome Res. 2016, 149, 282–289. [Google Scholar] [CrossRef]
  60. Giovannotti, M.; Trifonov, V.A.; Paoletti, A.; Kichigin, I.G.; O’Brien, P.C.M.; Kasai, F.; Giovagnoli, G.; Ng, B.L.; Ruggeri, P.; Cerioni, P.N.; et al. New insights into sex chromosome evolution in anole lizards (Reptilia, Dactyloidae). Chromosoma 2017, 126, 245–260. [Google Scholar] [CrossRef]
  61. Lisachov, A.P.; Tishakova, K.V.; Romanenko, S.A.; Molodtseva, A.S.; Prokopov, D.Y.; Pereira, J.C.; Ferguson-Smith, M.A.; Borodin, P.M.; Trifonov, V.A. Whole-Chromosome Fusions in the Karyotype Evolution of Sceloporus (Iguania, Reptilia) are More Intense in Sex Chromosomes Than Autosomes. bioRxiv 2020. [Google Scholar]
  62. Kupriyanova, L.; Kuksin, A.; Odierna, G. Karyotype, chromosome structure, reproductive modalities of three Southern Eurasian populations of the common lacertid lizard, Zootoca vivipara Jacquin, 1787. Acta Herpetol. 2008, 2, 99–106. [Google Scholar]
  63. Xu, L.; Auer, G.; Peona, V.; Suh, A.; Deng, Y.; Feng, S.; Zhang, G.; Blom, M.P.K.; Christidis, L.; Prost, S.; et al. Dynamic evolutionary history and gene content of sex chromosomes across diverse songbirds. Nat. Ecol. Evol. 2019, 3, 834–844. [Google Scholar] [CrossRef] [PubMed]
  64. Saxena, R.; Brown, L.G.; Hawkins, T.; Alagappan, R.K.; Skaletsky, H.; Reeve, M.P.; Reijo, R.; Rozen, S.; Dinulos, M.B.; Disteche, C.M.; et al. The DAZ gene cluster on the human Y chromosome arose from an autosomal gene that was transposed, repeatedly amplified and pruned. Nat. Genet. 1996, 14, 292–299. [Google Scholar] [CrossRef] [PubMed]
  65. Tsai, K.L.; Evans, J.M.; Noorai, R.E.; Starr-Moss, A.N.; Clark, L.A. Novel Y Chromosome Retrocopies in Canids Revealed through a Genome-Wide Association Study for Sex. Genes 2019, 10, 320. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Lisachov, A.P.; Makunin, A.I.; Giovannotti, M.; Pereira, J.C.; Druzhkova, A.S.; Caputo Barucchi, V.; Ferguson-Smith, M.A.; Trifonov, V.A. Genetic Content of the Neo-Sex Chromosomes in Ctenonotus and Norops (Squamata, Dactyloidae) and Degeneration of the Y Chromosome as Revealed by High-Throughput Sequencing of Individual Chromosomes. Cytogenet. Genome Res. 2019, 157, 115–122. [Google Scholar] [CrossRef]
  67. Giovannotti, M.; Rojo, V.; Nisi Cerioni, P.; González-Tizón, A.; Martínez-Lage, A.; Splendiani, A.; Naveira, H.; Ruggeri, P.; Arribas, Ó.; Olmo, E.; et al. Isolation and characterization of two satellite DNAs in some Iberian rock lizards (Squamata, Lacertidae). J. Exp. Zool. B Mol. Dev. Evol. 2014, 322, 13–26. [Google Scholar] [CrossRef] [PubMed]
  68. Rojo, V.; Martínez-Lage, A.; Giovannotti, M.; González-Tizón, A.M.; Nisi Cerioni, P.; Caputo Barucchi, V.; Galán, P.; Olmo, E.; Naveira, H. Evolutionary dynamics of two satellite DNA families in rock lizards of the genus Iberolacerta (Squamata, Lacertidae): Different histories but common traits. Chromosome Res. 2015, 23, 441–461. [Google Scholar] [CrossRef]
  69. Peona, V.; Blom, M.P.K.; Xu, L.; Burri, R.; Sullivan, S.; Bunikis, I.; Liachko, I.; Haryoko, T.; Jønsson, K.A.; Zhou, Q.; et al. Identifying the causes and consequences of assembly gaps using a multiplatform genome assembly of a bird-of-paradise. Mol. Ecol. Resour. 2021, 21, 263–286. [Google Scholar] [CrossRef]
  70. Ichimura, T.; Watanabe, S.; Sakamoto, Y.; Aoto, T.; Fujita, N.; Nakao, M. Transcriptional repression and heterochromatin formation by MBD1 and MCAF/AM family proteins. J. Biol. Chem. 2005, 280, 13928–13935. [Google Scholar] [CrossRef] [Green Version]
  71. Timms, R.T.; Tchasovnikarova, I.A.; Antrobus, R.; Dougan, G.; Lehner, P.J. ATF7IP-Mediated Stabilization of the Histone Methyltransferase SETDB1 Is Essential for Heterochromatin Formation by the HUSH Complex. Cell Rep. 2016, 17, 653–659. [Google Scholar] [CrossRef] [Green Version]
  72. Gardner, E.J.; Nizami, Z.F.; Talbot, C.C.; Gall, J.G. Stable intronic sequence RNA (sisRNA), a new class of noncoding RNA from the oocyte nucleus of Xenopus tropicalis. Genes Dev. 2012, 26, 2550–2559. [Google Scholar] [CrossRef] [Green Version]
  73. Liu, J.; Guo, X.; Chen, D.; Li, J.; Yue, B.; Zeng, X. Diversification and historical demography of the rapid racerunner (Eremias velox) in relation to geological history and Pleistocene climatic oscillations in arid Central Asia. Mol. Phylogenetics Evol. 2019, 130, 244–258. [Google Scholar] [CrossRef] [PubMed]
  74. Wang, C.; Tang, X.; Xin, Y.; Yue, F.; Yan, X.; Liu, B.; An, B.; Wang, X.; Chen, Q. Identification of Sex Chromosomes by Means of Comparative Genomic Hybridization in a Lizard, Eremias multiocellata. Zool. Sci. 2015, 32, 151–156. [Google Scholar] [CrossRef] [PubMed]
Genes 12 00779 g001 550
Figure 1. The distribution of coverage depth for the E. velox sequencing reads along the ATF7IP2 gene of P. muralis. (a) The reads derived from the W-specific DNA sample; (b) the reads derived from the total genomic DNA sample. The Y axes indicate coverage. The X axis indicates the coordinates in the chromosome 14 scaffold of the P. muralis genome assembly (GCF_004329235.1). Exons are indicated as solid bars, introns are indicated by thin lines. The image is obtained via the Integrative Genomics Viewer (https://www.igv.org, accessed on 7 August 2019).
Figure 1. The distribution of coverage depth for the E. velox sequencing reads along the ATF7IP2 gene of P. muralis. (a) The reads derived from the W-specific DNA sample; (b) the reads derived from the total genomic DNA sample. The Y axes indicate coverage. The X axis indicates the coordinates in the chromosome 14 scaffold of the P. muralis genome assembly (GCF_004329235.1). Exons are indicated as solid bars, introns are indicated by thin lines. The image is obtained via the Integrative Genomics Viewer (https://www.igv.org, accessed on 7 August 2019).
Genes 12 00779 g001
Genes 12 00779 g002 550
Figure 2. The reconstructed scheme of the ATF7IP2-derived repeat in the W chromosome of E. velox. (a) Sequence composition of the repeat unit. Pink indicates ATF7IP2-derived sequences, green indicates the POM/Squam1-SINE mobile element, blue indicates other SINE elements, and gray indicates sequences of unknown origin. (b) The reference-based scheme of the ATF7IP2-derived sequences (see Figure 1 for description), with their relative positions within the repeat unit indicated by the red and blue lines.
Figure 2. The reconstructed scheme of the ATF7IP2-derived repeat in the W chromosome of E. velox. (a) Sequence composition of the repeat unit. Pink indicates ATF7IP2-derived sequences, green indicates the POM/Squam1-SINE mobile element, blue indicates other SINE elements, and gray indicates sequences of unknown origin. (b) The reference-based scheme of the ATF7IP2-derived sequences (see Figure 1 for description), with their relative positions within the repeat unit indicated by the red and blue lines.
Genes 12 00779 g002
Genes 12 00779 g003 550
Figure 3. FISH with the PCR-amplified probe to the fragment of the ATF7IP2 gene (red fluorescence) and the W-specific microdissected probe (green fluorescence). Chromosomes are counterstained with DAPI. (a) Metaphase plate of female E. velox. Scale bar: 20 μm. (b) Lampbrush ZW bivalent of E. velox. Scale bar: 15 μm. Insert: agarose gel electrophoresis with the PCR-amplified fragment of the ATF7IP2 gene.
Figure 3. FISH with the PCR-amplified probe to the fragment of the ATF7IP2 gene (red fluorescence) and the W-specific microdissected probe (green fluorescence). Chromosomes are counterstained with DAPI. (a) Metaphase plate of female E. velox. Scale bar: 20 μm. (b) Lampbrush ZW bivalent of E. velox. Scale bar: 15 μm. Insert: agarose gel electrophoresis with the PCR-amplified fragment of the ATF7IP2 gene.
Genes 12 00779 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lisachov, A.; Andreyushkova, D.; Davletshina, G.; Prokopov, D.; Romanenko, S.; Galkina, S.; Saifitdinova, A.; Simonov, E.; Borodin, P.; Trifonov, V. Amplified Fragments of an Autosome-Borne Gene Constitute a Significant Component of the W Sex Chromosome of Eremias velox (Reptilia, Lacertidae). Genes 2021, 12, 779. https://doi.org/10.3390/genes12050779

AMA Style

Lisachov A, Andreyushkova D, Davletshina G, Prokopov D, Romanenko S, Galkina S, Saifitdinova A, Simonov E, Borodin P, Trifonov V. Amplified Fragments of an Autosome-Borne Gene Constitute a Significant Component of the W Sex Chromosome of Eremias velox (Reptilia, Lacertidae). Genes. 2021; 12(5):779. https://doi.org/10.3390/genes12050779

Chicago/Turabian Style

Lisachov, Artem, Daria Andreyushkova, Guzel Davletshina, Dmitry Prokopov, Svetlana Romanenko, Svetlana Galkina, Alsu Saifitdinova, Evgeniy Simonov, Pavel Borodin, and Vladimir Trifonov. 2021. "Amplified Fragments of an Autosome-Borne Gene Constitute a Significant Component of the W Sex Chromosome of Eremias velox (Reptilia, Lacertidae)" Genes 12, no. 5: 779. https://doi.org/10.3390/genes12050779

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop