Design and performance of a turbot (Scophthalmus maximus) oligo-microarray based on ESTs from immune tissues

Analysis of apolipoprotein multigene family in spotted sea bass (Lateolabrax maculatus) and their expression profiles in response to Vibrio harveyi infection

Apolipoproteins (Apos), which are the protein components of plasma lipoproteins, play important roles in lipid transport in vertebrates. It has been demonstrated that in teleosts, several Apos display antimicrobial activity and play crucial roles in innate immunity. Despite their importance, apo genes have not been systematically characterized in many aquaculture fish species. In our study, a complete set of 23 apo genes was identified and annotated from spotted sea bass (Lateolabrax maculatus). Phylogenetic and homology analyses provided evidence for their annotation and evolutionary relationships. To investigate their potential roles in the immune response, the expression patterns of 23 apo genes were determined in the liver and intestine by qRT-PCR after Vibrio harveyi infection. After infection, a total of 20 differentially expressed apo genes were observed, and their expression profiles varied among the genes and tissues. 5 apo genes (apoA1, apoA4a.1, apoC2, apoF and apoO) were dramatically induced or suppressed (log2 fold change >4, P < 0.05), suggesting their involvement in the immune response of spotted sea bass. Our study provides a valuable foundation for future studies aimed at uncovering the specific roles of each apo gene during bacterial infection in spotted sea bass and other teleost species.

Characterization of the European Sea Bass (Dicentrarchus labrax) Gonadal Transcriptome During Sexual Development

The European sea bass is one of the most important cultured fish in Europe and has a marked sexual growth dimorphism in favor of females. It is a gonochoristic species with polygenic sex determination, where a combination between still undifferentiated genetic factors and environmental temperature determines sex ratios. The molecular mechanisms responsible for gonadal sex differentiation are still unknown. Here, we sampled fish during the gonadal developmental period (110 to 350 days post fertilization, dpf), and performed a comprehensive transcriptomic study by using a species-specific microarray. This analysis uncovered sex-specific gonadal transcriptomic profiles at each stage of development, identifying larger number of differentially expressed genes in ovaries when compared to testis. The expression patterns of 54 reproduction-related genes were analyzed. We found that hsd17β10 is a reliable marker of early ovarian differentiation. Further, three genes, pdgfb, snx1, and nfy, not previously related to fish sex differentiation, were tightly associated with testis development in the sea bass. Regarding signaling pathways, lysine degradation, bladder cancer, and NOD-like receptor signaling were enriched for ovarian development while eight pathways including basal transcription factors and steroid biosynthesis were enriched for testis development. Analysis of the transcription factor abundance showed an earlier increase in females than in males. Our results show that, although many players in the sex differentiation pathways are conserved among species, there are peculiarities in gene expression worth exploring. The genes identified in this study illustrate the diversity of players involved in fish sex differentiation and can become potential biomarkers for the management of sex ratios in the European sea bass and perhaps other cultured species.

Long-term affected flat oyster (Ostrea edulis) haemocytes show differential gene expression profiles from naïve oysters in response to Bonamia ostreae

European flat oyster (Ostrea edulis) production has suffered a severe decline due to bonamiosis. The responsible parasite enters in oyster haemocytes, causing an acute inflammatory response frequently leading to death. We used an immune-enriched oligo-microarray to understand the haemocyte response to Bonamia ostreae by comparing expression profiles between naïve (NS) and long-term affected (AS) populations along a time series (1 d, 30 d, 90 d). AS showed a much higher response just after challenge, which might be indicative of selection for resistance. No regulated genes were detected at 30 d in both populations while a notable reactivation was observed at 90 d, suggesting parasite latency during infection. Genes related to extracellular matrix and protease inhibitors, up-regulated in AS, and those related to histones, down-regulated in NS, might play an important role along the infection. Twenty-four candidate genes related to resistance should be further validated for selection programs aimed to control bonamiosis.

Gene expression analysis of Ruditapes philippinarum haemocytes after experimental Perkinsus olseni zoospore challenge and infection in the wild

The production of Manila clam (Ruditapes philippinarum) is seriously threatened by the protistan parasite Perkinsus olseni. We characterized and compared gene expression of Manila clam haemocytes in response to P. olseni in a time-course (10 h, 24 h, 8 d) controlled laboratory challenge (LC), representing the first step of infection, and in a more complex infection in the wild (WI), using a validated oligo-microarray containing 11,232 transcripts, mostly annotated. Several immune-genes involved in NIK/NF-kappaB signalling, Toll-like receptor signalling and apoptosis were activated at LC-10 h. However, down-regulation of genes encoding lysozyme, histones, cathepsins and heat shock proteins indicated signals of immunodepression, which persisted at LC-24 h, when only down-regulated genes were detected. A rebound of haemocyte activity occurred at LC-8 d as shown by up-regulation of genes involved in cytoskeleton organization and cell survival. The WI study showed a more complex picture, and several immune-relevant processes including cytoskeleton organization, cell survival, apoptosis, encapsulation, cell redox- and lipid-homeostasis were activated, illustrating the main mechanism of host response. Our results provide useful information, including potential biomarkers, to develop strategies for controlling Manila clam perkinsosis.

