Molecular structure and chromosome distribution of three repetitive DNA families in Anemone hortensis L. (Ranunculaceae)

Mega-sized pericentromeric blocks of simple telomeric repeats and their variants reveal patterns of chromosome evolution in ancient Cycadales genomes

Simple telomeric repeats composed of six to seven iterating nucleotide units are important sequences typically found at the ends of chromosomes. Here we analyzed their abundance and homogeneity in 42 gymnosperm (29 newly sequenced), 29 angiosperm (one newly sequenced), and eight bryophytes using bioinformatics, conventional cytogenetic and molecular biology approaches to explore their diversity across land plants. We found more than 10 000-fold variation in the amounts of telomeric repeats among the investigated taxa. Repeat abundance was positively correlated with increasing intragenomic sequence heterogeneity and occurrence at non-telomeric positions, but there was no correlation with genome size. The highest abundance/heterogeneity was found in the gymnosperm genus Cycas (Cycadaceae), in which megabase-sized blocks of telomeric repeats (i.e., billions of copies) were identified. Fluorescent in situ hybridization experiments using variant-specific probes revealed canonical Arabidopsis-type telomeric TTTAGGG repeats at chromosome ends, while pericentromeric blocks comprised at least four major telomeric variants with decreasing abundance: TTTAGGG>TTCAGGG >TTTAAGG>TTCAAGG. Such a diversity of repeats was not found in the sister cycad family Zamiaceae or in any other species analyzed. Using immunocytochemistry, we showed that the pericentromeric blocks of telomeric repeats overlapped with histone H3 serine 10 phosphorylation signals. We show that species of Cycas have amplified their telomeric repeats in centromeric and telomeric positions on telocentric chromosomes to extraordinary high levels. The ancestral chromosome number reconstruction suggests their occurrence is unlikely to be the product of ancient Robertsonian chromosome fusions. We speculate as to how the observed chromosome dynamics may be associated with the diversification of cycads.

Structure and Methylation of 35S rDNA in Allopolyploids Anemone multifida (2n = 4x = 32, BBDD) and Anemone baldensis (2n = 6x = 48, AABBDD) and Their Parental Species Show Evidence of Nucleolar Dominance

Transcriptional silencing of 35S rDNA loci inherited from one parental species is occurring relatively frequently in allopolyploids. However, molecular mechanisms by which it is selected for transcriptional silencing remain unclear. We applied NGS, silver staining and bisulfite sequencing to study the structure, expression and methylation landscape of 35S rDNA in two allopolyploids of common origin, allotetraploid Anemone multifida (2n = 4x = 32, genome composition BBDD) and allohexaploid A. baldensis (2n = 6x = 48, AABBDD), and their genome donors, A. sylvestris (2n = 16, AA), A. cylindrica (2n = 16, BB) and A. parviflora (2n = 16, DD). The size of the recovered 35S rDNA units varied from 10,489 bp in A. cylindrica to 12,084 bp in A. sylvestris. Anemone showed an organization typical of most ribosomal 35S rDNA composed of NTS, ETS, rRNA genes, TTS and TIS with structural features of plant IGS sequences and all functional elements needed for rRNA gene activity. The NTS was more variable than the ETS and consisted of SRs which are highly variable among Anemone. Five to six CpG-rich islands were found within the ETS. CpG island located adjacent to the transcription initiation site (TIS) was highly variable regarding the sequence size and methylation level and exhibited in most of the species lower levels of methylation than CpG islands located adjacent to the 18S rRNA gene. Our results uncover hypomethylation of A. sylvestris- and A. parviflora-derived 35S rDNA units in allopolyploids A. multifida and A. baldensis. Hypomethylation of A. parviflora-derived 35S rDNA was more prominent in A. baldensis than in A. multifida. We showed that A. baldensis underwent coupled A. sylvestris-derived 35S rDNA array expansion and A. parviflora-derived 35S rDNA copy number decrease that was accompanied by lower methylation level of A. sylvestris-derived 35S rDNA units in comparison to A. parviflora-derived 35S rDNA units. These observations suggest that in A. baldensis nucleolar dominance is directed toward A. sylvestris-derived chromosomes. This work broadens our current knowledge of the 35S rDNA organization in Anemone and provides evidence of the progenitor-specific 35S rDNA methylation in nucleolar dominance.

The Pontastacus leptodactylus (Astacidae) Repeatome Provides Insight Into Genome Evolution and Reveals Remarkable Diversity of Satellite DNA

