Junctions between repetitive DNAs on the PSR chromosome of Nasonia vitripennis: association of palindromes with recombination

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.

Adjacent sequences disclose potential for intra-genomic dispersal of satellite DNA repeats and suggest a complex network with transposable elements

Background

Satellite DNA (satDNA) sequences are typically arranged as arrays of tandemly repeated monomers. Due to the similarity among monomers, their organizational pattern and abundance, satDNAs are hardly accessible to structural and functional studies and still represent the most obscure genome component. Although many satDNA arrays of diverse length and even single monomers exist in the genome, surprisingly little is known about transition from satDNAs to other sequences. Studying satDNA monomers at junctions and identifying DNA sequences adjacent to them can help to understand the processes that (re)distribute satDNAs and significance that evolution of these sequence elements might have in creating the genomic landscape.

Results

We explored sets of randomly selected satDNA-harboring genomic fragments in four mollusc species to examine satDNA transition sites, and the nature of adjacent sequences. All examined junctions are characterized by abrupt transitions from satDNAs to other sequences. Among them, junctions of only one examined satDNA mapped non-randomly (within the palindrome), indicating that well-defined sequence feature is not a necessary prerequisite in the junction formation. In the studied sample, satDNA flanking sequences can be roughly classified into two groups. The first group is composed of anonymous DNA sequences which occasionally include short segments of transposable elements (TEs) as well as segments of other satDNA sequences. In the second group, satDNA repeats and the array flanking sequences are identified as parts of TEs of the Helitron superfamily. There, some array flanking regions hold fragmented satDNA monomers alternating with anonymous sequences of comparable length as missing monomer parts, suggesting a process of sequence reorganization by a mechanism able to excise short monomer parts and replace them with unrelated sequences.

Conclusions

The observed architecture of satDNA transition sites can be explained as a result of insertion and/or recombination events involving short arrays of satDNA monomers and TEs, in combination with hypothetical transposition-related ability of satDNA monomers to be shuffled independently in the genome. We conclude that satDNAs and TEs can form a complex network of sequences which essentially share the propagation mechanisms and in synergy shape the genome.

Electronic supplementary material

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

Evolutionary dynamics and sites of illegitimate recombination revealed in the interspersion and sequence junctions of two nonhomologous satellite DNAs in cactophilic Drosophila species

Satellite DNA (satDNA) is a major component of genomes but relatively little is known about the fine-scale organization of unrelated satDNAs residing at the same chromosome location, and the sequence structure and dynamics of satDNA junctions. We studied the organization and sequence junctions of two nonhomologous satDNAs, pBuM and DBC-150, in three species from the neotropical Drosophila buzzatii cluster (repleta group). In situ hybridization to microchromosomes, interphase nuclei and extended DNA fibers showed frequent interspersion of the two satellites in D. gouveai, D. antonietae and, to a lesser extent, D. seriema. We isolated by PCR six pBuM x DBC-150 junctions: four are exclusive to D. gouveai and two are exclusive to D. antonietae. The six junction breakpoints occur at different positions within monomers, suggesting independent origin. Four junctions showed abrupt transitions between the two satellites, whereas two junctions showed a distinct 10 bp tandem duplication before the junction. Unlike pBuM, DBC-150 junction repeats are more variable than randomly cloned monomers and showed diagnostic features in common to a 3-monomer higher-order repeat seen in the sister species D. serido. The high levels of interspersion between pBuM and DBC-150 repeats suggest extensive rearrangements between the two satellites, maybe favored by specific features of the microchromosomes. Our interpretation is that the junctions evolved by multiples events of illegitimate recombination between nonhomologous satDNA repeats, with subsequent rounds of unequal crossing-over expanding the copy number of some of the junctions.

Satellite DNA junctions identify the potential origin of new repetitive elements in the beetle Tribolium madens

Two related satellite DNA families (satellite I and satellite II) with complex higher-order repeat (HOR) monomers represent major DNA components equilocated in the pericentromeric heterochromatin of all Tribolium madens chromosomes. Fragments obtained upon genomic DNA restriction revealed two subfamilies of satellite II monomers, and also identified regions of transition between satellite I and satellite II sequences. The two subfamilies differ not only in diagnostic nucleotides, but also in flipped orientation of constituent subunits. Hybrid genomic fragments comprise directly linked satellite I and satellite II monomers that cannot be distinguished from randomly cloned monomers of corresponding families. An exception is the most proximal satellite I monomer in the hybrid fragment named TMADhinf, which shows sequence divergence typical for repeats evolving at array ends, in zones of low homogenization efficiency. This pattern points to the extensive rearrangement processes generating abrupt transitions between satellite arrays combined with array maintenance by unequal crossover. Switching points between adjacent satellites as well as the edges of flipped subunits are localized within a short sequence segment, indicating a preferential site of recombination within satellite subunits. Multiple copies of TMADhinf junction fragment support the hypothesis that sites of evolutionary origin of novel satellite repeat (sub)families can be localized at array ends, in regions of enhanced sequence divergence.

