Before transcriptome assembly, adapters and low-quality reads (minimum read length, 35 bp) were removed using Cutadapt ( Because the available P. parnellii genome (4) is highly fragmented, we used the Trinity (version 2.0.6) de novo assembly package ( to assemble the P. parnellii transcriptome. Highly similar transcripts were collapsed using CD-Hit ( Annotation of the transcriptome was carried out using Blast+ against human complementary DNA sequences from Ensembl (version 86). We kept the best Blast hit with the lowest E value and highest percent identity, which associates each transcript to a human Ensembl gene.