Pontastacus leptodactylus is a native European crayfish species found in both freshwater and brackish environments. It has commercial importance for fisheries and aquaculture industries. Up till now, most studies concerning P. leptodactylus have focused onto gaining knowledge about its phylogeny and population genetics. However, little is known about the chromosomal evolution and genome organization of this species. Therefore, we performed clustering analysis of a low coverage genomic dataset to identify and characterize repetitive DNA in the P. leptodactylus genome. In addition, the karyogram of P. leptodactylus (2n = 180) is presented here for the first time consisting of 75 metacentric, 14 submetacentric, and a submetacentric/metacentric heteromorphic chromosome pair. We determined the genome size to be at ~18.7 gigabase pairs. Repetitive DNA represents about 54.85% of the genome. Satellite DNA repeats are the most abundant type of repetitive DNA, making up to ~28% of the total amount of repetitive elements, followed by the Ty3/Gypsy retroelements (~15%). Our study established a surprisingly high diversity of satellite repeats in P. leptodactylus. The genome of P. leptodactylus is by far the most satellite-rich genome discovered to date with 258 satellite families described. Of the five mapped satellite DNA families on chromosomes, PlSAT3-411 co-localizes with the AT-rich DAPI positive probable (peri)centromeric heterochromatin on all chromosomes, while PlSAT14-79 co-localizes with the AT-rich DAPI positive (peri)centromeric heterochromatin on one chromosome and is also located subterminally and intercalary on some chromosomes. PlSAT1-21 is located intercalary in the vicinity of the (peri)centromeric heterochromatin on some chromosomes, while PlSAT6-70 and PlSAT7-134 are located intercalary on some P. leptodactylus chromosomes. The FISH results reveal amplification of interstitial telomeric repeats (ITRs) in P. leptodactylus. The prevalence of repetitive elements, especially the satellite DNA repeats, may have provided a driving force for the evolution of the P. leptodactylus genome.

Bioinformatic and Molecular Analysis of Satellite Repeat Diversity in Vaccinium Genomes

Bioinformatic and molecular characterization of satellite repeats was performed to understand the impact of their diversification on Vaccinium genome evolution. Satellite repeat diversity was evaluated in four cultivated and wild species, including the diploid species Vaccinium myrtillus and Vaccinium uliginosum, as well as the tetraploid species Vaccinium corymbosum and Vaccinium arctostaphylos. We comparatively characterized six satellite repeat families using in total 76 clones with 180 monomers. We observed that the monomer units of VaccSat1, VaccSat2, VaccSat5, and VaccSat6 showed a higher order repeat (HOR) structure, likely originating from the organization of two adjacent subunits with differing similarity, length and size. Moreover, VaccSat1, VaccSat3, VaccSat6, and VaccSat7 were found to have sequence similarity to parts of transposable elements. We detected satellite-typical tandem organization for VaccSat1 and VaccSat2 in long arrays, while VaccSat5 and VaccSat6 distributed in multiple sites over all chromosomes of tetraploid V. corymbosum, presumably in long arrays. In contrast, very short arrays of VaccSat3 and VaccSat7 are dispersedly distributed over all chromosomes in the same species, likely as internal parts of transposable elements. We provide a comprehensive overview on satellite species specificity in Vaccinium, which are potentially useful as molecular markers to address the taxonomic complexity of the genus, and provide information for genome studies of this genus.

The Repetitive DNA Composition in the Natural Pesticide Producer Tanacetum cinerariifolium: Interindividual Variation of Subtelomeric Tandem Repeats

Dalmatian pyrethrum (Tanacetum cinerariifolium (Trevir.) Sch. Bip.), a plant species endemic to the east Adriatic coast, is used worldwide for production of the organic insecticide, pyrethrin. Most studies concerning Dalmatian pyrethrum have focused on its morphological and biochemical traits relevant for breeding. However, little is known about the chromosomal evolution and genome organization of this species. Our study aims are to identify, classify, and characterize repetitive DNA in the T. cinerariifolium genome using clustering analysis of a low coverage genomic dataset. Repetitive DNA represents about 71.63% of the genome. T. cinerariifolium exhibits linked 5S and 35S rDNA configuration (L-type). FISH reveals amplification of interstitial telomeric repeats (ITRs) in T. cinerariifolium. Of the three newly identified satellite DNA families, TcSAT1 and TcSAT2 are located subterminally on most of T. cinerariifolium chromosomes, while TcSAT3 family is located intercalary within the longer arm of two chromosome pairs. FISH reveals high levels of polymorphism of the TcSAT1 and TcSAT2 sites by comparative screening of 28 individuals. TcSAT2 is more variable than TcSAT1 regarding the number and position of FISH signals. Altogether, our data highlights the dynamic nature of DNA sequences associated with subtelomeres in T. cinerariifolium and suggests that subtelomeres represent one of the most dynamic and rapidly evolving regions in eukaryotic genomes.

