Supplemental files for study investigating the impact of Vibrio anguillarum infection on diploid and triploid Chinook salmon (Oncorhynchus tshawytscha) was investigated to identify if there was any significant immune regulation by microRNAs (miRNA). Small RNAs from hindgut (hg), head kidney (hk), and spleen (sp) were sequenced to determine if miRNA transcript abundance was altered due to ploidy and infection in nine-month old full-sibling diploids and triploids. MiARma-seq and bowtie1 were utilized to align reads to the Chinook salmon genome, and then read abundance associated with miRNA genes quantified using FeatureCounts. Read counts were then analyzed using DESeq2 to identify differentially expressed miRNA between ploidies as well as after infection in all three tissues. Finally, targets of differentially expressed miRNA were predicted using miRanda and pathway impacts predicted using KEGG analysis.
Files included: run_log.log - miARma-seq logfiles for hindgut (hg), head kidney (hk), and spleen (sp) stats.log - miARma-seq stats information for hindgut (hg), head kidney (hk), and spleen (sp) miarma.cadonic.known.tshawyt.sample.ini - initiation file for miARma-seq which includes all parameters that were utilized to analyze raw reads deseq2_cadonic_sample.R - sample of R code for DEseq2 used to identify differentially expressed miRNA genes Sasa_Ensembl_3UTR_clean.fasta - File of 3' untranslated region (UTR) for annotated genes in Ensembl biomart (Ensembl release 105) that was analyzed with miRanda (version v3.3a) to predict miRNA-mRNA relationships. Sequences with unavailable UTRs were cleaned by replacing empty spaces with N using "sed -e '/^[^>]/s/[^ATGC]/N/g' Sasa_Ensembl_3UTR.txt > Sasa_Ensembl_3UTR_clean.fasta". cadonic_2023_Bioinformatics_protocol_extended.docx - Extended explanation for the bioinformatic steps and justification for the workflow used in the study.
Files included: run_log.log - miARma-seq logfiles for hindgut (hg), head kidney (hk), and spleen (sp) stats.log - miARma-seq stats information for hindgut (hg), head kidney (hk), and spleen (sp) miarma.cadonic.known.tshawyt.sample.ini - initiation file for miARma-seq which includes all parameters that were utilized to analyze raw reads deseq2_cadonic_sample.R - sample of R code for DEseq2 used to identify differentially expressed miRNA genes Sasa_Ensembl_3UTR_clean.fasta - File of 3' untranslated region (UTR) for annotated genes in Ensembl biomart (Ensembl release 105) that was analyzed with miRanda (version v3.3a) to predict miRNA-mRNA relationships. Sequences with unavailable UTRs were cleaned by replacing empty spaces with N using "sed -e '/^[^>]/s/[^ATGC]/N/g' Sasa_Ensembl_3UTR.txt > Sasa_Ensembl_3UTR_clean.fasta". cadonic_2023_Bioinformatics_protocol_extended.docx - Extended explanation for the bioinformatic steps and justification for the workflow used in the study.