Sequence divergence and conservation in organizationally distinct subfamilies of Donax trunculus satellite DNA

Sequence characterization of the 5S ribosomal DNA and the internal transcribed spacer (ITS) region in four European Donax species (Bivalvia: Donacidae)

Background

The whole repeat unit of 5S rDNA and the internal transcribed spacer (ITS) of four European Donax species were analysed. After amplifying, cloning and sequencing several 5S and ITS units, their basic features and their variation were described. The phylogenetic usefulness of 5S and ITS sequences in the inference of evolutionary relationships among these wedge clams was also investigated.

Results

The length of the 5S repeat presented little variation among species, except D. trunculus that differed from the rest of the Donax species in 170–210 bp. The deduced coding region covered 120 bp, and showed recognizable internal control regions (ICRs) involved in the transcription. The length of non-transcribed spacer region (NTS) ranged from 157 bp to 165 bp in Donax trunculus and from 335 bp to 367 bp in the other three species. The conservation degree of transcriptional regulatory regions was analysed revealing a conserved TATA-like box in the upstream region. Regarding ITS sequences, the four Donax species showed slight size differences among clones due to the variation occurring in the ITS1 and ITS2, except Donax variegatus did not display size differences in the ITS2. The total length of the ITS sequence ranged between 814 and 1014 bp. Resulting phylogenetic trees display that the two ribosomal DNA regions provide well-resolved phylogenies where the four European Donax species form a single clade receiving high support in nodes. The topology obtained with 5S sequences was in agreement with Donax evolutionary relationships inferred from several sequences of different nature in previous studies.

Conclusions

This is not only a basic research work, where new data and new knowledge is provided about Donax species, but also have allowed the authentication of these wedge clams and offers future applications to provide other genetic resources.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0684-x) contains supplementary material, which is available to authorized users.

New insights into the genome repetitive fraction of the Antarctic bivalve Adamussium colbecki

Repetitive DNA represents the major component of the genome in both plant and animal species. It includes transposable elements (TEs), which are dispersed throughout the genome, and satellite DNAs (satDNAs), which are tandemly organized in long arrays. The study of the structure and organization of repetitive DNA contributes to our understanding of genome architecture and the mechanisms leading to its evolution. Molluscs represent one of the largest groups of invertebrates and include organisms with a wide variety of morphologies and lifestyles. To increase our knowledge of bivalves at the genome level, we analysed the Antarctic scallop Adamussium colbecki. The screening of the genomic library evidenced the presence of two novel satDNA elements and the CvA transposon. The interspecific investigation performed in this study demonstrated that one of the two satDNAs isolated in A. colbecki is widespread in polar molluscan species, indicating a possible link between repetitive DNA and abiotic factors. Moreover, the transcriptional activity of CvA and its presence in long-diverged bivalves suggests a possible role for this ancient element in shaping the genome architecture of this clade.

Are Pericentric Inversions Reorganizing Wedge Shell Genomes?

Wedge shells belonging to the Donacidae family are the dominant bivalves in exposed beaches in almost all areas of the world. Typically, two or more sympatric species of wedge shells differentially occupy intertidal, sublittoral, and offshore coastal waters in any given locality. A molecular cytogenetic analysis of two sympatric and closely related wedge shell species, Donax trunculus and Donax vittatus, was performed. Results showed that the karyotypes of these two species were both strikingly different and closely alike; whilst metacentric and submetacentric chromosome pairs were the main components of the karyotype of D. trunculus, 10–11 of the 19 chromosome pairs were telocentric in D. vittatus, most likely as a result of different pericentric inversions. GC-rich heterochromatic bands were present in both species. Furthermore, they showed coincidental 45S ribosomal RNA (rRNA), 5S rRNA and H3 histone gene clusters at conserved chromosomal locations, although D. trunculus had an additional 45S rDNA cluster. Intraspecific pericentric inversions were also detected in both D. trunculus and D. vittatus. The close genetic similarity of these two species together with the high degree of conservation of the 45S rRNA, 5S rRNA and H3 histone gene clusters, and GC-rich heterochromatic bands indicate that pericentric inversions contribute to the karyotype divergence in wedge shells.

Methylation profile of a satellite DNA constituting the intercalary G+C-rich heterochromatin of the cut trough shell Spisula subtruncata (Bivalvia, Mactridae)

Tandemly repeated DNAs usually constitute significant portions of eukaryotic genomes. In bivalves, however, repetitive DNAs are habitually not widespread. In our search for abundant repetitive DNAs in trough shells, we discovered a novel satellite DNA, SSUsat, which constitutes at least 1.3% of the genome of Spisula subtruncata. As foreseen by the satellite DNA library hypothesis, we confirmed that this satellite DNA is also present in two other Mactridae species, showing a highly conserved nucleotide sequence together with a dramatic diminution in the number of repeats. Predominantly located at the G + C-rich intercalary heterochromatin of S. subtruncata, SSUsat displays several DNA methylation peculiarities. The level of methylation of SSUsat is high (3.38%) in comparison with bivalve standards and triplicates the mean of the S. subtruncata genome (1.13%). Methylation affects not only the cytosines in CpG dinucleotides but also those in CHH and CHG trinucleotides, a feature common in plants but scarce and without any clear known relevance in animals. SSUsat segments enriched in methylated cytosines partly overlap those showing higher sequence conservation. The presence of a chromosome pair showing an accumulation of markedly under-methylated SSUsat monomers additionally indicates that the methylation processes that shape repetitive genome compartments are quite complex.

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.

