Machine learning approaches to investigate the impact of PCBs on the transcriptome of the common bottlenose dolphin (Tursiops truncatus)

Transcriptome profiling of blood from common bottlenose dolphins (Tursiops truncatus) in the northern Gulf of Mexico to enhance health assessment capabilities

Following the 2010 Deepwater Horizon disaster and subsequent unusual mortality event, adverse health impacts have been reported in bottlenose dolphins in Barataria Bay, LA including impaired stress response and reproductive, pulmonary, cardiac, and immune function. These conditions were primarily diagnosed through hands-on veterinary examinations and analysis of standard diagnostic panels. In human and veterinary medicine, gene expression profiling has been used to identify molecular mechanisms underlying toxic responses and disease states. Identification of molecular markers of exposure or disease may enable earlier detection of health effects or allow for health evaluation when the use of specialized methodologies is not feasible. To date this powerful tool has not been applied to augment the veterinary data collected concurrently during dolphin health assessments. This study examined transcriptomic profiles of blood from 76 dolphins sampled in health assessments during 2013–2018 in the waters near Barataria Bay, LA and Sarasota Bay, FL. Gene expression was analyzed in conjunction with the substantial suite of health data collected using principal component analysis, differential expression testing, over-representation analysis, and weighted gene co-expression network analysis. Broadly, transcript profiles of Barataria Bay dolphins indicated a shift in immune response, cytoskeletal alterations, and mitochondrial dysfunction, most pronounced in dolphins likely exposed to Deepwater Horizon oiling. While gene expression profiles in Barataria Bay dolphins were altered compared to Sarasota Bay for all years, profiles from 2013 exhibited the greatest alteration in gene expression. Differentially expressed transcripts included genes involved in immunity, inflammation, reproductive failure, and lung or cardiac dysfunction, all of which have been documented in dolphins from Barataria Bay following the Deepwater Horizon oil spill. The genes and pathways identified in this study may, with additional research and validation, prove useful as molecular markers of exposure or disease to assist wildlife veterinarians in evaluating the health of dolphins and other cetaceans.

A review of the application of machine learning in water quality evaluation

With the rapid increase in the volume of data on the aquatic environment, machine learning has become an important tool for data analysis, classification, and prediction. Unlike traditional models used in water-related research, data-driven models based on machine learning can efficiently solve more complex nonlinear problems. In water environment research, models and conclusions derived from machine learning have been applied to the construction, monitoring, simulation, evaluation, and optimization of various water treatment and management systems. Additionally, machine learning can provide solutions for water pollution control, water quality improvement, and watershed ecosystem security management. In this review, we describe the cases in which machine learning algorithms have been applied to evaluate the water quality in different water environments, such as surface water, groundwater, drinking water, sewage, and seawater. Furthermore, we propose possible future applications of machine learning approaches to water environments.

Comparative Analysis of Blood Transcriptome in the Yangtze Finless Porpoise (Neophocaena asiaeorientalis)

The Yangtze finless porpoise (Neophocaena asiaeorientalis) is the sole freshwater subspecies of Neophocaenaphocaenoides, and there is a lack of data on its transcriptome. In this study, we applied RNA-seq technology to assemble, de novo, a transcriptome and analyzed differential expressed genes (DEGs). About 6 Gb of clean data was generated for the Yangtze finless porpoise blood (n = 6) through de novo sequencing. In total, 151,211 unigenes were generated and a total of 119,039 of these unigenes (78.72%) were functionally annotated when searched for within the NCBI Nr, SwissProt, GO, COG, and KEGG databases. Diverse and extensive expressed gene catalogs were sampled for the Yangtze finless porpoise. DESeq2 was used to analyze the differential expression genes (DEGs) obtained from the assembled transcriptome. The results indicated that DEGs have close relationships with the Yangtze finless porpoise’s development, evolution and adaptation. Further, we found that genes involved in cetacean TAG synthesis might directly explain the molecular basis of cetacean blubber thickening. These transcriptome data will assist in understanding molecular mechanisms of Yangtze finless porpoise adaptation.

