Genomic characterization and expression analysis of four apolipoprotein A-IV paralogs in Senegalese sole (Solea senegalensis Kaup)

Establishment of a lipid metabolism disorder model in ApoEb mutant zebrafish

BACKGROUND AND AIMS

ApoEb is a zebrafish homologous to mammalian ApoE, whose deficiency would lead to lipid metabolism disorders (LMDs) like atherosclerosis. We attempted to knock out the zebrafish ApoEb, then establish a zebrafish model with LMD.

METHODS

ApoEb was knocked out using the CRISPR/Cas9 system, and the accumulation of lipids was confirmed by Oil Red O staining, confocal imaging, and lipid measurements. The lipid-lowering effects of simvastatin (SIM), ezetimibe (EZE) and Xuezhikang (XZK), an extract derived from red yeast rice, were evaluated through in vivo imaging in zebrafish larvae.

RESULTS

In the ApoEb mutant, significant vascular lipid deposition occurred, and lipid measurement performed in the whole-body homogenate of larvae and adult plasma showed significantly increased lipid levels. SIM, EZE and XZK apparently relieved hyperlipidemia in ApoEb mutants, and XZK had a significant inhibitory effect on the recruitment of neutrophils and macrophages.

CONCLUSIONS

In this study, an LMD model has been established in ApoEb mutant zebrafish. We suggest that this versatile model could be applied in studying hypercholesterolemia and related vascular pathology in the context of early atherosclerosis, as well as the physiological function of ApoE.

The effects of feeding and starvation on antioxidant defence, fatty acid composition and lipid peroxidation in reared Oncorhynchus mykiss fry

The effects of feeding and starvation have been studied with respect to oxidative stress and enzymatic antioxidant activities in the whole body of 4 cm rainbow trout fry Oncorhynchus mykiss (Walbaum 1792). The experiment was conducted for 28 days. The selected biomarkers for the study were determined, including non-enzymic scavengers glutathione (GSH), oxidized glutathione (GSSG) and malondialdehyde (MDA) contents and a number of enzymes are known to have major antioxidant activity, such as activities of süperoksit dismutaz (SOD), catalase (CAT), glutatyon peroksidaz (GSHpx), glutatyon Redüktaz (GR) and Glutatyon-S-Transferaz (GST). There is an endogenous cellular glutathione pool which consists of two forms of glutathione, i.e. the GSH and the GSSG. Oxidative damage was measured by the formation of MDA as an indication of lipid peroxidation. The activities of SOD in 14th and 28th day and the activity of CAT in 14th day were increased significantly during the 28 days of starvation. GSHpx and GR activities in starved fry decreased significantly in 28th day. GST activity in all starved fry showed the most significant increases the period of 28 days starving. The highest ΣSFA (Total Saturated Fatty Acid) content was obtained from 28 day starved fry. In starved fry, there was an apparent preference in utilization of C18:1n-9 than in the fed fry. In both starved and fed fry, C16:1n-7 was preferentially kept during the same period. Fry kept 28 days under starvation conditions exhausted C15:0, C17:0, C18:3n-6, C22:0, C24:0. They utilized less C20:5n-3 acid and conserved strongly C22:6n-3 acid. Concentrations of C20:5n-3, C22:5n-3, C22:6n-3 and total n-3 fatty acids significantly increased and C18:3n-3 significantly decreased in the whole body of starved fry during starvation period. A significant increase in the concentrations of C22:5n-3 and C22:6n-3 was determined in the fed fries in the last 2 weeks. Fat-soluble vitamins, cholesterol, stigmasterol and β-sitosterol levels were also determined in the same period of O. mykiss fry.

Re-evaluation of common carp (Cyprinus carpio L.) housekeeping genes for gene expression studies - considering duplicated genes

