Jianbo Yao

Characterization of the rainbow trout transcriptome using Sanger and 454-pyrosequencing approaches

BackgroundRainbow trout are important fish for aquaculture and recreational fisheries and serves as a model species for research investigations associated with carcinogenesis, comparative immunology, toxicology and evolutionary biology. However, to date there is no genome reference sequence to facilitate the development of molecular technologies that utilize high-throughput characterizations of gene expression and genetic variation. Alternatively, transcriptome sequencing is a rapid and efficient means for gene discovery and genetic marker development. Although a large number (258,973) of EST sequences are publicly available, the nature of rainbow trout duplicated genome hinders assembly and complicates annotation.ResultsHigh-throughput deep sequencing of the Swanson rainbow trout doubled-haploid transcriptome using 454-pyrosequencing technology yielded ~1.3 million reads with an average length of 344 bp, a total of 447 million bases. De novo assembly of the sequences yielded 151,847 Tentative Consensus (TC) sequences (average length of 662 bp) and 224,391 singletons. A combination assembly of both the 454-pyrosequencing ESTs and the pre-existing sequences resulted in 161,818 TCs (average length of 758 bp) and 261,071 singletons. Gene Ontology analysis of the combination assembly showed high similarities to transcriptomes of other fish species with known genome sequences.ConclusionThe 454 library significantly increased the suite of ESTs available for rainbow trout, allowing improved assembly and annotation of the transcriptome. Furthermore, the 454 sequencing enables functional genome research in rainbow trout, providing a wealth of sequence data to serve as a reference transcriptome for future studies including identification of paralogous sequences and/or allelic variation, digital gene expression and proteomic research.

Single nucleotide polymorphism discovery in rainbow trout by deep sequencing of a reduced representation library.

BackgroundTo enhance capabilities for genomic analyses in rainbow trout, such as genomic selection, a large suite of polymorphic markers that are amenable to high-throughput genotyping protocols must be identified. Expressed Sequence Tags (ESTs) have been used for single nucleotide polymorphism (SNP) discovery in salmonids. In those strategies, the salmonid semi-tetraploid genomes often led to assemblies of paralogous sequences and therefore resulted in a high rate of false positive SNP identification. Sequencing genomic DNA using primers identified from ESTs proved to be an effective but time consuming methodology of SNP identification in rainbow trout, therefore not suitable for high throughput SNP discovery. In this study, we employed a high-throughput strategy that used pyrosequencing technology to generate data from a reduced representation library constructed with genomic DNA pooled from 96 unrelated rainbow trout that represent the National Center for Cool and Cold Water Aquaculture (NCCCWA) broodstock population.ResultsThe reduced representation library consisted of 440 bp fragments resulting from complete digestion with the restriction enzyme Hae III; sequencing produced 2,000,000 reads providing an average 6 fold coverage of the estimated 150,000 unique genomic restriction fragments (300,000 fragment ends). Three independent data analyses identified 22,022 to 47,128 putative SNPs on 13,140 to 24,627 independent contigs. A set of 384 putative SNPs, randomly selected from the sets produced by the three analyses were genotyped on individual fish to determine the validation rate of putative SNPs among analyses, distinguish apparent SNPs that actually represent paralogous loci in the tetraploid genome, examine Mendelian segregation, and place the validated SNPs on the rainbow trout linkage map. Approximately 48% (183) of the putative SNPs were validated; 167 markers were successfully incorporated into the rainbow trout linkage map. In addition, 2% of the sequences from the validated markers were associated with rainbow trout transcripts.ConclusionThe use of reduced representation libraries and pyrosequencing technology proved to be an effective strategy for the discovery of a high number of putative SNPs in rainbow trout; however, modifications to the technique to decrease the false discovery rate resulting from the evolutionary recent genome duplication would be desirable.

Effect of starvation on global gene expression and proteolysis in rainbow trout (Oncorhynchus mykiss)

