The enigma of Y chromosome degeneration: TRAM, a novel retrotransposon is preferentially located on the Neo-Y chromosome of Drosophila miranda

Establishment of H3K9me3-dependent heterochromatin during embryogenesis in Drosophila miranda

Heterochromatin is a key architectural feature of eukaryotic genomes crucial for silencing of repetitive elements. During Drosophila embryonic cellularization, heterochromatin rapidly appears over repetitive sequences, but the molecular details of how heterochromatin is established are poorly understood. Here, we map the genome-wide distribution of H3K9me3-dependent heterochromatin in individual embryos of Drosophila miranda at precisely staged developmental time points. We find that canonical H3K9me3 enrichment is established prior to cellularization and matures into stable and broad heterochromatin domains through development. Intriguingly, initial nucleation sites of H3K9me3 enrichment appear as early as embryonic stage 3 over transposable elements (TEs) and progressively broaden, consistent with spreading to neighboring nucleosomes. The earliest nucleation sites are limited to specific regions of a small number of recently active retrotransposon families and often appear over promoter and 5' regions of LTR retrotransposons, while late nucleation sites develop broadly across the entirety of most TEs. Interestingly, early nucleating TEs are strongly associated with abundant maternal piRNAs and show early zygotic transcription. These results support a model of piRNA-associated co-transcriptional silencing while also suggesting additional mechanisms for site-restricted H3K9me3 nucleation at TEs in pre-cellular Drosophila embryos.

Sex chromosome repeats tip the balance towards speciation

Because sex chromosomes, by definition, carry genes that determine sex, mutations that alter their structural and functional stability can have immediate consequences for the individual by reducing fertility, but also for a species by altering the sex ratio. Moreover, the sex-specific segregation patterns of heteromorphic sex chromosomes make them havens for selfish genetic elements that not only create suboptimal sex ratios but can also foster sexual antagonism. Compensatory mutations to mitigate antagonism or return sex ratios to a Fisherian optimum can create hybrid incompatibility and establish reproductive barriers leading to species divergence. The destabilizing influence of these selfish elements is often manifest within populations as copy number variants (CNVs) in satellite repeats and transposable elements (TE) or as CNVs involving sex-determining genes, or genes essential to fertility and sex chromosome dosage compensation. This review catalogs several examples of well-studied sex chromosome CNVs in Drosophilids and mammals that underlie instances of meiotic drive, hybrid incompatibility and disruptions to sex differentiation and sex chromosome dosage compensation. While it is difficult to pinpoint a direct cause/effect relationship between these sex chromosome CNVs and speciation, it is easy to see how their effects in creating imbalances between the sexes, and the compensatory mutations to restore balance, can lead to lineage splitting and species formation.

Uncovering the evolutionary history of neo-XY sex chromosomes in the grasshopper Ronderosia bergii (Orthoptera, Melanoplinae) through satellite DNA analysis

BACKGROUND

Neo-sex chromosome systems arose independently multiple times in evolution, presenting the remarkable characteristic of repetitive DNAs accumulation. Among grasshoppers, occurrence of neo-XY was repeatedly noticed in Melanoplinae. Here we analyzed the most abundant tandem repeats of R. bergii (2n = 22, neo-XY♂) using deep Illumina sequencing and graph-based clustering in order to address the neo-sex chromosomes evolution.

RESULTS

The analyses revealed ten families of satDNAs comprising about ~1% of the male genome, which occupied mainly C-positive regions of autosomes. Regarding the sex chromosomes, satDNAs were recorded within centromeric or interstitial regions of the neo-X chromosome and four satDNAs occurred in the neo-Y, two of them being exclusive (Rber248 and Rber299). Using a combination of probes we uncovered five well-defined cytological variants for neo-Y, originated by multiple paracentric inversions and satDNA amplification, besides fragmented neo-Y. These neo-Y variants were distinct in frequency between embryos and adult males.

CONCLUSIONS

The genomic data together with cytogenetic mapping enabled us to better understand the neo-sex chromosome dynamics in grasshoppers, reinforcing differentiation of neo-X and neo-Y and revealing the occurrence of multiple additional rearrangements involved in the neo-Y evolution of R. bergii. We discussed the possible causes that led to differences in frequency for the neo-Y variants between embryos and adults. Finally we hypothesize about the role of DNA satellites in R. bergii as well as putative historical events involved in the evolution of the R. bergii neo-XY.

Satellite DNA: An Evolving Topic

Satellite DNA represents one of the most fascinating parts of the repetitive fraction of the eukaryotic genome. Since the discovery of highly repetitive tandem DNA in the 1960s, a lot of literature has extensively covered various topics related to the structure, organization, function, and evolution of such sequences. Today, with the advent of genomic tools, the study of satellite DNA has regained a great interest. Thus, Next-Generation Sequencing (NGS), together with high-throughput in silico analysis of the information contained in NGS reads, has revolutionized the analysis of the repetitive fraction of the eukaryotic genomes. The whole of the historical and current approaches to the topic gives us a broad view of the function and evolution of satellite DNA and its role in chromosomal evolution. Currently, we have extensive information on the molecular, chromosomal, biological, and population factors that affect the evolutionary fate of satellite DNA, knowledge that gives rise to a series of hypotheses that get on well with each other about the origin, spreading, and evolution of satellite DNA. In this paper, I review these hypotheses from a methodological, conceptual, and historical perspective and frame them in the context of chromosomal organization and evolution.

High-throughput analysis of the satellitome revealed enormous diversity of satellite DNAs in the neo-Y chromosome of the cricket Eneoptera surinamensis

Satellite DNAs (satDNAs) constitute large portion of eukaryote genomes, comprising non-protein-coding sequences tandemly repeated. They are mostly found in heterochromatic regions of chromosomes such as around centromere or near telomeres, in intercalary heterochromatin, and often in non-recombining segments of sex chromosomes. We examined the satellitome in the cricket Eneoptera surinamensis (2n = 9, neo-X1X2Y, males) to characterize the molecular evolution of its neo-sex chromosomes. To achieve this, we analyzed illumina reads using graph-based clustering and complementary analyses. We found an unusually high number of 45 families of satDNAs, ranging from 4 bp to 517 bp, accounting for about 14% of the genome and showing different modular structures and high diversity of arrays. FISH mapping revealed that satDNAs are located mostly in C-positive pericentromeric regions of the chromosomes. SatDNAs enrichment was also observed in the neo-sex chromosomes in comparison to autosomes. Especially astonishing accumulation of satDNAs loci was found in the highly differentiated neo-Y, including 39 satDNAs over-represented in this chromosome, which is the greatest satDNAs diversity yet reported for sex chromosomes. Our results suggest possible involvement of satDNAs in genome increasing and in molecular differentiation of the neo-sex chromosomes in this species, contributing to the understanding of sex chromosome composition and evolution in Orthoptera.

