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

Distinct structural variants and repeat landscape shape the genomes of the ancient grapes Aglianico and Falanghina

Mounting evidence recognizes structural variations (SVs) and repetitive DNA sequences as crucial players in shaping the existing grape phenotypic diversity at intra- and inter-species levels. To deepen our understanding on the abundance, diversity, and distribution of SVs and repetitive DNAs, including transposable elements (TEs) and tandemly repeated satellite DNA (satDNAs), we re-sequenced the genomes of the ancient grapes Aglianico and Falanghina. The analysis of large copy number variants (CNVs) detected candidate polymorphic genes that are involved in the enological features of these varieties. In a comparative analysis of Aglianico and Falanghina sequences with 21 publicly available genomes of cultivated grapes, we provided a genome-wide annotation of grape TEs at the lineage level. We disclosed that at least two main clusters of grape cultivars could be identified based on the TEs content. Multiple TEs families appeared either significantly enriched or depleted. In addition, in silico and cytological analyses provided evidence for a diverse chromosomal distribution of several satellite repeats between Aglianico, Falanghina, and other grapes. Overall, our data further improved our understanding of the intricate grape diversity held by two Italian traditional varieties, unveiling a pool of unique candidate genes never so far exploited in breeding for improved fruit quality.

Differential amplification of the subtelomeric satellite DNA JcSAT1 in the genus Jatropha L. (Euphorbiaceae)

Satellite DNAs (satDNAs) are highly repetitive sequences that occur in virtually all eukaryotic genomes and can undergo rapid copy number and nucleotide sequence variation among relatives. After chromosomal mapping of the satDNA JcSAT1, it was found a large accumulation at subtelomeres of Jatropha curcas (subgenus Curcas), but an absence of these monomers in J. integerrima (subgenus Jatropha). This fact suggests a dynamic scenario for this satellite repeat in Jatropha genomes. Here, we used a multitasking approach (sequence analysis, DNA blotting and chromosomal mapping) to investigate the molecular organization and chromosomal abundance and distribution of JcSAT1 in a broader group of species from the subgenus Jatropha (J. gossypiifolia, J. mollissima, J. podagrica, and J. multifida) in addition to J. curcas, with the aiming of understanding the evolution of this satDNA. Based on the analysis of BAC clone sequences of J. curcas, a large array (~ 30 kb) of 80 homogeneous monomers of JcSAT1 was identified in BAC 23J11. The monomer size was conserved (~ 358 bp) and contained a telomeric motif at the 5' end. PCR amplification coupled with a Southern blot revealed the presence of JcSAT1-like sequences in all species examined. However, a large set of genome copies was identified only in J. curcas, where a ladder-like pattern with multimers of different sizes was observed. In situ hybridization of BAC 23J11 confirmed the subtelomeric pattern for J. curcas, but showed no signals on chromosomes of species from the subgenus Jatropha. Our data indicate that JcSAT1 is a highly homogeneous satDNA that originated from a region near the telomeres and spread throughout the chromosomal subtermini, possibly due to frequent ectopic recombination between these regions. The abundance of JcSAT1 in the genome of J. curcas suggests that an amplification event occurred either at the base of the subgenus Curcas or at least in this species, although the repeat is shared by all species of the genus studied so far.

The nature and organization of satellite DNAs in Petunia hybrida, related, and ancestral genomes

Introduction

The garden petunia, Petunia hybrida (Solanaceae) is a fertile, diploid, annual hybrid species (2n=14) originating from P. axillaris and P. inflata 200 years ago. To understand the recent evolution of the P. hybrida genome, we examined tandemly repeated or satellite sequences using bioinformatic and molecular cytogenetic analysis.

Methods

Raw reads from available genomic assemblies and survey sequences of P. axillaris N (PaxiN), P. inflata S6, (PinfS6), P. hybrida (PhybR27) and the here sequenced P. parodii S7 (PparS7) were used for graph and k-mer based cluster analysis of TAREAN and RepeatExplorer. Analysis of repeat specific monomer lengths and sequence heterogeneity of the major tandem repeat families with more than 0.01% genome proportion were complemented by fluorescent in situ hybridization (FISH) using consensus sequences as probes to chromosomes of all four species.

Results

Seven repeat families, PSAT1, PSAT3, PSAT4, PSAT5 PSAT6, PSAT7 and PSAT8, shared high consensus sequence similarity and organisation between the four genomes. Additionally, many degenerate copies were present. FISH in P. hybrida and in the three wild petunias confirmed the bioinformatics data and gave corresponding signals on all or some chromosomes. PSAT1 is located at the ends of all chromosomes except the 45S rDNA bearing short arms of chromosomes II and III, and we classify it as a telomere associated sequence (TAS). It is the most abundant satellite repeat with over 300,000 copies, 0.2% of the genomes. PSAT3 and the variant PSAT7 are located adjacent to the centromere or mid-arm of one to three chromosome pairs. PSAT5 has a strong signal at the end of the short arm of chromosome III in P. axillaris and P.inflata, while in P. hybrida additional interstitial sites were present. PSAT6 is located at the centromeres of chromosomes II and III. PSAT4 and PSAT8 were found with only short arrays.