A High-Resolution Comparative Chromosome Map of Cricetus cricetus and Peromyscus eremicus Reveals the Involvement of Constitutive Heterochromatin in Breakpoint Regions

Compared to humans and other mammals, rodent genomes, specifically Muroidea species, underwent intense chromosome reshuffling in which many complex structural rearrangements occurred. This fact makes them preferential animal models for studying the process of karyotype evolution. Here, we present the first combined chromosome comparative maps between 2 Cricetidae species, Cricetus cricetus and Peromyscus eremicus, and the index species Mus musculus and Rattus norvegicus. Comparative chromosome painting was done using mouse and rat paint probes together with in silico analysis from the Ensembl genome browser database. Hereby, evolutionary events (inter- and intrachromosomal rearrangements) that occurred in C. cricetus and P. eremicus since the putative ancestral Muroidea genome could be inferred, and evolutionary breakpoint regions could be detected. A colocalization of constitutive heterochromatin and evolutionary breakpoint regions in each genome was observed. Our results suggest the involvement of constitutive heterochromatin in karyotype restructuring of these species, despite the different levels of conservation of the C. cricetus (derivative) and P. eremicus (conserved) genomes.

Tandem repeat-containing MITEs in the clam Donax trunculus

Two distinct classes of repetitive sequences, interspersed mobile elements and satellite DNAs, shape eukaryotic genomes and drive their evolution. Short arrays of tandem repeats can also be present within nonautonomous miniature inverted repeat transposable elements (MITEs). In the clam Donax trunculus, we characterized a composite, high copy number MITE, named DTC84. It is composed of a central region built of up to five core repeats linked to a microsatellite segment at one array end and flanked by sequences holding short inverted repeats. The modular composition and the conserved putative target site duplication sequence AA at the element termini are equivalent to the composition of several elements found in the cupped oyster Crassostrea virginica and in some insects. A unique feature of D. trunculus element is ordered array of core repeat variants, distinctive by diagnostic changes. Position of variants in the array is fixed, regardless of alterations in the core repeat copy number. Each repeat harbors a palindrome near the junction with the following unit, being a potential hotspot responsible for array length variations. As a consequence, variations in number of tandem repeats and variations in flanking sequences make every sequenced element unique. Core repeats may be thus considered as individual units within the MITE, with flanking sequences representing a "cassette" for internal repeats. Our results demonstrate that onset and spread of tandem repeats can be more intimately linked to processes of transposition than previously thought and suggest that genomes are shaped by interplays within a complex network of repetitive sequences.

Molecular and cytogenetic characterization of repetitive DNA in the Antarctic polyplacophoran Nuttallochiton mirandus

Two highly repeated DNAs, designated NmE1/NmE2 and NmE5, were identified by EcoRV digestion in the chiton Nuttallochiton mirandus (Mollusca: Polyplacophora). The comparison of the sequences obtained showed high similarity in 5' and 3' regions and the NmE5 sequence displayed an inserted sequence that might arise from a transposable element. Southern blotting analyses suggested a tandem organization of both satellite DNA families identified. Moreover, dot blot analyses, performed on several molluscan species, revealed a different degree of conservation of the repeated DNAs. Fluorescence in-situ hybridizations (FISH) on metaphase chromosomes showed that both satellite DNAs are located at centromeric regions.

Satellite DNAs between selfishness and functionality: structure, genomics and evolution of tandem repeats in centromeric (hetero)chromatin

Satellite DNAs (tandemly repeated, non-coding DNA sequences) stretch over almost all native centromeres and surrounding pericentromeric heterochromatin. Once considered as inert by-products of genome dynamics in heterochromatic regions, recent studies showed that satellite DNA evolution is interplay of stochastic events and selective pressure. This points to a functional significance of satellite sequences, which in (peri)centromeres may play some fundamental functional roles. First, specific interactions with DNA-binding proteins are proposed to complement sequence-independent epigenetic processes. The second role is achieved through RNAi mechanism, in which transcripts of satellite sequences initialize heterochromatin formation. In addition, satellite DNAs in (peri)centromeric regions affect chromosomal dynamics and genome plasticity. Paradoxically, while centromeric function is conserved through eukaryotes, the profile of satellite DNAs in this region is almost always species-specific. We argue that tandem repeats may be advantageous forms of DNA sequences in (peri)centromeres due to concerted evolution, which maintains high intra-array and intrapopulation sequence homogeneity of satellite arrays, while allowing rapid changes in nucleotide sequence and/or composition of satellite repeats. This feature may be crucial for long-term stability of DNA-protein interactions in centromeric regions.