The intragenomic polymorphism of a partially inverted repeat (PIR) in Gallus gallus domesticus, potential role of inverted repeats in satellite DNAs evolution

We report here the molecular characterization of the basic repeating unit of a novel repetitive family, partially inverted repeat (PIR), previously identified from chicken genome. This repetitive DNA family shares a close evolutionary relationship with XhoI/EcoRI repeats and chicken nuclear-membrane-associated (CNM) repeat. Sequence analyses reveal the 1430 bp basic repeating unit can be divided into two regions: the central region ( approximately 1000 bp) and the flanking region ( approximately 430 bp). Within the central region, a pair of repeats (86 bp) flanks the central core ( approximately 828 bp) in inversed orientation. Due to the tandem array feature shared by the repeating units, the inverted repeats fall between the central core and flanking region. Southern blot analyses further reveal the intragenomic polymorphism of PIR, and the molecular size of repeating units ranges from 1.1 kb to 1.6 kb. The identified monomer variants may result from multiple crossing-over events, implying the potential roles of inverted repeats in satellite DNAs variation.

NOR activity and repeat sequences of the paternal sex ratio chromosome of the parasitoid wasp Trichogramma kaykai

Part of the male population of the wasp Trichogramma kaykai carries a B chromosome that manipulates its host sex ratio in favour of males. The only known repeat on this paternal sex ratio (PSR) chromosome is the 45S rDNA, which includes here five different internal transcribed spacer 2 (ITS2) sequences. In this report, we describe that only part of these ITS2 sequences is transcribed. The absence of transcription of some ITS2 sequences might explain the presence of multiple ITS2 sequences on the PSR chromosome since homogenization of rDNA spacers is thought to occur only in transcribed regions. Analysis of the only other known tandem repeat in Trichogramma, the EcoRI repeat, showed that it is absent from the PSR chromosome, and that the T. kaykai EcoRI repeat has 98 and 77% DNA sequence homology with the T. deion and T. brassicae EcoRI repeats, respectively. The size of the PSR chromosome measures 9 Mbp and is equal to 3.9% of the haploid T. kaykai genome. Finally, fluorescent in situ hybridization with a pool of high and moderate repetitive T. kaykai DNA (C0t-50) revealed only a very few major tandem repeats on the Trichogramma genome and only 45S rDNA on the PSR chromosome.

Mapping of paternal-sex-ratio deletion chromosomes localizes multiple regions involved in expression and transmission

The paternal-sex-ratio (PSR) chromosome in the parasitic wasp Nasonia vitripennis is a submetacentric supernumerary (B chromosome). Males transmit PSR, but after fertilization it causes the loss of the paternal autosomes. Paternal genome loss caused by PSR results in the conversion of a female (diploid) zygote into a male (haploid) under haplodiploid sex determination. In this study, site-specific markers were developed to assay deletion derivatives of PSR. Both polymerase chain reaction and Southern hybridization were used to detect the presence/absence of 16 single-site markers on a set of 20 functional and nine nonfunctional deletion chromosomes. Based on the pattern of marker loss on the deletion chromosomes, the basic organization of PSR was revealed. Two sets of markers were deleted independently, apparently representing the two arms of the submetacentric chromosome. The presence or absence of specific regions was examined in relation to phenotypic characteristics of the deletion chromosomes; ability to cause paternal genome loss, and stability in mitotic cell divisions. Rather than identifying a single region on PSR as being responsible for PSR function, the results suggest that the retention of one of two chromosomal regions is sufficient for causing paternal genome loss. Furthermore, a region was identified that is tightly correlated with mitotic stability, as measured from chromosomal transmission rates. Functional chromosomes with short-arm deletions had high (approximately 100%) transmission rates, whereas functional chromosomes with long-arm deletions had low (approximately 85%) transmission rates.

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.

The molecular organisation of a B chromosome tandem repeat sequence from Brachycome dichromosomatica

A high copy, tandemly repeated, sequence (Bd49) specific to the B chromosome and located near the centromere in Brachycome dichromosomatica was used to identify lambda genomic clones from DNA of a 3B plant. Only one clone of those analysed was composed entirely of a tandem array of the B-specific repeat unit. In other clones, the Bd49 repeats were linked to, or interspersed with, sequences that are repetitious and distributed elsewhere on the A and B chromosomes. One such repetitious flanking sequence has similarity to retrotransposon sequences and a second is similar to chloroplast DNA sequences. Of the four separate junctions analysed of Bd49-like sequence with flanking sequence, three were associated with the same A/T-rich region in Bd49 and the fourth was close to a 25 bp imperfect dyadic sequence. No novel B-specific sequences were detected within the genomic clones.

Isolation and characterization of a retroelement from B chromosome (PSR) in the parasitic wasp Nasonia vitripennis

