A gene model was built based on the WS260 annotation. Tag counts for each gene were extracted from STAR aligned BAM files, anddifferential gene expression between N2 and mutant backgrounds was tested using DESeq2. A false discovery rate(FDR) < 0.01 and LFC > 0.5849 was used to define genes as upregulated, and FDR < 0.01 and LFC < -1 was used to define genes asdownregulated.