Microarray gene expression profiles from mature gonad tissues of Atlantic bluefin tuna, Thunnus thynnus in the Gulf of Mexico

Transcriptomics integrated with metabolomics reveals the effect of Lycium barbarum polysaccharide on apoptosis in Nile tilapia (Oreochromis niloticus)

Lycium barbarum polysaccharide (LBP) is one of the main active ingredients in the fruit of L. barbarum L. It has been used as herbal medicine for thousands of years in China. In this study, Nile tilapia (Oreochromis niloticus) was taken as the research object. After feeding tilapia with 5 different doses of LBP (0 mg/kg, 500 mg/kg, 1000 mg/kg, 1500 mg/kg, 2000 mg/kg) for 55 d, it was found that LBP could promote the growth of tilapia, and this effect was the strongest at Group 1500 mg/kg. Apoptosis analysis in the liver and spleen showed that dietary supplementation with 1000 mg/kg LBP had the best protective effect on the spleen and liver in tilapia. Combined transcriptomics and metabolomics of the spleen in tilapia at Group 0 mg/kg and 1000 mg/kg showed that the differentially expressed genes (DEGs) such as NT5C2L1, pmm1, FasL and the differentially metabolites such as xanthine, dGMP, guanine and glutamate were mainly concentrated in signaling pathways such as Purine metabolism and FoxO signaling pathway. In conclusion, LBP regulates the metabolic waste levels of tilapia mainly through Purine metabolism and the FoxO signaling pathway, thereby inhibiting cell apoptosis, improving the utilization of nutrients, and promoting the growth of tilapia. This study not only provides a theoretical basis for the application of LBP in aquatic animals but also provides useful information for the healthy development of the aquaculture.

Ancient DNA SNP-panel data suggests stability in bluefin tuna genetic diversity despite centuries of fluctuating catches in the eastern Atlantic and Mediterranean

Atlantic bluefin tuna (Thunnus thynnus; BFT) abundance was depleted in the late 20th and early 21st century due to overfishing. Historical catch records further indicate that the abundance of BFT in the Mediterranean has been fluctuating since at least the 16th century. Here we build upon previous work on ancient DNA of BFT in the Mediterranean by comparing contemporary (2009–2012) specimens with archival (1911–1926) and archaeological (2nd century BCE–15th century CE) specimens that represent population states prior to these two major periods of exploitation, respectively. We successfully genotyped and analysed 259 contemporary and 123 historical (91 archival and 32 archaeological) specimens at 92 SNP loci that were selected for their ability to differentiate contemporary populations or their association with core biological functions. We found no evidence of genetic bottlenecks, inbreeding or population restructuring between temporal sample groups that might explain what has driven catch fluctuations since the 16th century. We also detected a putative adaptive response, involving the cytoskeletal protein synemin which may be related to muscle stress. However, these results require further investigation with more extensive genome-wide data to rule out demographic changes due to overfishing, and other natural and anthropogenic factors, in addition to elucidating the adaptive drivers related to these.

Transcriptome-wide changes associated with the reproductive behaviour of male guppies exposed to 17α-ethinyl estradiol

