Nucleic acids movement and its relation to genome dynamics of repetitive DNA: Is cellular and intercellular movement of DNA and RNA molecules related to the evolutionary dynamic genome components?: Is cellular and intercellular movement of DNA and RNA molecules related to the evolutionary dynamic genome components?

There is growing evidence of evolutionary genome plasticity. The evolution of repetitive DNA elements, the major components of most eukaryotic genomes, involves the amplification of various classes of mobile genetic elements, the expansion of satellite DNA, the transfer of fragments or entire organellar genomes and may have connections with viruses. In addition to various repetitive DNA elements, a plethora of large and small RNAs migrate within and between cells during individual development as well as during evolution and contribute to changes of genome structure and function. Such migration of DNA and RNA molecules often results in horizontal gene transfer, thus shaping the whole genomic network of interconnected species. Here, we propose that a high evolutionary dynamism of repetitive genome components is often related to the migration/movement of DNA or RNA molecules. We speculate that the cytoplasm is probably an ideal compartment for such evolutionary experiments.

HiC-TE: a computational pipeline for Hi-C data analysis to study the role of repeat family interactions in the genome 3D organization.

MOTIVATION: The role of repetitive DNA in the 3D organization of the interphase nucleus is a subject of intensive study. In studies of 3D nucleus organization, mutual contacts of various loci can be identified by Hi-C sequencing. Typical analyses use binning of read pairs by location to reduce noise. We use binning by repeat families instead to make similar conclusions about repeat regions. RESULTS: To achieve this, we combined Hi-C data, reference genome data and tools for repeat analysis into a Nextflow pipeline identifying and quantifying the contacts of specific repeat families. As an output, our pipeline produces heatmaps showing contact frequency and circular diagrams visualizing repeat contact localization. Using our pipeline with tomato data, we revealed the preferential homotypic interactions of ribosomal DNA, centromeric satellites and some LTR retrotransposon families and, as expected, little contact between organellar and nuclear DNA elements. While the pipeline can be applied to any eukaryotic genome, results in plants provide better coverage, since the built-in TE-greedy-nester software only detects tandems and LTR retrotransposons. Other repeats can be fed via GFF3 files. This pipeline represents a novel and reproducible way to analyze the role of repetitive elements in the 3D organization of genomes. AVAILABILITY AND IMPLEMENTATION: https://gitlab.fi.muni.cz/lexa/hic-te/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Long-read sequencing technology indicates genome-wide effects of non-B DNA on polymerization speed and error rate.

DNA conformation may deviate from the classical B-form in ∼13% of the human genome. Non-B DNA regulates many cellular processes; however, its effects on DNA polymerization speed and accuracy have not been investigated genome-wide. Such an inquiry is critical for understanding neurological diseases and cancer genome instability. Here, we present the first simultaneous examination of DNA polymerization kinetics and errors in the human genome sequenced with Single-Molecule Real-Time (SMRT) technology. We show that polymerization speed differs between non-B and B-DNA: It decelerates at G-quadruplexes and fluctuates periodically at disease-causing tandem repeats. Analyzing polymerization kinetics profiles, we predict and validate experimentally non-B DNA formation for a novel motif. We demonstrate that several non-B motifs affect sequencing errors (e.g., G-quadruplexes increase error rates), and that sequencing errors are positively associated with polymerase slowdown. Finally, we show that highly divergent G4 motifs have pronounced polymerization slowdown and high sequencing error rates, suggesting similar mechanisms for sequencing errors and germline mutations.

TE-greedy-nester: structure-based detection of LTR retrotransposons and their nesting.

