Transcriptome analysis of northern elephant seal (Mirounga angustirostris) muscle tissue provides a novel molecular resource and physiological insights

Time Domains of Hypoxia Responses and -Omics Insights

The ability to respond rapidly to changes in oxygen tension is critical for many forms of life. Challenges to oxygen homeostasis, specifically in the contexts of evolutionary biology and biomedicine, provide important insights into mechanisms of hypoxia adaptation and tolerance. Here we synthesize findings across varying time domains of hypoxia in terms of oxygen delivery, ranging from early animal to modern human evolution and examine the potential impacts of environmental and clinical challenges through emerging multi-omics approaches. We discuss how diverse animal species have adapted to hypoxic environments, how humans vary in their responses to hypoxia (i.e., in the context of high-altitude exposure, cardiopulmonary disease, and sleep apnea), and how findings from each of these fields inform the other and lead to promising new directions in basic and clinical hypoxia research.

Elephant seal muscle cells adapt to sustained glucocorticoid exposure by shifting their metabolic phenotype

Elephant seals experience natural periods of prolonged food deprivation while breeding, molting, and undergoing postnatal development. Prolonged food deprivation in elephant seals increases circulating glucocorticoids without inducing muscle atrophy, but the cellular mechanisms that allow elephant seals to cope with such conditions remain elusive. We generated a cellular model and conducted transcriptomic, metabolic, and morphological analyses to study how seal cells adapt to sustained glucocorticoid exposure. Seal muscle progenitor cells differentiate into contractile myotubes with a distinctive morphology, gene expression profile, and metabolic phenotype. Exposure to dexamethasone at three ascending concentrations for 48 h modulated the expression of six clusters of genes related to structural constituents of muscle and pathways associated with energy metabolism and cell survival. Knockdown of the glucocorticoid receptor (GR) and downstream expression analyses corroborated that GR mediates the observed effects. Dexamethasone also decreased cellular respiration, shifted the metabolic phenotype toward glycolysis, and induced mitochondrial fission and dissociation of mitochondria-endoplasmic reticulum (ER) interactions without decreasing cell viability. Knockdown of DNA damage-inducible transcript 4 (DDIT4), a GR target involved in the dissociation of mitochondria-ER membranes, recovered respiration and modulated antioxidant gene expression in myotubes treated with dexamethasone. These results show that adaptation to sustained glucocorticoid exposure in elephant seal myotubes involves a metabolic shift toward glycolysis, which is supported by alterations in mitochondrial morphology and a reduction in mitochondria-ER interactions, resulting in decreased respiration without compromising cell survival.

Pinniped Ontogeny as a Window into the Comparative Physiology and Genomics of Hypoxia Tolerance

Diving physiology has received considerable scientific attention as it is a central element of the extreme phenotype of marine mammals. Many scientific discoveries have illuminated physiological mechanisms supporting diving, such as massive, internally bound oxygen stores and dramatic cardiovascular regulation. However, the cellular and molecular mechanisms that support the diving phenotype remain mostly unexplored as logistic and legal restrictions limit the extent of scientific manipulation possible. With next-generation sequencing (NGS) tools becoming more widespread and cost-effective, there are new opportunities to explore the diving phenotype. Genomic investigations come with their own challenges, particularly those including cross-species comparisons. Studying the regulatory pathways that underlie diving mammal ontogeny could provide a window into the comparative physiology of hypoxia tolerance. Specifically, in pinnipeds, which shift from terrestrial pups to elite diving adults, there is potential to characterize the transcriptional, epigenetic, and posttranslational differences between contrasting phenotypes while leveraging a common genome. Here we review the current literature detailing the maturation of the diving phenotype in pinnipeds, which has primarily been explored via biomarkers of metabolic capability including antioxidants, muscle fiber typing, and key aerobic and anaerobic metabolic enzymes. We also discuss how NGS tools have been leveraged to study phenotypic shifts within species through ontogeny, and how this approach may be applied to investigate the biochemical and physiological mechanisms that develop as pups become elite diving adults. We conclude with a specific example of the Antarctic Weddell seal by overlapping protein biomarkers with gene regulatory microRNA datasets.

Functional Studies with Primary Cells Provide a System for Genome-to-Phenome Investigations in Marine Mammals

