Contrasting evolution of a satellite DNA and its ancestral IGS rDNA in Phaseolus (Fabaceae)

Cytomolecular diversity among Vigna Savi (Leguminosae) subgenera

The genus Vigna (Leguminosae) comprises about 150 species grouped into five subgenera. The present study aimed to improve the understanding of karyotype diversity and evolution in Vigna, using new and previously published data through different cytogenetic and DNA content approaches. In the Vigna subgenera, we observed a random distribution of rDNA patterns. The 35S rDNA varied in position, from terminal to proximal, and in number, ranging from one (V. aconitifolia, V. subg. Ceratotropis) to seven pairs (V. unguiculata subsp. unguiculata, V. subg. Vigna). On the other hand, the number of 5S rDNA was conserved (one or two pairs), except for V. radiata (V. subg. Ceratotropis), which had three pairs. Genome size was relatively conserved within the genus, ranging from 1C = 0.43 to 0.70 pg in V. oblongifolia and V. unguiculata subsp. unguiculata, respectively, both belonging to V. subg. Vigna. However, we observed a positive correlation between DNA content and the number of 35S rDNA sites. In addition, data from chromosome-specific BAC-FISH suggest that the ancestral 35S rDNA locus is conserved on chromosome 6 within Vigna. Considering the rapid diversification in the number and position of rDNA sites, such conservation is surprising and suggests that additional sites may have spread out from this ancestral locus.

Chromosomal dynamics in Senna: comparative PLOP-FISH analysis of tandem repeats and flow cytometric nuclear genome size estimations

Introduction

Tandem repeats (TRs) occur abundantly in plant genomes. They play essential roles that affect genome organization and evolution by inducing or generating chromosomal rearrangements such as duplications, deletions, inversions, and translocations. These impact gene expression and chromosome structure and even contribute to the emergence of new species.

Method

We investigated the effects of TRs on speciation in Senna genus by performing a comparative analysis using fluorescence in situ hybridization (FISH) with S. tora-specific TR probes. We examined the chromosomal distribution of these TRs and compared the genome sizes of seven Senna species (estimated using flow cytometry) to better understand their evolutionary relationships.

Results

Two (StoTR03_159 and StoTR04_55) of the nine studied TRs were not detected in any of the seven Senna species, whereas the remaining seven were found in all or some species with patterns that were similar to or contrasted with those of S. tora. Of these studies species, only S. angulata showed significant genome rearrangements and dysploid karyotypes resembling those of S. tora. The genome sizes varied among these species and did not positively correlate with chromosome number. Notably, S. angulata had the fewest chromosomes (2n = 22) but a relatively large genome size.

Discussion

These findings reveal the dynamics of TRs and provide a cytogenetic depiction of chromosomal rearrangements during speciation in Senna. To further elucidate the dynamics of repeat sequences in Senna, future studies must include related species and extensive repeatomic studies, including those on transposable elements.

Study and Physical Mapping of the Species-Specific Tandem Repeat CS-237 Linked with 45S Ribosomal DNA Intergenic Spacer in Cannabis sativa L

Hemp (Cannabis sativa L.) is a valuable crop and model plant for studying sex chromosomes. The scientific interest in the plant has led to its whole genome sequencing and the determination of its cytogenetic characteristics. A range of cytogenetic markers (subtelomeric repeat CS-1, 5S rDNA, and 45S rDNA) has been mapped onto hemp's chromosomes by fluorescent in situ hybridization (FISH). In this study, another cytogenetic marker (the tandem repeat CS-237, with a 237 bp monomer) was found, studied, and localized on chromosomes by FISH. The signal distribution and karyotyping revealed that the CS-237 probe was localized in chromosome 6 with one hybridization site and in chromosome 8 with two hybridization sites, one of which colocalizes with the 45S rDNA probe (with which a nucleolus organizer region, NOR, was detected). A BLAST analysis of the genomic data and PCR experiments showed that the modified CS-237 monomers (delCS-237, 208 bp in size) were present in the intergenic spacers (IGSs) of hemp 45S rDNA monomers. Such a feature was firstly observed in Cannabaceae species. However, IGS-linked DNA repeats were found in several plant species of other families (Fabaceae, Solanaceae, and Asteraceae). This phenomenon is discussed in this article. The example of CS-237 may be useful for further studying the phenomenon as well as for the physical mapping of hemp chromosomes.