Contrasting the Chromosomal Organization of Repetitive DNAs in Two Gryllidae Crickets with Highly Divergent Karyotypes

A large percentage of eukaryotic genomes consist of repetitive DNA that plays an important role in the organization, size and evolution. In the case of crickets, chromosomal variability has been found using classical cytogenetics, but almost no information concerning the organization of their repetitive DNAs is available. To better understand the chromosomal organization and diversification of repetitive DNAs in crickets, we studied the chromosomes of two Gryllidae species with highly divergent karyotypes, i.e., 2n(♂) = 29,X0 (Gryllus assimilis) and 2n = 9, neo-X1X2Y (Eneoptera surinamensis). The analyses were performed using classical cytogenetic techniques, repetitive DNA mapping and genome-size estimation. Conserved characteristics were observed, such as the occurrence of a small number of clusters of rDNAs and U snDNAs, in contrast to the multiple clusters/dispersal of the H3 histone genes. The positions of U2 snDNA and 18S rDNA are also conserved, being intermingled within the largest autosome. The distribution and base-pair composition of the heterochromatin and repetitive DNA pools of these organisms differed, suggesting reorganization. Although the microsatellite arrays had a similar distribution pattern, being dispersed along entire chromosomes, as has been observed in some grasshopper species, a band-like pattern was also observed in the E. surinamensis chromosomes, putatively due to their amplification and clustering. In addition to these differences, the genome of E. surinamensis is approximately 2.5 times larger than that of G. assimilis, which we hypothesize is due to the amplification of repetitive DNAs. Finally, we discuss the possible involvement of repetitive DNAs in the differentiation of the neo-sex chromosomes of E. surinamensis, as has been reported in other eukaryotic groups. This study provided an opportunity to explore the evolutionary dynamics of repetitive DNAs in two non-model species and will contribute to the understanding of chromosomal evolution in a group about which little chromosomal and genomic information is known.

Achilles, a New Family of Transcriptionally Active Retrotransposons from the Olive Fruit Fly, with Y Chromosome Preferential Distribution

Sex chromosomes have many unusual features relative to autosomes. The in depth exploration of their structure will improve our understanding of their origin and divergence (degeneration) as well as the evolution of genetic sex determination pathways which, most often are attributed to them. In Tephritids, the structure of Y chromosome, where the male-determining factor M is localized, is largely unexplored and limited data concerning its sequence content and evolution are available. In order to get insight into the structure and organization of the Y chromosome of the major olive insect pest, the olive fly Bactrocera oleae, we characterized sequences from a Pulse Field Gel Electrophoresis (PFGE)-isolated Y chromosome. Here, we report the discovery of the first olive fly LTR retrotransposon with increased presence on the Y chromosome. The element belongs to the BEL-Pao superfamily, however, its sequence comparison with the other members of the superfamily suggests that it constitutes a new family that we termed Achilles. Its ~7.5 kb sequence consists of the 5'LTR, the 5'non-coding sequence and the open reading frame (ORF), which encodes the polyprotein Gag-Pol. In situ hybridization to the B. oleae polytene chromosomes showed that Achilles is distributed in discrete bands dispersed on all five autosomes, in all centromeric regions and in the granular heterochromatic network corresponding to the mitotic sex chromosomes. The between sexes comparison revealed a variation in Achilles copy number, with male flies possessing 5-10 copies more than female (CI range: 18-38 and 12-33 copies respectively per genome). The examination of its transcriptional activity demonstrated the presence of at least one intact active copy in the genome, showing a differential level of expression between sexes as well as during embryonic development. The higher expression was detected in male germline tissues (testes). Moreover, the presence of Achilles-like elements in different species of the Tephritidae family suggests an ancient origin of this element.

Tracking the evolution of sex chromosome systems in Melanoplinae grasshoppers through chromosomal mapping of repetitive DNA sequences

BACKGROUND

The accumulation of repetitive DNA during sex chromosome differentiation is a common feature of many eukaryotes and becomes more evident after recombination has been restricted or abolished. The accumulated repetitive sequences include multigene families, microsatellites, satellite DNAs and mobile elements, all of which are important for the structural remodeling of heterochromatin. In grasshoppers, derived sex chromosome systems, such as neo-XY♂/XX♀ and neo-X1X2Y♂/X1X1X2X2♀, are frequently observed in the Melanoplinae subfamily. However, no studies concerning the evolution of sex chromosomes in Melanoplinae have addressed the role of the repetitive DNA sequences. To further investigate the evolution of sex chromosomes in grasshoppers, we used classical cytogenetic and FISH analyses to examine the repetitive DNA sequences in six phylogenetically related Melanoplinae species with X0♂/XX♀, neo-XY♂/XX♀ and neo-X1X2Y♂/X1X1X2X2♀ sex chromosome systems.

RESULTS

Our data indicate a non-spreading of heterochromatic blocks and pool of repetitive DNAs (C0t-1 DNA) in the sex chromosomes; however, the spreading of multigene families among the neo-sex chromosomes of Eurotettix and Dichromatos was remarkable, particularly for 5S rDNA. In autosomes, FISH mapping of multigene families revealed distinct patterns of chromosomal organization at the intra- and intergenomic levels.

CONCLUSIONS

These results suggest a common origin and subsequent differential accumulation of repetitive DNAs in the sex chromosomes of Dichromatos and an independent origin of the sex chromosomes of the neo-XY and neo-X1X2Y systems. Our data indicate a possible role for repetitive DNAs in the diversification of sex chromosome systems in grasshoppers.

BEL/Pao retrotransposons in metazoan genomes

BACKGROUND

