Molecular characterization and gene expression of six trypsinogens in the flatfish Senegalese sole (Solea senegalensis Kaup) during larval development and in tissues

Linking Pedigree Information to the Gene Expression Phenotype to Understand Differential Family Survival Mechanisms in Highly Fecund Fish: A Case Study in the Larviculture of Pacific Bluefin Tuna

Mass spawning in fish culture often brings about a marked variance in family size, which can cause a reduction in effective population sizes in seed production for stock enhancement. This study reports an example of combined pedigree information and gene expression phenotypes to understand differential family survival mechanisms in early stages of Pacific bluefin tuna, Thunnus orientalis, in a mass culture tank. Initially, parentage was determined using the partial mitochondrial DNA control region sequence and 11 microsatellite loci at 1, 10, 15, and 40 days post-hatch (DPH). A dramatic proportional change in the families was observed at around 15 DPH; therefore, transcriptome analysis was conducted for the 15 DPH larvae using a previously developed oligonucleotide microarray. This analysis successfully addressed the family-specific gene expression phenotypes with 5739 differentially expressed genes and highlighted the importance of expression levels of gastric-function-related genes at the developmental stage for subsequent survival. This strategy demonstrated herein can be broadly applicable to species of interest in aquaculture to comprehend the molecular mechanism of parental effects on offspring survival, which will contribute to the optimization of breeding technologies.

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.

Trypsin gene expression in adults and larvae of tropical gar Atractosteus tropicus

Trypsin gene (try) expression levels were quantified in different organs of wild and captive tropical gar (Atractosteus tropicus) adults, and changes in expression during initial ontogeny of the species were determined. RNA was extracted from the pancreas, and cDNA was synthesized and later amplified by endpoint PCR using oligonucleotides designed from different try sequences of fish registered in GenBank. Subsequently, specific oligonucleotides were designed from the partial sequences. Gene expression was measured after RNA extraction and synthesis of the cDNA of 11 organs (liver, pancreas, stomach, esophagus, intestine, pyloric caeca, brain, muscle, gills, gonad, and kidney) of captive and wild adults. Likewise, samples of A. tropicus larvae were taken on days 0 (embryo), 5, 10, 15, 20, 25, and 30 days after hatching (DAH), the RNA was extracted, and the synthesis of cDNA was carried out to measure real-time gene expression (qPCR). The results showed that the highest relative try expression occurred mainly in the esophagus, liver, stomach, and pancreas of both wild and captive adult fish; however, captive organisms had a higher try expression level than wild fish. Although try expression during initial ontogeny was high in embryos (0 DAH), it did not reach the maximum value until 15 DAH. It was concluded that try expression levels in captive adults are due to the high protein content in the balanced feed (trout diet). The highest try expression level during larviculture was detected at 15 DAH, which indicates that A. tropicus larvae have a mature digestive system and can efficiently hydrolyze proteins from feed at this developmental stage.

Peptidylarginine deiminase and deiminated proteins are detected throughout early halibut ontogeny - Complement components C3 and C4 are post-translationally deiminated in halibut (Hippoglossus hippoglossus L.)

Post-translational protein deimination is mediated by peptidylarginine deiminases (PADs), which are calcium dependent enzymes conserved throughout phylogeny with physiological and pathophysiological roles. Protein deimination occurs via the conversion of protein arginine into citrulline, leading to structural and functional changes in target proteins. In a continuous series of early halibut development from 37 to 1050° d, PAD, total deiminated proteins and deiminated histone H3 showed variation in temporal and spatial detection in various organs including yolksac, muscle, skin, liver, brain, eye, spinal cord, chondrocytes, heart, intestines, kidney and pancreas throughout early ontogeny. For the first time in any species, deimination of complement components C3 and C4 is shown in halibut serum, indicating a novel mechanism of complement regulation in immune responses and homeostasis. Proteomic analysis of deiminated target proteins in halibut serum further identified complement components C5, C7, C8 C9 and C1 inhibitor, as well as various other immunogenic, metabolic, cytoskeletal and nuclear proteins. Post-translational deimination may facilitate protein moonlighting, an evolutionary conserved phenomenon, allowing one polypeptide chain to carry out various functions to meet functional requirements for diverse roles in immune defences and tissue remodelling.

