Gene expression changes with tumor disease and leech parasitism in the juvenile green sea turtle skin transcriptome

Emerging infectious diseases are a major threat to biodiversity in the 21st century. Fibropapillomatosis (FP) is an epithelial tumor disease that affects immature and adult marine turtles worldwide, particularly green turtles (Chelonia mydas). We know little about the host factors contributing to FP susceptibility, in part because transcriptomic studies that compare transcript expression in turtles with and without FP are lacking. Here, we performed RNA-Seq on healthy skin tissue from immature C. mydas in the Indian River Lagoon, Florida, USA, comparing turtles (1) with and without FP and (2) with and without leech parasites, a putative vector of FP. We assembled a de novo C. mydas skin transcriptome to identify transcripts with significant differential expression (DE) across FP and leech categories. Significant DE transcripts were found across FP and leech comparisons, including 10 of the same transcripts with DE across both comparisons. Leech-positive individuals significantly upregulated different immune and viral interaction transcripts than did leech-negative individuals, including viral interaction transcripts associated with herpesvirus interactions. This finding strengthens the role of marine leeches as mechanical vectors of Chelonid herpesvirus 5 (ChHV5) which has been implicated as a causative agent of FP. FP-positive turtles upregulated several tumor progression and suppression transcripts relative to FP-negative turtles, which had no significant DE tumor progression transcripts. FP-positive turtles also upregulated significantly more protein interaction transcripts than FP-negative turtles. DE transcripts across leech comparisons showed no functional enrichment, whereas DE transcripts across FP comparisons showed some GO terms were enriched in FP-positive and FP negative turtles. Notably, only FP-negative turtles were enriched for GO terms involved in acquired and inflammatory immune gene regulation. Overall, our DE transcripts included several candidate genes that may play important roles in C. mydas resistance to or recovery from FP, highlighting that transcriptomics provides a promising venue to understand this impactful disease. Continued investigation of C. mydas responses to FP and leech affliction is imperative for species persistence and the conservation of marine ecosystems worldwide due to the essential role of sea turtles in ecosystem function and stability.