Long terminal repeat (LTR) retrotransposons are a widespread kind of transposable element present in eukaryotic genomes. They are a major factor in genome evolution due to their ability to create large scale mutations and genome rearrangements. Compared to other transposable elements, little attention has been paid to elements belonging to the metazoan BEL/Pao subclass of LTR retrotransposons. No comprehensive characterization of these elements is available so far. The aim of this study was to describe all BEL/Pao elements in a set of 62 sequenced metazoan genomes, and to analyze their phylogenetic relationship.

RESULTS

We identified a total of 7,861 BEL/Pao elements in 53 of our 62 study genomes. We identified BEL/Pao elements in 20 genomes where such elements had not been found so far. Our analysis shows that BEL/Pao elements are the second-most abundant class of LTR retrotransposons in the genomes we study, more abundant than Ty1/Copia elements, and second only to Ty3/Gypsy elements. They occur in multiple phyla, including basal metazoan phyla, suggesting that BEL/Pao elements arose early in animal evolution. We confirm findings from previous studies that BEL/Pao elements do not occur in mammals. The elements we found can be grouped into more than 1725 families, 1623 of which are new, previously unknown families. These families fall into seven superfamilies, only five of which have been characterized so far. One new superfamily is a major subdivision of the Pao superfamily which we propose to call Dan, because it is restricted to the genome of the zebrafish Danio rerio. The other new superfamily comprises 83 elements and is restricted to lower aquatic eumetazoans. We propose to call this superfamily Flow. BEL/Pao elements do not show any signs of recent horizontal gene transfer between distantly related species.

CONCLUSIONS

In sum, our analysis identifies thousands of new BEL/Pao elements and provides new insights into their distribution, abundance, and evolution.

Applying mobile genetic elements for genome analysis and evolution

Transposable elements (TEs) are ubiquitous components of all living organisms, and in the course of their coexistence with their respective host genomes, these parasitc DNAs have played important roles in the evolution of complex genetic networks. The interaction between mobile DNAs and their host genomes are quite diverse, ranging from modifications of gene structure and regulation to alterations in general genome architecture. Thus during evolutionary time these elements can be regarded as natural molecular tools in shaping the organization, structure, and function of eukaryotic genes and genomes. Based on their intrinsic properties and features, mobile DNAs are widely applied at present as a technical "toolbox," essential for studying a diverse spectrum of biological questions. In this review, we aim to summarize both the evolutionary impact of TEs on genome evolution and their valuable and diverse methodological applications as molecular tools.

Molecular cytogenetic characterization of Rumex papillaris, a dioecious plant with an XX/XY(1)Y (2) sex chromosome system

Rumex papillaris Boiss, & Reut., an Iberian endemic, belongs to the section Acetosa of the genus Rumex whose main representative is R. acetosa L., a species intensively studied in relation to sex-chromosome evolution. Here, we characterize cytogenetically the chromosomal complement of R. papillaris in an effort to enhance future comparative genomic approaches and to better our understanding of sex chromosome structure in plants. Rumex papillaris, as is common in this group, is a dioecious species characterized by the presence of a multiple sex chromosome system (with females 2n = 12 + XX and males 2n = 12 + XY(1)Y(2)). Except for the X chromosome both Y chromosomes are the longest in the karyotype and appear heterochromatic due to the accumulation of at least two satellite DNA families, RAE180 and RAYSI. Each chromosome of pair VI has an additional major heterochromatin block at the distal region of the short arm. These supernumerary heterochromatic blocks are occupied by RAE730 satellite DNA family. The Y-related RAE180 family is also present in an additional minor autosomal locus. Our comparative study of the chromosomal organization of the different satellite-DNA sequences in XX/XY and XX/XY(1)Y(2) Rumex species demonstrates that of active mechanisms of heterochromatin amplification occurred and were accompanied by chromosomal rearrangements giving rise to the multiple XX/XY(1)Y(2) chromosome systems observed in Rumex. Additionally, Y(1) and Y(2) chromosomes have undergone further rearrangements leading to differential patterns of Y-heterochromatin distribution between Rumex species with multiple sex chromosome systems.

Abundance and chromosomal distribution of six Drosophila buzzatii transposons: BuT1, BuT2, BuT3, BuT4, BuT5, and BuT6

The abundance and chromosomal distribution of six class-II transposable elements (TEs) of Drosophila buzzatii have been analyzed by Southern blotting and in situ hybridization. These six transposons had been previously found at the breakpoints of inversions 2j and 2q ( 7 ) of D. buzzatii. These two polymorphic inversions were generated by an ectopic recombination event between two copies of Galileo, a Foldback element. The four breakpoints became hotspots for TE insertions after the generation of the inversion and the transposons analyzed in this work are considered to be secondary invaders of these regions. Insertions of the six transposons are present in the euchromatin but show an increased density in the pericentromeric euchromatin-heterochromatin transition region and the dot chromosome. They are also more abundant in the inverted segments of chromosome 2 rearrangements. We further observed that the accumulation of TE insertions varies between elements and is correlated between dot, proximal regions, and inverted segments. These observations fully agree with previous data in Drosophila melanogaster and support recombination rate as the chief force explaining the chromosomal distribution of TEs.

The Y chromosome as a target for acquired and amplified genetic material in evolution

The special properties of the Y chromosome stem form the fact that it is a non-recombining degenerate derivative of the X chromosome. The absence of homologous recombination between the X and the Y chromosome leads to gradual degeneration of various Y chromosome genes on an evolutionary timescale. The absence of recombination, however, also favors the accumulation of transposable elements on the Y chromosome during its evolution, as seen with both Drosophila and mammalian Y chromosomes. Alongside these processes, the acquisition and amplification of autosomal male benefit genes occur. This review will focus on recent studies that reveal the autosome-acquired genes on the Y chromosome of both Drosophila and humans. The evolution of the acquired and amplified genes on the Y chromosome is also discussed. Molecular and comparative analyses of Y-linked repeats in the Drosophila melanogaster genome demonstrate that there was a period of their degeneration followed by a period of their integration into RNAi silencing, which was beneficial for male fertility. Finally, the function of non-coding RNA produced by amplified Y chromosome genetic elements will be discussed.

The origin and evolution of the variability in a Y-specific satellite-DNA of Rumex acetosa and its relatives