Discussion

These results demonstrate that (i) repeat families occupy distinct niches within chromosomes, (ii) they differ in the copy number, cluster organization and homogenization events, and that (iii) the recent genome hybridization in breeding P. hybrida preserved the chromosomal position of repeats but affected the copy number of repetitive DNA.

Living on the edge: morphological, karyological and genetic diversity studies of the Hungarian Plantago maxima populations and established ex situ collection

BACKGROUND

The analysis of genetic diversity of protected plant species can greatly support conservation efforts. Plantago maxima Juss. ex Jacq. is a perennial species distributed along the Eurasian steppe. The westernmost range edge of the species' distribution is located in the Pannonian basin, in Hungary where it is represented by a few, fragmented and highly endangered populations. We studied population diversity of all Hungarian range edge, natural populations, and one established ex situ population. One population from the centre of distribution (Kazakhstan) was implemented in the cpDNA haplotype study to compare the peripheral vs. central populations. We performed morphometric trait-based analysis, chromosome studies (morphometric analyses and FISH) and genetic diversity evaluations using inter simple sequence repeats (ISSR) and cpDNA trnL-trnF to evaluate differences between the in situ and ex situ populations as well as central vs. peripheral populations.

RESULTS

Our results showed no obvious morphological differences among the in situ and ex situ populations in the period between 2018 and 2020. One ex situ subpopulation develops flowers three years in a row from 2019, which is a favourable indicator of the introduction success. Hungarian populations are exclusively diploids (2n = 2x = 12). The karyogram consists of 5 metacentric and 1 acrocentric chromosome pair. Plantago maxima has one 35S and two 5S rDNA loci, located on the acrocentric chromosome pair. Eight variable ISSR primers yielded 100 fragments, of which 74.6% were polymorphic (mean H_e = 0.220). A high level of genetic variation within population was observed (92%) while the genetic differentiation among the populations was only 8%. STRUCTURE analysis revealed that the largest Kunpeszér population separated from the rest of the Hungarian populations, indicating a high rate of admixture among the other ones. Based on the trnL-trnF sequence analysis the Hungarian populations represent a single haplotype, which can indicate a reduced diversity due to isolation and recent population decline. By contrast, Kazakh population represents a distinct haplotype compared to the Hungarian samples.

CONCLUSIONS

The present study draws the attention to the high conservation value of the Plantago maxima populations from the westernmost range edge of the species' distribution.

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.

Development of Microsatellite Markers for Tanacetum cinerariifolium (Trevis.) Sch. Bip., a Plant with a Large and Highly Repetitive Genome

Dalmatian pyrethrum (Tanacetum cinerariifolium (Trevis.) Sch. Bip.) is an outcrossing plant species (2n = 18) endemic to the eastern Adriatic coast and source of the natural insecticide pyrethrin. Due to the high repeatability and large genome (1C-value = 9.58 pg) our previous attempts to develop microsatellite markers using the traditional method were unsuccessful. Now we have used Illumina paired-end whole genome sequencing and developed a specific procedure to obtain useful microsatellite markers. A total of 796,130,142 high-quality reads (approx. 12.5× coverage) were assembled into 6,909,675 contigs using two approaches (de novo assembly and joining of overlapped pair-end reads). A total of 31,380 contigs contained one or more microsatellite sequences, of which di-(59.7%) and trinucleotide (25.9%) repeats were the most abundant. Contigs containing microsatellites were filtered according to various criteria to achieve better yield of functional markers. After two rounds of testing, 17 microsatellite markers were developed and characterized in one natural population. Twelve loci were selected for preliminary genetic diversity analysis of three natural populations. Neighbor-joining tree, based on the proportion of shared alleles distances, grouped individuals into clusters according to population affiliation. The availability of codominant SSR markers will allow analysis of genetic diversity and structure of natural Dalmatian pyrethrum populations as well as identification of breeding lines and cultivars.

Interstitial Arabidopsis-Type Telomeric Repeats in Asteraceae