MOTIVATION: Transposable elements (TEs) in eukaryotes often get inserted into one another, forming sequences that become a complex mixture of full-length elements and their fragments. The reconstruction of full-length elements and the order in which they have been inserted is important for genome and transposon evolution studies. However, the accumulation of mutations and genome rearrangements over evolutionary time makes this process error-prone and decreases the efficiency of software aiming to recover all nested full-length TEs. RESULTS: We created software that uses a greedy recursive algorithm to mine increasingly fragmented copies of full-length LTR retrotransposons in assembled genomes and other sequence data. The software called TE-greedy-nester considers not only sequence similarity but also the structure of elements. This new tool was tested on a set of natural and synthetic sequences and its accuracy was compared to similar software. We found TE-greedy-nester to be superior in a number of parameters, namely computation time and full-length TE recovery in highly nested regions. AVAILABILITY AND IMPLEMENTATION: http://gitlab.fi.muni.cz/lexa/nested. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Overexpression of a flower-specific aerolysin-like protein from the dioecious plant Rumex acetosa alters flower development and induces male sterility in transgenic tobacco.

Sex determination in Rumex acetosa, a dioecious plant with a complex XY1 Y2 sex chromosome system (females are XX and males are XY1 Y2 ), is not controlled by an active Y chromosome but depends on the ratio between the number of X chromosomes and autosomes. To gain insight into the molecular mechanisms of sex determination, we generated a subtracted cDNA library enriched in genes specifically or predominantly expressed in female floral buds in early stages of development, when sex determination mechanisms come into play. In the present paper, we report the molecular and functional characterization of FEM32, a gene encoding a protein that shares a common architecture with proteins in different plants, animals, bacteria and fungi of the aerolysin superfamily; many of these function as ß pore-forming toxins. The expression analysis, assessed by northern blot, RT-PCR and in situ hybridization, demonstrates that this gene is specifically expressed in flowers in both early and late stages of development, although its transcripts accumulate much more in female flowers than in male flowers. The ectopic expression of FEM32 under both the constitutive promoter 35S and the flower-specific promoter AP3 in transgenic tobacco showed no obvious alteration in vegetative development but was able to alter floral organ growth and pollen fertility. The 35S::FEM32 and AP3::FEM32 transgenic lines showed a reduction in stamen development and pollen viability, as well as a diminution in fruit set, fruit development and seed production. Compared with other floral organs, pistil development was, however, enhanced in plants overexpressing FEM32. According to these effects, it is likely that FEM32 functions in Rumex by arresting stamen and pollen development during female flower development. The aerolysin-like pore-forming proteins of eukaryotes are mainly involved in defence mechanisms against bacteria, fungi and insects and are also involved in apoptosis and programmed cell death (PCD), a mechanism that could explain the role of FEM32 in Rumex sex determination.

Quadruplex DNA in long terminal repeats in maize LTR retrotransposons inhibits the expression of a reporter gene in yeast.

BACKGROUND: Many studies have shown that guanine-rich DNA sequences form quadruplex structures (G4) in vitro but there is scarce evidence of guanine quadruplexes in vivo. The majority of potential quadruplex-forming sequences (PQS) are located in transposable elements (TEs), especially close to promoters within long terminal repeats of plant LTR retrotransposons. RESULTS: In order to test the potential effect of G4s on retrotransposon expression, we cloned the long terminal repeats of selected maize LTR retrotransposons upstream of the lacZ reporter gene and measured its transcription and translation in yeast. We found that G4s had an inhibitory effect on translation in vivo since "mutants" (where guanines were replaced by adenines in PQS) showed higher expression levels than wild-types. In parallel, we confirmed by circular dichroism measurements that the selected sequences can indeed adopt G4 conformation in vitro. Analysis of RNA-Seq of polyA RNA in maize seedlings grown in the presence of a G4-stabilizing ligand (NMM) showed both inhibitory as well as stimulatory effects on the transcription of LTR retrotransposons. CONCLUSIONS: Our results demonstrate that quadruplex DNA located within long terminal repeats of LTR retrotransposons can be formed in vivo and that it plays a regulatory role in the LTR retrotransposon life-cycle, thus also affecting genome dynamics.

The slowdown of Y chromosome expansion in dioecious Silene latifolia due to DNA loss and male-specific silencing of retrotransposons.

