Satellite DNA in Plants: More than Just Rubbish

TAREAN: a computational tool for identification and characterization of satellite DNA from unassembled short reads

Satellite DNA is one of the major classes of repetitive DNA, characterized by tandemly arranged repeat copies that form contiguous arrays up to megabases in length. This type of genomic organization makes satellite DNA difficult to assemble, which hampers characterization of satellite sequences by computational analysis of genomic contigs. Here, we present tandem repeat analyzer (TAREAN), a novel computational pipeline that circumvents this problem by detecting satellite repeats directly from unassembled short reads. The pipeline first employs graph-based sequence clustering to identify groups of reads that represent repetitive elements. Putative satellite repeats are subsequently detected by the presence of circular structures in their cluster graphs. Consensus sequences of repeat monomers are then reconstructed from the most frequent k-mers obtained by decomposing read sequences from corresponding clusters. The pipeline performance was successfully validated by analyzing low-pass genome sequencing data from five plant species where satellite DNA was previously experimentally characterized. Moreover, novel satellite repeats were predicted for the genome of Vicia faba and three of these repeats were verified by detecting their sequences on metaphase chromosomes using fluorescence in situ hybridization.

TAREAN: a computational tool for identification and characterization of satellite DNA from unassembled short reads

Satellite DNA is one of the major classes of repetitive DNA, characterized by tandemly arranged repeat copies that form contiguous arrays up to megabases in length. This type of genomic organization makes satellite DNA difficult to assemble, which hampers characterization of satellite sequences by computational analysis of genomic contigs. Here, we present tandem repeat analyzer (TAREAN), a novel computational pipeline that circumvents this problem by detecting satellite repeats directly from unassembled short reads. The pipeline first employs graph-based sequence clustering to identify groups of reads that represent repetitive elements. Putative satellite repeats are subsequently detected by the presence of circular structures in their cluster graphs. Consensus sequences of repeat monomers are then reconstructed from the most frequent k-mers obtained by decomposing read sequences from corresponding clusters. The pipeline performance was successfully validated by analyzing low-pass genome sequencing data from five plant species where satellite DNA was previously experimentally characterized. Moreover, novel satellite repeats were predicted for the genome of Vicia faba and three of these repeats were verified by detecting their sequences on metaphase chromosomes using fluorescence in situ hybridization.

The Peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity

BACKGROUND

The availability of the peach genome sequence has fostered relevant research in peach and related Prunus species enabling the identification of genes underlying important horticultural traits as well as the development of advanced tools for genetic and genomic analyses. The first release of the peach genome (Peach v1.0) represented a high-quality WGS (Whole Genome Shotgun) chromosome-scale assembly with high contiguity (contig L50 214.2 kb), large portions of mapped sequences (96%) and high base accuracy (99.96%). The aim of this work was to improve the quality of the first assembly by increasing the portion of mapped and oriented sequences, correcting misassemblies and improving the contiguity and base accuracy using high-throughput linkage mapping and deep resequencing approaches.

RESULTS

Four linkage maps with 3,576 molecular markers were used to improve the portion of mapped and oriented sequences (from 96.0% and 85.6% of Peach v1.0 to 99.2% and 98.2% of v2.0, respectively) and enabled a more detailed identification of discernible misassemblies (10.4 Mb in total). The deep resequencing approach fixed 859 homozygous SNPs (Single Nucleotide Polymorphisms) and 1347 homozygous indels. Moreover, the assembled NGS contigs enabled the closing of 212 gaps with an improvement in the contig L50 of 19.2%.

CONCLUSIONS

The improved high quality peach genome assembly (Peach v2.0) represents a valuable tool for the analysis of the genetic diversity, domestication, and as a vehicle for genetic improvement of peach and related Prunus species. Moreover, the important phylogenetic position of peach and the absence of recent whole genome duplication (WGD) events make peach a pivotal species for comparative genomics studies aiming at elucidating plant speciation and diversification processes.

