Screening of intracellular parasites
To explore the presence of relevant intracellular parasites that could explain their speciation process (60), we collapsed unique reads from the butterfly-free sequences and removed from the dataset low complexity sequences using Prinseq (61). Afterwards, we applied kraken2 to assign reads against standard contaminant databases (bacteria, archaea, fungi, protozoa and viral)(62).