Integrating genomic resources of flatfish (Pleuronectiformes) to boost aquaculture production

Flatfish have a high market acceptance thus representing a profitable aquaculture production. The main farmed species is the turbot (Scophthalmus maximus) followed by Japanese flounder (Paralichthys olivaceous) and tongue sole (Cynoglossus semilaevis), but other species like Atlantic halibut (Hippoglossus hippoglossus), Senegalese sole (Solea senegalensis) and common sole (Solea solea) also register an important production and are very promising for farming. Important genomic resources are available for most of these species including whole genome sequencing projects, genetic maps and transcriptomes. In this work, we integrate all available genomic information of these species within a common framework, taking as reference the whole assembled genomes of turbot and tongue sole (>210× coverage). New insights related to the genetic basis of productive traits and new data useful to understand the evolutionary origin and diversification of this group were obtained. Despite a general 1:1 chromosome syntenic relationship between species, the comparison of turbot and tongue sole genomes showed huge intrachromosomic reorganizations. The integration of available mapping information supported specific chromosome fusions along flatfish evolution and facilitated the comparison between species of previously reported genetic associations for productive traits. When comparing transcriptomic resources of the six species, a common set of ~2500 othologues and ~150 common miRNAs were identified, and specific sets of putative missing genes were detected in flatfish transcriptomes, likely reflecting their evolutionary diversification.

Construction of an Ostrea edulis database from genomic and expressed sequence tags (ESTs) obtained from Bonamia ostreae infected haemocytes: Development of an immune-enriched oligo-microarray

The flat oyster, Ostrea edulis, is one of the main farmed oysters, not only in Europe but also in the United States and Canada. Bonamiosis due to the parasite Bonamia ostreae has been associated with high mortality episodes in this species. This parasite is an intracellular protozoan that infects haemocytes, the main cells involved in oyster defence. Due to the economical and ecological importance of flat oyster, genomic data are badly needed for genetic improvement of the species, but they are still very scarce. The objective of this study is to develop a sequence database, OedulisDB, with new genomic and transcriptomic resources, providing new data and convenient tools to improve our knowledge of the oyster's immune mechanisms. Transcriptomic and genomic sequences were obtained using 454 pyrosequencing and compiled into an O. edulis database, OedulisDB, consisting of two sets of 10,318 and 7159 unique sequences that represent the oyster's genome (WG) and de novo haemocyte transcriptome (HT), respectively. The flat oyster transcriptome was obtained from two strains (naïve and tolerant) challenged with B. ostreae, and from their corresponding non-challenged controls. Approximately 78.5% of 5619 HT unique sequences were successfully annotated by Blast search using public databases. A total of 984 sequences were identified as being related to immune response and several key immune genes were identified for the first time in flat oyster. Additionally, transcriptome information was used to design and validate the first oligo-microarray in flat oyster enriched with immune sequences from haemocytes. Our transcriptomic and genomic sequencing and subsequent annotation have largely increased the scarce resources available for this economically important species and have enabled us to develop an OedulisDB database and accompanying tools for gene expression analysis. This study represents the first attempt to characterize in depth the O. edulis haemocyte transcriptome in response to B. ostreae through massively sequencing and has aided to improve our knowledge of the immune mechanisms of flat oyster. The validated oligo-microarray and the establishment of a reference transcriptome will be useful for large-scale gene expression studies in this species.

Vaccine-induced modulation of gene expression in turbot peritoneal cells. A microarray approach

We used a microarray approach to examine changes in gene expression in turbot peritoneal cells after injection of the fish with vaccines containing the ciliate parasite Philasterides dicentrarchi as antigen and one of the following adjuvants: chitosan-PVMMA microspheres, Freund́s complete adjuvant, aluminium hydroxide gel or Matrix-Q (Isconova, Sweden). We identified 374 genes that were differentially expressed in all groups of fish. Forty-two genes related to tight junctions and focal adhesions and/or actin cytoskeleton were differentially expressed in free peritoneal cells. The profound changes in gene expression related to cell adherence and cytoskeleton may be associated with cell migration and also with the formation of cell-vaccine masses and their attachment to the peritoneal wall. Thirty-five genes related to apoptosis were differentially expressed. Although most of the proteins coded by these genes have a proapoptotic effect, others are antiapoptotic, indicating that both types of signals occur in peritoneal leukocytes of vaccinated fish. Interestingly, many of the genes related to lymphocytes and lymphocyte activity were downregulated in the groups injected with vaccine. We also observed decreased expression of genes related to antigen presentation, suggesting that macrophages (which were abundant in the peritoneal cavity after vaccination) did not express these during the early inflammatory response in the peritoneal cavity. Finally, several genes that participate in the inflammatory response were differentially expressed, and most participated in resolution of inflammation, indicating that an M2 macrophage response is generated in the peritoneal cavity of fish one day post vaccination.

