Genomic resources for flatfish research and their applications

Genome sequencing and analysis of black flounder (Paralichthys orbignyanus) reveals new insights into Pleuronectiformes genomic size and structure

Black flounder (Paralichthys orbignyanus, Pleuronectiformes) is a commercially significant marine fish with promising aquaculture potential in Argentina. Despite extensive studies on Black flounder aquaculture, its limited genetic information available hampers the crucial role genetics plays in the development of this activity. In this study, we first employed Illumina sequencing technology to sequence the entire genome of Black flounder. Utilizing two independent libraries-one from a female and another from a male-with 150 bp paired-end reads, a mean insert length of 350 bp, and over 35 X-fold coverage, we achieved assemblies resulting in a genome size of ~ 538 Mbp. Analysis of the assemblies revealed that more than 98% of the core genes were present, with more than 78% of them having more than 50% coverage. This indicates a somehow complete and accurate genome at the coding sequence level. This genome contains 25,231 protein-coding genes, 445 tRNAs, 3 rRNAs, and more than 1,500 non-coding RNAs of other types. Black flounder, along with pufferfishes, seahorses, pipefishes, and anabantid fish, displays a smaller genome compared to most other teleost groups. In vertebrates, the number of transposable elements (TEs) is often correlated with genome size. However, it remains unclear whether the sizes of introns and exons also play a role in determining genome size. Hence, to elucidate the potential factors contributing to this reduced genome size, we conducted a comparative genomic analysis between Black flounder and other teleost orders to determine if the small genomic size could be explained by repetitive elements or gene features, including the whole genome genes and introns sizes. We show that the smaller genome size of flounders can be attributed to several factors, including changes in the number of repetitive elements, and decreased gene size, particularly due to lower amount of very large and small introns. Thus, these components appear to be involved in the genome reduction in Black flounder. Despite these insights, the full implications and potential benefits of genome reduction in Black flounder for reproduction and aquaculture remain incompletely understood, necessitating further research.

Insights into the Mitochondrial Genetic Makeup and Miocene Colonization of Primitive Flatfishes (Pleuronectiformes: Psettodidae) in the East Atlantic and Indo-West Pacific Ocean

The mitogenomic evolution of the Psettodes flatfishes is still poorly known from their range distribution in eastern Atlantic and Indo-West Pacific Oceans. The study delves into the matrilineal evolutionary pathway of these primitive flatfishes, with a specific focus on the complete mitogenome of the Psettodes belcheri species, as determined through next-generation sequencing. The mitogenome in question spans a length of 16,747 base pairs and comprises a total of 37 genes, including 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region. Notably, the mitogenome of P. belcheri exhibits a bias towards AT base pairs, with a composition of 54.15%, mirroring a similar bias observed in its close relative, Psettodes erumei, which showcases percentages of 53.07% and 53.61%. Most of the protein-coding genes commence with an ATG initiation codon, except for Cytochrome c oxidase I (COI), which initiates with a GTG codon. Additionally, four protein-coding genes commence with a TAA termination codon, while seven others exhibit incomplete termination codons. Furthermore, two protein-coding genes, namely NAD1 and NAD6, terminate with AGG and TAG stop codons, respectively. In the mitogenome of P. belcheri, the majority of transfer RNAs demonstrate the classical cloverleaf secondary structures, except for tRNA-serine, which lacks a DHU stem. Comparative analysis of conserved blocks within the control regions of two Psettodidae species unveiled that the CSB-II block extended to a length of 51 base pairs, surpassing the other blocks and encompassing highly variable sites. A comprehensive phylogenetic analysis using mitochondrial genomes (13 concatenated PCGs) categorized various Pleuronectiformes species, highlighting the basal position of the Psettodidae family and showed monophyletic clustering of Psettodes species. The approximate divergence time (35-10 MYA) between P. belcheri and P. erumei was estimated, providing insights into their separation and colonization during the early Miocene. The TimeTree analysis also estimated the divergence of two suborders, Psettodoidei and Pleuronectoidei, during the late Paleocene to early Eocene (56.87 MYA). The distribution patterns of Psettodes flatfishes were influenced by ocean currents and environmental conditions, contributing to their ecological speciation. In the face of climate change and anthropogenic activities, the conservation implications of Psettodes flatfishes are emphasized, underscoring the need for regulated harvesting and adaptive management strategies to ensure their survival in changing marine ecosystems. Overall, this study contributes to understanding the evolutionary history, genetic diversity, and conservation needs of Psettodes flatfishes globally. However, the multifaceted exploration of mitogenome and larger-scale genomic data of Psettodes flatfish will provide invaluable insights into their genetic characterization, evolutionary history, environmental adaptation, and conservation in the eastern Atlantic and Indo-West Pacific Oceans.