BACKGROUND: The rise and fall of the Y chromosome was demonstrated in animals but plants often possess the large evolutionarily young Y chromosome that is thought has expanded recently. Break-even points dividing expansion and shrinkage phase of plant Y chromosome evolution are still to be determined. To assess the size dynamics of the Y chromosome, we studied intraspecific genome size variation and genome composition of male and female individuals in a dioecious plant Silene latifolia, a well-established model for sex-chromosomes evolution. RESULTS: Our genome size data are the first to demonstrate that regardless of intraspecific genome size variation, Y chromosome has retained its size in S. latifolia. Bioinformatics study of genome composition showed that constancy of Y chromosome size was caused by Y chromosome DNA loss and the female-specific proliferation of recently active dominant retrotransposons. We show that several families of retrotransposons have contributed to genome size variation but not to Y chromosome size change. CONCLUSIONS: Our results suggest that the large Y chromosome of S. latifolia has slowed down or stopped its expansion. Female-specific proliferation of retrotransposons, enlarging the genome with exception of the Y chromosome, was probably caused by silencing of highly active retrotransposons in males and represents an adaptive mechanism to suppress degenerative processes in the haploid stage. Sex specific silencing of transposons might be widespread in plants but hidden in traditional hermaphroditic model plants.

Sex and the flower - developmental aspects of sex chromosome evolution.

Background: The evolution of dioecious plants is occasionally accompanied by the establishment of sex chromosomes: both XY and ZW systems have been found in plants. Structural studies of sex chromosomes are now being followed up by functional studies that are gradually shedding light on the specific genetic and epigenetic processes that shape the development of separate sexes in plants. Scope: This review describes sex determination diversity in plants and the genetic background of dioecy, summarizes recent progress in the investigation of both classical and emerging model dioecious plants and discusses novel findings. The advantages of interspecies hybrids in studies focused on sex determination and the role of epigenetic processes in sexual development are also overviewed. Conclusions: We integrate the genic, genomic and epigenetic levels of sex determination and stress the impact of sex chromosome evolution on structural and functional aspects of plant sexual development. We also discuss the impact of dioecy and sex chromosomes on genome structure and expression.

From Mendel's discovery on pea to today's plant genetics and breeding : Commemorating the 150th anniversary of the reading of Mendel's discovery.

KEY MESSAGE: This work discusses several selected topics of plant genetics and breeding in relation to the 150th anniversary of the seminal work of Gregor Johann Mendel. In 2015, we celebrated the 150th anniversary of the presentation of the seminal work of Gregor Johann Mendel. While Darwin's theory of evolution was based on differential survival and differential reproductive success, Mendel's theory of heredity relies on equality and stability throughout all stages of the life cycle. Darwin's concepts were continuous variation and "soft" heredity; Mendel espoused discontinuous variation and "hard" heredity. Thus, the combination of Mendelian genetics with Darwin's theory of natural selection was the process that resulted in the modern synthesis of evolutionary biology. Although biology, genetics, and genomics have been revolutionized in recent years, modern genetics will forever rely on simple principles founded on pea breeding using seven single gene characters. Purposeful use of mutants to study gene function is one of the essential tools of modern genetics. Today, over 100 plant species genomes have been sequenced. Mapping populations and their use in segregation of molecular markers and marker-trait association to map and isolate genes, were developed on the basis of Mendel's work. Genome-wide or genomic selection is a recent approach for the development of improved breeding lines. The analysis of complex traits has been enhanced by high-throughput phenotyping and developments in statistical and modeling methods for the analysis of phenotypic data. Introgression of novel alleles from landraces and wild relatives widens genetic diversity and improves traits; transgenic methodologies allow for the introduction of novel genes from diverse sources, and gene editing approaches offer possibilities to manipulate gene in a precise manner.

Transposable elements and G-quadruplexes.