BackgroundFast, efficiently growing animals have increased protein synthesis and/or reduced protein degradation relative to slow, inefficiently growing animals. Consequently, minimizing the energetic cost of protein turnover is a strategic goal for enhancing animal growth. Characterization of gene expression profiles associated with protein turnover would allow us to identify genes that could potentially be used as molecular biomarkers to select for germplasm with improved protein accretion.ResultsWe evaluated changes in hepatic global gene expression in response to 3-week starvation in rainbow trout (Oncorhynchus mykiss). Microarray analysis revealed a coordinated, down-regulated expression of protein biosynthesis genes in starved fish. In addition, the expression of genes involved in lipid metabolism/transport, aerobic respiration, blood functions and immune response were decreased in response to starvation. However, the microarray approach did not show a significant increase of gene expression in protein catabolic pathways. Further studies, using real-time PCR and enzyme activity assays, were performed to investigate the expression of genes involved in the major proteolytic pathways including calpains, the multi-catalytic proteasome and cathepsins. Starvation reduced mRNA expression of the calpain inhibitor, calpastatin long isoform (CAST-L), with a subsequent increase in the calpain catalytic activity. In addition, starvation caused a slight but significant increase in 20S proteasome activity without affecting mRNA levels of the proteasome genes. Neither the mRNA levels nor the activities of cathepsin D and L were affected by starvation.ConclusionThese results suggest a significant role of calpain and 20S proteasome pathways in protein mobilization as a source of energy during fasting and a potential association of the CAST-L gene with fish protein accretion.

Identification of Novel Bovine Cumulus Cell Molecular Markers Predictive of Oocyte Competence: Functional and Diagnostic Implications

Abstract The present study was undertaken to discover molecular markers in bovine cumulus cells predictive of oocyte competence and to elucidate their functional significance. Differences in RNA transcript abundance in cumulus cells harvested from oocytes of adult versus prepubertal animals (a model of poor oocyte quality) were identified by microarray analysis. Four genes of interest encoding for the lysosomal cysteine proteinases cathepsins B, S, K, and Z and displaying greater transcript abundance in cumulus cells surrounding oocytes harvested from prepubertal animals were chosen for further investigation. Greater mRNA abundance for such genes in cumulus cells of prepubertal oocytes was confirmed by real-time RT-PCR. Elevated transcript abundance for cathepsins B, S, and Z also was observed in cumulus cells surrounding adult metaphase II oocytes that developed to the blastocyst stage at a low percentage following parthenogenetic activation versus those that developed at a high percentage. Functional significance of cumulus cell cathepsin expression to oocyte competence was confirmed by treatment of cumulus-oocyte complexes during in vitro oocyte maturation with a cell-permeable cysteine proteinase (cathepsin) inhibitor. Inhibitor treatment decreased apoptotic nuclei in the cumulus layer and enhanced development of parthenogenetically activated and in vitro-fertilized adult oocytes to the blastocyst stage. Stimulatory effects of inhibitor treatment during meiotic maturation on subsequent embryonic development were not observed when oocytes were matured in the absence of cumulus cells. The present results support a functional role for cumulus cell cathepsins in compromised oocyte competence and suggest that cumulus cell cathepsin mRNA abundance may be predictive of oocyte quality.

RNA-Seq Identifies SNP Markers for Growth Traits in Rainbow Trout

Fast growth is an important and highly desired trait, which affects the profitability of food animal production, with feed costs accounting for the largest proportion of production costs. Traditional phenotype-based selection is typically used to select for growth traits; however, genetic improvement is slow over generations. Single nucleotide polymorphisms (SNPs) explain 90% of the genetic differences between individuals; therefore, they are most suitable for genetic evaluation and strategies that employ molecular genetics for selective breeding. SNPs found within or near a coding sequence are of particular interest because they are more likely to alter the biological function of a protein. We aimed to use SNPs to identify markers and genes associated with genetic variation in growth. RNA-Seq whole-transcriptome analysis of pooled cDNA samples from a population of rainbow trout selected for improved growth versus unselected genetic cohorts (10 fish from 1 full-sib family each) identified SNP markers associated with growth-rate. The allelic imbalances (the ratio between the allele frequencies of the fast growing sample and that of the slow growing sample) were considered at scores >5.0 as an amplification and <0.2 as loss of heterozygosity. A subset of SNPs (n = 54) were validated and evaluated for association with growth traits in 778 individuals of a three-generation parent/offspring panel representing 40 families. Twenty-two SNP markers and one mitochondrial haplotype were significantly associated with growth traits. Polymorphism of 48 of the markers was confirmed in other commercially important aquaculture stocks. Many markers were clustered into genes of metabolic energy production pathways and are suitable candidates for genetic selection. The study demonstrates that RNA-Seq at low sequence coverage of divergent populations is a fast and effective means of identifying SNPs, with allelic imbalances between phenotypes. This technique is suitable for marker development in non-model species lacking complete and well-annotated genome reference sequences.