Purification and characterization of a sarcoplasmic serine proteinase from threadfin bream Nemipterus virgatus muscle

A sarcoplasmic serine proteinase (SSP) was purified from threadfin bream (Nemipterus virgatus) belly muscle by ammonium sulfate precipitation and a series of chromatographies including Q-Sepharose, Phenyl Sepharose and Superdex 200. The SSP was purified 1967 folds with a yield of 4.8%. The molecular weight of the SSP was estimated to be 43.5 kDa and 22.5 kDa on SDS-PAGE under non-reducing and reducing conditions, respectively. The N-terminal amino acid sequence of the two protein bands were determined as IVGGYEXQPYSQAHQVSLNSGY and corresponded. It is suggested that the SSP exists as a homodimer. Optimum pH and temperature were 9.5 and 50 °C, using Boc-Val-Pro-Arg-MCA as a substrate. Substrate specificity and effects of inhibitors indicated that the SSP was a trypsin-like serine proteinase. The SSP was responsible for hydrolyzing myosin heavy chain (MHC) and inducing modori phenomenon in the threadfin bream surimi gel. Thus, the SSP was considered as a modori-inducing proteinase.

The yellowtail (Seriola quinqueradiata) genome and transcriptome atlas of the digestive tract

Seriola quinqueradiata (yellowtail) is the most widely farmed and economically important fish in aquaculture in Japan. In this study, we used the genome of haploid yellowtail fish larvae for de novo assembly of whole-genome sequences, and built a high-quality draft genome for the yellowtail. The total length of the assembled sequences was 627.3 Mb, consisting of 1,394 scaffold sequences (>2 kb) with an N50 length of 1.43 Mb. A total of 27,693 protein-coding genes were predicted for the draft genome, and among these, 25,832 predicted genes (93.3%) were functionally annotated. Given our lack of knowledge of the yellowtail digestive system, and using the annotated draft genome as a reference, we conducted an RNA-Seq analysis of its three digestive organs (stomach, intestine and rectum). The RNA-Seq results highlighted the importance of certain genes in encoding proteolytic enzymes necessary for digestion and absorption in the yellowtail gastrointestinal tract, and this finding will accelerate development of formulated feeds for this species. Since this study offers comprehensive annotation of predicted protein-coding genes, it has potential broad application to our understanding of yellowtail biology and aquaculture.

Dietary administration of the probiotic Shewanella putrefaciens Pdp11 promotes transcriptional changes of genes involved in growth and immunity in Solea senegalensis larvae

Senegalese sole (Solea senegalensis) has been proposed as a high-potential species for aquaculture diversification in Southern Europe. It has been demonstrated that a proper feeding regimen during the first life stages influences larval growth and survival, as well as fry and juvenile quality. The bacterial strain Shewanella putrefaciens Pdp11 (SpPdp11) has shown very good probiotic properties in Senegalese sole, but information is scarce about its effect in the earliest stages of sole development. Thus, the aim of this study was to investigate the effect of SpPdp11, bioencapsulated in live diet, administered during metamorphosis (10-21 dph) or from the first exogenous feeding of Senegalese sole (2-21 dph). To evaluate the persistence of the probiotic effect, we sampled sole specimens from metamorphosis until the end of weaning (from 23 to 73 dph). This study demonstrated that probiotic administration from the first exogenous feeding produced beneficial effects on Senegalese sole larval development, given that specimens fed this diet exhibited higher and less dispersed weight, as well as increases in both total protein concentration and alkaline phosphatase activity, and in non-specific immune response. Moreover, real-time PCR documented changes in the expression of a set of genes involved in central metabolic functions including genes related to growth, genes coding for proteases (including several digestive enzymes), and genes implicated in the response to stress and in immunity. Overall, these results support the application of SpPdp11 in the first life stages of S. senegalensis as an effective tool with the clear potential to benefit sole aquaculture.

