Transcriptomic Analysis of Early Stages of Intestinal Regeneration in Holothuria Glaberrima David J
Total Page:16
File Type:pdf, Size:1020Kb
Transcriptomic Analysis of Early Stages of Intestinal Regeneration in Holothuria glaberrima David J. Quispe-Parra1, Joshua G. Medina-Feliciano1, Sebastián Cruz-González1, Humberto Ortiz-Zuazaga2, José E. García-Arrarás1* 1University of Puerto Rico, Biology Department, San Juan, 00925, Puerto Rico. 2University of Puerto Rico, Department of Computer Sciences, San Juan, 00925, Puerto Rico. Table S1. Results of transcriptome assessment with BUSCO Parameter BUSCO result Core genes queried 978 Complete core genes detected 99.1% Complete single copy core genes 27.2% Complete duplicated core genes 71.9% Fragmented core genes detected 0.4% Missing core genes 0.5% Table S2. Transcriptome length statistics and composition assessments with gVolante Parameter Result Number of sequences 491 436 Total length (nt) 408 930 895 Longest sequence (nt) 34 610 Shortest sequence (nt) 200 Mean sequence length (nt) 832 N50 sequence length (nt) 1 691 Table S3. RNA-seq data read statistic values Accession Quantity of reads Sample Mapped reads (SRA) Before Filtering After Filtering SRR12564573 NormalA 89 424 588.00 88 112 720.00 94.18% SRR12564572 NormalB 74 105 862.00 72 983 434.00 93.79% SRR12564570 NormalC 60 626 094.00 59 767 896.00 91.55% SRR12564564 Day1A 32 499 910.00 31 721 976.00 90.67% SRR12564563 Day1B 35 164 514.00 34 237 960.00 91.37% SRR12564571 Day1C 43 519 280.00 42 179 220.00 91.53% SRR12564567 Day3A 64 094 536.00 62 995 792.00 90.26% SRR12564566 Day3B 71 654 462.00 70 158 078.00 91.07% SRR12564565 Day3C 71 259 560.00 70 418 332.00 90.75% SRR12564569 Day3D* 107 688 184.00 106 083 898.00 - SRR12564568 Day3E* 108 372 334.00 106 767 522.00 - *Samples used for the assembly but not for differential expression analysis Table S4. Genes and primers used for RT-qPCR Gene Primers F: AGGTCAAAATGTCCGGCCTCTGT VBP R: TCACGCTGGCGGTGTCAAAACT F: ACGGCGAAGGTAGCATCCGTTT FoxA R: ATGGCAATCCAGCAAGCGCCAA F: CCTGCAATAATCCGGTGAGT Wnt6 R: CGCGAATGTAAATGTCATGG F: CGCTCTCCATTTCAGAGGAC Sox4 R: GTCTTTCTGCCTCCTCAACG F: TGCCTGAGCTTTCTGTTCCT Tap26 R: CCCTCAAAGAATGGTGGAAA F: CAATGGTTGTTGCTGGAGTCTTT R: NADH CGCAGAAGTAGCCGCGAATAT F: CGTTGAAATACCGTCCATCC Tenascin-N R: AACGAAGGATACGCTGGAGA Table S5. Mean Ct values obtained from RT-qPCR Gene Stage Ct Mean Ct SD Normal 28.636 1.268 VBP Day 1 28.263 0.640 Day 3 27.958 0.738 Normal 23.611 0.371 FoxA Day 1 27.643 3.124 Day 3 27.824 3.201 Wnt6 Normal 28.607 1.399 Day 1 25.253 1.965 Day 3 27.842 3.703 Normal 26.038 0.784 Sox4 Day 1 25.424 1.143 Day 3 29.242 1.014 Normal 27.634 1.319 Tap26 Day 1 26.727 0.483 Day 3 29.512 1.311 Normal 18.843 0.271 NADH Day 1 20.829 0.355 Day 3 22.583 0.341 Normal 21.025 3.217 Tenascin Day 3 30.992 1.519 Table S6. Comparison of expressed transcripts in RNA-seq with previous studies RT- Northern Blot/In RNA-seq PCR/qPCR situ Gene Day 1 Day 3 Day 3 Days 2-3 (Log2FC) (Log2FC) Myc 2.72 ↑ 2.05 ↑ ↑ NA β-catenin 0.88 0.46NS -NS NA Melanotransferrin 3.55 ↑ 3.7 ↑ ↑ NA Serum amyloid A 2.39 ↑ 3.2 ↑ NA ↑ protein WNT9 2.69 ↑ 2.88 ↑ NA ↑ BMP 3.53 ↑ 3.31 ↑ NA ↑ survivin -0.90NS -0.31NS -NS - ↑: Upregulated; ↓: Downregulated; -: No significant change; NS: Not significant (Pvalue > 0.05) Table S7. Comparison of expressed transcripts in RNA-seq with microarray study Gene RNA-seq Microarray Day 1 Day 3 Day 3 (Log2FC) (Log2FC) Tensc-R 3.26↑ 2.03 NS ↑ Actin-1 6.12↑ 7.47↑ ↑ Actin-2 1.71↑ 1.56↑ ↑ Collagen alfa-1 -1.54↓ -0.67 NS ↑ Laminin alpha1 -3.67↓ -2.79↓ ↑ Stromelysin-3 -4.16↓ -2.65↓ ↑ (MMP-11) Hox9 1.03↑ -0.83 NS ↑ Krueppel like -1.09↓ 0.12 NS ↓ Gelsolin -0.79 NS -2.69↓ ↓ ↑: Upregulated; ↓: Downregulated; -: No significant change; NS: Not significant (Pvalue > 0.05) Table S8. Differentially expressed transcripts between day 3 vs 1 comparison Adjusted p- ID Description Log2FC value Transcript_223195 15-hydroxyprostaglandin dehydrogenase 2.13 3.58E-02 [NAD(+)] Transcript_298217 4-aminobutyrate aminotransferase, 2.98 2.81E-05 mitochondrial Transcript_334800 40S ribosomal protein S15 6.73 1.89E-04 Transcript_254071 actin, cytoskeletal 3B|actin, cytoskeletal 3 3.00 1.19E-02 Transcript_432551 actin, muscle 3.65 5.45E-03 Transcript_383802 adenylyltransferase and sulfurtransferase 3.21 3.38E-02 MOCS3-like Transcript_308851 alcohol dehydrogenase class-3 2.33 1.59E-02 Transcript_425811 calreticulin 2.45 9.78E-08 Transcript_298293 collagen alpha-1(XII) chain|tenascin-X 2.08 7.99E-05 Transcript_196327 creatine kinase, flagellar 5.16 4.69E-02 Transcript_380459 cryptochrome-2 2.03 3.75E-03 Transcript_226806 cytochrome P450 3A24 2.25 1.83E-02 Transcript_262552 cytochrome P450 3A9 4.03 3.31E-02 Transcript_175314 cytochrome P450 4V2 2.00 4.85E-02 Transcript_172290 cytosol aminopeptidase 2.21 1.97E-05 Transcript_184414 dehydrogenase/reductase SDR family 2.25 1.39E-03 member 7 Transcript_175770 deleted in malignant brain tumors 1 protein 2.65 2.63E-02 Transcript_461577 deoxyribodipyrimidine photo-lyase 2.68 2.66E-03 Transcript_432963 DNA replication licensing factor mcm5 2.57 4.81E-07 Transcript_359454 dynein light chain LC6, flagellar outer arm 2.35 6.25E-04 Transcript_347359 ectonucleotide 3.16 4.27E-03 pyrophosphatase/phosphodiesterase family member 7-like Transcript_273238 epidermal growth factor-like protein 7 2.49 1.23E-05 Transcript_276903 ERI1 exoribonuclease 2 2.86 3.80E-03 Transcript_365063 extracellular transglutaminase 2.62 5.99E-09 Transcript_350713 fibrillin-1 2.08 9.82E-03 Transcript_356914 ficolin-2-like 3.74 1.04E-02 Transcript_179574 ficolin-2-like 3.05 3.76E-02 Transcript_221519 formin-J 2.02 1.31E-04 Transcript_358088 fucose mutarotase 2.46 1.86E-08 Transcript_351559 G patch domain-containing protein 3 4.02 5.91E-03 Transcript_185057 GDH/6PGL endoplasmic bifunctional 2.07 2.55E-06 protein|hexose-6-phosphate dehydrogenase (glucose 1-dehydrogenase) Transcript_402303 GTP 3',8-cyclase, mitochondrial|cyclic 2.32 8.19E-03 pyranopterin monophosphate synthase, mitochondrial|molybdenum cofactor biosynthesis protein 1 Transcript_438891 H2.0-like homeobox protein|Homeobox 3.45 1.75E-11 domain-containing protein Transcript_251717 hairy/enhancer-of-split related with YRPW 3.17 2.39E-02 motif protein 1 Transcript_456276 heparan sulfate glucosamine 3-O- 3.31 1.15E-02 sulfotransferase 1 Transcript_438353 homeobox protein EMX1 3.86 8.19E-03 Transcript_290553 isoamyl acetate-hydrolyzing esterase 1 2.80 5.95E-05 homolog Transcript_409717 isochorismatase domain-containing protein 2.95 4.66E-03 2|isochorismatase domain-containing protein 2, mitochondrial Transcript_258900 KRP170 2.69 1.74E-02 Transcript_189831 L-gulonolactone oxidase 2.45 6.22E-03 Transcript_297492 laminin subunit alpha-2|laminin subunit 2.21 1.74E-05 alpha-1 Transcript_305145 LOW QUALITY PROTEIN: DNA 2.94 6.14E-07 replication licensing factor mcm2 Transcript_381204 LOW QUALITY PROTEIN: transcription 3.54 6.83E-09 factor Sox-10 Transcript_285150 LOW QUALITY PROTEIN: zygotic DNA 2.29 2.08E-07 replication licensing factor mcm6-B Transcript_294196 low-density lipoprotein receptor-related 2.14 2.64E-04 protein 1 Transcript_328691 MAM and LDL-receptor class A domain- 2.14 7.82E-04 containing protein 1 Transcript_257273 methylsterol monooxygenase 1 3.02 3.79E-03 Transcript_324570 microfibril-associated glycoprotein 2.48 3.32E-02 4|tenascin-N Transcript_263515 polycystic kidney disease protein 1-like 2 2.32 1.61E-02 Transcript_269460 probable D-lactate dehydrogenase, 3.55 2.01E-02 mitochondrial|lactate dehydrogenase D Transcript_236563 protein FAM166B-like 17.42 4.96E-04 Transcript_437843 protein lin-52 homolog 2.10 2.16E-03 Transcript_375637 protein Wnt-6 2.25 2.75E-03 Transcript_197666 putative aminopeptidase W07G4.4 2.06 2.59E-04 Transcript_311247 putative hydroxypyruvate isomerase 2.52 5.99E-04 Transcript_443748 pyridine nucleotide-disulfide 2.41 1.73E-02 oxidoreductase domain-containing protein 2|pyridine nucleotide-disulphide oxidoreductase domain 2 Transcript_297199 ryncolin-1-like 5.18 3.37E-02 Transcript_406209 sulfite oxidase 3.85 2.99E-07 Transcript_54760 testis-specific serine/threonine-protein 5.30 3.42E-02 kinase 4-like Transcript_251414 transcription factor Sp5 2.40 4.19E-03 Transcript_384474 transmembrane protein KIAA1109 2.08 1.48E-05 Transcript_266841 tuftelin 2.47 1.33E-02 Transcript_368226 ubiquitin 6.29 1.18E-03 Transcript_270219 ubiquitin-40S ribosomal protein S27a 6.47 4.00E-04 Transcript_437026 uncharacterized protein K02A2.6-like 3.43 2.72E-05 Transcript_190019 uncharacterized protein K02A2.6- 2.32 1.79E-02 like|Reverse transcriptase domain- containing protein Transcript_291904 uncharacterized protein 2.14 3.79E-04 LOC105438010|Glycoside hydrolase family 31 domain containing protein Transcript_283678 uncharacterized protein 4.75 9.63E-03 LOC583353|neurotrypsin Transcript_26623 uncharacterized protein LOC592324 7.59 2.53E-06 Transcript_442525 uncharacterized protein 2.01 2.72E-04 LOC753842|Protease inhibitor I35 (TIMP) domain containing protein Transcript_270528 uncharacterized protein 2.01 7.71E-05 LOC757055|glyoxylate/hydroxypyruvate reductase A HPR2 Transcript_263687 valacyclovir hydrolase-like 2.34 1.86E-03 Transcript_428635 zygotic DNA replication licensing factor 2.64 3.93E-04 mcm3 Transcript_167684 Transcript_294311 - -3.44 2.27E-02 Transcript_251847 - -2.32 7.18E-05 Transcript_254514 - -2.01 2.89E-02 Transcript_395263 14-3-3 family protein artA|14-3-3 protein 3 -5.90 2.83E-03 Transcript_285870 26S proteasome non-ATPase regulatory -3.03 1.24E-09 subunit 5 Transcript_42693 26S proteasome regulatory subunit 6A- -2.16 6.22E-06 B|26S protease regulatory subunit 6A- B|26S protease regulatory subunit 6A-like Transcript_48454 40S ribosomal protein