Gilthead sea bream (Sparus auratus) and European sea bass (Dicentrarchus labrax) expressed sequence tags: Characterization, tissue-specific expression and gene markers

Skeletal Morphogenesis and Anomalies in Gilthead Seabream: A Comprehensive Review

The gilthead seabream, one of the most important species in Mediterranean aquaculture, with an increasing status of exploitation in terms of production volume and aquafarming technologies, has become an important research topic over the years. The accumulation of knowledge from several studies conducted during recent decades on their functional and biological characteristics has significantly improved their aquacultural aspects, namely their reproductive success, survival, and growth. Despite the remarkable progress in the aquaculture industry, hatchery conditions are still far from ideal, resulting in frequent abnormalities at the beginning of intensive culture, entailing significant economic losses. Those deformities are induced during the embryonic and post-embryonic periods of life, and their development is still poorly understood. In the present review, we created a comprehensive synthesis that covers the various aspects of skeletal morphogenesis and anomalies in the gilthead seabream, highlighting the genetic, environmental, and nutritional factors contributing to bone deformities and emphasized the potential of the gilthead seabream as a model organism for understanding bone morphogenesis in both aquaculture and translational biological research. This review article addresses the existing lack in the literature regarding gilthead seabream bone deformities, as there are currently no comprehensive reviews on this subject.

Single Nucleotide Polymorphism Markers with Applications in Conservation and Exploitation of Aquatic Natural Populations

Simple Summary

In recent decades, societies, states and local authorities have become increasingly aware that for effective long-term management and protection of aquatic ecosystems and populations, it is necessary to take into account the genetic changes occurring in these populations. One type of high-resolution molecular marker suitable for studying the neutral and adaptive genetic diversity of populations is single nucleotide polymorphism (SNP). This review is an attempt to show the benefits of using SNPs to recognize natural populations of aquatic animals and detect the threats to them from accidentally or intentionally released farm animals, fishery and global climate changes. It is postulated that conservation actions should protect not only pristine natural populations that are endangered or overfished, but also protect populations of non-threatened species from unnecessarily released semi-domesticated animals. The enhancement of natural populations with farmed material usually reduces their genetic diversity. Experimental size-selective catches of artificially created populations have caused evolutionary changes in the life cycles of fishes. However, fishery-induced evolution in natural populations is difficult to observe. The negative measurable effects on populations can be expected when the number of breeding individuals is reduced below 100, which occurs very rarely in the sea and more often in fragmented freshwater streams, ponds and seasonal rivers.

Abstract

An increasing number of aquatic species have been studied for genetic polymorphism, which extends the knowledge on their natural populations. One type of high-resolution molecular marker suitable for studying the genetic diversity of large numbers of individuals is single nucleotide polymorphism (SNP). This review is an attempt to show the range of applications of SNPs in studies of natural populations of aquatic animals. In recent years, SNPs have been used in the genetic analysis of wild and enhanced fish and invertebrate populations in natural habitats, exploited migratory species in the oceans, migratory anadromous and freshwater fish and demersal species. SNPs have been used for the identification of species and their hybrids in natural environments, to study the genetic consequences of restocking for conservation purposes and the negative effects on natural populations of fish accidentally escaping from culture. SNPs are very useful for identifying genomic regions correlated with phenotypic variants relevant for wildlife protection, management and aquaculture. Experimental size-selective catches of populations created in tanks have caused evolutionary changes in life cycles of fishes. The research results have been discussed to clarify whether the fish populations in natural conditions can undergo changes due to selective harvesting targeting the fastest-growing fishes.

The European sea bass: a key marine fish model in the wild and in aquaculture

The European sea bass (Dicentrarchus labrax L.) is a marine fish of key economic and cultural importance in Europe. It is now more an aquaculture than a fisheries species (>96% of the production in 2016), although modern rearing techniques date back only from the late 1980s. It also has high interest for evolutionary studies, as it is composed of two semispecies (Atlantic and Mediterranean lineages) that have come into secondary contact following the last glaciation. Based on quantitative genetics studies of most traits of interest over the past 10–15 years, selective breeding programs are now applied to this species, which is at the beginning of its domestication process. The availability of a good quality reference genome has accelerated the development of new genomic resources, including SNP arrays that will enable genomic selection to improve genetic gain. There is a need to improve feed efficiency, both for economic and environmental reasons, but this will require novel phenotyping approaches. Further developments will likely focus on the understanding of genotype‐by‐environment interactions, which will be important both for efficient breeding of farmed stocks and for improving knowledge of the evolution of natural populations. At the interface between both, the domestication process must be better understood to improve production and also to fully evaluate the possible impact of aquaculture escapees on wild populations. The latter is an important question for all large‐scale aquaculture productions.

