Dauer juvenile recovery transcriptome of two contrasting EMS mutants of the entomopathogenic nematode Heterorhabditis bacteriophora

The entomopathogenic nematode Heterorhabditis bacteriophora, symbiotically associated with enterobacteria of the genus Photorhabdus, is a biological control agent against many insect pests. Dauer Juveniles (DJ) of this nematode are produced in industrial-scale bioreactors up to 100 m3 in liquid culture processes lasting approximately 11 days. A high DJ yield (> 200,000 DJ·mL−1) determines the success of the process. To start the mass production, a DJ inoculum proceeding from a previous monoxenic culture is added to pre-cultured (24 h) Photorhabdus bacteria. Within minutes after contact with the bacteria, DJ are expected to perceive signals that trigger their further development (DJ recovery) to reproductive hermaphrodites. A rapid, synchronized, and high DJ recovery is a key factor for an efficient culture process. In case of low percentage of DJ recovery, the final DJ yield is drastically reduced, and the amount of non-desired stages (males and non-fertilized females) hinders the DJ harvest. In a preliminary work, a huge DJ recovery phenotypic variability in H. bacteriophora ethyl methanesulphonate (EMS) mutants was determined. In the present study, two EMS-mutant lines (M31 and M88) with high and low recovery phenotypes were analyzed concerning their differences in gene expression during the first hours of contact with Photorhabdus supernatant containing food signals triggering recovery. A snapshot (RNA-seq analysis) of their transcriptome was captured at 0.5, 1, 3 and 6 h after exposure. Transcripts (3060) with significant regulation changes were identified in the two lines. To analyze the RNA-seq data over time, we (1) divided the expression profiles into clusters of similar regulation, (2) identified over and under-represented gene ontology categories for each cluster, (3) identified Caenorhabditis elegans homologous genes with recovery-related function, and (4) combined the information with available single nucleotide polymorphism (SNP) data. We observed that the expression dynamics of the contrasting mutants (M31 and M88) differ the most within the first 3 h after Photorhabdus supernatant exposure, and during this time, genes related to changes in the DJ cuticle and molting are more active in the high-recovery line (M31). Comparing the gene expression of DJ exposed to the insect food signal in the haemolymph, genes related to host immunosuppressive factors were not found in DJ upon bacterial supernatant exposure. No link between the position of SNPs associated with high recovery and changes in gene expression was determined for genes with high differential expression. Concerning specific transcripts, nine H. bacteriophora gene models with differential expression are provided as candidate genes for further studies.