A blubber gene expression index for evaluating stress in marine mammals

Blubber gene expression and cortisol concentrations reveal changing physiological stress in a Southern ocean sentinel species

The health of migratory eastern Australian humpback whales (Megaptera novaeangliae) can reflect the condition of their remote polar foraging environments. This study used gene expression (LEP, LEPR, ADIQ, AhR, TNF-α, HSP-70), blubber hormone concentrations (cortisol, testosterone), and photogrammetric body condition to assess this sentinel species during a period of unprecedented changes to anthropogenic activity and natural processes. The results revealed higher cortisol concentrations in 2020 compared to 2021, suggesting a decline in physiological stress between years. Additionally, metabolic transcripts LEPR, and AhR, which is also linked to xenobiotic metabolism, were upregulated during the 2020 southbound migration. These differences suggest that one or more environmental stressors were reduced between 2020 and 2021, with upregulated AhR possibly indicating a Southern Ocean pollutant declined between the years. This research confirms a Southern Ocean-wide decrease in whale stress during the study period and informs efforts to identify key stressors on Antarctic marine ecosystems.

A novel RT-qPCR health assay reveals differential expression of stress and immunoregulatory genes between the seasonal migrations of humpback whales (Megaptera novaeangliae)

Health information is essential for the conservation management of whale species. However, assessing the health of free-ranging whales is challenging as samples are primarily limited to skin and blubber tissue. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) offers a method to measure health from blubber RNA, providing insights into energetic status, stress and immune activity. To identify changes in health, natural differences in baseline gene expression linked to an individual's sex, reproductive status and life-history stage must first be quantified. This study aimed to establish baseline gene expression indices of health in migrating humpback whales (Megaptera novaeangliae). To do this, we developed an assay to quantify seven health-related gene transcripts (Leptin, Leptin Receptor, Adiponectin, Aryl Hydrocarbon Receptor, Tumour Necrosis Factor-α, Interleukin-6, Heat Shock Protein-70) and used Bayesian mixed effect models to assess differential baseline expression based on sex, lactation status and migration stage (northbound to and southbound from the annual breeding grounds). Results showed no significant contribution of sex to differential baseline expression. However, lactating individuals exhibited downregulated AhR and HSP-70 compared to non-lactating conspecifics. Additionally, southbound individuals demonstrated significantly upregulated HSP-70 and downregulated TNF-alpha, suggesting a relationship between these inflammation-linked transcripts and migratory fasting. Our results suggest that baseline differences due to migratory stage and lactation status should be considered in health applications of this assay. Future monitoring efforts can use our baseline measurements to better understand how gene expression is tied to population-level impacts, such as reduced prey availability or migratory stressors.

Adrenal response to ACTH challenge alters thyroid and immune function and varies with body reserves in molting adult female northern elephant seals

Intrinsic stressors associated with life-history stages may alter the responsiveness of the hypothalamic-pituitary-adrenal axis and responses to extrinsic stressors. We administered adrenocorticotropic hormone (ACTH) to 24 free-ranging adult female northern elephant seals (NESs) at two life-history stages: early and late in their molting period and measured a suite of endocrine, immune, and metabolite responses. Our objective was to evaluate the impact of extended, high-energy fasting on adrenal responsiveness. Animals were blood sampled every 30 min for 120 min post-ACTH injection, then blood was sampled 24 h later. In response to ACTH injection, cortisol levels increased 8- to 10-fold and remained highly elevated compared with baseline at 24 h. Aldosterone levels increased 6- to 9-fold before returning to baseline at 24 h. The magnitude of cortisol and aldosterone release were strongly associated, and both were greater after extended fasting. We observed an inverse relationship between fat mass and the magnitude of cortisol and aldosterone responses, suggesting that body reserves influenced adrenal responsiveness. Sustained elevation in cortisol was associated with alterations in thyroid hormones; both tT3 and tT4 concentrations were suppressed at 24 h, while rT3 increased. Immune cytokine IL-1β was also suppressed after 24 h of cortisol elevation, and numerous acute and sustained impacts on substrate metabolism were evident. Our data suggest that female NESs are more sensitive to stress after the molt fast and that acute stress events can have important impacts on metabolism and immune function. These findings highlight the importance of considering life-history context when assessing the impacts of anthropogenic stressors on wildlife.

Comprehensive molecular and morphological resolution of blubber stratification in a deep-diving, fasting-adapted seal

