Large-scale sequence analyses of Atlantic cod

ExTaxsI: an exploration tool of biodiversity molecular data

BACKGROUND

The increasing availability of multi-omics data is leading to regularly revised estimates of existing biodiversity data. In particular, the molecular data enable novel species to be characterized and the information linked to those already observed to be increased with new genomics data. For this reason, the management and visualization of existing molecular data, and their related metadata, through the implementation of easy-to-use IT tools have become a key point to design future research. The more users are able to access biodiversity-related information, the greater the ability of the scientific community to expand its knowledge in this area.

RESULTS

In this article we focus on the development of ExTaxsI (Exploring Taxonomy Information), an IT tool that can retrieve biodiversity data stored in NCBI databases and provide a simple and explorable visualization. We use 3 case studies to show how an efficient organization of the available data can lead to obtaining new information that is fundamental as a starting point for new research. Using this approach highlights the limits in the distribution of data availability, a key factor to consider in the experimental design phase of broad-spectrum studies such as metagenomics.

CONCLUSIONS

ExTaxsI can easily retrieve molecular data and its metadata with an explorable visualization, with the aim of helping researchers to improve experimental designs and highlight the main gaps in the coverage of available data.

Fine-scale population differences in Atlantic cod reproductive success: A potential mechanism for ecological speciation in a marine fish

Successful resource-management and conservation outcomes ideally depend on matching the spatial scales of population demography, local adaptation, and threat mitigation. For marine fish with high dispersal capabilities, this remains a fundamental challenge. Based on daily parentage assignments of more than 4,000 offspring, we document fine-scaled temporal differences in individual reproductive success for two spatially adjacent (<10 km) populations of a broadcast-spawning marine fish. Distinguished by differences in genetics and life history, Atlantic cod (Gadus morhua) from inner- and outer-fjord populations were allowed to compete for mating and reproductive opportunities. After accounting for phenotypic variability in several traits, reproductive success of outer-fjord cod was significantly lower than that of inner-fjord cod. This finding, given that genomically different cod ecotypes inhabit inner- and outer-fjord waters, raises the intriguing hypothesis that the populations might be diverging because of ecological speciation. Individual reproductive success, skewed within both sexes (more so among males), was positively affected by body size, which also influenced the timing of reproduction, larger individuals spawning later among females but earlier among males. Our work suggests that spatial mismatches between management and biological units exist in marine fishes and that studies of reproductive interactions between putative populations or ecotypes can provide an informative basis on which determination of the scale of local adaptation can be ascertained.

A Single Vibrionales 16S rRNA Oligotype Dominates the Intestinal Microbiome in Two Geographically Separated Atlantic cod Populations

Atlantic cod (Gadus morhua) provides an interesting species for the study of host-microbe interactions because it lacks the MHC II complex that is involved in the presentation of extracellular pathogens. Nonetheless, little is known about the diversity of its microbiome in natural populations. Here, we use high-throughput sequencing of the 16S rRNA V4 region, amplified with the primer design of the Earth Microbiome Project (EMP), to investigate the microbial composition in gut content and mucosa of 22 adult individuals from two coastal populations in Norway, located 470 km apart. We identify a core microbiome of 23 OTUs (97% sequence similarity) in all individuals that comprises 93% of the total number of reads. The most abundant orders are classified as Vibrionales, Fusobacteriales, Clostridiales, and Bacteroidales. While mucosal samples show significantly lower diversity than gut content samples, no differences in OTU community composition are observed between the two geographically separated populations. All specimens share a limited number of abundant OTUs. Moreover, the most abundant OTU consists of a single oligotype (order Vibrionales, genus Photobacterium) that represents nearly 50% of the reads in both locations. Our results suggest that these microbiomes comprise a limited number of species or that the EMP V4 primers do not yield sufficient resolution to confidently separate these communities. Our study contributes to a growing body of literature that shows limited spatial differentiation of the intestinal microbiomes in marine fish based on 16S rRNA sequencing, highlighting the need for multi-gene approaches to provide more insight into the diversity of these communities.

Mitochondrial genome variation of Atlantic cod

OBJECTIVE

The objective of this study was to analyse intraspecific sequence variation of Atlantic cod mitochondrial DNA, based on a comprehensive collection of completely sequenced mitochondrial genomes.

RESULTS

We determined the complete mitochondrial DNA sequence of 124 cod specimens from the eastern and western part of the species' distribution range in the North Atlantic Ocean. All specimens harboured a unique mitochondrial DNA haplotype. Nine hundred and fifty-two polymorphic sites were identified, including 109 non-synonymous sites within protein coding regions. Eighteen variable sites were identified as indels, exclusively distributed in structural RNA genes and non-coding regions. Phylogeographic analyses based on 156 available cod mitochondrial genomes did not reveal a clear structure. There was a lack of mitochondrial genetic differentiation between two ecotypes of cod in the eastern North Atlantic, but eastern and western cod were differentiated and mitochondrial genome diversity was higher in the eastern than the western Atlantic, suggesting deviating population histories. The geographic distribution of mitochondrial genome variation seems to be governed by demographic processes and gene flow among ecotypes that are otherwise characterized by localized genomic divergence associated with chromosomal inversions.

Discovery of microRNAs associated with the antiviral immune response of Atlantic cod macrophages

