The application of LTR retrotransposons as molecular markers in plants

InpactorDB: A Plant LTR Retrotransposon Reference Library

LTR retrotransposons (LTR-RT) are major components of plant genomes. These transposable elements participate in the structure and evolution of genes and genomes through their mobility and their copy number amplification. For example, they are commonly used as evolutionary markers in genetic, genomic, and cytogenetic approaches. However, the plant research community is faced with the near absence of free availability of full-length, curated, and lineage-level classified LTR retrotransposon reference sequences. In this chapter, we will introduce InpactorDB, an LTR retrotransposon sequence database of 181 plant species representing 98 plant families for a total of 67,241 non-redundant elements. We will introduce how to use newly sequenced genomes to identify and classify LTR-RTs in a similar way with a standardized procedure using the Inpactor tool. InpactorDB is freely available at https://inpactordb.github.io .

Transposable Elements in the Revealing of Polymorphism-Based Differences in the Seeds of Flax Varieties Grown in Remediated Chernobyl Area

The nuclear reactor accident in Chernobyl, Ukraine, resulted in effects both locally and farther away. Most of the contaminated areas were the agricultural fields and forests. Experimental fields were established near Chernobyl—radioactively contaminated fields localized 5 km from Chernobyl Nuclear Power Plant as well as the remediated soil that is localized directly in the Chernobyl town. Two flax varieties growing under chronic exposition to ionizing radiation were used for this study—the local Ukrainian variety Kyivskyi and a commercial variety Bethune. The screening of the length polymorphism generated by transposable elements insertions were performed. All known types of common flax transposon, retrotransposons and iPBS approach were used. In the iPBS multiplex analyze, for the Kyivskyi variety, a unique addition was found in the seeds from the radioactive contaminated field and for the Bethune variety, a total of five amplicon additions were obtained and one deletion. For the TRIM Cassandra fingerprints, two amplicon additions were generated in the seeds from radioactive contaminated fields for the Bethune variety. In summary, the obtained data represent the genetic diversity between control and irradiated subgroups of flax seeds from Chernobyl area and the presence of activated transposable elements due to the irradiation stress.

Neonicotinoid insecticides exert diverse cytotoxic and genotoxic effects on cultivated sunflower

Contamination with neonicotinoids is a global problem affecting environment and target and non-target organisms including plants. The present study explored the potential genotoxic and cytotoxic effects of the insecticides Actara 25 WD and Nuprid 200 SL containing the active substances thiamethoxam (TMX) and imidacloprid (IMI), respectively, on cultivated sunflower (Helianthus annuus L.). The half maximal effective concentration (½EC₅₀) of the tested substances was calculated using a dose-response inhibition analysis of the growth of plant roots relative to the corresponding controls. Application of approximately ½EC₅₀ or higher TMX doses significantly increased the antioxidant activity in sunflower leaves, whereas IMI led to a significant decrease in root antioxidant capacity, indicating organ-specific insecticide effects on sunflower plants. Even low doses (½EC₅₀) of the studied neonicotinoids led to irregularities in mitotic phases and abnormalities in the cytokinesis and chromosome segregation, such as bridges, laggards, stickiness, and C-mitosis. Genotoxic effects manifested by a dose-independent induction of primary DNA damages and retrotransposon dynamics were also observed. The used set of physiological, biochemical, and genetic traits provides new information about the organ-specific effects of neonicotinoids in sunflower plants and elaborates on the complexity of mechanisms underpinning these effects that include DNA damages, cytokinesis defects, and genome instability.

Basic concepts and methodologies of DNA marker systems in plant molecular breeding