The Genomics of Plant Satellite DNA

The twenty-first century began with a certain indifference to the research of satellite DNA (satDNA). Neither genome sequencing projects were able to accurately encompass the study of satDNA nor classic methodologies were able to go further in undertaking a better comprehensive study of the whole set of satDNA sequences of a genome. Nonetheless, knowledge of satDNA has progressively advanced during this century with the advent of new analytical techniques. The enormous advantages that genome-wide approaches have brought to its analysis have now stimulated a renewed interest in the study of satDNA. At this point, we can look back and try to assess more accurately many of the key questions that were left unsolved in the past about this enigmatic and important component of the genome. I review here the understanding gathered on plant satDNAs over the last few decades with an eye on the near future.

Subgenome Discrimination in Brassica and Raphanus Allopolyploids Using Microsatellites

Intergeneric crosses between Brassica species and Raphanus sativus have produced crops with prominent shoot and root systems of Brassica and R. sativus, respectively. It is necessary to discriminate donor genomes when studying cytogenetic stability in distant crosses to identify homologous chromosome pairing, and microsatellite repeats have been used to discriminate subgenomes in allopolyploids. To identify genome-specific microsatellites, we explored the microsatellite content in three Brassica species (B. rapa, AA, B. oleracea, CC, and B. nigra, BB) and R. sativus (RR) genomes, and validated their genome specificity by fluorescence in situ hybridization. We identified three microsatellites showing A, C, and B/R genome specificity. ACBR_msat14 and ACBR_msat20 were detected in the A and C chromosomes, respectively, and ACBR_msat01 was detected in B and R genomes. However, we did not find a microsatellite that discriminated the B and R genomes. The localization of ACBR_msat20 in the 45S rDNA array in ×Brassicoraphanus 977 corroborated the association of the 45S rDNA array with genome rearrangement. Along with the rDNA and telomeric repeat probes, these microsatellites enabled the easy identification of homologous chromosomes. These data demonstrate the utility of microsatellites as probes in identifying subgenomes within closely related Brassica and Raphanus species for the analysis of genetic stability of new synthetic polyploids of these genomes.

Comparative FISH analysis of Senna tora tandem repeats revealed insights into the chromosome dynamics in Senna

BACKGROUND

DNA tandem repeats (TRs) are often abundant and occupy discrete regions in eukaryotic genomes. These TRs often cause or generate chromosomal rearrangements, which, in turn, drive chromosome evolution and speciation. Tracing the chromosomal distribution of TRs could therefore provide insights into the chromosome dynamics and speciation among closely related taxa. The basic chromosome number in the genus Senna is 2n = 28, but dysploid species like Senna tora have also been observed.

OBJECTIVE

To understand the dynamics of these TRs and their impact on S. tora dysploidization.

METHODS

We performed a comparative fluorescence in situ hybridization (FISH) analysis among nine closely related Senna species and compared the chromosomal distribution of these repeats from a cytotaxonomic perspective by using the ITS1-5.8S-ITS2 sequence to infer phylogenetic relationships.

RESULTS

Of the nine S. tora TRs, two did not show any FISH signal whereas seven TRs showed similar and contrasting patterns to other Senna species. StoTR01_86, which was localized in the pericentromeric regions in all S. tora, but not at the nucleolar organizer region (NOR) site, was colocalized at the NOR site in all species except in S. siamea. StoTR02_7_tel was mostly localized at chromosome termini, but some species had an interstitial telomeric repeat in a few chromosomes. StoTR05_180 was distributed in the subtelomeric region in most species and was highly amplified in the pericentromeric region in some species. StoTR06_159 was either absent or colocalized in the NOR site in some species, and StoIGS_463, which was localized at the NOR site in S. tora, was either absent or localized at the subtelomeric or pericentromeric regions in other species.

