#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

De novo identification of satellite DNAs in the sequenced genomes of Drosophila virilis and D. americana using the RepeatExplorer and TAREAN pipelines


Authors: Bráulio S. M. L. Silva aff001;  Pedro Heringer aff001;  Guilherme B. Dias aff001;  Marta Svartman aff001;  Gustavo C. S. Kuhn aff001
Authors place of work: Departamento de Genética, Ecologia e Evolução, Universidade Federal de Minas Gerais, Belo Horizonte, Minas Gerais, Brasil aff001
Published in the journal: PLoS ONE 14(12)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0223466

Summary

Satellite DNAs are among the most abundant repetitive DNAs found in eukaryote genomes, where they participate in a variety of biological roles, from being components of important chromosome structures to gene regulation. Experimental methodologies used before the genomic era were insufficient, too laborious and time-consuming to recover the collection of all satDNAs from a genome. Today, the availability of whole sequenced genomes combined with the development of specific bioinformatic tools are expected to foster the identification of virtually all the “satellitome” of a particular species. While whole genome assemblies are important to obtain a global view of genome organization, most of them are incomplete and lack repetitive regions. We applied short-read sequencing and similarity clustering in order to perform a de novo identification of the most abundant satellite families in two Drosophila species from the virilis group: Drosophila virilis and D. americana, using the Tandem Repeat Analyzer (TAREAN) and RepeatExplorer pipelines. These species were chosen because they have been used as models to understand satDNA biology since the early 70’s. We combined the computational approach with data from the literature and chromosome mapping to obtain an overview of the major tandem repeat sequences of these species. The fact that all of the abundant tandem repeats (TRs) we detected were previously identified in the literature allowed us to evaluate the efficiency of TAREAN in correctly identifying true satDNAs. Our results indicate that raw sequencing reads can be efficiently used to detect satDNAs, but that abundant tandem repeats present in dispersed arrays or associated with transposable elements are frequent false positives. We demonstrate that TAREAN with its parent method RepeatExplorer may be used as resources to detect tandem repeats associated with transposable elements and also to reveal families of dispersed tandem repeats.

Keywords:

Transposable elements – Genome analysis – Drosophila – Drosophila melanogaster – Clustering algorithms – Invertebrate genomics – Tandem repeats – Polytene chromosomes

Introduction

The genome of eukaryotes encloses a variety of repetitive DNA sequences which comprises most of the nuclear DNA of several organisms, including animals, plants and insects [1,2]. Among them are the satellite DNAs (satDNAs), usually defined as abundant, tandemly repeated noncoding DNA sequences, forming large arrays (hundreds of kilobases up to megabases), typically located in the heterochromatic regions of the chromosomes [3,4], although short arrays may additionally be present in the euchromatin [5,6].

The collection of satDNAs in the genome, also known as the “satellitome”, usually represents a significant fraction (>30%) of several animal and plant genomes. Other classes of noncoding tandem repeats include the microsatellites, with repeat units less than 10 bp long, array sizes around 100 bp and scattered throughout the genome; and the minisatellites, with repeats between 10 to 100 bp long, forming up to kb-size arrays, located at several euchromatic regions, with a high density at terminal chromosome regions [3,4]. Therefore, the best criteria to distinguish satellites from micro and minisatellites are long array sizes and preferential accumulation at heterochromatin for the former.

SatDNAs do not encode proteins, but they may play important functional roles in the chromosomes, most notably related to chromatin modulation and the establishment of centromeres [79]. They are among the fastest evolving components of the genome (although some conserved satellites have also been reported) [1012], and such behavior combined to their abundance and structural role have major implications for the evolution and diversification of genomes and species [8,13].

Since the discovery of satDNAs in the early 60’s, species from the genus Drosophila have been used as a model to address several aspects of satDNA biology, such as their origin, organization, variation, evolution and function (e.g. [7,1418]).

Currently, several Drosophila genomes have been sequenced by next-generation technologies and new bioinformatic tools have been designed for the identification of repetitive DNAs from this vast source of genomic resources [19]. Among them, the RepeatExplorer software [20] has been successfully used for de novo identification of repetitive DNAs directly from unassembled short sequence reads, and the recently implemented Tandem Repeat Analyzer (TAREAN) pipeline [21] was introduced to specifically identify putative satDNAs. Such a combination between sequenced genomes and bioinformatic tools is now expected to foster the identification of the full “satellitome” of any given species (e.g. [2226]). Despite the availability of all such resources, only a few Drosophila species had their satDNA landscape determined with these new approaches [23].

In the genus Drosophila genome sizes vary between ~130 Mb to ~400 Mb, but most analyzed species have genome with around 180–200 Mb, such as D. melanogaster [27,28]. The satDNA content also varies across species, from ~2% in D. buzzatii [23] to ~60% in D. nasutoides [29]. Some studies suggest a positive correlation between genome size and the amount of satDNAs in Drosophila [28,30,31].

The genome size of D. virilis (virilis group), with ~400 Mb, is among the largest reported for Drosophila. Accordingly, the estimated satDNA in this species is also high (>40%) [28,32]. Previous studies using CsCl density gradients revealed that three evolutionary related satDNAs with 7 bp long repeat units and only one mutation difference, named satellite1 (5’ ACAAACT 3’), satellite2 (5’ ATAAACT 3’) and satellite3 (5’ ACAAATT 3’) together represent ~40% of its genome [32,33]. These satellites mapped predominantly to the heterochromatic regions of all chromosomes except the Y. Another satDNA identified in this species, but using genomic DNA digestion with restriction endonucleases, was named pvB370, and consists of 370 bp long repeat units [34] predominantly located at sub-telomeric regions and, to a lesser extent, along some discrete euchromatic loci [35]. Other abundant TRs have been identified in the D. virilis genome, such as the 220TR and 154TR families, which belong to the internal structure of transposable elements [16,36], the 225 bp family, present in the intergenic spacer of ribosomal genes, and the less characterized 172 bp family [37]. A recent study reported additional tandem repeats less than 20 bp long but at low abundance [18].