A linkage map of transcribed single nucleotide polymorphisms in rohu (Labeo rohita) and QTL associated with resistance to Aeromonas hydrophila

BACKGROUND

Production of carp dominates world aquaculture. More than 1.1 million tonnes of rohu carp, Labeo rohita (Hamilton), were produced in 2010. Aeromonas hydrophila is a bacterial pathogen causing aeromoniasis in rohu, and is a major problem for carp production worldwide. There is a need to better understand the genetic mechanisms affecting resistance to this disease, and to develop tools that can be used with selective breeding to improve resistance. Here we use a 6 K SNP array to genotype 21 full-sibling families of L. rohita that were experimentally challenged intra-peritoneally with a virulent strain of A. hydrophila to scan the genome for quantitative trait loci associated with disease resistance.

RESULTS

In all, 3193 SNPs were found to be informative and were used to create a linkage map and to scan for QTL affecting resistance to A. hydrophila. The linkage map consisted of 25 linkage groups, corresponding to the number of haploid chromosomes in L. rohita. Male and female linkage maps were similar in terms of order, coverage (1384 and 1393 cM, respectively) and average interval distances (1.32 and 1.35 cM, respectively). Forty-one percent of the SNPs were annotated with gene identity using BLAST (cut off E-score of 0.001). Twenty-one SNPs mapping to ten linkage groups showed significant associations with the traits hours of survival and dead or alive (P <0.05 after Bonferroni correction). Of the SNPs showing significant or suggestive associations with the traits, several were homologous to genes of known immune function or were in close linkage to such genes. Genes of interest included heat shock proteins (70, 60, 105 and "small heat shock proteins"), mucin (5b precursor and 2), lectin (receptor and CD22), tributyltin-binding protein, major histocompatibility loci (I and II), complement protein component c7-1, perforin 1, ubiquitin (ligase, factor e4b isoform 2 and conjugation enzyme e2 c), proteasome subunit, T-cell antigen receptor and lymphocyte specific protein tyrosine kinase.

CONCLUSIONS

A panel of markers has been identified that will be validated for use with both genomic and marker-assisted selection to improve resistance of L. rohita to A. hydrophila.

Regional environmental pressure influences population differentiation in turbot (Scophthalmus maximus)

Unravelling the factors shaping the genetic structure of mobile marine species is challenging due to the high potential for gene flow. However, genetic inference can be greatly enhanced by increasing the genomic, geographical or environmental resolution of population genetic studies. Here, we investigated the population structure of turbot (Scophthalmus maximus) by screening 17 random and gene-linked markers in 999 individuals at 290 geographical locations throughout the northeast Atlantic Ocean. A seascape genetics approach with the inclusion of high-resolution oceanographical data was used to quantify the association of genetic variation with spatial, temporal and environmental parameters. Neutral loci identified three subgroups: an Atlantic group, a Baltic Sea group and one on the Irish Shelf. The inclusion of loci putatively under selection suggested an additional break in the North Sea, subdividing southern from northern Atlantic individuals. Environmental and spatial seascape variables correlated marginally with neutral genetic variation, but explained significant proportions (respectively, 8.7% and 10.3%) of adaptive genetic variation. Environmental variables associated with outlier allele frequencies included salinity, temperature, bottom shear stress, dissolved oxygen concentration and depth of the pycnocline. Furthermore, levels of explained adaptive genetic variation differed markedly between basins (3% vs. 12% in the North and Baltic Sea, respectively). We suggest that stable environmental selection pressure contributes to relatively strong local adaptation in the Baltic Sea. Our seascape genetic approach using a large number of sampling locations and associated oceanographical data proved useful for the identification of population units as the basis of management decisions.

Molecular analysis of endocrine disruption in hornyhead turbot at wastewater outfalls in southern california using a second generation multi-species microarray