Transcription of histones H1 and H2B is regulated by several immune stimuli in gilthead seabream and European sea bass

Histones (H1 to H4) are the primary proteins which mediate the folding of DNA into chromatin; however, and in addition to this function, histones have been also related to antimicrobial peptides (AMPs) activity in vertebrates, in fact, mammalian H1 is mobilized as part as the anti-viral immune response. In fish, histones with AMP activity have been isolated and characterized mainly from skin and gonads. One of most threatening pathogens for wild and cultured fish species nowadays is nodavirus (NNV), which target tissues are the brain and retina, but it is also able to colonize the gonad and display vertical transmission. Taking all this into account we have identified the h1 and h2b coding sequences in European sea bass (Dicentrarchus labrax) and gilthead seabream (Sparus aurata) fish species and studied their pattern of expression under naïve conditions and NNV in vivo infection. The data obtained prompted us to study their role on the immune response of gonad and head-kidney leucocytes upon viral (NNV), bacteria (Vibrio anguillarum or Photobacterium damselae), pathogen-associated molecular patterns (PAMPs) or mitogens stimulation. The h1 and h2b genes are expressed in a wide range of tissues and their expression is modify by infection or other immune stimuli, but further studies will be needed to determine the significance of these changes. These results suggest that h1 expression is related to the immune response against NNV in the brain, while h2b transcription seems to be more important in the head-kidney. Moreover, the potential role of histones as anti-viral agents is suggested and further characterization is in progress.

Advances in European sea bass genomics and future perspectives

Only recently available sequenced and annotated teleost fish genomes were restricted to a few model species, none of which were for aquaculture. The application of marker assisted selection for improved production traits had been largely restricted to the salmon industry and genetic and Quantitative Trait Loci (QTL) maps were available for only a few species. With the advent of next generation sequencing the landscape is rapidly changing and today the genomes of several aquaculture species have been sequenced. The European sea bass, Dicentrarchus labrax, is a good example of a commercially important aquaculture species in Europe for which in the last decade a wealth of genomic resources, including a chromosomal scale genome assembly, physical and linkage maps as well as relevant QTL have been generated. The current challenge is to stimulate the uptake of the resources by the industry so that the full potential of this scientific endeavor can be exploited and produce benefits for producers and the public alike.

Efficient isolation of polymorphic microsatellites from high-throughput sequence data based on number of repeats

Transcriptome data are a good resource to develop microsatellites due to their potential in targeting candidate genes. However, developing microsatellites can be a time-consuming enterprise due to the numerous primer pairs to be tested. Therefore, the use of methodologies that make it efficient to identify polymorphic microsatellites is desirable. Here we used a 62,038 contigs transcriptome assembly, obtained from pyrosequencing a peacock blenny (Salaria pavo) multi-tissue cDNA library, to mine for microsatellites and in silico evaluation of their polymorphism. A total of 4190 microsatellites were identified in 3670 unique unigenes, and from these microsatellites, in silico polymorphism was detected in 733. We selected microsatellites based either on their in silico polymorphism and annotation results or based only on their number of repeats. Using these two approaches, 28 microsatellites were successfully amplified in twenty-six individuals, and all but 2 were found to be polymorphic, being the first genetic markers for this species. Our results showed that the strategy of selection based on number of repeats is more efficient in obtaining polymorphic microsatellites than the strategy of in silico polymorphism (allelic richness was 8.2±3.85 and 4.56±2.45 respectively). This study demonstrates that combining the knowledge of number of repeats with other predictors of variability, for example in silico microsatellite polymorphism, improves the rates of polymorphism, yielding microsatellites with higher allelic richness, and decreases the number of monomorphic microsatellites obtained.

Deep sequencing for de novo construction of a marine fish (Sparus aurata) transcriptome database with a large coverage of protein-coding transcripts

Background

The gilthead sea bream (Sparus aurata) is the main fish species cultured in the Mediterranean area and constitutes an interesting model of research. Nevertheless, transcriptomic and genomic data are still scarce for this highly valuable species. A transcriptome database was constructed by de novo assembly of gilthead sea bream sequences derived from public repositories of mRNA and collections of expressed sequence tags together with new high-quality reads from five cDNA 454 normalized libraries of skeletal muscle (1), intestine (1), head kidney (2) and blood (1).

Results

