Transcriptomic Profiling of Fibropapillomatosis in Green Sea Turtles (Chelonia mydas) From South Texas

Blood transcriptome analysis of common kestrel nestlings living in urban and non-urban environments

Urbanisation is one of the main anthropogenic forms of land cover affecting an ever-increasing number of wild animals and their habitats. Physiological plasticity represents an important process through which animals can adjust to the novel conditions of anthropogenic environments. Relying on the analysis of gene expression, it is possible to identify the molecular responses to the habitat conditions and infer possible environmental factors that affect the organismal physiology. We have quantified for the first time the blood transcriptome of common kestrel (Falco tinnunculus) nestlings living in urban sites and compared it to the transcriptome of kestrel nestlings inhabiting rural and natural environments. We found mild differences in the expression of genes among sites, indicating adaptability or acclimation of the birds to the urban habitat. We identified 58 differentially expressed genes between urban and natural kestrels, and 12 differentially expressed genes between urban and rural kestrels. The most striking differences among sites involved inflammatory-immunological, metabolic, apoptosis, DNA repair and development genes. In particular, we found that (i) urban kestrel nestlings had higher expression of genes linked to inflammation, repair of DNA damage, or apoptosis than natural kestrel nestlings, and (ii) natural and rural kestrel nestlings had higher expression of genes linked to the development and activation of immune cells, type I interferon response, or major histocompatibility complex than urban kestrel nestlings. Finally, the KEGG enrichment analysis identified the insulin signalling as the main pathway that differed between natural and urban kestrel nestlings. This is one of a limited number of studies on vertebrates that revealed habitat-associated differences in the transcriptome. It paves the way for further in-depth studies on the links between physiological variation and habitat structure at different spatial and temporal scales.

Immunity in Sea Turtles: Review of a Host-Pathogen Arms Race Millions of Years in the Running

The immune system of sea turtles is not completely understood. Sea turtles (as reptiles) bridge a unique evolutionary gap, being ectothermic vertebrates like fish and amphibians and amniotes like birds and mammals. Turtles are ectotherms; thus, their immune system is influenced by environmental conditions like temperature and season. We aim to review the turtle immune system and note what studies have investigated sea turtles and the effect of the environment on the immune response. Turtles rely heavily on the nonspecific innate response rather than the specific adaptive response. Turtles' innate immune effectors include antimicrobial peptides, complement, and nonspecific leukocytes. The antiviral defense is understudied in terms of the diversity of pathogen receptors and interferon function. Turtles also mount adaptive responses to pathogens. Lymphoid structures responsible for lymphocyte activation and maturation are either missing in reptiles or function is affected by season. Turtles are a marker of health for their marine environment, and their immune system is commonly dysregulated because of disease or contaminants. Fibropapillomatosis (FP) is a tumorous disease that afflicts sea turtles and is thought to be caused by a virus and an environmental factor. We aim, by exploring the current understanding of the immune system in turtles, to aid the investigation of environmental factors that contribute to the pathogenesis of this disease and provide options for immunotherapy.

Microbial composition of carapace, feces, and water column in captive juvenile green sea turtles with carapacial ulcers

Introduction

Green sea turtles are endangered marine reptiles. Carapacial ulcers will develop on juvenile green sea turtles during artificial rescue, seriously affecting their health and potentially leading to death.

Methods

To determine the pathogens causing ulcerative carapacial disease, we performed 16S and ITS high-throughput sequencing, and microbial diversity analysis on samples from carapacial ulcers, healthy carapaces, feces, and seawater of juvenile green sea turtles.

Results

Our analysis showed that changes in microbial diversity of green sea turtle feces and seawater were not significantly associated with ulcerative carapacial disease.

Discussion

Psychrobacter sp. is the dominant species in the carapacial ulcers of green sea turtles. The bacterium is present in both healthy turtles and seawater where carapacial ulcers did not occur and decreasing seawater temperatures are likely responsible for the infection of juvenile green turtles with Psychrobacter sp. This is the first study on carapacial ulcers in captive juvenile green sea turtles. Our research provides theoretical guidance for the prevention and control of carapacial ulcers in captive juvenile green sea turtles.

Adaptive evolution of major histocompatibility complex class I immune genes and disease associations in coastal juvenile sea turtles