Quantitative real-time PCR is one of the most widely used techniques for measuring changes in the expression of target transcripts due to its sensitivity, specificity, and cost-effectiveness. However, the essential step that determines appropriate and correct data interpretation is the selection of proper endogenous control genes. Identifying useful reference genes with stable expression is critical for accurate normalization and precise results. Functional divergence of duplicated genes in tetraploid species, like common carp, can complicate the choice for a proper reference gene. In the present study, we determined the expression stability of duplicated genes of 40s, b2m, ef1α, gapdh, g6pd, and odc1 in different tissues of common carp (Cyprinus carpio L.). Gene expression analysis comprised healthy control fish, fish under bacterial and parasitic infections, and across the early stage of common carp development. Obtained data were compared with the actb gene, which is used widely as a reference in RT-qPCR analysis. The application of the three different algorithms - geNorm, NormFinder, BestKeeper, allowed comparative evaluation of the expression stability of the tested genes. Subsequently, the RefFinder, a web-based tool, was used to rank the examined housekeeping genes comprehensively. We demonstrate variable transcription stability levels in the examined mRNAs as well as differences in expression between paralog gene copies. The 40s, b2m, ef1α and actb genes showed the most stable expression across all physiological conditions and tissues. The gapdh, odc1, and g6pd gene variants demonstrated lower stability. Differences in expression patterns between duplicated genes underline the possibility of functional divergence between them. This aspect should be considered in polyploid species before selecting the reference gene(s). Our study also points on the importance of choice for a reference gene (paralog) when expressing newly identified genes and the spatial expression profile is performed. SUBJECTS: Aquaculture, Molecular Biology, Fish Science.

Development of whole-genome multiplex assays and construction of an integrated genetic map using SSR markers in Senegalese sole

The Senegalese sole (Solea senegalensis) is an economically important flatfish species. In this study, a genome draft was analyzed to identify microsatellite (SSR) markers for whole-genome genotyping. A subset of 224 contigs containing SSRs were preselected and validated by using a de novo female hybrid assembly. Overall, the SSR density in the genome was 886.7 markers per megabase of genomic sequences and the dinucleotide motif was the most abundant (52.4%). In silico comparison identified a set of 108 SSRs (with di-, tetra- or pentanucleotide motifs) widely distributed in the genome and suitable for primer design. A total of 106 markers were structured in thirteen multiplex PCR assays (with up to 10-plex) and the amplification conditions were optimized with a high-quality score. Main genetic diversity statistics and genotyping reliability were assessed. A subset of 40 high polymorphic markers were selected to optimize four supermultiplex PCRs (with up to 11-plex) for pedigree analysis. Theoretical exclusion probabilities and real parentage allocation tests using parent–offspring information confirmed their robustness and effectiveness for parental assignment. These new SSR markers were combined with previously published SSRs (in total 229 makers) to construct a new and improved integrated genetic map containing 21 linkage groups that matched with the expected number of chromosomes. Synteny analysis with respect to C. semilaevis provided new clues on chromosome evolution in flatfish and the formation of metacentric and submetacentric chromosomes in Senegalese sole.

An improved de novo assembling and polishing of Solea senegalensis transcriptome shed light on retinoic acid signalling in larvae

Senegalese sole is an economically important flatfish species in aquaculture and an attractive model to decipher the molecular mechanisms governing the severe transformations occurring during metamorphosis, where retinoic acid seems to play a key role in tissue remodeling. In this study, a robust sole transcriptome was envisaged by reducing the number of assembled libraries (27 out of 111 available), fine-tuning a new automated and reproducible set of workflows for de novo assembling based on several assemblers, and removing low confidence transcripts after mapping onto a sole female genome draft. From a total of 96 resulting assemblies, two "raw" transcriptomes, one containing only Illumina reads and another with Illumina and GS-FLX reads, were selected to provide SOLSEv5.0, the most informative transcriptome with low redundancy and devoid of most single-exon transcripts. It included both Illumina and GS-FLX reads and consisted of 51,348 transcripts of which 22,684 code for 17,429 different proteins described in databases, where 9527 were predicted as complete proteins. SOLSEv5.0 was used as reference for the study of retinoic acid (RA) signalling in sole larvae using drug treatments (DEAB, a RA synthesis blocker, and TTNPB, a RA-receptor agonist) for 24 and 48 h. Differential expression and functional interpretation were facilitated by an updated version of DEGenes Hunter. Acute exposure of both drugs triggered an intense, specific and transient response at 24 h but with hardly observable differences after 48 h at least in the DEAB treatments. Activation of RA signalling by TTNPB specifically increased the expression of genes in pathways related to RA degradation, retinol storage, carotenoid metabolism, homeostatic response and visual cycle, and also modified the expression of transcripts related to morphogenesis and collagen fibril organisation. In contrast, DEAB mainly decreased genes related to retinal production, impairing phototransduction signalling in the retina. A total of 755 transcripts mainly related to lipid metabolism, lipid transport and lipid homeostasis were altered in response to both treatments, indicating non-specific drug responses associated with intestinal absorption. These results indicate that a new assembling and transcript sieving were both necessary to provide a reliable transcriptome to identify the many aspects of RA action during sole development that are of relevance for sole aquaculture.

Assessment of Growth, Lipid Metabolism and Gene Expression Responses in Senegalese Sole Larvae Fed With Low Dietary Phospholipid Levels