Although many pharmaceutical compounds (and their metabolites) can induce harmful impacts at the molecular, physiological and behavioural levels, their underlying mechanistic associations have remained largely unexplored. Here, we utilized RNA-Seq to build a whole brain transcriptome profile to examine the impact of a common endocrine disrupting pharmaceutical (17α-ethinyl estradiol, EE2) on reproductive behaviour in wild guppies (Poecilia reticulata). Specifically, we annotated 16,791 coding transcripts in whole brain tissue in relation to the courtship behaviour (i.e. sigmoid display) of EE2 exposed (at environmentally relevant concentration of 8 ng/L for 28-days) and unexposed guppies. Further, we obtained 10,960 assembled transcripts matching in the non-coding orthologous genomes. Behavioural responses were assessed using a standard mate choice experiment, which allowed us to disentangle chemical cues from visual cues. We found that a high proportion of the RNAseq reads aligned back to our de novo assembled transcriptome with 80.59% mapping rate. Behavioural experiments showed that when males were presented only with female visual cues, there was a significant interaction between male treatment and female treatment in the time spent in the preference zone. This is one of the first studies to show that transcriptome-wide changes are associated with the reproductive behaviour of fish: EE2 exposed male guppies that performed high levels of courtship had a gene profile that deviated the most from the other treatment groups, while both non-courting EE2 and control males had similar gene signatures. Using Gene Ontology pathway analysis, our study shows that EE2-exposed males had gene transcripts enriched for pathways associated with altered immunity, starvation, altered metabolism and spermatogenesis. Our study demonstrates that multiple gene networks orchestrate courting behaviour, emphasizing the importance of investigating impacts of pharmaceuticals on gene networks instead of single genes.

Insights on the seasonal variations of reproductive features in the Eastern Atlantic Bluefin Tuna

The Atlantic Bluefin Tuna (ABFT, Thunnus thynnus) is one of the most intensely exploited fisheries resources in the world. In spite of the years of studies on ABFT, basic aspects of its reproductive biology remain uncertain. To gain insight regarding the seasonal changes of the reproductive characteristics of the eastern stock of ABFT, blood and tissue samples were collected from mature specimens caught in the Mediterranean basin during the reproductive (May-June) and non-reproductive season (Oct-Nov). Histological analysis of the gonads of May-June samples indicated that there were females which were actively spawning (contained post-ovulatory follicles) and females that were not actively spawning that had previtellogenic and fully vitellogenic oocytes. In males, testis were at early or late stage of spermatogenesis during the reproductive season. In Oct-Nov, ovaries contained mostly previtellogenic oocytes as well as β and α atretic follicles while the testis predominantly contained spermatogonia and few cysts with spermatocytes and spermatozoa. Gonadosomatic index (GSI) in females was highest among the actively spawning individuals while in males GSI was higher in early and late spermatogenic individuals compared to those that were spent. Plasma sex steroids levels varied with the reproductive season. In females, estradiol (E₂), was higher in May-June while testosterone (T) and progesterone (P) did not vary. In males, E₂ and T were higher in May-June while P levels were similar at the two sampling points. Circulating follicle stimulating hormone (FSH) was higher in Oct-Nov than in May-June both in males and females. Vitellogenin (VTG) was detected in plasma from both males and females during the reproductive season with levels in females significantly higher than in males. VTG was undetected in Oct-Nov samples. Since choriogenesis is an important event during follicle growth, the expression of three genes involved in vitelline envelope formation and hardening was measured and results showed significantly higher levels in ovaries in fish caught in May-June with respect to those sampled in Oct-Nov. In addition, a set of genes encoding for ion channels that are responsible for oocyte hydration and buoyancy, as well as sperm viability, were characterized at the two time points, and these were found to be more highly expressed in females during the reproductive season. Finally, the expression level of three mRNAs encoding for different lipid-binding proteins was analyzed with significantly higher levels detected in males, suggesting sex-specific expression. Our findings provide additional information on the reproductive biology of ABFT, particularly on biomarkers for the assessment of the state of maturation of the gonad, highlighting gender-specific signals and seasonal differences.

Suitable reference gene for quantitative real-time PCR analysis of gene expression in gonadal tissues of minnow Puntius sophore under high-temperature stress

BACKGROUND