CONCLUSIONS

These data suggest that TRs play important roles in S. tora dysploidy and suggest the involvement of 45S rDNA intergenic spacers in "carrying" repeats during genome reshuffling.

Chromosomal Mapping of Tandem Repeats Revealed Massive Chromosomal Rearrangements and Insights Into Senna tora Dysploidy

Tandem repeats can occupy a large portion of plant genomes and can either cause or result from chromosomal rearrangements, which are important drivers of dysploidy-mediated karyotype evolution and speciation. To understand the contribution of tandem repeats in shaping the extant Senna tora dysploid karyotype, we analyzed the composition and abundance of tandem repeats in the S. tora genome and compared the chromosomal distribution of these repeats between S. tora and a closely related euploid, Senna occidentalis. Using a read clustering algorithm, we identified the major S. tora tandem repeats and visualized their chromosomal distribution by fluorescence in situ hybridization. We identified eight independent repeats covering ~85 Mb or ~12% of the S. tora genome. The unit lengths and copy numbers had ranges of 7-5,833 bp and 325-2.89 × 10⁶, respectively. Three short duplicated sequences were found in the 45S rDNA intergenic spacer, one of which was also detected at an extra-NOR locus. The canonical plant telomeric repeat (TTTAGGG)_n was also detected as very intense signals in numerous pericentromeric and interstitial loci. StoTR05_180, which showed subtelomeric distribution in Senna occidentalis, was predominantly pericentromeric in S. tora. The unusual chromosomal distribution of tandem repeats in S. tora not only enabled easy identification of individual chromosomes but also revealed the massive chromosomal rearrangements that have likely played important roles in shaping its dysploid karyotype.

Sequence Composition Underlying Centromeric and Heterochromatic Genome Compartments of the Pacific Oyster Crassostrea gigas

Segments of the genome enriched in repetitive sequences still present a challenge and are omitted in genome assemblies. For that reason, the exact composition of DNA sequences underlying the heterochromatic regions and the active centromeres are still unexplored for many organisms. The centromere is a crucial region of eukaryotic chromosomes responsible for the accurate segregation of genetic material. The typical landmark of centromere chromatin is the rapidly-evolving variant of the histone H3, CenH3, while DNA sequences packed in constitutive heterochromatin are associated with H3K9me3-modified histones. In the Pacific oyster Crassostrea gigas we identified its centromere histone variant, Cg-CenH3, that shows stage-specific distribution in gonadal cells. In order to investigate the DNA composition of genomic regions associated with the two specific chromatin types, we employed chromatin immunoprecipitation followed by high-throughput next-generation sequencing of the Cg-CenH3- and H3K9me3-associated sequences. CenH3-associated sequences were assigned to six groups of repetitive elements, while H3K9me3-associated-ones were assigned only to three. Those associated with CenH3 indicate the lack of uniformity in the chromosomal distribution of sequences building the centromeres, being also in the same time dispersed throughout the genome. The heterochromatin of C. gigas exhibited general paucity and limited chromosomal localization as predicted, with H3K9me3-associated sequences being predominantly constituted of DNA transposons.

Diversity of repetitive sequences within compact genomes of Phaseolus L. beans and allied genera Cajanus L. and Vigna Savi