Blubber is a modified subcutaneous adipose tissue in marine mammals that provides energy storage, thermoregulation, hydrodynamic locomotion, and buoyancy. Blubber displays vertical stratification by lipid content, fatty acid composition, and vascularization, leading to the assumption that deeper blubber layers are metabolically active, while superficial layers are mainly structural and thermoregulatory. However, few studies have examined functional stratification of marine mammal blubber directly, especially in pinnipeds. We characterized morphological and transcriptional differences across blubber layers in the northern elephant seal, a deep-diving and fasting-adapted phocid. We collected blubber from seals early in their fasting period and divided blubber cores into three similarly sized portions. We hypothesized that the innermost blubber portion would have higher 1) heterogeneity in adipocyte size, 2) microvascular density, and 3) expression of genes associated with metabolism and hormone signaling than outer blubber. We found that adipocyte area and variance increased from outermost (skin-adjacent) to innermost (muscle-adjacent) blubber layers, suggesting that inner blubber has a higher capacity for lipid storage and turnover than outer blubber. Inner blubber had a higher proportion of CD144+ endothelial cells, suggesting higher microvascular density. In contrast, outer blubber had a higher proportion of CD4+ immune cells than inner blubber, suggesting higher capacity for response to tissue injury. Transcriptome analysis identified 61 genes that were differentially expressed between inner and outer blubber layers, many of which have not been studied previously in marine mammals. Based on known functions of these genes in other mammals, we suggest that inner blubber has potentially higher 1) adipogenic capacity, 2) cellular diversity, and 3) metabolic and neuroendocrine signaling activity, while outer blubber may have higher 1) extracellular matrix synthesis activity and 2) responsiveness to pathogens and cell stressors. We further characterized expression of nine genes of interest identified by transcriptomics and two adipokines with higher precision across blubber layers using targeted assays. Our study provides functional insights into stratification of blubber in marine mammals and a molecular key, including CD144, CD4, HMGCS2, GABRG2, HCAR2, and COL1A2, for distinguishing blubber layers for physiological and functional studies in seals.

Ontogeny of Carbon Monoxide-Related Gene Expression in a Deep-Diving Marine Mammal

Marine mammals such as northern elephant seals (NES) routinely experience hypoxemia and ischemia-reperfusion events to many tissues during deep dives with no apparent adverse effects. Adaptations to diving include increased antioxidants and elevated oxygen storage capacity associated with high hemoprotein content in blood and muscle. The natural turnover of heme by heme oxygenase enzymes (encoded by HMOX1 and HMOX2) produces endogenous carbon monoxide (CO), which is present at high levels in NES blood and has been shown to have cytoprotective effects in laboratory systems exposed to hypoxia. To understand how pathways associated with endogenous CO production and signaling change across ontogeny in diving mammals, we measured muscle CO and baseline expression of 17 CO-related genes in skeletal muscle and whole blood of three age classes of NES. Muscle CO levels approached those of animals exposed to high exogenous CO, increased with age, and were significantly correlated with gene expression levels. Muscle expression of genes associated with CO production and antioxidant defenses (HMOX1, BVR, GPX3, PRDX1) increased with age and was highest in adult females, while that of genes associated with protection from lipid peroxidation (GPX4, PRDX6, PRDX1, SIRT1) was highest in adult males. In contrast, muscle expression of mitochondrial biogenesis regulators (PGC1A, ESRRA, ESRRG) was highest in pups, while genes associated with inflammation (HMOX2, NRF2, IL1B) did not vary with age or sex. Blood expression of genes involved in regulation of inflammation (IL1B, NRF2, BVR, IL10) was highest in pups, while HMOX1, HMOX2 and pro-inflammatory markers (TLR4, CCL4, PRDX1, TNFA) did not vary with age. We propose that ontogenetic upregulation of baseline HMOX1 expression in skeletal muscle of NES may, in part, underlie increases in CO levels and expression of genes encoding antioxidant enzymes. HMOX2, in turn, may play a role in regulating inflammation related to ischemia and reperfusion in muscle and circulating immune cells. Our data suggest putative ontogenetic mechanisms that may enable phocid pups to transition to a deep-diving lifestyle, including high baseline expression of genes associated with mitochondrial biogenesis and immune system activation during postnatal development and increased expression of genes associated with protection from lipid peroxidation in adulthood.

Species and population specific gene expression in blood transcriptomes of marine turtles

BACKGROUND

Transcriptomic data has demonstrated utility to advance the study of physiological diversity and organisms' responses to environmental stressors. However, a lack of genomic resources and challenges associated with collecting high-quality RNA can limit its application for many wild populations. Minimally invasive blood sampling combined with de novo transcriptomic approaches has great potential to alleviate these barriers. Here, we advance these goals for marine turtles by generating high quality de novo blood transcriptome assemblies to characterize functional diversity and compare global transcriptional profiles between tissues, species, and foraging aggregations.

RESULTS

We generated high quality blood transcriptome assemblies for hawksbill (Eretmochelys imbricata), loggerhead (Caretta caretta), green (Chelonia mydas), and leatherback (Dermochelys coriacea) turtles. The functional diversity in assembled blood transcriptomes was comparable to those from more traditionally sampled tissues. A total of 31.3% of orthogroups identified were present in all four species, representing a core set of conserved genes expressed in blood and shared across marine turtle species. We observed strong species-specific expression of these genes, as well as distinct transcriptomic profiles between green turtle foraging aggregations that inhabit areas of greater or lesser anthropogenic disturbance.

CONCLUSIONS

Obtaining global gene expression data through non-lethal, minimally invasive sampling can greatly expand the applications of RNA-sequencing in protected long-lived species such as marine turtles. The distinct differences in gene expression signatures between species and foraging aggregations provide insight into the functional genomics underlying the diversity in this ancient vertebrate lineage. The transcriptomic resources generated here can be used in further studies examining the evolutionary ecology and anthropogenic impacts on marine turtles.