Differential expression analysis can be a powerful method to investigate expressed differences between closely related species. Our ambition is to highlight differentially expressed nuclear genes to explain the hybrid incompatibilities among the Amazonian, Centro-American and Orinocan semispecies of Drosophila paulistorum. RNA sequencing (RNA-seq) establishes the foundation of the study where we first evaluate the influence of two distinct alignment references. We discover the benefits of concatenating a de novo assembly instead of using the genome reference of a close relative. The bioinformatic pipeline handles the interesting inclusion of D. melanogaster and D. willistoni, where their contribution assists in the search for previously studied speciation genes. Among the down- and upregulated subsets we can see a diverse mix of general biological processes such as regulatory functions and transcriptional factors. In the end we uncover potential indications to why the Amazonian seems to be the least compatible semispecie to produce hybrids. This study provides a competitive working frame for comparative RNA-seq studies between closely related species.