Omics and cytokine discovery in fish: Presenting the Yellowtail kingfish (Seriola lalandi) as a case study

A continued programme of research is essential to overcome production bottlenecks in any aquacultured fish species. Since the introduction of genetic and molecular techniques, the quality of immune research undertaken in fish has greatly improved. Thousands of species specific cytokine genes have been discovered, which can be used to conduct more sensitive studies to understand how fish physiology is affected by aquaculture environments or disease. Newly available transcriptomic technologies, make it increasingly easier to study the immunogenetics of farmed species for which little data exists. This paper reviews how the application of transcriptomic procedures such as RNA Sequencing (RNA-Seq) can advance fish research. As a case study, we present some preliminary findings using RNA-Seq to identify cytokine related genes in Seriola lalandi. These will allow in-depth investigations to understand the immune responses of these fish in response to environmental change or disease and help in the development of therapeutic approaches.

Comprehensive transcriptomic analysis of the process of gonadal sex differentiation in the turbot (Scophthalmus maximus)

The turbot is a flatfish with a ZW/ZZ sex determination system but with a still unknown sex determining gene(s), and with a marked sexual growth dimorphism in favor of females. To better understand sexual development in turbot we sampled young turbot encompassing the whole process of gonadal differentiation and conducted a comprehensive transcriptomic study on its sex differentiation using a validated custom oligomicroarray. Also, the expression profiles of 18 canonical reproduction-related genes were studied along gonad development. The expression levels of gonadal aromatase cyp19a1a alone at three months of age allowed the accurate and early identification of sex before the first signs of histological differentiation. A total of 56 differentially expressed genes (DEG) that had not previously been related to sex differentiation in fish were identified within the first three months of age, of which 44 were associated with ovarian differentiation (e.g., cd98, gpd1 and cry2), and 12 with testicular differentiation (e.g., ace, capn8 and nxph1). To identify putative sex determining genes, ∼4.000 DEG in juvenile gonads were mapped and their positions compared with that of previously identified sex- and growth-related quantitative trait loci (QTL). Although no genes mapped to the previously identified sex-related QTLs, two genes (foxl2 and 17βhsd) of the canonical reproduction-related genes mapped to growth-QTLs in linkage group (LG) 15 and LG6, respectively, suggesting that these genes are related to the growth dimorphism in this species.

First haploid genetic map based on microsatellite markers in Senegalese sole (Solea senegalensis, Kaup 1858)

The Senegalese sole (Solea senegalensis, Kaup 1858) is a flatfish species of great value for aquaculture. In this study, we develop the first linkage map in this species based on microsatellite markers characterized from genomic DNA libraries and EST databases of Senegalese sole and from other flatfish species. Three reference gynogenetic families were obtained by chromosome-manipulation techniques: two haploid gynogenetics, used to assign and order microsatellites to linkage groups and another diploid gynogenetic family, used for estimating marker-centromere distances. The consensus map consists of 129 microsatellites distributed in 27 linkage groups (LG), with an average density of 4.7 markers per LG and comprising 1,004 centimorgans (cM). Additionally, 15 markers remained unlinked. Through half-tetrad analysis, we were able to estimate the centromere distance for 81 markers belonging to 24 LG, representing an average of 3 markers per LG. Comparative mapping was performed between flatfish species LG and model fish species chromosomes (stickleback, Tetraodon, medaka, fugu and zebrafish). The usefulness of microsatellite markers and the genetic map as tools for comparative mapping and evolution studies is discussed.