In this paper, we analyze a satellite-DNA family, the RAYSI family, which is specific of the Y chromosomes of Rumex acetosa, a dioecious plant species with a multiple sex-chromosome system in which the females are XX and the males are XY(1)Y(2). Here, we demonstrate that this satellite DNA is common to other relatives of R. acetosa, including Rumex papillaris, Rumex intermedius, Rumex thyrsoides and Rumex tuberosus that are also dioecious species with a multiple system of sex chromosomes. This satellite-DNA family is absent from the genomes of other dioecious Rumex species having an XX/XY sex-chromosome system. Our data confirm recent molecular phylogenies that support a unique origin for all dioecious species of Rumex and two separate lineages for species with single or complex sex-chromosome systems. Our data also support an accelerated degeneration of Y-chromosome in XX/XY(1)Y(2) species by the accumulation of satellite-DNA sequences. On the other hand, the particular non-recombining nature of the Y chromosomes of R. acetosa and their closest relatives lead to a particular mode of evolution of RAYSI sequences. Thus, mechanisms leading to the suppression of recombination between the Y chromosomes reduced the rate of concerted evolution and gave rise to the apparition of different RAYSI subfamilies. Thus, R. acetosa and R. intermedius have two subfamilies (the RAYSI-S and RAYSI-J subfamilies and the INT-A and INT-B subfamilies, respectively), while R. papillaris only has one, the RAYSI-J subfamily. The RAYSI-S and RAYSI-J subfamilies of R. acetosa differ in 83 fixed diagnostic sites and several diagnostic deletions while the INT-A and the INT-B of R. intermedius differ in 27 fixed diagnostic sites. Pairwise comparisons between RAYSI-S and RAYSI-J sequences or between INT-A and INT-B sequences revealed these sites to be shared mutations detectable in repeats of the same variant in same positions. Evolutionary comparisons suggest that the subfamily RAYSI-J has appeared in the common ancestor of R. acetosa and R. papillaris, in which RAYSI-J has replaced totally (R. papillaris) or almost totally the ancestral sequence (R. acetosa). This scenario assumes that RAYSI-S sequences should be considered ancestral sequences and that a secondary event of subfamily subdivision should be occurring in R. intermedius, with their RAYSI subfamilies more closely related to one another than with other RAYSI sequences. Our analysis suggests that the different subfamilies diverged by a gradual and cohesive way probably mediated by sister-chromatid interchanges while their expansion or contraction in number might be explained by alternating cycles of sudden mechanisms of amplification or elimination.

Y chromosomes: born to be destroyed

Suppression of recombination is the prerequisite for stable genetically determined sex systems. A consequence of suppression of recombination is the strong bias in the distribution of transposable elements (TEs), mostly retrotransposons. Our results and those from others indicate that the major force driving the degeneration of Y chromosomes are retrotransposons in remodelling former euchromatic chromosome structures into heterochromatic ones. We put forward the following hypotheses. (1) A massive accumulation of retrotransposons occurs early in non-recombining regions. (2) Heterochromatic nucleation centres are formed as a genomic defence mechanism against invasive parasitic elements. The newly established nucleation centres become epigenetically inherited. (3) Spreading of heterochromatin from the nucleation centres into flanking regions induces, in an adaptive process, transcriptional gene silencing of neighbourhood genes that could either be still intact or in an already eroded condition. (4) Constitutive silenced genes are not under the same genetic selection pressure as active genes. They are more exposed to the decay process. (5) Gene dosage balance is re-established by the parallel evolution of dosage compensation mechanisms.

Partial deletions of the W chromosome due to reciprocal translocation in the silkworm Bombyx mori

In the silkworm, Bombyx mori (female, ZW; male, ZZ), femaleness is determined by the presence of a single W chromosome, irrespective of the number of autosomes or Z chromosomes. The W chromosome is devoid of functional genes, except the putative female-determining gene (Fem). However, there are strains in which chromosomal fragments containing autosomal markers have been translocated on to W. In this study, we analysed the W chromosomal regions of the Zebra-W strain (T(W;3)Ze chromosome) and the Black-egg-W strain (T(W;10)+(w-2) chromosome) at the molecular level. Initially, we undertook a project to identify W-specific RAPD markers, in addition to the three already established W-specific RAPD markers (W-Kabuki, W-Samurai and W-Kamikaze). Following the screening of 3648 arbitrary 10-mer primers, we obtained nine W-specific RAPD marker sequences (W-Bonsai, W-Mikan, W-Musashi, W-Rikishi, W-Sakura, W-Sasuke, W-Yukemuri-L, W-Yukemuri-S and BMC1-Kabuki), almost all of which contained the border regions of retrotransposons, namely portions of nested retrotransposons. We confirmed the presence of eleven out of twelve W-specific RAPD markers in the normal W chromosomes of twenty-five silkworm strains maintained in Japan. These results indicate that the W chromosomes of the strains in Japan are almost identical in type. The Zebra-W strain (T(W;3)Ze chromosome) lacked the W-Samurai and W-Mikan RAPD markers and the Black-egg-W strain (T(W;10)+(w-2) chromosome) lacked the W-Mikan RAPD marker. These results strongly indicate that the regions containing the W-Samurai and W-Mikan RAPD markers or the W-Mikan RAPD marker were deleted in the T(W;3)Ze and T(W;10)+(w-2) chromosomes, respectively, due to reciprocal translocation between the W chromosome and the autosome. This deletion apparently does not affect the expression of Fem; therefore, this deleted region of the W chromosome does not contain the putative Fem gene.

Reduced rates of sequence evolution of Y-linked satellite DNA in Rumex (Polygonaceae)

One characteristic of sex chromosomes is the accumulation of a set of different types of repetitive DNA sequences in the Y chromosomes. However, little is known about how this occurs or about how the absence of recombination affects the subsequent evolutionary fate of the repetitive sequences in the Y chromosome. Here we compare the evolutionary pathways leading to the appearance of three different families of satellite-DNA sequences within the genomes of Rumex acetosa and R. papillaris, two dioecious plant species with a complex XX/XY(1)Y(2) sex-chromosome system. We have found that two of these families, one autosomic (the RAE730 family) and one Y-linked (the RAYSI family), arose independently from the ancestral duplication of the same 120-bp repeat unit. Conversely, a comparative analysis of the three satellite-DNA families reveals no evolutionary relationships between these two and the third, RAE180, also located in the Y chromosomes. However, we have demonstrated that, regardless of the mechanisms that gave rise to these families, satellite-DNA sequences have different evolutionary fates according to their location in different types of chromosomes. Specifically, those in the Y chromosomes have evolved at half the rate of those in the autosomes, our results supporting the hypothesis that satellite DNAs in nonrecombining Y chromosomes undergo lower rates of sequence evolution and homogenization than do satellite DNAs in autosomes.