A significant part of eukaryotic genomes is formed by transposable elements (TEs) containing not only genes but also regulatory sequences. Some of the regulatory sequences located within TEs can form secondary structures like hairpins or three-stranded (triplex DNA) and four-stranded (quadruplex DNA) conformations. This review focuses on recent evidence showing that G-quadruplex-forming sequences in particular are often present in specific parts of TEs in plants and humans. We discuss the potential role of these structures in the TE life cycle as well as the impact of G-quadruplexes on replication, transcription, translation, chromatin status, and recombination. The aim of this review is to emphasize that TEs may serve as vehicles for the genomic spread of G-quadruplexes. These non-canonical DNA structures and their conformational switches may constitute another regulatory system that, together with small and long non-coding RNA molecules and proteins, contribute to the complex cellular network resulting in the large diversity of eukaryotes.

Impact of repetitive DNA on sex chromosome evolution in plants.

Structurally and functionally diverged sex chromosomes have evolved in many animals as well as in some plants. Sex chromosomes represent a specific genomic region(s) with locally suppressed recombination. As a consequence, repetitive sequences involving transposable elements, tandem repeats (satellites and microsatellites), and organellar DNA accumulate on the Y (W) chromosomes. In this paper, we review the main types of repetitive elements, their gathering on the Y chromosome, and discuss new findings showing that not only accumulation of various repeats in non-recombining regions but also opposite processes form Y chromosome. The aim of this review is also to discuss the mechanisms of repetitive DNA spread involving (retro) transposition, DNA polymerase slippage or unequal crossing-over, as well as modes of repeat removal by ectopic recombination. The intensity of these processes differs in non-recombining region(s) of sex chromosomes when compared to the recombining parts of genome. We also speculate about the relationship between heterochromatinization and the formation of heteromorphic sex chromosomes.

Quadruplex-forming sequences occupy discrete regions inside plant LTR retrotransposons.

Retrotransposons with long terminal repeats (LTR) form a significant proportion of eukaryotic genomes, especially in plants. They have gag and pol genes and several regulatory regions necessary for transcription and reverse transcription. We searched for potential quadruplex-forming sequences (PQSs) and potential triplex-forming sequences (PTSs) in 18 377 full-length LTR retrotransposons collected from 21 plant species. We found that PQSs were often located in LTRs, both upstream and downstream of promoters from which the whole retrotransposon is transcribed. Upstream-located guanine PQSs were dominant in the minus DNA strand, whereas downstream-located guanine PQSs prevailed in the plus strand, indicating their role both at transcriptional and post-transcriptional levels. Our circular dichroism spectroscopy measurements confirmed that these PQSs readily adopted guanine quadruplex structures-some of them were paralell-stranded, while others were anti-parallel-stranded. The PQS often formed doublets at a mutual distance of up to 400 bp. PTSs were most abundant in 3'UTR (but were also present in 5'UTR). We discuss the potential role of quadruplexes and triplexes as the regulators of various processes participating in LTR retrotransposon life cycle and as potential recombination sites during post-insertional retrotransposon-based genome rearrangements.

Fully automated pipeline for detection of sex linked genes using RNA-Seq data.