The high throughput and low cost of current whole-genome sequencing technologies have made it possible to obtain genome assemblies for a wide range of organisms. However, de novo whole-genome shotgun strategies are still largely unable to fully recover highly repetitive regions such as centromeres and pericentromeric regions and, as a result, satDNAs are usually misrepresented or absent from such assemblies [19]. One way of circumventing the assembly bottleneck is to directly identify repeats from raw sequencing reads. One of such approaches is implemented in the RepeatExplorer pipeline, already used in a wide range of plant and animal species [22,38,39]. RepeatExplorer performs similarity-based clustering of raw short sequencing reads and partial consensus assembly, allowing for repeat identification even from small samples of genome coverage. A recent development of RepeatExplorer includes the TAREAN pipeline for the specific detection of tandem repeats by searching for circular structures in directed read clusters [21].

In the present study, we aimed to test the ability of TAREAN to correctly identify and estimate the abundance of satDNAs in D. virilis. To refine and expand our knowledge of the identified putative satDNAs, in some cases we mapped them in mitotic and polytene chromosomes using fluorescent in situ hybridization (FISH) technique.

There are several examples showing that satDNA abundance may vary widely even across closely related species [10,40]. For example, one species may present few repeats in the genome (therefore not being identified as a satellite), while a closely related species presents thousands, reaching a satDNA status. For this reason, we also added to our study D. americana, a species belonging from the virilis group, but separated from D. virilis by ~4.1 Myr [41].

Material and methods

RepeatExplorer and TAREAN analyses

The in silico identification of putative satDNAs was performed using the RepeatExplorer and TAREAN pipelines [20,21] implemented in the Galaxy platform [42]. These algorithms were developed to identify and characterize repetitive DNA elements from unassembled short read sequences. We used the publicly available Drosophila virilis strain 160 (SRX669289), Drosophila americana strain H5 (ERX1035147) and Drosophila americana strain W11 (ERX1035149) [43] Illumina paired-end sequences. The sequences were obtained through the “European Nucleotide Archive” (EBI) database and their quality scores measured with the “FASTQC” tool. We used “FASTQ Groomer” (Sanger & Illumina 1.8 +) to convert all the sequences to a single fastqsanger format. We removed adapters and excluded any reads with more than 5% of its sequence in low quality bases (Phred cutoff < 10) using the “Preprocessing of fastq paired-reads” tool included in the RepeatExplorer Galaxy instance. The interlaced filtered paired-end reads were used as input data for the RepeatExplorer clustering and Tandem Repeat Analyzer tools with the following settings: “sample size = 2,000,000—select taxon and protein domain database version (REXdb): Metazoa version 3.0—select queue: extra-long and slow”. For the TAREAN analyses we also used the “perform cluster merging” tool for reducing the redundancy of the results.

The results were provided in a HTML archive report and all the data were downloaded in a single archive for further investigation. We analyzed clusters representing >0.5% of the genome of Drosophila virilis strain 160.

Clusters with tandem repeats identified by TAREAN are denoted as putative high or low confidence satellites. These estimates are denoted according to the “Connected component index (C)” and “Pair completeness index (P)”. The C index indicates clusters formed by tandemly repeated genomic sequences, while the P index measures the ratio between complete read pairs in the cluster and the number of broken pairs, that is directly related to the length of continuous tandem arrays [21].

Fluorescent probe construction

We extracted total genomic DNA from a pool of 20 adult Drosophila virilis (strain 15010–1051.51 from Santiago, Chile) and D. americana (strain H5 from Mississipi, United States of America) with the Wizard® Genomic DNA Purification Kit (Promega Corporation). For primer’s design, we used the consensus sequences from each satDNA identified by RepeatExplorer/TAREAN and multiple sequence alignments by selecting the most conserved nucleotide regions. Satellite DNAs were PCR amplified with the following primers forward (F) and reverse (R):

  • Sat1_F (ACAAACTACAAACTACAAACTACAAACTACAAACT), Sat1_R (AGTTTGTAGTTTGTAGTTTGTAGTTTGTAGTTTGT), 172TR_F (ATTTATGGGCTGGGAAGCTTTGACGTATG), 172TR_R (CGGTCAAATCTCATCCGATTTTCATGAGG), 225TR_F (GCGACACCACTCCCTATATAGG), 225TR_R (CGCGCAAGGCATGTCATATG), pvB370_F (TAGTAGGGATCCGTACAAATTCAA), pvB370_R (GTACGGATCCCTACTAATAATTGGCAT).

All primers were used to amplify the target sequences from genomic DNA, with the exception of Sat1 in which the amplification process was conducted by forward and reverse primers self-annealing without genomic DNA. The PCR products were excised from agarose gels and ligated into pGEM-T vector plasmids (Promega) with T4 DNA ligase (Promega). For cloning, the plasmids were multiplied into E.coli cells and then eluted with the PureLink Quick Plasmid Miniprep Kit (Invitrogen). To ensure the presence of the inserts, the final samples were Sanger sequenced in an ABI3130 and later analyzed in the Chromas software (Technelysium). Clones with satDNA inserts were later prepared as probes for FISH.