Retrotransposable elements on the W chromosome of the silkworm, Bombyx mori

The sex chromosomes of the silkworm, Bombyxmori, are designated ZW(XY) for females and ZZ(XX) for males. The W chromosome of B. mori does not recombine with the Z chromosome and autosomes and no genes for morphological characters have been mapped to the W chromosome as yet. Furthermore, femaleness is determined by the presence of a single W chromosome, regardless of the number of autosomes or Z chromosomes. To understand these interesting features of the W chromosome, it is necessary to analyze the W chromosome at the molecular biology level. Initially to isolate DNA sequences specific for the W chromosome as randomly amplified polymorphic DNA (RAPD) markers, we compared the genomic DNAs between males and females by PCR with arbitrary 10-mer primers. To the present, we have identified 12 W-specific RAPD markers, and with the exception of one RAPD marker, all of the deduced amino acid sequences of these W-specific RAPD markers show similarity to previously reported amino acid sequences of retrotransposable elements from various organisms. After constructing a genomic DNA lambda phage library of B. mori we obtained two lambda phage clones, one containing the W-Kabuki RAPD sequence and one containing the W-Samurai RAPD sequence and found that these DNA sequences comprised nested structures of many retrotransposable elements. To further analyze the W chromosome, we obtained 14 W-specific bacterial artificial chromosome (BAC) clones from three BAC libraries and subjected these clones to shotgun sequencing. The resulting assembly of sequences did not produce a single contiguous sequence due to the presence of many retrotransposable elements. Therefore, we coupled PCR with shotgun sequencing. Through these analyses, we found that many long terminal repeat (LTR) and non-LTR retrotransposons, retroposons, DNA transposons and their derivatives, have accumulated on the W chromosome as strata. These results strongly indicate that retrotransposable elements are the main structural component of the W chromosome.

Retroelements: tools for sex chromosome evolution

Many eukaryotic taxa inherit a heteromorphic sex chromosome pair. It is a generally accepted hypothesis that the sex chromosome pair is derived from a pair of homologous autosomes that has developed after the occurrence of a sex differentiator in an evolutionary process into two structurally and functionally different partners. In most of the analyzed systems the occurrence of the dominant sex differentiator is paralleled by the suppression of recombination within and close by that region. The recombinational isolation can spread in an evolutionary selection process from neighboring regions finally over the whole chromosome. Suppression of recombination strongly biases the distribution of retrotransposons in the genome. Our results and that from others indicate that the major force driving the evolution of Y chromosomes are retrotransposons, remodeling euchromatic chromosome structures into heterochromatic ones. In our model, intact or already eroded retrotransposons become trapped due to their inherent transposition mechanisms in non-recombining regions. The massive accumulation of retrotransposons interferes strongly with the activity of genes. We hypothesize that Y chromosome degeneration is a stepwise evolutionary process: (1) Massive accumulation of retrotransposons occurs in the non-recombining regions. (2) Heterochromatic nucleation centers are formed as a consequence of genomic defense against invasive parasitic elements; the established nucleation centers become epigenetically inherited. (3) Spreading of heterochromatin from the nucleation centers into flanking regions induces in an adaptive process gene silencing of neighbored genes that could either be still intact or in an already eroded condition, e.g., showing point mutations, deletions, insertions; the retroelements should be subjects to the same forces of deterioration as the genes themselves. (4) Constitutive silenced genes are not committed to the same genetic selection pressure as active genes and therefore more exposed to the decay process. (5) Gene dosage balance is reestablished by the parallel evolution of dosage compensation mechanisms. The evolving secondary sex chromosomes, neo-X and neo-Y, of Drosophila miranda are revealed to be a unique and potent model system to catch the evolutionary Y deterioration process in progress.

Molecular characterization and chromosomal distribution of Galileo, Kepler and Newton, three foldback transposable elements of the Drosophila buzzatii species complex

Galileo is a foldback transposable element that has been implicated in the generation of two polymorphic chromosomal inversions in Drosophila buzzatii. Analysis of the inversion breakpoints led to the discovery of two additional elements, called Kepler and Newton, sharing sequence and structural similarities with Galileo. Here, we describe in detail the molecular structure of these three elements, on the basis of the 13 copies found at the inversion breakpoints plus 10 additional copies isolated during this work. Similarly to the foldback elements described in other organisms, these elements have long inverted terminal repeats, which in the case of Galileo possess a complex structure and display a high degree of internal variability between copies. A phylogenetic tree built with their shared sequences shows that the three elements are closely related and diverged approximately 10 million years ago. We have also analyzed the abundance and chromosomal distribution of these elements in D. buzzatii and other species of the repleta group by Southern analysis and in situ hybridization. Overall, the results suggest that these foldback elements are present in all the buzzatti complex species and may have played an important role in shaping their genomes. In addition, we show that recombination rate is the main factor determining the chromosomal distribution of these elements.

Detachment analysis of the translocated W chromosome shows that the female-specific randomly amplified polymorphic DNA (RAPD) marker, female-218, is derived from the second chromosome fragment region of the translocated W chromosome of the sex-limited p(B) silkworm (Bombyx mori ) strain

The sex chromosomes of the silkworm, Bombyx mori, are designated ZW for the female and ZZ for the male. We previously characterized a female-specific randomly amplified polymorphic DNA (RAPD) marker, designated Female-218, from the translocation-bearing W chromosomes. These W chromosomes contain a region of the second chromosome, which carries visible larval markers of the p loci. We used strain TWPB in which female larvae have black skin due to the p(B) gene (T(W;2)p(B), +p/+p) while male larvae have whitish skin (+p/+p). To determine whether the Female-218 RAPD marker is derived from the "W region" or a "second chromosome fragment", we induced a detachment of the translocated W chromosome, T(W;2)p(B), by treating the eggs with hot water at an early developmental stage. After hot water treatment, we obtained 27 white female larvae out of 4850 female larvae. The Female-218 RAPD marker was not amplified in 26 out of 27 white female larvae, and was amplified from one white female larva. Moreover, we obtained 11 black male larvae out of 5377 male larvae. Eight out of 11 black male larvae became adult moths, and the Female-218 RAPD marker was amplified from all eight male moths. Examination of the genetic relationship between the Female-218 RAPD marker and the second chromosome fragment of the translocated W chromosome strongly indicates that the Female-218 RAPD marker is amplified from the region of second chromosome fragment of the T(W;2)p(B) chromosome.

