Identification of a novel allele of peroxisome proliferator-activated receptor gamma (PPARG) and its association with resistance to Aeromonas salmonicida in Atlantic salmon (Salmo salar)

The Peroxisome Proliferator-Activated Receptors of Ray-Finned Fish: Unique Structures, Elusive Functions

The Actinopterygian and specifically the Teleostean peroxisome proliferator-activated receptors (PPARs) present an impressive variability and complexity in their structures, both at the gene and protein levels. These structural differences may also reflect functional divergence from their mammalian homologs, or even between fish species. This review, taking advantage of the data generated from the whole-genome sequencing of several fish species, highlights the differences in the primary structure of the receptors, while discussing results from the literature pertaining to the functions of fish PPARs and their activation by natural and synthetic compounds.

Dietary Exposure to Environmentally Relevant Levels of Chemical Contaminants Reduces Growth and Survival in Juvenile Chinook Salmon

Chemical pollution can degrade aquatic ecosystems. Chinook salmon in contaminated habitats are vulnerable to health impacts from toxic exposures. Few studies have been conducted on adverse health outcomes associated with current levels and mixtures of contaminants. Fewer still address effects specific to the juvenile life-stage of salmonids. The present study evaluated contaminant-related effects from dietary exposure to environmentally relevant concentrations and mixture profiles in juvenile Chinook salmon from industrialized waterways in the U.S. Pacific Northwest using two end points: growth assessment and disease susceptibility. The dose and chemical proportions were reconstituted based on environmental sampling and analysis using the stomach contents of juvenile Chinook salmon recently collected from contaminated, industrialized waterways. Groups of fish were fed a mixture with fixed proportions of 10 polychlorinated biphenyls (PCBs), 3 dichlorodiphenyltrichloroethanes (DDTs), and 13 polycyclic aromatic hydrocarbons (PAHs) at five concentrations for 35 days. These contaminant compounds were selected because of elevated concentrations and the widespread presence in sediments throughout industrialized waterways. Fork length and otolith microstructural growth indicators were significantly reduced in fish fed environmentally relevant concentrations of these contaminants. In addition, contaminant-exposed Chinook salmon were more susceptible to disease during controlled challenges with the pathogen Aeromonas salmonicida. Our results indicate that dietary exposure to contaminants impairs growth and immune function in juvenile Chinook salmon, thereby highlighting that current environmental exposure to chemicals of potential management concern threatens the viability of exposed salmon.

Gill and Liver Transcript Expression Changes Associated With Gill Damage in Atlantic Salmon (Salmo salar)

Gill damage represents a significant challenge in the teleost fish aquaculture industry globally, due to the gill’s involvement in several vital functions and direct contact with the surrounding environment. To examine the local and systemic effects accompanying gill damage (which is likely to negatively affect gill function) of Atlantic salmon, we performed a field sampling to collect gill and liver tissue after several environmental insults (e.g., harmful algal blooms). Before sampling, gills were visually inspected and gill damage was scored; gill scores were assigned from pristine [gill score 0 (GS0)] to severely damaged gills (GS3). Using a 44K salmonid microarray platform, we aimed to compare the transcriptomes of pristine and moderately damaged (i.e., GS2) gill tissue. Rank Products analysis (5% percentage of false-positives) identified 254 and 34 upregulated and downregulated probes, respectively, in GS2 compared with GS0. Differentially expressed probes represented genes associated with functions including gill remodeling, wound healing, and stress and immune responses. We performed gill and liver qPCR for all four gill damage scores using microarray-identified and other damage-associated biomarker genes. Transcripts related to wound healing (e.g., neb and klhl41b) were significantly upregulated in GS2 compared with GS0 in the gills. Also, transcripts associated with immune and stress-relevant pathways were dysregulated (e.g., downregulation of snaclec 1-like and upregulation of igkv3) in GS2 compared with GS0 gills. The livers of salmon with moderate gill damage (i.e., GS2) showed significant upregulation of transcripts related to wound healing (i.e., chtop), apoptosis (e.g., bnip3l), blood coagulation (e.g., f2 and serpind1b), transcription regulation (i.e., pparg), and stress-responses (e.g., cyp3a27) compared with livers of GS0 fish. We performed principal component analysis (PCA) using transcript levels for gill and liver separately. The gill PCA showed that PC1 significantly separated GS2 from all other gill scores. The genes contributing most to this separation were pgam2, des, neb, tnnt2, and myom1. The liver PCA showed that PC1 significantly separated GS2 from GS0; levels of hsp70, cyp3a27, pparg, chtop, and serpind1b were the highest contributors to this separation. Also, hepatic acute phase biomarkers (e.g., serpind1b and f2) were positively correlated to each other and to gill damage. Gill damage-responsive biomarker genes and associated qPCR assays arising from this study will be valuable in future research aimed at developing therapeutic diets to improve farmed salmon welfare.