Fluorescent in situ hybridization (FISH)

The metaphase and polytene chromosomes were obtained from neuroblasts and salivary glands of third instar larvae of D. virilis (strain 15010–1051.51) and D. americana (strain H5), according to [44,45]. Probe labeling and FISH experiment conditions were conducted according to [16]. The satDNA probes were immunodetected with antidigoxigenin-Rhodamine and avidin-FITC (Roche Applied Science).

We used DAPI “4,6-diamidino-2-phenylindole” (Roche) in “SlowFade” antifade reagent (Invitrogen) for DNA counterstaining. The analyses were conducted under an Axio Imager A2 epifluorescence microscope equipped with the AxiocamMRm camera (Zeiss). Images were captured with Axiovision (Zeiss) and edited in Adobe Photoshop.

Results

Identification of putative satDNAs in D. virilis and D. americana

The most abundant putative satDNAs (covering >0.5% of the genome) identified by the RepeatExplorer and TAREAN pipelines are shown in Table 1 (see S1 Fig for histogram summary analyses and S4S15 Figs for detailed data from each cluster retrieved). All of the six identified tandem repeat families (Sat1, 154TR, pvB370, 172TR, 225TR, 36TR) are shared by both species and have been previously identified.

Tab. 1. Putative satellite DNAs in D. virilis strain 160 and D. americana strain H5 identified by TAREAN and Repeat Explorer.
Putative satellite DNAs in <i>D</i>. <i>virilis</i> strain 160 and <i>D</i>. <i>americana</i> strain H5 identified by TAREAN and Repeat Explorer.

Although the total abundance of these six tandem repeats is similar (~17%) in the two species, there are differences in the estimated proportion occupied by each putative satDNA between the species.

In order to check if these differences are predominantly inter-specific, we used RepeatExplorer and TAREAN to compare the abundances of each tandem repeat between two D. americana strains (H5 and W11), which were sequenced using the same sequencing platform and methods. Our analysis showed that differences in repeat proportion among D. americana strains are somewhat comparable with the ones observed between D. virilis and D. americana (Fig 1). These results indicate that comparisons of tandem repeat abundance between taxa using RepeatExplorer and TAREAN should be taken with caution as significant differences can also be observed among lineages within the same species. Interestingly, repeat abundance variations between lineages in these species were also detected by [46].

Genome content for six putative satellite DNAs in two <i>Drosophila americana</i> strains (H5 and W11) according to RepeatExplorer and TAREAN analyses.
Fig. 1. Genome content for six putative satellite DNAs in two Drosophila americana strains (H5 and W11) according to RepeatExplorer and TAREAN analyses.

To further characterize the tandem repeat families identified in silico, we constructed DNA probes using the consensus sequences generated by TAREAN from three families and used them to verify their localization in metaphase and polytene chromosomes. In the following sections we describe our in silico and FISH analyses for each identified family, comparing the results with previous studies and discussing if TAREAN correctly identified and distinguished satDNAs from other classes of tandem repeats. The tandem repeat families are described below in order of their abundance (higher to lower) as revealed for D. virilis strain 160.

Sat1

The most abundant tandem repeat identified by TAREAN in D. virilis and D. americana is composed by a 7 bp long repeat corresponding to the previous described satellite I [33]. In D. virilis, our FISH experiments in metaphase chromosomes showed this satDNA occupying the pericentromeric region of all autosomes except the small dot chromosomes, and in the X and Y chromosomes (Fig 2A). However, the hybridization in polytene chromosomes revealed that Sat1 also localizes in the pericentromeric region of the dot chromosome (Fig 3A). [32] showed a similar hybridization pattern, although their results did not consistently demonstrate Sat1 signals in the dot and Y chromosomes.

Mapping of Sat1, 172TR and 225TR by FISH on metaphase chromosomes.
Fig. 2. Mapping of Sat1, 172TR and 225TR by FISH on metaphase chromosomes.
(A) Drosophila virilis and (B) Drosophila americana. Upper panel: Sat1 (green) and 172TR (red). Lower panel: 225TR (red). The mitotic chromosomes of D. virilis were identified by their sizes combined with the hybridization signals on polytene chromosomes (see Fig 3).
Mapping of Sat1, 172TR and 225TR by FISH on polytene chromosomes.
Fig. 3. Mapping of Sat1, 172TR and 225TR by FISH on polytene chromosomes.
(A, C) Drosophila virilis and (B, D) Drosophila americana. (A, B) Sat1 (green) and 172TR (red). (C, D) 225TR (red). Scale bars represent 10μm. The polytene chromosome arms were identified according to [47].

In D. americana, Sat1 signals were detected in the pericentromeric region of all autosomes in metaphase chromosomes, except the dot (Fig 2B), while in polytene chromosomes, Sat1 signals were also observed in the dot chromosomes (Fig 3B). However, differently to what was observed in D. virilis, our Sat1 hybridizations in the D. americana polytene dot chromosomes did not give enough information about the precise location of this satDNA, although it also appears to occupy a portion of the pericentromeric region. As another difference from D. virilis, Sat1 sequences appear to be absent from the Y chromosome in D. americana (Fig 2B). Our FISH results corroborate the smaller genomic fraction occupied by this satDNA in D. americana (~9% against ~12% in D. virilis), revealed by the in silico analysis (Figs 2, 3A and 3B). These new findings in D. americana and D. virilis also agree with recent results from [46].