Evolutionary course of CsRn1 long-terminal-repeat retrotransposon and its heterogeneous integrations into the genome of the liver fluke, Clonorchis sinensis

The evolutionary course of the CsRn1 long-terminal-repeat (LTR) retrotransposon was predicted by conducting a phylogenetic analysis with its paralog LTR sequences. Based on the clustering patterns in the phylogenetic tree, multiple CsRn1 copies could be grouped into four subsets, which were shown to have different integration times. Their differential sequence divergences and heterogeneous integration patterns strongly suggested that these subsets appeared sequentially in the genome of C. sinensis. Members of recently expanding subset showed the lowest level of divergence in their LTR and reverse transcriptase gene sequences. They were also shown to be highly polymorphic among individual genomes of the trematode. The CsRn1 element exhibited a preference for repetitive, agenic chromosomal regions in terms of selecting integration targets. Our results suggested that CsRn1 might induce a considerable degree of intergenomic variation and, thereby, have influenced the evolution of the C. sinensis genome.

Hybridization, transgressive segregation and evolution of new genetic systems inDrosophila

Introgressive hybridization facilitates incorporation of genes from one species into the gene pool of another. Studies on long-term effects of introgressive hybridization in animal systems are sparse. Drosophila nasuta (2n = 8) and D. albomicans (2n = 6)-a pair of allopatric, morphologically almost identical, cross-fertile members of the nasuta subgroup of the immigrans species group-constitute an excellent system to analyse the impact of hybridization followed by transgressive segregation of parental characters in the hybrid progeny. Hybrid populations of D. nasuta and D. albomicans maintained for over 500 generations in the laboratory constitute new recombinant hybrid genomes, here termed cytoraces. The impact of hybridization, followed by introgression and transgressive segregation, on chromosomal constitution and karyotypes, some fitness parameters, isozymes, components of mating behaviour and mating preference reveals a complex pattern of interracial divergence among parental species and cytoraces. This assemblage of characters in different combinations in a laboratory hybrid zone allows us to study the emergence of new genetic systems. Here, we summarize results from our ongoing studies comparing these hybrid cytoraces with the parental species, and discuss the implications of these findings for our understanding of the evolution of new genetic systems.

An Ac -like Transposable Element Family With Transcriptionally Active Y-Linked Copies in the White Campion, Silene latifolia

An RFLP genomic subtraction was used to isolate male-specific sequences in the species Silene latifolia. One isolated fragment, SLP2, shares similarity to a portion of the Activator (Ac) transposase from Zea mays and to related proteins from other plant species. Southern blot analysis of male and female S. latifolia genomic DNA shows that SLP2 belongs to a low-copy-number repeat family with two Y-linked copies. Screening of a S. latifolia male genomic library using SLP2 as a probe led to the isolation of five clones, which were partially sequenced. One clone contains two large open reading frames that can be joined into a sequence encoding a putative protein of 682 amino acids by removing a short intron. Database searches and phylogenetic analysis show that this protein belongs to the hAT superfamily of transposases, closest to Tag2 (Arabidopsis thaliana), and contains all of the defined domains critical for the activity of these transposases. PCR with genomic and cDNA templates from S. latifolia male, female, and hermaphrodite individuals revealed that one of the Y-linked copies is transcriptionally active and alternatively spliced. This is the first report of a transcriptionally active transposable element (TE) family in S. latifolia and the first DNA transposon residing on a plant Y chromosome. The potential activity and regulation of this TE family and its use for Y chromosome gene discovery is discussed.

Molecular characterization of the genomic region linked with apomixis in Pennisetum/Cenchrus

Apomixis is defined as asexual reproduction through seeds, although this outcome can be achieved by multiple pathways. Since little is known about the molecular control of these pathways, how they might intersect is also a mystery. Two of these pathways in the grass family, diplospory and apospory, are receiving attention from molecular biologists. Apospory in Pennisetum/Cenchrus, two genera of panicoid grasses, results in the formation of four-nucleate embryo sacs that lack antipodals. Sexual reproduction frequently aborts so that the resulting seed is composed of (1) a parthenogenetically derived embryo that is genetically identical to the mother and (2) endosperm formed through pseudogamy. The transmission of apomixis is associated with the transfer of a linkage block on a single chromosome. This linkage block contains repetitive sequences as well as hemizygous, low-copy DNA sequences. Fluorescence in situ hybridization has demonstrated that these DNA regions occur on only a single chromosome, but not its homologs, in the polyploid apomicts studied. Features of the apomixis-associated region resemble those of other chromosomal segments isolated from recombination and replete with "selfish" DNAs.

A survey of the DNA sequences surrounding the Bari1 repeats in the pericentromeric h39 region of Drosophila melanogaster

In Drosophila melanogaster, clustered copies of the Bari1 transposon are only present in the pericentromeric h39 region of the second chromosome, where other clusters of repetitive elements, either found organized in large tandem arrays only in the h39 region (Responder, PortoI), or both in the h39 region and in other heterochromatic regions (Hoppel), are also observed. The topological relationship among the repetitive sequences of the h39 region and the nature of the sequences separating its large repeat clusters are at present largely unknown. To get new insights on the sequence composition of the heterochromatin and on the forces governing its origin and maintenance, we have cloned and analyzed part of the DNA sequences flanking the h39 Bari1 repeats. In a region spanning 3 and 9 kb, respectively, from the ends of a Bari1 array we found only single copies of the PortoI and Hoppel transposable elements, and five copies of a variant form of the Responder repeats. No large tandem arrays of any repeated element were present. In addition, a highly conserved 596 bp sequence, that may have a functional role, is present on both sides of the Bari1 repeats. We suggest that the current organization of the h39 heterochromatin implies some topological or functional constraint that prevents the formation of further arrays of repetitive elements in the region.

Identification and characterization of the ARIADNE gene family in Arabidopsis. A group of putative E3 ligases