BACKGROUND: Sex chromosomes present a genomic region which to some extent, differs between the genders of a single species. Reliable high-throughput methods for detection of sex chromosomes specific markers are needed, especially in species where genome information is limited. Next generation sequencing (NGS) opens the door for identification of unique sequences or searching for nucleotide polymorphisms between datasets. A combination of classical genetic segregation analysis along with RNA-Seq data can present an ideal tool to map and identify sex chromosome-specific expressed markers. To address this challenge, we established genetic cross of dioecious plant Rumex acetosa and generated RNA-Seq data from both parental generation and male and female offspring. RESULTS: We present a pipeline for detection of sex linked genes based on nucleotide polymorphism analysis. In our approach, tracking of nucleotide polymorphisms is carried out using a cross of preferably distant populations. For this reason, only 4 datasets are needed - reads from high-throughput sequencing platforms for parent generation (mother and father) and F1 generation (male and female progeny). Our pipeline uses custom scripts together with external assembly, mapping and variant calling software. Given the resource-intensive nature of the computation, servers with high capacity are a requirement. Therefore, in order to keep this pipeline easily accessible and reproducible, we implemented it in Galaxy - an open, web-based platform for data-intensive biomedical research. Our tools are present in the Galaxy Tool Shed, from which they can be installed to any local Galaxy instance. As an output of the pipeline, user gets a FASTA file with candidate transcriptionally active sex-linked genes, sorted by their relevance. At the same time, a BAM file with identified genes and alignment of reads is also provided. Thus, polymorphisms following segregation pattern can be easily visualized, which significantly enhances primer design and subsequent steps of wet-lab verification. CONCLUSIONS: Our pipeline presents a simple and freely accessible software tool for identification of sex chromosome linked genes in species without an existing reference genome. Based on combination of genetic crosses and RNA-Seq data, we have designed a high-throughput, cost-effective approach for a broad community of scientists focused on sex chromosome structure and evolution.

Guanine quadruplexes are formed by specific regions of human transposable elements.

BACKGROUND: Transposable elements form a significant proportion of eukaryotic genomes. Recently, Lexa et al. (Nucleic Acids Res 42:968-978, 2014) reported that plant long terminal repeat (LTR) retrotransposons often contain potential quadruplex sequences (PQSs) in their LTRs and experimentally confirmed their ability to adopt four-stranded DNA conformations. RESULTS: Here, we searched for PQSs in human retrotransposons and found that PQSs are specifically localized in the 3'-UTR of LINE-1 elements, in LTRs of HERV elements and are strongly accumulated in specific regions of SVA elements. Circular dichroism spectroscopy confirmed that most PQSs had adopted monomolecular or bimolecular guanine quadruplex structures. Evolutionarily young SVA elements contained more PQSs than older elements and their propensity to form quadruplex DNA was higher. Full-length L1 elements contained more PQSs than truncated elements; the highest proportion of PQSs was found inside transpositionally active L1 elements (PA2 and HS families). CONCLUSIONS: Conservation of quadruplexes at specific positions of transposable elements implies their importance in their life cycle. The increasing quadruplex presence in evolutionarily young LINE-1 and SVA families makes these elements important contributors toward present genome-wide quadruplex distribution.

Identification of a novel retrotransposon with sex chromosome-specific distribution in Silene latifolia.

Silene latifolia is a dioecious plant species with chromosomal sex determination. Although the evolution of sex chromosomes in S. latifolia has been the subject of numerous studies, a global view of X chromosome structure in this species is still missing. Here, we combine X chromosome microdissection and BAC library screening to isolate new X chromosome-linked sequences. Out of 8 identified BAC clones, only BAC 86M14 showed an X-preferential signal after FISH experiments. Further analysis revealed the existence of the Athila retroelement which is enriched in the X chromosome and nearly absent in the Y chromosome. Based on previous data, the Athila retroelement belongs to the CL3 group of most repetitive sequences in the S. latifolia genome. Structural, transcriptomics and phylogenetic analyses revealed that Athila CL3 represents an old clade in the Athila lineage. We propose a mechanism responsible for Athila CL3 distribution in the S. latifolia genome.

Possible mechanisms responsible for absence of a retrotransposon family on a plant Y chromosome.