Multi-Omics Approaches and Radiation on Lipid Metabolism in Toothed Whales

Lipid synthesis pathways of toothed whales have evolved since their movement from the terrestrial to marine environment. The synthesis and function of these endogenous lipids and affecting factors are still little understood. In this review, we focused on different omics approaches and techniques to investigate lipid metabolism and radiation impacts on lipids in toothed whales. The selected literature was screened, and capacities, possibilities, and future approaches for identifying unusual lipid synthesis pathways by omics were evaluated. Omics approaches were categorized into the four major disciplines: lipidomics, transcriptomics, genomics, and proteomics. Genomics and transcriptomics can together identify genes related to unique lipid synthesis. As lipids interact with proteins in the animal body, lipidomics, and proteomics can correlate by creating lipid-binding proteome maps to elucidate metabolism pathways. In lipidomics studies, recent mass spectroscopic methods can address lipid profiles; however, the determination of structures of lipids are challenging. As an environmental stress, the acoustic radiation has a significant effect on the alteration of lipid profiles. Radiation studies in different omics approaches revealed the necessity of multi-omics applications. This review concluded that a combination of many of the omics areas may elucidate the metabolism of lipids and possible hazards on lipids in toothed whales by radiation.

Transcriptome analysis of cadmium exposure in kidney fibroblast cells of the North Atlantic Right Whale (Eubalaena glacialis)

An 8X15k oligonucleotide microarray was developed consisting of 2334 Eubalaena glacialis probes and 2166 Tursiops truncatus probes and used to measure the effects, at transcriptomic level, of cadmium exposure in right whale kidney fibroblast cells. Cells were exposed to three concentrations (1 μM, 0.1 μM, and 0.01 μM) of cadmium chloride (CdCl₂) for three exposure times (1, 4, and 24 h). Cells exposed to 1 μM CdCl₂ for 4 h and 24 h showed upregulated genes involved in protection from metal toxicity and oxidative stress, protein renaturation, apoptosis inhibition, as well as several regulators of cellular processes. Downregulated genes represented a suite of functions including cell proliferation, transcription regulation, actin polymerization, and stress fiber synthesis. The collection of differentially expressed genes in this study support proposed mechanisms of cadmium-induced apoptosis such as ubiquitin proteasome system disruption, Ca²⁺ homeostasis interference, mitochondrial membrane potential collapse, reactive oxygen species (ROS) production, and cell cycle arrest. The results also have confirmed the right whale microarray as a reproducible tool in measuring differentiated gene expression that could be a valuable asset for transcriptome analysis of other baleen whales and potential health assessment protocols.

Skin distress associated with xenobiotics exposure: An epigenetic study in the Mediterranean fin whale (Balaenoptera physalus)

The phenotypic plasticity of many organisms is mediated in part by epigenetics, the heritable changes in gene activity that occur without any alterations to DNA sequence. A major mechanism in epigenetics is the DNA methylation (DNAm). Hypo- and hyper-methylation are generalized responses to control gene expression however recent studies have demonstrated that classes of contaminants could mark specific DNAm signatures, that could usefully signal prior environmental exposure. We collected skin and blubber from 6 free-ranging fin whale (Balaenoptera physalus) individuals sampled as a part of a previous published study in the northern Mediterranean Sea. Genomic DNA extracted from the skin of the fin whales and levels of contaminants measured in the blubber of the same individuals were used for DNAm profiling through reduced representation bisulfite sequencing (RRBS). We tested the hypothesis that differences in the methylation patterns could be related to environmental exposure to contaminants and load in the whale tissues. The aims of this study were to determine the DNAm profiles of the methylation contexts (CpGs and non-CpGs) of differently contaminated groups using the RRBS, and to identify potential contaminant exposure related genes. Amount and proportion of methylcytosines in CpG and non-CpG regions (CHH and CHG) was very similar across the 6 samples. The proportion of methylcytosines sites in CpG was n = 32,682, the highest among all the sequence contexts (n = 3216 in CHH; n = 1743 in CHG). The majority of the methylcytosine occurred in the intron regions, followed by exon and promoter regions in CpG, CHH and CHG. Gene Ontology results indicated that DNAm affected genes that take place in cell differentiation and function in cutaneous, vascular and nervous systems. The identification of cellular response pathways allows a better understanding of the organism biological reaction to a specific environmental challenge and the development of sensitive tools based on the predictive responses. Eco-epigenetics analyses have an extraordinary potential to address growing issues on pollution biomonitoring, ecotoxicity assessment, conservation and management planning.

