Characterization of full-length sequenced cDNA inserts (FLIcs) from Atlantic salmon (Salmo salar)

In silico evaluation of the potential allergenicity of a fungal biomass from Rhizomucor pusillus for use as a novel food ingredient

The world's hunger for novel food ingredients drives the development of safe, sustainable, and nutritious novel food products. For foods containing novel proteins, potential allergenicity of the proteins is a key safety consideration. One such product is a fungal biomass obtained from the fermentation of Rhizomucor pusillus. The annotated whole genome sequence of this strain was subjected to sequence homology searches against the AllergenOnline database (sliding 80-amino acid windows and full sequence searches). In a stepwise manner, proteins were designated as potentially allergenic and were further compared to proteins from commonly consumed foods and from humans. From the sliding 80-mer searches, 356 proteins met the conservative >35% Codex Alimentarius threshold, 72 of which shared ≥50% identity over the full sequence. Although matches were identified between R. pusillus proteins and proteins from allergenic food sources, the matches were limited to minor allergens from these sources, and they shared a greater degree of sequence homology with those from commonly consumed foods and human proteins. Based on the in silico analysis and a literature review for the source organism, the risk of allergenic cross-reactivity of R. pusillus is low.

Molecular characterization of adenosine monophosphate deaminase 1 and the correlation analysis between its mRNA expression levels and inosine monophosphate content in large yellow croaker (Larimichthys crocea)

Declining flesh quality has drawn considerable attention in the farmed large yellow croaker (LYC; Larimichthys crocea) industry. Inosine monophosphate (IMP) is the primary flavor substance in aquatic animals. Adenosine monophosphate deaminase 1 (AMPD1) plays a critical role in IMP formation by catalyzing the deamination of AMP to IMP in the purine nucleotide cycle. To further evaluate the correlation between ampd1 mRNA expression levels and IMP content in the LYC muscle tissue, the relevant open reading frame (ORF) of L. crocea (Lcampd1) was cloned, and the IMP content and Lcampd1 mRNA expression in the muscles of LYCs of different sizes were examined. The ORF cDNA of Lcampd1 was 2211 bp in length and encoded a polypeptide of 736 amino acids (AAs). The deduced protein, LcAMPD1, possesses conserved AMPD active regions (SLSTDDP) and shows high homology with AMPD proteins of other teleost fishes. The genomic DNA sequence of Lcampd1 exhibits a high degree of evolutionary conservation in terms of structural organization among species. Phylogenetic analysis of the deduced AA sequence revealed that teleost fish and mammalian AMPD1 were separate from each other and formed a cluster with AMPD3, suggesting that AMPD1 and AMPD3 arose by duplication of a common primordial gene. In healthy LYC, Lcampd1 mRNA was expressed only in the muscle tissue. The IMP content in the muscle of LYCs with different average body weights was measured by high-performance liquid chromatography; the results showed that the IMP content in the muscle of LYCs with greater body weight was significantly higher than that in LYC with lower body weight. Moreover, a similar trend in Lcampd1 expression was observed in these muscle tissues. The Pearson correlation analysis further showed that the Lcampd1 mRNA expression was positively correlated with IMP content in the muscles of different-sized LYCs. These results suggest the potential function of Lcampd1 in determining the IMP content in LYC and provide a theoretical basis for flesh quality improvement, as well as a scientific basis for the development of the molecular breeding of LYC.

Hsp70/Hsp90 Organising Protein (Hop): Coordinating Much More than Chaperones

The Hsp70/Hsp90 organising protein (Hop, also known as stress-inducible protein 1/STI1/STIP1) has received considerable attention for diverse cellular functions in both healthy and diseased states. There is extensive evidence that intracellular Hop is a co-chaperone of the major chaperones Hsp70 and Hsp90, playing an important role in the productive folding of Hsp90 client proteins, although recent evidence suggests that eukaryotic Hop is regulatory within chaperone complexes rather than essential. Consequently, Hop is implicated in many key signalling pathways, including aberrant pathways leading to cancer. Hop is also secreted, and it is now well established that Hop interacts with the prion protein, PrP^C, to mediate multiple signalling events. The intracellular and extracellular forms of Hop most likely represent two different isoforms, although the molecular determinants of these divergent functions are yet to be identified. There is also a growing body of research that reports the involvement of Hop in cellular activities that appear independent of either chaperones or PrP^C. While the various cellular functions of Hop have been described, its biological function remains elusive. However, recent knockout studies in mammals suggest that Hop has an important role in embryonic development. This review provides a critical overview of the latest molecular, cellular and biological research on Hop, critically evaluating its function in healthy systems and how this function is adapted in diseased states.

MicroSalmon: A Comprehensive, Searchable Resource of Predicted MicroRNA Targets and 3'UTR Cis-Regulatory Elements in the Full-Length Sequenced Atlantic Salmon Transcriptome