A genetic linkage map of sole (Solea solea): a tool for evolutionary and comparative analyses of exploited (flat)fishes

Linkage maps based on markers derived from genes are essential evolutionary tools for commercial marine fish to help identify genomic regions associated with complex traits and subject to selective forces at play during exploitation or selective breeding. Additionally, they allow the use of genomic information from other related species for which more detailed information is available. Sole (solea solea L.) is a commercially important flatfish species in the North Sea, subject to overexploitation and showing evidence of fisheries-induced evolutionary changes in growth- and maturation-related traits. Sole would definitely benefit from a linkage map to better understand how evolution has shaped its genome structure. This study presents a linkage map of sole based on 423 single nucleotide polymorphisms derived from expressed sequence tags and 8 neutral microsatellite markers. The total map length is 1233.8 cM and consists of 38 linkage groups with a size varying between 0 to 92.1 cM. Being derived from expressed sequence tags allowed us to align the map with the genome of four model fish species, namely medaka (Oryzias latipes), Nile tilapia (Oreochromis niloticus), three-spined stickleback (Gasterosteus aculeatus) and green spotted pufferfish (Tetraodon nigroviridis). This comparison revealed multiple conserved syntenic regions with all four species, and suggested that the linkage groups represent 21 putative sole chromosomes. The map was also compared to the linkage map of turbot (Scophthalmus maximus), another commercially important flatfish species and closely related to sole. For all putative sole chromosomes (except one) a turbot homolog was detected, confirming the even higher degree of synteny between these two flatfish species.

Microarray analysis of hepatic gene expression in juvenile Japanese flounder Paralichthys olivaceus fed diets supplemented with fish or vegetable oils

Gene expression profiling was performed in Japanese flounder Paralichthys olivaceus fed diets supplemented with fish oil (FO), linseed oil (LO), or olive oil (OO) for 6 weeks. The LO and OO groups showed significantly retarded growth, lower feed intake, lower protein efficiency ratio, and lower hepatosomatic index (P < 0.05). Liver fatty acid composition reflected the dietary fatty acid composition. Microarray analysis revealed that dietary n - 3 highly unsaturated fatty acid (HUFA) deficiency affected 169 transcripts. In the LO group, 57 genes were up-regulated and 38 genes were down-regulated, whereas in the OO group nine genes were up-regulated and 87 genes were down-regulated. Analysis of the functional annotations suggested that dietary n - 3 HUFA affected genes involved in signal transduction (23.2 %), cellular processes (21.1 %), metabolism (including glucose, lipid, and nucleobase; 15.5 %), transport (11.3 %), regulation of transcription (10.5 %), and immune response (4.2 %). Several genes encoding serine/threonine kinases such as protein kinase C and calmodulin-dependent kinase and nuclear hormone receptors such as vitamin D receptor, retinoic acid receptor, and receptors for cytokines (bone morphogenic protein and transforming growth factor β) were affected. Among 169 transcripts, 22 genes were affected in both LO and OO groups. The present study identified several genes involved in n - 3 HUFA deficiency-sensitive pathways, which will be useful for selective breeding of flounder strains able to adapt to n - 3 HUFA deficiency.

Stage-specific gene expression during spermatogenesis in the Senegalese sole (Solea senegalensis), a fish with semi-cystic type of spermatogenesis, as assessed by laser capture microdissection and absolute quantitative PCR