High ambient temperature is known to affect fish gonadal development and physiology in a variety of ways depending on the severity and duration of exposure; however, the underlying molecular mechanisms are poorly understood. Gonadal gene expression influence the gonadal development, physiology and the quality of egg/sperm produced in teleosts and the mechanistic understanding of spatio-temporal changes in the gonadal gene expression could be instrumental in controlling the fate of egg/sperm and the quality of seed produced. Real time-quantititative polymerase chain reaction (RT-qCR), is a high throughput, sensitive and reproducible methodology used for understanding gene expression patterns by measuring the relative abundance of mRNA transcripts. However, its accuracy relies upon a suitable reference gene whose expression levels remain stable across various experimental conditions. In the present study, we evaluated the suitability of ten potential reference genes to be used as internal controls in RT-qPCR analysis in gonadal tissues (ovary and testis) of minnow Puntius sophore exposed to high temperature stress for different time periods (7 days, 60 days). Expression analysis of ten different constitutively expressed genes viz. 18S ribosomal RNA (18S rRNA), beta actin (βactin), β-2 microglobulin (b2mg), eukaryotic elongation factor-1 (eef1), glyceraldehyde-3phosphate dehydrogenase (gapdh), glucose-6-phosphate dehydrogenase (g6pd), ribosomal binding protein L13 (rpl13), tubulin (tub), tata box binding protein (tbp), ubiquitin (ubi) was carried out by using RT-qPCR and the stability in their expressions were evaluated by using four different algorithms; namely, delta Ct, BestKeeper, geNorm and NormFinder.

RESULTS

In ovary, eef1 was found to be the most suitable reference gene in all the algorithms used. In testis, b2mg was found to be the most suitable reference gene in delta Ct, BestKeeper, NormFinder analysis while tbp and eef1 were found to be the most suitable reference genes in geNorm analysis.

CONCLUSIONS

In conclusion, eef1 and b2mg were found to be the most suitable reference genes in ovary and testis, respectively, of Puntius sophore exposed to high temperature stress, and could be used as internal controls for gene expression analysis in gonadal tissues of Puntius sophore.

Characterisation of the transcriptome of male and female wild-type guppy brains with RNA-Seq and consequences of exposure to the pharmaceutical pollutant, 17α-ethinyl estradiol

Waterways are increasingly being contaminated by chemical compounds that can disrupt the endocrinology of organisms. One such compound is 17α-ethinyl estradiol (EE2), a synthetic estrogen used in the contraceptive pill. Despite considerable research interest in the effects of EE2 on reproduction and gene expression, surprisingly, only a few studies have capitalised on technologies, such as next-generation sequencing (NGS), to uncover the molecular pathways related to EE2 exposure. Accordingly, using high-throughput sequencing technologies, the aim of our study was to explore the effects of EE2 on brain transcriptome in wild-type male and female guppy (Poecilia reticulata). We conducted two sets of experiments, where fish were exposed to EE2 (measured concentrations: 8ng/L and 38ng/L) in a flow-through system for 21days. The effects on the brain transcriptome on both males and females were assessed using Illumina sequencing (MiSeq and HiSeq) platform followed by bioinformatics analysis (edgeR, DESeq2). Here, we report that exposure to EE2 caused both up- and downregulation of specific transcript abundances, and affected transcript abundance in a sex-specific manner. Specifically, we found 773 transcripts, of which 60 were male-specific, 61 female-specific and 285 treatment-specific. EE2 affected expression of 165 transcripts in males, with 88 downregulated and 77 upregulated, while in females, 120 transcripts were affected with 62 downregulated and 58 upregulated. Finally, RT-qPCR validation demonstrated that expression of transcripts related to transposable elements, neuroserpin and heat shock protein were significantly affected by EE2-exposure. Our study is the first to report brain transcriptome libraries for guppies exposed to EE2. Not only does our study provide a valuable resource, it offers insights into the mechanisms underlying the feminizing effects on the brains of organisms exposed to environmentally realistic concentrations of EE2.

Sex Change in Clownfish: Molecular Insights from Transcriptome Analysis