Complete 3′UTRs unambiguously assigned to specific mRNA isoforms from the Atlantic salmon full-length (FL) transcriptome were collected into a 3′UTRome. miRNA response elements (MREs) and other cis-regulatory motifs were subsequently predicted and assigned to 3′UTRs of all FL-transcripts. The MicroSalmon GitHub repository provides all results. RNAHybrid and sRNAtoolbox tools predicted the MREs. UTRscan and the Teiresias algorithm predicted other 3′UTR cis-acting motifs, both known vertebrate motifs and putative novel motifs. MicroSalmon provides search programs to retrieve all FL-transcripts targeted by a miRNA (median number 1487), all miRNAs targeting an FL-transcript (median number 27), and other cis-acting motifs. As thousands of FL-transcripts may be targets of each miRNA, additional experimental strategies are necessary to reduce the likely true and relevant targets to a number that may be functionally validated. Low-complexity motifs known to affect mRNA decay in vertebrates were over-represented. Many of these were enriched in the terminal end, while purine- or pyrimidine-rich motifs with unknown functions were enriched immediately downstream of the stop codon. Furthermore, several novel complex motifs were over-represented, indicating conservation and putative function. In conclusion, MicroSalmon is an extensive and useful, searchable resource for study of Atlantic salmon transcript regulation by miRNAs and cis-acting 3′UTR motifs.

A de novo Full-Length mRNA Transcriptome Generated From Hybrid-Corrected PacBio Long-Reads Improves the Transcript Annotation and Identifies Thousands of Novel Splice Variants in Atlantic Salmon

Atlantic salmon (Salmo salar) is a major species produced in world aquaculture and an important vertebrate model organism for studying the process of rediploidization following whole genome duplication events (Ss4R, 80 mya). The current Salmo salar transcriptome is largely generated from genome sequence based in silico predictions supported by ESTs and short-read sequencing data. However, recent progress in long-read sequencing technologies now allows for full-length transcript sequencing from single RNA-molecules. This study provides a de novo full-length mRNA transcriptome from liver, head-kidney and gill materials. A pipeline was developed based on Iso-seq sequencing of long-reads on the PacBio platform (HQ reads) followed by error-correction of the HQ reads by short-reads from the Illumina platform. The pipeline successfully processed more than 1.5 million long-reads and more than 900 million short-reads into error-corrected HQ reads. A surprisingly high percentage (32%) represented expressed interspersed repeats, while the remaining were processed into 71 461 full-length mRNAs from 23 071 loci. Each transcript was supported by several single-molecule long-read sequences and at least three short-reads, assuring a high sequence accuracy. On average, each gene was represented by three isoforms. Comparisons to the current Atlantic salmon transcripts in the RefSeq database showed that the long-read transcriptome validated 25% of all known transcripts, while the remaining full-length transcripts were novel isoforms, but few were transcripts from novel genes. A comparison to the current genome assembly indicates that the long-read transcriptome may aid in improving transcript annotation as well as provide long-read linkage information useful for improving the genome assembly. More than 80% of transcripts were assigned GO terms and thousands of transcripts were from genes or splice-variants expressed in an organ-specific manner demonstrating that hybrid error-corrected long-read transcriptomes may be applied to study genes and splice-variants expressed in certain organs or conditions (e.g., challenge materials). In conclusion, this is the single largest contribution of full-length mRNAs in Atlantic salmon. The results will be of great value to salmon genomics research, and the pipeline outlined may be applied to generate additional de novo transcriptomes in Atlantic Salmon or applied for similar projects in other species.

Effects of environmental enrichment on forebrain neural plasticity and survival success of stocked Atlantic salmon

Fish reared for stocking programmes are severely stimulus deprived compared with their wild conspecifics raised under natural conditions. This leads to reduced behavioural plasticity and low post-release survival of stocked fish. Environmental enrichment can have positive effects on important life skills, such as predator avoidance and foraging behaviour, but the neural mechanisms underpinning these behavioural changes are still largely unknown. In this study, juvenile Atlantic salmon (Salmo salar) were reared in an enriched hatchery environment for 7 weeks, after which neurobiological characteristics and post-release survival were compared with those of fish reared under normal hatchery conditions. Using in situ hybridization and qPCR, we quantified the expression of brain-derived neurotrophic factor (bdnf) and the neural activity marker cfos in telencephalic subregions associated with relational memory, emotional learning and stress reactivity. Aside from lower expression of bdnf in the Dlv (a region associated with relational memory) of enriched salmon, we observed no other significant effects of enrichment in the studied regions. Exposure to an enriched environment increased post-release survival during a 5 month residence in a natural river by 51%. Thus, we demonstrate that environmental enrichment can improve stocking success of Atlantic salmon parr and that environmental enrichment is associated with changes in bdnf expression in the fish's hippocampus-equivalent structure.

miRNAs associated with immune response in teleost fish