Characterizing polymorphism at the major histocompatibility complex (MHC) genes is key to understanding the vertebrate immune response to disease. Despite being globally afflicted by the infectious tumour disease fibropapillomatosis (FP), immunogenetic variation in sea turtles is minimally explored. We sequenced the α 1 peptide-binding region of MHC class I genes (162 bp) from 268 juvenile green (Chelonia mydas) and 88 loggerhead (Caretta caretta) sea turtles in Florida, USA. We recovered extensive variation (116 alleles) and trans-species polymorphism. Supertyping analysis uncovered three functional MHC supertypes corresponding to the three well-supported clades in the phylogeny. We found significant evidence of positive selection at seven amino acid sites in the class I exon. Random forest modelling and risk ratio analysis of Ch. mydas alleles uncovered one allele weakly associated with smooth FP tumour texture, which may be associated with disease outcome. Our study represents the first characterization of MHC class I diversity in Ch. mydas and the largest sample of sea turtles used to date in any study of adaptive genetic variation, revealing tremendous genetic variation and high adaptive potential to viral pathogen threats. The novel associations we identified between MHC diversity and FP outcomes in sea turtles further highlight the importance of evaluating genetic predictors of disease, including MHC and other functional markers.

De novo assembly and functional annotation of blood transcriptome of loggerhead turtle, and in silico characterization of peroxiredoxins and thioredoxins

The aim of this study was to generate and analyze the atlas of the loggerhead turtle blood transcriptome by RNA-seq, as well as identify and characterize thioredoxin (Tnxs) and peroxiredoxin (Prdxs) antioxidant enzymes of the greatest interest in the control of peroxide levels and other biological functions. The transcriptome of loggerhead turtle was sequenced using the Illumina Hiseq 2000 platform and de novo assembly was performed using the Trinity pipeline. The assembly comprised 515,597 contigs with an N50 of 2,631 bp. Contigs were analyzed with CD-Hit obtaining 374,545 unigenes, of which 165,676 had ORFs encoding putative proteins longer than 100 amino acids. A total of 52,147 (31.5%) of these transcripts had significant homology matches in at least one of the five databases used. From the enrichment of GO terms, 180 proteins with antioxidant activity were identified, among these 28 Prdxs and 50 putative Tnxs. The putative proteins of loggerhead turtles encoded by the genes Prdx1, Prdx3, Prdx5, Prdx6, Txn and Txnip were predicted and characterized in silico. When comparing Prdxs and Txns of loggerhead turtle with homologous human proteins, they showed 18 (9%), 52 (18%) 94 (43%), 36 (16%), 35 (33%) and 74 (19%) amino acid mutations respectively. However, they showed high conservation in active sites and structural motifs (98%), with few specific modifications. Of these, Prdx1, Prdx3, Prdx5, Prdx6, Txn and Txnip presented 0, 25, 18, three, six and two deleterious changes. This study provides a high quality blood transcriptome and functional annotation of loggerhead sea turtles.

Fibropapillomatosis and Chelonid Alphaherpesvirus 5 Infection in Kemp's Ridley Sea Turtles (Lepidochelys kempii)

Simple Summary

The Kemp’s ridley sea turtle is an endangered species that is susceptible to a tumor disease called fibropapillomatosis (FP) and its associated virus, chelonid alphaherpesvirus 5 (ChHV5). The goal of our study was to describe FP in Kemp’s ridley turtles, including estimated disease prevalence and pathologyg, and case demographics and outcomes, to better understand the risk posed by FP to Kemp’s ridley population recovery. During 2006–2020, we identified 22 cases of Kemp’s ridley turtles with FP, including 12 adult turtles, a reproductively valuable age class. Molecular diagnostics were used to identify ChHV5 DNA in blood (7.8%) and tumor (91.7%) samples collected from free-ranging Kemp’s ridley turtles. Genomic sequencing was conducted to identify ChHV5 variants in tumor samples collected from Kemp’s ridley turtles with FP. Along with case data, phylogenetic analysis of resultant sequences suggests increasing, spatiotemporal spread of ChHV5 infections and FP among Kemp’s ridley turtles in coastal areas, including the Gulf of Mexico and the southwestern Atlantic Ocean, where they share habitat with green sea turtles (in which FP is enzootic). This is concerning because FP has an uncertain pathogenesis, is potentially related to anthropogenic environmental degradation, and can cause suffering and/or death in severely afflicted turtles.