Sequential hermaphroditism is a unique reproductive strategy among teleosts that is displayed mainly in fish species living in the coral reef environment. The reproductive biology of hermaphrodites has long been intriguing; however, very little is known about the molecular pathways underlying their sex change. Here, we provide the first de novo transcriptome analyses of a hermaphrodite teleost´s undergoing sex change in its natural environment. Our study has examined relative gene expression across multiple groups-rather than just two contrasting conditions- and has allowed us to explore the differential expression patterns throughout the whole process. Our analysis has highlighted the rapid and complex genomic response of the brain associated with sex change, which is subsequently transmitted to the gonads, identifying a large number of candidate genes, some well-known and some novel, involved in the process. The present study provides strong evidence of the importance of the sex steroidogenic machinery during sex change in clownfish, with the aromatase gene playing a central role, both in the brain and the gonad. This work constitutes the first genome-wide study in a social sex-changing species and provides insights into the genetic mechanism governing social sex change and gonadal restructuring in protandrous hermaphrodites.

Transcriptome analysis reveals differentially expressed genes associated with germ cell and gonad development in the Southern bluefin tuna (Thunnus maccoyii)

BACKGROUND

Controlling and managing the breeding of bluefin tuna (Thunnus spp.) in captivity is an imperative step towards obtaining a sustainable supply of these fish in aquaculture production systems. Germ cell transplantation (GCT) is an innovative technology for the production of inter-species surrogates, by transplanting undifferentiated germ cells derived from a donor species into larvae of a host species. The transplanted surrogates will then grow and mature to produce donor-derived seed, thus providing a simpler alternative to maintaining large-bodied broodstock such as the bluefin tuna. Implementation of GCT for new species requires the development of molecular tools to follow the fate of the transplanted germ cells. These tools are based on key reproductive and germ cell-specific genes. RNA-Sequencing (RNA-Seq) provides a rapid, cost-effective method for high throughput gene identification in non-model species. This study utilized RNA-Seq to identify key genes expressed in the gonads of Southern bluefin tuna (Thunnus maccoyii, SBT) and their specific expression patterns in male and female gonad cells.

RESULTS

Key genes involved in the reproductive molecular pathway and specifically, germ cell development in gonads, were identified using analysis of RNA-Seq transcriptomes of male and female SBT gonad cells. Expression profiles of transcripts from ovary and testis cells were compared, as well as testis germ cell-enriched fraction prepared with Percoll gradient, as used in GCT studies. Ovary cells demonstrated over-expression of genes related to stem cell maintenance, while in testis cells, transcripts encoding for reproduction-associated receptors, sex steroids and hormone synthesis and signaling genes were over-expressed. Within the testis cells, the Percoll-enriched fraction showed over-expression of genes that are related to post-meiosis germ cell populations.

CONCLUSIONS

Gonad development and germ cell related genes were identified from SBT gonads and their expression patterns in ovary and testis cells were determined. These expression patterns correlate with the reproductive developmental stage of the sampled fish. The majority of the genes described in this study were sequenced for the first time in T. maccoyii. The wealth of SBT gonadal and germ cell-related gene sequences made publicly available by this study provides an extensive resource for further GCT and reproductive molecular biology studies of this commercially valuable fish.

Comprehensive transcriptomic analysis of the process of gonadal sex differentiation in the turbot (Scophthalmus maximus)

The turbot is a flatfish with a ZW/ZZ sex determination system but with a still unknown sex determining gene(s), and with a marked sexual growth dimorphism in favor of females. To better understand sexual development in turbot we sampled young turbot encompassing the whole process of gonadal differentiation and conducted a comprehensive transcriptomic study on its sex differentiation using a validated custom oligomicroarray. Also, the expression profiles of 18 canonical reproduction-related genes were studied along gonad development. The expression levels of gonadal aromatase cyp19a1a alone at three months of age allowed the accurate and early identification of sex before the first signs of histological differentiation. A total of 56 differentially expressed genes (DEG) that had not previously been related to sex differentiation in fish were identified within the first three months of age, of which 44 were associated with ovarian differentiation (e.g., cd98, gpd1 and cry2), and 12 with testicular differentiation (e.g., ace, capn8 and nxph1). To identify putative sex determining genes, ∼4.000 DEG in juvenile gonads were mapped and their positions compared with that of previously identified sex- and growth-related quantitative trait loci (QTL). Although no genes mapped to the previously identified sex-related QTLs, two genes (foxl2 and 17βhsd) of the canonical reproduction-related genes mapped to growth-QTLs in linkage group (LG) 15 and LG6, respectively, suggesting that these genes are related to the growth dimorphism in this species.

