Elucidating the molecular bases of epigenetic inheritance in non-model invertebrates: the case of the root-knot nematode Meloidogyne incognita

Chromatin Landscape Dynamics in the Early Development of the Plant Parasitic Nematode Meloidogyne incognita

In model organisms, epigenome dynamics underlies a plethora of biological processes. The role of epigenetic modifications in development and parasitism in nematode pests remains unknown. The root-knot nematode Meloidogyne incognita adapts rapidly to unfavorable conditions, despite its asexual reproduction. However, the mechanisms underlying this remarkable plasticity and their potential impact on gene expression remain unknown. This study provides the first insight into contribution of epigenetic mechanisms to this plasticity, by studying histone modifications in M. incognita. The distribution of five histone modifications revealed the existence of strong epigenetic signatures, similar to those found in the model nematode Caenorhabditis elegans. We investigated their impact on chromatin structure and their distribution relative to transposable elements (TE) loci. We assessed the influence of the chromatin landscape on gene expression at two developmental stages: eggs, and pre-parasitic juveniles. H3K4me3 histone modification was strongly correlated with high levels of expression for protein-coding genes implicated in stage-specific processes during M. incognita development. We provided new insights in the dynamic regulation of parasitism genes kept under histone modifications silencing. In this pioneering study, we establish a comprehensive framework for the importance of epigenetic mechanisms in the regulation of the genome expression and its stability in plant-parasitic nematodes.

The Centromere Histone Is Conserved and Associated with Tandem Repeats Sharing a Conserved 19-bp Box in the Holocentromere of Meloidogyne Nematodes

Although centromeres have conserved function, centromere-specific histone H3 (CenH3) and centromeric DNA evolve rapidly. The centromere drive model explains this phenomenon as a consequence of the conflict between fast-evolving DNA and CenH3, suggesting asymmetry in female meiosis as a crucial factor. We characterized evolution of the CenH3 protein in three closely related, polyploid mitotic parthenogenetic species of the Meloidogyne incognita group, and in the distantly related meiotic parthenogen Meloidogyne hapla. We identified duplication of the CenH3 gene in a putative sexual ancestral Meloidogyne. We found that one CenH3 (αCenH3) remained conserved in all extant species, including in distant Meloidogyne hapla, whereas the other evolved rapidly and under positive selection into four different CenH3 variants. This pattern of CenH3 evolution in Meloidogyne species suggests the subspecialization of CenH3s in ancestral sexual species. Immunofluorescence performed on mitotic Meloidogyne incognita revealed a dominant role of αCenH3 on its centromere, whereas the other CenH3s have lost their function in mitosis. The observed αCenH3 chromosome distribution disclosed cluster-like centromeric organization. The ChIP-Seq analysis revealed that in M. incognita αCenH3-associated DNA dominantly comprises tandem repeats, composed of divergent monomers which share a completely conserved 19-bp long box. Conserved αCenH3-associated DNA is also confirmed in the related mitotic Meloidogyne incognita group species suggesting preservation of both centromere protein and DNA constituents. We hypothesize that the absence of centromere drive in mitosis might allow for CenH3 and its associated DNA to achieve an equilibrium in which they can persist for long periods of time.

Genome assembly and annotation of Meloidogyne enterolobii, an emerging parthenogenetic root-knot nematode

Root-knot nematodes (genus Meloidogyne) are plant parasites causing huge economic loss in the agricultural industry and affecting severely numerous developing countries. Control methods against these plant pests are sparse, the preferred one being the deployment of plant cultivars bearing resistance genes against Meloidogyne species. However, M. enterolobii is not controlled by the resistance genes deployed in the crop plants cultivated in Europe. The recent identification of this species in Europe is thus a major concern. Here, we sequenced the genome of M. enterolobii using short and long-read technologies. The genome assembly spans 240 Mbp with contig N50 size of 143 kbp, enabling high-quality annotations of 59,773 coding genes, 4,068 non-coding genes, and 10,944 transposable elements (spanning 8.7% of the genome). We validated the genome size by flow cytometry and the structure, quality and completeness by bioinformatics metrics. This ensemble of resources will fuel future projects aiming at pinpointing the genome singularities, the origin, diversity, and adaptive potential of this emerging plant pest.

Measurement(s)	genome • sequence_assembly • sequence feature annotation
Technology Type(s)	DNA sequencing assay • sequence assembly process • sequence annotation
Sample Characteristic - Organism	Meloidogyne enterolobii

Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.12410363

