Rainbow smelt (Osmerus mordax) genomic library and EST resources

A Y-linked anti-Müllerian hormone type-II receptor is the sex-determining gene in ayu, Plecoglossus altivelis

Whole-genome duplication and genome compaction are thought to have played important roles in teleost fish evolution. Ayu (or sweetfish), Plecoglossus altivelis, belongs to the superorder Stomiati, order Osmeriformes. Stomiati is phylogenetically classified as sister taxa of Neoteleostei. Thus, ayu holds an important position in the fish tree of life. Although ayu is economically important for the food industry and recreational fishing in Japan, few genomic resources are available for this species. To address this problem, we produced a draft genome sequence of ayu by whole-genome shotgun sequencing and constructed linkage maps using a genotyping-by-sequencing approach. Syntenic analyses of ayu and other teleost fish provided information about chromosomal rearrangements during the divergence of Stomiati, Protacanthopterygii and Neoteleostei. The size of the ayu genome indicates that genome compaction occurred after the divergence of the family Osmeridae. Ayu has an XX/XY sex-determination system for which we identified sex-associated loci by a genome-wide association study by genotyping-by-sequencing and whole-genome resequencing using wild populations. Genome-wide association mapping using wild ayu populations revealed three sex-linked scaffolds (total, 2.03 Mb). Comparison of whole-genome resequencing mapping coverage between males and females identified male-specific regions in sex-linked scaffolds. A duplicate copy of the anti-Müllerian hormone type-II receptor gene (amhr2bY) was found within these male-specific regions, distinct from the autosomal copy of amhr2. Expression of the Y-linked amhr2 gene was male-specific in sox9b-positive somatic cells surrounding germ cells in undifferentiated gonads, whereas autosomal amhr2 transcripts were detected in somatic cells in sexually undifferentiated gonads of both genetic males and females. Loss-of-function mutation for amhr2bY induced male to female sex reversal. Taken together with the known role of Amh and Amhr2 in sex differentiation, these results indicate that the paralog of amhr2 on the ayu Y chromosome determines genetic sex, and the male-specific amh-amhr2 pathway is critical for testicular differentiation in ayu.

Draft Genome of the Korean smelt Hypomesus nipponensis and its transcriptomic responses to heat stress in the liver and muscle

Pond smelt (Hypomesus nipponensis) is a cold-freshwater fish species and a winter economic aquaculture resource in South Korea. Because of its high susceptibility to abnormal water temperature from global warming, a large number of smelt die in hot summers. Here, we present the first draft genome of H. nipponensis and transcriptomic changes in molecular mechanisms or intracellular responses under heat stress. We combined Illumina and PacBio sequencing technologies to generate the draft genome of H. nipponensis. Based on the reference genome, we conducted transcriptome analysis of liver and muscle tissues under normal (NT, 5°C) vs. warm (HT, 23°C) conditions to identify heat stress–induced genes and gene categories. We observed a total of 1987 contigs with N₅₀ of 0.46 Mbp, with the largest contig (3.03 Mbp) in the assembled genome. A total of 20,644 protein-coding genes were predicted, and 19,224 genes were functionally annotated: 15,955 genes for Gene Ontology terms and 11,560 genes for KEGG Orthology. We conducted the lost and gained genes analysis compared with three species that: human, zebrafish, and salmon. In the lost genes analysis, we detected that smelt lost 4461 (22.16%), 2825 (10.62%), and 1499 (3.09%) genes compare with above three species, respectively. In the gained genes analysis, we observed that smelt gained 1133 (5.49%), 1670 (8.09%), and 229 (1.11%) genes compared with the above species, respectively. From transcriptome analysis, a total of 297 and 331 differentially expressed genes (DEGs) with a false discovery rate <0.05 were identified in the liver and muscle tissues, respectively. Gene enrichment analysis of DEGs indicates that upregulated genes were significantly enriched for lipid biosynthetic process (GO:0008610, P < 0.001) and regulation of apoptotic process (GO:0042981, P < 0.01), and genes were downregulated by immune responses such as myeloid cell differentiation (GO:0030099, P < 0.001) in the liver under heat stress. In muscle tissue, upregulated genes were enriched for hypoxia (GO:0001666, P < 0.05), transcription regulator activity (GO:0140110, P < 0.001), and calcium-release channel activity (GO:0015278, P < 0.01), and genes were downregulated for a nicotinamide nucleotide biosynthetic process (GO:0019359, P < 0.01). The results of KEGG pathway analysis were similar to that of gene enrichment analysis. The draft genome and transcriptomic of H. nipponensis will be a useful genetic resource for functional and evolutionary studies. Our findings will improve understanding of molecular mechanisms and heat responses and be useful for predicting survival of the smelt and its closely related species under global warming.