The Preservation of PPARγ Genome Duplicates in Some Teleost Lineages: Insights into Lipid Metabolism and Xenobiotic Exploitation

Three peroxisome proliferator-activated receptor paralogues (PPARα, -β and -γ) are currently recognized in vertebrate genomes. PPARγ is known to modulate nutrition, adipogenesis and immunity in vertebrates. Natural ligands of PPARγ have been proposed; however, the receptor also binds synthetic ligands such as endocrine disruptors. Two paralogues of PPARα and PPARβ have been documented in teleost species, a consequence of the 3R WGD. Recently, two PPARγ paralogue genes were also identified in Astyanax mexicanus. We aimed to determine whether the presence of two PPARγ paralogues is prevalent in other teleost genomes, through genomic and phylogenetic analysis. Our results showed that besides Characiformes, two PPARγ paralogous genes were also identified in other teleost taxa, coinciding with the teleost-specific, whole-genome duplication and with the retention of both genes prior to the separation of the Clupeocephala. To functionally characterize these genes, we used the European sardine (Sardina pilchardus) as a model. PPARγA and PPARγB display a different tissue distribution, despite the similarity of their functional profiles: they are unresponsive to tested fatty acids and other human PPARγ ligands yet yield a transcriptional response in the presence of tributyltin (TBT). This observation puts forward the relevance of comparative analysis to decipher alternative binding architectures and broadens the disruptive potential of man-made chemicals for aquatic species.

Vaccine-Induced Protection Against Furunculosis Involves Pre-emptive Priming of Humoral Immunity in Arctic Charr

With respect to salmonid aquaculture, one of the most important bacterial pathogens due to high mortality and antibiotic usage is the causative agent of typical furunculosis, Aeromonas salmonicida spp. salmonicida (Asal). In Atlantic salmon, Salmo salar, the host response during infections with Asal is well-documented, with furunculosis outbreaks resulting in significant mortality in commercial settings. However, less is known about the host-pathogen interactions in the emerging aquaculture species, Arctic charr Salvelinus alpinus. Furthermore, there is no data on the efficacy or response of this species after vaccination with commonly administered vaccines against furunculosis. To this end, we examined the immunological response of S. alpinus during infection with Asal, with or without administration of vaccines (Forte Micro®, Forte Micro® + Renogen®, Elanco Animal Health). Artic charr (vaccinated or unvaccinated) were i.p.-injected with a virulent strain of Asal (10⁶ CFUs/mL) and tissues were collected pre-infection/post-vaccination, 8, and 29 days post-infection. Unvaccinated Arctic charr were susceptible to Asal with 72% mortalities observed after 31 days. However, there was 72–82% protection in fish vaccinated with either the single or dual-vaccine, respectively. Protection in vaccinated fish was concordant with significantly higher serum IgM concentrations, and following RNA sequencing and transcriptome assembly, differential expression analysis revealed several patterns and pathways associated with the improved survival of vaccinated fish. Most striking was the dramatically higher basal expression of complement/coagulation factors, acute phase-proteins, and iron hemostasis proteins in pre-challenged, vaccinated fish. Remarkably, following Asal infection, this response was abrogated and instead the transcriptome was characterized by a lack of immune-stimulation compared to that of unvaccinated fish. Furthermore, where pathways of actin assembly and FcγR-mediated phagocytosis were significantly differentially regulated in unvaccinated fish, vaccinated fish showed either the opposite regulation (ForteMicro®), or no impact at all (ForteMicro®Renogen®). The present data indicates that vaccine-induced protection against Asal relies on the pre-activation and immediate control of humoral immune parameters that is coincident with reduced activation of apoptotic (e.g., NF-κB) and actin-associated pathways.

Transcriptomics analysis reveals candidate genes and pathways for susceptibility or resistance to Singapore grouper iridovirus in orange-spotted grouper (Epinephelus coioides)