Defining the role of omics in assessing ecosystem health: Perspectives from the Canadian environmental monitoring program

Scientific reviews and studies continue to describe omics technologies as the next generation of tools for environmental monitoring, while cautioning that there are limitations and obstacles to overcome. However, omics has not yet transitioned into national environmental monitoring programs designed to assess ecosystem health. Using the example of the Canadian Environmental Effects Monitoring (EEM) program, the authors describe the steps that would be required for omics technologies to be included in such an established program. These steps include baseline collection of omics endpoints across different species and sites to generate a range of what is biologically normal within a particular ecosystem. Natural individual variability in the omes is not adequately characterized and is often not measured in the field, but is a key component to an environmental monitoring program, to determine the critical effect size or action threshold for management. Omics endpoints must develop a level of standardization, consistency, and rigor that will allow interpretation of the relevance of changes across broader scales. To date, population-level consequences of routinely measured endpoints such as reduced gonad size or intersex in fish is not entirely clear, and the significance of genome-wide molecular, proteome, or metabolic changes on organism or population health is further removed from the levels of ecological change traditionally managed. The present review is not intended to dismiss the idea that omics will play a future role in large-scale environmental monitoring studies, but rather outlines the necessary actions for its inclusion in regulatory monitoring programs focused on assessing ecosystem health.

Development and validation of a mixed-tissue oligonucleotide DNA microarray for Atlantic bluefin tuna, Thunnus thynnus (Linnaeus, 1758)

BACKGROUND

The largest of the tuna species, Atlantic bluefin tuna (Thunnus thynnus), inhabits the North Atlantic Ocean and the Mediterranean Sea and is considered to be an endangered species, largely a consequence of overfishing. T. thynnus aquaculture, referred to as fattening or farming, is a capture based activity dependent on yearly renewal from the wild. Thus, the development of aquaculture practices independent of wild resources can provide an important contribution towards ensuring security and sustainability of this species in the longer-term. The development of such practices is today greatly assisted by large scale transcriptomic studies.

RESULTS

We have used pyrosequencing technology to sequence a mixed-tissue normalised cDNA library, derived from adult T. thynnus. A total of 976,904 raw sequence reads were assembled into 33,105 unique transcripts having a mean length of 893 bases and an N50 of 870. Of these, 33.4% showed similarity to known proteins or gene transcripts and 86.6% of them were matched to the congeneric Pacific bluefin tuna (Thunnus orientalis) genome, compared to 70.3% for the more distantly related Nile tilapia (Oreochromis niloticus) genome. Transcript sequences were used to develop a novel 15 K Agilent oligonucleotide DNA microarray for T. thynnus and comparative tissue gene expression profiles were inferred for gill, heart, liver, ovaries and testes. Functional contrasts were strongest between gills and ovaries. Gills were particularly associated with immune system, signal transduction and cell communication, while ovaries displayed signatures of glycan biosynthesis, nucleotide metabolism, transcription, translation, replication and repair.

CONCLUSIONS

Sequence data generated from a novel mixed-tissue T. thynnus cDNA library provide an important transcriptomic resource that can be further employed for study of various aspects of T. thynnus ecology and genomics, with strong applications in aquaculture. Tissue-specific gene expression profiles inferred through the use of novel oligo-microarray can serve in the design of new and more focused transcriptomic studies for future research of tuna physiology and assessment of the welfare in a production environment.