Genome structure and content of the rice root-knot nematode (Meloidogyne graminicola)

Discovered in the 1960s, Meloidogyne graminicola is a root-knot nematode species considered as a major threat to rice production. Yet, its origin, genomic structure, and intraspecific diversity are poorly understood. So far, such studies have been limited by the unavailability of a sufficiently complete and well-assembled genome. In this study, using a combination of Oxford Nanopore Technologies and Illumina sequencing data, we generated a highly contiguous reference genome (283 scaffolds with an N50 length of 294 kb, totaling 41.5 Mb). The completeness scores of our assembly are among the highest currently published for Meloidogyne genomes. We predicted 10,284 protein-coding genes spanning 75.5% of the genome. Among them, 67 are identified as possibly originating from horizontal gene transfers (mostly from bacteria), which supposedly contribute to nematode infection, nutrient processing, and plant defense manipulation. Besides, we detected 575 canonical transposable elements (TEs) belonging to seven orders and spanning 2.61% of the genome. These TEs might promote genomic plasticity putatively related to the evolution of M. graminicola parasitism. This high-quality genome assembly constitutes a major improvement regarding previously available versions and represents a valuable molecular resource for future phylogenomic studies of Meloidogyne species. In particular, this will foster comparative genomic studies to trace back the evolutionary history of M. graminicola and its closest relatives.

Profiling microRNAs through development of the parasitic nematode Haemonchus identifies nematode-specific miRNAs that suppress larval development

Parasitic nematodes transition between dramatically different free-living and parasitic stages, with correctly timed development and migration crucial to successful completion of their lifecycle. However little is known of the mechanisms controlling these transitions. microRNAs (miRNAs) negatively regulate gene expression post-transcriptionally and regulate development of diverse organisms. Here we used microarrays to determine the expression profile of miRNAs through development and in gut tissue of the pathogenic nematode Haemonchus contortus. Two miRNAs, mir-228 and mir-235, were enriched in infective L3 larvae, an arrested stage analogous to Caenorhabditis elegans dauer larvae. We hypothesized that these miRNAs may suppress development and maintain arrest. Consistent with this, inhibitors of these miRNAs promoted H. contortus development from L3 to L4 stage, while genetic deletion of C. elegans homologous miRNAs reduced dauer arrest. Epistasis studies with C. elegans daf-2 mutants showed that mir-228 and mir-235 synergise with FOXO transcription factor DAF-16 in the insulin signaling pathway. Target prediction suggests that these miRNAs suppress metabolic and transcription factor activity required for development. Our results provide novel insight into the expression and functions of specific miRNAs in regulating nematode development and identify miRNAs and their target genes as potential therapeutic targets to limit parasite survival within the host.

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

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4686-x) contains supplementary material, which is available to authorized users.

Hybridization and polyploidy enable genomic plasticity without sex in the most devastating plant-parasitic nematodes

Root-knot nematodes (genus Meloidogyne) exhibit a diversity of reproductive modes ranging from obligatory sexual to fully asexual reproduction. Intriguingly, the most widespread and devastating species to global agriculture are those that reproduce asexually, without meiosis. To disentangle this surprising parasitic success despite the absence of sex and genetic exchanges, we have sequenced and assembled the genomes of three obligatory ameiotic and asexual Meloidogyne. We have compared them to those of relatives able to perform meiosis and sexual reproduction. We show that the genomes of ameiotic asexual Meloidogyne are large, polyploid and made of duplicated regions with a high within-species average nucleotide divergence of ~8%. Phylogenomic analysis of the genes present in these duplicated regions suggests that they originated from multiple hybridization events and are thus homoeologs. We found that up to 22% of homoeologous gene pairs were under positive selection and these genes covered a wide spectrum of predicted functional categories. To biologically assess functional divergence, we compared expression patterns of homoeologous gene pairs across developmental life stages using an RNAseq approach in the most economically important asexually-reproducing nematode. We showed that >60% of homoeologous gene pairs display diverged expression patterns. These results suggest a substantial functional impact of the genome structure. Contrasting with high within-species nuclear genome divergence, mitochondrial genome divergence between the three ameiotic asexuals was very low, signifying that these putative hybrids share a recent common maternal ancestor. Transposable elements (TE) cover a ~1.7 times higher proportion of the genomes of the ameiotic asexual Meloidogyne compared to the sexual relative and might also participate in their plasticity. The intriguing parasitic success of asexually-reproducing Meloidogyne species could be partly explained by their TE-rich composite genomes, resulting from allopolyploidization events, and promoting plasticity and functional divergence between gene copies in the absence of sex and meiosis.