Abstract

Fibropapillomatosis (FP), a debilitating, infectious neoplastic disease, is rarely reported in endangered Kemp’s ridley sea turtles (Lepidochelys kempii). With this study, we describe FP and the associated chelonid alphaherpesvirus 5 (ChHV5) in Kemp’s ridley turtles encountered in the United States during 2006–2020. Analysis of 22 case reports of Kemp’s ridley turtles with FP revealed that while the disease was mild in most cases, 54.5% were adult turtles, a reproductively valuable age class whose survival is a priority for population recovery. Of 51 blood samples from tumor-free turtles and 12 tumor samples from turtles with FP, 7.8% and 91.7%, respectively, tested positive for ChHV5 DNA via quantitative polymerase chain reaction (qPCR). Viral genome shotgun sequencing and phylogenetic analysis of six tumor samples show that ChHV5 sequences in Kemp’s ridley turtles encountered in the Gulf of Mexico and northwestern Atlantic cluster with ChHV5 sequences identified in green (Chelonia mydas) and loggerhead (Caretta caretta) sea turtles from Hawaii, the southwestern Atlantic Ocean, and the Caribbean. Results suggest an interspecific, spatiotemporal spread of FP among Kemp’s ridley turtles in regions where the disease is enzootic. Although FP is currently uncommon in this species, it remains a health concern due to its uncertain pathogenesis and potential relationship with habitat degradation.

Evolutionary Comparisons of Chelonid Alphaherpesvirus 5 (ChHV5) Genomes from Fibropapillomatosis-Afflicted Green (Chelonia mydas), Olive Ridley (Lepidochelys olivacea) and Kemp's Ridley (Lepidochelys kempii) Sea Turtles

The spreading global sea turtle fibropapillomatosis (FP) epizootic is threatening some of Earth's ancient reptiles, adding to the plethora of threats faced by these keystone species. Understanding this neoplastic disease and its likely aetiological pathogen, chelonid alphaherpesvirus 5 (ChHV5), is crucial to understand how the disease impacts sea turtle populations and species and the future trajectory of disease incidence. We generated 20 ChHV5 genomes, from three sea turtle species, to better understand the viral variant diversity and gene evolution of this oncogenic virus. We revealed previously underappreciated genetic diversity within this virus (with an average of 2035 single nucleotide polymorphisms (SNPs), 1.54% of the ChHV5 genome) and identified genes under the strongest evolutionary pressure. Furthermore, we investigated the phylogeny of ChHV5 at both genome and gene level, confirming the propensity of the virus to be interspecific, with related variants able to infect multiple sea turtle species. Finally, we revealed unexpected intra-host diversity, with up to 0.15% of the viral genome varying between ChHV5 genomes isolated from different tumours concurrently arising within the same individual. These findings offer important insights into ChHV5 biology and provide genomic resources for this oncogenic virus.

Environmental DNA monitoring of oncogenic viral shedding and genomic profiling of sea turtle fibropapillomatosis reveals unusual viral dynamics

Pathogen-induced cancers account for 15% of human tumors and are a growing concern for endangered wildlife. Fibropapillomatosis is an expanding virally and environmentally co-induced sea turtle tumor epizootic. Chelonid herpesvirus 5 (ChHV5) is implicated as a causative virus, but its transmission method and specific role in oncogenesis and progression is unclear. We applied environmental (e)DNA-based viral monitoring to assess viral shedding as a direct means of transmission, and the relationship between tumor burden, surgical resection and ChHV5 shedding. To elucidate the abundance and transcriptional status of ChHV5 across early, established, regrowth and internal tumors we conducted genomics and transcriptomics. We determined that ChHV5 is shed into the water column, representing a likely transmission route, and revealed novel temporal shedding dynamics and tumor burden correlations. ChHV5 was more abundant in the water column than in marine leeches. We also revealed that ChHV5 is latent in fibropapillomatosis, including early stage, regrowth and internal tumors; higher viral transcription is not indicative of poor patient outcome, and high ChHV5 loads predominantly arise from latent virus. These results expand our knowledge of the cellular and shedding dynamics of ChHV5 and can provide insights into temporal transmission dynamics and viral oncogenesis not readily investigable in tumors of terrestrial species.