Small ocean temperature increases elicit stage-dependent changes in DNA methylation and gene expression in a fish, the European sea bass

In natural fish populations, temperature increases can result in shifts in important phenotypic traits. DNA methylation is an epigenetic mechanism mediating phenotypic changes. However, whether temperature increases of the magnitude predicted by the latest global warming models can affect DNA methylation is unknown. Here, we exposed European sea bass to moderate temperature increases in different periods within the first two months of age. We show that increases of even 2 °C in larvae significantly changed global DNA methylation and the expression of ecologically-relevant genes related to DNA methylation, stress response, muscle and organ formation, while 4 °C had no effect on juveniles. Furthermore, DNA methylation changes were more marked in larvae previously acclimated to a different temperature. The expression of most genes was also affected by temperature in the larvae but not in juveniles. In conclusion, this work constitutes the first study of DNA methylation in fish showing that temperature increases of the magnitude predicted by the latest global warming models result in stage-dependent alterations in global DNA methylation and gene expression levels. This study, therefore, provides insights on the possible consequences of climate change in fish mediated by genome-wide epigenetic modifications.

Express Sequence Tag Analysis - Identification of Anseriformes Trypsin Genes from Full-Length cDNA Library of the Duck (Anas platyrhynchos) and Characterization of Their Structure and Function

Trypsins are key proteins important in animal protein digestion by breaking down the peptide bonds on the carboxyl side of lysine and arginine residues, hence it has been used widely in various biotechnological processes. In the current study, a full-length cDNA library with capacity of 5·10(5) CFU/ml from the duck (Anas platyrhynchos) was constructed. Using express sequence tag (EST) sequencing, genes coding two trypsins were identified and two full-length trypsin cDNAs were then obtained by rapid-amplification of cDNA end (RACE)-PCR. Using Blast, they were classified into the trypsin I and II subfamilies, but both encoded a signal peptide, an activation peptide, and a 223-a.a. mature protein located in the C-terminus. The two deduced mature proteins were designated as trypsin-IAP and trypsin-IIAP, and their theoretical isoelectric points (pI) and molecular weights (MW) were 7.99/23466.4 Da and 4.65/24066.0 Da, respectively. Molecular characterizations of genes were further performed by detailed bioinformatics analysis. Phylogenetic analysis revealed that trypsin-IIAP has an evolution pattern distinct from trypsin-IAP, suggesting its evolutionary advantage. Then the duck trypsin-IIAP was expressed in an Escherichia coli system, and its kinetic parameters were measured. The three dimensional structures of trypsin-IAP and trypsin-IIAP were predicted by homology modeling, and the conserved residues required for functionality were identified. Two loops controlling the specificity of the trypsin and the substrate-binding pocket represented in the model are almost identical in primary sequences and backbone tertiary structures of the trypsin families.

Elucidation of different cold-adapted Atlantic cod (Gadus morhua) trypsin X isoenzymes

Trypsins from Atlantic cod (Gadus morhua), consisting of several isoenzymes, are highly active cold-adapted serine proteases. These trypsins are isolated for biomedical use in an eco-friendly manner from underutilized seafood by-products. Our group has explored the biochemical properties of trypsins and their high potential in biomedicine. For broader utilization of cod trypsins, further characterization of biochemical properties of the individual cod trypsin isoenzymes is of importance. For that purpose, a benzamidine purified trypsin isolate from Atlantic cod was analyzed. Anion exchange chromatography revealed eight peaks containing proteins around 24kDa with tryptic activity. Based on mass spectrometric analysis, one isoenzyme gave the best match to cod trypsin I and six isoenzymes gave the best match to cod trypsin X. Amino terminal sequencing of two of these six trypsin isoenzymes showed identity to cod trypsin X. Three sequence variants of trypsin X were identified by cDNA analysis demonstrating that various forms of this enzyme exist. One trypsin X isoenzyme was selected for further characterization based on abundance and stability. Stepwise increase in catalytic efficiency (k_cat/K_m) of this trypsin X isoenzyme was obtained with substrates containing one to three amino acid residues. The study demonstrates that the catalytic efficiency of this trypsin X isoenzyme is comparable to that of cod trypsin I, the most abundant and highly active isoenzyme in the benzamidine cod trypsin isolate. Differences in pH stability and sensitivity to inhibitors of the trypsin X isoenzyme compared to cod trypsin I were detected that may be important for practical use.