Mitochondrial Haplotype-based Identification of Root-knot Nematodes (Meloidogyne spp.) on Cut Foliage Crops in Florida

Florida accounts for more than 75% of the national cut foliage production. Unfortunately, root-knot nematodes (RKN) (Meloidogyne spp.) are a serious problem on these crops, rendering many farms unproductive. Currently, information on the Meloidogyne spp. occurring on most commonly cultivated cut foliage crops in Florida, and tools for their rapid identification are lacking. The objectives of this study were to (i) identify specific RKN infecting common ornamental cut foliage crops in Florida and (ii) evaluate the feasibility of using the mtDNA haplotype as a molecular diagnostic tool for rapid identification of large samples of RKN. A total of 200 Meloidogyne females were collected from cut foliage plant roots. Meloidogyne spp. were identified by PCR and RFLP of mitochondrial DNA. PCR and RFLP of mitochondrial DNA were effective in discriminating the Meloidogyne spp. present. Meloidogyne incognita is the most dominant RKN on cut foliage crops in Florida and must be a high target for making management decisions. Other Meloidogyne spp. identified include M. javanica, M. hapla, Meloidogyne sp. 1, and Meloidogyne sp. 2. The results for this study demonstrate the usefulness of the mtDNA haplotype-based designation as a valuable molecular tool for identification of Meloidogyne spp.

Epigenetic Inheritance and Its Role in Evolutionary Biology: Re-Evaluation and New Perspectives

Epigenetics increasingly occupies a pivotal position in our understanding of inheritance, natural selection and, perhaps, even evolution. A survey of the PubMed database, however, reveals that the great majority (>93%) of epigenetic papers have an intra-, rather than an inter-generational focus, primarily on mechanisms and disease. Approximately ~1% of epigenetic papers even mention the nexus of epigenetics, natural selection and evolution. Yet, when environments are dynamic (e.g., climate change effects), there may be an “epigenetic advantage” to phenotypic switching by epigenetic inheritance, rather than by gene mutation. An epigenetically-inherited trait can arise simultaneously in many individuals, as opposed to a single individual with a gene mutation. Moreover, a transient epigenetically-modified phenotype can be quickly “sunsetted”, with individuals reverting to the original phenotype. Thus, epigenetic phenotype switching is dynamic and temporary and can help bridge periods of environmental stress. Epigenetic inheritance likely contributes to evolution both directly and indirectly. While there is as yet incomplete evidence of direct permanent incorporation of a complex epigenetic phenotype into the genome, doubtlessly, the presence of epigenetic markers and the phenotypes they create (which may sort quite separately from the genotype within a population) will influence natural selection and, so, drive the collective genotype of a population.

Distribution of CpG Motifs in Upstream Gene Domains in a Reef Coral and Sea Anemone: Implications for Epigenetics in Cnidarians

Coral reefs are under assault from stressors including global warming, ocean acidification, and urbanization. Knowing how these factors impact the future fate of reefs requires delineating stress responses across ecological, organismal and cellular scales. Recent advances in coral reef biology have integrated molecular processes with ecological fitness and have identified putative suites of temperature acclimation genes in a Scleractinian coral Acropora hyacinthus. We wondered what unique characteristics of these genes determined their coordinate expression in response to temperature acclimation, and whether or not other corals and cnidarians would likewise possess these features. Here, we focus on cytosine methylation as an epigenetic DNA modification that is responsive to environmental stressors. We identify common conserved patterns of cytosine-guanosine dinucleotide (CpG) motif frequencies in upstream promoter domains of different functional gene groups in two cnidarian genomes: a coral (Acropora digitifera) and an anemone (Nematostella vectensis). Our analyses show that CpG motif frequencies are prominent in the promoter domains of functional genes associated with environmental adaptation, particularly those identified in A. hyacinthus. Densities of CpG sites in upstream promoter domains near the transcriptional start site (TSS) are 1.38x higher than genomic background levels upstream of -2000 bp from the TSS. The increase in CpG usage suggests selection to allow for DNA methylation events to occur more frequently within 1 kb of the TSS. In addition, observed shifts in CpG densities among functional groups of genes suggests a potential role for epigenetic DNA methylation within promoter domains to impact functional gene expression responses in A. digitifera and N. vectensis. Identifying promoter epigenetic sequence motifs among genes within specific functional groups establishes an approach to describe integrated cellular responses to environmental stress in reef corals and potential roles of epigenetics on survival and fitness in the face of global climate change.