Spermatogenesis is a complex process where hormonal signals regulate the interaction of different cell types in a tight spatial and temporal fashion. The Senegalese sole (Solea senegalensis) is a marine flatfish that, in contrast to many fish, exhibits a semi-cystic, asynchronous pattern of spermatogenesis progression. This pattern is characterized by the release of spermatids into the tubule lumen, where they transform into spermatozoa. In this study, we used laser capture microdissection (LCM) to isolate cells from cysts containing spermatogonia, spermatocytes, spermatids or spermatozoa in order to investigate developmental patterns of gene expression. Furthermore, we also analyzed the stage-specific expression of the same set of genes throughout spermatogenesis (early-mid, late and maturing spermatogenic stages) in tissue fragments of the Senegalese sole testis. Genes analyzed by absolute qPCR in cysts isolated by LCM and stage-specific testis samples included genes involved in steroid synthesis and action (3β-hsd, 17β-hsd, 20β-hsd, star, star-like, progesterone receptor), gonadotropin action (fshr, lhr), the kisspeptin system (kiss2, kiss2r) and other genes important for the production of mature gametes (zona pellucida 2.2, claudin and clusterin). Our results show that, in general, steroidogenesis-related genes tended to increase with spermatogenesis progression and that 3β-hsd and 20β-hsd were expressed in germ cells but 17β-hsd was not. Our results also show that fshr is expressed in most testicular cell types, including germ cells. In contrast, lhr is expressed only in late spermatogenesis and is not expressed in any of the germ cell types examined, indicating that, in contrast to fshr, lhr may be primarily expressed in non-germinal cells (e.g. Leydig cells). Furthermore, kisspeptin and its receptor were expressed in all germ cell types examined and, as expected, gamete maturation-related genes were more expressed in mature stages. These results illustrate that key factors that participate in the hormonal regulation of spermatogenesis in the Senegalese sole testis show complex cell type- and stage-specific patterns of gene expression.

Population genetics provides new insights into biomarker prevalence in dab (Limanda limanda L.): a key marine biomonitoring species

Bioindicators are species for which some quantifiable aspect of its biology, a biomarker, is assumed to be sensitive to ecosystem health. However, there is frequently a lack of information on the underlying genetic and environmental drivers shaping the spatiotemporal variance in prevalence of the biomarkers employed. Here, we explore the relative role of potential variables influencing the spatiotemporal prevalence of biomarkers in dab, Limanda limanda, a species used as a bioindicator of marine contaminants. Firstly, the spatiotemporal genetic structure of dab around UK waters (39 samples across 15 sites for four years: 2005–2008) is evaluated with 16 microsatellites. Two temporally stable groups are identified corresponding to the North and Irish Seas (average between basin = 0.007; = 0.022). Secondly, we examine the association between biomarker prevalence and several variables, including genetic structuring, age and contaminant exposure. Genetic structure had significant interactive effects, together with age and some contaminants, in the prevalence of some of the biomarkers considered, namely hyperpigmentation and liver lesions. The integration of these data sets enhanced our understanding of the relationship between biomarker prevalence, exposure to contaminants and population-specific response, thereby yielding more informative predictive models of response and prospects for environmental remediation.

Exploring the larval transcriptome of the common sole (Solea solea L.)

BACKGROUND

The common sole (Solea solea) is a promising candidate for European aquaculture; however, the limited knowledge of the physiological mechanisms underlying larval development in this species has hampered the establishment of successful flatfish aquaculture. Although the fact that genomic tools and resources are available for some flatfish species, common sole genomics remains a mostly unexplored field. Here, we report, for the first time, the sequencing and characterisation of the transcriptome of S. solea and its application for the study of molecular mechanisms underlying physiological and morphological changes during larval-to-juvenile transition.

RESULTS

The S. solea transcriptome was generated from whole larvae and adult tissues using the Roche 454 platform. The assembly process produced a set of 22,223 Isotigs with an average size of 726 nt, 29 contigs and a total of 203,692 singletons. Of the assembled sequences, 75.2% were annotated with at least one known transcript/protein; these transcripts were then used to develop a custom oligo-DNA microarray. A total of 14,674 oligonucleotide probes (60 nt), representing 12,836 transcripts, were in situ synthesised onto the array using Agilent non-contact ink-jet technology. The microarray platform was used to investigate the gene expression profiles of sole larvae from hatching to the juvenile form. Genes involved in the ontogenesis of the visual system are up-regulated during the early stages of larval development, while muscle development and anaerobic energy pathways increase in expression over time. The gene expression profiles of key transcripts of the thyroid hormones (TH) cascade and the temporal regulation of the GH/IGF1 (growth hormone/insulin-like growth factor I) system suggest a pivotal role of these pathways in fish growth and initiation of metamorphosis. Pre-metamorphic larvae display a distinctive transcriptomic landscape compared to previous and later stages. Our findings highlighted the up-regulation of gene pathways involved in the development of the gastrointestinal system as well as biological processes related to folic acid and retinol metabolism. Additional evidence led to the formation of the hypothesis that molecular mechanisms of cell motility and ECM adhesion may play a role in tissue rearrangement during common sole metamorphosis.

