Fig. 4

A Comparative analyses of MelSat DNAs expression in M. incognita (isolates 28 and 27), B M. javanica and M. arenaria using whole transcriptome data from [28]. We mapped the MelSat consensus sequences on the available RNA-seq data using Bowtie2 [41]. The quantification step included read counts and scaled read counts, and scaling method applied was log2RPKM (Reads Per Kilobase Million) (Additional file 1: Table S3). Expression profile was shown as logarithmic transformation of RPKM. Gray boxes show the absence of particular satRNA in the transcriptome. Red lines indicate the level of GAPDH gene expression. C The proportion of MalSat transcripts and genome abundance for each MelSat was calculated as a ratio of MalSat transcripts in respect to the total transcripts and copy number of MelSat in respect to genome sizes and length of satDNA repeat unit for each MelSat. Light blue boxes indicate MelSat DNAs with positive ratio in all species where MelSat were found. The expression analyses were also performed for 4 M. incognita housekeeping genes (Additional file 1: Table S4)