The Peach v2.0 release: high-resolution linkage mapping and deep resequencing improve chromosome-scale assembly and contiguity

BACKGROUND

The availability of the peach genome sequence has fostered relevant research in peach and related Prunus species enabling the identification of genes underlying important horticultural traits as well as the development of advanced tools for genetic and genomic analyses. The first release of the peach genome (Peach v1.0) represented a high-quality WGS (Whole Genome Shotgun) chromosome-scale assembly with high contiguity (contig L50 214.2 kb), large portions of mapped sequences (96%) and high base accuracy (99.96%). The aim of this work was to improve the quality of the first assembly by increasing the portion of mapped and oriented sequences, correcting misassemblies and improving the contiguity and base accuracy using high-throughput linkage mapping and deep resequencing approaches.

RESULTS

Four linkage maps with 3,576 molecular markers were used to improve the portion of mapped and oriented sequences (from 96.0% and 85.6% of Peach v1.0 to 99.2% and 98.2% of v2.0, respectively) and enabled a more detailed identification of discernible misassemblies (10.4 Mb in total). The deep resequencing approach fixed 859 homozygous SNPs (Single Nucleotide Polymorphisms) and 1347 homozygous indels. Moreover, the assembled NGS contigs enabled the closing of 212 gaps with an improvement in the contig L50 of 19.2%.

CONCLUSIONS

The improved high quality peach genome assembly (Peach v2.0) represents a valuable tool for the analysis of the genetic diversity, domestication, and as a vehicle for genetic improvement of peach and related Prunus species. Moreover, the important phylogenetic position of peach and the absence of recent whole genome duplication (WGD) events make peach a pivotal species for comparative genomics studies aiming at elucidating plant speciation and diversification processes.

Evolutionary dynamics of satellite DNA repeats from Phaseolus beans

Common bean (Phaseolus vulgaris) subtelomeres are highly enriched for khipu, the main satellite DNA identified so far in this genome. Here, we comparatively investigate khipu genomic organization in Phaseolus species from different clades. Additionally, we identified and characterized another satellite repeat, named jumper, associated to khipu. A mixture of P. vulgaris khipu clones hybridized in situ confirmed the presence of khipu-like sequences on subterminal chromosome regions in all Phaseolus species, with differences in the number and intensity of signals between species and when species-specific clones were used. Khipu is present as multimers of ∼500 bp and sequence analyses of cloned fragments revealed close relationship among khipu repeats. The new repeat, named jumper, is a 170-bp satellite sequence present in all Phaseolus species and inserted into the nontranscribed spacer (NTS) of the 5S rDNA in the P. vulgaris genome. Nevertheless, jumper was found as a high-copy repeat at subtelomeres and/or pericentromeres in the Phaseolus microcarpus lineage only. Our data argue for khipu as an important subtelomeric satellite DNA in the genus and for a complex satellite repeat composition of P. microcarpus subtelomeres, which also contain jumper. Furthermore, the differential amplification of these repeats in subtelomeres or pericentromeres reinforces the presence of a dynamic satellite DNA library in Phaseolus.

Tandem repeats of Allium fistulosum associated with major chromosomal landmarks