Some transposable elements (TEs) show extraordinary variance in abundance along sex chromosomes but the mechanisms responsible for this variance are unknown. Here, we studied Ogre long terminal repeat (LTR) retrotransposons in Silene latifolia, a dioecious plant with evolutionarily young heteromorphic sex chromosomes. Ogre elements are ubiquitous in the S. latifolia genome but surprisingly absent on the Y chromosome. Bacterial artificial chromosome (BAC) library analysis and fluorescence in situ hybridization (FISH) were used to determine Ogre structure and chromosomal localization. Next generation sequencing (NGS) data were analysed to assess the transcription level and abundance of small RNAs. Methylation of Ogres was determined by bisulphite sequencing. Phylogenetic analysis was used to determine mobilization time and selection forces acting on Ogre elements. We characterized three Ogre families ubiquitous in the S. latifolia genome. One family is nearly absent on the Y chromosome despite all the families having similar structures and spreading mechanisms. We showed that Ogre retrotransposons evolved before sex chromosomes appeared but were mobilized after formation of the Y chromosome. Our data suggest that the absence of one Ogre family on the Y chromosome may be caused by 24-nucleotide (24-nt) small RNA-mediated silencing leading to female-specific spreading. Our findings highlight epigenetic silencing mechanisms as potentially crucial factors in sex-specific spreading of some TEs, but other possible mechanisms are also discussed.

Telomeric retrotransposons show propensity to form G-quadruplexes in various eukaryotic species.

BACKGROUND: Canonical telomeres (telomerase-synthetised) are readily forming G-quadruplexes (G4) on the G-rich strand. However, there are examples of non-canonical telomeres among eukaryotes where telomeric tandem repeats are invaded by specific retrotransposons. Drosophila melanogaster represents an extreme example with telomeres composed solely by three retrotransposons-Het-A, TAHRE and TART (HTT). Even though non-canonical telomeres often show strand biased G-distribution, the evidence for the G4-forming potential is limited. RESULTS: Using circular dichroism spectroscopy and UV absorption melting assay we have verified in vitro G4-formation in the HTT elements of D. melanogaster. Namely 3 in Het-A, 8 in TART and 2 in TAHRE. All the G4s are asymmetrically distributed as in canonical telomeres. Bioinformatic analysis showed that asymmetric distribution of potential quadruplex sequences (PQS) is common in telomeric retrotransposons in other Drosophila species. Most of the PQS are located in the gag gene where PQS density correlates with higher DNA sequence conservation and codon selection favoring G4-forming potential. The importance of G4s in non-canonical telomeres is further supported by analysis of telomere-associated retrotransposons from various eukaryotic species including green algae, Diplomonadida, fungi, insects and vertebrates. Virtually all analyzed telomere-associated retrotransposons contained PQS, frequently with asymmetric strand distribution. Comparison with non-telomeric elements showed independent selection of PQS-rich elements from four distinct LINE clades. CONCLUSION: Our findings of strand-biased G4-forming motifs in telomere-associated retrotransposons from various eukaryotic species support the G4-formation as one of the prerequisites for the recruitment of specific retrotransposons to chromosome ends and call for further experimental studies.

Variation in G-quadruplex sequence and topology differentially impacts human DNA polymerase fidelity.

G-quadruplexes (G4s), a type of non-B DNA, play important roles in a wide range of molecular processes, including replication, transcription, and translation. Genome integrity relies on efficient and accurate DNA synthesis, and is compromised by various stressors, to which non-B DNA structures such as G4s can be particularly vulnerable. However, the impact of G4 structures on DNA polymerase fidelity is largely unknown. Using an in vitro forward mutation assay, we investigated the fidelity of human DNA polymerases delta (δ4, four-subunit), eta (η), and kappa (κ) during synthesis of G4 motifs representing those in the human genome. The motifs differ in sequence, topology, and stability, features that may affect DNA polymerase errors. Polymerase error rate hierarchy (δ4 < κ < Î·) is largely maintained during G4 synthesis. Importantly, we observed unique polymerase error signatures during synthesis of VEGF G4 motifs, stable G4s which form parallel topologies. These statistically significant errors occurred within, immediately flanking, and encompassing the G4 motif. For pol δ4, the errors were deletions, insertions and complex errors within the G4 or encompassing the G4 motif and surrounding sequence. For pol η, the errors occurred in 3' sequences flanking the G4 motif. For pol κ, the errors were frameshift mutations within G-tracts of the G4. Because these error signatures were not observed during synthesis of an antiparallel G4 and, to a lesser extent, a hybrid G4, we suggest that G4 topology and/or stability could influence polymerase fidelity. Using in silico analyses, we show that most polymerase errors are predicted to have minimal effects on predicted G4 stability. Our results provide a unique view of G4s not previously elucidated, showing that G4 motif heterogeneity differentially influences polymerase fidelity within the motif and flanking sequences. Thus, our study advances the understanding of how DNA polymerase errors contribute to G4 mutagenesis.