The concepts, methodologies and applications of some of the major molecular or DNA markers commonly used in plant science have been presented. The general principles of molecular marker techniques have been elucidated with detailed explanation of some notable basic concepts associated with marker applications: marker polymorphism, dominant or co-dominant mode of inheritance, agronomic trait-marker linkage, genetic mutations and variation. The molecular marker methods that have been extensively reviewed are RFLP, RAPD, SCAR, AFLP, SSR, CpSSR, ISSR, RAMP, SAMPL, SRAP, SSCP, CAPS, SNP, DArT, EST, and STS. In addition, the practicality of the retrotransposon-based marker methods, IRAP, REMAP, RBIP, and IPBS, have been discussed. Moreover, some salient characteristics of DNA markers have been compared and the various marker systems classified as PCR- or non-PCR-based, dominantly or co-dominantly inherited, locus specific or non-specific as well as at the levels of marker polymorphism and efficiency of marker reproducibility. Furthermore, the principles and methods of the following DNA markers have been highlighted: Penta-primer amplification refractory mutation system (PARMS), Conserved DNA-Derived Polymorphism (CDDP), P450-based analogue (PBA) markers, Tubulin-Based Polymorphism (TBP), Inter-SINE amplified polymorphism (ISAP), Sequence specific amplified polymorphism (S-SAP), Intron length polymorphisms (ILPs), Inter small RNA polymorphism (iSNAP), Direct amplification of length polymorphisms (DALP), Promoter anchored amplified polymorphism (PAAP), Target region amplification polymorphism (TRAP), Conserved region amplification polymorphism (CoRAP), Start Codon Targeted (SCoT) Polymorphism, and Directed Amplification of Minisatellite DNA (DAMD). Some molecular marker applications that have been recently employed to achieve various objectives in plant research have also been outlined. This review will serve as a useful reference resource for plant breeders and other scientists, as well as technicians and students who require basic know-how in the use of molecular or DNA marker technologies.

Development of molecular markers based on LTR retrotransposon in the Cleistogenes songorica genome

Long terminal repeat retrotransposons (LTR-RTs) contribute a large fraction of many sequenced plant genomes and play important roles in genomic diversity and phenotypic variations. LTR-RTs are abundantly distributed in plant genomes, facilitating the development of markers based on LTR-RTs for a variety of genotyping purposes. Whole-genome analysis of LTR-RTs was performed in Cleistogenes songorica. A total of 299,079 LTR-RTs were identified and classified as Gypsy type, Copia type, or other type. LTR-RTs were widely distributed in the genome, enriched in the heterochromatic region of the chromosome, and negatively correlated with gene distribution. However, approximately one-fifth of genes were still interrupted by LTR-RTs, and these genes are annotated. Furthermore, four types of primer pairs (PPs) were designed, namely, retrotransposon-based insertion polymorphisms, inter-retrotransposon amplified polymorphisms, insertion site-based polymorphisms, and retrotransposon-microsatellite amplified polymorphisms. A total of 350 PPs were screened in 23 accessions of the genus Cleistogenes, of which 80 PPs showed polymorphism, and 72 PPs showed transferability among Gramineae and non-Gramineae species. In addition, a comparative analysis of homologous LTR-RTs was performed with other related grasses. Taken together, the study will serve as a valuable resource for genotyping applications for C. songorica and related grasses.

Sequence-Specific Amplified Polymorphism (SSAP) and Sequence Characterized Amplified Region (SCAR) Markers in Zea mays

Transposable elements (TEs) are mobile, recurring DNA sequences scattered throughout genome and have a large impact on genome structure and function. Several genetic marker techniques were developed to exploit their ubiquitous nature. Sequence-specific amplified polymorphism (SSAP) is a TE-based genetic marker system that has been used in various purposes such as measuring genetic relatedness between species, deciphering the population structures, molecular tagging for agronomic development in marker-assisted breeding (MAS). In addition to SSAP, sequence characterized amplified region (SCAR) from the SSAP markers provides an added advantage in identifying qualitative traits. Once developed SCAR markers are efficient, fast, and reliable method for genetic evaluations. These methods can be useful especially for the crops which have no genetic sequence information. With improved discriminatory ability they offer access to dynamic and polymorphic regions of genome. These techniques can be useful in breeding programs to improve or develop high yielding crops.

Transposon-Based Tagging In Silico Using FastPCR Software

Retrotransposons are ubiquitous, generally dispersed components of eukaryotic genomes. These properties, together with their "copy and paste" lifecycle that generates insertional polymorphism without need for excision, makes them widely useful as a molecular-genetic tags. Various tagging systems have been developed that exploit the sequence conservation of retrotransposon components, such as those found in their long terminal repeats (LTRs). To detect polymorphisms for retrotransposon insertions, marker systems generally rely on PCR amplification between the termini and some component of flanking genomic DNA. As complements to various "wet lab" protocols for retrotransposon tagging, in silico bioinformatics approaches are useful for predicting likely outcomes from unsequenced accessions on the basis of reference genomes. In this chapter, we describe protocols for in silico retrotransposon-based fingerprinting techniques using the FastPCR software as an integrated tools environment for in silico PCR primer design and analysis.