Sentinel fish hornyhead turbot (Pleuronichthysverticalis) captured near wastewater outfalls are used for monitoring exposure to industrial and agricultural chemicals of ~ 20 million people living in coastal Southern California. Although analyses of hormones in blood and organ morphology and histology are useful for assessing contaminant exposure, there is a need for quantitative and sensitive molecular measurements, since contaminants of emerging concern are known to produce subtle effects. We developed a second generation multi-species microarray with expanded content and sensitivity to investigate endocrine disruption in turbot captured near wastewater outfalls in San Diego, Orange County and Los Angeles California. Analysis of expression of genes involved in hormone [e.g., estrogen, androgen, thyroid] responses and xenobiotic metabolism in turbot livers was correlated with a series of phenotypic end points. Molecular analyses of turbot livers uncovered altered expression of vitellogenin and zona pellucida protein, indicating exposure to one or more estrogenic chemicals, as well as, alterations in cytochrome P450 (CYP) 1A, CYP3A and glutathione S-transferase-α indicating induction of the detoxification response. Molecular responses indicative of exposure to endocrine disruptors were observed in field-caught hornyhead turbot captured in Southern California demonstrating the utility of molecular methods for monitoring environmental chemicals in wastewater outfalls. Moreover, this approach can be adapted to monitor other sites for contaminants of emerging concern in other fish species for which there are few available gene sequences.

A combined strategy involving Sanger and 454 pyrosequencing increases genomic resources to aid in the management of reproduction, disease control and genetic selection in the turbot (Scophthalmus maximus)

Background

Genomic resources for plant and animal species that are under exploitation primarily for human consumption are increasingly important, among other things, for understanding physiological processes and for establishing adequate genetic selection programs. Current available techniques for high-throughput sequencing have been implemented in a number of species, including fish, to obtain a proper description of the transcriptome. The objective of this study was to generate a comprehensive transcriptomic database in turbot, a highly priced farmed fish species in Europe, with potential expansion to other areas of the world, for which there are unsolved production bottlenecks, to understand better reproductive- and immune-related functions. This information is essential to implement marker assisted selection programs useful for the turbot industry.

Results

Expressed sequence tags were generated by Sanger sequencing of cDNA libraries from different immune-related tissues after several parasitic challenges. The resulting database (“Turbot 2 database”) was enlarged with sequences generated from a 454 sequencing run of brain-hypophysis-gonadal axis-derived RNA obtained from turbot at different development stages. The assembly of Sanger and 454 sequences generated 52,427 consensus sequences (“Turbot 3 database”), of which 23,661 were successfully annotated. A total of 1,410 sequences were confirmed to be related to reproduction and key genes involved in sex differentiation and maturation were identified for the first time in turbot (AR, AMH, SRY-related genes, CYP19A, ZPGs, STAR FSHR, etc.). Similarly, 2,241 sequences were related to the immune system and several novel key immune genes were identified (BCL, TRAF, NCK, CD28 and TOLLIP, among others). The number of genes of many relevant reproduction- and immune-related pathways present in the database was 50–90% of the total gene count of each pathway. In addition, 1,237 microsatellites and 7,362 single nucleotide polymorphisms (SNPs) were also compiled. Further, 2,976 putative natural antisense transcripts (NATs) including microRNAs were also identified.

Conclusions

The combined sequencing strategies employed here significantly increased the turbot genomic resources available, including 34,400 novel sequences. The generated database contains a larger number of genes relevant for reproduction- and immune-associated studies, with an excellent coverage of most genes present in many relevant physiological pathways. This database also allowed the identification of many microsatellites and SNP markers that will be very useful for population and genome screening and a valuable aid in marker assisted selection programs.

Microarray analysis of the inflammatory and immune responses in head kidney turbot leucocytes treated with resveratrol

A DNA oligo-microarray enriched in genes and involved in inflammatory and immune responses was used to evaluate the effects of resveratrol on gene expression in turbot head kidney leucocytes. Leucocytes were cultured for 3, 6 and 24 h, in the presence or absence of resveratrol, or were stimulated with the membrane fraction of the parasite Philasterides dicentrarchi or with the membrane plus resveratrol. Gene expression changed considerably in control cells, and several of the regulated genes were related to inflammatory and immune responses and to the cytoskeleton. Similar changes in gene expression occurred in control cells and in cells stimulated with P. dicentrarchi membrane fraction. Treatment with resveratrol induced changes in the expression (mostly down-regulation) of several genes involved in immune responses and inflammation. Thus, the down-regulation of the transcription factor PU.1, pentraxin-multidomain protein, heme oxygenase 1, S100 calcium-binding protein A-16 (S100A16) and the signal transducer and activator of transcription 4 was observed after all three incubation times. The down-regulation of the suppressor of cytokine signalling 3a, LPS-induced tumour necrosis alpha, hepcidin, metallothionein, TLR8 and the calcium dependent lectin A was observed after 3 and 6 h. Resveratrol also decreased the expression of CCL20, IL-8, apolipoprotein E and glutathione S-transferase after incubation for 6 and 24 h, and of TNF-α after incubation for 3 and 24 h. Resveratrol also induced strong regulation of several cytoskeleton-related genes. The use of the turbot oligo-microarray enabled us to discover genes whose expression was not previously suspected of being modulated by this polyphenol.

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