MicroRNAs (miRNAs) have been identified as important post transcriptional regulators of gene expression. In higher vertebrates, a subset of miRNAs has been identified as important regulators of a number of key genes in immune system gene networks, and this paper review recent studies on miRNAs associated with immune response in teleost fish. Challenge studies conducted in several species have identified differently expressed miRNAs associated with viral or bacterial infection. The results from these studies point out several miRNAs that are likely to have evolutionary conserved functions that are related to immune response in teleost fish. Changed expression levels of mature miRNAs from the five miRNA genes miRNA-462, miRNA-731, miRNA-146, miRNA-181 and miRNA-223 are observed following viral as well as bacterial infection in several teleost fish. Furthermore, significant changes in expression of mature miRNAs from the five genes miRNA-21, miRNA-155, miRNA-1388, miRNA-99 and miRNA-100 are observed in multiple studies of virus infected fish while changes in expression of mature miRNA from the three genes miRNA-122, miRNA-192 and miRNA-451 are observed in several studies of fish with bacterial infections. Interestingly, some of these genes are not present in higher vertebrates. The function of the evolutionary conserved miRNAs responding to infection depends on the target gene(s) they regulate. A few target genes have been identified while a large number of target genes have been predicted by in silico analysis. The results suggest that many of the targets are genes from the host's immune response gene networks. We propose a model with expected temporal changes in miRNA expression if they target immune response activators/effector genes or immune response inhibitors, respectively. The best way to understand the function of a miRNA is to identify its target gene(s), but as the amount of genome resources for teleost fish is limited, with less well characterized genomes and transcriptomes, identifying the true target genes of the miRNAs associated with the immune response is a challenge. Identifying such target genes by applying new methods and approaches will likely be the next important step to understand the function of the miRNAs associated with immune response in teleost fish.

EGAR, A Food Protein-Derived Tetrapeptide, Reduces Seizure Activity in Pentylenetetrazole-Induced Epilepsy Models Through α-Amino-3-Hydroxy-5-Methyl-4-Isoxazole Propionate Receptors

A primary pathogeny of epilepsy is excessive activation of α-amino-3-hydroxy-5-methyl-4-isoxazole propionate receptors (AMPARs). To find potential molecules to inhibit AMPARs, high-throughput screening was performed in a library of tetrapeptides in silico. Computational results suggest that some tetrapeptides bind stably to the AMPAR. We aligned these sequences of tetrapeptide candidates with those from in vitro digestion of the trout skin protein. Among salmon-derived products, Glu-Gly-Ala-Arg (EGAR) showed a high biological affinity toward AMPAR when tested in silico. Accordingly, natural EGAR was hypothesized to have anticonvulsant activity, and in vitro experiments showed that EGAR selectively inhibited AMPAR-mediated synaptic transmission without affecting the electrophysiological properties of hippocampal pyramidal neurons. In addition, EGAR reduced neuronal spiking in an in vitro seizure model. Moreover, the ability of EGAR to reduce seizures was evaluated in a rodent epilepsy model. Briefer and less severe seizures versus controls were shown after mice were treated with EGAR. In conclusion, the promising experimental results suggest that EGAR inhibitor against AMPARs may be a target for antiepilepsy pharmaceuticals. Epilepsy is a common brain disorder characterized by the occurrence of recurring, unprovoked seizures. Twenty to 30 % of persons with epilepsy do not achieve adequate seizure control with any drug. Here we provide a possibility in which a natural and edible tetrapeptide, EGAR, can act as an antiepileptic agent. We have combined computation with in vitro experiments to show how EGAR modulates epilepsy. We also used an animal model of epilepsy to prove that EGAR can inhibit seizures in vivo. This study suggests EGAR as a potential pharmaceutical for the treatment of epilepsy.

Hsp70/Hsp90 organising protein (hop): beyond interactions with chaperones and prion proteins

The Hsp70/Hsp90 organising protein (Hop), also known as stress-inducible protein 1 (STI1), has received considerable attention for diverse cellular functions in both healthy and diseased states. There is extensive evidence that intracellular Hop is a co-chaperone of the major chaperones Hsp70 and Hsp90, playing an important role in the productive folding of Hsp90 client proteins. Consequently, Hop is implicated in a number of key signalling pathways, including aberrant pathways leading to cancer. However, Hop is also secreted and it is now well established that Hop also serves as a receptor for the prion protein, PrP(C). The intracellular and extracellular forms of Hop most likely represent two different isoforms, although the molecular determinants of these divergent functions are yet to be identified. There is also a growing body of research that reports the involvement of Hop in cellular activities that appear independent of either chaperones or PrP(C). While Hop has been shown to have various cellular functions, its biological function remains elusive. However, recent knockout studies in mammals suggest that Hop has an important role in embryonic development. This review provides a critical overview of the latest molecular, cellular and biological research on Hop, critically evaluating its function in healthy systems and how this function is adapted in diseases states.

The activity of carp muscle-specific creatine kinase at low temperature is enhanced by decreased hydrophobicity of residue 268