Two potential fish glycerol-3-phosphate phosphatases

Winter-acclimated rainbow smelt (Osmerus mordax Mitchill) produce high levels of glycerol as an antifreeze. A common pathway to glycerol involves the enzyme glycerol-3-phosphate phosphatase (GPP), but no GPP has yet been identified in fish or any other animal. Here, two phosphatases assembled from existing EST libraries (from winter-acclimated smelt and cold-acclimated smelt hepatocytes) were found to resemble a glycerol-associated phosphatase from a glycerol-producing alga, Dunaliella salina, and a recently discovered GPP from a bacterium, Mycobacterium tuberculosis. Recombinant proteins were generated and were found to have GPP activity on the order of a few μMol Pi/mg enzyme/min. The two enzymes have acidic pH optima (~5.5) similar to that previously determined for GPP activity in liver tissue, with about 1/3 of their peak activities at neutral pH. The two enzymes appear to account for the GPP activity of smelt liver, but due to their reduced activities at neutral pH, their contributions to glycerol production in vivo remain unclear. Similar enzymes may be active in a glycerol-producing insect, Dendroctonus ponderosae.

Cloning and characterization of aquaglyceroporin genes from rainbow smelt (Osmerus mordax) and transcript expression in response to cold temperature

Aquaglyceroporins (GLPs) are integral membrane proteins that facilitate passive movement of water, glycerol and urea across cellular membranes. In this study, GLP-encoding genes were characterized in rainbow smelt (Osmerus mordax mordax), an anadromous teleost that accumulates high glycerol and modest urea levels in plasma and tissues as an adaptive cryoprotectant mechanism in sub-zero temperatures. We report the gene and promoter sequences for two aqp10b paralogs (aqp10ba, aqp10bb) that are 82% identical at the predicted amino acid level, and aqp9b. Aqp10bb and aqp9b have the 6 exon structure common to vertebrate GLPs. Aqp10ba has 8 exons; there are two additional exons at the 5' end, and the promoter sequence is different from aqp10bb. Molecular phylogenetic analysis suggests that the aqp10b paralogs arose from a gene duplication event specific to the smelt lineage. Smelt GLP transcripts are ubiquitously expressed; however, aqp10ba transcripts were highest in kidney, aqp10bb transcripts were highest in kidney, intestine, pyloric caeca and brain, and aqp9b transcripts were highest in spleen, liver, red blood cells and kidney. In cold-temperature challenge experiments, plasma glycerol and urea levels were significantly higher in cold- compared to warm-acclimated smelt; however, GLP transcript levels were generally either significantly lower or remained constant. The exception was significantly higher aqp10ba transcript levels in kidney. High aqp10ba transcripts in smelt kidney that increase significantly in response to cold temperature in congruence with plasma urea suggest that this gene duplicate may have evolved to allow the re-absorption of urea to concomitantly conserve nitrogen and prevent freezing.

Origins and functional diversification of salinity-responsive Na(+) , K(+) ATPase α1 paralogs in salmonids