We evaluated the expression profiles of turbot in spleen, liver, and head kidney across five temporal points of the Philasterides dicentrarchi infection process using an 8x15K 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 221 differentially expressed (DE) genes (8.1% of the whole microarray), 113 in spleen, 83 in liver, and 90 in head kidney, in at least 1 of the 5 temporal points sampled for each organ. Most of these genes could be annotated (83.0%) and functionally categorized using GO 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. A high proportion of DE genes were organ-specific (70.6%), although their associated GO functions showed notable similarities in the three organs. The most striking difference in functional distribution was observed between the up- and downregulated gene groups. Upregulated genes were mostly associated to immune functions, while downregulated ones mainly involved metabolism-related genes. Genetic response appeared clustered in a few groups of genes with similar expression profiles along the temporal series. The information obtained will aid to understand the turbot immune response and will specifically be valuable to develop strategies of defense to P. dicentrarchi to achieve more resistant broodstocks for turbot industry.

Microarray-based identification of differentially expressed genes in families of turbot (Scophthalmus maximus) after infection with viral haemorrhagic septicaemia virus (VHSV)

Viral haemorrhagic septicaemia virus (VHSV) is one of the major threats to the development of the aquaculture industry worldwide. The present study was aimed to identify genes differentially expressed in several turbot (Scophthalmus maximus) families showing different mortality rates after VHSV. The expression analysis was conducted through genome-wide expression profiling with an oligo-microarray in the head kidney. A significant proportion of the variation in the gene expression profiles seemed to be explained by the genetic background, indicating that the mechanisms by which particular species and/or populations can resist a pathogen(s) are complex and multifactorial. Before the experimental infections, fish from resistant families (low mortality rates after VHSV infection) showed high expression of different antimicrobial peptides, suggesting that their pre-immune state may be stronger than fish of susceptible families (high mortality rates after VHSV infection). After infection, fish from both high- and low-mortality families showed an up-modulation of the interferon-induced Mx2 gene, the IL-8 gene and the VHSV-induced protein 5 gene compared with control groups. Low levels of several molecules secreted in the mucus were observed in high-mortality families, but different genes involved in viral entrance into target cells were down-regulated in low-mortality families. Moreover, these families also showed a strong down-modulation of marker genes related to VHSV target organs, including biochemical markers of renal dysfunction and myocardial injury. In general, the expression of different genes involved in the metabolism of sugars, lipids and proteins were decreased in both low- and high-mortality families after infection. The present study serves as an initial screen for genes of interest and provides an extensive overview of the genetic basis underlying the differences between families that are resistant or susceptible to VHSV infection.

Advances in genomics for flatfish aquaculture

Fish aquaculture is considered to be one of the most sustainable sources of protein for humans. Many different species are cultured worldwide, but among them, marine flatfishes comprise a group of teleosts of high commercial interest because of their highly prized white flesh. However, the aquaculture of these fishes is seriously hampered by the scarce knowledge on their biology. In recent years, various experimental 'omics' approaches have been applied to farmed flatfishes to increment the genomic resources available. These tools are beginning to identify genetic markers associated with traits of commercial interest, and to unravel the molecular basis of different physiological processes. This article summarizes recent advances in flatfish genomics research in Europe. We focus on the new generation sequencing technologies, which can produce a massive amount of DNA sequencing data, and discuss their potentials and applications for de novo genome sequencing and transcriptome analysis. The relevance of these methods in nutrigenomics and foodomics approaches for the production of healthy animals, as well as high quality and safety products for the consumer, is also briefly discussed.

An Expressed Sequence Tag (EST)-enriched genetic map of turbot (Scophthalmus maximus): a useful framework for comparative genomics across model and farmed teleosts

Background

The turbot (Scophthalmus maximus) is a relevant species in European aquaculture. The small turbot genome provides a source for genomics strategies to use in order to understand the genetic basis of productive traits, particularly those related to sex, growth and pathogen resistance. Genetic maps represent essential genomic screening tools allowing to localize quantitative trait loci (QTL) and to identify candidate genes through comparative mapping. This information is the backbone to develop marker-assisted selection (MAS) programs in aquaculture. Expressed sequenced tag (EST) resources have largely increased in turbot, thus supplying numerous type I markers suitable for extending the previous linkage map, which was mostly based on anonymous loci. The aim of this study was to construct a higher-resolution turbot genetic map using EST-linked markers, which will turn out to be useful for comparative mapping studies.

Results