Abstract The muscle-specific forms of creatine kinase in rabbit (RM-CK) and carp (M1-CK) exhibit different temperature-dependent functional properties. Replacing the glycine at residue 268 of RM-CK with asparagine increases the enzyme's activity at 10°C. In this study, we investigated how hydrophobicity of residue 268 affects the biochemical properties of RM-CK and M1-CK at low temperature. We generated three mutants of both RM-CK and M1-CK: Asp268, Lys268, and Leu268. The secondary structures of these mutants were similar, as revealed by their circular dichroism spectra. Similar to the Asn268 mutants, the Asp268 and Lys268 mutants of RM-CK and M1-CK exhibited higher specific activities at 10°C and pH 8.0. However, no such effect was observed for the RM-CK and M1-CK Leu268 mutants. While in the presence of cryoprotectant (sucrose or trehalose), the activities of wild-type RM-CK and M1-CK mutant enzymes with a hydrophobic residue at 268 were higher, and the effect was more profound at pH 8.0. It may be inferred that water molecules affect protein conformation around residue 268, thereby influencing protein stability at low temperature.

Differential transcriptome analysis between Paulownia fortunei and its synthesized autopolyploid

Paulownia fortunei is an ecologically and economically important tree species that is widely used as timber and chemical pulp. Its autotetraploid, which carries a number of valuable traits, was successfully induced with colchicine. To identify differences in gene expression between P. fortunei and its synthesized autotetraploid, we performed transcriptome sequencing using an Illumina Genome Analyzer IIx (GAIIx). About 94.8 million reads were generated and assembled into 383,056 transcripts, including 18,984 transcripts with a complete open reading frame. A conducted Basic Local Alignment Search Tool (BLAST) search indicated that 16,004 complete transcripts had significant hits in the National Center for Biotechnology Information (NCBI) non-redundant database. The complete transcripts were given functional assignments using three public protein databases. One thousand one hundred fifty eight differentially expressed complete transcripts were screened through a digital abundance analysis, including transcripts involved in energy metabolism and epigenetic regulation. Finally, the expression levels of several transcripts were confirmed by quantitative real-time PCR. Our results suggested that polyploidization caused epigenetic-related changes, which subsequently resulted in gene expression variation between diploid and autotetraploid P. fortunei. This might be the main mechanism affected by the polyploidization. Our results represent an extensive survey of the P. fortunei transcriptome and will facilitate subsequent functional genomics research in P. fortunei. Moreover, the gene expression profiles of P. fortunei and its autopolyploid will provide a valuable resource for the study of polyploidization.

The Atlantic Bonito (Sarda sarda, Bloch 1793) transcriptome and detection of differential expression during larvae development

The Atlantic bonito (Sarda sarda, Bloch 1793) belongs to the important marine fish species with a wide geographical distribution covering the Atlantic Ocean, the Mediterranean and its bordering seas. Aquaculture practices for this species are still in their infancies and scientific studies are seldom undertaken, mainly because of difficulties in sampling. Thus for small tuna species like the Atlantic bonito only little is known about its biology and regarding the molecular background even less information is available. In the production of marine fish it is known that the most critical period is the larval stages, as high growth rates as well as significant developmental changes take place. In this study we have investigated the transcriptome of the Atlantic bonito of five larvae stages applying Illumina sequencing technology. For non-model species like aquaculture species, transcriptome analysis of RNA samples from individuals using Illumina sequencing technology is technically efficient and cost effective. In the present study a total number of 169,326,711 paired-end reads with a read length of 100 base pairs were generated resulting in a reference transcriptome of 68,220 contigs with an average length of 2054 base pairs. For differential expression analyses single end reads were obtained from different developmental stages and mapped to the constructed reference transcriptome. Differential expression analyses revealed in total 18,657 differentially expressed transcripts and were assigned to five distinguished groups. Each of the five clusters shows stage specific gene expression. We present for the first time in the Atlantic bonito an extensive RNA-Seq based characterization of its transcriptome as well as significant information on differential expression among five developmental larvae stages. The generated transcripts, including SNP and microsatellite information for candidate molecular markers and gene expression information will be a valuable resource for future genetic and molecular studies.

Molecular responses to toxicological stressors: profiling microRNAs in wild Atlantic salmon (Salmo salar) exposed to acidic aluminum-rich water

Atlantic salmon (Salmo salar) is among the most sensitive organisms toward acidic, aluminum exposure. Main documented responses to this type of stress are a combination of hypoxia and loss of blood plasma ions. Physiological responses to stress in fish are often grouped into primary, secondary and tertiary responses, where the above mentioned effects are secondary responses, while primary responses include endocrine changes as measurable levels of catecholamines and corticosteroids. In this study we have exposed young (14 months) Atlantic salmon to acid/Al water (pH ≈ 5.6, Al(i) ≈ 80 μg L⁻¹) for 3 days, and obtained clear and consistent decrease of Na⁺ and Cl⁻ ions, and increases of glucose in blood plasma, hematocrit and P(CO₂) in blood. We did not measure plasma cortisol (primary response compound), but analyzed effects on microRNA level (miRNA) in muscle tissue, as this may represent initial markers of primary stress responses. miRNAs regulate diverse biological processes, many are evolutionarily conserved, and hundreds have been identified in various animals, although only in a few fish species. We used a novel high-throughput sequencing (RNA-Seq) method to identify miRNAs in Atlantic salmon and specific miRNAs as potential early markers for stress. A total of 18 miRNAs were significantly differentially expressed (FDR<0.1) in exposed compared to control fish, four down-regulated and 14 up-regulated. An unsupervised hierarchical clustering of significant miRNAs revealed two clusters representing exposed and non-exposed individuals. Utilizing the genome of the zebrafish and bioinformatic tools, we identified 224 unique miRNAs in the Atlantic salmon samples sequenced. Additional laboratory studies focusing on function, stress dose-responses and temporal expression of the identified miRNAs will facilitate their use as initial markers for stress responses.