MicroRNAs (miRNAs) are known to play important immunoregulatory roles in teleosts, although miRNAs involved in the antiviral immune response of Atlantic cod (Gadus morhua) were previously uncharacterised. Using deep sequencing and qPCR, the present study was conducted to identify miRNAs responsive to the viral mimic, polyriboinosinic polyribocytidylic acid (pIC) in Atlantic cod macrophages. Macrophage samples isolated from Atlantic cod (n=3) and treated with pIC or phosphate buffered saline (PBS control) for 24 and 72h were used for miRNA profiling. Following deep sequencing, DESeq2 analyses identified four (miR-731-3p, miR-125b-3-3p, miR-150-3p and miR-462-3p) and two (miR-2188-3p and miR-462-3p) significantly differentially expressed miRNAs at 24 and 72h post-stimulation (HPS), respectively. Sequencing-identified miRNAs were subjected to qPCR validation using a larger number of biological replicates (n=6) exposed to pIC or PBS over time (i.e. 12, 24, 48 and 72 HPS). As in sequencing, miR-731-3p, miR-462-3p and miR-2188-3p showed significant up-regulation by pIC. The sequencing results were not qPCR-validated for miR-125b-3-3p and miR-150-3p as up- and down-regulated miRNAs at 24 HPS, respectively; however, qPCR results showed significant up-regulation in response to pIC stimulation at later time points (i.e. 48 and/or 72 HPS). We also used qPCR to assess the expression of other miRNAs that were previously shown as immune responsive in other vertebrates. qPCR results at 48 and/or 72 HPS revealed that miR-128-3-5p, miR-214-1-5p and miR-451-3p were induced by pIC, whereas miR-30b-3p and miR-199-1-3p expression were repressed in response to pIC. The present study identified ten pIC-stimulated miRNAs, suggesting them as important in antiviral immune responses of Atlantic cod macrophages. Some pIC-responsive miRNAs identified in this study were predicted to target putative immune-related genes of Atlantic cod (e.g. miR-30b-3p targeting herc4), although the regulatory functions of these miRNAs need to be validated by future studies.

Study of the plasma proteome of Atlantic cod (Gadus morhua): Effect of exposure to two PAHs and their corresponding diols

Occurrence of polycyclic aromatic hydrocarbon (PAH) contamination in the marine environment represents a risk to marine life and humans. In this study, plasma samples from Atlantic cod (Gadus morhua) were analysed by shotgun mass spectrometry to investigate the plasma proteome in response to exposure to single PAHs (naphthalene or chrysene) and their corresponding metabolites (dihydrodiols). In total, 369 proteins were identified and ranked according to their relative abundance. The levels of 12 proteins were found significantly altered in PAH exposed fish and are proposed as new biomarker candidates. Eleven proteins were upregulated, primarily immunoglobulin components, and one protein was downregulated (antifreeze protein type IV.) The uniformity of the upregulated proteins suggests a triggered immune response in the exposed fish. Overall, the results provide valuable knowledge for future studies of the Atlantic cod plasma proteome and generate grounds for establishing new plasma protein biomarkers for environmental monitoring of PAH related exposure.

Identification of an albumin-like protein in plasma of Atlantic cod (Gadus morhua) and its biomarker potential for PAH contamination

Increased research efforts are currently focusing on Atlantic cod (Gadus morhua) and its significance for monitoring the contaminant situation in marine environments. Polycyclic aromatic hydrocarbons (PAHs) are well known toxic and carcinogenic compounds, thus continuous monitoring is required to ensure ecosystem sustainability and human food safety. A sensitive biomarker of PAH exposure in humans is the detection of PAH metabolites bound to albumin in blood. The potential of a similar PAH-albumin biomarker in Atlantic cod was therefore investigated by a desktop bioinformatic study followed by liquid chromatography mass spectrometry/mass spectrometry analysis of plasma from 16 fish. For the first time, an albumin-like protein in plasma of Atlantic cod is described, and the biomarker potential based on PAH-albumin adduct detection is discussed. Due to the detected low abundance of the albumin-like protein, it was found unlikely to be applicable as a new biomarker tool for evaluation of PAH exposure.

Low plastic ingestion rate in Atlantic cod (Gadus morhua) from Newfoundland destined for human consumption collected through citizen science methods

Marine microplastics are a contaminant of concern because their small size allows ingestion by a wide range of marine life. Using citizen science during the Newfoundland recreational cod fishery, we sampled 205 Atlantic cod (Gadus morhua) destined for human consumption and found that 5 had eaten plastic, an ingestion prevalence rate of 2.4%. This ingestion rate for Atlantic cod is the second lowest recorded rate in the reviewed published literature (the lowest is 1.4%), and the lowest for any fish in the North Atlantic. This is the first report for plastic ingestion in fish in Newfoundland, Canada, a province dependent on fish for sustenance and livelihoods.

The Ontogeny and Brain Distribution Dynamics of the Appetite Regulators NPY, CART and pOX in Larval Atlantic Cod (Gadus morhua L.)