Initial analysis of tandemly repetitive sequences in the genome of Zhikong scallop (Chlamys farreri Jones et Preston)

Tandemly repetitive sequences are widespread in all eukaryotic genomes, but data on tandem repeats are limited in Zhikong scallop (Chlamys farreri). In the present study, paired-end sequencing of 2016 individual fosmid clones resulted in 3646 sequences. A total of 2,286,986 bp of genomic sequences were generated, representing approximately 1.84 per thousand of the Zhikong scallop genome. Using tandem repeats finder (TRF) software, a total of 2500 tandem repeats were found, including 313 satellites, 1816 minisatellites and 371 microsatellites. The cumulative length of tandem repeats was 552,558 bp, accounting for 24.16% of total length. Specifically, the length of microsatellites, minisatellites and satellites was 9425, 336,001 and 207,132 bp, accounting for 1.71, 60.81 and 37.49% of the length of tandem repeats, and 0.41, 14.69 and 9.06% of total length, respectively. The detailed information on the characteristic of all repeat units was also represented, which will provide a useful resource for physical mapping and better utilization of the existing genomic information in Zhikong scallop.

Repetitive DNA, molecular cytogenetics and genome organization in the King scallop (Pecten maximus)

We studied the structure, organization and relationship of repetitive DNA sequences in the genome of the scallop, Pecten maximus, a bivalve that is important both commercially and in marine ecology. Recombinant DNA libraries were constructed after partial digestion of genomic DNA from scallop with PstI and ApaI restriction enzymes. Clones containing repetitive DNA were selected by hybridisation to labelled DNA from scallop, oyster and mussel; colonies showing strong hybridisation only to scallop were selected for analysis and sequencing. Six non-homologous tandemly repeated sequences were identified in the sequences, and Southern hybridisation with all repeat families to genomic DNA digests showed characteristic ladders of hybridised bands. Three families had monomer lengths around 40 bp while three had repeats characteristic of the length wrapping around one (170 bp), or two (326 bp) nucleosomes. In situ hybridisation to interphase nuclei showed each family had characteristic numbers of clusters indicating contrasting arrangements. Two of the repeats had unusual repetitions of bases within their sequence, which may relate to the nature of microsatellites reported in bivalves. The study of these rapidly evolving sequences is valuable to understand an important source of genomic diversity, has the potential to provide useful markers for population studies and gives a route to identify mechanisms of DNA sequence evolution.

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.

High conservation of the differentially amplified MPA2 satellite DNA family in parthenogenetic root-knot nematodes

Sequence variability and distribution of a newly characterized MPA2 satellite DNA family are described in five root-knot nematode species of the genus Meloidogyne, the mitotic parthenogens M. paranaensis, M. incognita, M. arenaria and M. javanica, and the meiotic/mitotic M. hapla (isolates A and B, respectively). The lack of distinctive mutations and the considerable contribution (40.8%) of ancestral changes disclose an ancient satellite DNA which existed in the common ancestor of extant parthenogenetic species in the same or similar form and remained preserved for a period of at least 43 My. Nonuniformly distributed polymorphic sites along the satellite monomer suggest differences in constraints acting on particular sequence segments. Sequence diversity is clearly unaffected by significant differences in genomic abundance of the MPA2 satellite DNA in the examined species. Observed results suggest that the dynamics of this satellite DNA family might be in the first instance a consequence of characteristics of its nucleotide sequence and possible constraints imposed on it. Under conditions of mitotic and meiotic parthenogenesis, slow accumulation of mutations and slow replacement of old MPA2 sequence variants with new ones may be equivalent to the dynamics of some satellite DNA sequences conserved for extremely long evolutionary periods in sexual species.

Conserved and variable domains in satellite DNAs of mitotic parthenogenetic root-knot nematode species

Two satellite DNAs have been characterized in the mitotic parthenogenetic root-knot nematodes Meloidogyne javanica and M. paranaensis, agriculturally important phytoparasitic species. The satellite repeat variants cloned from M. javanica could not be resolved from those described earlier in M. arenaria [Castagnone-Sereno, P., Leroy, F., Abad, P., 2000. Cloning and characterization of an extremely conserved satellite DNA family from the root-knot nematode Meloidogyne arenaria. Genome 43, 346-353] and are therefore classified as a single satellite named MARJA. However, this satellite shows 34.3% sequence divergence in comparison with the MPA1 satellite characterized in M. paranaensis, and monomer variants of both satellites are clearly distinguished by homogenized nucleotide substitutions. Nucleotide variability analysis revealed in one segment of the satellite monomer domains of high and low variability, conserved both within and between monomer variants of the two satellites. Intersatellite conservation of these domains indicates evolution of satellite sequence under different constraints, probably due to some functional interactions. In addition, high intrasatellite homogeneity, presence of ancestral mutations in groups of MARJA monomers in both M. javanica and M. arenaria and highly homogenized divergent positions in comparison with the MPA1 indicate similar sequence dynamics in mitotic parthenogenetic taxa to that observed in amphimictic species.