In this study, the transcriptional response of grouper to Singapore grouper iridovirus (SGIV) stimulation was characterized using RNA sequencing. Transcriptome sequencing of three test groups in the grouper was performed using the Illumina MiSeq platform. The three test groups were a control group, which was injected with PBS buffer; a high-susceptible (HS) group, which died shortly after the SGIV injection; and a high-resistance (HR) group, which survived the SGIV injection. In total, 38,253 unigenes were generated. When the HS group was compared with the control group, 885 unigenes were upregulated and 487 unigenes were downregulated. When the HR and control groups were compared, 1114 unigenes were upregulated and 420 were downregulated, and when the HR and HS groups were compared, 1010 unigenes were upregulated and 375 were downregulated. In the KEGG analysis, two immune-related pathways, the p53 and peroxisome proliferator-activated receptor pathways, were detected with highly significant enrichment. In addition, 7465 microsatellites and 22,1569 candidate single nucleotide polymorphisms were identified from our transcriptome data. The results suggested several pathways that are associated with traits of disease susceptibility or disease resistance, and provided extensive information about novel gene sequences, gene expression profiles, and genetic markers. This may contribute to vaccine research and a breeding program against SGIV infection in grouper.

Cloning and expression characterization of peroxisome proliferator-activated receptors (PPARs) with their agonists, dietary lipids, and ambient salinity in rabbitfish Siganus canaliculatus

Rabbitfish Siganus canaliculatus is the first marine teleost reported to have the ability of biosynthesizing C_20-22 long-chain polyunsaturated fatty acids (LC-PUFA) from C₁₈ precursors, and thus provides a model for studying the regulatory mechanisms of LC-PUFA biosynthesis in teleosts. To investigate the possible roles of peroxisome proliferator-activated receptors (PPARs), critical transcription factors involved in the regulation of lipid metabolism, in the regulation of LC-PUFA biosynthesis in rabbitfish, the PPAR genes were cloned and their expression characterization with PPAR agonists, dietary lipid resource, and ambient salinity were examined. Three cDNA sequences respectively encoding 477, 516 and 519 amino acids of PPARα, PPARβ, and PPARγ isoforms were obtained. PPARα exhibited a wide tissue expression with its highest levels in the heart and brain; PPARβ was predominantly expressed in the gills, while PPARγ was highly expressed in the intestine and gills. In rabbitfish primary hepatocytes, both the PPAR agonists 2-bromopalmitate (2-Bro) and fenofibrate (FF) increased the expression of PPARγ, SREBP1c and Elovl5, whereas FF depressed the expression of Δ6/Δ5 Fad. Moreover, a higher hepatic PPARβ expression was observed in fish fed diets with vegetable oils (VO) than that with fish oil (FO), in the former the expression of PPARα, PPARβ, and PPARγ were increased at the low ambient salinity (10ppt), where an increasing expression of Δ5/Δ6 Fad, Δ4 Fad and Elovl5 genes was previously reported. These results suggest that PPARs might be involved in the upregulation of LC-PUFA biosynthesis with dietary VO and low ambient salinity in rabbitfish.

The Role of Peroxisome Proliferator-Activated Receptor Gamma (PPARG) in Adipogenesis: Applying Knowledge from the Fish Aquaculture Industry to Biomedical Research

The tropical freshwater zebrafish has recently emerged as a valuable model organism for the study of adipose tissue biology and obesity-related disease. The strengths of the zebrafish model system are its wealth of genetic mutants, transgenic tools, and amenability to high-resolution imaging of cell dynamics within live animals. However, zebrafish adipose research is at a nascent stage and many gaps exist in our understanding of zebrafish adipose physiology and metabolism. By contrast, adipose research within other, closely related, teleost species has a rich and extensive history, owing to the economic importance of these fish as a food source. Here, we compare and contrast knowledge on peroxisome proliferator-activated receptor gamma (PPARG)-mediated adipogenesis derived from both biomedical and aquaculture literatures. We first concentrate on the biomedical literature to (i) briefly review PPARG-mediated adipogenesis in mammals, before (ii) reviewing Pparg-mediated adipogenesis in zebrafish. Finally, we (iii) mine the aquaculture literature to compare and contrast Pparg-mediated adipogenesis in aquaculturally relevant teleosts. Our goal is to highlight evolutionary similarities and differences in adipose biology that will inform our understanding of the role of adipose tissue in obesity and related disease.

Peroxisome proliferator-activated receptor gamma (PPARγ) in brown trout: Interference of estrogenic and androgenic inputs in primary hepatocytes