154TR

The genomic distribution of 154TR has been recently studied in detail in D. virilis and D. americana using FISH in metaphase and polytene chromosomes [36]. This sequence was independently identified in silico by [37] and [48]. The 154TR was characterized as a tandem repeat derived from a Helitron transposable element [37], which was studied in detail and classified as a family named DINE-TR1 [36]. DINE-TR1 elements containing 154TR homologous sequences were found in several Acalyptratae species, mostly within the Drosophila genus, although long arrays (> 10 copies) of 154TR were only detected in three species (D. virilis, D. americana and D. biarmipes) [36].

FISH in metaphase and polytene chromosomes revealed that 154TR is located in the distal pericentromeric region (β-heterochromatin) and many euchromatic loci of all autosomes and the X chromosome of D. virilis and D. americana. In addition, this tandem repeat covers a large portion of the Y chromosome in both species. In D. virilis, 154TR signals are very abundant in the centromeric heterochromatin of chromosome 5 and are also found in a discrete region within the pericentromeric region of the X chromosome [36].

Our results from the TAREAN analysis classified 154TR as a putative satellite with low confidence in both species (Table 1). We suggest that this result is probably a consequence of 154TR being both tandemly repeated, like a satDNA, and dispersed, like a transposable element. In this case, even though the connected component index (C) of 154TR is high, its relatively low pair completeness index (P) contributes to its classification as a putative satellite with low confidence by TAREAN (Table 1). We suggest that 154TR is not a satDNA and thus, should be classified as a highly abundant dispersed tandem repeat.

pvB370

The pvB370 satellite was first described by [34], who also identified this family as deriving from the direct terminal repeats of pDv transposable elements [49]. In a following study, [35] showed that in D. virilis and D. americana, pvB370 is located at several euchromatic loci and at the telomeric region of all chromosomes.

Because pvB370 was previously mapped in the chromosomes of D. virilis and D. americana using FISH, we did not conduct a throughout analysis on both species. However, because pvB370 seems to display a euchromatic distribution [35] similar to the one we observed for 172TR (Fig 3A and 3B) we hybridized both pvB370 and 172TR probes concomitantly in D. americana polytene chromosomes. Our results showed little or no overlap between pvB370 and 172TR, although many arrays from the two families are very close (at least a few kbp) to each other (Fig 4).

Chromosome location of 172TR and pvB370 by FISH on polytene chromosomes of <i>Drosophila americana</i>.
Fig. 4. Chromosome location of 172TR and pvB370 by FISH on polytene chromosomes of Drosophila americana.
There is little or no overlap between these tandem repeats. Red (172TR) and green (pvB370). Scale bar represents 10μm.

172TR

The 172TR family corresponds to the 172 bp tandem repeats previously identified in silico by [37]. Our FISH results in the metaphase and polytene chromosomes of D. virilis revealed that 172TR is distributed throughout the arms of autosomes 3, 4 and 5, in several loci at the X chromosome and in at least two loci in chromosome 2, including the subtelomeric region (Figs 2A and 3A). Most of the arrays are located at distal chromosome regions. No hybridization signals were detected in the dot and Y chromosomes.

The FISH results in D. americana showed 172TR signals at multiple loci along all autosomes, except the dot, and more equally distributed in both distal and proximal regions of chromosome arms (Figs 2B and 3B). Similarly to D. virilis, no hybridization signal was detected in the Y chromosome (Fig 2B). The FISH data (Figs 2, 3A and 3B) clearly showed a higher number of 172TR loci in D. americana compared to D. virilis, a result that is consistent with the higher overall abundance of 172TR repeats in D. americana predicted by the in silico analysis (Table 1).

225TR

The putative satDNA detected in our in silico analyses as 225TR was previously identified as a component of intergenic spacers (IGS) of ribosomal genes from D. virilis located at the chromocenter and nucleolus regions of polytene chromosomes [37]. Our FISH experiments in polytene chromosomes confirmed these results in D. virilis (Fig 3C), additionally showing that in D. americana this family displays the same pattern of localization (Fig 3D).

In addition, we also performed FISH with a 225TR probe in metaphase chromosomes of both species for the first time, that revealed its location in the pericentromeric region of the X chromosome and in the pericentromeric and telomeric regions of chromosome Y (Fig 2A and 2B). This result is in accordance with previous studies showing the location of these IGS sequences in the sex chromosomes of Drosophila [50].

Although the TAREAN pipeline failed to detect the 225TR in D. americana, RepeatExplorer revealed the presence of this family. This indicates a possible limitation of TAREAN in detecting less abundant tandem repeats in comparison with RepeatExplorer. Moreover, TAREAN only retrieves clusters with highly circular structures, and therefore excludes 225TR repeats that are associated with linear structures (S12 and S13 Figs). These observations indicate that, although 225TR is an abundant tandem repeat, it does not have all the typical features of a satDNA.

36TR

A previous work made by [49] identified the presence of 36 bp tandem repeats inside the pDv transposable element and a subsequent work by [34] showed that array size variation exists among different pDv copies in D. virilis. This TR was not retrieved by the TAREAN pipeline but we found it in high abundance (~0.73% in D. virilis and ~0.48% in D. americana; Table 1) among the results from RepeatExplorer, that further classified this TR as a low confidence satDNA. Interestingly, the RepeatExplorer pipeline revealed that the cluster corresponding to this 36 bp tandem repeat has a high number of shared reads with the pvB370 cluster (S2 and S3 Figs). In this case, the link between 36TR and pvB370 clusters is explained by their co-occurrence as complete (36 bp) and partial (pvB370) sequences within the pDv transposable element [34]. This result shows that the RepeatExplorer pipeline is able to detect putative relationships between distinct repetitive sequences.