Evidence for degeneration of the Y chromosome in the dioecious plant Silene latifolia.

The human Y--probably because of its nonrecombining nature--has lost 97% of its genes since X and Y chromosomes started to diverge [1, 2]. There are clear signs of degeneration in the Drosophila miranda neoY chromosome (an autosome fused to the Y chromosome), with neoY genes showing faster protein evolution [3-6], accumulation of unpreferred codons [6], more insertions of transposable elements [5, 7], and lower levels of expression [8] than neoX genes. In the many other taxa with sex chromosomes, Y degeneration has hardly been studied. In plants, many genes are expressed in pollen [9], and strong pollen selection may oppose the degeneration of plant Y chromosomes [10]. Silene latifolia is a dioecious plant with young heteromorphic sex chromosomes [11, 12]. Here we test whether the S. latifolia Y chromosome is undergoing genetic degeneration by analyzing seven sex-linked genes. S. latifolia Y-linked genes tend to evolve faster at the protein level than their X-linked homologs, and they have lower expression levels. Several Y gene introns have increased in length, with evidence for transposable-element accumulation. We detect signs of degeneration in most of the Y-linked gene sequences analyzed, similar to those of animal Y-linked and neo-Y chromosome genes.

Microsatellite distribution on sex chromosomes at different stages of heteromorphism and heterochromatinization in two lizard species (Squamata: Eublepharidae: Coleonyx elegans and lacertidae: Eremias velox).

BACKGROUND: The accumulation of repetitive sequences such as microsatellites during the differentiation of sex chromosomes has not been studied in most squamate reptiles (lizards, amphisbaenians and snakes), a group which has a large diversity of sex determining systems. It is known that the Bkm repeats containing tandem arrays of GATA tetranucleotides are highly accumulated on the degenerated W chromosomes in advanced snakes. Similar, potentially homologous, repetitive sequences were found on sex chromosomes in other vertebrates. Using FISH with probes containing all possible mono-, di-, and tri-nucleotide sequences and GATA, we studied the genome distribution of microsatellite repeats on sex chromosomes in two lizard species (the gecko Coleonyx elegans and the lacertid Eremias velox) with independently evolved sex chromosomes. The gecko possesses heteromorphic euchromatic sex chromosomes, while sex chromosomes in the lacertid are homomorphic and the W chromosome is highly heterochromatic. Our aim was to test whether microsatellite distribution on sex chromosomes corresponds to the stage of their heteromorphism or heterochromatinization. Moreover, because the lizards lie phylogenetically between snakes and other vertebrates with the Bkm-related repeats on sex chromosomes, the knowledge of their repetitive sequence is informative for the determination of conserved versus convergently evolved repetitive sequences across vertebrate lineages. RESULTS: Heteromorphic sex chromosomes of C. elegans do not show any sign of microsatellite accumulation. On the other hand, in E. velox, certain microsatellite sequences are extensively accumulated over the whole length or parts of the W chromosome, while others, including GATA, are absent on this heterochromatinized sex chromosome. CONCLUSION: The accumulation of microsatellite repeats corresponds to the stage of heterochromatinization of sex chromosomes rather than to their heteromorphism. The lack of GATA repeats on the sex chromosomes of both lizards suggests that the Bkm-related repeats on sex chromosomes in snakes and other vertebrates evolved convergently. The comparison of microsatellite sequences accumulated on sex chromosomes in E. velox and in other eukaryotic organisms suggests that historical contingency, not characteristics of particular sequences, plays a major role in the determination of which microsatellite sequence is accumulated on the sex chromosomes in a particular lineage.