Predose and Postdose Blood Gene Expression Profiles Identify the Individuals Susceptible to Acetaminophen-Induced Liver Injury in Rats

The extent of drug-induced liver injury (DILI) can vary greatly between different individuals. Thus, it is crucial to identify susceptible population to DILI. The aim of this study was to determine whether transcriptomics analysis of predose and postdose rat blood would allow prediction of susceptible individuals to DILI using the widely applied analgesic acetaminophen (APAP) as a model drug. Based on ranking in alanine aminotransferase levels, five most susceptible and five most resistant rats were identified as two sub-groups after APAP treatment. Predose and postdose gene expression profiles of blood samples from these rats were determined by microarray analysis. The expression of 158 genes innately differed in the susceptible rats from the resistant rats in predose data. In order to identify more reliable biomarkers related to drug responses for detecting individuals susceptibility to APAP-induced liver injury (AILI), the changes of these genes' expression posterior to APAP treatment were detected. Through the further screening method based on the trends of gene expression between the two sub-groups before and after drug treatment, 10 genes were identified as potential predose biomarkers to distinguish between the susceptible and resistant rats. Among them, four genes, Incenp, Rpgrip1, Sbf1, and Mmp12, were found to be reproducibly in real-time PCR with an independent set of animals. They were all innately higher expressed in resistant rats to AILI, which are closely related to cell proliferation and tissue repair functions. It indicated that rats with higher ability of cell proliferation and tissue repair prior to drug treatment might be more resistant to AILI. In this study, we demonstrated that combination of predose and postdose gene expression profiles in blood might identify the drug related inter-individual variation in DILI, which is a novel and important methodology for identifying susceptible population to DILI.

A functional genomics tool for the Pacific bluefin tuna: Development of a 44K oligonucleotide microarray from whole-genome sequencing data for global transcriptome analysis

Bluefin tunas are one of the most important fishery resources worldwide. Because of high market values, bluefin tuna farming has been rapidly growing during recent years. At present, the most common form of the tuna farming is based on the stocking of wild-caught fish. Therefore, concerns have been raised about the negative impact of the tuna farming on wild stocks. Recently, the Pacific bluefin tuna (PBT), Thunnus orientalis, has succeeded in completing the reproduction cycle under aquaculture conditions, but production bottlenecks remain to be solved because of very little biological information on bluefin tunas. Functional genomics approaches promise to rapidly increase our knowledge on biological processes in the bluefin tuna. Here, we describe the development of the first 44K PBT oligonucleotide microarray (oligo-array), based on whole-genome shotgun (WGS) sequencing and large-scale expressed sequence tags (ESTs) data. In addition, we also introduce an initial 44K PBT oligo-array experiment using in vitro grown peripheral blood leukocytes (PBLs) stimulated with immunostimulants such as lipopolysaccharide (LPS: a cell wall component of Gram-negative bacteria) or polyinosinic:polycytidylic acid (poly I:C: a synthetic mimic of viral infection). This pilot 44K PBT oligo-array analysis successfully addressed distinct immune processes between LPS- and poly I:C- stimulated PBLs. Thus, we expect that this oligo-array will provide an excellent opportunity to analyze global gene expression profiles for a better understanding of diseases and stress, as well as for reproduction, development and influence of nutrition on tuna aquaculture production.

Expression Comparison of Oil Biosynthesis Genes in Oil Palm Mesocarp Tissue Using Custom Array

Gene expression changes that occur during mesocarp development are a major research focus in oil palm research due to the economic importance of this tissue and the relatively rapid increase in lipid content to very high levels at fruit ripeness. Here, we report the development of a transcriptome-based 105,000-probe oil palm mesocarp microarray. The expression of genes involved in fatty acid (FA) and triacylglycerol (TAG) assembly, along with the tricarboxylic acid cycle (TCA) and glycolysis pathway at 16 Weeks After Anthesis (WAA) exhibited significantly higher signals compared to those obtained from a cross-species hybridization to the Arabidopsis (p-value < 0.01), and rice (p-value < 0.01) arrays. The oil palm microarray data also showed comparable correlation of expression (r² = 0.569, p < 0.01) throughout mesocarp development to transcriptome (RNA sequencing) data, and improved correlation over quantitative real-time PCR (qPCR) (r² = 0.721, p < 0.01) of the same RNA samples. The results confirm the advantage of the custom microarray over commercially available arrays derived from model species. We demonstrate the utility of this custom microarray to gain a better understanding of gene expression patterns in the oil palm mesocarp that may lead to increasing future oil yield.