Repetitive sequences are ubiquitous and fast-evolving elements responsible for size variation and large-scale organization of plant genomes. Within tribe Phaseoleae (Fabaceae), some genera, such as Phaseolus, Vigna, and Cajanus, show small genome and mostly stable chromosome number. Here, we applied a combined computational and cytological approach to study the organization and diversification of repetitive elements in some species of these genera. Sequences were classified in terms of type and repetitiveness and the most abundant were mapped to chromosomes. We identified long terminal repeat (LTR) retrotransposons, especially Ogre and Chromovirus elements, making up most of genomes, other than P. acutifolius and Vigna species. Satellite DNAs (SatDNAs) were less representative, but highly diverse among species, showing a clear phylogenetic relationship. In situ localization revealed preferential location at pericentromeres and centromeres for both types of sequences, suggesting a heterogeneous composition, especially for centromeres. Few elements showed subterminal accumulation. Copy number variation among chromosomes within and among species was observed for all nine identified SatDNAs. Altogether, our data pointed two main elements (Ty3/Gypsy retrotransponsons and SatDNAs) to the diversification on the repetitive landscape in Phaseoleae, with a typical set of repeats in each species. The high turnover of these sequences, however, did not affect total genome size.

Pilot satellitome analysis of the model plant, Physcomitrellapatens, revealed a transcribed and high-copy IGS related tandem repeat