Similar to many marine teleost species, Atlantic cod undergo remarkable physiological changes during the early life stages with concurrent and profound changes in feeding biology and ecology. In contrast to the digestive system, very little is known about the ontogeny and the localization of the centers that control appetite and feed ingestion in the developing brain of fish. We examined the expression patterns of three appetite regulating factors (orexigenic: neuropeptide Y, NPY; prepro-orexin, pOX and anorexigenic: cocaine- and amphetamine-regulated transcript, CART) in discrete brain regions of developing Atlantic cod using chromogenic and double fluorescent in situ hybridization. Differential temporal and spatial expression patterns for each appetite regulator were found from first feeding (4 days post hatch; dph) to juvenile stage (76 dph). Neurons expressing NPY mRNA were detected in the telencephalon (highest expression), diencephalon, and optic tectum from 4 dph onward. CART mRNA expression had a wider distribution along the anterior-posterior brain axis, including both telencephalon and diencephalon from 4 dph. From 46 dph, CART transcripts were also detected in the olfactory bulb, region of the nucleus of medial longitudinal fascicle, optic tectum and midbrain tegmentum. At 4 and 20 dph, pOX mRNA expression was exclusively found in the preoptic region, but extended to the hypothalamus at 46 and 76 dph. Co-expression of both CART and pOX genes were also observed in several hypothalamic neurons throughout larval development. Our results show that both orexigenic and anorexigenic factors are present in the telencephalon, diencephalon and mesencephalon in cod larvae. The telencephalon mostly contains key factors of hunger control (NPY), while the diencephalon, and particularly the hypothalamus may have a more complex role in modulating the multifunctional control of appetite in this species. As the larvae develop, the overall progression in temporal and spatial complexity of NPY, CART and pOX mRNAs expression might be correlated to the maturation of appetite control regulation. These observations suggest that teleost larvae continue to develop the regulatory networks underlying appetite control after onset of exogenous feeding.

Intraspecific Variation in Mitogenomes of Five Crassostrea Species Provides Insight into Oyster Diversification and Speciation

A large number of Crassostrea oysters are found in Asia-Pacific. While analyses of interspecific variation have helped to establish historical relationships among these species, studies on intraspecific variation are necessary to understand their recent evolutionary history and current forces driving population biology. We resequenced 18 and analyzed 31 mitogenomes of five Crassostrea species from China: Crassostrea gigas, Crassostrea angulata, Crassostrea sikamea, Crassostrea ariakensis, and Crassostrea hongkongensis. Our analysis finds abundant insertions, deletions, and single-nucleotide polymorphisms in all species. Intraspecific variation varies greatly among species with polymorphic sites ranging from 54 to 293 and nucleotide diversity ranging from 0.00106 to 0.00683. In all measurements, C. hongkongensis that has the narrowest geographic distribution exhibits the least sequence diversity; C. ariakensis that has the widest distribution shows the highest diversity, and species with intermediate distribution show intermediate levels of diversity. Low sequence diversity in C. hongkongensis may reflect recent bottlenecks that are probably exacerbated by human transplantation. High diversity in C. ariakensis is likely due to divergence of northern and southern China populations that have been separated without gene flow. The significant differences in mitogenome diversity suggest that the five sister species of Crassostrea have experienced different evolutionary forces since their divergence. The recent divergence of two C. ariakensis populations and the C. gigas/angulata species complex provides evidence for continued diversification and speciation of Crassostrea species along China's coast, which are shaped by unknown mechanisms in a north-south divide.

Transcriptome profiling reveals that feeding wild zooplankton to larval Atlantic cod (Gadus morhua) influences suites of genes involved in oxidation-reduction, mitosis, and selenium homeostasis

BACKGROUND

Larval nutrition and growth are key issues for wild and cultured cod. While it was shown previously that larval cod fed wild zooplankton grow faster than those fed only rotifers, the mechanisms involved in this enhanced growth are not completely understood. We used microarrays to identify larval cod transcripts that respond to feeding with small amounts of wild zooplankton (5-10 % of live prey items). The larval transcriptome was compared between 3 treatment groups [fed rotifers (RA), rotifers with protein hydrolysate (RA-PH), or rotifers with zooplankton (RA-Zoo)] at 9-10 mm length [26-30 days post-hatch (dph)] to identify a robust suite of zooplankton-responsive genes (i.e. differentially expressed between RA-Zoo and both other groups).

RESULTS

The microarray experiment identified 147 significantly up-regulated and 156 significantly down-regulated features in RA-Zoo compared with both RA and RA-PH. Gene ontology terms overrepresented in the RA-Zoo responsive gene set included "response to selenium ion" and several related to cell division and oxidation-reduction. Ten selenoprotein-encoding genes, and 2 genes involved in thyroid hormone generation, were up-regulated in RA-Zoo compared with both other groups. Hierarchical clustering of RA-Zoo responsive genes involved in oxidation-reduction and selenium homeostasis demonstrated that only the zooplankton treatment had a considerable and consistent impact on the expression of these genes. Fourteen microarray-identified genes were selected for QPCR involving 9-13 mm larvae, and 13 of these were validated as differentially expressed between RA-Zoo and both other groups at ~9 mm. In contrast, in age-matched (34-35 dph; ~11 mm RA and RA-PH, ~13 mm RA-Zoo) and size-matched (~13 mm) older larvae, only 2 and 3 genes, respectively, showed the same direction of RA-Zoo-responsive change as in ~9 mm larvae.

CONCLUSIONS