Maternally derived trypsin may have multiple functions in the early development of turbot (Scopthalmus maximus)

Trypsin is an important serine protease that is considered to be involved in digestion of protein in teleost fish. Nevertheless, studies on trypsin/trypsinogen in fish embryos are very limited. In this study, the trypsinogen of turbot (Scophthalmus maximus) (tTG) was identified and the expression patterns and activity of trypsinogen/trypsin were investigated. The results showed that the tTG mRNA was evenly distributed in the oocytes and was also expressed along the yolk periphery in early embryos. At later embryo stages and 1 days after hatching (dph), the tTG mRNA concentrated at the alimentary tract and head. Quantitative expression analysis showed that the tTG transcripts decreased after fertilization until the gastrula stage, then increased with the embryo and larvae development. This result was also confirmed by the specific activity analysis of trypsin and in-situ-hybridization (ISH). All of the results indicated that tTG in early embryo stages was maternally derived and expressed by itself after gastrula stages. Additionally, location of tTG mRNA in embryos and larvae was investigated; we considered that trypsin may have multiple functions during the embryo development process. Based on our results regarding trypsinogen in embryos and early development, we concluded that the trypsin/trypsinogen in turbot embryos was inherited from a maternal source and we suggested that trypsin in early development has multiple functions in the process of development.

De novo assembly, characterization and functional annotation of Senegalese sole (Solea senegalensis) and common sole (Solea solea) transcriptomes: integration in a database and design of a microarray

Background

Senegalese sole (Solea senegalensis) and common sole (S. solea) are two economically and evolutionary important flatfish species both in fisheries and aquaculture. Although some genomic resources and tools were recently described in these species, further sequencing efforts are required to establish a complete transcriptome, and to identify new molecular markers. Moreover, the comparative analysis of transcriptomes will be useful to understand flatfish evolution.

Results

A comprehensive characterization of the transcriptome for each species was carried out using a large set of Illumina data (more than 1,800 millions reads for each sole species) and 454 reads (more than 5 millions reads only in S. senegalensis), providing coverages ranging from 1,384x to 2,543x. After a de novo assembly, 45,063 and 38,402 different transcripts were obtained, comprising 18,738 and 22,683 full-length cDNAs in S. senegalensis and S. solea, respectively. A reference transcriptome with the longest unique transcripts and putative non-redundant new transcripts was established for each species. A subset of 11,953 reference transcripts was qualified as highly reliable orthologs (>97% identity) between both species. A small subset of putative species-specific, lineage-specific and flatfish-specific transcripts were also identified. Furthermore, transcriptome data permitted the identification of single nucleotide polymorphisms and simple-sequence repeats confirmed by FISH to be used in further genetic and expression studies. Moreover, evidences on the retention of crystallins crybb1, crybb1-like and crybb3 in the two species of soles are also presented. Transcriptome information was applied to the design of a microarray tool in S. senegalensis that was successfully tested and validated by qPCR. Finally, transcriptomic data were hosted and structured at SoleaDB.

Conclusions

Transcriptomes and molecular markers identified in this study represent a valuable source for future genomic studies in these economically important species. Orthology analysis provided new clues regarding sole genome evolution indicating a divergent evolution of crystallins in flatfish. The design of a microarray and establishment of a reference transcriptome will be useful for large-scale gene expression studies. Moreover, the integration of transcriptomic data in the SoleaDB will facilitate the management of genomic information in these important species.