Genome-Wide Identification and Development of LTR Retrotransposon-Based Molecular Markers for the Melilotus Genus

Melilotus is an important genus of legumes with industrial and medicinal value, partly due to the production of coumarin. To explore the genetic diversity and population structure of Melilotus, 40 accessions were analyzed using long terminal repeat (LTR) retrotransposon-based markers. A total of 585,894,349 bp of LTR retrotransposon sequences, accounting for 55.28% of the Melilotus genome, were identified using bioinformatics tools. A total of 181,040 LTR retrotransposons were identified and classified as Gypsy, Copia, or another type. A total of 350 pairs of primers were designed for assessing polymorphisms in 15 Melilotus albus accessions. Overall, 47 polymorphic primer pairs were screened for their availability and transferability in 18 Melilotus species. All the primer pairs were transferable, and 292 alleles were detected at 47 LTR retrotransposon loci. The average polymorphism information content (PIC) value was 0.66, which indicated that these markers were highly informative. Based on unweighted pair group method with arithmetic mean (UPGMA) dendrogram cluster analysis, the 18 Melilotus species were classified into three clusters. This study provides important data for future breeding programs and for implementing genetic improvements in the Melilotus genus.

Development and Deployment of High-Throughput Retrotransposon-Based Markers Reveal Genetic Diversity and Population Structure of Asian Bamboo

Bamboo, a non-timber grass species, known for exceptionally fast growth is a commercially viable crop. Long terminal repeat (LTR) retrotransposons, the main class I mobile genetic elements in plant genomes, are highly abundant (46%) in bamboo, contributing to genome diversity. They play significant roles in the regulation of gene expression, chromosome size and structure as well as in genome integrity. Due to their random insertion behavior, interspaces of retrotransposons can vary significantly among bamboo genotypes. Capitalizing this feature, inter-retrotransposon amplified polymorphism (IRAP) is a high-throughput marker system to study the genetic diversity of plant species. To date, there are no transposon based markers reported from the bamboo genome and particularly using IRAP markers on genetic diversity. Phyllostachys genus of Asian bamboo is the largest of the Bambusoideae subfamily, with great economic importance. We report structure-based analysis of bamboo genome for the LTR-retrotransposon superfamilies, Ty3-gypsy and Ty1-copia, which revealed a total of 98,850 retrotransposons with intact LTR sequences at both the ends. Grouped into 64,281 clusters/scaffold using CD-HIT-EST software, only 13 clusters of retroelements were found with more than 30 LTR sequences and with at least one copy having all intact protein domains such as gag and polyprotein. A total of 16 IRAP primers were synthesized, based on the high copy numbers of conserved LTR sequences. A study using these IRAP markers on genetic diversity and population structure of 58 Asian bamboo accessions belonging to the genus Phyllostachys revealed 3340 amplicons with an average of 98% polymorphism. The bamboo accessions were collected from nine different provinces of China, as well as from Italy and America. A three phased approach using hierarchical clustering, principal components and a model based population structure divided the bamboo accessions into four sub-populations, PhSP1, PhSP2, PhSP3 and PhSP4. All the three analyses produced significant sub-population wise consensus. Further, all the sub-populations revealed admixture of alleles. The analysis of molecular variance (AMOVA) among the sub-populations revealed high intra-population genetic variation (75%) than inter-population. The results suggest that Phyllostachys bamboos are not well evolutionarily diversified, although geographic speciation could have occurred at a limited level. This study highlights the usability of IRAP markers in determining the inter-species variability of Asian bamboos.

Retrotransposon-based genetic variation of Poa annua populations from contrasting climate conditions

Background

Poa annua L. is an example of a plant characterized by abundant, worldwide distribution from polar to equatorial regions. Due to its high plasticity and extraordinary expansiveness, P. annua is considered an invasive species capable of occupying and surviving in a wide range of habitats including pioneer zones, areas intensively transformed by human activities, remote subarctic meadows and even the Antarctic Peninsula region.

Methods

In the present study, we evaluated the utility of inter-primer binding site (iPBS) markers for assessing the genetic variation of P. annua populations representing contrasting environments from the worldwide range of this species. The electrophoretic patterns of polymerase chain reaction products obtained for each individual were used to estimate the genetic diversity and differentiation between populations.