Phospholipids (PL) are essential molecules for larval growth and development. In this study, growth, lipid metabolism and gene expression responses associated with different dietary PL levels in pelagic sole larvae were evaluated. In a first trial, the long-term effects on growth and survival of two experimental microdiets (MD) containing high (High-PL) or low (Low-PL) PL levels were tested and compared to a diet based on live prey (rotifers). The MD were supplied from 3 to 10 days post-hatch (dph) and Artemia from day 8 to 29 dph. High-PL fed larvae had higher dry mass (1.2-fold) than Low-PL fed larvae at 8 dph and both MD were smaller (2.9-fold) than larvae fed live preys. However, a compensatory growth (33% between 8 and 20 dph) occurred when MD were substituted by Artemia and by the end of the trial no significant differences in mass or survival occurred between the dietary treatments. In a second trial, growth, lipid metabolism and gene expression profiles of larvae fed with MD up to 8 dph were analyzed. Growth data confirmed that mass of larvae fed with High-PL was higher (1.3-fold) than the those fed Low-PL and they had lower levels of triacylglycerol (2.8-fold), cholesterol (1.2-fold) and cetoleic acid (1.7-fold). Histological analysis indicated an excess of lipid vacuoles in larvae fed with Low-PL and the expression analysis revealed a coordinated response to enhance lipid mobilization since the expression of genes involved in PL intermediate synthesis, PL remodeling as well as eight apolipoprotein was up-regulated. The down-regulation of apolipoprotein apob2 in larvae fed with Low-PL indicated a specific regulation by PL levels. The present work provides insight into the responses associated with dietary PL in early fish larvae, which will be of use for future studies aimed as designing effective larval sole diets.

Red and White Chinook Salmon (Oncorhynchus tshawytscha): Differences in the Transcriptome Profile of Muscle, Liver, and Pylorus

Astaxanthin (Ax), the main carotenoid responsible for the distinct red flesh color in salmonids (Oncorhynchus, Salvelinus, Salmo, and Parahucho), is added to the diet of farmed fish at a substantial cost. Despite the great economical value for the salmon industry, the key molecular mechanisms involved in the regulation of muscle coloration are poorly understood. Chinook salmon (Oncorhynchus tshawytscha) represent an ideal model to study flesh coloration because they exhibit a distinct color polymorphism responsible for two color morphs, white and red flesh pigmented fish. This study was designed to identify the molecular basis for the development of red and white coloration of fish reared under the same experimental conditions and to better understand the absorption mechanism of Ax in salmonids. Pyloric caeca, liver, and muscle of both groups (n = 6 each) were selected as the most likely critical target organs to be involved respectively in the intestinal uptake, metabolism, and retention of Ax. Difference in the transcriptome profile of each tissue using next-generation sequencing technology was conducted. Ten KEGG pathways were significantly enriched for differentially expressed genes between red and white salmon pylorus tissue, while none for the transcriptome profile in the other two tissues. Differential expressed gene (DE) analyses showed that there were relatively few differences in muscle (31 DE genes, p < 0.05) and liver (43 DE genes, p < 0.05) of white and red Chinook salmon compared approximately 1125 DE genes characterized in the pylorus tissue, with several linked to Ax binding ability, absorption, and metabolism.

Tissue distribution of transcription for 29 lipid metabolism-related genes in Takifugu rubripes, a marine teleost storing lipid predominantly in liver

The tissue distribution pattern of lipid is highly diverse among different fish species. Tiger puffer has a special lipid storage pattern, storing lipid predominantly in liver. In order to better understand the lipid physiology in fish storing lipid in liver, the present study preliminarily investigated the tissue distribution of transcription for 29 lipid metabolism-related genes in tiger puffer, which are involved in lipogenesis, fatty acid oxidation, biosynthesis and hydrolysis of glycerides, lipid transport, and relevant transcription regulation. Samples of eight tissues, brain, eye, heart, spleen, liver, intestine, skin, and muscle, from fifteen juvenile tiger puffer were used in the qRT-PCR analysis. The intestine and brain had high transcription of lipogenic genes, whereas the liver and muscle had low expression levels. The intestine also had the highest transcription level of most apolipoproteins and lipid metabolism-related transcription factors. The transcription of fatty acid β-oxidation-related genes was low in the muscle. The peroxisomal fatty acid oxidation may dominate over mitochondrial β-oxidation in the liver and intestine of tiger puffer, and the MAG pathway probably predominates over the G3P pathway in re-acylation of absorbed lipids in the intestine. The intracellular glyceridases were highly transcribed in the brain, eye, and heart. In conclusion, in tiger puffer, the intestine could be a center of lipid metabolism whereas the liver is more likely a pure storage organ for lipid. The lipid metabolism in the muscle could also be inactive, possibly due to the very low level of intramuscular lipid.