CONCLUSIONS

Next-generation sequencing provided a good representation of the sole transcriptome, and the combination of different approaches led to the annotation of a high number of transcripts. The construction of a microarray platform for the characterisation of the larval sole transcriptome permitted the definition of the main processes involved in organogenesis and larval growth.

Multi-organ histological observations on juvenile Senegalese soles exposed to low concentrations of waterborne cadmium

A histopathological screening was performed on juvenile Senegalese soles exposed to environmentally realistic concentrations of waterborne Cd (0.5, 5 and 10 μg L(-1)) for 28 days. The severity and dissemination of histopathological changes were variable and limited to the kidney, liver, spleen, gills and skin goblet cells. Contradicting available literature that refers the liver as the most affected organ upon acute exposure and the kidney following chronic exposure, the liver was the most impacted organ (even at the lowest concentration), in a trend that could relate to the duration of exposure and Cd concentration. The most noticeable hepatic alterations related to inflammation, although hepatocellular alterations like lipidosis and eosinophilic foci also occurred. The trunk kidney of exposed fish endured moderate inflammation, apoptosis and necrosis, however, without a clear time-dependent effect. The spleen of fish subjected to the highest concentrations revealed diffuse necrotic foci accompanied by melanomacrophage intrusion. The gills, albeit the most important apical uptake organ of dissolved toxicants, sustained only moderate damage, from epithelial hyperplasia and pavement cell detachment to the potentially more severe chloride cell alterations. In the skin, an increase in goblet cell size occurred, most notoriously correlated to Cd concentration at earlier stages of exposure. The results show that a metal-naïve juvenile fish can endure deleterious effects when exposed to low, ecologically relevant, concentrations of a common toxic metal and that the pattern of Cd-induced histopathological alterations can be complex and linked to organ-specific responses and metal translocation within the organism.

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-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.

Hepatic proteome changes in Solea senegalensis exposed to contaminated estuarine sediments: a laboratory and in situ survey

Assessing toxicity of contaminated estuarine sediments poses a challenge to ecotoxicologists due to the complex geochemical nature of sediments and to the combination of multiple classes of toxicants. Juvenile Senegalese soles were exposed for 14 days in the laboratory and in situ (field) to sediments from three sites (a reference plus two contaminated) of a Portuguese estuary. Sediment characterization confirmed the combination of metals, polycyclic aromatic hydrocarbons and organochlorines in the two contaminated sediments. Changes in liver cytosolic protein regulation patterns were determined by a combination of two-dimensional electrophoresis with de novo sequencing by tandem mass spectrometry. From the forty-one cytosolic proteins found to be deregulated, nineteen were able to be identified, taking part in multiple cellular processes such as anti-oxidative defence, energy production, proteolysis and contaminant catabolism (especially oxidoreductase enzymes). Besides a clear distinction between animals exposed to the reference and contaminated sediments, differences were also observed between laboratory- and in situ-tested fish. Soles exposed in the laboratory to the contaminated sediments failed to induce, or even markedly down-regulated, many proteins, with the exception of a peroxiredoxin (an anti-oxidant enzyme) and a few others, when compared to reference fish. In situ exposure to the contaminated sediments revealed significant up-regulation of basal metabolism-related enzymes, comparatively to the reference condition. Down-regulation of basal metabolism enzymes, related to energy production and gene transcription, in fish exposed in the laboratory to the contaminated sediments, may be linked to sediment-bound contaminants and likely compromised the organisms' ability to deploy adequate responses against insult.

Exploitation of a turbot (Scophthalmus maximus L.) immune-related expressed sequence tag (EST) database for microsatellite screening and validation