A consensus gene-enriched genetic map of the turbot was constructed using 463 SNP and microsatellite markers in nine reference families. This map contains 438 markers, 180 EST-linked, clustered at 24 linkage groups. Linkage and comparative genomics evidences suggested additional linkage group fusions toward the consolidation of turbot map according to karyotype information. The linkage map showed a total length of 1402.7 cM with low average intermarker distance (3.7 cM; ~2 Mb). A global 1.6:1 female-to-male recombination frequency (RF) ratio was observed, although largely variable among linkage groups and chromosome regions. Comparative sequence analysis revealed large macrosyntenic patterns against model teleost genomes, significant hits decreasing from stickleback (54%) to zebrafish (20%). Comparative mapping supported particular chromosome rearrangements within Acanthopterygii and aided to assign unallocated markers to specific turbot linkage groups.

Conclusions

The new gene-enriched high-resolution turbot map represents a useful genomic tool for QTL identification, positional cloning strategies, and future genome assembling. This map showed large synteny conservation against model teleost genomes. Comparative genomics and data mining from landmarks will provide straightforward access to candidate genes, which will be the basis for genetic breeding programs and evolutionary studies in this species.

High-throughput sequence analysis of turbot (Scophthalmus maximus) transcriptome using 454-pyrosequencing for the discovery of antiviral immune genes

Background

Turbot (Scophthalmus maximus L.) is an important aquacultural resource both in Europe and Asia. However, there is little information on gene sequences available in public databases. Currently, one of the main problems affecting the culture of this flatfish is mortality due to several pathogens, especially viral diseases which are not treatable. In order to identify new genes involved in immune defense, we conducted 454-pyrosequencing of the turbot transcriptome after different immune stimulations.

Methodology/Principal Findings

Turbot were injected with viral stimuli to increase the expression level of immune-related genes. High-throughput deep sequencing using 454-pyrosequencing technology yielded 915,256 high-quality reads. These sequences were assembled into 55,404 contigs that were subjected to annotation steps. Intriguingly, 55.16% of the deduced protein was not significantly similar to any sequences in the databases used for the annotation and only 0.85% of the BLASTx top-hits matched S. maximus protein sequences. This relatively low level of annotation is possibly due to the limited information for this specie and other flatfish in the database. These results suggest the identification of a large number of new genes in turbot and in fish in general. A more detailed analysis showed the presence of putative members of several innate and specific immune pathways.

Conclusions/Significance

To our knowledge, this study is the first transcriptome analysis using 454-pyrosequencing for turbot. Previously, there were only 12,471 EST and less of 1,500 nucleotide sequences for S. maximus in NCBI database. Our results provide a rich source of data (55,404 contigs and 181,845 singletons) for discovering and identifying new genes, which will serve as a basis for microarray construction, gene expression characterization and for identification of genetic markers to be used in several applications. Immune stimulation in turbot was very effective, obtaining an enormous variety of sequences belonging to genes involved in the defense mechanisms.

Gene expression profiling of Atlantic cod (Gadus morhua) embryogenesis using microarray

Atlantic cod (Gadus morhua) is a fish species of high importance, as a key species in a range of Northern ecosystems, in fisheries, and as an emerging species in aquaculture. So far, little is known about the transcriptional activity during early developmental stages of Atlantic cod. Hence, we decided to use a cDNA microarray covering 7,000 genes to analyze the temporal activity of the transcriptome during cod embryogenesis. Twelve different embryonic time points were selected, covering major developmental stages and processes such as maternally derived mRNA, blastula, gastrula, segmentation, hatching, and first-feeding larval stage. The microarray analysis revealed a highly dynamic transcriptional profile, showing for the first time the differential expression of thousands of known and unknown genes during Atlantic cod embryogenesis. These initial findings will serve as an important baseline for future in-depth studies of candidate genes involved in development, reproductive control, disease resistance, growth, nutrient digestion, and metabolism.

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.

QTL detection for Aeromonas salmonicida resistance related traits in turbot (Scophthalmus maximus)

Background

Interactions between fish and pathogens, that may be harmless under natural conditions, often result in serious diseases in aquaculture systems. This is especially important due to the fact that the strains used in aquaculture are derived from wild strains that may not have had enough time to adapt to new disease pressures. The turbot is one of the most promising European aquaculture species. Furunculosis, caused by the bacterium Aeromonas salmonicida, produces important losses to turbot industry. An appealing solution is to achieve more robust broodstock, which can prevent or diminish the devastating effects of epizooties. Genomics strategies have been developed in turbot to look for candidate genes for resistance to furunculosis and a genetic map with appropriate density to screen for genomic associations has been also constructed. In the present study, a genome scan for QTL affecting resistance and survival to A. salmonicida in four turbot families was carried out. The objectives were to identify consistent QTL using different statistical approaches (linear regression and maximum likelihood) and to locate the tightest associated markers for their application in genetic breeding strategies.

Results