Discussion

Here we performed de novo identification of the most abundant tandem repeat families in D. virilis and D. americana. These species were chosen because they have larger genomes compared to other Drosophila species and because they have been used as models to understand satDNA biology since the early 70’s. In order to do that, we combined the RepeatExplorer and TAREAN results with data from the literature and, in some cases, with new chromosome mapping data obtained by us using FISH in metaphase and polytene chromosomes.

Because all of the repeats identified herein had been previously detected by other methods, we were able to test if the TAREAN pipeline could correctly classify them as satDNAs or not.

TAREAN identified the heptanucleotide Sat1 as a satDNA with high confidence, which agrees with all attributes known for this family and the satDNA definition (i.e. high copy-number, long-arrays, predominant heterochromatic location) [32,33]. Sat1 was identified as the most abundant tandem repeat in both D. virilis and D. americana, which is also in accordance with previous work [32,33]. However, the other two less abundant heptanucleotide satellites, Sat2 and Sat3, were not detected by TAREAN. As these three satellites differ from each other by a single nucleotide substitution, they were likely all included in the Sat1 cluster by TAREAN. This clustering of variants appears to be a relevant disadvantage that might influence the identification of not only the heptanucleotide satDNA family but other short repeat families with similar features (e.g. short monomer size and high sequence similarity). Therefore, to analyze these type of sequences in detail, it might be advisable to also use tools that are more appropriate for this aim, for example, the software k-Seek [51]. It is also worth mentioning that the heptanucleotide satDNA genomic fractions revealed by TAREAN (~12% for D. virilis strain 160 and ~9% for D. americana strain H5) are significantly below the previously estimated of >40% genomic fraction, based on density gradient ultracentrifugation methods [32,52]. Although TAREAN may not be ideally suitable to quantify satellites with short repeat units [21], it is worth mentioning that [46] have recently demonstrated that Illumina sequence reads containing the heptanucleotide satellites from D. virilis tend to be highly enriched for low quality scores. Furthermore, the use of raw reads from different sequencing platforms did not allowed the recovery of simple satellites at the predicted ~40% genomic fraction indicated by previous works [46]. The difference between these estimates (12% to 40%) may reflect an intrinsic bias in current sequencing methods. A second possibility, which does not reject the first is the existence of real differences in satDNA content between different strains of the same species.

TAREAN classified the 154TR, pvB370 and 36TR families as putative satellites in D. virilis and D. americana. With the exception of pvB370 in D. americana, which was classified with high confidence, all remaining repeats had low confidence calls from TAREAN (Table 1). These tandem repeats are known to be abundant and associated with transposable elements (as integral parts or evolutionarily related), suggesting that RepeatExplorer and TAREAN could be used as resources to detect tandem repeats associated with transposable elements. In the case of 154TR, pvB370 and 36TR, the relationship could be checked directly in the RepeatExplorer pipeline by identifying clusters of tandem repeats sharing a high number of reads with clusters associated to transposable elements (see S2 and S3 Figs), or indirectly in the TAREAN pipeline, by investigating the tandem repeats classified as putative satellites with low confidence (or lower values of satellite probability). The rationale behind this last procedure is that identified families with a ‘low satellite score’ may represent repetitive DNAs with intermediate features, being both highly dispersed and tandemly repeated. One situation in which this scenario is expected is the case where tandem repeats belonging to the terminal or internal portions of transposable elements underwent array expansion [36,53]. Nonetheless, some highly dispersed tandem repeats are not necessarily associated with transposable elements, which is the case of 172TR shown here and the 1.688 satDNA from D. melanogaster [5].

It is interesting to note that, in D. virilis and D. americana, the families 172TR, pvB370 and 154TR were either classified as putative satellites with low confidence, or with high confidence but associated with a relatively low satellite probability (Table 1). Because all these three families were found distributed along the euchromatic regions of chromosomes, we suggest that a low 'satellite score' in the TAREAN pipeline is a good predictor of dispersed tandem repeats. As mentioned above, although there is no indication of a relationship between the 172TR family with any known transposable element, its lower satellite score from the in silico analysis correctly predicts the dispersed array distribution observed in polytene chromosomes (Fig 3A and 3B).

In conclusion, six abundant putative satDNAs were identified in D. virilis and D. americana by TAREAN and RepeatExplorer: Sat1, 154TR, pvB370, 172TR, 225TR and 36TR. All of them have been previously characterized to a higher or lesser extent in previous works, but using different methodologies. The main advantage of TAREAN and RepeatExplorer in comparison with previous methods aiming to identify satDNAs in D. virilis refers to their relative lack of bias compared to the in silico digestion applied by [37], that identifies only tandem repeats presenting restriction sites, and the k-Seek method [51] applied by [18] that specifically identifies short tandem repeats with less than 20 bp.

