Adaptive divergence in a high gene flow environment: Hsc70 variation in the European flounder (Platichthys flesus L.)

Little is known about local adaptations in marine fishes since population genetic surveys in these species have typically not applied genetic markers subject to selection. In this study, we used a candidate gene approach to investigate adaptive population divergence in the European flounder (Platichthys flesus L.) throughout the northeastern Atlantic. We contrasted patterns of genetic variation in a presumably neutral microsatellite baseline to patterns from a heat-shock cognate protein gene, Hsc70. Using two different neutrality tests we found that the microsatellite data set most likely represented a neutral baseline. In contrast, Hsc70 strongly deviated from neutral expectations. Importantly, when estimating standardized levels of population divergence (F(ST)'), we also found a large discrepancy in the patterns of structuring in the two data sets. Thus, samples grouped according to geographical or historical proximity with regards to microsatellites, but according to environmental similarities with regards to Hsc70. The differences between the data sets were particularly pronounced in pairwise comparisons involving populations in the western and central Baltic Sea. For instance, the genetic differentiation between geographically close Baltic Sea and North Sea populations was found to be 0.02 and 0.45 for microsatellites and Hsc70 respectively. Our results strongly suggest adaptive population divergence and indicate local adaptations at the DNA level in a background of high levels of gene flow, typically found in many marine fish species. Furthermore, this study highlights the usefulness of the candidate gene approach for demonstrating local selection in non-model organisms such as most marine fishes.

Evolutionary mechanisms shaping the genetic population structure of marine fishes; lessons from the European flounder (Platichthys flesusL.)

A number of evolutionary mechanisms have been suggested for generating low but significant genetic structuring among marine fish populations. We used nine microsatellite loci and recently developed methods in landscape genetics and coalescence-based estimation of historical gene flow and effective population sizes to assess temporal and spatial dynamics of the population structure in European flounder (Platichthys flesus L.). We collected 1062 flounders from 13 localities in the northeast Atlantic and Baltic Seas and found temporally stable and highly significant genetic differentiation among samples covering a large part of the species' range (global F(ST) = 0.024, P < 0.0001). In addition to historical processes, a number of contemporary acting evolutionary mechanisms were associated with genetic structuring. Physical forces, such as oceanographic and bathymetric barriers, were most likely related with the extreme isolation of the island population at the Faroe Islands. A sharp genetic break was associated with a change in life history from pelagic to benthic spawners in the Baltic Sea. Partial Mantel tests showed that geographical distance per se was not related with genetic structuring among Atlantic and western Baltic Sea samples. Alternative factors, such as dispersal potential and/or environmental gradients, could be important for generating genetic divergence in this region. The results show that the magnitude and scale of structuring generated by a specific mechanism depend critically on its interplay with other evolutionary mechanisms, highlighting the importance of investigating species with wide geographical and ecological distributions to increase our understanding of evolution in the marine environment.

Complete nucleotide sequence of Japanese flounder (Paralichthys olivaceus) mitochondrial genome: structural properties and cue for resolving teleostean relationships

We cloned and sequenced the complete mitochondrial genome of Japanese flounder (Paralichthys olivaceus). A circular 17,090 bp mitochondrial genome from the flounder contains 37 structural genes as in other vertebrates so far reported. This is the first report of the complete mitochondrial sequence from a higher teleostean fish (Acanthopterygii). The organization including gene order is quite similar to that of other teleostean fishes as well as placental mammals. The putative control region of the Japanese flounder mitochondrial genome contains a length variable region of about a 74 bp tandem repeat cluster. As a preliminary study we adopted the maximum likelihood and neighbor-joining inference methods to examine phylogenetic relationships among teleostean and related fishes. Comparisons of amino acid sequences of protein-coding genes and nucleotide sequences of tRNA genes resolved some middle to deep branches among some teleostean fishes. The flounder mitochondrial genome does not show an indication of evolutionary rate difference among teleosts leading to difficulty in phylogenetic analyses, and our data is useful for future evolutionary studies dealing with higher teleostean fishes.

Hepatic transcriptomic profiles of European flounder (Platichthys flesus) from field sites and computational approaches to predict site from stress gene responses following exposure to model toxicants

Genomic technologies offer opportunities to gain a more global assessment of the health status of an organism through an understanding of the functional pathways that are responding to pollutant exposure. We have developed a 13,000 clone cDNA toxicogenomics microarray for Platichthys flesus, the European flounder (EU-GENIPOL Project). We aimed to distinguish the origins of flounder taken from six sampling sites of different pollution status in Northern Europe according to their hepatic gene expression profile using bioinformatic approaches. To determine which gene expression differences may relate to pollutant impact, we have completed complementary laboratory exposures of flounder to selected toxicants and determined the associated gene expression profiles. Using multivariate variable selection coupled with a statistical modelling procedure (GALGO) we can predict geographical site but the accuracy is limited to specific sites. The search space for a combination of genes that effectively predicts class membership is very large, however, by combining the signatures derived from acute laboratory exposure to individual chemicals to limit the search space, a very accurate model for classification of all the different environmental sites was achieved. The final model utilised the expression profiles of 16 clones and validation with a qPCR array comprising these genes correctly assigned the site of origin for fish obtained from three of the sites in an independent sampling. These data would imply that the gene expression fingerprints obtained with these arrays are primarily attributable to variations in chemical pollutant responses at the different sites, indicating their potential utility in environmental impact assessment.

Indirect enzyme-linked immunosorbent assay for the identification of sole (Solea solea), European plaice (Pleuronectes platessa), flounder (Platichthys flesus), and Greenland halibut (Reinhardtius hippoglossoides)

Polyclonal antibodies produced against soluble muscle protein extracts from sole (Solea solea), European plaice (Pleuronectes platessa), flounder (Platichthys flesus), and Greenland halibut (Reinhardtius hippoglossoides) were used in an indirect enzyme-linked immunosorbent assay for the specific identification of fillets from these flatfish species. The assay was performed in two different formats: microtiter plates and immunostick tubes. Immunorecognition of antibodies adsorbed to their specific fish samples was made with goat antirabbit immunoglobulins conjugated to the enzyme horseradish peroxidase. Subsequent enzymatic conversion of the substrate allowed unequivocal identification of all flatfish species studied.

Characterization of novel precursor miRNAs using next generation sequencing and prediction of miRNA targets in Atlantic halibut

BACKGROUND

microRNAs (miRNAs) are implicated in regulation of many cellular processes. miRNAs are processed to their mature functional form in a step-wise manner by multiple proteins and cofactors in the nucleus and cytoplasm. Many miRNAs are conserved across vertebrates. Mature miRNAs have recently been characterized in Atlantic halibut (Hippoglossus hippoglossus L.). The aim of this study was to identify and characterize precursor miRNA (pre-miRNAs) and miRNA targets in this non-model flatfish. Discovery of miRNA precursor forms and targets in non-model organisms is difficult because of limited source information available. Therefore, we have developed a methodology to overcome this limitation.

METHODS

Genomic DNA and small transcriptome of Atlantic halibut were sequenced using Roche 454 pyrosequencing and SOLiD next generation sequencing (NGS), respectively. Identified pre- miRNAs were further validated with reverse-transcription PCR. miRNA targets were identified using miRanda and RNAhybrid target prediction tools using sequences from public databases. Some of miRNA targets were also identified using RACE-PCR. miRNA binding sites were validated with luciferase assay using the RTS34st cell line.

RESULTS

We obtained more than 1.3 M and 92 M sequence reads from 454 genomic DNA sequencing and SOLiD small RNA sequencing, respectively. We identified 34 known and 9 novel pre-miRNAs. We predicted a number of miRNA target genes involved in various biological pathways. miR-24 binding to kisspeptin 1 receptor-2 (kiss1-r2) was confirmed using luciferase assay.

CONCLUSION

This study demonstrates that identification of conserved and novel pre-miRNAs in a non-model vertebrate lacking substantial genomic resources can be performed by combining different next generation sequencing technologies. Our results indicate a wide conservation of miRNA precursors and involvement of miRNA in multiple regulatory pathways, and provide resources for further research on miRNA in non-model animals.

Morphological evidence of correlational selection and ecological segregation between dextral and sinistral forms in a polymorphic flatfish, Platichthys stellatus

Phenotypic polymorphisms in natural systems are often maintained by ecological selection, but only if niche segregation between morphs exists. Polymorphism for eyed-side direction is rare among the approximately 700 species of flatfish (Pleuronectiformes), and the evolutionary mechanisms that maintain it are unknown. Platichthys stellatus (starry flounder) is a polymorphic pleuronectid flatfish exhibiting large, clinal variation in proportion of left-eyed (sinistral) morphs, from 50% in California to 100% in Japan. Here I examined multiple traits related to swimming and foraging performance between sinistral and dextral morphs of P. stellatus from 12 sites to investigate if the two morphs differ in ways that may affect function and ecology. Direction of body asymmetry was correlated with several other characters: on an average, dextral morphs had longer, wider caudal peduncles, shorter snouts and fewer gill rakers than sinistral morphs. Although the differences were small in magnitude, they were consistent in direction across samples, implying that dextral and sinistral starry flounder may be targeting different prey types. Morphological differences between morphs were greatest in samples where the chances of competitive interactions between them were the greatest. These results suggest that the two morphs are not ecologically identical, may represent a rare example of divergent selection maintaining polymorphism of asymmetric forms, and that correlational selection between body asymmetry and other characters may be driven by competitive interactions between sinistral and dextral flatfish. This study is one of very few that demonstrates the ecological significance of direction in a species with polymorphic asymmetric forms.

Inductive expression and characterization analysis of Paralichthys olivaceus pigment epithelium-derived factor in a virally infected cell line

Pigment epithelium-derived factor (PEDF) is acknowledged to be a non-inhibitory member of the serine protease inhibitor (serpin) superfamily, with antiangiogenesis, and neuroprotective and immunoregulatory function, mainly in the tissues of nervous system. Here, A PEDF gene homolog, Paralichthys olivaceus PEDF (PoPEDF), was isolated from flounder embryonic cells (FEC) treated with UV-inactivated Grass carp hemorrhage virus (GCHV) and subsequently identified as a differentially expressed gene. The full length of PoPEDF cDNA is 1803bp with an open reading frame of 1212bp encoding a 403-amino-acid protein. This deduced protein contains an N-terminal signal peptide, a glycosylation site, a consensus serpin motif, and a 34-mer and a 44-mer fragment, all of which are very conserved in the PEDF family. PoPEDF gene exhibits a conserved exon-intron arrangement with 8 exons and 7 introns. This conserved evolutionary relationship was further confirmed by a phylogenetic analysis, where fish PEDFs and mammalian members formed a well-supported clade. Constitutive expression of PoPEDF was widely detected in many tissues. In response to UV-inactivated GCHV or poly(I:C), PEDF mRNA was upregulated in FEC cells with time. This is the first report on the transcriptional induction of PEDF in virally infected cells.

Molecular Cloning, Promoter Analysis and Expression Profiles of the sox3 Gene in Japanese Flounder, Paralichthys olivaceus

Sox3, which belongs to the SoxB1 subgroup, plays major roles in neural and gonadal development. In the present study, Japanese flounder Paralichthys olivaceus sox3 gene (Posox3) and its promoter sequence were isolated and characterized. The deduced PoSox3 protein contained 298 amino acids with a characteristic HMG-box domain. Alignment and phylogenetic analyses indicated that PoSox3 shares highly identical sequence with Sox3 homologues from different species. The promoter region of Posox3 has many potential transcription factor (TF) binding sites. The expression profiles of Posox3 in different developmental stages and diverse adult tissues were analyzed by quantitative real-time RT-PCR (qRT-PCR). Posox3 mRNA was maternally inherited, and maintained at a considerably high expression level between the blastula stage and the hatching stage during embryonic development. Posox3 was abundantly expressed in the adult brain and showed sexually dimorphic expression pattern. In situ hybridization (ISH) was carried out to investigate the cellular distribution of Posox3 in the ovary, and results showed the uniform distribution of Posox3 throughout the cytoplasm of oogonia and stage I–III oocytes. These results indicate that Posox3 has potentially vital roles in embryonic and neural development and may be involved in the oogenesis process. Our work provides a fundamental understanding of the structure and potential functions of Sox3 in Paralichthys olivaceus.

Sequence and organization of the complete mitochondrial genomes of spotted halibut (Verasper variegatus) and barfin flounder (Verasper moseri)

In this work, the mitochondrial genomes for spotted halibut (Verasper variegatus) and barfin flounder (Verasper moseri) were completely sequenced. The entire mitochondrial genome sequences of the spotted halibut and barfin flounder were 17,273 and 17,588 bp in length, respectively. The organization of the two mitochondrial genomes was similar to those reported from other fish mitochondrial genomes containing 37 genes (2 rRNAs, 22 tRNAs and 13 protein-coding genes) and two non-coding regions (control region (CR) and WANCY region). In the CR, the termination associated sequence (ETAS), six central conserved block (CSB-A,B,C,D,E,F), three conserved sequence blocks (CSB1-3) and a region of 61-bp tandem repeat cluster at the end of CSB-3 were identified by similarity comparison with fishes and other vertebrates. The tandem repeat sequences show polymorphism among the different individuals of the two species. The complete mitochondrial genomes of spotted halibut and barfin flounder should be useful for evolutionary studies of flatfishes and other vertebrate species.

Locus number estimation of MHC class II B in stone flounder and Japanese flounder

Members of major histocompatibility complex (MHC) family are important in immune systems. Great efforts have been made to reveal their complicated gene structures. But many existing studies focus on partial sequences of MHC genes. In this study, by gene cloning and sequencing, we identified cDNA sequences and DNA sequences of the MHC class II B in two flatfishes, stone flounder (Kareius bicoloratus) and homozygous diploid Japanese flounder (Paralichthys olivaceus). Eleven cDNA sequences were acquired from eight stone flounder individuals, and most of the polymorphic sites distributed in exons 2 and 3. Twenty-eight alleles were identified from the DNA fragments in these eight individuals. It could be deduced from their Bayesian inference phylogenetic tree that at least four loci of MHC class II B exist in stone flounder. The detailed whole-length DNA sequences in one individual were analyzed, revealing that the intron length varied among different loci. Four different cDNA sequences were identified from one homozygous diploid Japanese flounder individual, implying the existence of at least four loci. Comparison of the cDNA sequences to the DNA sequence confirmed that six exons existed in this gene of Japanese flounder, which was a common feature shared by Pleuronectiformes fishes. Our results proved the multi-locus feature of MHC class II B. The sequences we obtained would provide detailed and systematic data for further research.

An immune responsive complement factor D/adipsin and kallikrein-like serine protease (PoDAK) from the olive flounder Paralichthys olivaceus

The cDNA encoding of a complement factor D/adipsin and kallikrein-like serine protease, designated PoDAK, was isolated from the olive flounder Paralichthys olivaceus. PoDAK cDNA encodes a polypeptide with 277 amino acids containing conserved catalytic triad residues of serine proteases. The amino acid sequence of PoDAK showed high similarity to the kallikrein-like protein of medaka, mammalian adipsin/complement factor D and tissue kallikrein homolog, KT-14 of trout, complement factor D of zebrafish, and shared 31.6-36.8% homology with complement factor D/adipsin known from other species, including mammals. Phylogenetic analysis revealed that PoDAK clustered with the kallikrein-like protein of medaka and mammalian adipsin/complement factor D and tissue kallikrein homolog KT-14 of trout. The expression of PoDAK mRNA was high in the gills and heart, moderate in muscle, liver, intestine, stomach, kidney, and spleen of healthy flounder, and increased in the kidney, liver, and spleen of flounder challenged by the viral hemorrhagic septicemia virus (VHSV) or Streptococcus iniae. In situ hybridization confirmed that PoDAK mRNA is localized in the kidney and heart of individuals infected with VHSV. Further investigations are needed to clarify the function of PoDAK in vivo and in vitro.

Identification and Characterization of a PRDM14 Homolog in Japanese Flounder (Paralichthys olivaceus)

PRDM14 is a PR (PRDI-BF1-RIZ1 homologous) domain protein with six zinc fingers and essential roles in genome-wide epigenetic reprogramming. This protein is required for the establishment of germ cells and the maintenance of the embryonic stem cell ground state. In this study, we cloned the full-length cDNA and genomic DNA of the Paralichthys olivaceus prdm14 (Po-prdm14) gene and isolated the 5' regulatory region of Po-prdm14 by whole-genome sequencing. Peptide sequence alignment, gene structure analysis, and phylogenetic analysis revealed that Po-PRDM14 was homologous to mammalian PRDM14. Results of real-time quantitative polymerase chain reaction amplification (RT-qPCR) and in situ hybridization (ISH) in embryos demonstrated that Po-prdm14 was highly expressed between the morula and late gastrula stages, with its expression peaking in the early gastrula stage. Relatively low expression of Po-prdm14 was observed in the other developmental stages. ISH of gonadal tissues revealed that the transcripts were located in the nucleus of the oocytes in the ovaries but only in the spermatogonia and not the spermatocytes in the testes. We also presume that the Po-prdm14 transcription factor binding sites and their conserved binding region among vertebrates. The combined results suggest that Po-PRDM14 has a conserved function in teleosts and mammals.

PoCRIP1, Paralichthys olivaceus cysteine-rich intestinal protein 1: molecular characterization, expression analysis upon Edwardsiella tarda challenge and a possible role in the immune regulation

Cysteine-rich intestinal protein (CRIP) is a LIM domain protein containing a zinc-finger motif and plays a role in the regulation of the inflammatory immune response. In the present study, we isolated a CRIP1 cDNA, designated PoCRIP1, from an olive flounder Paralichthys olivaceus intestine cDNA library by EST analysis. The PoCRIP cDNA consists of 421 bp with a polyadenylation signal sequence, AATAAA, and a poly(A) tail; it encodes a polypeptide of 76 amino acids containing a double zinc-finger motif (Cys(2)HisCys and Cys(4) sequences). The deduced amino acid sequence of PoCRIP1 showed 75.3-94.7% homology with CRIP1s of other species, including mammals. The PoCRIP1 transcript was highly expressed in the intestine and pyloric ceca and moderately expressed in the gill, heart, kidney, liver, muscle, spleen, skin, and stomach of normal conditioned flounder. Inducible expression of the PoCRIP1 transcript was observed in flounder challenged with Edwardsiella tarda, an economically important pathogen for aquaculture of flounder. Over-expression of PoCRIP1 augmented p65-driven flounder IL-6 promoter activity in HINAE cells. These results suggest that PoCRIP1 may function in the immune response of the flounder through the regulation of cytokine expression.

Growth differentiation factor 15, a novel acute phase response gene in Japanese flounder, Paralichthys olivaceus

The acute phase response, an important aspect of innate immunity, leads to the production of acute phase proteins (APPs) in the liver which would consequently help restore homeostasis to the body. Here, we identified a novel cytokine, growth differentiation factor 15 (GDF15) from Japanese flounder. Three out of the 384 EST sequences derived from liver of Japanese flounder treated with formalin-killed Edwardsiella tarda showed significant homology with GDF of various species. After obtaining the full-length cDNA, the deduced amino acid sequence exhibited low identity (<30%) with GDF15s of higher vertebrates. The predicted ORF of JFGDF15 revealed a signaling peptide at the N terminal, a TGFbeta propeptide domain and a TGFbeta domain. The mature peptide domain of JFGDF15 contains an RXXR motif, a furin cleavage site, required for the release of the mature peptide and conserved amino acids, which are signature features of TGFbeta superfamily proteins. JFGDF15 mRNA transcripts were detected in fish, 6h post-injection with PBS. The transcripts were highly up-regulated in liver at 6h post-injection with formalin-killed E. tarda. Moreover, up-regulation of the transcripts was also observed at 12h post-injection with formalin-killed Streptococcus iniae.

Structure analysis of mtDNA control region of spotted halibut (Verasper variegatus) and its related species

Spotted halibut (Verasper variegatus) is the only species of Genus Verasper in China. The fish was naturally distributed in Yellow Sea and Bohai Sea in northern China and Kyushu in Japan and in Korean sea area. Using PCR product direct sequencing, mitochondrial control region sequences of 24 individuals of spotted halibut was confirmed and analyzed. 4 control region haplotypes, resulting from length heteroplamy of the tandem repeat region, was obtained from these 24 fish. Sequence analysis demonstrated that there were four similar structures in the control region, i.e., extended terminal associated sequences (ETAS), central conserved sequence block (CSB), conserved sequence block (CSB), and repeat region, in V. moseri, Limanda ferruginea, Reinhardtius hippoglossoides, Heppoglossoides platessoides, Paralichthys olivaceous, Solea solea, S. senegalensis, and S. lascari. By comparing with other vertebrates, we found that there were similar repeated sequences immediately after the CSB-3 in all of the anuran species.

Role of Endogenous Cathepsin L in Muscle Protein Degradation in Olive Flounder (Paralichthys olivaceus) Surimi Gel

We investigated the effect of endogenous cathepsin L on surimi gel produced from olive flounder (Paralichthys olivaceus). The amino acid sequences of six proteins predicted or identified as cathepsin L were obtained from the olive flounder genome database, and a phylogenetic analysis was conducted. Next, cathepsin L activity toward N-α-benzyloxycarbonyl-l-phenylalanyl-l-arginine-(7-amino-4-methylcoumarin) (Z-F-R-AMC) was detected in crude olive flounder extract and a crude enzyme preparation. A considerable decrease in the level of myosin heavy chain (MHC) in surimi occurred during autolysis at 60 °C. In contrast, the levels of actin, troponin-T, and tropomyosin decreased only slightly. To prevent protein degradation by cathepsin L, a protease inhibitor was added to surimi. In the presence of 1.0% protease inhibitor, the autolysis of olive flounder surimi at 60 °C was inhibited by 12.2%; the degree of inhibition increased to 44.2% as the inhibitor concentration increased to 3.0%. In addition, the deformation and hardness of modori gel increased as the inhibitor concentration increased to 2.0%. Therefore, cathepsin L plays an important role in protein degradation in surimi, and the quality of surimi gel could be enhanced by inhibiting its activity.

Molecular cloning and expression of CCAAT/enhancer binding proteins in Japanese flounder Paralichthys olivaceus

The findings of this study represent the first report, to the authors' knowledge, of CCAAT/enhancer binding protein (C/EBP) cDNA sequence in a fish species. C/EBP epsilon of Japanese flounder was 1861 bp in length (ORF of 822 bp) encoding for 274 amino acids, with a calculated molecular weight of 30 kDa. Japanese flounder C/EBP beta was found to be 1561 bp in length (ORF of 1041 bp), encoding for 347 amino acids and a calculated molecular weight of 39 kDa. These genes were expressed in various fish organs, tissues and secretions. C/EBP epsilon was detected by Northern blot from total RNA of head and posterior kidney, heart and spleen. However, RT-PCR also detected C/EBP epsilon in brain, spleen and peritoneal cavity fluid and peripheral blood leucocyte cDNA. C/EBP beta was detected by Northern blot analysis in the head and posterior kidney, spleen, intestine, liver, brain, heart, gill and testis and further found by RT-PCR to be detected in mucus, peritoneal cavity fluid, peripheral blood leucocytes and eye cDNA. Phylogenetic analysis placed the Japanese flounder C/EBP beta within the same cluster as previously reported C/EBP beta sequences. However, Japanese flounder C/EBP epsilon sequence data were not found to cluster with the three reported mammalian C/EBP epsilon sequences currently available. Understanding C/EBP transcriptional gene control in commercially important fish species may lead to a better control of disease.

New record of Engyprosopon macrolepis (Regan, 1908) from Taiwan (Pleuronectiformes: Bothidae)

Specimens of the largescale dwarf flounder, Engyprosopon macrolepis (Regan, 1908) (Pleuronectiformes: Bothidae), were collected from the fishery landing of Chi-kan, Peng-hu Islands, western Taiwan. Although this species has been recorded from western Pacific Ocean including the Philippines and Japan, it has never been recorded from Taiwanese waters. Therefore, our specimens represent the first record for Taiwanese waters and fill a gap in the distribution of this species between the two countries. A detailed description of the specimens, including coloration of the fresh condition is provided and compared to the data of specimens collected in other regions. With the addition of E. macrolepis, a total of 11 species of Engyprosopon are recorded from Taiwanese waters, and we suggest that more species are expected to occur off Taiwan.