Tandem repeats are often associated with important chromosomal landmarks, such as centromeres, telomeres, subtelomeric, and other heterochromatic regions, and can be good candidates for molecular cytogenetic markers. Tandem repeats present in many plant species demonstrate dramatic differences in unit length, proportion in the genome, and chromosomal organization. Members of genus Allium with their large genomes represent a challenging task for current genetics. Using the next generation sequencing data, molecular, and cytogenetic methods, we discovered two tandemly organized repeats in the Allium fistulosum genome (2n = 2C = 16), HAT58 and CAT36. Together, these repeats comprise 0.25% of the bunching onion genome with 160,000 copies/1 C of HAT58 and 93,000 copies/1 C of CAT36. Fluorescent in situ hybridization (FISH) and C-banding showed that HAT58 and CAT36 associated with the interstitial and pericentromeric heterochromatin of the A. fistulosum chromosomes 5, 6, 7, and 8. FISH with HAT58 and CAT36 performed on A. cepa (2n = 2C = 16) and A. wakegi (2n = 2C = 16), a natural allodiploid hybrid between A. fistulosum and A. cepa, revealed that these repeats are species specific and produced specific hybridization patterns only on A. fistulosum chromosomes. Thus, the markers can be used in interspecific breeding programs for monitoring of alien genetic material. We applied Non-denaturing FISH that allowed detection of the repeat bearing chromosomes within 3 h. A polymorphism of the HAT58 chromosome location was observed. This finding suggests that the rapid evolution of the HAT58 repeat is still ongoing.

Transcription of highly repetitive tandemly organized DNA in amphibians and birds: A historical overview and modern concepts

Tandemly organized highly repetitive DNA sequences are crucial structural and functional elements of eukaryotic genomes. Despite extensive evidence, satellite DNA remains an enigmatic part of the eukaryotic genome, with biological role and significance of tandem repeat transcripts remaining rather obscure. Data on tandem repeats transcription in amphibian and avian model organisms is fragmentary despite their genomes being thoroughly characterized. Review systematically covers historical and modern data on transcription of amphibian and avian satellite DNA in somatic cells and during meiosis when chromosomes acquire special lampbrush form. We highlight how transcription of tandemly repetitive DNA sequences is organized in interphase nucleus and on lampbrush chromosomes. We offer LTR-activation hypotheses of widespread satellite DNA transcription initiation during oogenesis. Recent explanations are provided for the significance of high-yield production of non-coding RNA derived from tandemly organized highly repetitive DNA. In many cases the data on the transcription of satellite DNA can be extrapolated from lampbrush chromosomes to interphase chromosomes. Lampbrush chromosomes with applied novel technical approaches such as superresolution imaging, chromosome microdissection followed by high-throughput sequencing, dynamic observation in life-like conditions provide amazing opportunities for investigation mechanisms of the satellite DNA transcription.

Chromatin-associated transcripts of tandemly repetitive DNA sequences revealed by RNA-FISH

Tandemly repetitive DNA sequences, also named satellite repeats, are major DNA components of heterochromatin and are often organized as long arrays in the pericentromeric, centromeric, and subtelomeric regions of eukaryotic chromosomes. An increasing amount of evidence indicates that transcripts derived from some satellite repeats play important roles in various biological functions. We used a RNA-fluorescence in situ hybridization (RNA-FISH) technique to investigate the transcription of the four well-characterized satellite repeats of maize (Zea mays), including the 180-bp knob repeat, the telomeric (TTTAGGG)_n repeat, the 156-bp centromeric repeat CentC, and a 350-bp subtelomeric repeat. Although few transcripts derived from these four repeats were found in the expressed sequence tag and RNA-seq databases, RNA-FISH consistently detected the transcripts from three of the four repeats on interphase nuclei, suggesting that the transcripts from the three repeats are largely integrated into chromatin. The transcripts from the knob and telomeric repeats were mapped to the related DNA loci. In contrast, the transcripts from the CentC repeats were mainly localized to the nucleolus, although nucleoplasmic CentC transcripts were also detectable. The nucleolus and nuclear RNAs appeared to be important for the nuclear localization for at least one centromeric protein, Mis12. We demonstrate that RNA-FISH is a powerful tool to assess the level of transcription as well as to physically map the nuclear locations of the transcripts derived from satellite repeats.

The NnCenH3 protein and centromeric DNA sequence profiles of Nelumbo nucifera Gaertn. (sacred lotus) reveal the DNA structures and dynamics of centromeres in basal eudicots

