DNA Modification Patterns within the Transposable Elements of the Fig (Ficus carica L.) Genome

Transposable element activity can be harmful to the host’s genome integrity, but it can also provide selective advantages. One strategy to cope with transposons is epigenetic control through DNA base modifications. We report the non-canonic DNA modification dynamics of fig (Ficus carica L.) by exploiting high-quality genome reference and related N4-methylcytosine (4mC) and N6-methyladenine (6mA) data. Overall, 1.49% of transposon nucleotides showed either 4mC or 6mA modifications: the 4mC/6mA ratio was similar in Class I and Class II transposons, with a prevalence of 4mC, which is comparable to coding genes. Different percentages of 4mC or 6mA were observed among LTR-retrotransposon lineages and sub-lineages. Furthermore, both the Copia and Gypsy retroelements showed higher modification rates in the LTR and coding regions compared with their neighbour regions. Finally, the unconventional methylation of retrotransposons is unrelated to the number of close genes, suggesting that the 4mC and 6mA frequency in LTR-retrotransposons should not be related to transcriptional repression in the adjacency of the element. In conclusion, this study highlighted unconventional DNA modification patterns in fig transposable elements. Further investigations will focus on functional implications, in regards to how modified retroelements affect the expression of neighbouring genes, and whether these epigenetic markers can spread from repeats to genes, shaping the plant phenotype.