Tandem repeats of telomeric-like motifs at intra-chromosomal regions, known as interstitial telomeric repeats (ITR), have drawn attention as potential markers of structural changes, which might convey information about evolutionary relationships if preserved through time. Building on our previous work that reported outstanding ITR polymorphisms in the genus Anacyclus, we undertook a survey across 132 Asteraceae species, focusing on the six most speciose subfamilies and considering all the ITR data published to date. The goal was to assess whether the presence, site number, and chromosomal location of ITRs convey any phylogenetic signal. We conducted fluorescent in situ hybridization (FISH) using an Arabidopsis-type telomeric sequence as a probe on karyotypes obtained from mitotic chromosomes. FISH signals of ITR sites were detected in species of subfamilies Asteroideae, Carduoideae, Cichorioideae, Gymnarhenoideae, and Mutisioideae, but not in Barnadesioideae. Although six small subfamilies have not yet been sampled, altogether, our results suggest that the dynamics of ITR formation in Asteraceae cannot accurately trace the complex karyological evolution that occurred since the early diversification of this family. Thus, ITRs do not convey a reliable signal at deep or shallow phylogenetic levels and cannot help to delimitate taxonomic categories, a conclusion that might also hold for other important families such as Fabaceae.

Interstitial Telomeric-like Repeats (ITR) in Seed Plants as Assessed by Molecular Cytogenetic Techniques: A Review

The discovery of telomeric repeats in interstitial regions of plant chromosomes (ITRs) through molecular cytogenetic techniques was achieved several decades ago. However, the information is scattered and has not been critically evaluated from an evolutionary perspective. Based on the analysis of currently available data, it is shown that ITRs are widespread in major evolutionary lineages sampled. However, their presence has been detected in only 45.6% of the analysed families, 26.7% of the sampled genera, and in 23.8% of the studied species. The number of ITR sites greatly varies among congeneric species and higher taxonomic units, and range from one to 72 signals. ITR signals mostly occurs as homozygous loci in most species, however, odd numbers of ITR sites reflecting a hemizygous state have been reported in both gymnosperm and angiosperm groups. Overall, the presence of ITRs appears to be poor predictors of phylogenetic and taxonomic relatedness at most hierarchical levels. The presence of ITRs and the number of sites are not significantly associated to the number of chromosomes. The longitudinal distribution of ITR sites along the chromosome arms indicates that more than half of the ITR presences are between proximal and terminal locations (49.5%), followed by proximal (29.0%) and centromeric (21.5%) arm regions. Intraspecific variation concerning ITR site number, chromosomal locations, and the differential presence on homologous chromosome pairs has been reported in unrelated groups, even at the population level. This hypervariability and dynamism may have likely been overlooked in many lineages due to the very low sample sizes often used in cytogenetic studies.

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.

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.

Characterization and Dynamics of Repeatomes in Closely Related Species of Hieracium (Asteraceae) and Their Synthetic and Apomictic Hybrids

The repetitive content of the plant genome (repeatome) often represents its largest fraction and is frequently correlated with its size. Transposable elements (TEs), the main component of the repeatome, are an important driver in the genome diversification due to their fast-evolving nature. Hybridization and polyploidization events are hypothesized to induce massive bursts of TEs resulting, among other effects, in an increase of copy number and genome size. Little is known about the repeatome dynamics following hybridization and polyploidization in plants that reproduce by apomixis (asexual reproduction via seeds). To address this, we analyzed the repeatomes of two diploid parental species, Hieracium intybaceum and H. prenanthoides (sexual), their diploid F1 synthetic and their natural triploid hybrids (H. pallidiflorum and H. picroides, apomictic). Using low-coverage next-generation sequencing (NGS) and a graph-based clustering approach, we detected high overall similarity across all major repeatome categories between the parental species, despite their large phylogenetic distance. Medium and highly abundant repetitive elements comprise ∼70% of Hieracium genomes; most prevalent were Ty3/Gypsy chromovirus Tekay and Ty1/Copia Maximus-SIRE elements. No TE bursts were detected, neither in synthetic nor in natural hybrids, as TE abundance generally followed theoretical expectations based on parental genome dosage. Slight over- and under-representation of TE cluster abundances reflected individual differences in genome size. However, in comparative analyses, apomicts displayed an overabundance of pararetrovirus clusters not observed in synthetic hybrids. Substantial deviations were detected in rDNAs and satellite repeats, but these patterns were sample specific. rDNA and satellite repeats (three of them were newly developed as cytogenetic markers) were localized on chromosomes by fluorescence in situ hybridization (FISH). In a few cases, low-abundant repeats (5S rDNA and certain satellites) showed some discrepancy between NGS data and FISH results, which is due partly to the bias of low-coverage sequencing and partly to low amounts of the satellite repeats or their sequence divergence. Overall, satellite DNA (including rDNA) was markedly affected by hybridization, but independent of the ploidy or reproductive mode of the progeny, whereas bursts of TEs did not play an important role in the evolutionary history of Hieracium.

Evolutionary Dynamics of the Repetitive DNA in the Karyotypes of Pipa carvalhoi and Xenopus tropicalis (Anura, Pipidae)

