The quality of the raw RNA-seq data was assessed with FastQC v0.11.5. Paired-end reads were subsequently trimmed and de novo assembled into transcripts with Trinity v2.4.0 (34, 35). Post-assembly quality control and taxonomic partitioning were done with BlobTools (36, 37) and phyloFlash (38). We did not find any substantial contamination. Transcript abundance was estimated with the RSEM (RNA-Seq by Expectation Maximization) method using the normalized expression values TPM (transcripts per kilobase million) and FPKM (fragments per kilobase million). Candidate coding regions were identified with TransDecoder v5.1.0. Completeness of the transcriptome was assessed with BUSCO v3.0.2b (39, 40) and compared to other available flatworm transcriptomes (fig. S13) (41). The transcriptome was annotated using the Trinotate v3.1.1 annotation pipeline (42), which includes blastx and blastp homology searches to the BLAST+/SwissProt databases (data files S1 and S2), protein domain identification with HMMER/Pfam, signal peptide and transmembrane domain prediction with signalP/tmHMM, and functional annotation against eggNOG/GO/KEGG databases.

Note: The content above has been extracted from a research article, so it may not display correctly.

Please log in to submit your questions online.
Your question will be posted on the Bio-101 website. We will send your questions to the authors of this protocol and Bio-protocol community members who are experienced with this method. you will be informed using the email address associated with your Bio-protocol account.

We use cookies on this site to enhance your user experience. By using our website, you are agreeing to allow the storage of cookies on your computer.