The modulation of genes involved in selenium binding, redox homeostasis, and thyroid hormone generation in ~9 mm RA-Zoo larvae in this study may be in response to the relatively high levels of selenium, iodine, and LC-PUFA (potentially causing oxidative stress) in zooplankton. Nonetheless, only a subset of zooplankton-responsive genes in ~9 mm larvae remained so in older larvae, suggesting that the observed transcriptome changes are largely involved in initiating the period of growth enhancement.

Characterization and expression analyses of five interferon regulatory factor transcripts (Irf4a, Irf4b, Irf7, Irf8, Irf10) in Atlantic cod (Gadus morhua)

The interferon regulatory factor (IRF) family of genes encodes a group of transcription factors that have important roles not only in regulating the expression of Type I interferons (IFNs) and other genes in the IFN pathway, but also in growth, development and the regulation of oncogenesis. In this study, several IRF family members (Irf4a, Irf4b, Irf7, Irf8, Irf10) in Atlantic cod (Gadus morhua) were characterized at the cDNA and putative amino acid levels, allowing for phylogenetic analysis of these proteins in teleost fish, as well as the development of gene-specific primers used in RT-PCR and quantitative PCR (QPCR) analyses. Two Atlantic cod Irf10 splice variants were identified and their presence confirmed by sequencing of the Irf10 genomic region. RT-PCR showed that Irf7, Irf8 and both Irf10 transcripts were expressed in all 15 cod tissues tested, while Irf4a and Irf4b were absent in some tissues. QPCR analysis of spleen expression expanded upon this, and upon previous work. All IRF transcripts in the study were responsive to stimulation by the viral mimic poly(I:C), and all except Irf4a were responsive to exposure to formalin-killed Aeromonas salmonicida (ASAL). These IRF genes, thus, are likely important in the cod immune response to both viral and bacterial infections. Increased temperature (10 °C to 16 °C) was also observed to modulate the antibacterial responses of all IRF transcripts, and the antiviral responses of Irf4b and Irf10-v2. This research supports earlier studies which reported that elevated temperature modulates the expression of many immune genes in Atlantic cod.

An evolutionary preserved intergenic spacer in gadiform mitogenomes generates a long noncoding RNA

Background

Vertebrate mitogenomes are economically organized and usually lack intergenic sequences other than the control region. Intergenic spacers located between the tRNA^Thr and tRNA^Pro genes (“T-P spacers”) have been observed in several taxa, including gadiform species, but information about their biological roles and putative functions is still lacking.

Results

Sequence characterization of the complete European hake Merluccius merluccius mitogenome identified a complex T-P spacer ranging in size from 223–532 bp. Further analyses of 32 gadiform species, representing 8 families and 28 genera, revealed the evolutionary preserved presence of T-P spacers across all taxa. Molecular complexity of the T-P spacers was found to be coherent with the phylogenetic relationships, supporting a common ancestral origin and gain of function during codfish evolution. Intraspecific variation of T-P spacer sequences was assessed in 225 Atlantic cod specimens and revealed 26 haplotypes. Pyrosequencing data representing the mito-transcriptome poly (A) fraction in Atlantic cod identified an abundant H-strand specific long noncoding RNA of about 375 nt. The T-P spacer corresponded to the 5’ part of this transcript, which terminated within the control region in a tail-to-tail configuration with the L-strand specific transcript (the 7S RNA).

Conclusions

The T-P spacer is inferred to be evolutionary preserved in gadiform mitogenomes due to gain of function through a long noncoding RNA. We suggest that the T-P spacer adds stability to the H-strand specific long noncoding RNA by forming stable hairpin structures and additional protein binding sites.

Analysing migrations of Atlantic cod Gadus morhua in the north-east Atlantic Ocean: then, now and the future

The application of data storage tags bears the potential for a quantum leap in the research on fish migrations, because not only first-capture and recapture positions are known, but at least theoretically, the migration path during the period at large can be reconstructed. Position, however, cannot be measured directly but has to be estimated using the available data on light, temperature, pressure and salinity. The reconstructed locations based on advanced estimation techniques have been termed geolocations. Examples are discussed which illustrate the applicability of geolocations in individual path descriptions, separation of reproductively isolated populations, timing and areas of spawning, tidal transport and use of protected areas. The examples are based on archival tag data from the North Sea, the Baltic Sea, the Barents Sea and Faroese and Icelandic Waters. Besides presenting the state-of-the-art geolocations for cod Gadus morhua in the north-east Atlantic Ocean, the major aim of this review is to raise awareness of gaps in knowledge and to identify ideas for new research.

Protein genes in repetitive sequence-antifreeze glycoproteins in Atlantic cod genome

Background

Highly repetitive sequences are the bane of genome sequence assembly, and the short read lengths produced by current next generation sequencing technologies further exacerbates this obstacle. An adopted practice is to exclude repetitive sequences in genome data assembly, as the majority of repeats lack protein-coding genes. However, this could result in the exclusion of important genotypes in newly sequenced non-model species. The absence of the antifreeze glycoproteins (AFGP) gene family in the recently sequenced Atlantic cod genome serves as an example.

Results