While Sat1 (identified by TAREAN as a satDNA with high confidence) is in fact a family that matches all features typically attributed for satDNAs, the classification of the other families as satDNAs (identified as a satDNA with low confidence on at least one species) is more controversial. The 154TR, pvB370 and 36TR families are associated with the internal structure of TEs, thus being distributed along the chromosome arms with different degrees of dispersion. The 225TR belongs to the IGS of ribosomal genes. In contrast, the 172TR family is an abundant tandem repeat but with exclusive euchromatic location, where they apparently do not to reach satDNA-like long arrays. Based on the repeat unit length of 172TR (172 bp), this family cannot be considered as a micro or minisatellite. In this context, it would be interesting to further investigate these five families (154TR, pvB370, 172TR, 225TR and 36TR) using long-read sequencing technologies, since they are expected to provide more detailed information about their copy number and array sizes.

Supporting information

S1 Fig [pdf]
TAREAN histogram summary analyses of (A) (strain 160) and (B) (strain H5).

S2 Fig [pdf]
pvB370 and 36TR supercluster analysis in strain 160.

S3 Fig [pdf]
pvB370 and 36TR supercluster analysis in strain H5.

S4 Fig [pdf]
Sat1 cluster analysis in strain 160.

S5 Fig [pdf]
Sat1 cluster analysis in strain H5.

S6 Fig [pdf]
154TR cluster analysis in strain 160.

S7 Fig [pdf]
154TR cluster analysis in strain H5.

S8 Fig [pdf]
pvB370 cluster analysis in strain 160.

S9 Fig [pdf]
pvB370 cluster analysis in strain H5.

S10 Fig [pdf]
172TR cluster analysis in strain 160.

S11 Fig [pdf]
172TR cluster analysis in strain H5.

S12 Fig [pdf]
225TR cluster analysis in strain 160.

S13 Fig [pdf]
225TR cluster analysis in strain H5.

S14 Fig [pdf]
36TR cluster analysis in strain 160.

S15 Fig [pdf]
36TR cluster analysis in Drosophila americana strain H5.


Zdroje

1. de Koning APJ, Gu WJ, Castoe TA, Batzer MA, Pollock DD. Repetitive elements may comprise over two thirds of the human genome. PLoS genetics. 2011;7(12). https://doi.org/10.1371/journal.pgen.1002384 22144907

2. Biscotti MA, Olmo E, Heslop-Harrison JS. Repetitive DNA in eukaryotic genomes. Chromosome Research. 2015;23(3):415–20. https://doi.org/10.1007/s10577-015-9499-z 26514350

3. Tautz D. Notes on the fefinition and nomenclature of tandemly repetitive DNA sequences. Exs. 1993;67:21–8. https://doi.org/10.1007/978-3-0348-8583-6_2 8400689

4. Charlesworth B, Sniegowski P, Stephan W. The evolutionary dynamics of repetitive DNA in eukaryotes. Nature. 1994;371(6494):215–20. https://doi.org/10.1038/371215a0 8078581

5. Kuhn GCS, Kuttler H, Moreira O, Heslop-Harrison JS. The 1.688 repetitive DNA of Drosophila: concerted evolution at different genomic scales and association with genes. Mol Biol Evol. 2012;29(1):7–11. https://doi.org/10.1093/molbev/msr173 21712468

6. Pavlek M, Gelfand Y, Plohl M, Mestrovic N. Genome-wide analysis of tandem repeats in Tribolium castaneum genome reveals abundant and highly dynamic tandem repeat families with satellite DNA features in euchromatic chromosomal arms. DNA Research. 2015;22(6):387–401. https://doi.org/10.1093/dnares/dsv021 26428853

7. Rosic S, Kohler F, Erhardt S. Repetitive centromeric satellite RNA is essential for kinetochore formation and cell division (vol 207, pg 335, 2014). J Cell Biol. 2014;207(5):673-. https://doi.org/10.1083/jcb.201404097 25365994

8. Kursel LE, Malik HS. The cellular mechanisms and consequences of centromere drive. Curr Opin Cell Biol. 2018;52:58–65. https://doi.org/10.1016/j.ceb.2018.01.011 29454259

9. Bracewell R, Chatla K, Nalley MJ, Bachtrog D. Dynamic turnover of centromeres drives karyotype evolution in Drosophila. BioRxiv [PrePrint]. 2019:733527. [posted 2019 Aug 27] https://www.biorxiv.org/content/10.1101/733527v1.full. https://doi.org/10.1101/733527 31524597

10. Kuhn GCS, Sene FM, Moreira-Filho O, Schwarzacher T, Heslop-Harrison JS. Sequence analysis, chromosomal distribution and long-range organization show that rapid turnover of new and old pBuM satellite DNA repeats leads to different patterns of variation in seven species of the Drosophila buzzatii cluster. Chromosome Research. 2008;16(2):307–24. https://doi.org/10.1007/s10577-007-1195-1 18266060

11. Plohl M, Meštrović N, Mravinac B. Satellite DNA evolution. Repetitive DNA: Karger Publishers; 2012. p. 126–52. https://doi.org/10.1159/000337122

12. Garrido-Ramos MA. Satellite DNA: an evolving topic. Genes-Basel. 2017;8(9). https://doi.org/10.3390/genes8090230 28926993

13. Ferree PM, Barbash DA. Species-specific heterochromatin prevents mitotic chromosome segregation to cause hybrid lethality in Drosophila. Plos Biol. 2009;7(10). https://doi.org/10.1371/journal.pbio.1000234 19859525

14. Strachan T, Webb D, Dover GA. Transition stages of molecular drive in multiple-copy DNA families in Drosophila. Embo J. 1985;4(7):1701–8. https://doi.org/10.1002/j.1460-2075.1985.tb03839.x 16453627