Discovery and characterization of miRNA genes in Atlantic salmon (Salmo salar) by use of a deep sequencing approach

Background

MicroRNAs (miRNAs) are an abundant class of endogenous small RNA molecules that downregulate gene expression at the posttranscriptional level. They play important roles in multiple biological processes by regulating genes that control developmental timing, growth, stem cell division and apoptosis by binding to the mRNA of target genes. Despite the position Atlantic salmon (Salmo salar) has as an economically important domesticated animal, there has been little research on miRNAs in this species. Knowledge about miRNAs and their target genes may be used to control health and to improve performance of economically important traits. However, before their biological function can be unravelled they must be identified and annotated. The aims of this study were to identify and characterize miRNA genes in Atlantic salmon by deep sequencing analysis of small RNA libraries from nine different tissues.

Results

A total of 180 distinct mature miRNAs belonging to 106 families of evolutionary conserved miRNAs, and 13 distinct novel mature miRNAs were discovered and characterized. The mature miRNAs corresponded to 521 putative precursor sequences located at unique genome locations. About 40% of these precursors were part of gene clusters, and the majority of the Salmo salar gene clusters discovered were conserved across species. Comparison of expression levels in samples from different tissues applying DESeq indicated that there were tissue specific expression differences in three conserved and one novel miRNA. Ssa-miR 736 was detected in heart tissue only, while two other clustered miRNAs (ssa-miR 212 and132) seems to be at a higher expression level in brain tissue. These observations correlate well with their expected functions as regulators of signal pathways in cardiac and neuronal cells, respectively. Ssa-miR 8163 is one of the novel miRNAs discovered and its function remains unknown. However, differential expression analysis using DESeq suggests that this miRNA is enriched in liver tissue and the precursor was mapped to intron 7 of the transferrin gene.

Conclusions

The identification and annotation of evolutionary conserved and novel Salmo salar miRNAs as well as the characterization of miRNA gene clusters provide biological knowledge that will greatly facilitate further functional studies on miRNAs in this species.

Validation of reference genes for expression studies during craniofacial development in arctic charr

Arctic charr (Salvelinus alpinus) is a highly polymorphic species and in Lake Thingvallavatn, Iceland, four phenotypic morphs have evolved. These differences in morphology, especially in craniofacial structures are already apparent during embryonic development, indicating that genes important in the formation of the craniofacial features are expressed differentially between the morphs. In order to generate tools to examine these expression differences in Arctic charr, the aim of the present study was to identify reference genes for quantitative real-time PCR (qPCR). The specific aim was to select reference genes which are able to detect very small expression differences among different morphs. We selected twelve candidate reference genes from the literature, identified corresponding charr sequences using data derived from transcriptome sequencing (RNA-seq) and examined their expression using qPCR. Many of the candidate reference genes were found to be stably expressed, yet their quality-rank as reference genes varied considerably depending on the type of analysis used. In addition to commonly used software for reference gene validation, we used classical statistics to evaluate expression profiles avoiding a bias for reference genes with similar expression patterns (co-regulation). Based on these analyses we chose three reference genes, ACTB, UB2L3 and IF5A1 for further evaluation. Their consistency was assessed in an expression study of three known craniofacially expressed genes, sparc (or osteonectin), matrix metalloprotease 2 (mmp2) and sox9 (sex-determining region Y box 9 protein) using qPCR in embryo heads derived from four charr groups at three developmental time points. The three reference genes were found to be very suitable for studying expression differences between the morphotypes, enabling robust detection of small relative expression changes during charr development. Further, the results showed that sparc and mmp2 are differentially expressed in embryos of different Arctic charr morphotypes.

Expression of two zona pellucida genes is regulated by 17α-ethinylestradiol in adult rare minnow Gobiocypris rarus