The Atlantic cod (Gadus morhua) genome was assembled entirely from Roche 454 short reads, demonstrating the feasibility of this approach. However, a well-known major adaptive trait, the AFGP, essential for survival in frigid Arctic marine habitats was absent in the annotated genome. To assess whether this resulted from population difference, we performed Southern blot analysis of genomic DNA from multiple individuals from the North East Arctic cod population that the sequenced cod belonged, and verified that the AFGP genotype is indeed present. We searched the raw assemblies of the Atlantic cod using our G. morhua AFGP gene, and located partial AFGP coding sequences in two sequence scaffolds. We found these two scaffolds constitute a partial genomic AFGP locus through comparative sequence analyses with our newly assembled genomic AFGP locus of the related polar cod, Boreogadus saida. By examining the sequence assembly and annotation methodologies used for the Atlantic cod genome, we deduced the primary cause of the absence of the AFGP gene family from the annotated genome was the removal of all repetitive Roche 454 short reads before sequence assembly, which would exclude most of the highly repetitive AFGP coding sequences. Secondarily, the model teleost genomes used in projection annotation of the Atlantic cod genome have no antifreeze trait, perpetuating the unawareness that the AFGP gene family is missing.

Conclusions

We recovered some of the missing AFGP coding sequences and reconstructed a partial AFGP locus in the Atlantic cod genome, bringing to light that not all repetitive sequences lack protein coding information. Also, reliance on genomes of model organisms as reference for annotating protein-coding gene content of a newly sequenced non-model species could lead to omission of novel genetic traits.

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

Background

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

Results

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

Conclusions

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

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.

Characterization and comparative profiling of MiRNA transcriptomes in bighead carp and silver carp

MicroRNAs (miRNAs) are small non-coding RNA molecules that are processed from large ‘hairpin’ precursors and function as post-transcriptional regulators of target genes. Although many individual miRNAs have recently been extensively studied, there has been very little research on miRNA transcriptomes in teleost fishes. By using high throughput sequencing technology, we have identified 167 and 166 conserved miRNAs (belonging to 108 families) in bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix), respectively. We compared the expression patterns of conserved miRNAs by means of hierarchical clustering analysis and log2 ratio. Results indicated that there is not a strong correlation between sequence conservation and expression conservation, most of these miRNAs have similar expression patterns. However, high expression differences were also identified for several individual miRNAs. Several miRNA* sequences were also found in our dataset and some of them may have regulatory functions. Two computational strategies were used to identify novel miRNAs from un-annotated data in the two carps. A first strategy based on zebrafish genome, identified 8 and 22 novel miRNAs in bighead carp and silver carp, respectively. We postulate that these miRNAs should also exist in the zebrafish, but the methodologies used have not allowed for their detection. In the second strategy we obtained several carp-specific miRNAs, 31 in bighead carp and 32 in silver carp, which showed low expression. Gain and loss of family members were observed in several miRNA families, which suggests that duplication of animal miRNA genes may occur through evolutionary processes which are similar to the protein-coding genes.

The genome sequence of Atlantic cod reveals a unique immune system

Atlantic cod (Gadus morhua) is a large, cold-adapted teleost that sustains long-standing commercial fisheries and incipient aquaculture. Here we present the genome sequence of Atlantic cod, showing evidence for complex thermal adaptations in its haemoglobin gene cluster and an unusual immune architecture compared to other sequenced vertebrates. The genome assembly was obtained exclusively by 454 sequencing of shotgun and paired-end libraries, and automated annotation identified 22,154 genes. The major histocompatibility complex (MHC) II is a conserved feature of the adaptive immune system of jawed vertebrates, but we show that Atlantic cod has lost the genes for MHC II, CD4 and invariant chain (Ii) that are essential for the function of this pathway. Nevertheless, Atlantic cod is not exceptionally susceptible to disease under natural conditions. We find a highly expanded number of MHC I genes and a unique composition of its Toll-like receptor (TLR) families. This indicates how the Atlantic cod immune system has evolved compensatory mechanisms in both adaptive and innate immunity in the absence of MHC II. These observations affect fundamental assumptions about the evolution of the adaptive immune system and its components in vertebrates.

Development and experimental validation of a 20K Atlantic cod (Gadus morhua) oligonucleotide microarray based on a collection of over 150,000 ESTs

The collapse of Atlantic cod (Gadus morhua) wild populations strongly impacted the Atlantic cod fishery and led to the development of cod aquaculture. In order to improve aquaculture and broodstock quality, we need to gain knowledge of genes and pathways involved in Atlantic cod responses to pathogens and other stressors. The Atlantic Cod Genomics and Broodstock Development Project has generated over 150,000 expressed sequence tags from 42 cDNA libraries representing various tissues, developmental stages, and stimuli. We used this resource to develop an Atlantic cod oligonucleotide microarray containing 20,000 unique probes. Selection of sequences from the full range of cDNA libraries enables application of the microarray for a broad spectrum of Atlantic cod functional genomics studies. We included sequences that were highly abundant in suppression subtractive hybridization (SSH) libraries, which were enriched for transcripts responsive to pathogens or other stressors. These sequences represent genes that potentially play an important role in stress and/or immune responses, making the microarray particularly useful for studies of Atlantic cod gene expression responses to immune stimuli and other stressors. To demonstrate its value, we used the microarray to analyze the Atlantic cod spleen response to stimulation with formalin-killed, atypical Aeromonas salmonicida, resulting in a gene expression profile that indicates a strong innate immune response. These results were further validated by quantitative PCR analysis and comparison to results from previous analysis of an SSH library. This study shows that the Atlantic cod 20K oligonucleotide microarray is a valuable new tool for Atlantic cod functional genomics research.