The large amphibian genomes contain numerous repetitive DNA components that have played an important role in the karyotypic diversification of this vertebrate group. Hypotheses based on the presumable primitive karyotype (2n = 20) of the anurans of the family Pipidae suggest that they have evolved principally through intrachromosomal rearrangements. Pipa is the only South American pipid, while all the other genera are found in Africa. The divergence of the South American lineages from the African ones occurred at least 136 million years ago and is thought to have had a strong biogeographic component. Here, we tested the potential of the repetitive DNA to enable a better understanding of the differentiation of the karyotype among the family Pipidae and to expand our capacity to interpret the chromosomal evolution in this frog family. Our results indicate a long history of conservation in the chromosome bearing the H3 histone locus, corroborating inferences on the chromosomal homologies between the species in pairs 6, 8, and 9. The chromosomal distribution of the microsatellite motifs also provides useful markers for comparative genomics at the chromosome level between Pipa carvalhoi and Xenopus tropicalis, contributing new insights into the evolution of the karyotypes of these species. We detected similar patterns in the distribution and abundance of the microsatellite arrangements, which reflect the shared organization in the terminal/subterminal region of the chromosomes between these two species. By contrast, the microsatellite probes detected a differential arrangement of the repetitive DNA among the chromosomes of the two species, allowing longitudinal differentiation of pairs that are identical in size and morphology, such as pairs 1, 2, 4, and 5. We also found evidence of the distinctive composition of the repetitive motifs of the centromeric region between the species analyzed in the present study, with a clear enrichment of the (CA) and (GA) microsatellite motifs in P. carvalhoi. Finally, microsatellite enrichment in the pericentromeric region of chromosome pairs 6, 8, and 9 in the P. carvalhoi karyotype, together with interstitial telomeric sequences (ITS), validate the hypothesis that pericentromeric inversions occurred during the chromosomal evolution of P. carvalhoi and reinforce the role of the repetitive DNA in the remodeling of the karyotype architecture of the Pipidae.

Identification of tandem repeat families from long-read sequences of Humulus lupulus

Hop (Humulus lupulus L.) is known for its use as a bittering agent in beer and has a rich history of cultivation, beginning in Europe and now spanning the globe. There are five wild varieties worldwide, which may have been introgressed with cultivated varieties. As a dioecious species, its obligate outcrossing, non-Mendelian inheritance, and genomic structural variability have confounded directed breeding efforts. Consequently, understanding the hop genome represents a considerable challenge, requiring additional resources. In order to facilitate investigations into the transmission genetics of hop, we report here a tandem repeat discovery pipeline developed using k-mer filtering and dot plot analysis of PacBio long-read sequences from the hop cultivar Apollo. From this we identified 17 new and distinct tandem repeat sequence families, which represent candidates for FISH probe development. For two of these candidates, HuluTR120 and HuluTR225, we produced oligonucleotide FISH probes from conserved regions of and demonstrated their utility by staining meiotic chromosomes from wild hop, var. neomexicanus to address, for example, questions about hop transmission genetics. Collectively, these tandem repeat sequence families represent new resources suitable for development of additional cytogenomic tools for hop research.

The Utility of Graph Clustering of 5S Ribosomal DNA Homoeologs in Plant Allopolyploids, Homoploid Hybrids, and Cryptic Introgressants

INTRODUCTION

Ribosomal DNA (rDNA) loci have been widely used for identification of allopolyploids and hybrids, although few of these studies employed high-throughput sequencing data. Here we use graph clustering implemented in the RepeatExplorer (RE) pipeline to analyze homoeologous 5S rDNA arrays at the genomic level searching for hybridogenic origin of species. Data were obtained from more than 80 plant species, including several well-defined allopolyploids and homoploid hybrids of different evolutionary ages and from widely dispersed taxonomic groups.

RESULTS

(i) Diploids show simple circular-shaped graphs of their 5S rDNA clusters. In contrast, most allopolyploids and other interspecific hybrids exhibit more complex graphs composed of two or more interconnected loops representing intergenic spacers (IGS). (ii) There was a relationship between graph complexity and locus numbers. (iii) The sequences and lengths of the 5S rDNA units reconstituted in silico from k-mers were congruent with those experimentally determined. (iv) Three-genomic comparative cluster analysis of reads from allopolyploids and progenitor diploids allowed identification of homoeologous 5S rRNA gene families even in relatively ancient (c. 1 Myr) Gossypium and Brachypodium allopolyploids which already exhibit uniparental partial loss of rDNA repeats. (v) Finally, species harboring introgressed genomes exhibit exceptionally complex graph structures.

CONCLUSION

We found that the cluster graph shapes and graph parameters (k-mer coverage scores and connected component index) well-reflect the organization and intragenomic homogeneity of 5S rDNA repeats. We propose that the analysis of 5S rDNA cluster graphs computed by the RE pipeline together with the cytogenetic analysis might be a reliable approach for the determination of the hybrid or allopolyploid plant species parentage and may also be useful for detecting historical introgression events.