Sequencing of the new 454 normalized libraries produced 2,945,914 high-quality reads and the de novo global assembly yielded 125,263 unique sequences with an average length of 727 nt. Blast analysis directed to protein and nucleotide databases annotated 63,880 sequences encoding for 21,384 gene descriptions, that were curated for redundancies and frameshifting at the homopolymer regions of open reading frames, and hosted at http://www.nutrigroup-iats.org/seabreamdb. Among the annotated gene descriptions, 16,177 were mapped in the Ingenuity Pathway Analysis (IPA) database, and 10,899 were eligible for functional analysis with a representation in 341 out of 372 IPA canonical pathways. The high representation of randomly selected stickleback transcripts by Blast search in the nucleotide gilthead sea bream database evidenced its high coverage of protein-coding transcripts.

Conclusions

The newly assembled gilthead sea bream transcriptome represents a progress in genomic resources for this species, as it probably contains more than 75% of actively transcribed genes, constituting a valuable tool to assist studies on functional genomics and future genome projects.

Genomics in eels--towards aquaculture and biology

Freshwater eels (genus Anguilla), especially the species inhabiting the temperate areas such as the European, American and Japanese eels, are important aquaculture species. Although artificial reproduction has been attempted since the 1930s and large numbers of studies have been conducted, it has not yet fully succeeded. Problems in eel artificial breeding are highly diverse, for instance, lack of basic information about reproduction in nature, no appropriate food for larvae, high mortality, and high individual variation in adults in response to maturation induction. Over the last decade, genomic data have been obtained for a variety of aquatic organisms. Recent technological advances in sequencing and computation now enable the accumulation of genomic information even for non-model species. The draft genome of the European eel Anguilla anguilla has been recently determined using Illumina technology and transcriptomic data based on next generation sequencing have been emerging. Extensive genomic information will facilitate many aspects of the artificial reproduction of eels. Here, we review the progress in genome-wide studies of eels, including additional analysis of the European eel genome data, and discuss future directions and implications of genomic data for aquaculture.

Transcriptomic characterization of the larval stage in gilthead seabream (Sparus aurata) by 454 pyrosequencing

Gilthead seabream (Sparus aurata) is a teleost belonging to the family Sparidae with a high economical relevance in the Mediterranean countries. Although genomic tools have been developed in this species in order to investigate its physiology at the molecular level and consequently its culture, genomic information on post-embryonic development is still scarce. In this study, we have investigated the transcriptome of a marine teleost during the larval stage (from hatching to 60 days after hatching) by the use of 454 pyrosequencing technology. We obtained a total of 68,289 assembled contigs, representing putative transcripts, belonging to 54,606 different clusters. Comparison against all S. aurata expressed sequenced tags (ESTs) from the NCBI database revealed that up to 34,722 contigs, belonging to about 61% of gene clusters, are sequences previously not described. Contigs were annotated through an iterative Blast pipeline by comparison against databases such as NCBI RefSeq from Danio rerio, SwissProt or NCBI teleost ESTs. Our results indicate that we have enriched the number of annotated sequences for this species by more than 50% compared with previously existing databases for the gilthead seabream. Gene Ontology analysis of these novel sequences revealed that there is a statistically significant number of transcripts with key roles in larval development, differentiation, morphology, and growth. Finally, all information has been made available online through user-friendly interfaces such as GBrowse and a Blast server with a graphical frontend.

Fast skeletal muscle transcriptome of the gilthead sea bream (Sparus aurata) determined by next generation sequencing

Background

The gilthead sea bream (Sparus aurata L.) occurs around the Mediterranean and along Eastern Atlantic coasts from Great Britain to Senegal. It is tolerant of a wide range of temperatures and salinities and is often found in brackish coastal lagoons and estuarine areas, particularly early in its life cycle. Gilthead sea bream are extensively cultivated in the Mediterranean with an annual production of 125,000 metric tonnes. Here we present a de novo assembly of the fast skeletal muscle transcriptome of gilthead sea bream using 454 reads and identify gene paralogues, splice variants and microsatellite repeats. An annotated transcriptome of the skeletal muscle will facilitate understanding of the genetic and molecular basis of traits linked to production in this economically important species.

Results