Cloning and analysis of fetal ovary microRNAs in cattle.

Ovarian folliculogenesis and early embryogenesis are complex processes, which require tightly regulated expression and interaction of a multitude of genes. Small endogenous RNA molecules, termed microRNAs (miRNAs), are involved in the regulation of gene expression during folliculogenesis and early embryonic development. To identify miRNAs in bovine oocytes/ovaries, a bovine fetal ovary miRNA library was constructed. Sequence analysis of random clones from the library identified 679 miRNA sequences, which represent 58 distinct bovine miRNAs. Of these distinct miRNAs, 42 are known bovine miRNAs present in the miRBase database and the remaining 16 miRNAs include 15 new bovine miRNAs that are homologous to miRNAs identified in other species, and one novel miRNA, which does not match any miRNAs in the database. The precursor sequences for 14 of the new 15 miRNAs as well as the novel miRNA were identified from the bovine genome database and their hairpin structures were predicted. Expression analysis of the 58 miRNAs in fetal ovaries in comparison to somatic tissue pools identified 8 miRNAs predominantly expressed in fetal ovaries. Further analysis of the eight miRNAs in germinal vesicle (GV) stage oocytes identified two miRNAs (bta-mir424 and bta-mir-10b), that are highly abundant in GV oocytes. Both miRNAs show similar expression patterns during oocyte maturation and preimplantation development of bovine embryos, being abundant in GV and MII stage oocytes, as well as in early stage embryos (until 16-cell stage). The amount of the novel miRNA is relatively small in oocytes and early cleavage embryos but greater in blastocysts, suggesting a role of this miRNA in blastocyst cell differentiation.

Cloning and characterization of microRNAs from rainbow trout (Oncorhynchus mykiss): Their expression during early embryonic development

BackgroundCurrent literature and our previous results on expression patterns of oocyte-specific genes and transcription factors suggest a global but highly regulated maternal mRNA degradation at the time of embryonic genome activation (EGA). MicroRNAs (miRNAs) are small, non-coding regulatory RNAs (19–23 nucleotides) that regulate gene expression by guiding target mRNA cleavage or translational inhibition. These regulatory RNAs are potentially involved in the degradation of maternally inherited mRNAs during early embryogenesis.ResultsTo identify miRNAs that might be important for early embryogenesis in rainbow trout, we constructed a miRNA library from a pool of unfertilized eggs and early stage embryos. Sequence analysis of random clones from the library identified 14 miRNAs, 4 of which are novel to rainbow trout. Real-time PCR was used to measure the expression of all cloned miRNAs during embryonic development. Four distinct expression patterns were observed and some miRNAs showed up-regulated expression during EGA. Analysis of tissue distribution of these miRNAs showed that some are present ubiquitously, while others are differentially expressed among different tissues. We also analyzed the expression patterns of Dicer, the enzyme required for the processing of miRNAs and Stat3, a transcription factor involved in activating the transcription of miR-21. Dicer is abundantly expressed during EGA and Stat3 is up-regulated before the onset of EGA.ConclusionThis study led to the discovery of 14 rainbow trout miRNAs. Our data support the notion that Dicer processes miRNAs and Stat3 induces expression of miR-21 and possibly other miRNAs during EGA. These miRNAs in turn guide maternal mRNAs for degradation, which is required for normal embryonic development.

A microRNA repertoire for functional genome research in rainbow trout (Oncorhynchus mykiss).