Satellite DNA (satDNA) constitutes a substantial part of eukaryotic genomes. In the last decade, it has been shown that satDNA is not an inert part of the genome and its function extends beyond the nuclear membrane. However, the number of model plant species suitable for studying the novel horizons of satDNA functionality is low. Here, we explored the satellitome of the model "basal" plant, Physcomitrellapatens (Hedwig, 1801) Bruch & Schimper, 1849 (moss), which has a number of advantages for deep functional and evolutionary research. Using a newly developed pyTanFinder pipeline (https://github.com/Kirovez/pyTanFinder) coupled with fluorescence in situ hybridization (FISH), we identified five high copy number tandem repeats (TRs) occupying a long DNA array in the moss genome. The nuclear organization study revealed that two TRs had distinct locations in the moss genome, concentrating in the heterochromatin and knob-rDNA like chromatin bodies. Further genomic, epigenetic and transcriptomic analysis showed that one TR, named PpNATR76, was located in the intergenic spacer (IGS) region and transcribed into long non-coding RNAs (lncRNAs). Several specific features of PpNATR76 lncRNAs make them very similar with the recently discovered human lncRNAs, raising a number of questions for future studies. This work provides new resources for functional studies of satellitome in plants using the model organism P.patens, and describes a list of tandem repeats for further analysis.

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.

Satellite DNA in Plants: More than Just Rubbish

For decades, satellite DNAs have been the hidden part of genomes. Initially considered as junk DNA, there is currently an increasing appreciation of the functional significance of satellite DNA repeats and of their sequences. Satellite DNA families accumulate in the heterochromatin in different parts of the eukaryotic chromosomes, mainly in pericentromeric and subtelomeric regions, but they also span the functional centromere. Tandem repeat sequences may spread from subtelomeric to interstitial loci, leading to the formation of chromosome-specific loci or to the accumulation in equilocal sites in different chromosomes. They also appear as the main components of the heterochromatin in the sex-specific region of sex chromosomes. Satellite DNA, required for chromosome organization, also plays a role in pairing and segregation. Some satellite repeats are transcribed and can participate in the formation and maintenance of heterochromatin structure and in the modulation of gene expression. In addition to the identification of the different satellite DNA families, their characteristics and location, we are interested in determining their impact on the genomes, by identifying the mechanisms leading to their appearance and amplification as well as in understanding how they change over time, the factors affecting these changes, and the influence exerted by the evolutionary history of the organisms. On the other hand, satellite DNA sequences are rapidly evolving sequences that may cause reproductive barriers between organisms and promote speciation. The accumulation of experimental data collected in recent years and the emergence of new approaches based on next-generation sequencing and high-throughput genome analysis are opening new perspectives that are changing our understanding of satellite DNA. This review examines recent data to provide a timely update on the overall information gathered about this part of the genome, focusing on the advances in the knowledge of its origin, its evolution, and its potential functional roles.

Differential amplification of satellite PaB6 in chromosomally hypervariable Prospero autumnale complex (Hyacinthaceae)

BACKGROUND AND AIMS

Chromosomal evolution, including numerical and structural changes, is a major force in plant diversification and speciation. This study addresses genomic changes associated with the extensive chromosomal variation of the Mediterranean Prospero autumnale complex (Hyacinthaceae), which includes four diploid cytotypes each with a unique combination of chromosome number (x = 5, 6, 7), rDNA loci and genome size.

METHODS

A new satellite repeat PaB6 has previously been identified, and monomers were reconstructed from next-generation sequencing (NGS) data of P. autumnale cytotype B(6)B(6) (2n = 12). Monomers of all other Prospero cytotypes and species were sequenced to check for lineage-specific mutations. Copy number, restriction patterns and methylation levels of PaB6 were analysed using Southern blotting. PaB6 was localized on chromosomes using fluorescence in situ hybridization (FISH).

KEY RESULTS

The monomer of PaB6 is 249 bp long, contains several intact and truncated vertebrate-type telomeric repeats and is highly methylated. PaB6 is exceptional because of its high copy number and unprecedented variation among diploid cytotypes, ranging from 10(4) to 10(6) copies per 1C. PaB6 is always located in pericentromeric regions of several to all chromosomes. Additionally, two lineages of cytotype B(7)B(7) (x = 7), possessing either a single or duplicated 5S rDNA locus, differ in PaB6 copy number; the ancestral condition of a single locus is associated with higher PaB6 copy numbers.

CONCLUSIONS

Although present in all Prospero species, PaB6 has undergone differential amplification only in chromosomally variable P. autumnale, particularly in cytotypes B(6)B(6) and B(5)B(5). These arose via independent chromosomal fusions from x = 7 to x = 6 and 5, respectively, accompanied by genome size increases. The copy numbers of satellite DNA PaB6 are among the highest in angiosperms, and changes of PaB6 are exceptionally dynamic in this group of closely related cytotypes of a single species. The evolution of the PaB6 copy numbers is discussed, and it is suggested that PaB6 represents a recent and highly dynamic system originating from a small pool of ancestral repeats.

Next generation sequencing analysis reveals a relationship between rDNA unit diversity and locus number in Nicotiana diploids

BACKGROUND

Tandemly arranged nuclear ribosomal DNA (rDNA), encoding 18S, 5.8S and 26S ribosomal RNA (rRNA), exhibit concerted evolution, a pattern thought to result from the homogenisation of rDNA arrays. However rDNA homogeneity at the single nucleotide polymorphism (SNP) level has not been detailed in organisms with more than a few hundred copies of the rDNA unit. Here we study rDNA complexity in species with arrays consisting of thousands of units.

METHODS

We examined homogeneity of genic (18S) and non-coding internally transcribed spacer (ITS1) regions of rDNA using Roche 454 and/or Illumina platforms in four angiosperm species, Nicotiana sylvestris, N. tomentosiformis, N. otophora and N. kawakamii. We compared the data with Southern blot hybridisation revealing the structure of intergenic spacer (IGS) sequences and with the number and distribution of rDNA loci.

RESULTS AND CONCLUSIONS

In all four species the intragenomic homogeneity of the 18S gene was high; a single ribotype makes up over 90% of the genes. However greater variation was observed in the ITS1 region, particularly in species with two or more rDNA loci, where >55% of rDNA units were a single ribotype, with the second most abundant variant accounted for >18% of units. IGS heterogeneity was high in all species. The increased number of ribotypes in ITS1 compared with 18S sequences may reflect rounds of incomplete homogenisation with strong selection for functional genic regions and relaxed selection on ITS1 variants. The relationship between the number of ITS1 ribotypes and the number of rDNA loci leads us to propose that rDNA evolution and complexity is influenced by locus number and/or amplification of orphaned rDNA units at new chromosomal locations.