# Genome-wide expert annotation of the epigenetic machinery of the plant-parasitic nematodes Meloidogyne spp., with a focus on the asexually reproducing species

Loris Pratx, Corinne Rancurel, Martine da Rocha, Etienne Danchin, Philippe

Castagnone-Sereno, Pierre Abad, Laetitia Zurletto

► To cite this version:

Loris Pratx, Corinne Rancurel, Martine da Rocha, Etienne Danchin, Philippe Castagnone-Sereno, et al.. Genome-wide expert annotation of the epigenetic machinery of the plant-parasitic nematodes Meloidogyne spp., with a focus on the asexually reproducing species. BMC Genomics, 2018, 19, pp.1-21. 10.1186/s12864-018-4686-x. hal-02626482

HAL Id: hal-02626482

https://hal.inrae.fr/hal-02626482

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
RESEARCH ARTICLE

Open Access

<figure><img src="image_2.png" /></figure>

CrossMark

# Genome-wide expert annotation of the epigenetic machinery of the plant-parasitic nematodes *Meloidogyne* spp., with a focus on the asexually reproducing species

Loris Pratx<fnref n="1" />,<fnref n="2" />, Corinne Rancurel<fnref n="1" />,<fnref n="2" />, Martine Da Rocha<fnref n="1" />,<fnref n="2" />, Etienne G. J. Danchin<fnref n="1" />,<fnref n="2" />, Philippe Castagnone-Sereno<fnref n="1" />,<fnref n="2" />, Pierre Abad<fnref n="1" />,<fnref n="2" /> and Laetitia Perfus-Barbeoch<fnref n="1" />,<fnref n="2" />*

**Abstract**

**Background:** The renewed interest in epigenetics has led to the understanding that both the environment and individual lifestyle can directly interact with the epigenome to influence its dynamics. Epigenetic phenomena are mediated by DNA methylation, stable chromatin modifications and non-coding RNA-associated gene silencing involving specific proteins called epigenetic factors.

Multiple organisms, ranging from plants to yeast and mammals, have been used as model systems to study epigenetics. The interactions between parasites and their hosts are models of choice to study these mechanisms because the selective pressures are strong and the evolution is fast.

The asexually reproducing root-knot nematodes (RKN) offer different advantages to study the processes and mechanisms involved in epigenetic regulation. RKN genomes sequencing and annotation have identified numerous genes, however, which of those are involved in the adaption to an environment and potentially relevant to the evolution of plant-parasitism is yet to be discovered.

**Results:** Here, we used a functional comparative annotation strategy combining orthology data, mining of curated genomics as well as protein domain databases and phylogenetic reconstructions.

Overall, we show that (i) neither RKN, nor the model nematode *Caenorhabditis elegans* possess any DNA methyltransferases (DNMT) (ii) RKN do not possess the complete machinery for DNA methylation on the 6th position of adenine (6mA) (iii) histone (de)acetylation and (de)methylation pathways are conserved between *C. elegans* and RKN, and the corresponding genes are amplified in asexually reproducing RKN (iv) some specific non-coding RNA families found in plant-parasitic nematodes are dissimilar from those in *C. elegans*. In the asexually reproducing RKN *Meloidogyne incognita*, expression data from various developmental stages supported the putative role of these proteins in epigenetic regulations.

**Conclusions:** Our results refine previous predictions on the epigenetic machinery of model species and constitute the most comprehensive description of epigenetic factors relevant to the plant-parasitic lifestyle and/or asexual mode of reproduction of RKN. Providing an atlas of epigenetic factors in RKN is an informative resource that will enable researchers to explore their potential role in adaptation of these parasites to their environment.

**Keywords:** Adaptation, Epigenetics, Functional comparative annotation, Orthology, Chromatin, Histone modifications, Root-knot nematodes, Meloidogyne, Plant parasitism, Asexual reproduction

<figure><img src="image_3.png" /></figure>

BMC

© The Author(s). 2018 **Open Access** This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

<footnote>* Correspondence: laetitia zurletto@inra.fr</footnote>

<footnote n="1">¹Université Côte d'Azur, INRA, ISA, Sophia Antipolis, France</footnote>

<footnote n="2">²Institut Sophia Agrobiotech, 400, route des chappes, BP 167 - 06903 Sophia Antipolis Cedex, France</footnote>
## Background

Epigenetic modifications are heritable yet metastable and cannot be explained by changes in nucleotide sequence [1]. In eukaryotes, packaging of DNA into chromatin has profound effects on cellular processes that utilize DNA as template, including transcription, replication, recombination and repair [2]. The nucleosome is the fundamental unit of chromatin and it is composed of an octamer of histones around which 147 base pairs of DNA are wrapped [3]. Subsequent compaction leads to higher order structures including the formation of very dense arrays of nucleosomes observed in heterochromatin [4]. Despite being tightly packed, the chromatin appears to be highly plastic thanks to several factors that influence both its local and global architecture. Two mechanisms of epigenetic regulations are generally considered: methylation of DNA and post-translational modifications of histones [5, 6]. Besides these major processes, a growing body of evidence indicates that regulatory non-coding RNAs play an important role in epigenetic control [7]. Epigenetics research conducted so far has raised new evidence about how environmental factors can impact the mechanisms through which biological processes and functions are regulated (for review, [8, 9]). In that respect, a pivotal role of epigenetic mechanisms has been shown in controlling various strategies of pathogens to hijack host cell pathways [10, 11].

The root-knot nematode (RKN) *Meloidogyne incognita* is a plant parasite of major agricultural importance which reproduces exclusively by mitotic parthenogenesis [12]. Despite its clonal mode of reproduction, this nematode can adapt to unfavorable conditions. For instance, avirulent *M. incognita* strains controlled by a resistance gene are able to overcome plant resistance and become virulent on these plants [13]. Considering that the proportion of individuals with this phenotype rises gradually over generations and does not follow a Mendelian transmission [14, 15], it is possible that epigenetic changes could represent one of the main driving forces of evolution in this species [16].

Here, we study the epigenetic determinants possibly involved in the phenotypic plasticity of this parasite based on four sequenced RKN genomes [17–19], and information on core epigenetic proteins available in public databases [20]. So far, a structured source for such information is missing for RKN but also for the model nematode *C. elegans*.

We provide a manually curated annotation giving information about more than 3500 candidate epigenetic regulators in 20 species, including four RKN, five other nematodes and a set of model species of interest.

For *M. incognita* genes, we included expression data across several developmental stages. Such combination of functional annotation on RKN epigenetic factors is

relevant to cover all possible mechanisms that could underlie the success of RKN pathogens.

## Methods

### Proteomes used for functional comparative annotation