Directed sequencing and annotation of three Dicentrarchus labrax L. chromosomes by applying Sanger- and pyrosequencing technologies on pooled DNA of comparatively mapped BAC clones

Dicentrarchus labrax is one of the major marine aquaculture species in the European Union. In this study, we have developed a directed-sequencing strategy to sequence three sea bass chromosomes and compared results with other teleosts. Three BAC DNA pools were created from sea bass BAC clones that mapped to stickleback chromosomes/groups V, XVII and XXI. The pools were sequenced to 17-39x coverage by pyrosequencing. Data assembly was supported by Sanger reads and mate pair data and resulted in superscaffolds of 13.2 Mb, 17.5 Mb and 13.7 Mb respectively. Annotation features of the superscaffolds include 1477 genes. We analyzed size change of exon, intron and intergenic sequence between teleost species and deduced a simple model for the evolution of genome composition in teleost lineage. Combination of second generation sequencing technologies, Sanger sequencing and genome partitioning strategies allows "high-quality draft assemblies" of chromosome-sized superscaffolds, which are crucial for the prediction and annotation of complete genes.

Characterization of mitochondrial mRNAs in codfish reveals unique features compared to mammals

Expression and processing of mitochondrial gene transcripts are fundamental to mitochondrial function, but information from early vertebrates like teleost fishes is essentially lacking. We have analyzed mitogenome sequences of ten codfishes (family Gadidae), and provide complete sequences from three new species (Saithe, Pollack and Blue whiting). Characterization of the mitochondrial mRNAs in Saithe and Atlantic cod identified a set of ten poly(A) transcripts, and six UAA stop codons are generated by posttranscriptional polyadenylation. Structural assessment of poly(A) sites is consistent with an RNaseP cleavage activity 5' of tRNA acceptor-like stems. COI, ND5 and ND6 mRNAs were found to harbor 3' UTRs with antisense potential extending into neighboring gene regions. While the 3' UTR of COI mRNA is complementary to the tRNA(Ser UCN) and highly similar to that detected in human mitochondria, the ND5 and ND6 3' UTRs appear more heterogenic. Deep sequencing confirms expression of all mitochondrial mRNAs and rRNAs, and provides information about the precise 5' ends in mature transcripts. Our study supports an overall evolutionary conservation in mitochondrial RNA processing events among vertebrates, but reveals some unique 5' and 3' end characteristics in codfish mRNAs with implications to antisense regulation of gene expression.

A novel nucleo-cytoplasmic hybrid clone formed via androgenesis in polyploid gibel carp

BACKGROUND

Unisexual vertebrates have been demonstrated to reproduce by gynogenesis, hybridogenesis, parthenogenesis, or kleptogenesis, however, it is uncertain how the reproduction mode contributes to the clonal diversity. Recently, polyploid gibel carp has been revealed to possess coexisting dual modes of unisexual gynogenesis and sexual reproduction and to have numerous various clones. Using sexual reproduction mating between clone D female and clone A male and subsequent 7 generation multiplying of unisexual gynogenesis, we have created a novel clone strain with more than several hundred millions of individuals. Here, we attempt to identify genetic background of the novel clone and to explore the significant implication for clonal diversity contribution.

METHODS

Several nuclear genome markers and one cytoplasmic marker, the mitochondrial genome sequence, were used to identify the genetic organization of the randomly sampled individuals from different generations of the novel clone.

RESULTS

Chromosome number, Cot-1 repetitive DNA banded karyotype, microsatellite patterns, AFLP profiles and transferrin alleles uniformly indicated that nuclear genome of the novel clone is identical to that of clone A, and significantly different from that of clone D. However, the cytoplasmic marker, its complete mtDNA genome sequence, is same to that of clone D, and different from that of clone A.

CONCLUSIONS

The present data indicate that the novel clone is a nucleo-cytoplasmic hybrid between the known clones A and D, because it originates from the offspring of gonochoristic sexual reproduction mating between clone D female and clone A male, and contains an entire nuclear genome from the paternal clone A and a mtDNA genome (cytoplasm) from the maternal clone D. It is suggested to arise via androgenesis by a mechanism of ploidy doubling of clone A sperm in clone D ooplasm through inhibiting the first mitotic division. Significantly, the selected nucleo-cytoplasmic hybrid female still maintains its gynogenetic ability. Based on the present and previous findings, we discuss the association of rapid genetic changes and high genetic diversity with various ploidy levels and multiple reproduction modes in several unisexual and sexual complexes of vertebrates and even other invertebrates.

Gene expression analysis for the identification of selection and local adaptation in fishes

In recent years, variation in gene expression has been recognized as an important component of environmental adaptation in multiple model species, including a few fish species. There is, however, still little known about the genetic basis of adaptation in gene expression resulting from variation in the aquatic environment (e.g. temperature, salinity and oxygen) and the physiological effect and costs of such differences in gene expression. This review presents and discusses progress and pitfalls of applying gene expression analyses to fishes and suggests simple frameworks to get started with gene expression analysis. It is emphasized that well-planned gene expression studies can serve as an important tool for the identification of selection in local populations of fishes, even for non-traditional model species where limited genomic information is available. Recent studies focusing on gene expression variation among natural fish populations are reviewed, highlighting the latest applications that combine genetic evidence from neutral markers and gene expression data.