The Salmoniform whole-genome duplication is hypothesized to have facilitated the evolution of anadromy, but little is known about the contribution of paralogs from this event to the physiological performance traits required for anadromy, such as salinity tolerance. Here, we determined when two candidate, salinity-responsive paralogs of the Na(+) , K(+) ATPase α subunit (α1a and α1b) evolved and studied their evolutionary trajectories and tissue-specific expression patterns. We found that these paralogs arose during a small-scale duplication event prior to the Salmoniform, but after the teleost, whole-genome duplication. The 'freshwater paralog' (α1a) is primarily expressed in the gills of Salmoniformes and an unduplicated freshwater sister species (Esox lucius) and experienced positive selection in the freshwater ancestor of Salmoniformes and Esociformes. Contrary to our predictions, the 'saltwater paralog' (α1b), which is more widely expressed than α1a, did not experience positive selection during the evolution of anadromy in the Coregoninae and Salmonine. To determine whether parallel mutations in Na(+) , K(+) ATPase α1 may contribute to salinity tolerance in other fishes, we studied independently evolved salinity-responsive Na(+) , K(+) ATPase α1 paralogs in Anabas testudineus and Oreochromis mossambicus. We found that a quarter of the mutations occurring between salmonid α1a and α1b in functionally important sites also evolved in parallel in at least one of these species. Together, these data argue that paralogs contributing to salinity tolerance evolved prior to the Salmoniform whole-genome duplication and that strong selection and/or functional constraints have led to parallel evolution in salinity-responsive Na(+) , K(+) ATPase α1 paralogs in fishes.

Glycerol synthesis in freeze-resistant rainbow smelt: towards the characterization of a key enzyme glycerol-3-phosphatase

Rainbow smelt (Osmerus mordax) synthesize high amounts of glycerol in winter as a cryoprotectant through the direct dephosphorylation of glycerol-3-phosphate by a phosphatase, glycerol-3-phosphatase (G3Pase). Such a protein is well described in a few species including fungi, bacteria and plants but never studied beyond tissue homogenates in any animal species. Purification, identification and characterization of this enzyme is thus crucial for a better comprehension of the biochemical adaptation in rainbow smelt in response to low temperature and more generally of the biochemical mechanisms involved in glycerol synthesis in animals. This work presents the first attempt to purify G3Pase from smelt liver, the main site of glycerol synthesis for the whole animal. A partial purification was performed, and some characteristics of the protein determined, including optimal pH, K(m) and cation requirements. Smelt G3Pase is most likely a low molecular weight, Mg⁺-dependent and cytosolic phosphatase.

Smelt was the likely beneficiary of an antifreeze gene laterally transferred between fishes

Background

Type II antifreeze protein (AFP) from the rainbow smelt, Osmerus mordax, is a calcium-dependent C-type lectin homolog, similar to the AFPs from herring and sea raven. While C-type lectins are ubiquitous, type II AFPs are only found in a few species in three widely separated branches of teleost fishes. Furthermore, several other non-homologous AFPs are found in intervening species. We have previously postulated that this sporadic distribution has resulted from lateral gene transfer. The alternative hypothesis, that the AFP evolved from a lectin present in a shared ancestor and that this gene was lost in most species, is not favored because both the exon and intron sequences are highly conserved.

Results

Here we have sequenced and annotated a 160 kb smelt BAC clone containing a centrally-located AFP gene along with 14 other genes. Quantitative PCR indicates that there is but a single copy of this gene within the smelt genome, which is atypical for fish AFP genes. The corresponding syntenic region has been identified and searched in a number of other species and found to be devoid of lectin or AFP sequences. Unlike the introns of the AFP gene, the intronic sequences of the flanking genes are not conserved between species. As well, the rate and pattern of mutation in the AFP gene are radically different from those seen in other smelt and herring genes.

Conclusions

These results provide stand-alone support for an example of lateral gene transfer between vertebrate species. They should further inform the debate about genetically modified organisms by showing that gene transfer between ‘higher’ eukaryotes can occur naturally. Analysis of the syntenic regions from several fishes strongly suggests that the smelt acquired the AFP gene from the herring.

Identification and validation of differentially expressed transcripts in a hepatocyte model of cold-induced glycerol production in rainbow smelt (Osmerus mordax)