Significant QTL for resistance were identified by the linear regression method in three linkage groups (LGs 4, 6 and 9) and for survival in two LGs (6 and 9). The maximum likelihood methodology identified QTL in three LGs (5, 6 and 9) for both traits. Significant association between disease traits and genotypes was detected for several markers, some of them explaining up to 17% of the phenotypic variance. We also identified candidate genes located in the detected QTL using data from previously mapped markers.

Conclusions

Several regions controlling resistance to A. salmonicida in turbot have been detected. The observed concordance between different statistical methods at particular linkage groups gives consistency to our results. The detected associated markers could be useful for genetic breeding strategies. A finer mapping will be necessary at the detected QTL intervals to narrow associations and around the closely associated markers to look for candidate genes through comparative genomics or positional cloning strategies. The identification of associated variants at specific genes will be essential, together with the QTL associations detected in this study, for future marker assisted selection programs.

Comparative expression analysis in mature gonads, liver and brain of turbot (Scophthalmus maximus) by cDNA-AFLPS

Turbot is one of the most important farmed fish in Europe. This species exhibits a considerable sexual dimorphism in growth and sexual maturity that makes the all-female production recommended for turbot farming. Our knowledge about the genetic basis of sex determination and the molecular regulation of gonad differentiation in this species is still limited. Our goal was to identify and compare gene expression and functions between testes and ovaries in adults in order to ascertain the relationship between the genes that could be involved in the gonad differentiation or related to the sex determination system. The identification of differentially expressed sex related genes is an initial step towards understanding the molecular mechanisms of gonad differentiation. For this, we carried out a transcriptome analysis based on cDNA-AFLP technique which allowed us to obtain an initial frame on sex-specific gene expression that will facilitate further analysis especially along the critical gonad differentiating period. With the aim of widening the study on sex-biased gene expression we reproduced the same experiments in two somatic tissues: liver and brain. We have selected the liver because it is the most analyzed one regarding sexual dimorphic gene expression and due to its importance in steroid hormones metabolism and the brain because the functional relationship between brain and gonad is documented. We found slight but important differences between sexes which deserve further investigation.

The search for alternative aqueous film forming foams (AFFF) with a low environmental impact: physiological and transcriptomic effects of two Forafac(®) fluorosurfactants in turbot

Fluorosurfactants are the key components in aqueous film forming foams (AFFF). They provide these fire fighting agents with the required low surface tension and they enable film formation on top of lighter fuels to prevent burn back. Development of effective and environmentally acceptable PFOS alternatives is one of the most important priorities in the fire fighting foam industry. DuPont™ offers the fluorosurfactant mixtures Forafac(®)1157 and Forafac(®)1157N for the formulation of AFFFs which are alternatives to the persistent and toxic perfluorooctane sulphonate (PFOS). Ecotoxicological testing of these inadequately documented mixtures is necessary to include them in AFFF hazard and risk assessment. Juvenile turbot (Scophthalmus maximus) were exposed for 14 days to 0.1; 0.5 and 1.5mg/L of the fluorosurfactant mixtures used in Forafac(®)1157 and Forafac(®)1157N. In an initial transcriptomics experiment, microarray analysis revealed differentially expressed transcripts of genes which were mainly involved in digestion and in the immune system. This discovery-driven screening approach offered the basis for new hypotheses that were tested in two subsequent experiments in which food intake, energy reserves, growth and a set of haematological parameters were examined. Additionally, effects of the two mixtures were compared to those of PFOS. Based on the results of this study, the mode of action of Forafac(®)1157N was the activation of the acute phase reaction resulting in increased leukocyte concentrations and the inhibition of growth due to the high energetic cost of toxicant exposure. For Forafac(®)1157, evidences of immunosuppression were found on the transcriptional level and the altered differential leukocyte profiles indicated that stress was induced in these fish. However, food intake, energy reserves and growth were not compromised, even at high exposure concentrations, which was in contrast to the effects seen after PFOS exposure. Taking into account that Forafac(®)1157 appeared to be less toxic than PFOS, this mixture could be considered as a more environmentally acceptable PFOS alternative for the use in AFFFs.

Lipovitellin is a non-self recognition receptor with opsonic activity

Lipovitellin (Lv), a glycolipoprotein, is a major component of the egg yolk, which is usually regarded as an energy reserve of nutrients essential for growth and development. We have purified Lv from ovulated eggs of the rosy barb Puntius conchonius by two-step chromatography and characterized it by staining with periodic acid/Schiff reagent and Sudan black B, amino acid composition analysis, and peptide mass fingerprinting. The results of ligand and bacterial binding assays, an enzyme-linked immunosorbent assay (ELISA), and the phagocytosis test revealed, for the first time, that the purified native form of P. conchonius Lv acts as a pattern recognition molecule with multiple specificities capable of identifying pathogen-associated molecular patterns (PAMPs), including those of lipopolysaccharide, lipoteichoic acid, and peptidoglycan, rather than self components and that it can bind Gram-negative and -positive bacteria, such as Escherichia coli and Staphylococcus aureus. These tests also showed that the P. conchonius Lv functions as an opsonin capable of enhancing macrophage phagocytosis. Taken together, these characteristics suggest that in developing embryos/larvae of P. conchonius, the native form of Lv may be physiologically involved in the sensing of invading pathogens via interaction with PAMPs and in the recruitment of the primitive macrophages that appear in early embryos to phagocytose and digest the pathogens, thereby protecting them from pathogenic attacks.