Around 2.7 million reads of mRNA sequence data were generated from the fast myotomal of adult fish (~2 kg) and juvenile fish (~0.09 kg) that had been either fed to satiation, fasted for 3-5d or transferred to low (11°C) or high (33°C) temperatures for 3-5d. Newbler v2.5 assembly resulted in 43,461 isotigs >100 bp. The number of sequences annotated by searching protein and gene ontology databases was 10,465. The average coverage of the annotated isotigs was x40 containing 5655 unique gene IDs and 785 full-length cDNAs coding for proteins containing 58–1536 amino acids. The v2.5 assembly was found to be of good quality based on validation using 200 full-length cDNAs from GenBank. Annotated isotigs from the reference transcriptome were attributable to 344 KEGG pathway maps. We identified 26 gene paralogues (20 of them teleost-specific) and 43 splice variants, of which 12 had functional domains missing that were likely to affect their biological function. Many key transcription factors, signaling molecules and structural proteins necessary for myogenesis and muscle growth have been identified. Physiological status affected the number of reads that mapped to isotigs, reflecting changes in gene expression between treatments.

Conclusions

We have produced a comprehensive fast skeletal muscle transcriptome for the gilthead sea bream, which will provide a resource for SNP discovery in genes with a large effect on production traits of commercial interest and for expression studies of growth and adaptation.

Molecular profiling of the gilthead sea bream (Sparus aurata L.) response to chronic exposure to the myxosporean parasite Enteromyxum leei

The aim of the present work was to investigate the transcriptome response of gilthead sea bream (Sparus aurata) after challenge with the myxosporean Enteromyxum leei, a wide-spread enteric parasite causing heavy economic losses in Mediterranean sparid farms. This parasite causes severe desquamative enteritis which usually leads to death of the fish, and there are no preventative or curative treatments for this enteromyxosis. After 113 days of exposure to parasite-contaminated effluent, fish were classified into three cohorts: control fish not exposed to parasite, those that were exposed and infected, and those that were exposed but not infected. In order to detect target genes that may be candidates for infective status or resistance, a cDNA microarray containing 18,490 cDNA clones enriched in genes differentially expressed after infection was hybridised with head kidney and intestine samples. In infected fish, 371 and 373 genes were differentially regulated at the >1.5-fold level in intestine and head kidney respectively, whereas in non-infected fish 175 and 501 genes were differentially regulated in these tissues, respectively. A global marked gene down-regulation was evident in infected fish, mainly in genes involved in the immune and acute phase response particularly complement and mannose binding lectin. Microarray analysis demonstrated a complex interplay between host and/or parasite derived proteases and protease inhibitors, apoptosis, cell proliferation and antioxidant defence genes in exposed fish. In the head kidney of non-infected fish a marked depression of genes involved in the acute phase response was evident. By contrast, in the intestine of non-infected fish, interferon-stimulated and MHC class II genes involved in antigen processing and presentation were up-regulated, possibly indicating that an active immune response at the local level is important to avoid infection with or proliferation of the parasite.

In silico mining and characterization of simple sequence repeats from gilthead sea bream (Sparus aurata) expressed sequence tags (EST-SSRs); PCR amplification, polymorphism evaluation and multiplexing and cross-species assays

We screened for simple sequence repeats (SSRs) found in ESTs derived from an EST-database development project ('Marine Genomics Europe' Network of Excellence). Different motifs of di-, tri-, tetra-, penta- and hexanucleotide SSRs were evaluated for variation in length and position in the expressed sequences, relative abundance and distribution in gilthead sea bream (Sparus aurata). We found 899 ESTs that harbor 997 SSRs (4.94%). On average, one SSR was found per 2.95 kb of EST sequence and the dinucleotide SSRs are the most abundant accounting for 47.6% of the total number. EST-SSRs were used as template for primer design. 664 primer pairs could be successfully identified and a subset of 206 pairs of primers was synthesized, PCR-tested and visualized on ethidium bromide stained agarose gels. The main objective was to further assess the potential of EST-SSRs as informative markers and investigate their cross-species amplification in sixteen teleost fish species: seven sparid species and nine other species from different families. Approximately 78% of the primer pairs gave PCR products of expected size in gilthead sea bream, and as expected, the rate of successful amplification of sea bream EST-SSRs was higher in sparids, lower in other perciforms and even lower in species of the Clupeiform and Gadiform orders. We finally determined the polymorphism and the heterozygosity of 63 markers in a wild gilthead sea bream population; fifty-eight loci were found to be polymorphic with the expected heterozygosity and the number of alleles ranging from 0.089 to 0.946 and from 2 to 27, respectively. These tools and markers are expected to enhance the available genetic linkage map in gilthead sea bream, to assist comparative mapping and genome analyses for this species and further with other model fish species and finally to help advance genetic analysis for cultivated and wild populations and accelerate breeding programs.