ARIADNE (ARI) proteins were recently identified in fruitfly (Drosophila melanogaster), mouse, and man because of their specific interaction with the ubiquitin-conjugating (E2) enzymes UbcD10, UbcM4, UbcH7, and UbcH8. They are characterized by specific motifs and protein structures that they share with PARKIN, and there is increasing evidence that ARI/PARKIN proteins function as E2-dependent ubiquitin-protein ligases. On the basis of homology and motif searches, 16 AtARI genes were identified in Arabidopsis. Analysis of the position of exons/introns and their chromosomal localization indicates that the AtARI gene family expanded via larger and smaller genome duplications. We present evidence that retroposition of processed mRNA may have also contributed to enlarging this gene family. Phylogenetic analyses divides the AtARI proteins into three subgroups. Two groups are absent in yeast, invertebrates, and vertebrates and may therefore represent new plant-specific subfamilies. Examination of the predicted protein sequences revealed that the ARI proteins share an additional leucine-rich region at the N terminus that is highly conserved in all phyla analyzed. Furthermore, conserved consensus signals for casein kinase II-dependent phosphorylation and for nuclear localization were identified. The in silico-based analyses were complemented with experimental data to quantify expression levels. Using real-time polymerase chain reaction, we show that the ARI genes are differentially transcribed. AtARI1 is highly expressed in all organs, whereas no transcripts could be detected for AtARI11, AtARI13, and AtARI14. AtARI12 and AtARI16 are expressed in an organ-specific manner in the roots and siliques, respectively.

Nested retrotransposons on the W chromosome of the wild silkworm Bombyx mandarina

The W chromosome of the silkworms Bombyx mori or B. mandarina is recombinationally isolated from the Z chromosome and the autosomes. We previously characterized a female-specific randomly amplified polymorphic DNA (RAPD), designated W-Yamato, derived from the W chromosome of the wild silkworm Bombyx mandarina. To further analyse the W chromosome of B. mandarina, we obtained a lambda phage clone that contains the W-Yamato RAPD sequence and sequenced the 16.7 kb DNA insert. We found that this DNA comprises a nested structure of at least seven elements: six retrotransposons and one transposable element-like sequence. The transposable element-like sequence is inserted into a micropia-like retrotransposon (Karate). The Karate and the non-long terminal repeat (non-LTR) retrotransposon BMC1 are inserted into a 412-like retrotransposon (Judo). Furthermore, this Judo, and two non-LTR retrotransposons (Kurosawa and Kendo) are inserted into a Pao-like retrotransposon (Yamato). These results indicate that the retrotransposons inserted into the W chromosome are not efficiently removed but accumulate gradually as strata without recombination.

Variability of genetic sex determination in poeciliid fishes

Poeciliids are one of the best-studied groups of fishes with respect to sex determination. They present an amazing variety of mechanisms, which span from simple XX-XY or ZZ-ZW systems to polyfactorial sex determination. The gonosomes of poeciliids generally are homomorphic, but very early stages of sex chromosome differentiation have been occasionally detected in some species. In the platyfish Xiphophorus maculatus, gene loci involved in melanoma formation, in different pigmentation patterns and in sexual maturity are closely linked to the sex-determining locus in the subtelomeric region of the X- and Y- chromosomes. The majority of traits encoded by these loci are highly polymorphic. This phenomenon might be explained by the high level of genomic plasticity apparently affecting the sex-determining region, where frequent rearrangements such as duplications, deletions, amplifications, and transpositions frequently occur. We propose that the high plasticity of the sex-determining region might explain the variability of sex determination in Xiphophorus and other poeciliids.

Jule from the fish Xiphophorus is the first complete vertebrate Ty3/Gypsy retrotransposon from the Mag family

Jule is the second complete long-terminal-repeat (LTR) Ty3/Gypsy retrotransposon identified to date in vertebrates. Jule, first isolated from the poeciliid fish Xiphophorus maculatus, is 4.8 kb in length, is flanked by two 202-bp LTRs, and encodes Gag (structural core protein) and Pol (protease, reverse transcriptase, RNase H, and integrase, in that order) but no envelope. There are three to four copies of Jule per haploid genome in X. maculatus. Two of them are located in a subtelomeric region of the sex chromosomes, where they are associated with the Xmrk receptor tyrosine kinase genes, of which oncogenic versions are responsible for the formation of hereditary melanoma in Xiphophorus. One almost intact copy of Jule was found in the first intron of the X-chromosomal allele of the Xmrk proto-oncogene, and a second, more corrupted copy is present only 56 nt downstream of the polyadenylation signal of the Xmrk oncogene. Jule-related elements were detected by Southern blot hybridization with less than 10 copies per haploid genome in numerous other poeciliids, as well as in more divergent fishes, including the medakafish Oryzias latipes and the tilapia Oreochromis niloticus. Database searches also identified Jule-related sequences in the zebrafish Danio rerio and in both genome project pufferfishes, Fugu rubripes and Tetraodon nigroviridis. Phylogenetic analysis revealed that Jule is the first member of the Mag family of Ty3/Gypsy retrotransposons described to date in vertebrates. This family includes the silkworm Mag and sea urchin SURL retrotransposons, as well as sequences from the nematode Caenorhabditis elegans. Additional related elements were identified in the genomes of the malaria mosquito Anopheles gambiae and the nematode Ascaris lumbricoides. Phylogeny of Mag-related elements suggested that the Mag family of retrotransposons is polyphyletic and is constituted of several ancient lineages that diverged before their host genomes more than 600 MYA.

Perspective: transposable elements, parasitic DNA, and genome evolution

The nature of the role played by mobile elements in host genome evolution is reassessed considering numerous recent developments in many areas of biology. It is argued that easy popular appellations such as "selfish DNA" and "junk DNA" may be either inaccurate or misleading and that a more enlightened view of the transposable element-host relationship encompasses a continuum from extreme parasitism to mutualism. Transposable elements are potent, broad spectrum, endogenous mutators that are subject to the influence of chance as well as selection at several levels of biological organization. Of particular interest are transposable element traits that early evolve neutrally at the host level but at a later stage of evolution are co-opted for new host functions.

The Origin of Interspecific Genomic Incompatibility via Gene Duplication