Transcriptome sequencing and microarray development for the Manila clam, Ruditapes philippinarum: genomic tools for environmental monitoring

BACKGROUND

The Manila clam, Ruditapes philippinarum, is one of the major aquaculture species in the world and a potential sentinel organism for monitoring the status of marine ecosystems. However, genomic resources for R. philippinarum are still extremely limited. Global analysis of gene expression profiles is increasingly used to evaluate the biological effects of various environmental stressors on aquatic animals under either artificial conditions or in the wild. Here, we report on the development of a transcriptomic platform for global gene expression profiling in the Manila clam.

RESULTS

A normalized cDNA library representing a mixture of adult tissues was sequenced using a ultra high-throughput sequencing technology (Roche 454). A database consisting of 32,606 unique transcripts was constructed, 9,747 (30%) of which could be annotated by similarity. An oligo-DNA microarray platform was designed and applied to profile gene expression of digestive gland and gills. Functional annotation of differentially expressed genes between different tissues was performed by enrichment analysis. Expression of Natural Antisense Transcripts (NAT) analysis was also performed and bi-directional transcription appears a common phenomenon in the R. philippinarum transcriptome. A preliminary study on clam samples collected in a highly polluted area of the Venice Lagoon demonstrated the applicability of genomic tools to environmental monitoring.

CONCLUSIONS

The transcriptomic platform developed for the Manila clam confirmed the high level of reproducibility of current microarray technology. Next-generation sequencing provided a good representation of the clam transcriptome. Despite the known limitations in transcript annotation and sequence coverage for non model species, sufficient information was obtained to identify a large set of genes potentially involved in cellular response to environmental stress.

Genomic resources for flatfish research and their applications

Flatfishes are a group of teleosts of high commercial and environmental interest, whose biology is still poorly understood. The recent rapid development of different 'omic' technologies is, however, enhancing the knowledge of the complex genetic control underlying different physiological processes of flatfishes. This review describes the different functional genomic approaches and resources currently available for flatfish research and summarizes different areas where microarray-based gene expression analysis has been applied. The increase in genome sequencing data has also allowed the construction of genetic linkage maps in different flatfish species; these maps are invaluable for investigating genome organization and identifying genetic traits of commercial interest. Despite the significant progress in this field, the genomic resources currently available for flatfish are still scarce. Further intensive research should be carried out to develop larger genomic sequence databases, high-density microarrays and, more detailed, complete linkage maps, using second-generation sequencing platforms. These tools will be crucial for further expanding the knowledge of flatfish physiology, and it is predicted that they will have important implications for wild fish population management, improved fish welfare and increased productivity in aquaculture.

Development of an oligo DNA microarray for the European sea bass and its application to expression profiling of jaw deformity

Background

The European sea bass (Dicentrarchus labrax) is a marine fish of great importance for fisheries and aquaculture. Functional genomics offers the possibility to discover the molecular mechanisms underlying productive traits in farmed fish, and a step towards the application of marker assisted selection methods in this species. To this end, we report here on the development of an oligo DNA microarray for D. labrax.

Results

A database consisting of 19,048 unique transcripts was constructed, of which 12,008 (63%) could be annotated by similarity and 4,692 received a GO functional annotation. Two non-overlapping 60mer probes were designed for each unique transcript and in-situ synthesized on glass slides using Agilent SurePrint™ technology. Probe design was positively completed for 19,035 target clusters; the oligo microarray was then applied to profile gene expression in mandibles and whole-heads of fish affected by prognathism, a skeletal malformation that strongly affects sea bass production. Statistical analysis identified 242 transcripts that are significantly down-regulated in deformed individuals compared to normal fish, with a significant enrichment in genes related to nervous system development and functioning. A set of genes spanning a wide dynamic range in gene expression level were selected for quantitative RT-PCR validation. Fold change correlation between microarray and qPCR data was always significant.

Conclusions

The microarray platform developed for the European sea bass has a high level of flexibility, reliability, and reproducibility. Despite the well known limitations in achieving a proper functional annotation in non-model species, sufficient information was obtained to identify biological processes that are significantly enriched among differentially expressed genes. New insights were obtained on putative mechanisms involved on mandibular prognathism, suggesting that bone/nervous system development might play a role in this phenomenon.