MicroRNAs (miRNAs) are small, highly conserved, non-coding RNAs that regulate gene expression of target mRNAs through cleavage or translational inhibition. miRNAs are most often identified through computational prediction from genome sequences. The rainbow trout genome sequence is not available yet, which does not allow miRNA prediction for this species which is of great economic interest for aquaculture and sport fisheries, and is a model research organism for studies related to carcinogenesis, toxicology, comparative immunology, disease ecology, physiology and nutrition. To identify miRNAs from rainbow trout, we constructed a miRNA library from a pool of nine somatic tissues. Analysis of the library identified 210 unique sequences representing 54 distinct miRNAs; 50 with conserved sequences matching previously identified miRNAs and four novel miRNAs. In addition, 13 miRNAs were computationally predicted from the rainbow trout transcriptome. Real-time PCR was used to measure miRNA expression patterns in adult somatic tissues and unfertilized eggs. The majority of the miRNAs showed characteristic tissue-specific expression patterns suggesting potential roles in maintaining tissue identity. Potential miRNA-target interactions were computationally predicted and single nucleotide polymorphisms (SNPs) were identified in the miRNAs and their target sites in the rainbow trout transcripts. The rainbow trout miRNAs identified and characterized in this study provide a new tool for functional genome research in salmonids. Tissue-specific miRNAs may serve as molecular markers, predictive of specific functional and diagnostic implications. The data on genetic polymorphisms in miRNA-target interactions is particularly useful for rainbow trout breeding programs.

Role of Importin Alpha8, a New Member of the Importin Alpha Family of Nuclear Transport Proteins, in Early Embryonic Development in Cattle

Abstract Nuclear proteins such as transcription and chromatin remodeling factors are required for initiation of transcription in early embryos before embryonic genome activation. The nuclear transport of these proteins is mediated by transport factors such as importins. Through analysis of expressed sequence tags from a bovine oocyte cDNA library, we identified a new member of the importin alpha family (named importin alpha8). The cloned cDNA for bovine importin alpha8 (KPNA7) is 1817 base pair in length, encoding a protein of 522 amino acids that contains a conserved importin beta-binding domain and seven armadillo motifs. The RT-PCR analysis revealed that KPNA7 mRNA is specifically expressed in ovaries and mature oocytes. Real-time PCR demonstrated that KPNA7 expression in germinal vesicle (GV) oocytes is 33 to 2396 times higher than that of other importin alpha genes and that KPNA7 mRNA is abundant in GV and metaphase II oocytes, as well as in early-stage embryos collected before embryonic genome activation, but is barely detectable in morula- and blastocyst-stage embryos. Similarly, expression of KPNA7 protein is very high in oocytes and early embryos but is low in blastocysts. A glutathione S-transferase pull-down assay revealed that KPNA7 has a strong binding affinity for the nuclear protein nucleoplasmin 2 relative to that of other importin alphas. RNA interference experiments demonstrated that knockdown of KPNA7 in early embryos results in a decreased proportion of embryos developing to the blastocyst stage. These results suggest that KPNA7 may have an important role in the transport of essential nuclear proteins required for early embryogenesis.

Molecular characterization of muscle atrophy and proteolysis associated with spawning in rainbow trout.

Severe muscle deterioration is a physiological response to the energetic demands of fish spawning. This response represents a suitable model to study mechanisms of muscle degradation in fish where typical tetrapod methods, such as muscle unloading, are not applicable. Enzyme activities and mRNA accumulations of genes in major proteolytic pathways, including cathepsins, calpains and the multi-catalytic proteasome, were measured in white muscles of rainbow trout during spawning and post-spawning seasons of gravid fish for comparisons to sterile fish. Fertile fish at spawning had less muscle tissue and less muscle protein compared to sterile fish and post-spawning fertile fish. Muscle deterioration of the fertile fish during spawning was associated with greater mRNA accumulation and elevated activity of cathepsin-L. Concurrently, muscle of spawning fish showed increased mRNA accumulations of cathepsin-D, the calpain regulatory subunit and the proteasome catalytic subunit alpha without corresponding increases in enzyme activities. In addition, elevated activity and increased mRNA accumulation of caspase-9, but not caspase-3, were observed in fertile fish during spawning. This study indicates that cathepsins mediate protein catabolism during spawning in rainbow trout and the catabolic process may involve activation of the apoptosis mediator, caspase-9, but not the apoptosis executioner, caspase-3.