Rainbow smelt ( Osmerus mordax ) avoid freezing by producing antifreeze protein (AFP) and accumulating glycerol. Glyceroneogenesis occurs in liver via a branch in glycolysis and gluconeogenesis and is activated by low temperature. Hepatocytes were isolated from the livers of fish acclimated to 8°C. Cells were incubated at warm (8°C; nonglycerol accumulating) or cold (0.4°C; glycerol accumulating) temperature over a 72-h time course. Reciprocal suppression subtractive hybridization libraries enriched for cold-responsive transcripts were constructed at 72 h. Microarray analyses using a 16K salmonid cDNA array were performed at 24, 48, and 72 h. Expression of type II AFP and 21 carbohydrate, amino acid, or lipid metabolism-related transcripts were validated using quantitative RT-PCR. Type II AFP transcript levels were not directly temperature related. In cold cells, levels of the glucose synthesis transcript were transiently higher. Increased glycerol production was not associated with increased phosphofructokinase or cytosolic glycerol-3-phosphate dehydrogenase transcript levels. Levels of transcripts (phosphoenolpyruvate carboxykinase, mitochondrial malate dehydrogenase, alanine aminotransferase, glutamate dehydrogenase, and aquaglyceroporin 9) associated with mobilization of amino acids to fuel glycerol accumulation were all transiently higher, suggesting a common regulatory mechanism. In cold compared with warm cells, pyruvate dehydrogenase kinase [an inhibitor of pyruvate dehydrogenase (PDH)] transcript levels were 20-fold higher. Potent inhibition of PDH would direct pyruvate and oxaloacetate derived from amino acids to glycerol, as opposed to oxidation via the citric acid cycle. Levels of a transcript potentially encoding glycerol-3-phosphatase, an enzyme not yet characterized in any vertebrate species, were higher following cold incubation. Finally, this study also presents the novel finding of increased glutamine synthetase transcript levels in response to low temperature.

Seasonal changes in hepatic gene expression reveal modulation of multiple processes in rainbow smelt (Osmerus mordax)

Rainbow smelt (Osmerus mordax) are freeze-resistant fish that accumulate glycerol and produce an antifreeze protein during winter. Quantitative reverse transcription PCR (qPCR) and subtractive hybridization studies have previously revealed five genes in rainbow smelt liver to be differentially regulated in winter in comparison with the fall when water temperatures are warmer. In order to further define the suite of processes that are regulated seasonally, we undertook a large-scale analysis of gene expression by hybridization of smelt cDNA to the salmonid 16K cGRASP microarray. In total, 69 genes were identified as up-regulated and 14 genes as down-regulated under winter conditions. A subset of these genes was examined for differential regulation by qPCR in the individual cDNA samples that were pooled for microarray analysis. Ten of the 15 genes tested showed significant change in the same direction as microarray results, whereas one showed significant change in the opposite direction. Fructose-bisphosphate aldolase B and the cytosolic NAD-dependent glycerol-3-phosphate dehydrogenase were among the most highly up-regulated genes, a result supporting a metabolic focus on glycerol synthesis during winter. Modulation of other processes, including endoplasmic reticulum stress, lipid metabolism and transport, and protein synthesis, was also suggested by the qPCR analysis of array-identified genes. The 15 genes were subsequently examined by qPCR for seasonal variation in expression over five sampling times between October and March, and ten showed significant variation in expression over the sampling period. Taken together, these results provide new understanding of the biochemical adaptations of vertebrates to an extremely low seasonal temperature.

Identification of immune genes of the Agamaki clam (Sinonovacula constricta) by sequencing and bioinformatic analysis of ESTs

The Agamaki clam (Sinonovacula constricta) is an economically important shellfish in Asia. However, genomic research on this species is still in its infancy, and genomic resources are largely unavailable. The objective of this study was to generate expressed sequence tags (ESTs) from a normalized liver complementary DNA library and to identify genes that function in immune defense. A total of 5,296 ESTs were sequenced, from which 540 contigs and 3,473 singletons were identified. BLAST homology analysis indicated that only 20.7% of these ESTs were homologues of known genes while the remaining 79.3% appeared to be novel sequences. Based on sequence similarities, 43 putative immune genes were identified that encode proteases and protease inhibitors, adhesive proteins, stress proteins, lysosomal enzymes, and signal transduction regulators. Our study thus provides both a large collection of novel transcripts and a detailed annotation of immune genes for an important bivalve species.

Generation and analysis of 10 000 ESTs from the half-smooth tongue sole Cynoglossus semilaevis and identification of microsatellite and SNP markers