Marine mammals exhibit some of the most dramatic physiological adaptations in their clade and offer unparalleled insights into the mechanisms driving convergent evolution on relatively short time scales. Some of these adaptations, such as extreme tolerance to hypoxia and prolonged food deprivation, are uncommon among most terrestrial mammals and challenge established metabolic principles of supply and demand balance. Non-targeted omics studies are starting to uncover the genetic foundations of such adaptations, but tools for testing functional significance in these animals are currently lacking. Cellular modeling with primary cells represents a powerful approach for elucidating the molecular etiology of physiological adaptation, a critical step in accelerating genome-to-phenome studies in organisms in which transgenesis is impossible (e.g., large-bodied, long-lived, fully aquatic, federally protected species). Gene perturbation studies in primary cells can directly evaluate whether specific mutations, gene loss, or duplication confer functional advantages such as hypoxia or stress tolerance in marine mammals. Here, we summarize how genetic and pharmacological manipulation approaches in primary cells have advanced mechanistic investigations in other non-traditional mammalian species, and highlight the need for such investigations in marine mammals. We also provide key considerations for isolating, culturing, and conducting experiments with marine mammal cells under conditions that mimic in vivo states. We propose that primary cell culture is a critical tool for conducting functional mechanistic studies (e.g., gene knockdown, over-expression, or editing) that can provide the missing link between genome- and organismal-level understanding of physiological adaptations in marine mammals.

Comprehensive Stress-Based De Novo Transcriptome Assembly and Annotation of Guar (Cyamopsis tetragonoloba (L.) Taub.): An Important Industrial and Forage Crop

The forage crop Guar (Cyamopsis tetragonoloba (L.) Taub.) has the ability to endure heat, drought, and mild salinity. A complete image on its genic architecture will promote our understanding about gene expression networks and different tolerance mechanisms at the molecular level. Therefore, whole mRNA sequence approach on the Guar plant was conducted to provide a snapshot of the mRNA information in the cell under salinity, heat, and drought stresses to be integrated with previous transcriptomic studies. RNA-Seq technology was employed to perform a 2 × 100 paired-end sequencing using an Illumina HiSeq 2500 platform for the transcriptome of leaves of C. tetragonoloba under normal, heat, drought, and salinity conditions. Trinity was used to achieve a de novo assembly followed by gene annotation, functional classification, metabolic pathway analysis, and identification of SSR markers. A total of 218.2 million paired-end raw reads (~44 Gbp) were generated. Of those, 193.5M paired-end reads of high quality were used to reconstruct a total of 161,058 transcripts (~266 Mbp) with N50 of 2552 bp and 61,508 putative genes. There were 6463 proteins having >90% full-length coverage against the Swiss-Prot database and 94% complete orthologs against Embryophyta. Approximately, 62.87% of transcripts were blasted, 50.46% mapped, and 43.50% annotated. A total of 4715 InterProScan families, 3441 domains, 74 repeats, and 490 sites were detected. Biological processes, molecular functions, and cellular components comprised 64.12%, 25.42%, and 10.4%, respectively. The transcriptome was associated with 985 enzymes and 156 KEGG pathways. A total of 27,066 SSRs were gained with an average frequency of one SSR/9.825 kb in the assembled transcripts. This resulting data will be helpful for the advanced analysis of Guar to multi-stress tolerance.

Blubber transcriptome responses to repeated ACTH administration in a marine mammal

Chronic physiological stress impacts animal fitness by catabolizing metabolic stores and suppressing reproduction. This can be especially deleterious for capital breeding carnivores such as marine mammals, with potential for ecosystem-wide effects. However, the impacts and indicators of chronic stress in animals are currently poorly understood. To identify downstream mediators of repeated stress responses in marine mammals, we administered adrenocorticotropic hormone (ACTH) once daily for four days to free-ranging juvenile northern elephant seals (Mirounga angustirostris) to stimulate endogenous corticosteroid release, and compared blubber tissue transcriptome responses to the first and fourth ACTH administrations. Gene expression profiles were distinct between blubber responses to single and repeated ACTH administration, despite similarities in circulating cortisol profiles. We identified 61 and 12 genes that were differentially expressed (DEGs) in response to the first ACTH and fourth administrations, respectively, 24 DEGs between the first and fourth pre-ACTH samples, and 12 DEGs between ACTH response samples from the first and fourth days. Annotated DEGs were associated with functions in redox and lipid homeostasis, suggesting potential negative impacts of repeated stress on capital breeding, diving mammals. DEGs identified in this study are potential markers of repeated stress in marine mammals, which may not be detectable by endocrine profiles alone.

Adipose transcriptome analysis provides novel insights into molecular regulation of prolonged fasting in northern elephant seal pups

