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.

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.

No abstract