Insights on Immune Function in Free-Ranging Green Sea Turtles (Chelonia mydas) with and without Fibropapillomatosis

Chelonid alphaherpesviruses 5 and 6 (ChHV5 and ChHV6) are viruses that affect wild sea turtle populations. ChHV5 is associated with the neoplastic disease fibropapillomatosis (FP), which affects green turtles (Chelonia mydas) in panzootic proportions. ChHV6 infection is associated with lung-eye-trachea disease (LETD), which has only been observed in maricultured sea turtles, although antibodies to ChHV6 have been detected in free-ranging turtles. To better understand herpesvirus prevalence and host immunity in various green turtle foraging aggregations in Florida, USA, our objectives were to compare measures of innate and adaptive immune function in relation to (1) FP tumor presence and severity, and (2) ChHV5 and ChHV6 infection status. Free-ranging, juvenile green turtles (N = 45) were captured and examined for external FP tumors in Florida's Big Bend, Indian River Lagoon, and Lake Worth Lagoon. Blood samples were collected upon capture and analyzed for ChHV5 and ChHV6 DNA, antibodies to ChHV5 and ChHV6, in vitro lymphocyte proliferation using a T-cell mitogen (concanavalin A), and natural killer cell activity. Despite an overall high FP prevalence (56%), ChHV5 DNA was only observed in one individual, whereas 20% of turtles tested positive for antibodies to ChHV5. ChHV6 DNA was not observed in any animals and only one turtle tested positive for ChHV6 antibodies. T-cell proliferation was not significantly related to FP presence, tumor burden, or ChHV5 seroprevalence; however, lymphocyte proliferation in response to concanavalin A was decreased in turtles with severe FP (N = 3). Lastly, green turtles with FP (N = 9) had significantly lower natural killer cell activity compared to FP-free turtles (N = 5). These results increase our understanding of immune system effects related to FP and provide evidence that immunosuppression occurs after the onset of FP disease.

The Concurrent Detection of Chelonid Alphaherpesvirus 5 and Chelonia mydas Papillomavirus 1 in Tumoured and Non-Tumoured Green Turtles

Simple Summary

Characterised by benign tumours, fibropapillomatosis is a debilitating disease that predominantly afflicts the endangered green turtle (Chelonia mydas). A growing body of evidence has associated these tumours with a herpesvirus. However, a recent study detected both herpesvirus and papillomavirus in these tumours. This result challenged the idea that the herpesvirus is the sole virus associated with this disease. The present study aimed to better understand the co-occurrence of these viruses in turtles with fibropapillomatosis (in both tumour samples and non-tumoured skin samples), in addition to samples from non-tumoured turtles. Both viruses were detected in all sample types, with the 43.5% of tumours containing both herpesvirus and papillomavirus. Tumour samples were found to contain the most herpesvirus while the highest amount of papillomavirus was detected in non-tumoured skin from turtles with tumours. Collectively, these results pivot the way we think about this disease; as an infectious disease where two separate viruses may be at play.

Abstract

Characterised by benign tumours, fibropapillomatosis (FP) is a debilitating disease that predominantly afflicts the endangered green turtle (Chelonia mydas). A growing body of histological and molecular evidence has associated FP tumours with Chelonid alphaherpesvirus 5 (ChHV5). However, a recent study which detected both ChHV5 and Chelonia mydas papillomavirus 1 (CmPV1) DNA in FP tumour tissues has challenged this hypothesis. The present study aimed to establish a probe-based qPCR to assess the wider prevalence of CmPV1 and co-occurrence with ChHV5 in 275 marine turtles foraging in waters adjacent to the east coast of Queensland, Australia: three categories: Group A (FP tumours), Group B (non-tumoured skin from FP turtles) and Group C (non-tumoured skin from turtles without FP). Concurrent detection of ChHV5 and CmPV1 DNA is reported for all three categories, where Group A had the highest rate (43.5%). ChHV5 viral loads in Group A were significantly higher than loads seen in Group B and C. This was not the case for CmPV1 where the loads in Group B were highest, followed by Group A. However, the mean CmPV1 load for Group A samples was not significantly different to the mean load reported from Group B or C samples. Collectively, these results pivot the way we think about FP; as an infectious disease where two separate viruses may be at play.