Peroxisome proliferator-activated receptor gamma (PPARγ) is a pivotal regulator of lipid and glucose metabolism in vertebrates. Here, we isolated and characterized for the first time the PPARγ gene from brown trout (Salmo trutta f. fario). Hormones have been reported to interfere with the regulatory function of PPARγ in various organisms, albeit with little focus on fish. Thus, primary hepatocytes isolated from juveniles of brown trout were exposed to 1, 10 and 50μM of ethinylestradiol (EE2) or testosterone (T). A significant (3 fold) decrease was obtained in response to 50μM of EE2 and to 10 and 50μM of T (13 and 14 folds), while a 3 fold increase was observed at 1μM of EE2. Therefore, trout PPARγ seems a target for natural/synthetic compounds with estrogenic or androgenic properties and so, we advocate considering PPARγ as another alert sensor gene when assessing the effects of sex-steroid endocrine disruptors.

Characterization of two novel gadd45a genes in hybrid tilapia and their responses to the infection of Streptococcus agalactiae

Diseases are one of the major challenges in tilapia aquaculture. Identification of DNA markers associated with disease resistance may facilitate the acceleration of the selection for disease resistance. Gadd45a (growth arrest and DNA damage 45 A), a stress-inducible gene in humans and mice, has not been studied in fish. We characterized the two prologues of Gadd45a genes in hybrid tilapia. Gadd45a1 and Gadd45a2 shared an identical gene structure and showed an amino acid sequence identity of 73.8%. Their expressions were detected in all 10 tissues examined, with the kidney and gill having high transcriptional expressions. The expression levels of Gadd45a1 were significantly lower than those of Gadd45a2 in all examined tissues. After a challenge with a bacterial pathogen Streptococcus agalactiae, the expressions of the two genes were up-regulated significantly in the spleen, kidney, liver and intestine. These findings suggest that the two Gadd45a genes play an important role in the resistance to S. agalactiae in tilapia. We identified 10 SNPs in the two genes. The SNP markers in the two Gadd45a genes could be used to examine whether they are associated with resistance to S. agalactiae.

Characterization of the duodenase-1 gene and its associations with resistance to Streptococuus agalactiae in hybrid tilapia (Oreochromis spp.)

Tilapia is a group of cultured teleost fishes whose production is threatened by some diseases. Identification of DNA markers associated with disease resistance in candidate genes may facilitate to accelerate the selection of disease resistance. The gene encoding a duodenase, which can trigger immune response, has not been studied in fish. We characterized the cDNA of duodenase-1 gene of hybrid tilapia. Its ORF is 759 bp, encoding a serine protease of 252 amino acids. This gene consisted of five exons and four introns. Its expression was detected in all 10 tissues examined, and it was highly expressed in the intestine and kidney. After a challenge with the bacterial pathogen, Streptococcus agalactiae, its expression was up-regulated significantly in the intestine, liver and spleen. We identified seven SNPs in the gene and found that four of them were significantly associated with the resistance to S. agalactiae (P < 0.05). The CGTCC haplotype, CAGTC/CGGTC and CGTCC/CGTCC diplotype were significantly associated with the resistance to S. agalactiae (P = 0.00, 0.04 and < 0.0001, respectively). In addition, one SNP was associated significantly with growth traits (P < 0.05). These findings suggest that the duodenase-1 gene plays an important role in the resistance to S. agalactiae in tilapia. The SNP markers in the duodenase-1 gene associated with resistance to the bacterial pathogen, may facilitate the selection of tilapia resistant to the bacterial disease.

Peroxisome proliferator-activated receptor gamma (PPARγ) in yellow catfish Pelteobagrus fulvidraco: molecular characterization, mRNA expression and transcriptional regulation by insulin in vivo and in vitro