Genetic architecture of sex determination in fish: applications to sex ratio control in aquaculture

Controlling the sex ratio is essential in finfish farming. A balanced sex ratio is usually good for broodstock management, since it enables to develop appropriate breeding schemes. However, in some species the production of monosex populations is desirable because the existence of sexual dimorphism, primarily in growth or first time of sexual maturation, but also in color or shape, can render one sex more valuable. The knowledge of the genetic architecture of sex determination (SD) is convenient for controlling sex ratio and for the implementation of breeding programs. Unlike mammals and birds, which show highly conserved master genes that control a conserved genetic network responsible for gonad differentiation (GD), a huge diversity of SD mechanisms has been reported in fish. Despite theory predictions, more than one gene is in many cases involved in fish SD and genetic differences have been observed in the GD network. Environmental factors also play a relevant role and epigenetic mechanisms are becoming increasingly recognized for the establishment and maintenance of the GD pathways. Although major genetic factors are frequently involved in fish SD, these observations strongly suggest that SD in this group resembles a complex trait. Accordingly, the application of quantitative genetics combined with genomic tools is desirable to address its study and in fact, when applied, it has frequently demonstrated a multigene trait interacting with environmental factors in model and cultured fish species. This scenario has notable implications for aquaculture and, depending upon the species, from chromosome manipulation or environmental control techniques up to classical selection or marker assisted selection programs, are being applied. In this review, we selected four relevant species or fish groups to illustrate this diversity and hence the technologies that can be used by the industry for the control of sex ratio: turbot and European sea bass, two reference species of the European aquaculture, and salmonids and tilapia, representing the fish for which there are well established breeding programs.

Effects of temperature acclimation on Pacific bluefin tuna (Thunnus orientalis) cardiac transcriptome

Little is known about the mechanisms underpinning thermal plasticity of vertebrate hearts. Bluefin tuna hearts offer a unique model to investigate processes underlying thermal acclimation. Their hearts, while supporting an endothermic physiology, operate at ambient temperature, and are presented with a thermal challenge when migrating to different thermal regimes. Here, we examined the molecular responses in atrial and ventricular tissues of Pacific bluefin tuna acclimated to 14°C, 20°C, and 25°C. Quantitative PCR studies showed an increase in sarcoplasmic reticulum Ca2+ ATPase gene expression with cold acclimation and an induction of Na+/Ca2+-exchanger gene at both cold and warm temperatures. These data provide evidence for thermal plasticity of excitation-contraction coupling gene expression in bluefin tunas and indicate an increased capacity for internal Ca2+ storage in cardiac myocytes at 14°C. Transcriptomic analysis showed profound changes in cardiac tissues with acclimation. A principal component analysis revealed that temperature effect was greatest on gene expression in warm-acclimated atrium. Overall data showed an increase in cardiac energy metabolism at 14°C, potentially compensating for cold temperature to optimize bluefin tuna performance in colder oceans. In contrast, metabolic enzyme activity and gene expression data suggest a decrease in ATP production at 25°C. Expression of genes involved in protein turnover and molecular chaperones was also decreased at 25°C. Expression of genes involved in oxidative stress response and programmed cell death suggest an increase in oxidative damage and apoptosis at 25°C, particularly in the atrium. These findings provide insights into molecular processes that may characterize cardiac phenotypes at upper thermal limits of teleosts.