Promise and progress in environmental genomics: a status report on the applications of gene expression-based microarray studies in ecologically relevant fish species

The advent of any new technology is typically met with great excitement. So it was a few years ago, when the combination of advances in sequencing technology and the development of microarray technology made measurements of global gene expression in ecologically relevant species possible. Many of the review papers published around that time promised that these new technologies would revolutionize environmental biology as they had revolutionized medicine and related fields. A few years have passed since these technological advancements have been made, and the use of microarray studies in non-model fish species has been adopted in many laboratories internationally. Has the relatively widespread adoption of this technology really revolutionized the fields of environmental biology, including ecotoxicology, aquaculture and ecology, as promised? Or have these studies merely become a novelty and a potential distraction for scientists addressing environmentally relevant questions? In this review, the promises made in early review papers, in particular about the advances that the use of microarrays would enable, are summarized; these claims are compared to the results of recent studies to determine whether the forecasted changes have materialized. Some applications, as discussed in the paper, have been realized and have led to advances in their field, others are still under development.

Genomic organization and gene expression of the multiple globins in Atlantic cod: conservation of globin-flanking genes in chordates infers the origin of the vertebrate globin clusters

Background

The vertebrate globin genes encoding the α- and β-subunits of the tetrameric hemoglobins are clustered at two unlinked loci. The highly conserved linear order of the genes flanking the hemoglobins provides a strong anchor for inferring common ancestry of the globin clusters. In fish, the number of α-β-linked globin genes varies considerably between different sublineages and seems to be related to prevailing physico-chemical conditions. Draft sequences of the Atlantic cod genome enabled us to determine the genomic organization of the globin repertoire in this marine species that copes with fluctuating environments of the temperate and Arctic regions.

Results

The Atlantic cod genome was shown to contain 14 globin genes, including nine hemoglobin genes organized in two unlinked clusters designated β5-α1-β1-α4 and β3-β4-α2-α3-β2. The diverged cod hemoglobin genes displayed different expression levels in adult fish, and tetrameric hemoglobins with or without a Root effect were predicted. The novel finding of maternally inherited hemoglobin mRNAs is consistent with a potential role played by fish hemoglobins in the non-specific immune response. In silico analysis of the six teleost genomes available showed that the two α-β globin clusters are flanked by paralogs of five duplicated genes, in agreement with the proposed teleost-specific duplication of the ancestral vertebrate globin cluster. Screening the genome of extant urochordate and cephalochordate species for conserved globin-flanking genes revealed linkage of RHBDF1, MPG and ARHGAP17 to globin genes in the tunicate Ciona intestinalis, while these genes together with LCMT are closely positioned in amphioxus (Branchiostoma floridae), but seem to be unlinked to the multiple globin genes identified in this species.

Conclusion

The plasticity of Atlantic cod to variable environmental conditions probably involves the expression of multiple globins with potentially different properties. The interspecific difference in number of fish hemoglobin genes contrasts with the highly conserved synteny of the flanking genes. The proximity of globin-flanking genes in the tunicate and amphioxus genomes resembles the RHBDF1-MPG-α-globin-ARHGAP17-LCMT linked genes in man and chicken. We hypothesize that the fusion of the three chordate linkage groups 3, 15 and 17 more than 800 MYA led to the ancestral vertebrate globin cluster during a geological period of increased atmospheric oxygen content.

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.

A brief review on the mechanisms of miRNA regulation

MicroRNAs (miRNAs) are a class of short, endogenously-initiated non-coding RNAs that post-transcriptionally control gene expression via either translational repression or mRNA degradation. It is becoming evident that miRNAs are playing significant roles in regulatory mechanisms operating in various organisms, including developmental timing and host-pathogen interactions as well as cell differentiation, proliferation, apoptosis and tumorigenesis. Likewise, as a regulatory element, miRNA itself is coordinatively modulated by multifarious effectors when carrying out basic functions, such as SNP, miRNA editing, methylation and circadian clock. This mini-review summarized the current understanding of interactions between miRNAs and their targets, including recent advancements in deciphering the regulatory mechanisms that control the biogenesis and functionality of miRNAs in various cellular processes.

RNA deep sequencing of the Atlantic cod transcriptome

The Atlantic cod (Gadus morhua) is an emerging aquaculture species. Efforts to develop and characterize its genomic recourses, including draft-grade genome sequencing, have been initiated by the research community. The transcriptome represents the whole complement of RNA transcripts in cells and tissues and reflects the expressed genes at various life stages, tissue types, physiological states, and environmental conditions. We are investigating the Atlantic cod transcriptome by Roche 454, Illumina GA, and ABI SOLiD deep sequencing platforms and corresponding bioinformatics. Both embryonic developmental stages and adult tissues are studied. Here we summarize our recent progress in the analyses of nuclear and mitochondrial polyA mRNAs, non-protein-coding intermediate RNAs, and regulatory microRNAs.

Toxicological application of primary hepatocyte cell cultures of Atlantic cod (Gadus morhua)--effects of BNF, PCDD and Cd

