- Research article
- Open Access
The repetitive component of the sunflower genome as shown by different procedures for assembling next generation sequencing reads
BMC Genomics volume 14, Article number: 686 (2013)
Next generation sequencing provides a powerful tool to study genome structure in species whose genomes are far from being completely sequenced. In this work we describe and compare different computational approaches to evaluate the repetitive component of the genome of sunflower, by using medium/low coverage Illumina or 454 libraries.
By varying sequencing technology (Illumina or 454), coverage (0.55 x-1.25 x), assemblers and assembly procedures, six different genomic databases were produced. The annotation of these databases showed that they were composed of different proportions of repetitive DNA families. The final assembly of the sequences belonging to the six databases produced a whole genome set of 283,800 contigs. The redundancy of each contig was estimated by mapping the whole genome set with a large Illumina read set and measuring the number of matched Illumina reads. The repetitive component amounted to 81% of the sunflower genome, that is composed mainly of numerous families of Gypsy and Copia retrotransposons. Also many families of non autonomous retrotransposons and DNA transposons (especially of the Helitron superfamily) were identified.
The results substantially matched those previously obtained by using a Sanger-sequenced shotgun library and a standard 454 whole-genome-shotgun approach, indicating the reliability of the proposed procedures also for other species. The repetitive sequences were collected to produce a database, SUNREP, that will be useful for the annotation of the sunflower genome sequence and for studying the genome evolution in dicotyledons.
Eukaryotic species show considerable variation in genome size. This is especially true in higher plants, whose genome size (1C) ranges from 63 Mbp in Genlisea margaretae to 150 Gbp in Paris japonica[1, 2]. Such differences have evolved mainly because of two processes: polyploidy and DNA amplification of transposons and related sequences. In eukaryotic genomes, the latter process has resulted in the accumulation of many repeated sequences – sequences that are similar or identical to sequences elsewhere in the genome, but whose number of copies is much higher than that possibly achieved through polyploidization. Differences in genome size among species largely depend on the size of this repetitive fraction. In fact, large genomes are filled with repetitive sequences, especially in plants . Some repeats appear to be non-functional, whereas others have played key roles in the evolution of species . For example, the mutagenic action of transposons provides substantial increases in genetic variability . Transposons also create novel functions, and alter the regulatory patterns of genes, resulting in phenotypic variation [5, 6].
The advent of next generation sequencing (NGS) represents a major advance for genetical and biological research, producing millions of genomic sequences at ever increasing speed and decreasing cost. Dozens of Gigabases of data can be sequenced in a week for the same cost as a few hundred kilobases of Sanger sequence (, updated at http://www.molecularecologist.com/next-gen-fieldguide-2013/). NGS technology has offered the opportunity to acquire genome-scale data for any organism [8–10].
In either reference guided or de novo assembly of NGS reads, a major computational task is to manage 'multi-reads', i.e. those reads that map to multiple locations and/or contain highly repeated k-mers [11, 12]. An algorithm for reference-guided assembly has three choices : 1) to ignore (hence discard) all multi-reads; 2) to perform the best match approach, in which only the best alignment is reported or, if equally good best match alignments occur, one at random or all of them are reported; 3) to report all alignments up to a maximum number. The first strategy restricts the analysis to unique regions in the genome, by discarding all repeats and limiting discovery of some biologically important variants. The other two strategies enable analyses of repetitive regions, with the best match approach providing a reasonable estimate of coverage  and reporting all possible alignments to avoid erroneous choices about read placement.
De novo assemblers belong to one of two classes, overlap-based and de Bruijn graph assemblers, that each create different types of graphs from the read data. The sequence assembly is then reconstructed by algorithms that traverse the graphs. Repeats cause branches in these graphs  and assemblers, making a guess as to which branch to follow, can create false joins and erroneous copy numbers. In a more conservative approach, the assembler breaks graphs at these branch points, producing an accurate but fragmented assembly.
The most common error of an assembler is the production of a chimaera by joining two repeats that are not close in the genome. To resolve chimaeras the first and most important tool is the use of paired-end reads. Because the distance between the paired reads is known, an assembler can use both the expected distance and the orientation of the reads to reconstruct the correct sequence.
Another strategy for handling repeats is to perform statistics on the depth of coverage for each contig. These statistics cannot show exactly how to assemble every repeat, but they do identify the repeats themselves. The assumption is that if a genome is sequenced, for example, to 50x coverage, the genome should be uniformly covered. This means that most contigs should also be covered at 50x. By contrast, a repetitive region will show deeper coverage, thereby allowing the algorithm to identify it as a repeat. In de novo sequencing of a genome, repeats are usually assembled after the assembly of unique regions, and assemblers use multiple paired-ends to link a repetitive contig to a unique one. However, when the objective of research is not to obtain a complete genome but merely single sequence families, the problem of assembling contigs into more extended ones is less stringent. In addition, in the case of repetitive DNA families, lower coverage can reduce the occurrence of multi-reads and hence can improve repeat assembly into contigs and repeat identification and reconstruction [15–18].
The sunflower genome is large (around 3,500 Mbp, ). The repetitive component has been recently characterized using a Sanger-sequenced small insert library . This library provided a first set of sequences (1,638) that were used to analyze the composition of sunflower genome in terms of types and abundance of repetitive elements. The fraction of repetitive sequences amounted to 62% of the sequences, while the putative functional genes accounted for 4%. The largest component of the repetitive fraction was represented by long terminal repeat (LTR) retrotransposons, especially of the Gypsy superfamily. Class II transposable elements were barely represented in that library.
A larger effort to characterize the repetitive component of the sunflower genome was then made by analyzing approximately 25% of the genome from 454 random sequence reads . In this study, the sunflower genome was shown to be composed of over 81% transposable elements, 77% of which were long terminal repeat (LTR) retrotransposons.
The retrotransposon component of the sunflower genome was also analysed in detail by assembling and analyzing bacterial artificial chromosome (BAC) clones [21, 22]. Buti et al.  analysed 3 BAC clones, identifying 18 full-length and 6 incomplete LTR retrotransposons. Among LTR-retrotransposons, non-autonomous elements (the so-called LARDs [23, 24]), which do not carry any protein-encoding sequence, were discovered for the first time in sunflower. The insertion time of intact retroelements was measured, based on the divergence of sister LTRs. All isolated elements were inserted relatively recently, especially those belonging to the Gypsy superfamily.
These results were confirmed and extended by Staton et al. . The LTR retrotransposon fraction was shown to be composed in large part by chromodomain-containing Gypsy LTR retrotransposons. The authors showed that there is a bias in the efficiency of homologous recombination in removing LTR retrotransposon DNA, and provided insight into the mechanisms associated with the composition of the transposable element fraction in the sunflower genome. They also showed that most intact LTR retrotransposons have likely inserted since the origin of this species, providing further evidence that biased LTR retrotransposon activity has played a major role in shaping the DNA landscape of the sunflower genome.
In other studies, retrotransposons of the sunflower have been shown to be conserved within the Helianthus genus  and transcriptionally active [20, 26–28]. Fluorescent in situ hybridization studies have suggested that the Gypsy and Copia superfamilies are most frequent in the heterochromatic regions close to centromeres and telomeres, respectively [29–31]. The genomic organization of Gypsy elements is conserved also in hybrid sunflower species derived from the common sunflower, despite large amplification of these elements in the genome of such species [31–33].
The aim of the present study is to verify the suitability of using different approaches of de novo assembling sequence reads obtained by NGS procedures (Illumina and 454) to gain a comprehensive characterization of the repetitive component of a plant species (Helianthus annuus), whose large-sized genome is being sequenced thanks to the efforts of an international sequencing consortium [21, 34]. The repeat structure of the sunflower genome obtained in this study is validated by comparison with those obtained using a sunflower Sanger-sequenced small insert library , Sanger- or 454-sequenced sunflower BAC clones [21, 22], and sunflower de novo assembled 454 reads .
Besides developing further resources needed to sequence the sunflower genome, this study highlights the extent to which the repetitive portion of a plant genome can be characterized using NGS, and describes the utility and concerns raised by NGS methods of surveying such sequences.
Comparison of different assembled sequence sets
By varying sequencing technology (Illumina or 454), coverage (0.55 ×-1.25 ×), assemblers and assembly procedures (with or without splitting of read packages), different genomic databases were produced. On the whole, it can be observed that for each of three packages of reads (Illumina, 454 large, and 454 small read packages) the split subpackages resulted in the production of a lower number of contigs (Table 1). However, contigs were far more repetitive than those produced by simple assembly of whole reads, as shown by higher values of average coverage (Table 1, Figure 1). In fact, sequences assembled from the split sets were from about three-fold (for 454 large package) to more than 27-fold (for Illumina reads) more repeated in the genome than those assembled from unsplit sets.
The annotations of the six sets of assembled sequences show large differences in functional composition (Figure 2). Differences were especially pronounced when the same set of reads was split into subpackages prior to assembly. Figure 2 shows that low redundancy sequences such as putative genes or non-LTR retrotransposons were more common when the assemblies were conducted with no preliminary splitting. In contrast, preliminary splitting resulted in the assembly of larger percentages of LTR-retrotransposons. This is especially true for Illumina reads (Figure 2), probably because of their shorter length compared to 454 reads. Also, contigs for which no significant similarity was found in the existing databases were more (and in certain cases, much more) frequent in the unsplit assembled read set (Figure 2).
The six assembled sequence sets (with and without splitting) from the Illumina, 454 large, and 454 small sets of reads were each assembled two by two (split and unsplit) and annotated. The functional composition of the three resulting assemblies was similar, except that the frequency of retrotransposons was higher in both 454 packages than in the Illumina read set. In addition, a larger frequency of unclassified sequences was obtained using the Illumina read set (Figure 2).
Because of the large differences in average coverage and functional composition among the six assembled sequence sets, a further assembly was performed to produce a comprehensive genomic sequence set for sunflower. A total of 283,800 sequences (including 54,427 supercontigs and 229,538 individual contigs) were obtained, representing a whole genome set of assembled sequences (WGSAS).
The reliability of this method to obtain accurate sequences was tested by comparison of these sequences to available, Sanger-sequenced ones. Twenty alignments between assembled contigs and real sunflower DNA sequences are shown as Additional file 1. Mismatches related to transitions/transversions represent only 7.1% of 12,727 aligned nucleotides, indels amount to 0.4%.
SUNREP, a database of sunflower repetitive sequences
The WGSAS was mapped with the large set of Illumina reads as above. The distribution of average coverage of the WGSAS is reported in Figure 3. The average coverage was used as a parameter by which the repetitive sequences could be discriminated from the others. In plants much of the genome may be repeated because of the polyploidy events that have occurred during their evolutionary history (, as an example). Therefore, we evaluated sequence redundancy in relation to the average coverage of five sunflower gene sequences that were considered as unique reference sequences. By mapping Illumina reads to the WGSAS to which the five genes were added, we obtained for those sequences an average coverage of 6.6. We conservatively identified as repeated sequences all of those contigs with an average coverage higher than five-fold the mean average coverage of the five reference sequences (6.6 × 5 = 33.0). By this method, we identified 47,924 repeated sequences that constitute a database of repetitive sequences of sunflower, hereafter called SUNREP. The remaining 235,876 sequences of WGSAS were classified as unique or low redundant.
The distribution of different sequence types in SUNREP is reported in Table 2. It can be observed that 11.50% of sequences included in SUNREP did not find any hits in the public databases used for annotation. Among the annotated sequence types, retrotransposons were by far the most represented in SUNREP. Of LTR-retrotransposons, sequences belonging to the Gypsy superfamily were 2.3-fold more represented than those belonging to the Copia superfamilies.
Interestingly, a large fraction of sequences showed similarity to LTR-retrotransposons, but the superfamily could not be determined. Such elements lack coding sequence, are non-autonomous and usually species-specific. They can be discovered only when long sequences are available because their identities are based on structural features and not on sequence similarity to retrotransposon coding domains. In this study, we identified these elements only by their sequence similarity to those first reported by Buti et al. . Non-LTR retrotransposons were poorly represented, as frequently observed in plant genomes.
Putative DNA transposons accounted for 1,334 sequences. A portion of these were classified as DNA transposons according to sequence similarity to the short domain of the transposase gene. All types of plant DNA transposons were putatively found in SUNREP, with a prevalence of MITEs and Helitrons.
SUNREP contigs showing sequence similarity to LTR-REs, non-LTR REs, and DNA transposons were also analysed using an all-by-all BLAST search to estimate the occurrence in SUNREP of similar sequences within those repeat classes, i.e. sequences that were assembled separately, even though sharing some sequence similarity. Each class of repeats was subdivided into families (i.e., composed by at least 2 SUNREP sequences) and singletons (i.e., sequences that did not share similarity). The distribution of such families according to the number of sequences for each of them is reported in Figure 4.
The most redundant family, belonging to the Gypsy repeat superfamily, included only 96 of the 47,924 sequences of SUNREP (0.20%). Only four Gypsy families were composed of more than 50 SUNREP sequences. Considering the 30 most numerous LTR-REs families, the vast majority belonged to the Gypsy superfamily (Figure 5). Among the 30 most numerous DNA transposons, the most common families belonged to the Helitron class, followed by putative MITEs (Figure 5). It should be noted that the number of sequences that belong to a family in SUNREP does not reflect the redundancy of that family in the genome.
Another class of repeats was classified as unknown, but showed similarity to sequences previously isolated and whose redundancy has already been measured in sunflower by molecular and bioinformatic procedures . These unknown repeats were largely represented in SUNREP. Many of these showed similarity to the most repeated sequence identified by Cavallini et al. , the so-called Contig 61, whose nature remains to be ascertained.
Finally, according to BLAST analysis, 483 sequences of SUNREP showed similarity to putative protein encoding sequences. Of these, 123 were classified as hypothetical proteins, without any further annotation, the others belonged to 199 different gene families, of which 5 were represented by at least 9 sequences in the database (Table 3).
The most redundant genes encode the NBS-LRR class of proteins, receptors that recognize highly variable pathogen effectors . Another redundant gene family encode DNAJ proteins. These proteins function in association with Hsp70 molecular chaperones to facilitate protein folding and play an active role in regulating normal cellular events like protein degradation, morphogenesis and cell cycle progression .
The third redundant gene family is very heterogeneous, encoding proteins with un-specified protein-kinase domains, that are involved in the transduction of signals to binding factors, to the centromeres, and to other effectors. Beside aspecific kinases, also serine/threonine/tyrosine kinases are encoded by a redundant gene family. The last redundant and heterogeneous gene family encodes F-box motif containing proteins.
Sunflower genome composition
In the case of a small insert library  or of a whole genomic shotgun sequence library , the composition of the sequence set directly reflects the composition of the sunflower genome. Conversely, in the case of a sequence set obtained by assembling Illumina and 454 reads, the simple composition of the set cannot offer a picture of the genome composition, because repeated sequences are assembled together and hence are underestimated. Consequently, we evaluated the composition of the sunflower genome by counting the number and percentage of reads that mapped to each sequence in the WGSAS. Mapping results are summarized in Table 4.
Based on their similarity to the sequences in the organellar database, we estimated that more than 11.6 millions of reads were of organellar DNA origin. Regarding other reads, around 42 million reads did not match any assembled sequence, indicating that the WGSAS does not cover the entire genome, as expected having assembled only a total of 2.66 × coverage. It is likely that much of the missing sequences were low copy-number regions in the genome and that the relatively low coverage used in our study did not allow assembly of such loci. Such low-copy sequences could be protein encoding genes or rare forms of repeats whose sequence was degenerated until becoming unique. Some of these unmapped reads also likely represent sequencing errors of some kind. On the other hand, it is also possible that stringent assembly procedures and shorter reads affecting alignment stringency have contributed to increase the number of unaligned reads.
Considering that Illumina reads in our experiments were sampled without bias for particular sequence types, the percentage of reads that matched to a sequence class indicated the proportion of that sequence class in the sunflower genome. So, it was estimated that the percentage of repetitive sequences in the H. annuus genome was very high, amounting at least to 52.13% (see Table 4), while unique or low redundant sequences (that should include the vast majority of protein-encoding genes) represented only 12.29% of the genome at least. The rest (35.58%) of the genomic reads did not match any contigs.
Sunflower genome composition was estimated also in terms of sequence types. The frequency of each repeat type was calculated based on mapping the WGSAS with the 4 x coverage of Illumina reads and counting the number of reads matching each sequence type. Such frequencies are reported in Table 5, adopting the nomenclature proposed by Wicker et al. .
It can be observed that retrotransposons (especially LTR-retrotransposons) were by far the most abundant class of sequences in the sunflower genome, accounting for at least 79.53% of the reads matching the WGSAS, while DNA transposons and non-LTR retrotransposons showed very low percentages (Table 5). Of LTR-retrotransposons, the vast majority belonged to the Gypsy superfamily, which is 2.56-fold represented compared to Copia superfamily. A large amount of the genome (9.90%) was apparently made up of LTR-retrotransposons of unknown superfamily. It is presumable that their frequency in the sunflower genome was underestimated, and will increase after the sunflower genome sequence becomes available.
In other analyses we mapped Illumina reads to a sample of 19 intact LTR-retrotransposons of sunflower, isolated by Buti et al. , to estimate the equilibrium between retrotransposon replication and retrotransposon loss. Illumina reads were mapped to these retrotransposons, keeping separated LTR sequences from the respective inter-LTR region (that is the encoding region for Gypsy and Copia retroelements and an apparently non-encoding sequence for unknown retroelements, respectively). The results of mapping are reported in Table 6. It can be observed that the ratios between LTR and inter-LTR average coverage ranged from 0.01 to 6.58. If all retrotransposons belonging to one and the same family were intact, i.e. composed by two LTRs and one inter-LTR region, the ratio should have been 2. For 9 out of 19 analysed LTR-REs the ratio was higher than 2, indicating the occurrence of solo-LTRs of that RE family in the genome. The other LTR-REs had a ratio ranging from 0.01 to 1.25, i.e. the inter-LTR region was more represented in the genome than the LTR. This result can be explained by the presence of different families that share, at least in part, the inter-LTR region. Interestingly, analysing separately Gypsy, Copia, and unknown elements, the mean ratio between LTR and inter-LTRs average coverage was higher than 2 only for Gypsy elements (Table 6).
Discussion and conclusion
In our experiments, different strategies were used for assembling original sunflower sequence reads and for obtaining contigs; i.e. different packages of reads (Illumina and 454) were subdivided into low-coverage subpackages prior to assembly.
Similar levels of sequence coverage have proven to be efficient in generating a considerable amount of biologically useful information and genomic resources in other species [15, 38]. By using low genome coverage, most of the assembled contigs do not represent specific genomic loci; instead, they are probably composed of reads derived from multiple copies of repetitive elements, thus representing “consensus” sequences of genomic repeats . Even though the exact sequence of this consensus does not necessarily occur in the genome, this representation of repetitive elements is sufficiently accurate to enable amplification of whole length repetitive elements by PCR . Indeed, our comparison of assembled contigs with available Sanger sequences demonstrates good correspondence between virtual and real sequences.
Our results clearly show that splitting the original packages of reads into a number of subpackages allowed us to assemble more contigs similar to repetitive sequences, although assembled contigs were fewer (in the case of Illumina reads, much fewer) than those obtained by assembling the sets of reads prior to splitting. The difference in number and redundancy of assembled sequences was more striking for Illumina reads than for 454 reads, probably due to the short length of Illumina reads. However, splitting the packages of reads did not apparently affect the mean length and the N50 of the assembled contigs.
The sequence sets obtained by using different pre-assembly approaches were different in sequence types and redundancy. Hence, the production of a WGSAS by further assembling of the different sets provided a more complete picture of sunflower genome composition. In addition, the analysis of redundancy based on mapping Illumina reads onto the WGSAS allowed us to quantify the redundancy of each contig.
Sunflower genome composition has already been ascertained using other methods, i.e. biochemical analyses , sequencing and analysis of a small insert library , whole genome 454 sequencing . All these analyses may have some potential weakness. Biochemical analyses  obviously do not consider DNA sequence but only denaturation and reassociation kinetics of DNA, so rare forms of repeats are excluded such as, for example, retrotransposon remnants. The Sanger-sequenced small insert library  comprised only 1,638 sequences, so conclusions are subject to sampling errors. The whole genome 454 sequence database  is based only on sequence similarity, however the number of 454 reads used (total coverage 0.23 ×) might be not sufficient to ensure accurate estimation of genome composition.
The present analyses showed that sunflower genome is mostly composed of LTR-retrotransposons (78.8%), similar to that already reported [21, 34]. It is known that the genome size is determined during evolution by an equilibrium between enlargement of the genome by insertion of REs and RE-mediated DNA removal [41–43]. DNA rearrangements, illegitimate recombination, and unequal homologous recombination drive DNA removal in plants by a number of mechanisms, such as the repair of double strand breaks (nonhomologous end-joining) and slipstrand mispairing [44–49].
The observed large number of retrotransposons indicates that such elements have been actively replicating during the evolution of this species. Recent studies have reported that sunflower LTR-REs are transcribed even at present [20, 26–28] and, in at least one case, RE transcription was shown to be followed by RE insertion .
Mapping Illumina reads to a set of 19 available intact LTR-retrotransposons suggested the occurrence of numerous solo-LTRs for 9 out of 19 REs, although the occurrence of REs sharing LTRs but having different internal regions cannot be ruled out and could lead to an overestimation of solo-LTR frequencies. Solo-LTRs are typically produced by illegitimate recombination. Our data suggest that massive amplification of these elements in the sunflower genome was partly counterbalanced by substantial DNA loss, especially related to Gypsy elements, although in other studies solo-LTRs have been found commonly for Copia elements as well [20, 21]. It is obvious that a very large number of intact retroelements are needed to validate this analysis.
Concerning the different RE superfamilies, the ratio between Gypsy and Copia retrotransposon frequencies amounted to 2.29, confirming the greater abundance of the former superfamily. This ratio is generally species-specific. Gypsy to Copia frequency ratio is even higher in papaya (5:1, ), Sorghum (4:1, ), and rice (3:1, ) than in the sunflower genome. In other cases, as in maize , poplar , and olive (Barghini, personal communication) a similar abundance of the two superfamilies was observed. Finally, in grapevine an opposite trend was found, with Copia elements two-fold more represented than Gypsy ones .
The large abundance of Gypsy elements compared to Copia can be explained by two hypotheses: Gypsy elements have been more active during sunflower evolution and/or they have been active more recently, so that are more easily recognizable by similarity searches, having been subjected to fewer mutations. Dating retrotransposon insertions in the sunflower genome indicate that Gypsy elements are generally younger than Copia, though some Copia elements are relatively young as well [21, 22].
Retrotransposon and DNA transposon sequences included in the redundant fraction of the WGSAS (SUNREP) were also assigned to different families within each superfamily, by an all-by-all BLAST search. The number of sequences composing each family was generally low, confirming that there are not prominent transposon families in this species [20, 21].
In a previous study , a different approach was used for determining the composition of different repeat types in terms of families, by using the graph-based method of Novak et al. . The families of LTR-retrotransposons and DNA transposons generally match the results reported in Staton et al. , with the exception of putative MITEs, that are more frequent than previously observed in other studies. Interestingly, the most frequent DNA transposon family belongs to the Helitron superfamily and is comprised of a number of sequences comparable to that of the most numerous LTR-RE subfamilies. Also the graph-based study included one Helitron subfamily among the 20 most redundant ones in the sunflower genome; all the others belonging to the LTR-RE class.
The results obtained by Staton et al.  and those reported in this study indicate that both the method by Novak et al.  and the all-by-all BLAST search (performed in our experiments) allow a precise estimation of repeat superfamilies and families. The first method allows information to be gained on repeat structure and provides putative consensus sequences of the repeat; all-by-all BLAST search (preceded by assembling all available sequences) can be applied to larger sets of reads.
Finally, mapping data indicated that a number of contigs showing similarity to putative protein encoding genes are to be considered as redundant. In many cases such contigs showed similarity to gene families already known to be repeated in plant genomes, such as NBS-LRR genes . In other cases, it is likely that functional domains, and not genes, are the cause of apparent redundancy. For example, F-box proteins are identified by the presence of protein interaction domains that bind ubiquitinilation targets and include a large variety of proteins . For other contigs, it might be that a gene (or a gene fragment) lies close to a repeated sequence and the redundancy of that contig is related to the repeated sequence and not to the gene sequence. Interestingly, our data indicated the occurrence, in the sunflower genome, of a relatively high number of putative Helitrons. These sequences, of transposable origin, are known to include DNA retro-transcribed on RNA transcripts  and might be responsible for the relatively high frequency of gene fragments in the redundant fraction of the WGSAS.
In conclusion, the results of our experiments show how different data partitioning and assembly approaches can be used to obtain valuable insights on genome composition using NGS technologies, either 454 or Illumina, or both technologies combined. Concerning sunflower, our data confirm the repeat structure of the genome and give new insights on different aspects of it. Moreover, they will facilitate the functional annotation of the H. annuus genome that is currently being sequenced and will be used for studies of intra- and interspecific variability related to H. annuus and its relatives.
Leaf tissue was sampled from a single individual from a highly inbred sunflower cultivar HA412-HO (PI 642777). Total genomic DNA was extracted using a CTAB procedure  and randomly sheared into fragments for sequencing.
For Illumina sequencing, library preparation followed the standard multi-step Illumina protocol . Ligated, size-selected fragments (~280 - 320 bp) were amplified through 18 cycles of polymerase chain reaction (PCR), using Phusion High-Fidelity PCR Master Mix (New England Biolabs, Ipswich, MA, USA) and standard Illumina primers. Resulting product concentrations were determined using a Nanodrop 1000 (ThermoFisher Scientific, Wilmington, DE, USA).
Enriched, cleaned product was diluted to 5 pM and submitted for sequencing on an Illumina GAII sequencer at the University of British Columbia (Vancouver). Denaturation, cluster generation, and subsequent sequencing followed the manufacturer’s recommendations. Image analysis, base-calling and error estimation were performed using the Illumina GA Pipeline version 1.5. Low quality bases, empty reads, and adapter sequences were removed using CLC-BIO Genomic Workbench, version 5.1 (CLC-BIO). A total of 28,236,626 random genomic sequences (mean read length 101 nt, total coverage 0.86 ×) were obtained. In a second experiment, 130,345,608 Illumina reads (mean read length 101 nt, total coverage 4 ×) were obtained.
In other experiments, standard 454 sequencing runs were performed on a Roche 454 GS FLX sequencer (Roche, http://www.roche.com) with XLR (Titanium) chemistry at The McGill University and Génome Québec Innovation Centre (Montreal, Canada). We obtained 22,666,169 WGS sequences, for a total of 1.8 × coverage, subdivided in two packages, a large one (mean read length 166.2 nt, coverage 1.25 ×), and a small one (mean read length 349.4 nt, coverage 0.55 ×).
In a first approach, an Illumina read set (genome coverage 0.86 ×), 454 large read set (1.25 ×), and 454 small read set (0.55 ×) were assembled separately. Each read set was assembled using CLC-BIO based on unambiguous overlaps. The resulting contigs were further assembled separately using Minimus2 software . CLC-BIO assembly parameters were: minimum contig length = 300; minimum distance = 200; max distance = 600 for Illumina package; minimum contig length = 300 for 454 large package; and minimum contig length = 300; minimum distance = 5,000; max distance = 15,000 for 454 small package. Minimus2 assembly parameters were REFCOUNT = 0 and MINID = 90. This second assembly produced three sets of supercontigs and single contigs. The resulting contigs were assessed as to their number, length, and N50 (Table 1).
In a second approach, each of the three read sets was split into low coverage subpackages. For the 454 large and small sets of reads, the split was performed to obtain 26 and 18 subpackages, with 0.048 and 0.031 × coverage, respectively. For the Illumina reads, we prepared 565 subpackages, each with less coverage (0.0015 ×) than those used for 454 reads, because preliminary experiments showed us that this level of coverage allows the largest recovery of repeated sequences (Barghini, personal communication). Each subpackage was individually assembled using CLC-BIO, then each group of subpackages was further assembled using Minimus2 with the above described parameters. This second assembly produced three additional sets of supercontigs and single contigs. Possible contaminants resembling organellar sequences were then removed from the six sets of contigs by masking them against an in-house sunflower organellar sequence database using RepeatMasker (http://www.repeatmasker.org/), and excluding all contigs showing at least 1% of their length similar to organellar sequences. The remaining contigs of the six sets were further assembled, using Minimus 2 with REFCOUNT = 0 and MINID = 90: in a first step, the six sets of contigs were assembled two by two (split and unsplit), to produce three sets of sequences; in a second step, the three sets were assembled into a unique set of assembled sequences. For all assemblies also, the resulting contigs were assessed for number, length, and N50 (Table 1).
Redundancy estimation of sequences
Relative redundancy of each sequence in the six sets of assembled sequences and in the WGSAS was estimated by mapping the sequences with a large Illumina sequence read set (total coverage 4 ×). Mapping was performed using CLC-BIO, which randomly places multi-reads, hence the number of mapped reads to a single sequence is only an indication of its redundancy. On the other hand, if all sequences of a repeat family or class are taken together, the total number of mapped reads (in respect to total genomic reads) indicates the effective redundancy of that family or class. To establish mapping parameters, sixty sequences were selected for which redundancy had been previously determined by slot blot and hybridization (; Giordani, personal communication). For these sequences, correlations were calculated between their known redundancy and their average coverage (the sum of the bases of the aligned part of all the reads divided by the length of the reference sequence) by using different parameters (mismatch cost, deletion cost, insertion cost, length fraction, similarity, Additional file 2). The parameters determining the largest correlation were selected to be used in the subsequent mapping of different sequence sets. The means and distributions of average coverage values for each contig of the six sets are reported in Table 1 and Figure 1, respectively.
In the case of the WGSAS, to evaluate the redundancy of DNA sequences, the same 4 × Illumina sequence read set was mapped onto the WGSAS plus one actin-encoding gene (FJ487620.1) and four unique gene sequences , encoding a lipid transfer protein (FR671365.1), a z-carotene desaturase (FR671183.1), an auxin-binding protein (FR671175.1), and an ABA-responsive C5 protein (FR671167.1). Then the average coverage was calculated for each gene sequence. We established an average coverage of 33.0 (i.e., five-fold that of the mean average coverage of the five above sequences, 6.6) as an arbitrary threshold for discriminating repeated sequences from unique or low-redundant ones.
The mapping procedure, using CLC-BIO and the same 4x Illumina sequence read set as above, was adopted for determining the relative redundancy of LTR and inter-LTR regions of a sample of sunflower retrotransposons.
Annotation of repetitive sequences
Sequences belonging to all assembled sequence sets were searched for homologies by using the NCBI BLAST with an E-value cut-off of 10-10 in the NR NCBI database (http://blast.ncbi.nlm.nih.gov/), containing all the non-redundant protein sequences, in the database RepBase  and in other available sunflower sequence sets [20, 22]. One region might match two or more homologous elements in the database. We then removed the redundant annotations by keeping only the region with the longest match and significant E-value. We also performed an analysis of SUNREP by using RepeatMasker to isolate microsatellites and low complexity sequences.
All sequences belonging to the WGSAS that were not included in SUNREP were annotated using RepeatMasker against SUNREP sequences, to isolate and annotate low copy remnants of repeated sequences, as transposons that had accumulated mutations.
After annotation, all SUNREP sequences similar to LTR-REs, non-LTR-REs, and DNA transposons were grouped into families according to their sequence similarity by performing an all-by-all search using BLAST with the following parameters: -r 4 -q −5 -e 1e-50. For each sequence, the most similar sequence (with a similarity score E-50 at least) in the database was recorded. Then, each sequence sharing similarity was attributed to the same family.
Greilhuber J, Borsch T, Müller K, Worberg A, Porembski S, Barthlott W: Smallest angiosperm genomes found in Lentibulariaceae, with chromosomes of bacterial size. Plant Biol. 2006, 8: 770-777. 10.1055/s-2006-924101.
Pellicer J, Fay MF, Leitch IJ: The largest eukaryotic genome of them all?. Bot J Linnean Soc. 2010, 164: 10-15. 10.1111/j.1095-8339.2010.01072.x.
Morgante M, De Paoli M, Radovic S: Transposable elements and the plant pan-genomes. Curr Opin Plant Biol. 2007, 10: 149-155. 10.1016/j.pbi.2007.02.001.
Britten RJ: Transposable element insertions have strongly affected human evolution. Proc Natl Acad Sci USA. 2010, 107: 19945-19948. 10.1073/pnas.1014330107.
Morgante M, Brunner S, Pea G, Fengler K, Zuccolo A, Rafalski A: Gene duplication and exon shuffling by helitron-like transposons generate intraspecies diversity in maize. Nat Genet. 2005, 37: 997-1002. 10.1038/ng1615.
Slotkin RK, Martienssen R: Transposable elements and the epigenetic regulation of the genome. Nat Rev Genet. 2007, 8: 272-285.
Glenn TC: Field guide to next-generation DNA sequencers. Mol Ecol Resour. 2011, 11: 759-769. 10.1111/j.1755-0998.2011.03024.x.
Hudson ME: Sequencing breakthroughs for genomic ecology and evolutionary biology. Mol Ecol Resour. 2008, 8: 3-17. 10.1111/j.1471-8286.2007.02019.x.
Rokas A, Abbot P: Harnessing genomics for evolutionary insights. Trends Ecol Evol. 2009, 24: 192-200. 10.1016/j.tree.2008.11.004.
Ekblom R, Galindo J: Applications of next generation sequencing in molecular ecology of non-model organisms. Heredity. 2011, 107: 1-15. 10.1038/hdy.2010.152.
Alkan C, Coe BP, Eichler EE: Genome structural variation discovery and genotyping. Nature Rev Genet. 2011, 12: 363-376. 10.1038/nrg2958.
Treangen TJ, Salzberg SL: Repetitive DNA and next-generation sequencing: computational challenges and solutions. Nature Rev Genet. 2012, 13: 36-46.
Simola DF, Kim J: Sniper: improved SNP discovery by multiply mapping deep sequenced reads. Genome Biol. 2011, 12: R55-10.1186/gb-2011-12-6-r55.
Pevzner PA, Tang H, Waterman MS: An Eulerian path approach to DNA fragment assembly. Proc Natl Acad Sci USA. 2001, 98: 9748-9753. 10.1073/pnas.171285098.
Rasmussen DA, Noor MAF: What can you do with 0.1× genome coverage? A case study based on a genome survey of the scuttle fly Megaselia scalaris (Phoridae). BMC Genomics. 2009, 10: 382-10.1186/1471-2164-10-382.
Lee RM, Thimmapuram J, Thinglum KA, Gong G, Hernandez AG, Wright CL, Kim RW, Mikel MA, Tranel PJ: Sampling the Waterhemp (Amaranthus tuberculatus) genome using pyrosequencing technology. Weed Sci. 2009, 57: 463-469. 10.1614/WS-09-021.1.
Swaminathan K, Alabady MS, Varala K, De Paoli E, Ho I, Rokhsar DS, Arumuganathan AK, Ming R, Green PJ, Meyers BC, Moose SP, Hudson ME: Genomic and small RNA sequencing of Miscanthus × giganteus shows the utility of sorghum as a reference genome sequence for Andropogoneae grasses. Genome Biol. 2010, 11: R12-10.1186/gb-2010-11-2-r12.
Straub SCK, Fishbein M, Livshultz T, Foster Z, Parks M, Weitemier K, Cronn RC, Liston A: Building a model: developing genomic resources for common milkweed (Asclepias syriaca) with low coverage genome sequencing. BMC Genomics. 2011, 12: 211-10.1186/1471-2164-12-211.
Baack EJ, Whitney KD, Rieseberg LH: Hybridization and genome size evolution: timing and magnitude of nuclear DNA content increases in Helianthus homoploid hybrid species. New Phytol. 2005, 167: 623-630. 10.1111/j.1469-8137.2005.01433.x.
Cavallini A, Natali L, Zuccolo A, Giordani T, Jurman I, Ferrillo V, Vitacolonna N, Sarri V, Cattonaro F, Ceccarelli M, Cionini PG, Morgante M: Analysis of transposons and repeat composition of the sunflower (Helianthus annuus L.) genome. Theor Appl Genet. 2010, 120: 491-508. 10.1007/s00122-009-1170-7.
Staton SE, Bakken BH, Blackman BK, Chapman MA, Kane NC, Tang S, Ungerer MC, Knapp SJ, Rieseberg LH, Burke JM: The sunflower (Helianthus annuus L.) genome reflects a recent history of biased accumulation of transposable elements. Plant J. 2012, 72: 142-153. 10.1111/j.1365-313X.2012.05072.x.
Buti M, Giordani T, Cattonaro F, Cossu RM, Pistelli L, Vukich M, Morgante M, Cavallini A, Natali L: Temporal dynamics in the evolution of the sunflower genome as revealed by sequencing and annotation of three large genomic regions. Theor Appl Genet. 2011, 123: 779-791. 10.1007/s00122-011-1626-4.
Kalendar R, Vicient CM, Peleg O, Anamthawat-Jonsson K, Bolshoy A, Schulman AH: Large retrotransposon derivatives: abundant, conserved but nonautonomous retroelements of barley and related genomes. Genetics. 2004, 166: 1437-1450. 10.1534/genetics.166.3.1437.
Wicker T, Sabot F, Hua-Van A, Bennetzen JL, Capy P, Chalhoub B, Flavell A, Leroy P, Morgante M, Panaud O, Paux E, SanMiguel P, Schulman AH: A unified classification system for eukaryotic transposable elements. Nature Rev Genet. 2007, 8: 973-982. 10.1038/nrg2165.
Vukich M, Schulman AH, Giordani T, Natali L, Kalendar R, Cavallini A: Genetic variability in sunflower (Helianthus annuus L.) and in the Helianthus genus as assessed by retrotransposon based molecular markers. Theor Appl Genet. 2009, 119: 1027-1038. 10.1007/s00122-009-1106-2.
Buti M, Giordani T, Vukich M, Gentzbittel L, Pistelli L, Cattonaro F, Morgante M, Cavallini A, Natali L: HACRE1, a recently inserted Copia-like retrotransposon of sunflower (Helianthus annuus L.). Genome. 2009, 11: 904-911.
Vukich M, Schulman AH, Giordani T, Natali L, Kalendar R, Cavallini A: Copia and Gypsy retrotransposons activity in sunflower (Helianthus annuus L.). BMC Plant Biol. 2009, 9: 150-10.1186/1471-2229-9-150.
Kawakami T, Dhakal P, Katterhenry AN, Heatherington CA, Ungerer MC: Transposable element proliferation and genome expansion are rare in contemporary sunflower hybrid populations despite widespread transcriptional activity of LTR retrotransposons. Genome Biol Evol. 2011, 3: 156-167. 10.1093/gbe/evr005.
Santini S, Cavallini A, Natali L, Minelli S, Maggini F, Cionini PG: Ty1/Copia- and Ty3/Gypsy-like DNA sequences in Helianthus species. Chromosoma. 2002, 111: 192-200. 10.1007/s00412-002-0196-2.
Natali L, Santini S, Giordani T, Minelli S, Maestrini P, Cionini PG, Cavallini A: Distribution of Ty3-Gypsy- and Ty1-Copia-like DNA sequences in the genus Helianthus and other Asteraceae. Genome. 2006, 49: 64-72. 10.1139/g05-058.
Staton SE, Ungerer MC, Moore RC: The genomic organization of Ty3/Gypsy-like retrotransposons in Helianthus (Asteraceae) homoploid hybrid species. Am J Bot. 2009, 96: 1646-1655. 10.3732/ajb.0800337.
Ungerer MC, Strakosh SC, Zhen Y: Genome expansion in three hybrid sunflower species is associated with retrotransposon proliferation. Curr Biol. 2006, 16: R872-R873. 10.1016/j.cub.2006.09.020.
Ungerer MC, Strakosh SC, Stimpson KM: Proliferation of Ty3/Gypsy-like retrotransposons in hybrid sunflower taxa inferred from phylogenetic data. BMC Biol. 2009, 7: 40-10.1186/1741-7007-7-40.
Kane NC, Gill N, King ME, Bowers JE, Berges H, Gouzy J, Bachlava E, Langlade NB, Lai Z, Stewart M, Burke JM, Vincourt P, Knapp SJ, Rieseberg LH: Progress towards a reference genome for sunflower. Botany. 2011, 89: 429-437. 10.1139/b11-032.
Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449: 463-467. 10.1038/nature06148.
McHale L, Tan X, Koehl P, Michelmore RW: Plant NBS-LRR proteins: adaptable guards. Genome Biol. 2006, 7: 212-
Frugis G, Mele G, Giannino D, Mariotti D: MsJ1, an alfalfa DnaJ-like gene, is tissue-specific and transcriptionally regulated during cell cycle. Plant Mol Biol. 1999, 40: 397-408. 10.1023/A:1006215231492.
Swaminathan K, Varala K, Hudson ME: Global repeat discovery and estimation of genomic copy number in a large, complex genome using a high-throughput 454 sequence survey. BMC Genomics. 2007, 8: 132-10.1186/1471-2164-8-132.
Novak P, Neumann P, Macas J: Graph-based clustering and characterization of repetitive sequences in next-generation sequencing data. BMC Bioinformatics. 2010, 11: 378-10.1186/1471-2105-11-378.
Cavallini A, Zolfino C, Cionini G, Cremonini R, Natali L, Sassoli O, Cionini PG: Nuclear DNA changes within Helianthus annuus L.: cytophotometric, karyological and biochemical analyses. Theor Appl Genet. 1986, 73: 20-26.
Devos KM, Brown JKM, Bennetzen JL: Genome size reduction through illegitimate recombination counteracts genome expansion in Arabidopsis. Genome Res. 2002, 12: 1075-1079. 10.1101/gr.132102.
Ma J, Devos KM, Bennetzen JL: Analyses of LTR retrotransposon structures reveal recent and rapid genomic DNA loss in rice. Genome Res. 2004, 14: 860-869. 10.1101/gr.1466204.
Grover C, Hawkins J, Wendel J: Phylogenetic insights into the pace and pattern of plant genome size evolution. Plant genomes. Genome dynamics vol 4 (Volff JN). 2008, Basel (CH): Karger, 57-68.
Kalendar R, Tanskanen J, Immonen S, Nevo E, Schulman AH: Genome evolution of wild barley (Hordeum spontaneum) by BARE-1 retrotransposon dynamics in response to sharp microclimatic divergence. Proc Natl Acad Sci USA. 2000, 97: 6603-6607. 10.1073/pnas.110587497.
Ma J, Bennetzen JL: Rapid recent growth and divergence of rice nuclear genomes. Proc Natl Acad Sci USA. 2004, 101: 12404-12410. 10.1073/pnas.0403715101.
Neumann P, Koblizkova A, Navratilova A, Macas J: Significant expansion of Vicia pannonica genome size mediated by amplification of a single type of giant retroelement. Genetics. 2006, 173: 1047-1056. 10.1534/genetics.106.056259.
Ammiraju JS, Zuccolo A, Yu Y, Song X, Piegu P, Chevalier F, Walling JG, Ma J, Talag J, Brar DS, SanMiguel PJ, Jiang N, Jackson SA, Panaud O, Wing RA: Evolutionary dynamics of an ancient retrotransposon family provides insights into evolution of genome size in the genus Oryza. Plant J. 2007, 52: 342-351. 10.1111/j.1365-313X.2007.03242.x.
Hawkins JS, Hu G, Rapp RA, Grafenberg JL, Wendel JF: Phylogenetic determination of the pace of transposable element proliferation in plants: Copia and LINE-like elements in Gossypium. Genome. 2008, 51: 11-18. 10.1139/G07-099.
Morse AM, Peterson DG, Islam-Faridi MN, Smith KE, Magbanua Z, Garcia SA, Kubisiak TL, Amerson HV, Carlson JE, Nelson CD, Davis JM: Evolution of genome size and complexity in Pinus. PLoS One. 2009, 4: e4332-10.1371/journal.pone.0004332.
Ming R, Hou S, Feng Y, Yu Q, Dionne-Laporte A, Albert H, Suzuki JY, Tripathi S, Moore PH, Gonsalves D: The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus). Nature. 2008, 452: 991-997. 10.1038/nature06856.
Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A, Schmutz J, Spannagl M, Tang H, Wang X, Wicker T: The Sorghum bicolor genome and the diversification of grasses. Nature. 2009, 457: 551-556. 10.1038/nature07723.
The International Rice Genome Sequencing Project: The map based sequence of the rice genome. Nature. 2005, 436: 793-800. 10.1038/nature03895.
Meyers BC, Tingey SV, Morgante M: Abundance, distribution, and transcriptional activity of repetitive elements in the maize genome. Genome Res. 2001, 11: 1660-1676. 10.1101/gr.188201.
Cossu RM, Buti M, Giordani T, Natali L, Cavallini A: A computational study of the dynamics of LTR retrotransposons in the Populus trichocarpa genome. Tree Genet Genomes. 2012, 8: 61-75. 10.1007/s11295-011-0421-3.
Jin J, Cardozo T, Lovering RC, Elledge SJ, Pagano M, Harper JW: Systematic analysis and nomenclature of mammalian F-box proteins. Genes Devel. 2004, 18: 2573-2580. 10.1101/gad.1255304.
Doyle JJ, Doyle JL: Isolation of plant DNA from fresh tissue. Focus. 1989, 12: 13-15.
Solexa Inc: Protocol for Whole Genome Sequencing using Solexa Technology. BioTechniques Protocol Guide. 2007, New York: BioTechniques, 2006-
Sommer DD, Delcher AL, Salzberg SL, Pop M: Minimus: a fast, lightweight genome assembler. BMC Bioinformatics. 2007, 8: 64-10.1186/1471-2105-8-64.
Giordani T, Buti M, Natali L, Pugliesi C, Cattonaro F, Morgante M, Cavallini A: An analysis of sequence variability in eight genes putatively involved in drought response in sunflower (Helianthus annuus L.). Theor Appl Genet. 2011, 122: 1039-1049. 10.1007/s00122-010-1509-0.
Jurka J: Repbase update: a database and an electronic journal of repetitive elements. Trends Genet. 2000, 16: 418-420. 10.1016/S0168-9525(00)02093-X.
The generation of Illumina and 454 sequence data was funded by Genome Canada- and Genome BC. This assembly and analyses were supported by PRIN-MIUR, Italy, Project “SUNREP: caratterizzazione molecolare della componente ripetitiva del genoma di girasole”. We thank Megan Stewart for isolating the DNA used for sequencing and Thuy Nguyen for processing the raw reads and making them available for analyses.
All whole-genome-shotgun sequences described are available upon request or on NCBI Sequence Read Archive under SRA Project number 64989 (http://www.ncbi.nlm.nih.gov/bioproject/64989). The SUNREP database (sequences and annotations) is available at the repository sequence page of the Department of Agriculture, Food, and Environment of the University of Pisa (http://www.agr.unipi.it/Sequence-Repository.358.0.html).
The authors declare that they have no competing interests.
Conceived and designed the study: LN AC; Generated the sequence data: NG NCK LR; Performed assembling, redundancy analyses and annotation of sequence data: RMC; Participated in performing assembling, redundancy analyses and annotation of sequence data: EB TG MB FM; Participated in the interpretation and discussion of results and contributed to the writing of the paper: LN RMC EB TG MB FM MM NG NCK LR AC; Principal investigator for the “SUNREP” PRIN-MIUR project and coordinator of the study: LN; Wrote the paper: LN, TG, NCK, LR, and AC; All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1:Alignment of a sample of 20 assembled sequences of the sunflower whole genome database (inbred line HA412-HO, numerical codes, below) to Sanger sequences from the small insert library. (PDF 140 KB)
Additional file 2:The 5 sets of CLC-BIO parameters used for mapping of the Illumina reads to 60 sunflower DNA sequences with known redundancy (; Giordani, personal communication) and correlation between known copy number and average coverage.(PDF 52 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Natali, L., Cossu, R.M., Barghini, E. et al. The repetitive component of the sunflower genome as shown by different procedures for assembling next generation sequencing reads. BMC Genomics 14, 686 (2013). https://doi.org/10.1186/1471-2164-14-686
- Genome structure
- Next Generation Sequencing
- Repetitive DNA