Zona pellucida (ZP) proteins are glycoproteins synthesized in liver, ovary or in both tissues in fish. In the present study, we aimed to determine the responsiveness of ZP2 and ZP3 to 17α-ethinylestradiol (EE2) in adult rare minnow Gobiocypris rarus. The full length of ZP3 cDNA was firstly characterized and its tissue distribution revealed that ZP3 mRNA was predominantly expressed in ovary of G. rarus. The gene expression profiles of ZP2, ZP3 and vitellogenin (VTG) were analyzed in gonad and liver of adult G. rarus exposed to EE2 at 1, 5, 25, and 125 ng/L for 3 and 6 days. The results show that ZP2 is more sensitive than ZP3 in gonads of both genders, and VTG in liver is extremely sensitive to EE2 in male fish. However, at lower concentrations (1 and 5 ng/L), the ZP2 in testes shows higher responsiveness to EE2 compared with VTG in rare minnow. The 5' flanking regions of ZP2 and ZP3 were isolated and the comparison of transcription factors in the regions of ZP2 and ZP3 suggested that the disparity for the responsiveness of ZP2 and ZP3 to EE2 could partly be a result of differential cis-elements such as oocyte-specific protein (Osp1) binding sites or/and sex-determining region Y (SRY) binding site.

Efficient assembly and annotation of the transcriptome of catfish by RNA-Seq analysis of a doubled haploid homozygote

Background

Upon the completion of whole genome sequencing, thorough genome annotation that associates genome sequences with biological meanings is essential. Genome annotation depends on the availability of transcript information as well as orthology information. In teleost fish, genome annotation is seriously hindered by genome duplication. Because of gene duplications, one cannot establish orthologies simply by homology comparisons. Rather intense phylogenetic analysis or structural analysis of orthologies is required for the identification of genes. To conduct phylogenetic analysis and orthology analysis, full-length transcripts are essential. Generation of large numbers of full-length transcripts using traditional transcript sequencing is very difficult and extremely costly.

Results

In this work, we took advantage of a doubled haploid catfish, which has two sets of identical chromosomes and in theory there should be no allelic variations. As such, transcript sequences generated from next-generation sequencing can be favorably assembled into full-length transcripts. Deep sequencing of the doubled haploid channel catfish transcriptome was performed using Illumina HiSeq 2000 platform, yielding over 300 million high-quality trimmed reads totaling 27 Gbp. Assembly of these reads generated 370,798 non-redundant transcript-derived contigs. Functional annotation of the assembly allowed identification of 25,144 unique protein-encoding genes. A total of 2,659 unique genes were identified as putative duplicated genes in the catfish genome because the assembly of the corresponding transcripts harbored PSVs or MSVs (in the form of pseudo-SNPs in the assembly). Of the 25,144 contigs with unique protein hits, around 20,000 contigs matched 50% length of reference proteins, and over 14,000 transcripts were identified as full-length with complete open reading frames. The characterization of consensus sequences surrounding start codon and the stop codon confirmed the correct assembly of the full-length transcripts.

Conclusions

The large set of transcripts assembled in this study is the most comprehensive set of genome resources ever developed from catfish, which will provide the much needed resources for functional genome research in catfish, serving as a reference transcriptome for genome annotation, analysis of gene duplication, gene family structures, and digital gene expression analysis. The putative set of duplicated genes provide a starting point for genome scale analysis of gene duplication in the catfish genome, and should be a valuable resource for comparative genome analysis, genome evolution, and genome function studies.

Sequencing and analysis of full-length cDNAs, 5'-ESTs and 3'-ESTs from a cartilaginous fish, the elephant shark (Callorhinchus milii)

Cartilaginous fishes are the most ancient group of living jawed vertebrates (gnathostomes) and are, therefore, an important reference group for understanding the evolution of vertebrates. The elephant shark (Callorhinchus milii), a holocephalan cartilaginous fish, has been identified as a model cartilaginous fish genome because of its compact genome (∼910 Mb) and a genome project has been initiated to obtain its whole genome sequence. In this study, we have generated and sequenced full-length enriched cDNA libraries of the elephant shark using the 'oligo-capping' method and Sanger sequencing. A total of 6,778 full-length protein-coding cDNA and 10,701 full-length noncoding cDNA were sequenced from six tissues (gills, intestine, kidney, liver, spleen, and testis) of the elephant shark. Analysis of their polyadenylation signals showed that polyadenylation usage in elephant shark is similar to that in mammals. Furthermore, both coding and noncoding transcripts of the elephant shark use the same proportion of canonical polyadenylation sites. Besides BLASTX searches, protein-coding transcripts were annotated by Gene Ontology, InterPro domain, and KEGG pathway analyses. By comparing elephant shark genes to bony vertebrate genes, we identified several ancient genes present in elephant shark but differentially lost in tetrapods or teleosts. Only ∼6% of elephant shark noncoding cDNA showed similarity to known noncoding RNAs (ncRNAs). The rest are either highly divergent ncRNAs or novel ncRNAs. In addition to full-length transcripts, 30,375 5'-ESTs and 41,317 3'-ESTs were sequenced and annotated. The clones and transcripts generated in this study are valuable resources for annotating transcription start sites, exon-intron boundaries, and UTRs of genes in the elephant shark genome, and for the functional characterization of protein sequences. These resources will also be useful for annotating genes in other cartilaginous fishes whose genomes have been targeted for whole genome sequencing.