Genomic and phylogenetic analysis of choriolysins, and biological activity of hatching liquid in the flatfish Senegalese sole

The hatching enzymes or choriolysins are key proteases in fish life cycle controlling the release of larvae to surrounding environment that have been suggested as target for novel biotechnological uses. Due to the large amounts of eggs released by the flatfish Solea senegalensis, during the spawning season, the hatching liquid properties and choriolysin-encoding genes were investigated in this species. A genomic analysis identified four putative genes referred to as SseHCEa, SseHCEb, SseLCE and SseHE. The phylogenetic analysis classified these paralogs into two clades, the clade I containing SseHCE paralogs and the clade II containing two well-supported subclades named as HE and LCE. The two SseHCE paralogs were intron-less and both genes were tandemly arrayed very close in the genome. The synteny and gene rearrangement identified in the flatfish lineage indicated that the duplication of these two paralogs occurred recently and they are under divergent evolution. The genes SseHE and SseLCE were structured in 8 exons and 7 introns and the synteny was conserved in teleosts. Expression studies confirmed that the four genes were expressed in the hatching gland cells and they migrate co-ordinately from the head to around the yolk sac close to the hatch with specific temporal and intensity expression profiles. Although the mRNA levels of the four genes peaked in the hours previous to larval hatching, the SseHCE and SseLCE paralogs kept a longer expression than SseHE after hatching. These expression patterns were consistent even when larvae were incubated at different temperatures that modified hatching times. The analysis of hatching-liquid using SDS-PAGE and zymography analyses of hatching liquid identified a major band of expected choriolysin size. The optimal pH for protease activity was 8.5 and inhibition assays using EDTA demonstrated that most of the activity in the hatching liquid was due to metalloproteases with Ca²⁺ ions acting as the most effective metal to restore the activity. All these data provide new clues about the choriolysin evolution and function in flatfish with impact in the aquaculture and the blue cosmetic industry.

Feed and immersion challenges with lymphocystis disease virus (LCDV) reveals specific mechanisms for horizontal transmission and immune response in senegalese sole post-larvae

The horizontal transmission of lymphocystis disease virus (LCDV) through contaminated water and feed (using artemia as vehicle) and the associated immune gene expression profiles in Senegalese sole post-larvae were investigated. All specimens analyzed were positive for LCDV DNA detection at 1-day post-challenge (1 dpc) with the highest viral levels in specimens infected through the immersion route. However, the percentage of LCDV-positive animals and number of viral DNA copies dropped progressively at 2 and 7 dpc. The histological analysis identified structural changes in the skin, muscle and gills of sole post-larvae LCDV-challenged by immersion. In situ hybridization confirmed a wide distribution of LCDV in the skin, gut, surrounding vessels in trunk muscle and head kidney in the immersion route, while the signals were restricted to the liver and lamina propria in the feeding treatment. Expression analysis using a set of 22 genes related to innate immune defense system demonstrated clear differences in the time-course response to LCDV as function of the infection route. Most antiviral defense genes, the proinflammatory cytokines, the complement c3, g-type lysozyme and T-cell markers cd4 and cd8a were rapidly induced in the feeding-infected post-larvae, and they were remained activated at 2 dpc. In contrast, in the immersion-infected post-larvae the induction of most defensive genes was delayed, with a low intensity at 2 dpc. All these data demonstrate that LCDV can horizontally infect Senegalese sole post-larvae through the water or feed although with different patterns of histopathological disorders, virus distribution and route-specific expression profiles.

Phylogeny, expression patterns and regulation of DNA Methyltransferases in early development of the flatfish, Solea senegalensis

Background

The identification of DNA methyltransferases (Dnmt) expression patterns during development and their regulation is important to understand the epigenetic mechanisms that modulate larval plasticity in marine fish. In this study, dnmt1 and dnmt3 paralogs were identified in the flatfish Solea senegalensis and expression patterns in early developmental stages and juveniles were determined. Additionally, the regulation of Dnmt transcription by a specific inhibitor (5-aza-2′-deoxycytidine) and temperature was evaluated.

Results