Contaminant Exposure Linked to Cellular and Endocrine Biomarkers in Southern California Bottlenose Dolphins

Cetaceans in the Southern California Bight (SCB) are exposed to high levels of halogenated organic contaminants (HOCs), which have previously been linked to impaired reproductive health and immune responses. We used a combination of molecular tools to examine the potential physiological impacts of HOC exposure in two bottlenose dolphin ( Tursiops truncatus) ecotypes in the SCB. We quantified 25 HOCs in the blubber of 22 biopsies collected from males between 2012 and 2016. We then analyzed genome-wide gene expression in skin using RNA-sequencing and measured blubber testosterone to compare HOC exposure with cellular and endocrine biomarkers. We found high levels of HOCs in both ecotypes with significantly higher total polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), tris(4-chlorophenyl)methanol (TCPMOH), and chlordane-related compounds in the coastal ecotype versus the offshore ecotype. We found evidence of PBDE bioaccumulation in both ecotypes, however, the pattern of bioaccumulation or endocrine disruption for other HOCs was different between the ecotypes, suggesting potential endocrine disruption in the coastal ecotype. We also observed correlations between HOCs and gene coexpression networks enriched for xenobiotic metabolism, hormone metabolism, and immune response that could indicate cellular effects from HOC exposure. By integrating measurements of HOC load with both transcriptome profiling and endocrine biomarkers, our approach provides insight into HOC exposure and potential impacts on wild cetacean health in southern California.

On the revolution of cetacean evolution

The order of Cetacea with 88 species including Odontoceti (or toothed whales) and Mysticeti (or baleen whales) is the most specialized and diversified group of mammals. The blue whale with a maximum recorded length of 29.9 m for 173 t of weight is the largest animal known to have ever existed, and any dolphin's brain is most powerful and complex than any other brain in the animal kingdom, second only to primate's. Nevertheless, Cetacea are mammals that re-entered the oceans only a little over 50 million years ago, a relatively short time on the evolutionary scale. During this time cetaceans and humans have developed marked morphological and behavioral differences, yet their genomes show a high level of similarity. This present review is focused on the description and significance of newly accessible cetacean genome tools and information, and their relevance in the study of the evolution of successful phenotypic adaptations associated to mammal's marine existence, and their applicability to the unresolved disease mechanisms in humans.

New Technologies for Monitoring Marine Mammal Health

Omics technologies have been developed in recent decades and applied to different subjects, although the greatest advancements have been achieved in human biology and disease. Genome sequencing and the exploration of its coding and noncoding regions are rapidly yielding meaningful answers to diverse questions, relating genome information to protein activity to environmental changes. In the past, marine mammal genetic and transcriptional studies have been restricted due to the lack of reference genomes. But the advance of high-throughput sequencing is revolutionizing the life sciences technologies. As long-lived organisms, at the top of the food chain, marine mammals play an important role in marine ecosystems and while their protected status is in favor of conservation of the species, it also complicates the researcher's approach to traditional measurements of health. Omics data generated by high-throughput technologies will represent an important key for improving the scientific basis for understanding both marine mammal and environment health.

De novo transcriptome assembly and RNA-Seq expression analysis in blood from beluga whales of Bristol Bay, AK