Three normalized cDNA libraries were constructed, two of which were constructed from reproductive tissues ovary and testis, and the other one from pooled immune tissues including head kidney, intestine, liver and spleen. A total of 10 542 clones were sequenced generating 10 128 expressed sequence tags (ESTs). Cluster analysis indicated a total of 5808 unique sequences including 1712 contigs and 4096 singletons. A total of 4249 (73%) of the unique ESTs had significant hits to the non-redundant protein database, 2253 of which were annotated using Gene Ontology (GO) terms. A total of 311 microsatellites (with 246 having sufficient flanking sequences for primer design) and 6294 putative SNPs were identified. These genome resources provide the material basis for future microarray development, marker validation and genetic linkage and QTL analysis.

Oogenesis in teleosts: how eggs are formed

One of the major objectives of the aquaculture industry is the production of a large number of viable eggs with high survival. Major achievements have been made in recent years in improving protocols for higher efficiency of egg production and viability of progeny. Main gaps remain, however, in understanding the dynamic processes associated with oogenesis, the formation of an egg, from the time that germ cells turn into oogonia, until the release of ova during spawning in teleosts. Recent studies on primordial germ-cells, yolk protein precursors and their processing within the developing oocyte, the deposition of vitamins in eggs, structure and function of egg envelopes and oocyte maturation processes, further reveal the complexity of oogenesis. Moreover, numerous circulating endocrine and locally-acting paracrine and autocrine factors regulate the various stages of oocyte development and maturation. Though it is clear that the major regulators during vitellogenesis and oocyte maturation are the pituitary gonadotropins (LH and FSH) and sex steroids, the picture emerging from recent studies is of complex hormonal cross-talk at all stages between the developing oocyte and its surrounding follicle layers to ensure coordination of the various processes that are involved in the production of a fertilizable egg. In this review we aim at highlighting recent advances on teleost fish oocyte differentiation, maturation and ovulation, including those involved in the degeneration and reabsorption of ovarian follicles (atresia). The role of blood-borne and local ovarian factors in the regulation of the key steps of development reveal new aspects associated with egg formation.

Stimulation of growth and changes in the hepatic transcriptome by 17β-estradiol in the yellow perch (Perca flavescens)

The effects of dietary 17β-estradiol (E2) on growth and liver transcriptomics were investigated in the yellow perch ( Perca flavescens). After a 3-mo treatment, E2 significantly stimulated an increase in length and weight of juvenile male and female perch relative to control animals. The increase was significantly greater in females compared with males. Separate, unnormalized cDNA libraries were constructed from equal quantities of RNA from 6 male and 6 female livers of E2-treated and control perch, and 3,546 and 3,719 expressed sequence tags (ESTs) were obtained, respectively. To characterize E2-regulated transcripts, EST frequencies between libraries were calculated within contiguous sequences that were assembled from the combined ESTs of both libraries. Frequencies were also determined in EST transcript groupings produced by aligning all of the ESTs from both libraries at the nucleotide level. From these analyses, there were 28 annotated transcripts that were regulated by 75% between libraries and for which there were at least 5 ESTs of the same transcript between libraries. Regulation of a subset ( 14 ) of these transcripts was confirmed by quantitative reverse transcription-polymerase chain reaction (QPCR). Transcripts that were upregulated by E2 included reproduction-related proteins, binding proteins, and proteases and protease inhibitors. While not part of the transcript frequency analysis, QPCR showed significant upregulation of estrogen receptor esr1 and of insulin-like growth factor I (IGF-I) in E2 livers. E2-downregulated transcripts represented a variety of functional categories including components of the respiratory chain, lipid transport and metabolism, glycolysis, amino acid and nitrogen metabolism, binding proteins, a hydrolytic enzyme, and a transcriptional regulator. In perch it appears that exogenous estrogen drastically shifts liver metabolism toward the production of lipoproteins and carbohydrate binding proteins, and that the growth-promoting action may involve an increase in hepatic IGF-I production.

Functional annotation and analysis of expressed sequence tags from the hepatopancreas of mitten crab (Eriocheir sinensis)