The physiological and cellular adaptations to extreme fasting in northern elephant seals ( Mirounga angustirostris, NES) are remarkable and may help to elucidate endocrine mechanisms that regulate lipid metabolism and energy homeostasis in mammals. Recent studies have highlighted the importance of thyroid hormones in the maintenance of a lipid-based metabolism during prolonged fasting in weaned NES pups. To identify additional molecular regulators of fasting, we used a transcriptomics approach to examine changes in global gene expression profiles before and after 6-8 wk of fasting in weaned NES pups. We produced a de novo assembly and identified 98 unique protein-coding genes that were differentially expressed between early and late fasting. Most of the downregulated genes were associated with lipid, carbohydrate, and protein metabolism. A number of downregulated genes were also associated with maintenance of the extracellular matrix, consistent with tissue remodeling during weight loss and the multifunctional nature of blubber tissue, which plays both metabolic and structural roles in marine mammals. Using this data set, we predict potential mechanisms by which NES pups sustain metabolism and regulate adipose stores throughout the fast, and provide a valuable resource for additional studies of extreme metabolic adaptations in mammals.

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.

Learning to starve: impacts of food limitation beyond the stress period

Starvation is common among wild animal populations, and many individuals experience repeated bouts of starvation over the course of their lives. Although much information has been gained through laboratory studies of acute starvation, little is known about how starvation affects an animal once food is again available (i.e. during the refeeding and recovery phases). Many animals exhibit a curious phenomenon – some seem to ‘get better’ at starving following exposure to one or more starvation events – by this we mean that they exhibit potentially adaptive responses, including reduced rates of mass loss, reduced metabolic rates, and lower costs of digestion. During subsequent refeedings they may also exhibit improved digestive efficiency and more rapid mass gain. Importantly, these responses can last until the next starvation bout or even be inherited and expressed in the subsequent generation. Currently, however, little is known about the molecular regulation and physiological mechanisms underlying these changes. Here, we identify areas of research that can fill in the most pressing knowledge gaps. In particular, we highlight how recently refined techniques (e.g. stable isotope tracers, quantitative magnetic resonance and thermal measurement) as well as next-generation sequencing approaches (e.g. RNA-seq, proteomics and holobiome sequencing) can address specific starvation-focused questions. We also describe outstanding unknowns ripe for future research regarding the timing and severity of starvation, and concerning the persistence of these responses and their interactions with other ecological stressors.

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.

De novo assembly of the ringed seal (Pusa hispida) blubber transcriptome: A tool that enables identification of molecular health indicators associated with PCB exposure

The ringed seal, Pusa hispida, is a keystone species in the Arctic marine ecosystem, and is proving a useful marine mammal for linking polychlorinated biphenyl (PCB) exposure to toxic injury. We report here the first de novo assembled transcriptome for the ringed seal (342,863 transcripts, of which 53% were annotated), which we then applied to a population of ringed seals exposed to a local PCB source in Arctic Labrador, Canada. We found an indication of energy metabolism imbalance in local ringed seals (n=4), and identified five significant gene transcript targets: plasminogen receptor (Plg-R(KT)), solute carrier family 25 member 43 receptor (Slc25a43), ankyrin repeat domain-containing protein 26-like receptor (Ankrd26), HIS30 (not yet annotated) and HIS16 (not yet annotated) that may represent indicators of PCB exposure and effects in marine mammals. The abundance profiles of these five gene targets were validated in blubber samples collected from 43 ringed seals using a qPCR assay. The mRNA transcript levels for all five gene targets, (Plg-R(KT), r²=0.43), (Slc25a43, r²=0.51), (Ankrd26, r²=0.43), (HIS30, r²=0.39) and (HIS16, r²=0.31) correlated with increasing levels of blubber PCBs. Results from the present study contribute to our understanding of PCB associated effects in marine mammals, and provide new tools for future molecular and toxicology work in pinnipeds.

Blubber transcriptome response to acute stress axis activation involves transient changes in adipogenesis and lipolysis in a fasting-adapted marine mammal

Stress can compromise an animal’s ability to conserve metabolic stores and participate in energy-demanding activities that are critical for fitness. Understanding how wild animals, especially those already experiencing physiological extremes (e.g. fasting), regulate stress responses is critical for evaluating the impacts of anthropogenic disturbance on physiology and fitness, key challenges for conservation. However, studies of stress in wildlife are often limited to baseline endocrine measurements and few have investigated stress effects in fasting-adapted species. We examined downstream molecular consequences of hypothalamic-pituitary-adrenal (HPA) axis activation by exogenous adrenocorticotropic hormone (ACTH) in blubber of northern elephant seals due to the ease of blubber sampling and its key role in metabolic regulation in marine mammals. We report the first phocid blubber transcriptome produced by RNAseq, containing over 140,000 annotated transcripts, including metabolic and adipocytokine genes of interest. The acute response of blubber to stress axis activation, measured 2 hours after ACTH administration, involved highly specific, transient (lasting <24 hours) induction of gene networks that promote lipolysis and adipogenesis in mammalian adipocytes. Differentially expressed genes included key adipogenesis factors which can be used as blubber-specific markers of acute stress in marine mammals of concern for which sampling of other tissues is not possible.

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.