Results

iPBS genotyping revealed a pattern of genetic variation differentiating the six studied P. annua populations characterized by their different climatic conditions. According to the analysis of molecular variance, the greatest genetic variation was recorded among populations, whereas 41.75% was observed between individuals within populations. The results of principal coordinates analysis (PCoA) and model-based clustering analysis showed a clear subdivision of analyzed populations. According to PCoA, populations from Siberia and the Kola Peninsula were the most different from each other and showed the lowest genetic variability. The application of STRUCTURE software confirmed the unique character of the population from the Kola Peninsula.

Discussion

The lowest variability of the Siberia population suggested that it was subjected to genetic drift. However, although demographic expansion was indicated by negative values of Fu’s F_S statistic and analysis of mismatch distribution, it was not followed by significant traces of a bottleneck or a founder effect. For the Antarctic population, the observed level of genetic variation was surprisingly high, despite the observed significant traces of bottleneck/founder effect following demographic expansion, and was similar to that observed in populations from Poland and the Balkans. For the Antarctic population, the multiple introduction events from different sources are considered to be responsible for such an observation. Moreover, the results of STRUCTURE and PCoA showed that the P. annua from Antarctica has the highest genetic similarity to populations from Europe.

Conclusions

The observed polymorphism should be considered as a consequence of the joint influence of external abiotic stress and the selection process. Environmental changes, due to their ability to induce transposon activation, lead to the acceleration of evolutionary processes through the production of genetic variability.

The landscape and structural diversity of LTR retrotransposons in Musa genome

Long terminal repeat retrotransposons represent a major component of plant genomes and act as drivers of genome evolution and diversity. Musa is an important fruit crop and also used as a starchy vegetable in many countries. BAC sequence analysis by dot plot was employed to investigate the LTR retrotransposons from Musa genomes. Fifty intact LTR retrotransposons from selected Musa BACs were identified by dot plot analysis and further BLASTN searches retrieved 153 intact copies, 61 truncated, and a great number of partial copies/remnants from GenBank database. LARD-like elements were also identified with several copies dispersed among the Musa genotypes. The predominant elements were the LTR retrotransposons Copia and Gypsy, while Caulimoviridae (pararetrovirus) were rare in the Musa genome. PCR amplification of reverse transcriptase (RT) sequences revealed their abundance in almost all tested Musa accessions and their ancient nature before the divergence of Musa species. The phylogenetic analysis based on RT sequences of Musa and other retrotransposons clustered them into Gypsy, Caulimoviridae, and Copia lineages. Most of the Musa-related elements clustered in their respective groups, while some grouped with other elements indicating homologous sequences. The present work will be helpful to understand the LTR retrotransposons landscape, giving a complete picture of the nature of the elements, their structural features, annotation, and evolutionary dynamics in the Musa genome.

Assessment of genetic diversity in Nordic timothy (Phleum pratense L.)

BACKGROUND

Timothy (Phleum pratense L.), a cool-season hexaploid perennial, is the most important forage grass species in Nordic countries. Earlier analyses of genetic diversity in a collection of 96 genebank accessions of timothy with SSR markers demonstrated high levels of diversity but could not resolve population structure. Therefore, we examined a subset of 51 accessions with REMAP markers, which are based on retrotransposons, and compared the diversity results with those obtained with SSR markers.

RESULTS

Using four primer combinations, 533 REMAP markers were analyzed, compared with 464 polymorphic alleles in the 13 SSR loci previously. The average marker index, which describes information obtained per experiment (per primer combination or locus) was over six times higher with REMAPs. Most of the variation found was within accessions, with somewhat less, 89 %, for REMAPs, than for SSR, with 93 %.

CONCLUSIONS

SSRs revealed differences in the level of diversity slightly better than REMAPs but neither marker type could reveal any clear clustering of accessions based on countries, vegetation zones, or different cultivar types. In our study, reliable evaluation of SSR allele dosages was not possible, so each allele had to be handled as a dominant marker. SSR and REMAP, which report from different mechanisms of generating genetic diversity and from different genomic regions, together indicate a lack of population structure. Taken together, this likely reflects the outcrossing and hexaploid nature of timothy rather than failures of either marker system.

Exploiting the Repetitive Fraction of the Wheat Genome for High-Throughput Single-Nucleotide Polymorphism Discovery and Genotyping

