Skin Transcriptomes of common bottlenose dolphins (Tursiops truncatus) from the northern Gulf of Mexico and southeastern U.S. Atlantic coasts.

Common bottlenose dolphins serve as sentinels for the health of their coastal environments as they are susceptible to health impacts from anthropogenic inputs through both direct exposure and food web magnification. Remote biopsy samples have been widely used to reveal contaminant burdens in free-ranging bottlenose dolphins, but do not address the health consequences of this exposure. To gain insight into whether remote biopsies can also identify health impacts associated with contaminant burdens, we employed RNA sequencing (RNA-seq) to interrogate the transcriptomes of remote skin biopsies from 116 bottlenose dolphins from the northern Gulf of Mexico and southeastern U.S. Atlantic coasts. Gene expression was analyzed using principal component analysis, differential expression testing, and gene co-expression networks, and the results correlated to season, location, and contaminant burden. Season had a significant impact, with over 60% of genes differentially expressed between spring/summer and winter months. Geographic location exhibited lesser effects on the transcriptome, with 23.5% of genes differentially expressed between the northern Gulf of Mexico and the southeastern U.S. Atlantic locations. Despite a large overlap between the seasonal and geographical gene sets, the pathways altered in the observed gene expression profiles were somewhat distinct. Co-regulated gene modules and differential expression analysis both identified epidermal development and cellular architecture pathways to be expressed at lower levels in animals from the northern Gulf of Mexico. Although contaminant burdens measured were not significantly different between regions, some correlation with contaminant loads in individuals was observed among co-expressed gene modules, but these did not include classical detoxification pathways. Instead, this study identified other, possibly downstream pathways, including those involved in cellular architecture, immune response, and oxidative stress, that may prove to be contaminant responsive markers in bottlenose dolphin skin.

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.

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.

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.

Transcriptomic responses of juvenile Pacific whiteleg shrimp, Litopenaeus vannamei, to hypoxia and hypercapnic hypoxia.

Estuarine crustaceans are often exposed to low dissolved O2 (hypoxia) accompanied by elevated CO2 (hypercapnia), which lowers water pH. Acclimatory responses to hypoxia have been widely characterized; responses to hypercapnia in combination with hypoxia (hypercapnic hypoxia) are less well known. Here we used oligonucleotide microarrays to characterize changes in global gene expression in the hepatopancreas of Pacific whiteleg shrimp, Litopenaeus vannamei, exposed to hypoxia or hypercapnic hypoxia for 4 or 24 h, compared with time-matched animals held in air-saturated water (normoxia). Unigenes whose expressions were significantly impacted by treatment and/or time were used to build artificial neural networks (ANNs) to identify genes with the greatest sensitivity in pairwise discriminations between treatments at each time point and between times for each treatment. ANN gene sets that discriminated hypoxia or hypercapnic hypoxia from normoxia shared functions of translation, mitochondrial energetics, and cellular defense. GO terms protein modification/phosphorylation/cellular protein metabolism and RNA processing/apoptosis/cell cycling occurred at highest frequency in discriminating hypercapnic hypoxia from hypoxia at 4 and 24 h, respectively. For 75.4% of the annotated ANN genes, exposure to hypercapnic hypoxia for 24 h reduced or reversed the transcriptional response to hypoxia alone. These results suggest that high CO2/low pH may interfere with transcriptionally based acclimation to hypoxia or elicit physiological or biochemical responses that relieve internal hypoxia. Whether these data reflect resilience or sensitivity of L. vannamei in the face of expanding hypoxic zones and rising levels of atmospheric CO2 may be important to understanding the survival of this and other estuarine species.

Transcriptome remodeling associated with chronological aging in the dinoflagellate, Karenia brevis.

The toxic dinoflagellate, Karenia brevis, forms dense blooms in the Gulf of Mexico that persist for many months in coastal waters, where they can cause extensive marine animal mortalities and human health impacts. The mechanisms that enable cell survival in high density, low growth blooms, and the mechanisms leading to often rapid bloom demise are not well understood. To gain an understanding of processes that underlie chronological aging in this dinoflagellate, a microarray study was carried out to identify changes in the global transcriptome that accompany the entry and maintenance of stationary phase up to the onset of cell death. The transcriptome of K. brevis was assayed using a custom 10,263 feature oligonucleotide microarray from mid-logarithmic growth to the onset of culture demise. A total of 2958 (29%) features were differentially expressed, with the mid-stationary phase timepoint demonstrating peak changes in expression. Gene ontology enrichment analyses identified a significant shift in transcripts involved in energy acquisition, ribosome biogenesis, gene expression, stress adaptation, calcium signaling, and putative brevetoxin biosynthesis. The extensive remodeling of the transcriptome observed in the transition into a quiescent non-dividing phase appears to be indicative of a global shift in the metabolic and signaling requirements and provides the basis from which to understand the process of chronological aging in a dinoflagellate.

Arsenic exposure to killifish during embryogenesis alters muscle development.

Epidemiological studies have correlated arsenic exposure in drinking water with adverse developmental outcomes such as stillbirths, spontaneous abortions, neonatal mortality, low birth weight, delays in the use of musculature, and altered locomotor activity. Killifish (Fundulus heteroclitus) were used as a model to help to determine the mechanisms by which arsenic could impact development. Killifish embryos were exposed to three different sodium arsenite concentrations and were collected at 32 h post-fertilization (hpf), 42 hpf, 168 hpf, or < 24 h post-hatch. A killifish oligo microarray was developed and used to examine gene expression changes between control and 25-ppm arsenic-exposed hatchlings. With artificial neural network analysis of the transcriptomic data, accurate prediction of each group (control vs. arsenic-exposed embryos) was obtained using a small subset of only 332 genes. The genes differentially expressed include those involved in cell cycle, development, ubiquitination, and the musculature. Several of the genes involved in cell cycle regulation and muscle formation, such as fetuin B, cyclin D-binding protein 1, and CapZ, were differentially expressed in the embryos in a time- and dose-dependent manner. Examining muscle structure in the hatchlings showed that arsenic exposure during embryogenesis significantly reduces the average muscle fiber size, which is coupled with a significant 2.1- and 1.6-fold upregulation of skeletal myosin light and heavy chains, respectively. These findings collectively indicate that arsenic exposure during embryogenesis can initiate molecular changes that appear to lead to aberrant muscle formation.