Mitochondrial coding genome analysis of tropical root-knot nematodes (Meloidogyne) supports haplotype based diagnostics and reveals evidence of recent reticulate evolution

The polyphagous parthenogenetic root-knot nematodes of the genus Meloidogyne are considered to be the most significant nematode pest in sub-tropical and tropical agriculture. Despite the crucial need for correct diagnosis, identification of these pathogens remains problematic. The traditionally used diagnostic strategies, including morphometrics, host-range tests, biochemical and molecular techniques, now appear to be unreliable due to the recently-suggested hybrid origin of root-knot nematodes. In order to determine a suitable barcode region for these pathogens nine quickly-evolving mitochondrial coding genes were screened. Resulting haplotype networks revealed closely related lineages indicating a recent speciation, an anthropogenic-aided distribution through agricultural practices, and evidence for reticulate evolution within M. arenaria. Nonetheless, nucleotide polymorphisms harbor enough variation to distinguish these closely-related lineages. Furthermore, completeness of lineage sorting was verified by screening 80 populations from widespread geographical origins and variable hosts. Importantly, our results indicate that mitochondrial haplotypes are strongly linked and consistent with traditional esterase isozyme patterns, suggesting that different parthenogenetic lineages can be reliably identified using mitochondrial haplotypes. The study indicates that the barcode region Nad5 can reliably identify the major lineages of tropical root-knot nematodes.

Host genetic resistance to root-knot nematodes, Meloidogyne spp., in Solanaceae: from genes to the field

Root-knot nematodes (RKNs) heavily damage most solanaceous crops worldwide. Fortunately, major resistance genes are available in a number of plant species, and their use provides a safe and economically relevant strategy for RKN control. From a structural point of view, these genes often harbour NBS-LRR motifs (i.e. a nucleotide binding site and a leucine rich repeat region near the carboxy terminus) and are organised in syntenic clusters in solanaceous genomes. Their introgression from wild to cultivated plants remains a challenge for breeders, although facilitated by marker-assisted selection. As shown with other pathosystems, the genetic background into which the resistance genes are introgressed is of prime importance to both the expression of the resistance and its durability, as exemplified by the recent discovery of quantitative trait loci conferring quantitative resistance to RKNs in pepper. The deployment of resistance genes at a large scale may result in the emergence and spread of virulent nematode populations able to overcome them, as already reported in tomato and pepper. Therefore, careful management of the resistance genes available in solanaceous crops is crucial to avoid significant reduction in the duration of RKN genetic control in the field. From that perspective, only rational management combining breeding and cultivation practices will allow the design and implementation of innovative, sustainable crop production systems that protect the resistance genes and maintain their durability.

COBRA-Seq: Sensitive and Quantitative Methylome Profiling

Combined Bisulfite Restriction Analysis (COBRA) quantifies DNA methylation at a specific locus. It does so via digestion of PCR amplicons produced from bisulfite-treated DNA, using a restriction enzyme that contains a cytosine within its recognition sequence, such as TaqI. Here, we introduce COBRA-seq, a genome wide reduced methylome method that requires minimal DNA input (0.1-1.0 mg) and can either use PCR or linear amplification to amplify the sequencing library. Variants of COBRA-seq can be used to explore CpG-depleted as well as CpG-rich regions in vertebrate DNA. The choice of enzyme influences enrichment for specific genomic features, such as CpG-rich promoters and CpG islands, or enrichment for less CpG dense regions such as enhancers. COBRA-seq coupled with linear amplification has the additional advantage of reduced PCR bias by producing full length fragments at high abundance. Unlike other reduced representative methylome methods, COBRA-seq has great flexibility in the choice of enzyme and can be multiplexed and tuned, to reduce sequencing costs and to interrogate different numbers of sites. Moreover, COBRA-seq is applicable to non-model organisms without the reference genome and compatible with the investigation of non-CpG methylation by using restriction enzymes containing CpA, CpT, and CpC in their recognition site.