Centromeres on eukaryotic chromosomes consist of large arrays of DNA repeats that undergo very rapid evolution. Nelumbo nucifera Gaertn. (sacred lotus) is a phylogenetic relict and an aquatic perennial basal eudicot. Studies concerning the centromeres of this basal eudicot species could provide ancient evolutionary perspectives. In this study, we characterized the centromeric marker protein NnCenH3 (sacred lotus centromere-specific histone H3 variant), and used a chromatin immunoprecipitation (ChIP)-based technique to recover the NnCenH3 nucleosome-associated sequences of sacred lotus. The properties of the centromere-binding protein and DNA sequences revealed notable divergence between sacred lotus and other flowering plants, including the following factors: (i) an NnCenH3 alternative splicing variant comprising only a partial centromere-targeting domain, (ii) active genes with low transcription levels in the NnCenH3 nucleosomal regions, and (iii) the prevalence of the Ty1/copia class of long terminal repeat (LTR) retrotransposons in the centromeres of sacred lotus chromosomes. In addition, the dynamic natures of the centromeric region showed that some of the centromeric repeat DNA sequences originated from telomeric repeats, and a pair of centromeres on the dicentric chromosome 1 was inactive in the metaphase cells of sacred lotus. Our characterization of the properties of centromeric DNA structure within the sacred lotus genome describes a centromeric profile in ancient basal eudicots and might provide evidence of the origins and evolution of centromeres. Furthermore, the identification of centromeric DNA sequences is of great significance for the assembly of the sacred lotus genome.

Distribution of 45S rDNA in Modern Rose Cultivars (Rosa hybrida), Rosa rugosa, and Their Interspecific Hybrids Revealed by Fluorescence in situ Hybridization

To elucidate the evolutionary dynamics of the location and number of rDNA loci in the process of polyploidization in the genus Rosa, we examined 45S rDNA sites in the chromosomes of 6 modern rose cultivars (R. hybrida), 5 R. rugosa cultivars, and 20 hybrid progenies by fluorescence in situ hybridization. Variation in the number of rDNA sites in parents and their interspecific hybrids was detected. As expected, 4 rDNA sites were observed in the genomes of 4 modern rose cultivars, while 3 hybridization sites were observed in the 2 others. Two expected rDNA sites were found in 2 R. rugosa cultivars, while in the other 3 R. rugosa cultivars 4 sites were present. Among the 20 R. hybrida × R. rugosa offspring, 13 carried the expected number of rDNA sites, and 1 had 6 hybridization sites, which exceeded the expected number by far. The other 6 offspring had either 2 or 3 hybridization sites, which was less than expected. Differences in the number of rDNA loci were observed in interspecific offspring, indicating that rDNA loci exhibit instability after distant hybridization events. Abnormal chromosome pairing may be the main factor explaining the variation in rDNA sites during polyploidization.

High-throughput analysis of the satellitome illuminates satellite DNA evolution

Satellite DNA (satDNA) is a major component yet the great unknown of eukaryote genomes and clearly underrepresented in genome sequencing projects. Here we show the high-throughput analysis of satellite DNA content in the migratory locust by means of the bioinformatic analysis of Illumina reads with the RepeatExplorer and RepeatMasker programs. This unveiled 62 satDNA families and we propose the term "satellitome" for the whole collection of different satDNA families in a genome. The finding that satDNAs were present in many contigs of the migratory locust draft genome indicates that they show many genomic locations invisible by fluorescent in situ hybridization (FISH). The cytological pattern of five satellites showing common descent (belonging to the SF3 superfamily) suggests that non-clustered satDNAs can become into clustered through local amplification at any of the many genomic loci resulting from previous dissemination of short satDNA arrays. The fact that all kinds of satDNA (micro- mini- and satellites) can show the non-clustered and clustered states suggests that all these elements are mostly similar, except for repeat length. Finally, the presence of VNTRs in bacteria, showing similar properties to non-clustered satDNAs in eukaryotes, suggests that this kind of tandem repeats show common properties in all living beings.