15. Bachmann L, Sperlich D. Gradual evolution of a specific satellite DNA family in Drosophila ambigua, D. tristis, and D. obscura. Mol Biol Evol. 1993;10(3):647–59. https://doi.org/10.1093/oxfordjournals.molbev.a040029 8336547

16. Dias GB, Svartman M, Delprat A, Ruiz A, Kuhn GCS. Tetris is a foldback transposon that provided the building blocks for an emerging satellite DNA of Drosophila virilis. Genome Biol Evol. 2014;6(6):1302–13. https://doi.org/10.1093/gbe/evu108 24858539

17. Khost DE, Eickbush DG, Larracuente AM. Single-molecule sequencing resolves the detailed structure of complex satellite DNA loci in Drosophila melanogaster. Genome research. 2017;27(5):709–21. https://doi.org/10.1101/gr.213512.116 28373483

18. Wei KHC, Lower SE, Caldas IV, Sless TJS, Barbash DA, Clark AG. Variable rates of simple satellite gains across the Drosophila phylogeny. Mol Biol Evol. 2018;35(4):925–41. https://doi.org/10.1093/molbev/msy005 29361128

19. Treangen TJ, Salzberg SL. Repetitive DNA and next-generation sequencing: computational challenges and solutions. Nature Reviews Genetics. 2012;13(1):36. https://doi.org/10.1038/nrg3117 22124482

20. Novak P, Neumann P, Pech J, Steinhaisl J, Macas J. RepeatExplorer: a Galaxy-based web server for genome-wide characterization of eukaryotic repetitive elements from next-generation sequence reads. Bioinformatics. 2013;29(6):792–3. https://doi.org/10.1093/bioinformatics/btt054 23376349

21. Novak P, Robledillo LA, Koblizkova A, Vrbova I, Neumann P, Macas J. TAREAN: a computational tool for identification and characterization of satellite DNA from unassembled short reads. Nucleic acids research. 2017;45(12). https://doi.org/10.1093/nar/gkx257 28402514

22. Ruiz-Ruano FJ, Lopez-Leon MD, Cabrero J, Camacho JPM. High-throughput analysis of the satellitome illuminates satellite DNA evolution. Scientific reports. 2016;6. https://doi.org/10.1038/srep28333 27385065

23. de Lima LG, Svartman M, Kuhn GCS. Dissecting the satellite DNA landscape in three cactophilic Drosophila sequenced genomes. G3-Genes Genom Genet. 2017;7(8):2831–43. https://doi.org/10.1534/g3.117.042093 28659292

24. Palacios-Gimenez OM, Dias GB, de Lima LG, Kuhn GCES, Ramos E, Martins C, et al. High-throughput analysis of the satellitome revealed enormous diversity of satellite DNAs in the neo-Y chromosome of the cricket Eneoptera surinamensis. Scientific reports. 2017;7. https://doi.org/10.1038/s41598-017-06822-8 28743997

25. Utsunomia R, Silva DMZD, Ruiz-Ruano FJ, Goes CAG, Melo S, Ramos LPE, et al. Satellitome landscape analysis of Megaleporinus macrocephalus (Teleostei, Anostomidae) reveals intense accumulation of satellite sequences on the heteromorphic sex chromosome. Scientific reports. 2019;9. https://doi.org/10.1038/s41598-019-42383-8 30971780

26. Liu Q, Li XY, Zhou XY, Li MZ, Zhang FJ, Schwarzacher T, et al. The repetitive DNA landscape in Avena (Poaceae): chromosome and genome evolution defined by major repeat classes in whole-genome sequence reads. Bmc Plant Biol. 2019;19. https://doi.org/10.1186/s12870-019-1769-z 31146681

27. Adams MD, Celniker SE, Holt RA, Evans CA, Gocayne JD, Amanatides PG, et al. The genome sequence of Drosophila melanogaster. Science. 2000;287(5461):2185–95. Epub 2000/03/25. https://doi.org/10.1126/science.287.5461.2185 10731132

28. Bosco G, Campbell P, Leiva-Neto JT, Markow TA. Analysis of Drosophila species genome size and satellite DNA content reveals significant differences among strains as well as between species. Genetics. 2007;177(3):1277–90. https://doi.org/10.1534/genetics.107.075069 18039867

29. Miklos G. Localized highly repetitive DNA sequences in vertebrate and invertebrate genomes. Molecular evolutionary genetics. 1985:241–321.

30. Gregory TR, Johnston JS. Genome size diversity in the family Drosophilidae. Heredity. 2008;101(3):228–38. https://doi.org/10.1038/hdy.2008.49 18523443

31. Craddock EM, Gall JG, Jonas M. Hawaiian Drosophila genomes: size variation and evolutionary expansions. Genetica. 2016;144(1):107–24. Epub 2016/01/23. https://doi.org/10.1007/s10709-016-9882-5 26790663

32. Gall JG, Cohen EH, Polan ML. Repetitive DNA sequences in Drosophila. Chromosoma. 1971;33(3):319-+. https://doi.org/10.1007/BF00284948 5088497

33. Gall JG, Atherton DD. Satellite DNA sequences in Drosophila virilis. J Mol Biol. 1974;85(4):633–64. https://doi.org/10.1016/0022-2836(74)90321-0 4854195

34. Heikkinen E, Launonen V, Muller E, Bachmann L. The pvB370 BamHI satellite DNA family of the Drosophila virilis group and its evolutionary relation to mobile dispersed genetic pDv elements. Journal of molecular evolution. 1995;41(5):604–14. https://doi.org/10.1007/BF00175819 7490775