Transposable elements (TEs) account for more than 80% of the wheat genome. Although they represent a major obstacle for genomic studies, TEs are also a source of polymorphism and consequently of molecular markers such as insertion site-based polymorphism (ISBP) markers. Insertion site-based polymorphisms have been found to be a great source of genome-specific single-nucleotide polymorphism (SNPs) in the hexaploid wheat ( L.) genome. Here, we report on the development of a high-throughput SNP discovery approach based on sequence capture of ISBP markers. By applying this approach to the reference sequence of chromosome 3B from hexaploid wheat, we designed 39,077 SNPs that are evenly distributed along the chromosome. We demonstrate that these SNPs can be efficiently scored with the KASPar (Kompetitive allele-specific polymerase chain reaction) genotyping technology. Finally, through genetic diversity and genome-wide association studies, we also demonstrate that ISBP-derived SNPs can be used in marker-assisted breeding programs.

Development of IRAP- and REMAP-derived SCAR markers for marker-assisted selection of the stripe rust resistance gene Yr15 derived from wild emmer wheat

Yr15 provides broad resistance to stripe rust, an important wheat disease. REMAP- and IRAP-derived co-dominant SCAR markers were developed and localize Yr15 to a 1.2 cM interval. They are reliable across many cultivars. Stripe rust [Pucinia striiformis f.sp. tritici (Pst)] is one of the most important fungal diseases of wheat, found on all continents and in over 60 countries. Wild emmer wheat (Triticum dicoccoides), which is the tetraploid progenitor of durum wheat, is a valuable source of novel stripe rust resistance genes for wheat breeding. T. dicoccoides accession G25 carries Yr15 on chromosome 1BS. Yr15 confers resistance to virtually all tested Pst isolates; it is effective in durum and bread wheat introgressions and their derivatives. Retrotransposons generate polymorphic insertions, which can be scored as Mendelian markers using techniques such as REMAP and IRAP. Six REMAP- and IRAP-derived SCAR markers were mapped using 1,256 F2 plants derived from crosses of the susceptible T. durum accession D447 (DW1) with its resistant BC3F9 and BC3F10 (B9 and B10) near isogenic lines, which carried Yr15 introgressed from G25. The nearest markers segregated 0.1 cM proximally and 1.1 cM distally to Yr15. These markers were also mapped and validated at the same position in another 500 independent F2 plants derived from crosses of B9 and B10 with the susceptible cultivar Langdon (LDN). SC2700 and SC790, defining Yr15 on an interval of 1.2 cM, were found to be reliable and robust co-dominant markers in a wide range of wheat lines and cultivars with and without Yr15. These markers are useful tags in marker-assisted wheat breeding programs that aim to incorporate Yr15 into elite wheat lines and cultivars for durable and broad-spectrum resistance to stripe rust.

A simple and efficient method to isolate LTR sequences of plant retrotransposon

Retrotransposons (RTNs) have important roles in the formation of plant genome size, structure, and evolution. Ubiquitous distributions, abundant copy numbers, high heterogeneities, and insertional polymorphisms of RTNs have made them as excellent sources for molecular markers development. However, the wide application of RTNs-based molecular markers is restricted by the scarcity of the LTR (long terminal repeat) sequences information. A new, simple, and efficient method to isolate LTR sequences of RTNs was presented based on the degenerate RNase H nested primers and PPT (polypurine tract) primer of RTNs in tree peony. This method combined the characteristics and advantages of high-efficiency thermal asymmetric interlaced PCR (hiTAIL-PCR), annealing control primer (ACP) system, and suppression PCR method. Nineteen LTR sequences were isolated using this new method in tree peony and the applicability of the LTR sequences based markers was validated by further SSAP analysis. The results showed that the new method is simple, of low-cost, and highly efficient, which is just conducted by three rounds of PCR and does not need any restriction enzymes and adapters, much less the hybridizations. This new method is rapid, economical, and cost- and time-saving, which could be easily used to isolate LTR sequences of RTNs.

Reverse transcriptase domain sequences from tree peony (Paeonia suffruticosa) long terminal repeat retrotransposons: sequence characterization and phylogenetic analysis