Electronic supplementary material

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

Molecular characterization of two trypsinogens in the orange-spotted grouper, Epinephelus coioides, and their expression in tissues during early development

In this study, we cloned two trypsinogens of the orange-spotted grouper, Epinephelus coioides, and analyzed their structure, expression, and activity. Full-length trypsinogen complementary (c)DNAs, named T1 and T2, were 900 and 875 nucleotides, and translated 242 and 244 deduced amino acid peptides, respectively. Both trypsinogens contained highly conserved residues essential for serine protease catalytic and conformational maintenance. Results from isoelectric and phylogenetic analyses suggested that both trypsinogens were grouped into trypsinogen group I. Both trypsinogens had similar expression patterns of negative relationship with body weight; expression was first detected at 1 day post-hatching (DPH) and exhibited steady-state expression during early development at 1-25 DPH. Both expression and activity levels significantly increased after 30 DPH due to metamorphosis. Grouper larval development is very slow with insignificant changes in total length and body weight before 8 DPH. The contribution of live food to an increase in the trypsin activity profile may explain their importance in food digestion and survival of larvae during early larval development.

Metamorphosis in teleosts

Teleosts are the largest and most diverse group of vertebrates, and many species undergo morphological, physiological, and behavioral transitions, "metamorphoses," as they progress between morphologically divergent life stages. The larval metamorphosis that generally occurs as teleosts mature from larva to juvenile involves the loss of embryo-specific features, the development of new adult features, major remodeling of different organ systems, and changes in physical proportions and overall phenotype. Yet, in contrast to anuran amphibians, for example, teleost metamorphosis can entail morphological change that is either sudden and profound, or relatively gradual and subtle. Here, we review the definition of metamorphosis in teleosts, the diversity of teleost metamorphic strategies and the transitions they involve, and what is known of their underlying endocrine and genetic bases. We suggest that teleost metamorphosis offers an outstanding opportunity for integrating our understanding of endocrine mechanisms, cellular processes of morphogenesis and differentiation, and the evolution of diverse morphologies and life histories.

Purification and characterization of trypsin from the pyloric ceca of orange-spotted grouper, Epinephelus coioides

Trypsin from the pyloric ceca of orange-spotted grouper, Epinephelus coioides, was purified by fractionation with ammonium sulfate, ionic exchange, and affinity chromatography. The protein was purified 161.85-fold with a yield of 4%. Purified trypsin had an apparent molecular weight of 24 kDa according to an SDS-PAGE analysis. Optimal profiles of temperature and pH of the enzyme were 50°C and 8-10, respectively, using Nα-benzoyl-L: -arginine ethyl ester as the substrate. The results of thermal and pH stability assays showed that the enzyme was stable at temperatures of up to 50°C and in the pH range of 6-8. Trypsin activity decreased with an increasing NaCl concentration (0-0.6 M). The activity of purified trypsin was effectively inhibited by a soybean trypsin inhibitor and N-p-tosyl-L: -lysine chloromethyl ketone, and was slightly inhibited by iodoacetic acid, ethylenediaminetetraacetic acid, 1-(L: -trans-epoxysuccinyl-leucylamino)-4-guanidinobutane, and pepstatin A. Protein identification of the purified protease showed that the sequences of two peptides, LGEHNI and NLDNDIML, were highly homologous to other fish trypsins. The measurement of trypsin activity in different tissues showed that the highest activity was detected in pyloric ceca, followed by anterior intestine, middle intestine, hind intestine and spleen, but very low activities were found in other tissues. An inverse relationship between the trypsin activity in four tissues of pyloric ceca, anterior intestine, middle intestine and hind intestine and fish body weight as a result of increased pepsin in stomach indicated grouper growth status was increased.

Effects of feeding with corn gluten meal on trypsin activity and mRNA expression in Fugu obscurus