35. Biessmann H, Zurovcova M, Yao JG, Lozovskaya E, Walter MF. A telomeric satellite in Drosophila virilis and its sibling species. Chromosoma. 2000;109(6):372–80. https://doi.org/10.1007/s004120000094 11072792

36. Dias GB, Heringer P, Svartman M, Kuhn GC. Helitrons shaping the genomic architecture of Drosophila: enrichment of DINE-TR1 in alpha and beta-heterochromatin, satellite DNA emergence, and piRNA expression. Chromosome research: an international journal on the molecular, supramolecular and evolutionary aspects of chromosome biology. 2015;23(3):597–613. Epub 2015/09/27. https://doi.org/10.1007/s10577-015-9480-x 26408292

37. Abdurashitov MA, Gonchar DA, Chernukhin VA, Tomilov VN, Tomilova JE, Schostak NG, et al. Medium-sized tandem repeats represent an abundant component of the Drosophila virilis genome. BMC genomics. 2013;14:771. Epub 2013/11/12. https://doi.org/10.1186/1471-2164-14-771 24209985

38. Garcia G, Rios N, Gutierrez V. Next-generation sequencing detects repetitive elements expansion in giant genomes of annual killifish genus Austrolebias (Cyprinodontiformes, Rivulidae). Genetica. 2015;143(3):353–60. Epub 2015/03/21. https://doi.org/10.1007/s10709-015-9834-5 25792372

39. Robledillo LÁ, Koblížková A, Novák P, Böttinger K, Vrbová I, Neumann P, et al. Satellite DNA in Vicia faba is characterized by remarkable diversity in its sequence composition, association with centromeres, and replication timing. Scientific reports. 2018;8(1):5838. https://doi.org/10.1038/s41598-018-24196-3 29643436

40. Ugarkovic D, Plohl M. Variation in satellite DNA profiles—causes and effects. Embo J. 2002;21(22):5955–9. https://doi.org/10.1093/emboj/cdf612 12426367

41. Morales-Hojas R, Reis M, Vieira CP, Vieira J. Resolving the phylogenetic relationships and evolutionary history of the Drosophila virilis group using multilocus data. Mol Phylogenet Evol. 2011;60(2):249–58. https://doi.org/10.1016/j.ympev.2011.04.022 21571080

42. Afgan E, Baker D, van den Beek M, Blankenberg D, Bouvier D, Cech M, et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2016 update. Nucleic acids research. 2016;44(W1):W3–W10. https://doi.org/10.1093/nar/gkw343 27137889

43. Fonseca NA, Morales-Hojas R, Reis M, Rocha H, Vieira CP, Nolte V, et al. Drosophila americana as a model species for comparative studies on the molecular basis of phenotypic variation. Genome Biol Evol. 2013;5(4):661–79. https://doi.org/10.1093/gbe/evt037 23493635

44. Baimai V. Chromosomal Polymorphisms of Constitutive Heterochromatin and inversions in Drosophila. Genetics. 1977;85(1):85–93. 838273

45. Ashburner M. Drosophila. A laboratory handbook: Cold spring harbor laboratory press; 1989. ISBN: 0879693215

46. Flynn JM, Long M, Wing RA, Clark AG. Evolutionary dynamics of abundant 7 bp satellites in the genome of Drosophila virilis. BioRxiv [PrePrint]. 2019:693077 [posted 2019 July 4] https://www.biorxiv.org/content/10.1101/693077v1.full. https://doi.org/10.1101/693077

47. Gubenko IS, Evgenev MB. Cytological and linkage maps of Drosophila virilis chromosomes. Genetica. 1984;65(2):127–39. https://doi.org/10.1007/BF00135277.

48. Melters DP, Bradnam KR, Young HA, Telis N, May MR, Ruby JG, et al. Comparative analysis of tandem repeats from hundreds of species reveals unique insights into centromere evolution. Genome biology. 2013;14(1). https://doi.org/10.1186/gb-2013-14-1-r10 23363705

49. Zelentsova ES, Vashakidze RP, Krayev AS, Evgenev MB. Dispersed repeats in Drosophila virilis: elements mobilized by interspecific hybridization. Chromosoma. 1986;93(6):469–76. https://doi.org/10.1007/BF00386786

50. Roy V, Monti-Dedieu L, Chaminade N, Siljak-Yakovlev S, Aulard S, Lemeunier F, et al. Evolution of the chromosomal location of rDNA genes in two Drosophila species subgroups: ananassae and melanogaster. Heredity. 2005;94(4):388. https://doi.org/10.1038/sj.hdy.6800612 15726113

51. Wei KH, Grenier JK, Barbash DA, Clark AG. Correlated variation and population differentiation in satellite DNA abundance among lines of Drosophila melanogaster. Proc Natl Acad Sci U S A. 2014;111(52):18793–8. https://doi.org/10.1073/pnas.1421951112 25512552

52. Cohen EH, Bowman SC. Detection and location of three simple sequence DNAs in polytene chromosomes from virilis group species of Drosophila. Chromosoma. 1979;73(3):327–55. Epub 1979/08/01. https://doi.org/10.1007/BF00288696 510073

53. Mestrovic N, Mravinac B, Pavlek M, Vojvoda-Zeljko T, Satovic E, Plohl M. Structural and functional liaisons between transposable elements and satellite DNAs. Chromosome Research. 2015;23(3):583–96. https://doi.org/10.1007/s10577-015-9483-7 26293606


Článek vyšel v časopise

PLOS One


2019 Číslo 12
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#