Five paralog genes of dnmt3, namely dnmt3aa, dnmt3ab, dnmt3ba, dnmt3bb.1 and dnmt3bb.2 and one gene for dnmt1 were identified. Phylogenetic analysis revealed that the dnmt gene family was highly conserved in teleosts and three fish-specific genes, dnmt3aa, dnmt3ba and dnmt3bb.2 have evolved. The spatio-temporal expression patterns of four dnmts (dnmt1, dnmt3aa, dnmt3ab and dnmt3bb.1) were different in early larval stages although all of them reduced expression with the age and were detected in neural organs and dnmt3aa appeared specific to somites. In juveniles, the four dnmt genes were expressed in brain and hematopoietic tissues such as kidney, spleen and gills. Treatment of sole embryos with 5-aza-2′-deoxycytidine down-regulated dntm1 and up-regulated dntm3aa. Moreover, in lecithotrophic larval stages, dnmt3aa and dnmt3ab were temperature sensitive and their expression was higher in larvae incubated at 16 °C relative to 20 °C.

Conclusion

Five dnmt3 and one dnmt1 paralog were identified in sole and their distinct developmental and tissue-specific expression patterns indicate that they may have different roles during development. The inhibitor 5-aza-2′-deoxycytidine modified the transcript abundance of dntm1 and dntm3aa in embryos, which suggests that a regulatory feedback mechanism exists for these genes. The impact of thermal regime on expression levels of dnmt3aa and dnmt3ab in lecithotrophic larval stages suggests that these paralogs might be involved in thermal programing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-017-0154-0) contains supplementary material, which is available to authorized users.

Terminal carbohydrates abundance, immune related enzymes, bactericidal activity and physico-chemical parameters of the Senegalese sole (Solea senegalensis, Kaup) skin mucus

Recently, interest in mucosal surfaces, more specifically fish skin and its secreted mucus, has greatly increased among immunologists. The abundance of terminal carbohydrates, several enzymes (proteases, lysozyme, peroxidase, alkaline phosphatase, esterases and ceruloplasmin), bactericidal activity against fish pathogenic and non-pathogenic bacteria and several physico-chemical parameters (protein concentration, pH, conductivity, redox potential, osmolarity, density and viscosity) in the skin mucus of Senegalese sole (Solea senegalensis, Kaup) have been evaluated. Present results evidence the abundance of N-acetylneuraminic acid, mannose, glucose and N-acetyl-galactosamine in skin mucus. The levels of lysozyme, proteases, esterases and alkaline phosphatase were very similar (from 20 to 30 Units mg^-1 protein). However, 93 Units mg^-1 protein were detected of ceruloplasmin and only 4'88 Units mg^-1 protein of peroxidase. Skin mucus of S. senegalensis showed high bactericidal activity against the tested pathogen bacteria but weak activity against non-pathogenic bacteria. Finally, a clear relationship between mucus density and temperature was detected, while viscosity showed a direct shear- and temperature-dependent behaviour. These results could be useful for better understanding the role of the skin mucus as a key component of the innate immune system, as well as, for elucidating possible relationships between biological and physico-chemical parameters and disease susceptibility.

Molecular characterization and developmental expression patterns of apolipoprotein A-I in Senegalese sole (Solea senegalensis Kaup)

The apolipoprotein A-I (ApoA-I) is an essential component of the high density lipoproteins (HDL). In this study, the cDNA and genomic sequences of this apolipoprotein were characterized for first time in Solea senegalensis. The predicted polypeptide revealed conserved structural features including ten repeats in the lipid-binding domain and some residues involved in cholesterol interaction and binding. The gene structure analysis identified four exons and three introns. Moreover, the synteny analysis revealed that apoA-I did not localize with other apolipoproteins indicating a divergent evolution with respect to the apoA-IV and apoE cluster. The phylogenetic analyses identified two distinct apoA-I paralogs in Ostariophysi (referred to as Ia and Ib) and only one (Ib) in Acanthopterygii. Whole-mount in situ hybridization located the apoA-I signal mainly in the yolk syncytial layer in lecitotrophic larval stages. Later at mouth opening, the mRNA signals were detected mainly in liver and intestine compatible with its role in the HDL formation. Moreover, a clear signal was detected in some regions of the brain, retina and neural cord suggesting a role in local regulation of cholesterol homeostasis. After metamorphosis, apoA-I was also detected in other tissues such as gills, head kidney and spleen suggesting a putative role in immunity. Expression analyses in larvae fed two diets with different triacylglycerol levels indicated that apoA-I mRNA levels were more associated to larval size and development than dietary lipid levels. Finally, qPCR analyses of immature and mature transcripts revealed distinct expression profiles suggesting a posttranscriptional regulatory mechanism.