Molecular characterization of estrogen receptor genes in loach Paramisgurnus dabryanus and their expression upon 17α-ethinylestradiol exposure in juveniles

The full-length cDNAs for estrogen receptor 1 (esr1), esr2a and esr2b were isolated and characterized from the loach (Paramisgurnus dabryanus, Cobitidae, cypriniformes). P. dabryanus Esr1, Esr2a and Esr2b share high amino acids identities with their counterparts of cyprinid species. Quantitative real-time PCR (qRT-PCR) was used to analyze the tissue distribution of esr mRNAs in one-year-old P. dabryanus. The mRNA expression of esr1 in female liver was extremely higher than that in other tissues. esr2a mRNA expression in female intestine and in male muscle was higher than that in other tissues. esr2b mRNA expression was the highest in both male and female intestine. Two-month-old P. dabryanus were exposed to 17α-ethinylestradiol (EE2) for 3weeks and the changes of esr mRNA expression in brain, gonad and liver were analyzed by qRT-PCR. Results showed that EE2 at 1, 5 and 25 ng/L significantly suppressed testicular esr1 mRNA expression in male. The ovarian esr2a mRNA expression was significantly up-regulated at 1 ng/L EE2. In female brain, esr1 mRNA expression was significantly down-regulated at 5 ng/L EE2. Both in males and females, EE2 exposure increased the hepatic esr1 mRNA expression in a concentration-dependent manner. The present study suggests that different esrs in different tissues have differential responsiveness to EE2 and the hepatic esr1 is a sensitive biomarker to EE2 at environmental concentrations in P. dabryanus juveniles. So, the loach P. dabryanus, a typical demersal fish, is a promising ecological model organism to detect estrogenic chemicals in the sediment of aquatic environment by using molecular biomarkers.

Detecting signatures of positive selection in partial sequences generated on a large scale: pitfalls, procedures and resources

Studying the actions of selection provides insight into adaptation, population divergence and gene function. Next-generation sequencing produces large amounts of partial sequences, potentially facilitating efforts to detect signatures of selection based on comparisons between synonymous (d(S)) and nonsynonymous (d(N)) substitutions, and single nucleotide polymorphism assays placed in selected genes would improve the ability to study adaptation in population surveys. However, sequences generated by these technologies are typically short. In nonmodel organisms that are a focus of evolutionary studies, the lack of a reference genome that facilitates the assembly of short sequences has limited surveys of positive selection in large numbers of genes. Here, we describe a series of steps to facilitate these surveys. We provide PERL scripts to assist data analysis, and describe the use of commonly available programs. We demonstrate these approaches in six salmon species, which have partially duplicated genomes. We recommend using multiway blast to optimize the number of alignments between partial coding sequences. Reading frames should be manually detected after alignment with sequences in Genbank using the BLASTX program. We encourage the use of a phylogenetic approach to separate orthologs from paralogs in duplicated genomes. Simple simulations on a gene known to have undergone selection in salmon species, transferrin, showed that the ability to detect selection in short sequences (<600 bp) depended on the proportion of codons under selection (1-2%) within that sequence. This relationship was less relevant in longer sequences. In this exploratory study, we detected 11 genes showing evidence of positive selection.

Targeted SNP discovery in Atlantic salmon (Salmo salar) genes using a 3'UTR-primed SNP detection approach

Background

Single nucleotide polymorphisms (SNPs) represent the most widespread type of DNA variation in vertebrates and may be used as genetic markers for a range of applications. This has led to an increased interest in identification of SNP markers in non-model species and farmed animals. The in silico SNP mining method used for discovery of most known SNPs in Atlantic salmon (Salmo salar) has applied a global (genome-wide) approach. In this study we present a targeted 3'UTR-primed SNP discovery strategy that utilizes sequence data from Salmo salar full length sequenced cDNAs (FLIcs). We compare the efficiency of this new strategy to the in silico SNP mining method when using both methods for targeted SNP discovery.

Results

The SNP discovery efficiency of the two methods was tested in a set of FLIc target genes. The 3'UTR-primed SNP discovery method detected novel SNPs in 35% of the target genes while the in silico SNP mining method detected novel SNPs in 15% of the target genes. Furthermore, the 3'UTR-primed SNP discovery strategy was the less labor intensive one and revealed a higher success rate than the in silico SNP mining method in the initial amplification step. When testing the methods we discovered 112 novel bi-allelic polymorphisms (type I markers) in 88 salmon genes [dbSNP: ss179319972-179320081, ss250608647-250608648], and three of the SNPs discovered were missense substitutions.

Conclusions

Full length insert cDNAs (FLIcs) are important genomic resources that have been developed in many farmed animals. The 3'UTR-primed SNP discovery strategy successfully utilized FLIc data to detect novel SNPs in the partially tetraploid Atlantic salmon. This strategy may therefore be useful for targeted SNP discovery in several species, and particularly useful in species that, like salmonids, have duplicated genomes.