Assessing the health of marine mammal sentinel species is crucial to understanding the impacts of environmental perturbations on marine ecosystems and human health. In Arctic regions, beluga whales, Delphinapterus leucas, are upper level predators that may serve as a sentinel species, potentially forecasting impacts on human health. While gene expression profiling from blood transcriptomes has widely been used to assess health status and environmental exposures in human and veterinary medicine, its use in wildlife has been limited due to the lack of available genomes and baseline data. To this end we constructed the first beluga whale blood transcriptome de novo from samples collected during annual health assessments of the healthy Bristol Bay, AK stock during 2012-2014 to establish baseline information on the content and variation of the beluga whale blood transcriptome. The Trinity transcriptome assembly from beluga was comprised of 91,325 transcripts that represented a wide array of cellular functions and processes and was extremely similar in content to the blood transcriptome of another cetacean, the bottlenose dolphin. Expression of hemoglobin transcripts was much lower in beluga (25.6% of TPM, transcripts per million) than has been observed in many other mammals. A T12A amino acid substitution in the HBB sequence of beluga whales, but not bottlenose dolphins, was identified and may play a role in low temperature adaptation. The beluga blood transcriptome was extremely stable between sex and year, with no apparent clustering of samples by principle components analysis and <4% of genes differentially expressed (EBseq, FDR<0.05). While the impacts of season, sexual maturity, disease, and geography on the beluga blood transcriptome must be established, the presence of transcripts involved in stress, detoxification, and immune functions indicate that blood gene expression analyses may provide information on health status and exposure. This study provides a wealth of transcriptomic data on beluga whales and provides a sizeable pool of preliminary data for comparison with other studies in beluga whale.

Marine mammal cell cultures: To obtain, to apply, and to preserve

The world's oceans today have become a place for the disposal of toxic waste, which leads to the degradation of marine mammal habitats and populations. Marine mammal cell cultures have proven to be a multifunctional tool for studying the peculiarities of the cell physiology and biochemistry of these animals as well as the destructive effects of anthropogenic and natural toxicants. This review describes the sources of marine mammal live tissues and the methods required for establishing cell cultures, their use, and long-term storage. Approaches to conserving rare animal species by applying cell biology methodologies are also discussed.

Brain transcriptomes of harbor seals demonstrate gene expression patterns of animals undergoing a metabolic disease and a viral infection

Diseases of marine mammals can be difficult to diagnose because of their life history and protected status. Stranded marine mammals have been a particularly useful resource to discover and comprehend the diseases that plague these top predators. Additionally, advancements in high-throughput sequencing (HTS) has contributed to the discovery of novel pathogens in marine mammals. In this study, we use a combination of HTS and stranded harbor seals (Phoca vitulina) to better understand a known and unknown brain disease. To do this, we used transcriptomics to evaluate brain tissues from seven neonatal harbor seals that expired from an unknown cause of death (UCD) and compared them to four neonatal harbor seals that had confirmed phocine herpesvirus (PhV-1) infections in the brain. Comparing the two disease states we found that UCD animals showed a significant abundance of fatty acid metabolic transcripts in their brain tissue, thus we speculate that a fatty acid metabolic dysregulation contributed to the death of these animals. Furthermore, we were able to describe the response of four young harbor seals with PhV-1 infections in the brain. PhV-1 infected animals showed a significant ability to mount an innate and adaptive immune response, especially to combat viral infections. Our data also suggests that PhV-1 can hijack host pathways for DNA packaging and exocytosis. This is the first study to use transcriptomics in marine mammals to understand host and viral interactions and assess the death of stranded marine mammals with an unknown disease. Furthermore, we show the value of applying transcriptomics on stranded marine mammals for disease characterization.

RNA-Seq analysis of seasonal and individual variation in blood transcriptomes of healthy managed bottlenose dolphins

Background