Tree peony is an important horticultural plant worldwide of great ornamental and medicinal value. Long terminal repeat retrotransposons (LTR-retrotransposons) are the major components of most plant genomes and can substantially impact the genome in many ways. It is therefore crucial to understand their sequence characteristics, genetic distribution and transcriptional activity; however, no information about them is available in tree peony. Ty1-copia-like reverse transcriptase sequences were amplified from tree peony genomic DNA by polymerase chain reaction (PCR) with degenerate oligonucleotide primers corresponding to highly conserved domains of the Ty1-copia-like retrotransposons in this study. PCR fragments of roughly 270 bp were isolated and cloned, and 33 sequences were obtained. According to alignment and phylogenetic analysis, all sequences were divided into six families. The observed difference in the degree of nucleotide sequence similarity is an indication for high level of sequence heterogeneity among these clones. Most of these sequences have a frame shift, a stop codon, or both. Dot-blot analysis revealed distribution of these sequences in all the studied tree peony species. However, different hybridization signals were detected among them, which is in agreement with previous systematics studies. Reverse transcriptase PCR (RT-PCR) indicated that Ty1-copia retrotransposons in tree peony were transcriptionally inactive. The results provide basic genetic and evolutionary information of tree peony genome, and will provide valuable information for the further utilization of retrotransposons in tree peony.

Genetic diversity of maize germplasm assessed by retrotransposon-based markers

Maize is one of the most important crops and also a model for grass genome research. Transposable elements comprise over 78% of the maize genome and their ability to generate new copies makes them good potential markers. Interretrotransposon-amplified polymorphism (IRAP) and retrotransposon microsatellite amplified polymorphism (REMAP) protocols were used for the first time in maize to study the genetic variability between maize cultivars. Ten PCR primers were selected based on a systematic analysis of the sequence conservation in the extremities of different high copy number transposable elements, whereas one primer was chosen based on a microsatellite sequence. Of the 16 primer combinations tested, 14 produced polymorphic bands. These markers were used to identify genetic similarity among 20 maize cultivars selected by their different kernel oil content. Genetic similarity analysis was performed based on the polymorphic band profiles and dendrograms were developed by the unweighted pair-group method with arithmetic averages. Clustering technique revealed that samples were grouped into three clusters that differed in their kernel oil content and size, and in their relative embryo size. In the current investigation, there is evidence that IRAP/REMAP may be useful as markers in maize.

Transposon fingerprinting using low coverage whole genome shotgun sequencing in cacao (Theobroma cacao L.) and related species

Background

Transposable elements (TEs) and other repetitive elements are a large and dynamically evolving part of eukaryotic genomes, especially in plants where they can account for a significant proportion of genome size. Their dynamic nature gives them the potential for use in identifying and characterizing crop germplasm. However, their repetitive nature makes them challenging to study using conventional methods of molecular biology. Next generation sequencing and new computational tools have greatly facilitated the investigation of TE variation within species and among closely related species.

Results

(i) We generated low-coverage Illumina whole genome shotgun sequencing reads for multiple individuals of cacao (Theobroma cacao) and related species. These reads were analysed using both an alignment/mapping approach and a de novo (graph based clustering) approach. (ii) A standard set of ultra-conserved orthologous sequences (UCOS) standardized TE data between samples and provided phylogenetic information on the relatedness of samples. (iii) The mapping approach proved highly effective within the reference species but underestimated TE abundance in interspecific comparisons relative to the de novo methods. (iv) Individual T. cacao accessions have unique patterns of TE abundance indicating that the TE composition of the genome is evolving actively within this species. (v) LTR/Gypsy elements are the most abundant, comprising c.10% of the genome. (vi) Within T. cacao the retroelement families show an order of magnitude greater sequence variability than the DNA transposon families. (vii) Theobroma grandiflorum has a similar TE composition to T. cacao, but the related genus Herrania is rather different, with LTRs making up a lower proportion of the genome, perhaps because of a massive presence (c. 20%) of distinctive low complexity satellite-like repeats in this genome.

Conclusions

(i) Short read alignment/mapping to reference TE contigs provides a simple and effective method of investigating intraspecific differences in TE composition. It is not appropriate for comparing repetitive elements across the species boundaries, for which de novo methods are more appropriate. (ii) Individual T. cacao accessions have unique spectra of TE composition indicating active evolution of TE abundance within this species. TE patterns could potentially be used as a “fingerprint” to identify and characterize cacao accessions.