Telomeres and genomic evolution

The terminal regions of eukaryotic chromosomes, composed of telomere repeat sequences and sub-telomeric sequences, represent some of the most variable and rapidly evolving regions of the genome. The sub-telomeric regions are characterized by segmentally duplicated repetitive DNA elements, interstitial telomere repeat sequences and families of variable genes. Sub-telomeric repeat sequence families are shared among multiple chromosome ends, often rendering detailed sequence characterization difficult. These regions are composed of constitutive heterochromatin and are subjected to high levels of meiotic recombination. Dysfunction within telomere repeat arrays, either due to disruption in the chromatin structure or because of telomere shortening, can lead to chromosomal fusion and the generation of large-scale genomic rearrangements across the genome. The dynamic nature of telomeric regions, therefore, provides functionally useful variation to create genetic diversity, but also provides a mechanism for rapid genomic evolution that can lead to reproductive isolation and speciation. This article is part of the theme issue 'Understanding diversity in telomere dynamics'.This article is part of the theme issue 'Understanding diversity in telomere dynamics'.

The impact of the Tekay chromoviral elements on genome organisation and evolution of Anemone s.l. (Ranunculaceae)

We studied the highly abundant chromoviral Tekay clade in species from three sister genera - Anemone, Pulsatilla and Hepatica (Ranunculaceae). With this clade, we performed a concomitant survey of its phylogenetic diversity, chromosomal organisation and transcriptional activity in Anemone s.l. in order to investigate dynamics of the Tekay elements at a finer scale than previously achieved in this or any other flowering clade. The phylogenetic tree built from Tekay sequences conformed to expected evolutionary relationships of the species; exceptions being A. nemorosa and A. sylvestris, which appeared more closely related that expected, and we invoke hybridisation events to explain the observed topology. The separation of elements into six clusters could be explained by episodic bursts of activity since divergence from a common ancestor at different points in their respective evolutionary histories. In Anemone s.l. the Tekay elements do not have a preferential position on chromosomes, i.e. they can have a: (i) centromeric/pericentromeric position; (ii) interstitial position in DAPI-positive AT-rich heterochromatic regions; can be (iii) dispersed throughout chromosomes; or even (iv) be absent from large heterochromatic blocks. Widespread transcriptional activity of the Tekay elements in Anemone s.l. taxa indicate that some copies of Tekay elements could still be active in this plant group, contributing to genome evolution and speciation within Anemone s.l. Identification of Tekay elements in Anemone s.l. provides valuable information for understanding how different localisation patterns might help to facilitate plant genome organisation in a structural and functional manner.

Towards a FISH-based karyotype of Rosa L. (Rosaceae)

The genus Rosa Linnaeus, 1753 has important economic value in ornamental sector and many breeding activities are going on supported by molecular studies. However, the cytogenetic studies of rose species are scarce and mainly focused on chromosome counting and chromosome morphology-based karyotyping. Due to the small size of the chromosomes and a high frequency of polyploidy in the genus, karyotyping is very challenging for rose species and requires FISH-based cytogenetic markers to be applied. Therefore, in this work the aim is to establish a FISH-based karyotype for Rosa wichurana (Crépin, 1888), a rose species with several benefits for advanced molecular cytogenetic studies of genus Rosa (Kirov et al. 2015a). It is shown that FISH signals from 5S, 45S and an Arabidopsis-type telomeric repeat are distributed on five (1, 2, 4, 5 and 7) of seven chromosome pairs. In addition, it is demonstrated that the interstitial telomeric repeat sequences (ITR) are located in the centromeric regions of four chromosome pairs. Using low hybridization stringency for ITR visualization, we showed that the number of ITR signals increases four times (1-4 signals). This study is the first to propose a FISH-based Rosa wichurana karyotype for the reliable identification of chromosomes. The possible origin of Rosa wichurana ITR loci is discussed.

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.

FISH using a gag-like fragment probe reveals a common Ty3-gypsy-like retrotransposon in genome of Coffea species