Fish primary hepatocyte cultures are commonly used for toxicological assessment of contaminants. So far no one has described a protocol on how to use Atlantic cod hepatocytes in bioassays. In this work we describe an experiment in which we were able to isolate intact liver cells from mature individuals. Hepatic cytochrome P450 1A (CYP1A) expression in the isolated cells was evaluated with in situ hybridization after intraperitoneal injection with the strong CYP1A inducer ss-naphthoflavone (BNF). Cod hepatocytes were further exposed to 1,2,3,7,8-polychlorinated dibenzo-p-dioxin (PCDD) and cadmium (Cd). Transcriptional responses of 11 genes were quantified (CYP1A, metallothionein (MT), aryl hydrocarbon receptor 2 (AhR2), UDP-glucuronosyltransferase (UGT), glutathione S-transferase (GST), vitellogenin B (VTGB), hypoxia-inducible factor 1 (HIF1), heme oxygenase 1 (HO-1), transferrin, glutathione peroxidase (GPx) and heat shock protein 70 (HSP70)). Immunohistochemisty evaluation clearly showed elevated CYP1A mRNA expression in primary hepatocytes isolated from BNF-exposed fish. The transcriptional results showed that PCDD exposure resulted in a 311-fold up-regulation of CYP1A and Cd a 1.82-fold increase of MT. Unexpectedly, AhR2 and UGT mRNA levels were not significantly up-regulated in PCDD-exposed cod hepatocytes. HO-1 and transferrin showed a dose-dependent transcriptional response to Cd exposure. Cd appears to act as an endocrine-disrupting metal in exposed primary Atlantic cod hepatocytes. This study demonstrates the use of Atlantic cod primary hepatocyte cultures in toxicological research.

Transcriptome sequencing and development of an expression microarray platform for the domestic ferret

BACKGROUND

The ferret (Mustela putorius furo) represents an attractive animal model for the study of respiratory diseases, including influenza. Despite its importance for biomedical research, the number of reagents for molecular and immunological analysis is restricted. We present here a parallel sequencing effort to produce an extensive EST (expressed sequence tags) dataset derived from a normalized ferret cDNA library made from mRNA from ferret blood, liver, lung, spleen and brain.

RESULTS

We produced more than 500000 sequence reads that were assembled into 16000 partial ferret genes. These genes were combined with the available ferret sequences in the GenBank to develop a ferret specific microarray platform. Using this array, we detected tissue specific expression patterns which were confirmed by quantitative real time PCR assays. We also present a set of 41 ferret genes with even transcription profiles across the tested tissues, indicating their usefulness as housekeeping genes.

CONCLUSION

The tools developed in this study allow for functional genomic analysis and make further development of reagents for the ferret model possible.

Atlantic cod (Gadus morhua) hemoglobin genes: multiplicity and polymorphism

BACKGROUND

Hemoglobin (Hb) polymorphism, assessed by protein gel electrophoresis, has been used almost exclusively to characterize the genetic structure of Atlantic cod (Gadus morhua) populations and to establish correlations with phenotypic traits such as Hb oxygen binding capacity, temperature tolerance and growth characteristics. The genetic system used to explain the results of gel electrophoresis entails the presence of one polymorphic locus with two major alleles (HbI-1; HbI-2). However, vertebrates have more than one gene encoding Hbs and recent studies have reported that more than one Hb gene is present in Atlantic cod. These observations prompted us to re-evaluate the number of Hb genes expressed in Atlantic cod, and to perform an in depth search for polymorphisms that might produce relevant phenotypes for breeding programs.

RESULTS

Analysis of Expressed Sequence Tags (ESTs) led to the identification of nine distinct Hb transcripts; four corresponding to the alpha Hb gene family and five to the beta Hb gene family. To gain insights about the Hb genes encoding these transcripts, genomic sequence data was generated from heterozygous (HbI-1/2) parents and fifteen progeny; five of each HbI type, i.e., HbI-1/1, HbI-1/2 and HbI-2/2. beta Hb genes displayed more polymorphism than alpha Hb genes. Two major allele types (beta1A and beta1B) that differ by two linked non-synonymous substitutions (Met55Val and Lys62Ala) were found in the beta1 Hb gene, and the distribution of these beta1A and beta1B alleles among individuals was congruent with that of the HbI-1 and HbI-2 alleles determined by protein gel electrophoresis. RT-PCR and Q-PCR analysis of the nine Hb genes indicates that all genes are expressed in adult fish, but their level of expression varies greatly; higher expression of almost all Hb genes was found in individuals displaying the HbI-2/2 electrophoretic type.

CONCLUSION

This study indicates that more Hb genes are present and expressed in adult Atlantic cod than previously documented. Our finding that nine Hb genes are expressed simultaneously in adult fish suggests that Atlantic cod, similarly to fish such as rainbow trout, carp, and goldfish, might be able to respond to environmental challenges such as chronic hypoxia or long-term changes in temperature by altering the level of expression of these genes. In this context, the role of the non-conservative substitution Lys62Ala found in the beta1 Hb gene, which appears to explain the occurrence of the HbI-1 and HbI-2 alleles described by gel electrophoresis, and which was found to be present in other fish such as eel, emerald rockcod, rainbow trout and moray, requires further investigation.