In this study, we identified and characterized 160 microsatellite loci from an expressed sequence tag (EST) database generated from immune-related organs of turbot (Scophthalmus maximus). A final set of 83 new polymorphic microsatellites were validated after the analysis of 40 individuals of Atlantic origin including both wild and farmed individuals. The allele number and the expected heterozygosity ranged from 2 to 18 and from 0.021 to 0.951, respectively. Evidences of null alleles at moderate-high frequencies were detected at six loci using population data. None of the analysed loci showed deviations from Mendelian segregation after the analysis of five full-sib families including approximately 92 individuals/family. The markers are used to consolidate the turbot genetic map, and because they are mostly EST-derived, they will be very useful for comparative genomic studies within flatfishes and with model fish species. Using an in silico approach, we detected significant homologies of microsatellite sequences with the EST databases of the flatfish species with highest genomic resources (Senegalese sole, Atlantic halibut, bastard halibut) in 31% of these turbot markers. The conservation of these microsatellites within Pleuronectiformes will pave the way for anchoring genetic maps of different species and identifying genomic regions related to productive traits.

Differential expression patterns of conserved miRNAs and isomiRs during Atlantic halibut development

BACKGROUND

MicroRNAs (miRNAs) play a major role in animal ontogenesis. Size variants of miRNAs, isomiRs, are observed along with the main miRNA types, but their origin and possible biological role are uncovered yet. Developmental profiles of miRNAs have been reported in few fish species only and, to our knowledge, differential expressions of isomiRs have not yet been shown during fish development. Atlantic halibut, Hippoglossus hippoglossus L., undergoes dramatic metamorphosis during early development from symmetrical pelagic larval stage to unsymmetrical flatfish. No data exist on role of miRNAs in halibut metamorphosis.

RESULTS

miRNA profiling using SOLiD deep sequencing technology revealed a total of 199 conserved, one novel antisense, and one miRNA* mature form. Digital expression profiles of selected miRNAs were validated using reverse transcription quantitative PCR. We found developmental transition-specific miRNA expression. Expression of some miRNA* exceeded the guide strand miRNA. We revealed that nucleotide truncations and/or additions at the 3' end of mature miRNAs resulted in size variants showing differential expression patterns during the development in a number of miRNA families. We confirmed the presence of isomiRs by cloning and Sanger sequencing. Also, we found inverse relationship between expression levels of sense/antisense miRNAs during halibut development.

CONCLUSION

Developmental transitions during early development of Atlantic halibut are associated with expression of certain miRNA types. IsomiRs are abundant and often show differential expression during the development.

Transcriptomic analyses in a benthic fish exposed to contaminated estuarine sediments through laboratory and in situ bioassays

The transcription of contaminant response-related genes was investigated in juvenile Senegalese soles exposed to sediments from three distinct sites (a reference plus two contaminated) of a Portuguese estuary (the Sado, W Portugal) through simultaneous 28-day laboratory and in situ bioassays. Transcription of cytochrome P450 1A (CYP1A), metallothionein 1 (MT1), glutathione peroxidase (GPx), catalase (CAT), caspase 3 (CASP3) and 90 kDa heat-shock protein alpha (HSP90AA) was surveyed in the liver by real-time PCR. CASP3 transcription analysis was complemented by surveying apoptosis through the TUNEL reaction. After 14 days of exposure, relative transcription was either reduced or decreased in fish exposed to the contaminated sediments, revealing a disturbance stress phase during which animals failed to respond to insult. After 28 days of exposure all genes' transcription responded to contamination but laboratory and in situ assays depicted distinct patterns of regulation. Although sediments revealed a combination of organic and inorganic toxicants, transcription of the CYP1A gene was consistently correlated to organic contaminants. Metallothionein regulation was found correlated to metallic and organic xenobiotic contamination in the laboratory and in situ, respectively. The transcription of oxidative stress-related genes can be a good indicator of general stress but caution is mandatory when interpreting the results since regulation may be influenced by multiple factors. As for MT1, HSP90 up-regulation has potential to be a good indicator for total contamination, as well as the CASP3 gene, even though hepatocyte apoptosis depicted values inconsistent with sediment contamination, showing that programmed cell death did not directly depend on caspase transcription alone.