Peroxisome proliferator-activated receptor gamma (PPARγ) is ligand-inducible transcription factor and has important roles in lipid metabolism, cell proliferation and inflammation. In the present study, yellow catfish Pelteobagrus fulvidraco PPARγ cDNA was isolated from liver by RT-PCR and RACE, and its molecular characterization and transcriptional regulation by insulin in vivo and in vitro were determined. The generation of PPARγ1 and PPARγ2 was due to alternative promoter of PPARγ gene. PPARγ1 and PPARγ2 mRNA covered 2426 bp and 2537 bp, respectively, with an open reading frame (ORF) of 1584 bp encoding 527 amino acid residues. Yellow catfish PPARγ gene was organized in a manner similar to that of their mammalian homologs, implying a modular organization of the protein's domains. A comparison between the yellow catfish PPARγ amino acid sequence and the correspondent sequences of several other species revealed the identity of 55-76.2%. Two PPARγ transcripts (PPARγ1 and PPARγ2) mRNAs were expressed in a wide range of tissues, but the abundance of each PPARγ mRNA showed the tissue- and developmental stage-dependent expression patterns. Intraperitoneal injection of insulin in vivo significantly stimulated the mRNA expression of total PPARγ and PPARγ1, but not PPARγ2 in the liver of yellow catfish. In contrast, incubation of hepatocytes with insulin in vitro increased the mRNA levels of PPARγ1, PPARγ2 and total PPARγ. To our knowledge, for the first time, the present study provides evidence that PPARγ1 and PPARγ2 are differentially expressed with and among tissues during different developmental stages and also regulated by insulin both in vivo and in vitro, which serves to increase our understanding on PPARγ physiological function in fish.

Altered gene expression patterns of innate and adaptive immunity pathways in transgenic rainbow trout harboring Cecropin P1 transgene

Background

We have recently developed several homozygous families of transgenic rainbow trout harbouring cecropin P1 transgene. These fish exhibit resistance characteristic to infection by Aeromonas salmonicida and infectious hematopoietic necrosis virus (IHNV). In our earlier studies we have reported that treatment of a rainbow trout macrophage cell line (RTS11) with a linear cationic α-helical antimicrobial peptide (e.g., cecropin B) resulted in elevated levels of expression of two pro-inflammatory relevant genes (e.g., IL-1β and COX-2). Therefore, we hypothesized that in addition to the direct antimicrobial activity of cecropin P1 in the disease resistant transgenic rainbow trout, this antimicrobial peptide may also affect the expression of immune relevant genes in the host. To confirm this hypothesis, we launched a study to determine the global gene expression profiles in three immune competent organs of cecropin P1 transgenic rainbow trout by using a 44k salmonid microarray.

Results

From the microarray data, a total of 2480 genes in the spleen, 3022 in the kidney, and 2102 in the liver were determined as differentially expressed genes (DEGs) in the cecropin P1 transgenic rainbow trout when compared to the non-transgenics. There were 478 DEGs in common among three tissues. Enrichment analyses conducted by two different bioinformatics tools revealed a tissue specific profile of functional pathway perturbation. Many of them were directly related to innate immune system such as phagocytosis, lysosomal processing, complement activation, antigen processing/presentation, and leukocyte migration. Perturbation of other biological functions that might contribute indirectly to host immunity was also observed.

Conclusions

The gene product of cecropin P1 transgene produced in the disease resistant transgenic rainbow trout not only can kill the pathogens directly but also exert multifaceted immunomodulatory properties to boost host immunity. The identified genes involved in different pathways related to immune function are valuable indicators associated with enhanced host immunity. These genes may serve as markers for selective breeding of rainbow trout or other aquaculture important fish species bearing traits of disease resistance.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-887) contains supplementary material, which is available to authorized users.

Gene expression profiles of the spleen, liver, and head kidney in turbot (Scophthalmus maximus) along the infection process with Aeromonas salmonicida using an immune-enriched oligo-microarray

We evaluated the expression profiles of turbot in the spleen, liver, and head kidney across five temporal points of the Aeromonas salmonicida infection process using an 8 × 15 K Agilent oligo-microarray. The microarray included 2,176 different fivefold replicated gene probes designed from a turbot 3' sequenced EST database. We were able to identify 471 differentially expressed (DE) genes (17.3% of the whole microarray), 223 in the spleen, 246 in the liver, and 125 in the head kidney, in at least one of the five temporal points sampled for each organ. Most of these genes could be annotated (83.0%) and functionally categorized using Gene Ontology terms (69.1%) after the additional sequencing of DE genes from the 5' end. Many DE genes were related to innate and acquired immune functions in accordance to previous studies with this pathogen in other fish species. A high proportion of DE genes were organ specific (77.1%), but their associated GO functions were rather similar in the three organs. The most striking difference in functional distribution was observed between the up- and down-regulated gene groups. Up-regulated genes were mostly associated to key immune functions while down-regulated ones mainly involved metabolism- and transport-related genes. Genetic response appeared clustered in groups of genes with similar expression profiles along the temporal series. The spleen showed the most clustering while the liver and head kidney displayed a higher diversification. The information obtained will aid to understand the turbot immune response and will specifically be valuable to develop strategies of defense to A. salmonicida to achieve more resistant broodstocks for turbot industry.