One of the great unsolved mysteries of evolutionary biology concerns the genetic mechanisms underlying the origin of genomic incompatibilities between species. Two prevailing thoughts are that such incompatibilities often result from epistatically interacting genes that act as loss-of-function alleles in hybrid backgrounds or from chromosomal rearrangements that result in mis-segregation during meiosis in hybrids. However, it is unclear how genes that cause a radical breakdown in hybrids arise without reducing fitness within species, and numerous cases of speciation appear to be unassociated with obvious chromosomal rearrangements. Here we suggest that duplicate genes, and more generally any kind of genomic redundancies, provide a powerful substrate for the origin of genomic incompatibilities in isolated populations. The divergent resolution of genomic redundancies, such that one population loses function from one copy while the second population loses function from a second copy at a different chromosomal location, leads to chromosomal repatterning such that gametes produced by hybrid individuals can be completely lacking in functional genes for a duplicate pair. Under this model, incompatibility factors accumulate with essentially no loss of fitness within populations as postulated under the Bateson-Dobzhansky-Muller (BDM) model of speciation and despite the fact that they arise from degenerative mutations. However, unlike the situation often envisioned under the BDM model, no change in the mode of gene action in hybrid backgrounds need be invoked. The plausibility of this model derives from a number of recent observations, including the fact that most genomes harbor substantial numbers of gene duplicates whose turnover is common and ongoing process and the fact that many genes have complex regulatory regions that facilitate their divergent resolution in sister taxa.

B-chromosome evolution

B chromosomes are extra chromosomes to the standard complement that occur in many organisms. They can originate in a number of ways including derivation from autosomes and sex chromosomes in intra- and interspecies crosses. Their subsequent molecular evolution resembles that of univalent sex chromosomes, which involves gene silencing, heterochromatinization and the accumulation of repetitive DNA and transposons. B-chromosome frequencies in populations result from a balance between their transmission rates and their effects on host fitness. Their long-term evolution is considered to be the outcome of selection on the host genome to eliminate B chromosomes or suppress their effects and on the B chromosome's ability to escape through the generation of new variants. Because B chromosomes interact with the standard chromosomes, they can play an important role in genome evolution and may be useful for studying molecular evolutionary processes.

Revising the selfish DNA hypothesis: new evidence on accumulation of transposable elements in heterochromatin

The bulk of the eukaryotic genome is composed of families of repetitive sequences that are genetically silent and exhibit various types of instability. Transposable elements (TEs) are particularly common in heterochromatic regions of the genome - a location where TEs might do less damage to their host. Recent advances suggest that the relationship between TEs and heterochromatin might not be quite so straightforward.

Retroposition of autosomal mRNA yielded testis-specific gene family on human Y chromosome

Most genes in the human NRY (non-recombining portion of the Y chromosome) can be assigned to one of two groups: X-homologous genes or testis-specific gene families with no obvious X-chromosomal homologues. The CDY genes have been localized to the human Y chromosome, and we report here that they are derivatives of a conventional single-copy gene, CDYL (CDY-like), located on human chromosome 13 and mouse chromosome 6. CDY genes retain CDYL exonic sequences but lack its introns. In mice, whose evolutionary lineage diverged before the appearance of the Y-linked derivatives, the autosomal Cdyl gene produces two transcripts; one is expressed ubiquitously and the other is expressed in testes only. In humans, autosomal CDYL produces only the ubiquitous transcript; the testis-specific transcript is the province of the Y-borne CDY genes. Our data indicate that CDY genes arose during primate evolution by retroposition of a CDYL mRNA and amplification of the retroposed gene. Retroposition contributed to the gene content of the human Y chromosome, together with two other molecular evolutionary processes: persistence of a subset of genes shared with the X chromosome and transposition of genomic DNA harbouring intact transcription units.

Identification of novel random amplified polymorphic DNAs (RAPDs) on the W chromosome of the domesticated silkworm, Bombyx mori, and the wild silkworm, B. mandarina, and their retrotransposable element-related nucleotide sequences

Genomic DNAs were compared between males and females of the domesticated silkworm, Bombyx mori, strains C108, C137, J137, p50, and WILD-W (constructed by crossing a wild silkworm, B. mandarina, female with a male of strain C108) by polymerase chain reaction (PCR) with 700 arbitrary 10-mer primers. Four female-specific RAPDs (W-Kabuki, W-Samurai, W-Kamikaze, and W-Yamato) were found. The sex chromosome formulas of B. mori and B. mandarina are ZW (XY) for the female and ZZ (XX) for the male. The four female-specific RAPDs are assumed to be derived from the W chromosome because the other chromosomes are shared by both sexes. A computer search for deduced amino acid sequences of these four RAPDs revealed that all of them showed homology to previously reported amino acid sequences encoded in known retrotransposable elements from various organisms.

Conspiracy of silence among repeated transgenes

Transgenic experiments in vertebrates often involve the insertion of tandem multiple-copy arrays at single sites. For many transgenes, expression is unpredictable from site to site, a phenomenon usually attributed to a repressive environment caused by nearby sequences. However, an alternative explanation comes from evidence that transgene repeat arrays in flies condense into heterochromatin, suggesting that low levels of expression in vertebrate transgene arrays might result from interactions between repeats within the array. A recent experiment using transgenic mouse lines demonstrates that reduction in copy number of silenced transgenes within an array leads to a striking increase in expression, demonstrating that silencing is intrinsic to the array, and is not attributable to position effects of nearby sequences. This work calls into question functions that have been attributed to vertebrate locus control regions and boundaries, and draws attention to the notion that repeat-induced gene silencing is a system for protection of eukaryotic genomes against threatening sequence elements.

Transposable element-host interactions: regulation of insertion and excision

Transposable elements propagate by inserting into new locations in the genomes of the hosts they inhabit. Their transposition might thus negatively affect the fitness of the host, suggesting the requirement for a tight control in the regulation of transposable element mobilization. The nature of this control depends on the structure of the transposable element. DNA elements encode a transposase that is necessary, and in most cases sufficient, for mobilization. In general, regulation of these elements depends on intrinsic factors with little direct input from the host. Retrotransposons require an RNA intermediate for transposition, and their frequency of mobilization is controlled at multiple steps by the host genome by regulating both their expression levels and their insertional specificity. As a result, a symbiotic relationship has developed between transposable elements and their host. Examples are now emerging showing that transposons can contribute significantly to the well being of the organisms they populate.