The genus Coffea possesses about 100 species, and the most economically important are Coffea canephora and Coffea arabica. The latter is predominantly self-compatible with 2n = 4x = 44, while the others of the genus are diploid with 2n = 2x = 22 and mostly self-incompatible. Studies using molecular markers have been useful to detect differences between genomes in Coffea; however, molecular and cytogenetic studies have produced only limited information on the karyotypes organization. We used DOP-PCR to isolate repetitive elements from genome of Coffea arabica var. typica. The pCa06 clone, containing a fragment of 775 bp length, was characterized by sequencing and used as a probe in chromosomes of C. arabica and six other species: C. canephora, Coffea eugenioides, Coffea kapakata, Coffea liberica var. dewevrei, Coffea racemosa, and Coffea stenophylla. This insert shows similarities with a gag protein of the Ty3-gypsy-like super-family. Dot blot and FISH analyses demonstrated that pCa06 is differentially accumulated between species and chromosomes. Signals appeared scattered and clustered on the chromosomes and were also associated with heterochromatic regions. While the literature shows that there is a high karyotype similarity between Coffea species, our results point out differences in the accumulation and dispersion of this Ty3-gypsy-like retrotransposon during karyotype differentiation of Coffea.

Evolution of the tetraploid Anemone multifida (2n = 32) and hexaploid A. baldensis (2n = 48) (Ranunculaceae) was accompanied by rDNA loci loss and intergenomic translocation: evidence for their common genome origin

BACKGROUND AND AIMS

In the genus Anemone two small groups of taxa occur with the highest ploidy levels 2n = 6x = 48, belonging to the closely related clades: the montane/alpine Baldensis clade and the more temperate Multifida clade. To understand the formation of polyploids within these groups, the evolution of allohexaploid A. baldensis (AABBDD, 2n = 6x = 48) from Europe and allotetraploid Anemone multifida (BBDD, 2n = 4x = 32) from America was analysed.

METHODS

Internal transcribed spacer and non-transcribed spacer sequences were used as molecular markers for phylogenetic analyses. Cytogenetic studies, including genomic in situ hybridization with genomic DNA of potential parental species as probe, fluorescence in situ hybridization with 5S and 18S rDNA as probes and 18S rDNA restriction analyses, were used to identify the parental origin of chromosomes and to study genomic changes following polyploidization.

KEY RESULTS

This study shows that A. multifida (BBDD, 2n= 4x = 32) and A. baldensis (AABBDD, 2n = 6x = 48) are allopolyploids originating from the crosses of diploid members of the Multifida (donor of the A and B subgenomes) and Baldensis groups (donor of the D subgenome). The A and B subgenomes are closely related to the genomes of A. sylvestris, A. virginiana and A. cylindrica, indicating that these species or their progeny might be the ancestral donors of the B subgenome of A. multifida and A and B subgenomes of A. baldensis. Both polyploids have undergone genomic changes such as interchromosomal translocation affecting B and D subgenomes and changes at rDNA sites. Anemone multifida has lost the 35S rDNA loci characteristic of the maternal donor (B subgenome) and maintained only the rDNA loci of the paternal donor (D subgenome).

CONCLUSIONS

It is proposed that A. multifida and A. baldensis probably had a common ancestor and their evolution was facilitated by vegetation changes during the Quaternary, resulting in their present disjunctive distribution.

Cytogenetic and phylogenetic studies of diploid and polyploid members of tribe Anemoninae (Ranunculaceae)

The ancestry, phylogenetic differentiation and systematic classification of the worldwide-distributed genus Anemone have been debated for many years. In this paper 11 Anemone, three Pulsatilla species and Hepatica nobilis were subjected to detailed karyotype analysis with the aim of obtaining new cytogenetic data that will contribute to karyotype evolutionary studies of the tribe Anemoninae. The results are interpreted in a phylogenetic context, established from the intergenic nontranscribed spacer (NTS) of 5S rDNA and internal transcribed spacer (ITS) of 35S rDNA. One to three 35S and one to three 5S rDNA loci are present in diploid and polyploid taxa. The 35S rDNA loci are located terminally on the short arm of acrocentric chromosomes, while for 5S rDNA there is no preferential chromosomal position as it exhibits terminal, subterminal, interstitial or pericentromeric positions, and is located either on acrocentric or metacentric chromosomes. The karyotype of hexaploid A. baldensis (2n = 6x = 48) is presented for the first time, and A. sylvestris is proposed as one of its putative parental species. Chromosome fusion/translocation is proposed as the key mechanism involved in reduction of the basic chromosome number from 8 in the Anemone subgenus to 7 in the Anemonidium subgenus. The cytogenetic data obtained are mainly supported by ITS and NTS phylogeny. Diversification of the genus Anemone was accompanied by a large reduction of heterochromatin, from the Mediterranean anemones that have large amounts of heterochromatin to the New World anemones without any detectable heterochromatic blocks.