Cathepsin D from the hepatopancreas of the cuttlefish (Sepia officinalis): purification and characterization

Cathepsin D from the hepatopancreas of cuttlefish ( Sepia officinalis ) was purified to homogeneity by precipitation with ammonium sulfate (30-60%, w/v), Sephadex G-100 gel filtration, Mono-S cation-exchange chromatography, Sephadex G-75 gel filtration, and Mono-S FPLC with a 54-fold increase in specific activity and 17% recovery. The molecular weight of the purified cathepsin D was estimated to be 37.5 kDa by sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). On the basis of the native-PAGE and hemoglobin zymography, the purified protease appeared as a single band. The optimum pH and temperature for the cathepsin D activity were pH 3.0 and 50 °C, respectively, using hemoglobin as a substrate. The purified enzyme was completely inhibited by pepstatin A; however, no inhibition was observed with phenylmethylsulfonyl fluoride and ethylenediaminetetraacetic acid. Moreover, the activity was strongly inhibited by SDS and molybdate and enhanced by ATP. The purified cathepsin D was activated by Mg(2+), Ni(2+), Zn(2+), Cu(2+), Cd(2+), Sr(2+), and Co(2+) ions, whereas it was not affected by Na(+), K(+), and Ca(2+) ions. The N-terminal amino acid sequence of the first 13 amino acids of the purified cathepsin D was APTPEPLSNYMDA. S. officinalis cathepsin D, which showed high homology with cathepsin D from marine vertebrates and invertebrates, had a Pro residue at position 6 and a Ser residue at position 8, where Thr and Lys are common in all marine vertebrates cathepsins D. S. officinalis cathepsin D showed high efficiency for the hydrolysis of myofibrillar proteins extracted from cuttlefish muscle.

Sequencing the genome of the Atlantic salmon (Salmo salar)

The International Collaboration to Sequence the Atlantic Salmon Genome (ICSASG) will produce a genome sequence that identifies and physically maps all genes in the Atlantic salmon genome and acts as a reference sequence for other salmonids.

Identification and characterization of full-length cDNAs in channel catfish (Ictalurus punctatus) and blue catfish (Ictalurus furcatus)

BACKGROUND

Genome annotation projects, gene functional studies, and phylogenetic analyses for a given organism all greatly benefit from access to a validated full-length cDNA resource. While increasingly common in model species, full-length cDNA resources in aquaculture species are scarce.

METHODOLOGY AND PRINCIPAL FINDINGS

Through in silico analysis of catfish (Ictalurus spp.) ESTs, a total of 10,037 channel catfish and 7,382 blue catfish cDNA clones were identified as potentially encoding full-length cDNAs. Of this set, a total of 1,169 channel catfish and 933 blue catfish full-length cDNA clones were selected for re-sequencing to provide additional coverage and ensure sequence accuracy. A total of 1,745 unique gene transcripts were identified from the full-length cDNA set, including 1,064 gene transcripts from channel catfish and 681 gene transcripts from blue catfish, with 416 transcripts shared between the two closely related species. Full-length sequence characteristics (ortholog conservation, UTR length, Kozak sequence, and conserved motifs) of the channel and blue catfish were examined in detail. Comparison of gene ontology composition between full-length cDNAs and all catfish ESTs revealed that the full-length cDNA set is representative of the gene diversity encoded in the catfish transcriptome.

CONCLUSIONS

This study describes the first catfish full-length cDNA set constructed from several cDNA libraries. The catfish full-length cDNA sequences, and data gleaned from sequence characteristics analysis, will be a valuable resource for ongoing catfish whole-genome sequencing and future gene-based studies of function and evolution in teleost fishes.

Salmo salar and Esox lucius full-length cDNA sequences reveal changes in evolutionary pressures on a post-tetraploidization genome

Background

Salmonids are one of the most intensely studied fish, in part due to their economic and environmental importance, and in part due to a recent whole genome duplication in the common ancestor of salmonids. This duplication greatly impacts species diversification, functional specialization, and adaptation. Extensive new genomic resources have recently become available for Atlantic salmon (Salmo salar), but documentation of allelic versus duplicate reference genes remains a major uncertainty in the complete characterization of its genome and its evolution.

Results

From existing expressed sequence tag (EST) resources and three new full-length cDNA libraries, 9,057 reference quality full-length gene insert clones were identified for Atlantic salmon. A further 1,365 reference full-length clones were annotated from 29,221 northern pike (Esox lucius) ESTs. Pairwise d_N/d_S comparisons within each of 408 sets of duplicated salmon genes using northern pike as a diploid out-group show asymmetric relaxation of selection on salmon duplicates.

Conclusions

9,057 full-length reference genes were characterized in S. salar and can be used to identify alleles and gene family members. Comparisons of duplicated genes show that while purifying selection is the predominant force acting on both duplicates, consistent with retention of functionality in both copies, some relaxation of pressure on gene duplicates can be identified. In addition, there is evidence that evolution has acted asymmetrically on paralogs, allowing one of the pair to diverge at a faster rate.