The blood transcriptome can reflect both systemic exposures and pathological changes in other organs of the body because immune cells recirculate through the blood, lymphoid tissues, and affected sites. In human and veterinary medicine, blood transcriptome analysis has been used successfully to identify markers of disease or pathological conditions, but can be confounded by large seasonal changes in expression. In comparison, the use of transcriptomic based analyses in wildlife has been limited. Here we report a longitudinal study of four managed bottlenose dolphins located in Waikoloa, Hawaii, serially sampled (approximately monthly) over the course of 1 year to establish baseline information on the content and variation of the dolphin blood transcriptome.

Results

Illumina based RNA-seq analyses were carried out using both the Ensembl dolphin genome and a de novo blood transcriptome as guides. Overall, the blood transcriptome encompassed a wide array of cellular functions and processes and was relatively stable within and between animals over the course of 1 year. Principal components analysis revealed moderate clustering by sex associated with the variation among global gene expression profiles (PC1, 22 % of variance). Limited seasonal change was observed, with < 2.5 % of genes differentially expressed between winter and summer months (FDR < 0.05). Among the differentially expressed genes, cosinor analysis identified seasonal rhythmicity for the observed changes in blood gene expression, consistent with studies in humans. While the proportion of seasonally variant genes in these dolphins is much smaller than that reported in humans, the majority of those identified in dolphins were also shown to vary with season in humans. Gene co-expression network analysis identified several gene modules with significant correlation to age, sex, or hematological parameters.

Conclusions

This longitudinal analysis of healthy managed dolphins establishes a preliminary baseline for blood transcriptome analysis in this species. Correlations with hematological parameters, distinct from muted seasonal effects, suggest that the otherwise relatively stable blood transcriptome may be a useful indicator of health and exposure. A robust database of gene expression in free-ranging and managed dolphins across seasons with known adverse health conditions or contaminant exposures will be needed to establish predictive gene expression profiles suitable for biomonitoring.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3020-8) contains supplementary material, which is available to authorized users.

Seasonal variation in the skin transcriptome of common bottlenose dolphins (Tursiops truncatus) from the northern Gulf of Mexico

As long-lived predators that integrate exposures across multiple trophic levels, cetaceans are recognized as sentinels for the health of marine ecosystems. Their utility as sentinels requires the establishment of baseline health parameters. Because cetaceans are protected, measurements obtained with minimal disruption to free ranging animals are highly desirable. In this study we investigated the utility of skin gene expression profiling to monitor health and contaminant exposure in common bottlenose dolphins (Tursiops truncatus). Remote integument biopsies were collected in the northern Gulf of Mexico prior to the Deepwater Horizon oil spill (May 2010) and during summer and winter for two years following oil contamination (2010-2011). A bottlenose dolphin microarray was used to characterize the skin transcriptomes of 94 individuals from three populations: Barataria Bay, Louisiana, Chandeleur Sound, Louisiana, and Mississippi Sound, Mississippi/Alabama. Skin transcriptomes did not differ significantly between populations. In contrast, season had a profound effect on gene expression, with nearly one-third of all genes on the array differing in expression between winter and the warmer seasons (moderated T-test; p<0.01, fold-change≥1.5). Persistent organic pollutants (POPs) in blubber changed concurrently, reaching >two-fold higher concentrations in summer compared to winter, due to a seasonal decrease in blubber thickness and loss of stored lipid. However, global gene expression did not correlate strongly with seasonally changing contaminant concentrations, most likely because the refractory, lipid-stored metabolites are not substrates for phase I or II xenobiotic detoxification pathways. Rather, processes related to cell proliferation, motility, and differentiation dominated the differences in expression in winter and the warmer seasons. More subtle differences were seen between spring and summer (1.5% of genes differentially expressed). However, two presumed oil-exposed animals from spring presented gene expression profiles more similar to the summer animals (presumed exposed) than to other spring animals. Seasonal effects have not previously been considered in studies assessing gene expression in cetaceans, but clearly must be taken into account when applying transcriptomic analyses to investigate their contaminant exposure or health status.