Genomic Methods Take the Plunge: Recent Advances in High-Throughput Sequencing of Marine Mammals

The dramatic increase in the application of genomic techniques to non-model organisms (NMOs) over the past decade has yielded numerous valuable contributions to evolutionary biology and ecology, many of which would not have been possible with traditional genetic markers. We review this recent progression with a particular focus on genomic studies of marine mammals, a group of taxa that represent key macroevolutionary transitions from terrestrial to marine environments and for which available genomic resources have recently undergone notable rapid growth. Genomic studies of NMOs utilize an expanding range of approaches, including whole genome sequencing, restriction site-associated DNA sequencing, array-based sequencing of single nucleotide polymorphisms and target sequence probes (e.g., exomes), and transcriptome sequencing. These approaches generate different types and quantities of data, and many can be applied with limited or no prior genomic resources, thus overcoming one traditional limitation of research on NMOs. Within marine mammals, such studies have thus far yielded significant contributions to the fields of phylogenomics and comparative genomics, as well as enabled investigations of fitness, demography, and population structure. Here we review the primary options for generating genomic data, introduce several emerging techniques, and discuss the suitability of each approach for different applications in the study of NMOs.

Transcriptome analysis of sika deer in China

Sika deer is of great commercial value because their antlers are used in tonics and alternative medicine and their meat is healthy and delicious. The goal of this study was to generate transcript sequences from sika deer for functional genomic analyses and to identify the transcripts that demonstrate tissue-specific, age-dependent differential expression patterns. These sequences could enhance our understanding of the molecular mechanisms underlying sika deer growth and development. In the present study, we performed de novo transcriptome assembly and profiling analysis across ten tissue types and four developmental stages (juvenile, adolescent, adult, and aged) of sika deer, using Illumina paired-end tag (PET) sequencing technology. A total of 1,752,253 contigs with an average length of 799 bp were generated, from which 1,348,618 unigenes with an average length of 590 bp were defined. Approximately 33.2 % of these (447,931 unigenes) were then annotated in public protein databases. Many sika deer tissue-specific, age-dependent unigenes were identified. The testes have the largest number of tissue-enriched unigenes, and some of them were prone to develop new functions for other tissues. Additionally, our transcriptome revealed that the juvenile-adolescent transition was the most complex and important stage of the sika deer life cycle. The present work represents the first multiple tissue transcriptome analysis of sika deer across four developmental stages. The generated data not only provide a functional genomics resource for future biological research on sika deer but also guide the selection and manipulation of genes controlling growth and development.

Muscle transcriptome response to ACTH administration in a free-ranging marine mammal

While much of our understanding of stress physiology is derived from biomedical studies, little is known about the downstream molecular consequences of adaptive stress responses in free-living animals. We examined molecular effectors of the stress hormones cortisol and aldosterone in the northern elephant seal, a free-ranging study system in which extreme physiological challenges and cortisol fluctuations are a routine part of life history. We stimulated the neuroendocrine stress axis by administering exogenous adrenocorticotropic hormone (ACTH) and examined the resultant effects by measuring corticosteroid hormones, metabolites, and gene expression before, during, and following administration. ACTH induced an elevation in cortisol, aldosterone, glucose, and fatty acids within 2 h, with complete recovery observed within 24 h of administration. The global transcriptional response of elephant seal muscle tissue to ACTH was evaluated by transcriptomics and involved upregulation of a highly coordinated network of conserved glucocorticoid (GC) target genes predicted to promote metabolic substrate availability without causing deleterious effects seen in laboratory animals. Transcriptional recovery from ACTH was characterized by downregulation of GC target genes and restoration of cell proliferation, metabolism, and tissue maintenance pathways within 24 h. Differentially expressed genes included several adipokines not previously described in muscle, reflecting unique metabolic physiology in fasting-adapted animals. This study represents one of the first transcriptome analyses of cellular responses to hypothalamic-pituitary-adrenal axis stimulation in a free-living marine mammal and suggests that compensatory, tissue-sparing mechanisms may enable marine mammals to maintain cortisol and aldosterone sensitivity while avoiding deleterious long-term consequences of stress.