For comparative annotation of epigenetic factors of the four RKN species (*M. incognita*, *M. javanica*, *M. arenaria*, *M. hapla*), we selected a total of 16 species from publicly available protein databases, including species with specific lifestyles (Table 1). Six model species, *C. elegans*, *Drosophila melanogaster*, *Homo sapiens*, *Saccharomyces cerevisae*, *Schizosaccharomyces pombe* and *Arabidopsis thaliana*, were selected because of their genome completeness, annotation quality and the existence of experimental characterization of epigenetic factors. We included 10 other species based on their phylogenetic position and/or lifestyles: two other plant-parasitic nematodes (*Globodera pallida*, *G. rostochiensis*), two animal-parasitic nematodes (*Trichinella spiralis*, *Ascaris suum*), and two plant-parasitic fungi with previous report of epigenetic regulations (*Botrytis cinerea*, *Leptosphaeria maculans* [21, 22]). Finally, we added four species which possess an epigenome of interest, known to be involved in phenotypical plasticity (*Acyrthosiphon pisum*, *Apis mellifera* [23, 24]) or in pathogenicity (*Schistosoma mansoni*, *Plasmodium falciparum* [25–28]). The 20 proteomes were all downloaded from Uniprot database [29], with the exception of *Meloidogyne* spp. available at "Meloidogyne genomic resources" website(https://meloidogyne.inra.fr) [19] and *Globodera* spp. downloaded from WormBase ParaSite [30, 31].

### Identification of core epigenetic factors

Epigenetic factors were identified thanks to a custom pipeline including four annotation steps (Fig. 1). The first step consisted in the identification of both (i) orthology links between the 20 proteomes by using OrthoMCL version 2.0 under standard parameters [32], and (ii) search of specific Pfam [33] protein domains assigned by Interproscan [34]. To explore and analyze homology clusters we used a web server called Family-Companion [35].

As a second step, we created a repertory of known epigenetic factors and associated protein domains from the 6 model species. Basically, we recovered all known chromatin factors based on (i) Swiss-Prot [29] and Wormbase [36] annotations, (ii) epigenetic factors-specific databases such as Histone database [37], Chro-moHub [38] and Histome [39], and (iii) literature [40–46]. This dataset of epigenetic factors constituted the reference dataset that we used for comparative genome annotation of RKN (Additional file 1: Table S1). Then, from this reference dataset, we looked for protein domains specifically associated with known epigenetic
<table><caption>**Table 1** Characteristics of the 20 species used for functional comparative annotation of core epigenetic factors</caption><thead><tr><th rowspan="3">Species</th><th colspan="10">Characteristics</th><th colspan="3">DNA methylation</th><th colspan="6">Non-coding RNA</th></tr><tr><th rowspan="2">Nematode</th><th rowspan="2">Model species</th><th colspan="2">Parasite</th><th colspan="3">RKN peculiarities</th><th colspan="2">5mC</th><th colspan="2">6mA</th><th rowspan="2">Methyl transfersases DNMT</th><th rowspan="2">Methyltransferase DAMT-1</th><th rowspan="2">Demethylase NMAD-1</th><th rowspan="2">Histone PTM</th><th colspan="2">Argonautes</th><th colspan="4">Non-coding RNA gene regulation</th></tr><tr><th>Animal</th><th>Plant</th><th>Reproductive mode</th><th>Genome structure</th><th>GO</th><th>PIWI</th><th>WAGO</th><th>DICER complex</th><th>RdRp</th><th>Drosha</th></tr></thead><tbody><tr><td>Meloidogyne incognita</td><td>X</td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>•</td><td></td><td>•</td><td>X</td><td></td><td>X</td><td>•</td><td>•</td><td>•</td></tr><tr><td>Meloidogyne javanica</td><td>X</td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>•</td><td></td><td>•</td><td></td><td>•</td><td>•</td><td>•</td><td>•</td><td>•</td></tr><tr><td>Meloidogyne arenaria</td><td>X</td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>•</td><td></td><td>•</td><td></td><td>•</td><td>•</td><td>•</td><td>•</td><td>•</td></tr><tr><td>Meloidogyne hapla</td><td>X</td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>•</td><td></td><td>•</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>•</td></tr><tr><td>Globodera pallida</td><td>X</td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>•</td><td></td><td>•</td><td></td><td>•</td><td>•</td><td>X</td><td>X</td><td>•</td></tr><tr><td>Globodera rostochiensis</td><td>X</td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>•</td><td></td><td>•</td><td></td><td>•</td><td>•</td><td>•</td><td>•</td><td>•</td></tr><tr><td>Ascaris suum</td><td>X</td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>•</td><td></td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>•</td></tr><tr><td>Trichinella spiralis</td><td>X</td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>•</td><td>•</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>•</td></tr><tr><td>Caenorhabditis elegans</td><td>X</td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td></tr><tr><td>Homo sapiens</td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td></tr><tr><td>Drosophila melanogaster</td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td></tr><tr><td>Saccharomyces cerevisiae S288c</td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>•</td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>X</td><td>X</td></tr><tr><td>Schizosaccharomyces pombe 972h-</td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>•</td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>X</td><td>X</td></tr><tr><td>Arabidopsis thaliana</td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td></tr><tr><td>Acyrthosiphon pisum</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>•</td><td>•</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td></tr><tr><td>Schistosoma mansoni</td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>•</td><td>•</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td></tr><tr><td>Plasmodium falciparum 3D7</td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Apis mellifera ligustica</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>X</td><td>•</td><td>•</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>•</td></tr><tr><td>Leptosphaeria maculans</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>•</td><td>•</td><td>•</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>•</td></tr><tr><td>Botrytis cinerea T4</td><td></td><td></td><td>X</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>•</td><td></td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>X</td><td>•</td></tr></tbody></table>

A cross indicates either if the characteristic is true or if data about epigenetic mechanism is available for one species. Blue circles indicate new core epigenetic factors annotated in this study

factors that could be used for functional annotation. These features were provided by Pfam database of curated protein families [47–49]. We only kept Pfam domains that could uniquely be linked to epigenetic factors (thereafter called “Pfam epigenetic factor list”; Additional file 2: Table S2). For instance, PCAF (P300/ CBP-associated factor) N-terminal domain (PF06466) and acetyltransferase 1 domain (PF00583) were both associated with Gcn5-related N-acetyltransferases (GNAT) family of histone acetyltransferase. However, while PCAF N-terminal domain (PF06466) was restricted to histone acetyltransferases, we found acetyltransferase 1 domain (PF00583) also carried by proteins that do not have an activity on histone (for instance choline acetyltransferases). In that case, only PCAF N-terminal domain (PF06466) was kept for functional annotation.

The goal of the third step was to identify epigenetic factors among OrthoMCL groups and/or Interproscan annotation, from step one. Proteins orthologous to at least one protein of the reference dataset and/or containing at least one Pfam domain of interest were annotated as putative epigenetic factors. Any protein not positioned into any OrthoMCL group was called “singleton” as it possessed no evident ortholog or in-paralog.

To validate the accuracy of the annotation, and to determine the closest orthologs between *C. elegans* and RKN, we performed phylogenetic constructions as the fourth step. Putative epigenetic factors protein sequences were aligned by MAFFT version 7.245 [50] with -auto option. The alignment was manually checked to remove misaligned sequences. We used trimAl version 1.2 [51] to remove gap-rich columns in the multiple alignments.

Maximum likelihood trees were built with RAxML version 8.1 [52] with an automatic detection of the fittest evolutionary model, and an estimated gamma distribution of the rates of evolution (PROTGAMMAUTO option). Rapid bootstrap replicates followed by a full ML analysis were conducted. We used the -autoMRE criterion to stop bootstrap replicated upon convergence. For each phylogeny, the best scoring ML tree with associated bootstrap support values was retrieved as final resulting ML topology. The trees were visualized with FigTree version 1.4.2 [53].

## Expression analysis

To enhance robustness of the expert annotation, we integrated an experimental transcript verification step. We focused on *M. incognita* because RNA-seq data were available on six developmental stages of this species in our lab [54]. For each *M. incognita* putative epigenetic factor, reads per kilobase of transcript per million mapped reads (RPKM) was calculated. We considered that putative epigenetic factor annotations were supported by existence of transcripts when RPKM≥ 5.

## Results

Epigenetic factors were identified thank to a comparative annotation based on 6 model species. We selected 716 proteins, to build the reference dataset of epigenetic factors. Among these 716 proteins, 288 were from *H. sapiens*, 116 from *C. elegans*, 64 from *D. melanogaster*, 49 from *S. cerevisiae*, 44 from *S. pombe* and 155 from *A. thaliana* (Additional file 1: Table S1 and Additional file 2: Table S2). Within this reference dataset, 29 Pfam domains were
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Pipeline for annotation of core epigenetic factors with a focus on root-knot nematodes (RKN). Step 1 consisted of the annotation, based on both OrthoMCL and Interproscan, of 4 RKN species together with 6 model species and 10 other species selected for their lifestyles and/or epigenetic interest. In parallel, step 2 consisted of creating a repertory of known epigenetic factors from public data on 6 model species. Step 3 consisted of the extraction of putative epigenetic factor from step 1 data, based on comparison with reference dataset created during step 2. Finally, step 4 consisted of the validation and characterization of epigenetic factors thanks to phylogenetic analysis</figcaption></figure>

found to be specific to epigenetic factors (i.e. not found outside epigenetic factors) (Additional file 2: Table S2).

## The methodology is effective for the identification of epigenetic factors

To assess the efficiency of this annotation methodology, we tested the pipeline (Fig. 1) on two species, *S. mansoni* and *A. pisum*, for which epigenetic factors have already been annotated [23–27].

In *S. mansoni*, 56 of the 67 previously identified epigenetic factors (84%) could be identified with this methodology (Additional file 3: Table S3). Among the 11 missing

*S. mansoni* proteins, one of them was positioned in the same OrthoMCL group as the human protein SNW1, which is not described to be a chromatin factor (Uniprot accession number: Q13573). Moreover, we were able to identify 36 proteins as new putative chromatin epigenetic factors in *S. mansoni*. In *A. pisum*, 82 of the 99 previously described epigenetic factors (83%) could be identified with our methodology (Additional file 4: Table S4). Four of the 17 missing *A. pisum* proteins were found as orthologs of non-chromatin factors. Moreover, we identified 150 proteins as new putative epigenetic factors of *A. pisum*. Taken together, these results show that this methodology allowed
identification of the majority of known epigenetic factors in two non-model species and identified new putative epi-genetic factors.

When applied to RKN, this annotation methodology identified putative core epigenetic factors orthologous to those of the model nematode *C. elegans* and 15 other species. For each family, numbers of core epigenetic factors were identified in *M. incognita*, and then were compared to those of *C. elegans* (Fig. 2; Additional file 5: Table S5). Annotations for each of the three epigenetic processes are detailed below.

**Cytosine methylation machinery is absent in RKN and in all tested nematodes with the exception of *T. spiralis***

In both plants and animals, methylation on the fifth carbon of cytosine (5mC) is common and has been extensively studied [55]. DNA methyltransferases (DNMT) were searched by orthology links with the 14 reference DNMT (Additional file 1: Table S1 and Additional file 2: Table S2) and on the basis of DNA methyltransferase Pfam domain presence (PF00145). No DNMT could be

identified in any plant-parasitic nematode (PPN). DNMT2 orthologs were identified in *A. suum* and *T. spiralis*, and DNMT3 was restricted to *T. spiralis*. Another putative DNMT (Uniprot reference: H2L057) was identified in *C. elegans* but could not be assigned to any DNMT family. H2L057 was identified by the presence of a DNA methyltransferase Pfam domain (PF00145) and was not associated to any orthoMCL group (“singleton”).

No DNMT1 was identified in any of the nine nematodes studied (*C. elegans*, *M. incognita*, *M. javanica*, *M. arenaria*, *M. hapla*, *G. pallida*, *G. rostochiensis*, *T. spiralis* and *A. suum*). However, according to the literature, DNMT1-like proteins were previously identified in nematodes, including *C. elegans*, *T. spiralis*, *A. suum* and *M. hapla* [16, 56]. We looked for these proteins in our dataset. *A. suum* DNMT1-like could not be found in Uniprot with the accession number given by the authors (GS_06989). *C. elegans* (Q9U1S4) and *T. spiralis* (E5SAG0) DNMT1-like proteins were not assigned to any OrthoMCL group (“singleton”) and did not possess the specific Pfam domain for DNMT (PF00145). *M.*

<figure><img src="image_5.png" /></figure>
hapla (Mhap1s0004g00613) DNMT1-like was positioned in an OrthoMCL group (GRP7217) containing 12 proteins from PPN and none of these proteins possessed the specific Pfam domain for DNMT (PF00145).

We built a phylogeny containing all DNMT, including the putative *C. elegans* DNMT (H2L057) and other nematodes DNMT1-like. Both the putative *C. elegans* DNMT and other nematodes DNMT-1 like grouped outside the known DNMT (Fig. 3a and b). Furthermore, the analysis of the protein domains carried by the DNMT1-like proteins showed they lack a DNA methylase domain (PF00145) and only contain a Pfam Zf-CXXC domain (PF02008) (Fig. 3c). Zf-CXXC domain is known to be associated to the recognition and binding to unmethylated cytosines in CpG. *H. sapiens* DNMT1 carries six protein domains, including a Zf-CXXC domain in position 646–691 and a DNA methylase domain in position 1140–1593. All nematode DNMT1-like could be aligned on a window from 635 to 693 in *H. sapiens* DNMT1, exactly on the Zf-CXXC domain position. However, the alignment of *Globodera* spp. and *C. elegans* DNMT1-like on the DNA methylase domain of *H. sapiens* DNMT1 never exceed 18 aminoacids while this domain size is 453 aminoacids, indicating that PPN, *Meloidogyne* spp. and *Globodera* spp., as well as *C. elegans*, do not possess any possibly functional DNMT. To date, in nematodes, only DNMT2 (*A. suum* and *T. spiralis*) and DNMT3 (*T. spiralis*) orthologs were found. Here, we provide an updated DNMT phylogeny in Fig. 3d.

## Adenine methylation (6mA) machinery is incomplete in RKN

Methylation of DNA on the sixth position of adenine (m6A) is involved in epigenetic transgenerational inheritance [57–59]. DNA N6-methyltransferases (DAMT) were searched by orthology links with the *C. elegans* reference DAMT-1 (Q09956; Table 2) and on the basis of S-adenosylmethionine-binding Pfam domain presence (PF05063). One ortholog of *C. elegans* DAMT-1 gene is found in *M. hapla* whereas multiple co-orthologs are found in the three asexually reproducing RKN species, *M. incognita*, *M. javanica*, *M. arenaria* (Table 2 and Additional file 6: Figure S1). One *M. incognita* ortholog, among the three putative Mi-DAMT-1 identified, was supported by transcriptomic data (Fig. 5). Then, we looked for 6mA demethylases by orthology links with the 5 *C. elegans* AlkB family members (Alkylated DNA repair protein B), together with the presence of AlkB family Pfam domain (PF13532). Proteins of the AlkB family catalyze the demethylation of methylated nucleotides from both DNA and RNA. However, until now, only NMAD-1 (Q8MNT9; Table 2) has been shown with 6mA demethylase activity in *C. elegans* [58]. No ortholog

of *C. elegans* NMAD-1 was found in the Meloidogyne spp. (Table 2 and Additional file 7: Figure S2).

## Histones and histone modifying-enzymes are conserved in RKN

Histones and histone modifying-enzymes were sought by either orthology links with the 631 reference proteins and/or on the basis of at least one of the 21 associated-Pfam domain presence (Additional file 1: Table S1 and Additional file 2: Table S2). This strategy led to the identification of 573 histones and candidate histone modifying-enzymes in *M. incognita*, 730 in *M. javanica*, 771 in *M. arenaria* and 177 in *M. hapla*. All known histones (H1, H2A, H2B, H3 and H4) and families of histone modifying-enzymes (Histone Acetyltransferases, HAT; Histone deacetylases, HDAC; Histone methyl-transferases, HMT; Histone demethylases, HDMT; Histone kinases, HK; Histone phosphatases, HP; Histone ubiquitinyl-transferases, HUT; Histone deubiquitinases, HDU) could be identified (Table 3; Additional file 8: Table S6). Then, we focused on histone acetylation and histone methylation and built phylogenies to associate every *Meloidogyne* spp. putative epigenetic factor to its closest ortholog in *C. elegans* (Fig. 4 and Additional file 9: Figure S3, Additional file 10: Figure S4, Additional file 11: Figure S5, Additional file 12: Figure S6, Additional file 13: Figure S7, Additional file 14: Figure S8, Additional file 15: Figure S9 and Additional file 16: Figure S10). For instance, the phylogenetic tree of GNAT family is shown in Fig. 4. This tree is representative of what was found in all families. Members of the GNAT family were found in all 20 species (Fig. 4a and b). A close up on the NAT10 lineage showed that when one gene is present in *C. elegans*, one ortholog is found in *M. hapla* whereas multiple co-orthologs are found in the three asexually reproducing RKN species, *M. incognita*, *M. javanica*, *M. arenaria* (Fig. 4b and c). In summary, this approach led to the identification of 65 protein lineages that contained at least one *C. elegans* protein. Among these 65 protein lineages, 48 (74%) were associated with at least one epigenetic factor of RKN (Table 3; Additional file 8: Table S6).

To assess the robustness of this annotation, we looked for transcriptional support of these putative histones and histone modifying enzymes based on *M. incognita* RNA-seq data. We found that most (42/65; 65%) of *C. elegans* lineages possessed at least one ortholog in *M. incognita* supported by transcriptomic data (Fig. 5). These results suggest that most of *C. elegans* histone modifying-pathways are conserved and expressed in *M. incognita*.

Furthermore, histone modifying enzymes without orthology to *C. elegans* were identified in RKN. Some of them were orthologs of other model species enzymes (Table 3; Additional file 8: Table S6). For instance, three *M. incognita* proteins were co-orthologs of human
<figure><img src="image_6.png" /><figcaption>**Fig. 3** DNA-methyltransferases (DNMT) and chromomethylases (CMT) annotation. Phylogenetic tree was built with a putative DNMT and CMT sequences from 6 model species and 9 nematodes without a priori, or b putative DNMT and CMT sequences from 6 model species and nematode DNMT1-like orthologs previously identified by Gao et al., 2012. c Alignment of nematode DNMT1-like against *H. sapiens* DNMT1. d Corrected phylogenetic tree for DNMT and CMT sequences from the 20 species. Nematode DNMT are indicated by arrows</figcaption></figure>

CARM1 (Histone-arginine methyltransferase 1). Candidate histone methyltransferases restricted to PPN (*Globodera* spp. and *Meloidogyne* spp.) were also identified based on the presence of the characteristic SET Pfam domain (PF00856)

and the lack of BLASTp hit against the NCBI nr database. They were called SET-PPN because they belong to the SET family of HMT and are, so far, specific to PPN. In *M. incognita*, these SET-PPN were supported by transcriptomic data
<table><caption>**Table 2** Summary of 6mA DNA methyltransferases and demethylases</caption><thead><tr><th rowspan="2">Type</th><th rowspan="2">Family</th><th rowspan="2">Lineage</th><th colspan="2">Ce</th><th colspan="6">RKN</th></tr><tr><th>Number</th><th>Sequences</th><th>Number</th><th>Mi<br>Transcriptional Support</th><th>Mj<br>Number</th><th>Ma<br>Number</th><th>Mh<br>Number</th></tr></thead><tbody><tr><td>Methyltransferase</td><td>MTA70 like</td><td>DAMT-1</td><td>1</td><td>DAMT-1 /Q09956</td><td>3</td><td>X</td><td>3</td><td>4</td><td>1</td></tr><tr><td rowspan="2">Demethylases</td><td rowspan="2">AlkB</td><td>NMAD-1</td><td>1</td><td>NMAD1 /Q8MNT9</td><td>0</td><td></td><td>0</td><td>0</td><td>0</td></tr><tr><td>AlkBH</td><td>3</td><td>alkb-8 /Q9U3P9<br>E0DBL0<br>Q7YWP5</td><td>7</td><td>X</td><td>11</td><td>10</td><td>2</td></tr></tbody></table>

In the model species *C. elegans* (Ce), the number, the names and the Uniprot accession numbers of genes are indicated. Number of orthologs in RKN (*M. incognita*, Mj; *M. javanica*, Mj; *M. arenaria*, Ma; *M. hapla*, Mh) is indicated for each lineage. For *M. incognita*, presence of RNA-seq transcriptional support is indicated by a cross. DAMT, DNA N6-methyl methyltransferase; NMAD, DNA N6-methyl adenine demethylase; AlkB, Alkylated DNA repair protein B; AlkBH, Alkylated DNA repair protein B homolog

(Table 3; Additional file 8: Table S6), suggesting they are functional.

### Small non-coding RNA epigenetic machinery exist in RKN but some components are missing

Small non-coding RNA (ncRNA) play a fundamental role in functional plasticity in eukaryotes [60]. In a general way, small ncRNAs serve as guide to Argonaute proteins (AGO) to regulate their respective targets for gene silencing. AGO proteins are characterized by the presence of PAZ (PF02170) and PIWI (PF02170) domains and can be classified into three clades ([46]; Table 3): (i) the AGO clade, which includes *A. thaliana* AGOs, human AGOs 1–4 and the *C. elegans* miRNA effectors ALG1/2; (ii) the PIWI clade, which includes the *C. elegans* PRG-1 and ERGO-1 (iii) an expanded family of worm-specific AGOs (WAGOs). To identify AGO proteins in RKN, we combined OrthoMCL groupings to the presence of Piwi domain (PF02171), PAZ domain (PF02170) and AGO1 domain (PF08699) [60]. To simplify the construction of the phylogeny, *M. incognita* and *M. hapla* were the only RKN tested because of the high similarity between *M. incognita* and the two other obligatory asexually reproducing RKN (*M. javanica*, *M. arenaria*). A total of 36 putative Argonautes was identified in *M. incognita* and 14 in *M. hapla* (Table 4). All of them could be associated to a specific non-coding RNA pathway, microRNA (miRNA) or small interfering RNA (siRNA), based on phylogeny (Fig. 6a). We could identify orthologs of *C. elegans* Argonautes involved in miRNA (ALG-1/2, HPO-24) and exogenous siRNA (RDE-1, ZK218.8) pathways (Fig. 6b). We could also identify all families of WAGO Argonautes (cytoplasmic WAGOs, nuclear WAGOs and the WAGO involved in self-recognition CSR-1; Additional file 17: Figure S11). Furthermore, WAGOs involved in self-recognition pathway seems amplified in RKN as we could identify 13 CSR-1-like in *M. incognita* and six in *M. hapla* (Table 4 and Additional file 17: Figure S11). All of the identified *M. incognita* Argonautes are expressed except for

cytoplasmic WAGO for which only half of the genes have transcriptional supports (Fig. 5). In contrast, we could not identify any ortholog of *C. elegans* Argonautes involved in endogenous siRNA (Argonautes triggering sperm-enriched siRNA, ALG-3/4) nor Argonautes involved in piwiRNA (piRNA) pathways (Table 4, Fig. 6b and Additional file 18: Figure S12).

To further investigate small ncRNA biogenesis pathways we looked for proteins from complexes that handle small ncRNA after export from the nucleus: the Drosha protein, Dicer protein, and some of the proteins binding to Drosha (PASH-1) or Dicer (DRH1/3). More especially, in *C. elegans* the ERI/DICER complex, composed of ERI-1/3/5, RRF-3, and DICER mediates RNAi processes [61]. Because siRNA can have several different origins, we also looked for more RNA-dependent RNA polymerase (RdRp) that can synthesize siRNA by copying out simple strand RNA [62]. With exception of ERI-3/5, when one gene is present in *C. elegans*, multiple co-orthologs are found in all four RKN (Table 5 and Fig. 7).

### Discussion

Three systems including DNA methylation, post-translational histone modifications and non-coding RNA-associated gene silencing are currently considered to initiate and sustain epigenetic changes [7]. Here we used bioinformatics-driven functional annotations and literature sources to identify epigenetic machinery genes of RKN (*Meloidogyne* spp.) with a particular interest on the asexually reproducing RKN *M. incognita*.

### Absence of DNMT1 in nematodes and presence of DNMT3 restricted to *T. spiralis*

In most cases, cystosine methylation promotes heterochromatin formation and gene silencing [63]. However, cytosine methylation is not heavily present among eukaryotes. Methylcytosines are only present in cryptic proportions in *Drosophila* and totally absent in *C. elegans* [64, 65]. More generally, in nematodes, no cytosine methylation has been reported except during *T. spiralis*
<table><caption>**Table 3** Summary of histones and histone modifying enzymes annotation</caption><thead><tr><th rowspan="2">Type</th><th rowspan="2">Family</th><th rowspan="2">Lineage</th><th colspan="4">Model species</th><th colspan="4">RKN</th></tr><tr><th colspan="2">Ce</th><th colspan="2">Hs</th><th colspan="2">Mi</th><th>Mj</th><th>Ma</th><th>Mh</th></tr><tr><th></th><th></th><th></th><th>Number</th><th>Sequences</th><th>Number</th><th>Sequences</th><th>Number</th><th>Transcriptional Support</th><th>Number</th><th>Number</th><th>Number</th></tr></thead><tbody><tr><td rowspan="5">Histone</td><td>Linker</td><td>H1</td><td>9</td><td>H1L-1/2/3/4/5/6/7/8/, HIS-24</td><td>11</td><td>H1, H1.0/1/2/3/4/5, H1F0, H1t, H1oo, H1x</td><td>2</td><td>X</td><td>7</td><td>10</td><td>2</td></tr><tr><td>Core</td><td>H2A</td><td>5</td><td>HIS-3/35, HTAS-1, HTZ-1, ZC155.2</td><td>18</td><td>H2A.1, H2A.J, H2A.V, H2A.Z, H2A1A/B/D/H/J, H2A2A/B/C, H2A3, H2AB1, H2AFB1, H2AG, H2AX, H2AW</td><td>21</td><td>X</td><td>28</td><td>23</td><td>5</td></tr><tr><td></td><td>H2B</td><td>5</td><td>HIS-4/11/39/41/48</td><td>22</td><td>H2B1A/B/C/D/E/F/G/H/I/K/L/M/N/O, H2B2C/D/E/F, H2B3B, H2BFM, H2BFS, H2BWT</td><td>10</td><td>X</td><td>12</td><td>16</td><td>1</td></tr><tr><td></td><td>H3</td><td>11</td><td>CPAR-1, HCP-3, HIS-2/40/69/70/71/72/73/74, F20D6.9</td><td>6</td><td>H3, H3.1/2/3, H3.1t, H3.3C</td><td>19</td><td>X</td><td>22</td><td>23</td><td>7</td></tr><tr><td></td><td>H4</td><td>1</td><td>HIS-1</td><td>2</td><td>H4, H4G</td><td>9</td><td>X</td><td>11</td><td>9</td><td>0</td></tr><tr><td rowspan="15">Histone acetyltransferase</td><td rowspan="15">GNAT</td><td>ELP3</td><td>1</td><td>ELPC-3</td><td>1</td><td>ELP3</td><td>0</td><td></td><td>0</td><td>0</td><td>0</td></tr><tr><td>F08F8.4</td><td>1</td><td>F08F8.4</td><td>2</td><td>ESCO1/2</td><td>2</td><td>X</td><td>4</td><td>4</td><td>1</td></tr><tr><td>HAT1</td><td>1</td><td>HAT-1</td><td>1</td><td>HAT1</td><td>3</td><td>X</td><td>4</td><td>4</td><td>1</td></tr><tr><td>NAA40</td><td>1</td><td>Y38A10A.7</td><td>1</td><td>NAA40</td><td>2</td><td>X</td><td>2</td><td>3</td><td>1</td></tr><tr><td>NAA50</td><td>1</td><td>F40F4.7</td><td>1</td><td>NAA50</td><td>2</td><td>X</td><td>3</td><td>3</td><td>1</td></tr><tr><td>NAA60</td><td>1</td><td>F30F8.10</td><td>1</td><td>NAA60</td><td>1</td><td></td><td>2</td><td>4</td><td>1</td></tr><tr><td>NAT10</td><td>1</td><td>NATH-10</td><td>1</td><td>NAT10</td><td>6</td><td>X</td><td>3</td><td>5</td><td>1</td></tr><tr><td>PCAF1</td><td>1</td><td>PCAF-1</td><td>1</td><td>KAT2A/B</td><td>3</td><td></td><td>4</td><td>4</td><td>1</td></tr><tr><td rowspan="4">MYS1</td><td>MYS1</td><td>1</td><td>MYS-1</td><td>1</td><td>KAT5</td><td>2</td><td>X</td><td>2</td><td>5</td><td>1</td></tr><tr><td>MYS2</td><td>1</td><td>MYS-2</td><td>1</td><td>KAT8</td><td>3</td><td>X</td><td>5</td><td>4</td><td>1</td></tr><tr><td>MYS4/LSY12 (LSY12)</td><td>1</td><td>LSY-12</td><td>2</td><td>KAT6A/B</td><td>3</td><td>X</td><td>6</td><td>3</td><td>1</td></tr><tr><td>MYS4/LSY12 (MYS4)</td><td>1</td><td>MYS-4</td><td>0</td><td>–</td><td>2</td><td>X</td><td>3</td><td>8</td><td>1</td></tr><tr><td>TF/NF</td><td>TAF1</td><td>1</td><td>TAF-1</td><td>2</td><td>TAF1/1 L</td><td>13</td><td>X</td><td>14</td><td>12</td><td>1</td></tr><tr><td rowspan="6">Histone deacetylase</td><td rowspan="2">Class I</td><td>Class1 (HDA1/3)</td><td>2</td><td>HDA-1/3</td><td>2</td><td>HDAC1/2</td><td>8</td><td>X</td><td>8</td><td>9</td><td>2</td></tr><tr><td>Class1 (HDA2)</td><td>1</td><td>HDA-2</td><td>1</td><td>HDAC3</td><td>4</td><td>X</td><td>3</td><td>3</td><td>1</td></tr><tr><td rowspan="3">Class II</td><td>Class II (HDA4)</td><td>1</td><td>HDA-4</td><td>4</td><td>HDAC4/5/7/9</td><td>4</td><td>X</td><td>3</td><td>4</td><td>1</td></tr><tr><td>Class II (HDAS/6/10)</td><td>1</td><td>HDA-6</td><td>0</td><td>–</td><td>0</td><td></td><td>0</td><td>0</td><td>0</td></tr><tr><td>Class II (HDAS/6/10)</td><td>3</td><td>HDA-5/10, F43G6.17</td><td>2</td><td>HDAC6/10</td><td>3</td><td>X</td><td>5</td><td>2</td><td>1</td></tr><tr><td>Class III</td><td>Class III (SIR2.1)</td><td>1</td><td>SIR-2.1</td><td>1</td><td>SIRT1</td><td>3</td><td></td><td>3</td><td>3</td><td>1</td></tr></tbody></table>
<table><caption>**Table 3** Summary of histones and histone modifying enzymes annotation (Continued)</caption><thead><tr><th rowspan="3">Type</th><th rowspan="3">Family</th><th rowspan="3">Lineage</th><th colspan="4">Model species</th><th colspan="4">RKN</th></tr><tr><th colspan="2">Ce</th><th colspan="2">Hs</th><th colspan="2">Mi</th><th rowspan="2">Mj</th><th rowspan="2">Ma</th><th rowspan="2">Mh</th></tr><tr><th>Number</th><th>Sequences</th><th>Number</th><th>Sequences</th><th>Number</th><th>Transcriptional Support</th></tr></thead><tbody><tr><td rowspan="10">Histone deacetylase</td><td rowspan="2">LSD</td><td>AMX1</td><td>1</td><td>AMX-1</td><td>1</td><td>KDM1B</td><td>0</td><td></td><td>0</td><td>0</td><td>0</td></tr><tr><td>LSD1/SPR5</td><td>2</td><td>LSD-1, SPR-5</td><td>1</td><td>KDM1A</td><td>4</td><td>X</td><td>6</td><td>6</td><td>2</td></tr><tr><td rowspan="4">JMJ</td><td>JMJC1</td><td>1</td><td>JMJC-1</td><td>1</td><td>NO66</td><td>6</td><td>X</td><td>7</td><td>8</td><td>2</td></tr><tr><td>JMJD1/JHDM1 (JHDM1)</td><td>1</td><td>JHDM-1</td><td>2</td><td>KDM2A/B</td><td>3</td><td>X</td><td>3</td><td>5</td><td>0</td></tr><tr><td>JMJD1/JHDM1 (JMJD1)</td><td>2</td><td>JMJD-1.1/1.2</td><td>3</td><td>KDMT7A, PHF2, PHF8</td><td>2</td><td>X</td><td>3</td><td>5</td><td>0</td></tr><tr><td>JMJD2</td><td>1</td><td>JMJD-2</td><td>5</td><td>KDM4A/B/C/D/E</td><td>3</td><td>X</td><td>9</td><td>10</td><td>1</td></tr><tr><td rowspan="2">RBR2</td><td rowspan="2"></td><td>RBR-2</td><td>4</td><td>KDM5A/B/C/D</td><td>4</td><td></td><td>6</td><td>5</td><td>1</td></tr><tr><td>JMJD3/UTX1 (JMJD3)</td><td>3</td><td>JMJD-3.1/3.2/3.3</td><td>0</td><td>–</td><td>0</td><td></td><td>0</td><td>0</td><td>0</td></tr><tr><td rowspan="2">JMJD4/PSR1 (JMJD4)</td><td rowspan="2"></td><td>JMJD3/UTX1 (UTX1)</td><td>1</td><td>UTX-1</td><td>3</td><td>KDM6A/B, UTY</td><td>2</td><td>X</td><td>2</td><td>4</td><td>1</td></tr><tr><td>JMJD4/PSR1 (JMJD4)</td><td>1</td><td>JMJD-4</td><td>0</td><td>–</td><td>0</td><td></td><td>0</td><td>0</td><td>0</td></tr><tr><td rowspan="2">JMJD5</td><td rowspan="2"></td><td>JMJD-5</td><td>1</td><td>KDM8</td><td>2</td><td>X</td><td>2</td><td>2</td><td>1</td></tr><tr><td>JMJD4/PSR1 (PSR1)</td><td>1</td><td>PSR-1</td><td>1</td><td>JMJD6</td><td>2</td><td>X</td><td>2</td><td>3</td><td>1</td></tr></tbody></table>

In the 2 model species (*C. elegans*, Ce; *H. sapiens*, Hs), the number and the names of genes are indicated. Number of orthologs in RKN (*M. incognita*, Mi; *M. javanica*, Mj; *M. arenaria*, Ma; *M. hapla*, Mh) is indicated for each lineage. For *M. incognita*, presence of RNA-seq transcriptional support is indicated by a cross. Substrates for the histone modifying enzymes of model species are indicated. Only the protein lineages tested by phylogeny are shown

development [56]. When present, cytosine methylation marks are established by DNMT. In mammals, DNMT1 and DNMT3 are respectively responsible for methylation maintenance and de novo methylations [66]. DNMT2 is a tRNA methyltransferase involved in RNA-mediated epigenetic inheritance [67, 68]. However, existence of a DNMT2 DNA methyltransferase activity remains unclear and is still under debate [69]. Although potential DNMT1 were identified in nematodes [56], these proteins only share the presence of a CXXC domain (PF02008) with others DNMT1. CXXC domains are involved in non-methylCpG binding and are present in a broad-range of chromatin-binding proteins, including histone modifying enzymes (example: MLL; the Mixed Lineage Leukemia gene) or transcription factors (example: CXXC1; CXXC-type zinc finger protein 1) [70]. For this reason, CXXC domain presence alone could not be used to formally identify DNMTs. Our results suggest that DNMT1 may actually be absent in all tested nematodes, and DNMT3 restricted to *T. spiralis*, the only nematode known to possess cytosine methylation. In mammals, DNMT1 is involved in methylation maintenance and is mostly expressed during adult life. DNMT3 establishes de novo methylation and is highly

expressed through development [66]. In *T. spiralis*, cytosine methylation is specifically associated with some precise developmental steps and is correlated with an increase of DNMT3 transcription. For this reason, DNMT3 presence could be sufficient to explain cytosine methylation in *T. spiralis*. Conversely, absence of DNMT3 suggests absence of cytosine methylation in other nematodes. This is consistent with the absence of methylcytosine in *C. elegans* [65] and the fact that no cytosine methylation was observed in *M. incognita* [16].

### There is no 6mA DNA demethylase counterpart in any RKN

6mA has been proposed to function as an epigenetic mark that carries heritable epigenetic information in eukaryotes, and especially in *C. elegans* [58]. Because *Meloidogyne* spp. possess no, or little, methylcytosine (5mC) DNA, we sought annotation for 6mA DNA machinery. While potential N6-adenine methyltransferases were identified in *Meloidogyne* spp., no demethylase was found. Based on the AlkB family phylogenetic tree and the identification of NMAD-1 ortholog in *G. palida*, a plant-parasitic nematode, we
a

<figure><img src="image_7.png" /></figure>

b

<table><thead><tr><th>Lineage</th><th>Hs</th><th>Dm</th><th>Sc</th><th>Sp</th><th>At</th><th>Ce</th><th>Mi</th><th>Mj</th><th>Ma</th><th>Mh</th><th>Gp</th><th>Gr</th><th>As</th><th>Ts</th><th>Ap</th><th>Am</th><th>Sm</th><th>Lm</th><th>Bc</th><th>Pf</th></tr></thead><tbody><tr><td>ELP3</td><td>ELP3</td><td>ELP3</td><td>ELP3</td><td>ELP3</td><td>HAG3</td><td>ELPC-3</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>1</td><td>0</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>ESCO1/2</td><td>ESCO1/2</td><td>-</td><td>-</td><td>-</td><td>-</td><td>F08F8.4</td><td>2</td><td>4</td><td>4</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>HAT1</td><td>HAT1</td><td>CG02051</td><td>HAT1</td><td>HAT1</td><td>HAG2</td><td>HAT-1</td><td>3</td><td>4</td><td>4</td><td>1</td><td>0</td><td>0</td><td>1</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>0</td><td>0</td></tr><tr><td>NAA40</td><td>NAA40</td><td>CG7593</td><td>HPA2/3, NAT4</td><td>SPCC825.0</td><td>At1g18335</td><td>Y38A10A.7</td><td>2</td><td>2</td><td>3</td><td>1</td><td>1</td><td>1</td><td>0</td><td>1</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td></tr><tr><td>NAA50</td><td>NAA450</td><td>SAN</td><td>-</td><td>-</td><td>F2I11_230</td><td>F40F4.7</td><td>2</td><td>3</td><td>3</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>NAA60</td><td>NAA60</td><td>NAA60</td><td>-</td><td>-</td><td>MCC1, F5E19_140</td><td>F30F8.10</td><td>1</td><td>2</td><td>4</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>NAT10</td><td>NAT10</td><td>NAT10</td><td>NAT10</td><td>NAT10</td><td>At1g10490, At3g57940</td><td>NATH-10</td><td>6</td><td>3</td><td>5</td><td>1</td><td>1</td><td>1</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>PCAF1</td><td>KAT2A/B</td><td>GCN5</td><td>GCN5</td><td>GCN5</td><td>HAG1</td><td>PCAF-1</td><td>3</td><td>4</td><td>4</td><td>1</td><td>2</td><td>1</td><td>1</td><td>1</td><td>3</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>IDM1</td><td>-</td><td>-</td><td>-</td><td>-</td><td>IDM1, At1g05380, At4g14920, At5g36670, At5g36740, At5g63900</td><td>-</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td><td>1</td><td>0</td></tr></tbody></table>

c

<figure><img src="image_8.png" /><figcaption>**Fig. 4** Representative example of histone modifying enzymes phylogenetic tree: GNAT histone acetyltransferase tree. GNAT histone acetyltransferases were identified in 20 species of interest. **a** Phylogenetic tree of all the 20 species GNAT histone acetyltransferases shows that each epigenetic factor protein fit into specific lineages (ELP3, ESCO1/2, HAT1, NAA40, NAA50, NAA60, NAT10, PCAF1 and IDM1). **b** Number of epigenetic factor protein for each lineage identified in the 20 species. **c** Close up of the NAT10 lineage</figcaption></figure>

hypothesized that such demethylase has been lost in RKN. In absence of protein to catalyze the demethylation of methylated DNA, it is unlikely that Meloidogyne spp. display 6mA genomic DNA, unless another protein has taken over the role of 6mA demethylation.

Histone (de)acetylation and (de)methylation machinery is conserved in RKN, and some families have expanded

The nucleosome is a structure composed by a DNA loop wrapped around two copies of each core histone (H2A, H2B, H3A, and H4), stabilized by linker histone (H1). Although core histones are extremely conserved through
<figure><img src="image_9.png" /><figcaption>**Fig. 5** *M. incognita* putative epigenetic factors supported by RNA-seq. Number of *M. incognita* putative epigenetic factors are classified according to *C. elegans* orthology. For each *C. elegans* epigenetic factor, the number of *M. incognita* orthologs is indicated. Dark histograms indicated that the *M. incognita* gene is expressed (RPKM≥ 5) whereas light histograms indicated that there is no RNA-seq support of gene expression (RPKM< 5)</figcaption></figure>

evolution [71], linker histone sequences differ greatly among species [72]. In *M. incognita*, 28 different core histone proteins were identified. This number of core histones is close to *C. elegans* which possesses 22 core histones. By contrast, linker histones were underrepresented in *M. incognita* (two *M. incognita* H1, six *C. elegans* H1) likely because divergence in linker histone sequences makes them difficult to identify.

HAT and HDAC generally possess a wide target range and could catalyze (de)acetylation on different histone positions, or on different proteins [73]. In most cases, histone hyperacetylation results in chromatin decompaction, and thus transcription of neighboring genes. Conversely, histone hypoacetylation promotes gene silencing.

Methylation of histones is more stable than acetylation and phosphorylation marks [74]. For this reason,
<table><caption>**Table 4** Argonautes annotation</caption><thead><tr><th rowspan="2">Type</th><th rowspan="2">Family</th><th rowspan="2">Lineage</th><th colspan="2">C. elegans</th><th colspan="3">RKN</th><th colspan="4">Other nematodes</th></tr><tr><th>Number</th><th>Sequences</th><th>Number</th><th>Transcriptional support</th><th>Mh</th><th>Number</th><th>Gp</th><th>Gr</th><th>As</th><th>Ts</th></tr></thead><tbody><tr><td rowspan="6">Argonaute</td><td rowspan="3">AGO</td><td>miRNA</td><td>3</td><td>ALG-1, ALG-2, HPO-24</td><td>2</td><td>X</td><td>1</td><td>1</td><td>1</td><td>2</td><td>4</td></tr><tr><td>sperm-enriched siRNA</td><td>2</td><td>ALG-3, ALG-4</td><td>0</td><td></td><td>0</td><td>0</td><td>0</td><td>1</td><td>86</td></tr><tr><td>exogenous siRNA</td><td>2</td><td>RDE-1, ZK218.8</td><td>2</td><td>X</td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td rowspan="2">PIWI</td><td>oogenesis-enriched siRNA</td><td>1</td><td>ERGO-1</td><td>0</td><td></td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>piRNA</td><td>1</td><td>PRG-1, (PRG-2)</td><td>0</td><td></td><td>0</td><td>0</td><td>0</td><td>0</td><td>3</td></tr><tr><td rowspan="2">WAGO</td><td rowspan="2">Cytoplasmic WAGO</td><td></td><td>8</td><td>WAGO-1, WAGO-2, WAGO-4, WAGO-5, SAGO-1, SAGO-2, PPW-1, PPW-2</td><td>15</td><td>X</td><td>4</td><td>6</td><td>9</td><td>1</td><td>0</td></tr><tr><td>Nuclear WAGO</td><td>5</td><td>HRDE-1, NRDE-3, WAGO-10, WAGO-11, C14B1.7</td><td>4</td><td>X</td><td>2</td><td>5</td><td>5</td><td>0</td><td>0</td></tr><tr><td colspan="2">Self recognition pathway</td><td></td><td>2</td><td>CSR-1, C04F12.1</td><td>13</td><td>X</td><td>6</td><td>7</td><td>9</td><td>1</td><td>0</td></tr></tbody></table>

In *C. elegans*, the number and the names of Argonautes are indicated. Number of orthologs in *M. incognita* (Mi), *M. hapla* (Mh) and other nematodes of interest (*G. pallida*, *Gp*; *G. rostochiensis*, *Gr*; *A. suum*, *As*; *T. spiralis*, *Ts*) is indicated for each lineage. For *M. incognita*, presence of RNA-seq transcriptional support is indicated by a cross

methylation could be considered as a “longer-term” mark which has been linked to trans-generational epigenetic heredity in *C. elegans* [75, 76]. Methylations on the positions H3K9, H3K36 and H4K20 are associated with gene inactivation while methylations on the positions H3K4 and H3K79 are associated with gene activation. These patterns of (in)activation are well conserved among evolution as they are systematically found by genome-wide approaches, in *A. thaliana* [77], *D. melanogaster* [78] and *C. elegans* [79]. All families of proteins involved in (de)acetylation and (de)methylation were found in RKN.

The fact that most *C. elegans* proteins involved in histone modifications exist in *Meloidogyne* spp. does not necessary mean that these proteins are functional. To address this question, we looked into *M. incognita* RNA-seq data for transcriptomic evidence. We found that among the 65 *C. elegans* proteins conserved in *M. incognita*, 42 possessed at least one *M. incognita* gene supported by transcriptomic data, indicating that, *C. elegans* histone (de)acetyl/(de)methylation pathways are present, conserved and probably functional in *M. incognita*.

A particular interest is transgenerational epigenetic inheritance phenomena observed in *C. elegans* and that could exist in *M. incognita*. For instance, *C. elegans* transgenerational inheritance of longevity [75] involves the H3K4-methyltransferase SET-2 which is found in four expressed genes in *M. incognita*. Another example

is the transgenerational inheritance of fertility [76]. In that case, transgenerational fertility defects were caused by loss of the H3K4-demethylase SPR-5 (four genes in *M. incognita*, among which two are supported by RNA-seq) and could be accelerated by losses of some H3K4/9-methyltransferases or H3K9-demethylases including MET-2, SET-30 and JMJD-2 which were identified in *M. incognita* (respectively: 3, 2 and 3 genes). Describing histone methylation patterns in *M. incognita* (especially H3K4 and H3K9 methylation) could be of particular interest using the ChIP-seq methodology we previously developed [14].

The majority of *C. elegans* proteins could be associated by phylogeny to one *M. hapla* protein whereas a higher number of genes for each histone modifying enzyme family was observed in asexually reproducing RKN (*M. incognita*, *M. javanica* and *M. arenaria*).

Such gene amplifications are in accordance with the duplicated genome structure of asexually reproducing RKN in comparison to the sexual species *M. hapla* [17, 19]. Indeed, most of *C. elegans* HAT/HDAC/HMT/HDMT associate with up to three *M. incognita* genes as expected by the duplicated (likely triploid) genome structure of *M. incognita*. However, four *C. elegans* epigenetic factors (NATH-10, TAF-1, SET-16, JMJC-1) appeared particularly amplified in *M. incognita* because at least six orthologs were identified for each of these 4 genes. NATH-10 is the *C. elegans* ortholog of the human
a

<figure><img src="image_10.png" /></figure>

b

<figure><img src="image_11.png" /></figure>

**Fig. 6** Argonaute phylogenetic tree. Putative Argonaute proteins from 7 nematodes (*C. elegans*, *M. incognita*, *M. hapla*, *G. pallida*, *G. rostochiensis*, *A. suum* and *T. spiralis*) were selected to build (a) whole Argonaute phylogenetic tree. From this whole Argonaute phylogenetic tree, four branches could be distinguished. One correspond to AGO and Piwi family of Argonautes (b) and three to WAGO (cytoplasmic WAGO, nuclear WAGO, self-recognition pathway WAGO) Argonautes (Additional file 17: Figure S11). *M. incognita* and *M. hapla* are colored in green. *C. elegans* proteins are colored in red. Other nematodes are colored in black

gene NATH-10, a GNAT histone acetyltransferase involved in diverse processes such as telomerase regulation [80], DNA damage response [81] and cellular division [82]. In *C. elegans*, NATH-10 is a positive regulator of vulval

induction. A mutation of NATH-10 increases the sperm production in hermaphrodite nematodes [83]. Six NATH-10 orthologs were found in *M. incognita* when only 1 NATH-10 ortholog was found in *M. hapla*. A
<table><caption>**Table 5** Selected components of small non-coding RNA biogenesis pathways</caption><thead><tr><th rowspan="2">Family</th><th rowspan="2">Lineage</th><th colspan="2">C. elegans</th><th colspan="6">RKN</th><th colspan="3">Other nematodes</th></tr><tr><th>Number</th><th>Sequences</th><th>Number</th><th>Transcriptional support</th><th>Mi</th><th>Mh</th><th>Gp</th><th>Gr</th><th>As</th><th>Ts</th></tr></thead><tbody><tr><td rowspan="4">ERI</td><td>ERI1</td><td>1</td><td>ERI1/O44406</td><td>3</td><td>X</td><td></td><td>1</td><td>1</td><td>1</td><td>1</td><td>0</td></tr><tr><td>ERI3</td><td>1</td><td>ERI3/Q9GZI7</td><td>0</td><td></td><td></td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>ERI5</td><td>1</td><td>ERI5/Q95XS0</td><td>0</td><td>X</td><td></td><td>0</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>DICER/ERI4</td><td>1</td><td>DCR1/P34529</td><td>2</td><td>X</td><td></td><td>2</td><td>2</td><td>1</td><td>0</td><td>4</td></tr><tr><td rowspan="4">RdRP</td><td>RRF3</td><td>1</td><td>RRF3/G5EE53</td><td>10</td><td>X</td><td></td><td>3</td><td>9</td><td>3</td><td>2</td><td>1</td></tr><tr><td>RdRP</td><td>3</td><td>EGO1/G5EBQ3<br>RRF1/G5ECM1<br>RRF2/G5EFA8</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>DRH</td><td>1</td><td>DRH1/G5EDI8</td><td>6</td><td>X</td><td></td><td>1</td><td>0</td><td>0</td><td>0</td><td>0</td></tr><tr><td>Drosha</td><td>1</td><td>DRH3/Q93413</td><td>16</td><td>X</td><td></td><td>9</td><td>6</td><td>6</td><td>2</td><td>0</td></tr><tr><td rowspan="2">DGCR8</td><td></td><td>1</td><td>DRSH1/O01326</td><td>6</td><td>X</td><td></td><td>1</td><td>1</td><td>1</td><td>0</td><td>1</td></tr><tr><td></td><td>1</td><td>PASH1/U4PRH5</td><td>3</td><td>X</td><td></td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr></tbody></table>

Number of orthologs in *M. incognita* (Mi), *M. hapla* (Mh) and other nematodes of interest (*G. pallida*, *Gp*; *G. rostochiensis*, *Gr*; *A. suum*, *As*; *T. spiralis*, *Ts*) is indicated for each lineage. For *M. incognita*, presence of RNA-seq transcriptional support is indicated by a cross. *ERI*, Endoribonuclease; *DCR*, Dicer; *RFF*, RNA-dependent RNA polymerase Family; *RdRp*, RNA-dependent RNA polymerase; *DRH*, Dicer-related helicase; *DRSH*, Drosha; *DGCR8*, Microprocessor complex subunit DGR8/PASH, Partner or Drosha

biological explanation to the apparent NATH-10 expansion in *M. incognita* could be functional redundancy of feminizing genes, in the context of *M. incognita* apomictic reproduction. However, only one of the six *M. incognita* NATH-10 genes was supported by transcriptomic data, suggesting either that RNA-seq do not cover *M. incognita* NATH-10 orthologs conditions of expression, or that pseudogeneization events may have occurred. TAF-1 is the main component of the transcription factor TFIID and possesses a histone acetyltransferase activity on histones H3 and H4, and a kinase activity. While TAF-1 regulates a small amount of genes [84, 85], these genes are of crucial importance. For instance, in mammals, TAF1 regulates apoptosis and cellular cycle [85, 86]; and in *C. elegans*, TAF-1 is required for the transcription of most of the embryonic genes [87]. TAF-1 genes are amplified in *M. incognita* (13 genes). Interestingly, nine of these 13 genes were supported by transcriptomic data suggesting that most of them possess a biological role. *C. elegans* SET-16 and its human orthologs KDMT2C/D are members of the MLL-like complex, involved in H3K4 trimethylation during development and hox genes regulation [88–90]. In *C. elegans* SET-16 is also specifically involved in vulval development. *C. elegans* SET-16 is characterized by the presence of four Pfam domains: SET (PF00856), FYRC (PF05965), zf-HC5HC2H (PF13771) and FYRN (PF05964) domains. SET-16 is amplified in *M. incognita* (10 genes), including eight genes supported by transcriptomic data. However, only three *M. incognita* SET-16 orthologs possess at least one of the SET-16 Pfam domains. This observation suggests that *M. incognita*

SET-16 genes are either truncated or wrongly-annotated, or that they could accomplish different functions.

In *C. elegans*, JMJC-1 is an HDMT from JMJ family involved in chromatin modifications and in stress-responses [91]. Six JMJC-1 orthologs were found in *M. incognita* and were supported by transcriptomic data. All but one *M. incognita* JMJC-1 are complete and possess the characteristic Cupin 4 domain (PF08007). The remaining one is truncated in 3' and cannot be extended because it has reached its scaffold end. In *Meloidogyne* spp., a wide variety of stresses have been linked to male differentiation [92–96]. Interestingly, it has been proposed that sex determination in *Meloidogyne* spp. could involve chromatin modification [97]. Because JMJC-1 is involved in both stress responses and chromatin modifications, and is amplified in asexually reproducing RKN (*M. incognita*, *M. javanica*, *M. arenaria*), it could be a candidate of choice to study in the context of sex determination. The human ortholog of JMJC-1, NO66, is known to possess a H3K4 and H3K36-specific demethylase activity [98]. For these reasons, the identification of H3K4 and/or H3K36 patterns in context of female/male development may be of interest to decipher.

### A set of HMT restricted to plant-parasitic nematodes

To refine comparative annotation of epigenetic factors in relation to the parasitic lifestyle of RKN, we included 5 other species with a parasitic lifestyle (*A. pisum*, *L. maculans*, *B. cinerea*, *P. falcipaum* and *S. mansoni*) trying to identify
<figure><img src="image_12.png" /></figure>

common epigenetic signatures. For instance, the same strategy succeeded to identify epigenetic factors involved in virulence of the parasite *L. maculans* based on orthology with HMT in the model fungi *Neurospora crassa* [21]. In *L. maculans*, effector production was shown to be regulated by the addition or deletion of chromatin marks [21].

In our study, a set of HMT was restricted to PPN (only present in *Meloidogyne* spp. and *Globodera* spp.) called thereafter PPN-SET. Because of the presence of a

functional HMT SET domain and transcriptional supports, these PPN-SET may probably possess a biological role in histone methylation, in PPN. However, their roles are unknown. Because these proteins were restricted to PPN, they could be involved in plant-parasitism. Interestingly, SET proteins are known to be involved in pathogenicity in bacteria [99]. However, our analysis should be reinforced by including additional nematode species exhibiting various lifestyles.
Although specific Argonautes involved in endogenous siRNA and piRNA processing are absent, CSR-1 is amplified in plant-parasitic nematodes

Three classes of small ncRNA are generally distinguished in animals: miRNA, siRNA and piRNA [100]. Argonautes involved in endogenous siRNA could not be identified in PPN. Endogenous siRNA that involve ERGO-1 are maternally inherited and required for zygote development [101]. These siRNA were proposed to be involved in the control of overexpressed genes that originates from gene expansion [102]. Endogenous siRNA that involve ALG-3/4 are necessary for spermatogenesis [103]. Although in *M. incognita* males fail to reproduce with females and do not contribute to the genome of the offspring, which can explain the loss of ALG-3/4, these Argonautes are also absent in either the facultative sexual species *M. hapla*, and in the obligatory sexual species *G. pallida* and *G. rostochiensis*. For this reason, absence of endogenous siRNA processing Argonautes suggests more probably a functional diversification of siRNA pathways in PPN. In addition, PRG-1/2 could not be identified in PPN. PRG-1/2 are involved in piRNA processing and act together with CSR-1 to distinguish self (CSR-1) and non-self (PRG-1/2) during gametogenesis [104]. This observation is in accordance with previous study that found that piwi RNAs are absent outside *C. elegans* clade in nematodes, and their function may have been replaced by pathways involving nematode and non-nematode specific RdRPs [105]. Intriguingly, CSR-1 are amplified and expressed in PPN: six genes in *M. hapla*, 13 genes in *M. incognita*, seven genes in *G. pallida* and nine genes in *G. rostochiensis* were found as *C. elegans* CSR-1 orthologs. By contrast, the animal parasitic-nematode *A. suum* only possesses one CSR-1 ortholog. Amplification of genes encoding small ncRNA machinery associated proteins, such RdRPs and DRHs, is also observed in the three asexually reproducing RKN species, *M. incognita*, *M. javanica*, *M. arenaria*. Altogether, the 13 CSR-1-like, 10 RdRP and 22 DRH genes in *M. incognita* represent potential interesting epigenetic regulators.

Conclusion

This analysis provides the first accurate, comprehensive and manually curated information about RKN proteins involved in epigenetic regulations. This analysis describes corresponding genes together with their expression levels in several developmental stages, from the asexually reproducing RKN of major agricultural importance, *M. incognita*. We believe that these functional annotations will be a valuable tool for researchers working in both the field of epigenetics, evolution, host-pathogen interaction and plant parasitism.

Additional files

**Additional file 1: Table S1.** Reference dataset of 716 known epigenetic factors in 6 model species. (XLSX 48 kb)

**Additional file 2: Table S2.** Protein domains (Pfam) specifically associated with epigenetic factors. For each of the 6 model species, number of proteins belonging to DNMT, histones, histone modifying enzymes (Histone Acetyltransferases, HAT; Histone deacetylases, HDAC; Histone methyltransferases, HMT; Histone demethylases, HDMT; Histone kinases, HK; Histone phosphatases, HP; Histone ubiquitinyl-transferases, HUT; Histone deubiquitinases, HDU), and Argonautes are indicated, as well as the specific protein domains (Pfam) that constitute the "Pfam epigenetic factor list". (XLSX 13 kb)

**Additional file 3: Table S3.** Epigenetic factors of *S. mansoni*. A cross indicates if the epigenetic factor has been previously identified in the literature and/or in the present study. (XLSX 15 kb)

**Additional file 4: Table S4.** Epigenetic factors of *A. pisum*. A cross indicates if the epigenetic factor has been previously identified in the literature and/or in the present study. (XLSX 25 kb)

**Additional file 5: Table S5.** Number of epigenetic factors, classified by process, identified in root-knot nematodes compared to *C. elegans*. Number of epigenetic factors for each families is detailed according to the method of identification: based on either sequence similarity (orthoMCL); or presence of specifically associated protein domain (Pfam); and after validation by phylogeny (validated). ND: Not Determined. (XLSX 15 kb)

**Additional file 6: Figure S1.** Phylogenetic tree of 6mA methyltransferases. (PPTX 81 kb)

**Additional file 7: Figure S2.** Phylogenetic tree of 6mA demethylases. (PPTX 118 kb)

**Additional file 8: Table S6.** Complete list of histones and histone modifying enzymes annotation. In the 6 model species (*C. elegans*, *H. sapiens*, *D. melanogaster*, *S. cerevisae*, *S. pombe* and *A. thaliana*), the number and the names of genes are indicated. Number of orthologs in RKN (*M. incognita*, *M. javanica*, *M. arenaria* and *M. hapla*) and other species of interest (*G. pallida*, *G. rostochiensis*, *A. suum*, *T. spiralis*, *A. pisum*, *A. mellifera*, *S. mansoni*, *B. cinerea*, *L. maculans* and *P. falciparum*) is indicated for each lineage. For *M. incognita*, presence of RNA-seq transcriptional support is indicated by a cross. Substrates for the histone modifying enzymes of model species are indicated. Protein lineages in grey were not tested by phylogeny. ND: Not determined. (XLSX 34 kb)

**Additional file 9: Figure S3.** Phylogenetic tree of MYST family Histone Acetyltransferases. (PPTX 102 kb)

**Additional file 10: Figure S4.** Phylogenetic tree of TF/NF family Histone Acetyltransferases. (PPTX 83 kb)

**Additional file 11: Figure S5.** Phylogenetic tree of Histone Deacetylases. (PPTX 115 kb)

**Additional file 12: Figure S6.** Phylogenetic tree of PRMT family Histone Methyltransferases. (PPTX 96 kb)

**Additional file 13: Figure S7.** Phylogenetic tree of Euchromatin SET Histone Methyltransferases 1. (PPTX 66 kb)

**Additional file 14: Figure S8.** Phylogenetic tree of Euchromatin SET Histone Methyltransferases 2. (PPTX 111 kb)

**Additional file 15: Figure S9.** Phylogenetic tree of Heterochromatin SET Histone Methyltransferases. (PPTX 103 kb)

**Additional file 16: Figure S10.** Phylogenetic tree of Histone Demethylases. (PPTX 119 kb)

**Additional file 17: Figure S11.** Phylogenetic tree of WAGO Argonautes. Putative Argonaute proteins from 7 nematodes (*C. elegans*, *M. incognita*, *M. hapla*, *G. pallida*, *G. rostochiensis*, *A. suum* and *T. spiralis*) were selected to build the whole argonaute tree (Fig. 3). From this whole Argonaute tree, three branches corresponded to WAGO (cytoplasmic WAGO, nuclear WAGO, self-recognition pathway WAGO) Argonautes. *M. incognita* and *M. hapla* are colored in green. *C. elegans* proteins are colored in red. (PPTX 165 kb)

**Additional file 18: Figure S12.** Phylogenetic tree of PIWI Argonautes. (PPTX 83 kb)
## Abbreviations

AlkB: Alkylated DNA repair protein B; CARM1: Histone-arginine methyltransferase 1; CBP/P300: CREB-binding protein and P300; CMT: Cromomethylase; DAMT: DNA N6-methyltransferase; DNMT: DNA-methyltransferase; DOT1: Disruptor of telomeric 1; DRH: Dicer related helicase; GNAT: Gcn5-related N-acetyltransferase; HAT: Histone acetyltransferase; HDAC: Histone deacetylase; HDMT: Histone demethylase; HK: Histone kinase; HMT: Histone methyltransferase; HP: Histone phosphatase; HUT: Histone ubiquitin-transferase; JMJ: Jumonji demethylase; LSD: Lysine-specific demethylase; miRNA: Micro RNA; MYST: MOZ/Ybf2/Sas3/Sas2/TIP60; ncRNA: non coding RNA; NMAD: DNA N6-methyl adenine demethylase; PCAF: P300/CBP-associated factor; piRNA: PIWI-interacting RNA; PPN: Plant-parasitic nematodes; PRMT: Protein arginin methyltransferase; PTMs: Post-translational modifications; RdRp: RNA-dependent RNA polymerase; RKN: Root-knot nematodes; RPKM: Reads per kilobase of transcript per million mapped reads; SET: Suppressor of variegation/enhancer of zeste/trithorax; siRNA: Small interfering RNA; WAGO: Worm-specific argonaute

## Acknowledgements

We are grateful to the Genotoul bioinformatics platform Toulouse Midi-Pyrenees (Bioinfo Genotoul) for providing computing resources; Chantal Castagnone for her assistance in the maintenance of experimental stocks; Bruno Favery for intensive reading of the manuscript and discussions; and we are especially thankful to Christoph Grunau and Gael Cristofari for sharing their expertise during Loris Pratx committee thesis meetings.

## Funding

L.P. was supported by the French Government (National Research Agency, ANR) through the "Investments for the Future" LABEX SIGNALIFE: program reference # ANR-11-LABX-0028-01. This project was supported by the ANR program ANR-13-JSV7-0006 - ASEXEVOL and a grant from the Division for Plant Health and Environment (SPE) of the French National Institute for Agricultural Research (INRA) to P.C. and L.P.B. The funding bodies did not have a role in the design of the study, data collection, analysis, interpretation of data, writing the manuscript, nor the decision to publish.

## Availability of data and materials

The nucleotide and protein sequences/annotations for *M. incognita* epigenetic factors together with the RNA-seq reads supporting the conclusions of this article are available at "Meloidogyne genomic resources" website: https://meloidogyne.inra.fr.

## Authors' contributions

LP and LPB conceived, designed and performed the experiments. EGJD participated to conception and design of phylogenetic analysis. MDR and CR performed the experiments. LP and LPB analyzed the data. LP, LPB, EGJD, PCS and PA contributed to the writing of the paper. LPB and PA supervised the project. All authors read and approved the final manuscript.

## Ethics approval and consent to participate

Not applicable.

## Competing interests

The authors declare that they have no competing interests.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Received: 11 April 2017 Accepted: 16 April 2018

Published online: 03 May 2018

## References

1. Russo V, Martienssen R, Riggs A. Epigenetic mechanisms of gene regulation. In: Riggs A, Martienssen R, Russo V, editors. Cold Spring Harbor laboratory press. New York: Cold spring harbor; 1996.

2. Kosak ST, Goudine M. Gene order and dynamic domains. Science. 2004; 306:644-7.

3. Luger K, Mäder AW, Richmond RK, Sargent DF, Richmond TJ. Crystal structure of the nucleosome core particle at 2.8 a resolution. Nature. 1997; 389:251-60.

4. Grewal SIS, Jia S. Heterochromatin revisited. Nat Rev Genet. 2007;8:35-46.

5. Jablonka E, Lamb MJ. The changing concept of epigenetics. Ann N Y Acad Sci. 2002;981:82-96.

6. Breiling A, Lyko F. Epigenetic regulatory functions of DNA modifications: 5-methylcytosine and beyond. Epigenetics Chromatin. 2015;8:24. https://doi.org/10.1186/s13072-015-0016-6. eCollection 2015. Review.

7. Liebers R, Rassoulzadegan M, Lyoko F. Epigenetic regulation by heritable RNA. PLoS Genet. 2014;10:e1004296.

8. Aguilera O, Fernández AF, Muñoz A, Fraga MF. Epigenetics and environment: a complex relationship. J Appl Physiol. 2010;109:243-51.

9. Hauser M-T, Aufsatz W, Jonak C, Luschnig C. Transgenerational epigenetic inheritance in plants. Biochim Biophys Acta. 1809;2011:459-68.

10. Bayer-Santos E, Marini MM, da Silveira JF. Non-coding RNAs in host-pathogen interactions: subversion of mammalian cell functions by protozoan parasites. Front Microbiol. 2017;8:1-8.

11. Zhu QH, Shan WX, Ayliffe MA, Wang MB. Epigenetic mechanisms: an emerging player in plant-microbe interactions. Mol Plant-Microbe Interact. 2016;29(3):187-96.

12. Castagnone-Sereno P, Danchin EGJ, Perfus-Barbeoch L, Abad P. Diversity and evolution of root-knot nematodes, genus Meloidogyne: new insights from the genomic era. Annu Rev Phytopathol. 2013;51:203-20.

13. Castagnone-Sereno P, Wajnberg E, Bongiovanni M, Leroy F, Dalmasso A. Genetic variation in Meloidogyne incognita virulence against the tomato mi resistance gene : evidence from isofemale line selection studies. Theor Appl Genet. 1994;88:749-53.

14. Bost SC, Triantaphyllou AC. Genetic basis of the epidemiologic effects of resistance to Meloidogyne incognita in the tomato cultivar small fry. J Nematol. 1982;14:540-4.

15. Jarquin-Barberena H, Dalmasso A, de Guiran G, Cardin MC. Acquired virulence in the plant parasitic nematode *Meloidogyne incognita*. *I*. Biological analysis of the phenomenon. *Rev Nématol*. 1991;14:299-303.

16. Perfus-Barbeoch L, Castagnone-Sereno P, Reichelt M, Fneich S, Roquis D, Pratx L, et al. Elucidating the molecular bases of epigenetic inheritance in non-model invertebrates: the case of the root-knot nematode Meloidogyne incognita. Front Physiol. 2014;5:211.

17. Abad P, Gouzy J, Aury J-M, Castagnone-Sereno P, Danchin EGJ, Deleury E, et al. Genome sequence of the metazoan plant-parasitic nematode Meloidogyne incognita. Nat Biotechnol. 2008;26:909-15.

18. Opperman CH, Bird DM, Williamson VM, Rokhsar DS, Burke M, Cohn J, et al. Sequence and genetic map of Meloidogyne hapla: a compact nematode genome for plant parasitism. Proc Natl Acad Sci U S A. 2008;105:14802-7.

19. Blanc-Mathieu R, Perfus-Barbeoch L, Aury JM, Da Rocha M, Gouzy J, Sallet E, Martin-Jimenez C, Bailly-Bechet M, Castagnone-Sereno P, Flot JF, Kozlowski DK, Cazareth J, Couloux A, Da Silva C, Guy J, Kim-Jo YJ, Rancurel C, Schiex T, Abad P, Wincker P, Danchin EGJ. Hybridization and polyploidy enable genomic plasticity without sex in the most devastating plant-parasitic nematodes. PLoS Genet. 2017;13(6):e1006777.

20. Medvedeva YA, Lennartsson A, Ehsani R, Kulakovskiy IV, Vorontsov IE, Panahandeh P, et al. EpiFactors: a comprehensive database of human epigenetic factors and complexes. Database. 2015;2015:bav067.

21. Soyer JL, El Ghalid M, Glaser N, Ollivier B, Linglin J, Grandaubert J, et al. Epigenetic control of effector gene expression in the plant pathogenic fungus Leptosphaeria maculans. Talbot NJ, editor. PLoS Genet. 2014;10:e1004227.

22. Bin Terhem R, van JAL K. Dual mating in Botrytis cinerea. In: Book of abstracts 10th international mycological congress; 2014. p. 403.

23. Rider SD, Srinivasan DG, Hilgarth RS. Chromatin-remodelling proteins of the pea aphid, Acyrthosiphon pisum (Harris). Insect Mol Biol. 2010;19:201-14.

24. Rasmussen EMK, Ardam GV. Cytosine modifications in the honey bee (Apis mellifera) worker genome. Front. Genet. 2015;6:1-5.

25. Anderson L, Pierce RJ, Verjovski-Almeida S. Schistosoma mansoni histones: from transcription to chromatin regulation; an in silico analysis. Mol Biochem Parasitol. 2012;183:105-14.

26. Mourão MM, Grunau C, LoVerde PT, Jones MK, Oliveira G. Recent advances in Schistosoma genomics. Parasite Immunol. 2012;34:151-62.

27. Bertin B, Oger F, Cornette J, Caby S, Noël C, Capron M, et al. Schistosoma mansoni CBP/p300 has a conserved domain structure and interacts functionally with the nuclear receptor SmFtz-F1. Mol Biochem Parasitol. 2006;146:180-91.
28. Ay F, Bunnik EM, Varoquaux N, Vert J-P, Noble WS, Le Roch KG. Multiple dimensions of epigenetic gene regulation in the malaria parasite plasmodium falciparum. *BioEssays*. 2015;37:182–94.

29. The UniProt Consortium. UniProt: a hub for protein information. *Nucleic Acids Res*. 2015;43:D204–12.

30. Howe KL, Bolt BJ, Cain S, Chan J, Chen WJ, Davis P, et al. WormBase 2016: expanding to enable helminth genomic research. *Nucleic Acids Res*. 2016; 44:D774–80.

31. Howe KL, Bolt BJ, Shafie M, Kersey P, Berriman M. WormBase ParaSite – a comprehensive resource for helminth genomics. *Mol Biochem Parasitol*. 2016; https://doi.org/10.1016/j.molbiopara.2016.11.005.

32. Li L, Jr CJS, Roos DS. OrthoMCL: identification of Ortholog groups for eukaryotic genomes. *Genome Res*. 2003;13:2178–89.

33. Finn RD, Coggill P, Eberhardt RY, Eddy SR, Mistry J, Mitchell AL, Potter SC, Punta M, Qureshi M, Sangrador-Vegas A, Salazar GA, Tate J, Bateman A. The Pfam protein families database: towards a more sustainable future. *Nucleic Acids Res*. 2016;44(D1):D279–D285. d.

34. Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. *Bioinformatics*. 2014;30:1236–40.

35. Cottet L, Rancurel C, Briand M, Carrere S. Family-companion: analyse, visualise, browse, query and share your homology clusters. *bioRxiv* 266742; doi: https://doi.org/10.1101/266742.

36. Harris TW, Baran J, Bieri T, Cabunoc A, Chan J, Chen WJ, et al. WormBase 2014: new views of curated biology. *Nucleic Acids Res*. 2014;42:789–93.

37. Mariño-Ramírez L, Levine KM, Morales M, Zhang S, Moreland RT, Baxevanis AD, et al. The histone database: an integrated resource for histones and histone fold-containing proteins. *Database*. 2011; https://doi.org/10.1093/database/bar048.

38. Liu L, Zhen XT, Denton E, Marsden BD, Schapira M. ChromoHub: a data hub for navigators of chromatin-mediated signalling. *Bioinformatics*. 2012;28:2205–6.

39. Khare SP, Habib F, Sharma R, Gadewal N, Gupta S, Galande S. Histome-a relational knowledgebase of human histone proteins and histone modifying enzymes. *Nucleic Acids Res*. 2012;40:D337–42.

40. Banerjee T, Chakravarti D. A peek into the complex realm of histone phosphorylation. *Mol Cell Biol*. 2011;31:4858–73.

41. Di Lorenzo A, Bedford MT. Histone arginine methylation. *FEBS Lett*. 2011; 585:2024–31. Federation of European Biochemical Societies

42. Finnegan EJ, Kovac KA. Plant DNA methyltransferases. *Plant Mol Biol*. 2000; 43:189–201.

43. Pikaard C, Mittelsten Scheid O. Epigenetic regulation in plants. *Cold Spring Harb Perspect Biol*. 2014;6:1–31.

44. Rossetto D, Avvakumov N, Côté J. Histone phosphorylation: a chromatin modification involved in diverse nuclear events. *Epigenetics*. 2012;7: 1098–108.

45. Simó-Riudalbas L, Esteller M. Targeting the histone orthography of cancer: drugs for writers, erasers and readers. *Br J Pharmacol*. 2014;172:2716–32.

46. Yigit E, Batista PJ, Bei Y, Pang KM, Chen CCG, Tolia NH, et al. Analysis of the C. Elegans Argonaute family reveals that distinct Argonautes act sequentially during RNAi. *Cell*. 2006;127:747–57.

47. Punta M, Coggill P, Eberhardt R, Mistry J, Tate J, Boursnell C, et al. The Pfam protein families databases. *Nucleic Acids Res*. 2012;40:D290–301. 30:1–12

48. Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. *Nucleic Acids Res*. 2014;42:D222–30.

49. Finn RD, Coggill P, Eberhardt RY, Eddy SR, Mistry J, Mitchell AL, et al. The Pfam protein families database: towards a more sustainable future. *Nucleic Acids Res*. 2016;44:D279–85.

50. Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. *Mol Biol Evol*. 2013;30:772–80.

51. Capella-Gutiérrez S, Silla-Martínez JM, Gabaldón T. trimAI: a tool for automated alignment trimming in large-scale phylogenetic analyses. *Bioinformatics*. 2009;25:1972–3.

52. Stamatakis A. RAxML version 8 a tool for phylogenetic analysis and post-analysis of large phylogenies. *Bioinformatics*. 2014;30:2010–1.

53. Rambaut A. FigTree, a graphical viewer of phylogenetic trees. 2007. http://tree.bio.ed.ac.uk/software/figtree. Accessed 9 July 2014.

54. Danchin EGJ, Arguel M-J, Campan-Fournier A, Perfus-Barbeoch L, Magliano M, Rosso M-N, et al. Identification of novel target genes for safer and more specific control of root-knot nematodes from a pan-genome mining. *PLoS Pathog*. 2013;9:e1003745.

55. Lee T, Zhai J, Meyers BC. Conservation and divergence in eukaryotic DNA methylation. *PNAS*. 2010;107(20):9027–8.

56. Gao F, Liu X, Wu X-P, Wang X-L, Gong D, Lu H, et al. Differential DNA methylation in discrete developmental stages of the parasitic nematode Trichinella spiralis. *Genome Biol*. 2012;13:R100.

57. Luo G-Z, Blanco MA, Greer EL, He C, Shi Y. DNA N6-Methyladenine: a new epigenetic mark in eukaryotes? *Nat Rev Mol Cell Biol*. 2015;16(12):705–10.

58. Greer EL, Blanco MA, Gu L, Sendinc E, Liu J, Aristizábal-Corrales D, Hsu CH, Aravind L, He C, Shi Y. DNA methylation on N6-adenine in *C. elegans*. *Cell*. 2015;161(4):868–78.

59. O'Brown ZK, Greer EL. N6-methyladenine: a conserved and dynamic DNA mark. *Adv Exp Med Biol*. 2016;945:213–46. https://doi.org/10.1007/978-3-319-43624-1_10.

60. Meister G. Argonaute proteins: functional insights and emerging roles. *Nat. Rev. Genet*. 2013;14:447–59.

61. Pavelec DM, Lachowiec J, Duchaine TF, Smith HE, Kennedy S. Requirement for the ERI/DICER complex in endogenous RNA interference and sperm development in Caenorhabditis elegans. *Genetics*. 2009;183(4):1283–95.

62. Billi AC, Fischer SE, Kim JK. Endogenous RNAi pathways in C. Elegans. *WormBook*. 2014;7:1–49.

63. Tariq M, Paszkowski J. DNA and histone methylation in plants. *Trends Genet*. 2004;20:244–51.

64. Salzberg A, Fisher O, Siman-Tov R, Ankri S. Identification of methylated sequences in genomic DNA of adult Drosophila melanogaster. *Biochem Biophys Res Commun*. 2004;322:465–9.

65. Suzuki MM, Bird A. DNA methylation landscapes: provocative insights from epigenomics. *Nat. Rev. Genet*. 2008;9:465–76.

66. Chen Z, Riggs AD. DNA methylation and demethylation in mammals. *J Biol Chem*. 2011;286:18347–53.

67. Kiani J, Grandjean V, Liebers R, Tuorto F, Ghanbarian H, Lyko F, et al. RNA-mediated epigenetic heredity requires the cytosine methyltransferase Dnmt2. *PLoS Genet*. 2013;9:e1003498.

68. Lyko F. The DNA methyltransferase family: a versatile toolkit for epigenetic regulation. *Nat Rev Genet*. 2018;19(2):81-92. https://doi.org/10.1038/nrg.2017.80. Epub 2017 Oct 16. Review.

69. Jeltsch A, Ehrenhofer-Murray A, Jurkowski T, Lyko F, Reuter G, Ankri S, et al. Mechanism and biological role of Dnmt2 in nucleic acid methylation. *RNA Biol*. 2016;6286:1–16.

70. Allen MD, Grummitt CG, Hilcenko C, Min SY, Tonkin LM, Johnson CM, et al. Solution structure of the nonmethyl-CpG-binding CXXC domain of the leukaemia-associated MLL histone methyltransferase. *EMBO J*. 2006;25: 4503–12.

71. Thatcher TH, Gorovsky M. A. Phylogenetic analysis of the core histones H2A, H2B, H3, and H4. *Nucleic Acids Res*. 1994;22:174–9.

72. Böhm L, Mitchell TC. Sequence conservation in the N-terminal domain of histone H1. *FEBS Lett*. 1985;193:1–4.

73. Kouzarides T. Chromatin modifications and their function. *Cell*. 2007;128: 693–705.

74. Sims RJ, Nishioka K, Reinberg D. Histone lysine methylation: a signature for chromatin function. *Trends Genet*. 2003;19:629–39.

75. Greer EL, Maures TJ, Ucar D, Hauswirth AG, Mancini E, Lim JP, et al. Transgenerational epigenetic inheritance of longevity in Caenorhabditis elegans. *Nature*. 2011;479:365–71.

76. Greer EL, Beese-Sims SE, Brookes E, Spadafora R, Zhu Y, Rothbart SB, et al. A histone methylation network regulates transgenerational epigenetic memory in C.Elegans. *Cell Rep*. 2014;7:113–26.

77. Roudier F, Ahmed I, Sarazin A, Mary-huard T, Cortijo S, Bouyer D, et al. Integrative epigenomic mapping defines four main chromatin states in Arabidopsis. *EMBO J*. 2011;30:1928–38.

78. Filion GJ, Van Bemmel JG, Braunschweig U, Talhou W, Kind J, Ward LD, et al. Systematic protein location mapping reveals five principal chromatin types in Drosophila cells. *Cell*. 2010;143:212–24.

79. Liu T, Rechtsteiner A, Egelhofer TA, Vielle A, Latorre I, Cheung M, et al. Broad chromosomal domains of histone modification patterns in C. Elegans. *Genome Res*. 2011;22:227–227.

80. Lv J, Liu H, Wang Q, Tang Z, Hou L, Zhang B. Molecular cloning of a novel human gene encoding histone acetyltransferase-like protein involved in transcriptional activation of hTERT. *Biochem Biophys Res Commun*. 2003; 311:506–13.

81. Liu H, Ling Y, Gong Y, Sun Y, Hou L, Zhang B. DNA damage induces N-acetyltransferase NAT10 gene expression through transcriptional activation. *Mol Cell Biochem*. 2007;300:249–58.
82. Shen Q, Zheng X, McNutt M a, Guang L, Sun Y, Wang J, et al. NAT10, a nucleolar protein, localizes to the midbody and regulates cytokinesis and acetylation of microtubules. Exp Cell Res. 2009;315:1653-67.

83. Duveau F, Félix M-A. Role of pleiotropy in the evolution of a cryptic developmental variation in Caenorhabditis elegans. PLoS Biol. 2012;10:e1001230.

84. Lee TI, Causton HC, Holstege FC, Shen WC, Hannett N, Jennings EG, et al. Redundant roles for the TFIID and SAGA complexes in global transcription. Nature. 2000;405:701-4.

85. O'Brien T, Tjian R. Different functional domains of TAFII250 modulate expression of distinct subsets of mammalian genes. Proc Natl Acad Sci U S A. 2000;97:2456-61.

86. Dunphy EL, Johnson T, Auerbach SS, Wang EH. Requirement for TAF(II)250 acetyltransferase activity in cell cycle progression. Mol Cell Biol. 2000;20:1134-9.

87. Walker AK, Shi Y, Blackwell TK. An extensive requirement for transcription factor IID-specific TAF-1 in Caenorhabditis elegans embryonic transcription. J Biol Chem. 2004;279:15339-47.

88. Andersen EC, Horvitz HR. Two C. Elegans histone methyltransferases repress lin-3 EGF transcription to inhibit vulval development. Development. 2007;134:2991-9.

89. Fisher K, Southall SM, Wilson JR, Poulin GB. Methylation and demethylation activities of a C. Elegans MLL-like complex attenuate RAS signalling. Dev. Biol. 2010;341:142-53.

90. Greer EL, Maures TJ, Hauswirth AG, Green EM, Leeman DS, Maro GS, et al. Members of the H3K4 trimethylation complex regulate lifespan in a germline-dependent manner in C. Elegans. Nature. 2010;466:383-7.

91. Kirienko NV, Fay DS. SLR-2 and JMJC-1 regulate an evolutionarily conserved stress-response network. EMBO J. 2010;29:727-39.

92. Moura RM, Davis EL, Luzzi BM, Boerma HR, Hussey RS. Post-Infectional development of Meloidogyne incognita on susceptible and resistant soybean genotypes. Nematropica. 1993;23:7-13.

93. Triantaphyllou AC. Environmental sex differentiation of nematodes in relation to pest management. Annu Rev Phytopathol. 1973;11:441-62.

94. Davide RG, Triantaphyllou AC. Influence of the environment on development and sex differenciation of root-knot nematodes. Nematologica. 1968;14:37-46.

95. Snyder DW, Opperman CH, Bird DM. A method for generating Meloidogyne incognita males. J Nematol. 2006;38:192-4.

96. Davide RG, Triantaphyllou AC. Influence of the environment on development and sex differentiation of root-knot nematodes. Nematologica. 1967;13:102-10.

97. Goldstein P, Triantaphyllou AC. Karyotype analysis of Meloidogyne hapla by 3-D reconstruction of synaptonemal complexes from electron microscopy of serial sections. Chromosoma. 1978;70:131-9.

98. Kooistra SM, Helin K. Molecular mechanisms and potential functions of histone demethylases. Nat. Rev. Mol. Cell Biol. 2012;13:297-311.

99. Pennini ME, Perrinet S, Dautry-Varsat A, Subtil A. Histone methylation by NUE, a novel nuclear effector of the intracellular pathogen chlamydia trachomatis. PLoS Pathog. 2010;6:1-12.

100. Ha M, Kim VN. Regulation of microRNA biogenesis. Nat Rev Mol Cell Biol. 2014;15:509-24.

101. Han T, Manoharan AP, Harkins TT, Bouffard P, Fitzpatrick C, Chu DS, et al. 26G endo-siRNAs regulate spermatogenic and zygotic gene expression in Caenorhabditis elegans. Proc Natl Acad Sci U S A. 2009;106:18674-9.

102. Vasale JJ, Gu W, Thivierge C, Batista PJ, Claycomb JM, Youngman EM, et al. Sequential rounds of RNA-dependent RNA transcription drive endogenous small-RNA biogenesis in the ERGO-1/Argonaute pathway. Proc Natl Acad Sci U S A. 2010;107:3582-7.

103. Conine CC, Batista PJ, Gu W, Claycomb JM, Chaves D a, Shirayama M, et al. Argonautes ALG-3 and ALG-4 are required for spermatogenesis-specific 26G-RNAs and thermotolerant sperm in Caenorhabditis elegans. Proc Natl Acad Sci U S A. 2010;107:3588-93.

104. Shirayama M, Seth M, Lee HC, Gu W, Ishidate T, Conte D, et al. PiRNAs initiate an epigenetic memory of nonself RNA in the C. Elegans germline. Cell. 2012;150:65-77.

105. Sarkies P, Selkirk ME, Jones JT, Blok V, Boothby T, Goldstein B, Hanelt B, Ardila-Garcia A, Fast NM, Schiffer PM, Kraus C, Taylor MJ, Koutsovoulos G, Blaxter ML, Miska EA. Ancient and novel small RNA pathways compensate for the loss of piRNAs in multiple independent nematode lineages. PLoS Biol. 2015;13(2):e1002061.

Ready to submit your research? Choose BMC and benefit from:

• fast, convenient online submission

• thorough peer review by experienced researchers in your field

• rapid publication on acceptance

• support for research data, including large and complex data types

• gold Open Access which fosters wider collaboration and increased citations

• maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

BMC