Due to its popularity as a traditional food, intensive harvesting of the mitten crab (Eriocheir sinensis) is common and has lead to an increase in disease incidence, resulting in catastrophic losses to crab aquaculture. The hepatopancreas of E. sinensis is not only an important digestive organ but also an indispensable immune organ. We constructed a nonnormalized cDNA library from the hepatopancreas of E. sinensis and acquired 3,297 high-quality expressed sequence tags representing 1,178 unigenes. More than half of these unigenes were novel genes for this species; the remaining had homologs in public databases, which is of great importance for future functional research. We also investigated the association of these genes with immune processes for insight into one of the main functions of the hepatopancreas besides metabolism. Despite the relatively low sampling scalar of our cDNA library, we were able to demonstrate several important properties of the hepatopancreatic transcriptome and identified numerous genes that were closely associated with immune responses. These results might serve as the basis for an in-depth genomics study of E. sinensis, including transcriptome analysis, physical mapping, and whole genome sequencing.

A salmonid EST genomic study: genes, duplications, phylogeny and microarrays

Background

Salmonids are of interest because of their relatively recent genome duplication, and their extensive use in wild fisheries and aquaculture. A comprehensive gene list and a comparison of genes in some of the different species provide valuable genomic information for one of the most widely studied groups of fish.

Results

298,304 expressed sequence tags (ESTs) from Atlantic salmon (69% of the total), 11,664 chinook, 10,813 sockeye, 10,051 brook trout, 10,975 grayling, 8,630 lake whitefish, and 3,624 northern pike ESTs were obtained in this study and have been deposited into the public databases. Contigs were built and putative full-length Atlantic salmon clones have been identified. A database containing ESTs, assemblies, consensus sequences, open reading frames, gene predictions and putative annotation is available. The overall similarity between Atlantic salmon ESTs and those of rainbow trout, chinook, sockeye, brook trout, grayling, lake whitefish, northern pike and rainbow smelt is 93.4, 94.2, 94.6, 94.4, 92.5, 91.7, 89.6, and 86.2% respectively. An analysis of 78 transcript sets show Salmo as a sister group to Oncorhynchus and Salvelinus within Salmoninae, and Thymallinae as a sister group to Salmoninae and Coregoninae within Salmonidae. Extensive gene duplication is consistent with a genome duplication in the common ancestor of salmonids. Using all of the available EST data, a new expanded salmonid cDNA microarray of 32,000 features was created. Cross-species hybridizations to this cDNA microarray indicate that this resource will be useful for studies of all 68 salmonid species.

Conclusion

An extensive collection and analysis of salmonid RNA putative transcripts indicate that Pacific salmon, Atlantic salmon and charr are 94–96% similar while the more distant whitefish, grayling, pike and smelt are 93, 92, 89 and 86% similar to salmon. The salmonid transcriptome reveals a complex history of gene duplication that is consistent with an ancestral salmonid genome duplication hypothesis. Genome resources, including a new 32 K microarray, provide valuable new tools to study salmonids.

Genomic resources for a commercial flatfish, the Senegalese sole (Solea senegalensis): EST sequencing, oligo microarray design, and development of the Soleamold bioinformatic platform

Background

The Senegalese sole, Solea senegalensis, is a highly prized flatfish of growing commercial interest for aquaculture in Southern Europe. However, despite the industrial production of Senegalese sole being hampered primarily by lack of information on the physiological mechanisms involved in reproduction, growth and immunity, very limited genomic information is available on this species.

Results

Sequencing of a S. senegalensis multi-tissue normalized cDNA library, from adult tissues (brain, stomach, intestine, liver, ovary, and testis), larval stages (pre-metamorphosis, metamorphosis), juvenile stages (post-metamorphosis, abnormal fish), and undifferentiated gonads, generated 10,185 expressed sequence tags (ESTs). Clones were sequenced from the 3'-end to identify isoform specific sequences. Assembly of the entire EST collection into contigs gave 5,208 unique sequences of which 1,769 (34%) had matches in GenBank, thus showing a low level of redundancy. The sequence of the 5,208 unigenes was used to design and validate an oligonucleotide microarray representing 5,087 unique Senegalese sole transcripts. Finally, a novel interactive bioinformatic platform, Soleamold, was developed for the Senegalese sole EST collection as well as microarray and ISH data.

Conclusion

New genomic resources have been developed for S. senegalensis, an economically important fish in aquaculture, which include a collection of expressed genes, an oligonucleotide microarray, and a publicly available bioinformatic platform that can be used to study gene expression in this species. These resources will help elucidate transcriptional regulation in wild and captive Senegalese sole for optimization of its production under intensive culture conditions.