The molecular regulation of the trypsin activity and the effects of corn gluten meal (CGM) on the trypsin activity and gene expression in Fugu obscurus were examined. The fish were fed with diets containing 0, 50, 100, 150, and 200 g kg⁻¹ CGM for 60 days. The trypsin activity and mRNA levels in the intestine, stomach, and hepatopancreas were measured before and after the experimental feeding. Enzyme activity was highest in the intestine, and mRNA expression was highest in the hepatopancreas prior to feeding the experimental diets. The effect of dietary CGM on the trypsin activity was dependent on the concentration. Enzyme activity was significantly lower in the hepatopancreas and stomach of the 200 g kg⁻¹ treatment group than in control group. In contrast, activity was enhanced in the stomach of the 100 g kg⁻¹ treatment group. Different amounts of CGM in the diets also resulted in significant changes in the trypsin mRNA levels in all three tissues. The trypsin mRNA expression in stomach, hepatopancreas, and intestine increased with increasing CGM from 0 to 10% and decreased with increasing CGM from 10 to 20%. The trypsin mRNA expression was highest in the 100 g kg⁻¹ treatment group, which was significantly higher than that of in the control. In addition, the changes in the trypsin mRNA expression levels were not reflected by changes in the trypsin activity. Our results suggest that the trypsin activity is regulated not only at the transcription level but also during translation and posttranslational modification.

Effect of different diets on proteolytic enzyme activity, trypsinogen gene expression and dietary carbon assimilation in Senegalese sole (Solea senegalensis) larvae

The effect of diet on larval growth, anionic trypsinogen gene expression (ssetryp1), and trypsin (EC 3.4.21.4) and chymotrypsin (EC 3.4.21.1) activities was assessed in Solea senegalensis. Changes in larval carbon stable isotope (δ(13)C) composition were used to estimate carbon assimilation. Diets were supplied for 20days to fish held in larval rearing tanks and consisted of live rotifers, Artemia sp. nauplii, rotifers followed by Artemia sp., rotifers co-fed with inert diet and inert diet alone. Growth was significantly faster in larvae fed only Artemia and those fed rotifers and Artemia (k=0.381-0.387day(-1)). Trypsin and chymotrypsin activities increased from 3 to 4days after hatching (DAH) in all dietary treatments, while ssetryp1 transcripts increased at 4-5 DAH only in larvae fed live prey. ssetryp1 gene expression was activated later in larvae fed only Artemia and this corresponded with Artemia δ(13)C values being reflected in larval tissue. Larval δ(13)C values also indicated greater selection and/or assimilation of rotifers in relation to the inert diet. Results demonstrate that during early larval development of sole, diet modulates ssetryp1 gene expression. The rapid and intense response to diets that promoted different growth and survival suggests the suitability of this biomarker as a nutritional status indicator in early sole larvae.

Molecular characterization of trypsinogens and development of trypsinogen gene expression and tryptic activities in grass carp (Ctenopharyngodon idellus) and topmouth culter (Culter alburnus)

This study examined the gene structures and expression of trypsinogens, as well as the trypsin activities of the grass carp Ctenopharyngodon idellus (herbivorous) and the topmouth culter Culter alburnus (carnivorous), which are commercially important freshwater species of the family Cyprinidae in China. Isolated full-length trypsinogen cDNA clones were 869 bp and 857 bp. The deduced amino acid sequences were 242aa and 247aa long, both containing the highly conserved residues essential for serine protease catalytic and conformational maintenance. The results from isoelectric and phylogenetic analyses suggest that grass carp trypsinogen is grouped with teleost trypsinogen group I, while topmouth culter trypsinogen is grouped with group II. The expression pattern of trypsinogen mRNA was similar between these two species, appearing 2 days post-hatching (dph) and reaching peaks at 11 and 23 dph. The trypsin-specific activities in both species were detected 2 dph and reached the major peaks at 8 dph, however the minor peaks were observed at 20 dph in the grass carp and 17 dph in the topmouth culter. The trypsin-specific activity was significantly higher in the grass carp than in the topmouth culter, which may be attributed to the nature of their different nutritional habits.