Molecular characterization of the paternal-sex-ratio (PSR) chromosome in Nasonia vitripennis (Hymenoptera: Pteromalidae) has led to the isolation of a dispersed repetitive element. The element is a LTR-containing retrotransposon which has been named NATE (NAsonia Transposable Element). NATE has direct terminal repeats and has an internal amino acid sequence similar to reverse transcriptases of other retroelements. Phylogenetic analysis indicates NATE is a member of the Gypsy/Ty3 group of retrotransposons, and represents the first isolated from Hymenoptera. Five closely related copies of NATE were isolated from the PSR chromosome, but cross-hybridizing elements were not detected on the autosomes of N. vitripennis. Strongly cross-hybridizing elements were, however, detected in two other Nasonia species. This observed distribution of NATE is interesting, because the supernumerary PSR chromosome may be derived from the genome of a sibling species of N. vitripennis.

Characterisation of repeated sequences from microdissected B chromosomes of Crepis capillaris

The B chromosome of Crepis capillaris was isolated from the standard chromosomes by microdissection, and the chromosomal DNA amplified using the degenerate oligonucleotide-primed polymerase chain reaction (DOP-PCR). The PCR product was cloned and a B-specific library created and characterised. Southern and in situ hybridisation analyses of the DOP-PCR product from microdissected B chromosomes confirmed that the B chromosome is composed mainly of sequences also present in the A chromosomes but lacks the main repeated DNA families located in the A-chromosomal heterochromatin. From 100 clones analysed, 12% of the generated B-chromosomal library was shown to be composed of dispersed repeats located in both the A and B chromosomes. No B-specific repeated sequence was detected. One of the most abundant repeated DNAs within the library, the family B134, was further characterised. Repeating units show a sequence similarity range from 69% to 90% and are characterised by their richness in (CA)n repeats. In situ hybridisation revealed that members of this family are dispersed throughout the A and B chromosomes but are more concentrated in the pericentromeric heterochromatin of the B, indicating that the molecular organization of B heterochromatin is different from that of the A chromosomes. Compared with the A chromosomes, the Bs contain about 20,000 copies per micron more of the B134 sequence. This indicates that B134 was amplified on the B chromosome after its origin. The B134 sequences in the B chromosomes have also diverged from those on the A chromosomes. Although the DNA composition of A and B chromosomes is similar, Bs are evolving separately from A chromosomes at the molecular level.

Organisation and origin of a B chromosome centromeric sequence from Brachycome dichromosomatica

Brachycome dichromosomatica is an Australian native daisy that has two pairs of A chromosomes and up to three B chromosomes in some populations. A putative B-specific tandem repeat DNA sequence (Bd49) was isolated previously. Here we describe further characterisation of this sequence and investigate its possible origin. Southern analysis showed that all individual B chromosomes examined have highly methylated tandem repeats of Bd49 but differences in banding pattern for distinct B isolates suggested that the sequence is in a state of flux. Using in situ hybridisation, the sequence was shown to be located at the centromeric region of the B chromosome. Southern analysis of genomic DNA with Bd49 demonstrated that multiple copies of the sequence exist in the genomes of B. eriogona, B. ciliaris, B. segmentosa and B. multifida (none of which have B chromosomes) whereas other species tested (including 0B plants of B. dichromosomatica and 0B and +B B. curvicarpa and B. dentata) have few or no copies. Genomic clones and Bd49-like sequences derived by the polymerase chain reaction (PCR) were obtained from five species but determination of phylogenetic relationships within the genus and inference as to the possible origin of the B chromosome were problematic because of extensive intragenomic heterogeneity of the sequences.

Hypervariability of ribosomal DNA at multiple chromosomal sites in lake trout (Salvelinus namaycush)

Variation in the intergenic spacer (IGS) of the ribosomal DNA (rDNA) of lake trout (Salvelinus namaycush) was examined. Digestion of genomic DNA with restriction enzymes showed that almost every individual had a unique combination of length variants with most of this variation occurring within rather than between populations. Sequence analysis of a 2.3 kilobase (kb) EcoRI-DraI fragment spanning the 3' end of the 28S coding region and approximately 1.8 kb of the IGS revealed two blocks of repetitive DNA. Putative transcriptional termination sites were found approximately 220 bases (b) downstream from the end of the 28S coding region. Comparison of the 2.3-kb fragments with two longer (3.1 kb) fragments showed that the major difference in length resulted from variation in the number of short (89 b) repeats located 3' to the putative terminator. Repeat units within a single nucleolus organizer region (NOR) appeared relatively homogeneous and genetic analysis found variants to be stably inherited. A comparison of the number of spacer-length variants with the number of NORs found that the number of length variants per individual was always less than the number of NORs. Examination of spacer variants in five populations showed that populations with more NORs had more spacer variants, indicating that variants are present at different rDNA sites on nonhomologous chromosomes.

The evolutionary dynamics of repetitive DNA in eukaryotes

Repetitive DNA sequences form a large portion of the genomes of eukaryotes. The 'selfish DNA' hypothesis proposes that they are maintained by their ability to replicate within the genome. The behaviour of repetitive sequences can result in mutations that cause genetic diseases, and confer significant fitness losses on the organism. Features of the organization of repetitive sequences in eukaryotic genomes, and their distribution in natural populations, reflect the evolutionary forces acting on selfish DNA.