- Research article
- Open Access
Transcriptome analysis of egg viability in rainbow trout, Oncorhynchus mykiss
BMC Genomicsvolume 20, Article number: 319 (2019)
Maternal transcripts are accumulated in the oocyte during oogenesis to provide for protein synthesis from oocyte maturation through early embryonic development, when nuclear transcription is silenced. The maternal mRNAs have short poly(A) tails after undergoing post-transcriptional processing necessary for stabilizing them for storage. The transcripts undergo cytoplasmic polyadenylation when they are to be translated. Transcriptome analyses comparing total mRNA and elongated poly(A) mRNA content among eggs of different quality can provide insight into molecular mechanisms affecting egg developmental competence in rainbow trout. The present study used RNA-seq to compare transcriptomes of unfertilized eggs of rainbow trout females yielding different eyeing rates, following rRNA removal and poly(A) retention for construction of the libraries.
The percentage of embryos to reach the 32-cell stage at 24 h post fertilization was significantly correlated to family eyeing rate, indicating that inviable embryos were developmentally compromised before zygotic genome activation. RNA sequencing identified 2 differentially expressed transcripts (DETs) from total mRNA sequencing comparing females with low-quality (< 5% eyeing), medium-quality (30–50% eyeing), and high-quality (> 80% eyeing) eggs. In contrast, RNA sequencing from poly(A) captured transcripts identified 945 DETs between low- and high-quality eggs, 1012 between low- and medium-quality eggs, and only 2 between medium- and high-quality eggs. The transcripts of mitochondrial genes were enriched with polyadenylated transcript sequencing and they were significantly reduced in low-quality eggs. Similarly, mitochondrial DNA was reduced in low-quality eggs compared with medium- and high-quality eggs. The functional gene analysis classified the 945 DETs between low- and high-quality eggs into 31 functional modules, many of which were related to ribosomal and mitochondrial functions. Other modules involved transcription, translation, cell division, apoptosis, and immune responses.
Our results indicate that differences in egg quality may be derived from differences in maternal nuclear transcript activation and cytoplasmic polyadenylation before ovulation, as opposed to accumulation and storage of maternal nuclear transcripts during oogenesis. Transcriptome comparisons suggest low-quality eggs suffered from impaired oxidative phosphorylation and translation. The DETs identified in this study provide insight into developmental competence in rainbow trout eggs.
Reliable production of high-quality eggs is essential for meeting production cycle demands for seed stock. Fertility is high in the rainbow trout industry when fish are maintained under optimal conditions. Nevertheless, quality of the eggs or ova can be affected by many intrinsic and extrinsic factors including the genetics, age and diet of brood fish [1,2,3,4,5,6]; pre-spawning exposure to stressors and photo-thermal cycles [7,8,9,10]; and postovulatory aging of the eggs [11,12,13]. Understanding mechanisms by which egg quality becomes compromised in response to suboptimal genetics, management, nutrition, and environmental conditions is critical to optimizing hatchery productivity.
The oocyte becomes transcriptionally inactive following oocyte growth and remains transcriptionally silent or greatly repressed until zygotic genome activation (ZGA) which usually takes place around the time of the mid-blastula transition (MBT) in most vertebrates. Therefore, the oocyte serves as a reservoir of biomolecules including proteins, lipids and RNAs deposited into the egg during oogenesis, for utilization from oocyte maturation through early embryonic development [14, 15]. Levels of certain proteins and lipids have been linked to egg viability in many fish species including rainbow trout . A relationship between the maternal transcriptome and developmental competence has also been supported in a variety of fishes using an assortment of molecular approaches, although few investigations involved next-generation sequencing [17, 18]. In rainbow trout, genes linked to decreased egg quality caused by postovulatory aging were identified using quantitative reverse transcription PCR [19,20,21] and transcripts associated with decreased egg quality in response to the use of photoperiod to shorten the time to spawning or hormone-induced ovulation were identified by microarray analyses [22, 23]. One microarray study compared the transcriptome of un-manipulated female rainbow trout that exhibited either 100% viability through eyeing or less than 64%, but confirmed only one differentially expressed transcript (DET) . Less than 200 DETs were identified among the studies with few overlapping transcripts between them. Although deep sequencing approaches have been used to identify some miRNAs and mitochondrial genome encoded small RNAs related to egg deterioration due to postovulatory aging [25, 26], global mRNA analysis techniques based on deep-sequencing technologies have not been applied to investigate a possible connection between mRNA content in unfertilized eggs and egg quality in fish.
Stored maternal transcripts generally have shortened polyadenylic acid (poly(A)) tails and are masked to inhibit both translation and degradation [27,28,29,30]. In the oocyte, new transcripts intended for sequestration are polyadenylated in the nucleus and then translocated to the cytoplasm where they are subsequently partially deadenylated for storage [30,31,32]. The stored maternal transcripts require cytoplasmic polyadenylation to allow translation. The cytoplasmic polyadenylation or deadenylation of stored transcripts is a critical control mechanism for translation during oocyte maturation and early embryonic development [28, 33, 34]. When comparing eggs of different quality, differences in expressed transcripts based on total mRNA content may indicate differences in accumulation or degradation of the transcripts throughout oogenesis and maturation, whereas differences in mRNAs with longer poly(A) tails may indicate differences and changes in the translational activity of the transcripts [35,36,37,38]. To further understand the relationship between the maternal transcriptome and egg quality, we used RNA-seq to compare mRNA transcriptomes of ovulated eggs from 20 individual females that produced eggs of disparate quality as determined by eyeing rate. We compared transcriptomes sequenced from libraries prepared following rRNA removal or by oligo(dT) capture of polyadenylated RNA. Oligo(dT) capture methodologies are not efficient at capturing transcripts with short poly(A) tails [35, 39, 40], and therefore, mRNA transcriptomes sequenced from libraries prepared by oligo(dT) capture of polyadenylated RNA should be enriched in activated transcripts with elongated poly(A) tails.
Results and discussion
Eyeing rate and early embryo viability
Viability was assessed at ~ 250 accumulated thermal units (ATUs) post fertilization, which we refer to as eyeing in the present manuscript. The 250 ATU mark is actually after retinal pigmentation but is often used by industry because mortality is generally very low after retinal pigmentation ; it is well after embryos are resistant to mechanical shock [42, 43], and is several days before hatching which allows time for the eggs to be sorted to remove dead and subviable eggs before shipment to production facilities. Eyeing rates for the families in the selective breeding program was typical of rainbow trout aquaculture operations , exhibiting a mean eyeing rate of 79.3% (Fig. 1). Only 16 of the 192 families evaluated had less than 50% eyeing and were classified as subfertile. Furthermore, only six families exhibited an eyeing rate of less than 10%, whereas eyeing rates were greater than 30% for all others. Sperm used to fertilize each of the subfertile families also yielded families with eyeing rates over 78%, substantiating the eggs and not the sperm as the cause of the subfertility. Visual inspection of eggs collected before fertilization did not show obvious signs of the eggs being compromised in ways that would allow for their being discarded by hatchery personnel.
Embryo cleavage was assessed for the 20 selected families at about 24 h post fertilization (Table 1). On average about 80% of the embryos in the high-quality families reached at least the 32-cell stage at this time. The percentage of embryos reaching the 32-cell stage was similar to the eyeing rates for most families, suggesting that most of the embryos that would not survive had died or had very delayed development by the 32-cell stage. The percentage of embryos reaching the 32-cell stage was significantly correlated to eyeing rate (R = 0.85; P < 0.001). Few embryos failed to reach the 8-cell stage by 24 h post fertilization. Together these data support most of the non-viable embryos were fertilized but failed developmentally before the 32-cell stage, which should be before the major wave of ZGA. Although the timing of the major wave of ZGA has not been characterized in rainbow trout, in general, more cleavage divisions are completed before ZGA in animals that develop more slowly . The ZGA for most fish species investigated, all of which develop more rapidly than rainbow trout, have been shown to begin during the MBT at about cell cycle 10 or ~ 1000 cells [45,46,47,48,49]. An exception is the medaka (Oryzias latipes) in which the ZGA begins at about the 64-cell stage; before the MBT . A previous study on the same rainbow trout broodstock population as in the present study reported most embryo mortality in subfertile families took place by the second cleavage interval ; earlier than in our study and well before ZGA. Although the timing of embryonic mortality cannot specify a cause of the mortality, most of the embryonic mortality reported for this population of rainbow trout has taken place while the embryos were dependent upon maternal transcripts and before many of the other forms of stored biomolecules would be either required or exhausted by the embryos, leaving aberrations in levels or activation of maternal transcripts as possible contributors to subfertility.
Mapping of sequencing reads following rRNA removal and poly(a) retention
RNA-seq analysis of the eggs of 20 females generated 31 to 58, and 43 to 68 million reads from the libraries constructed with the Illumina® TruSeq® Stranded Total RNA Library Prep Kit with Ribo-Zero Gold (rRNA removal) and TruSeq® Stranded mRNA Sample Prep Kit (poly(A) retention) respectively (Table 2). Similar percentages of the reads were mapped to rRNA gene sequences using the two approaches, with an average of 5.4% derived from the libraries constructed by rRNA removal and 4.7% derived from the libraries constructed by poly(A) retention kits respectively. Whereas the percentage of reads mapped to the nuclear transcriptome was half as great with poly(A) retention compared with rRNA removal, the percentage of reads mapped to mitochondrial RNA was more than 10-fold greater for poly(A) retention compared with rRNA removal. This reversal in trends is consistent with a higher proportion of the mitochondrial mRNA transcripts being polyadenylated than the proportion of the nuclear mRNA transcripts that are polyadenylated or efficiently captured by the poly(A) retention procedures.
Mitochondrial mRNA transcripts are polyadenylated with a tail of approximately 50 nucleotides in vertebrates as part of transcript processing and therefore most would be expected to be captured by the poly(A) retention [52, 53]. Similarly, nuclear mRNA transcripts are also polyadenylated as part of processing. In most cells the majority of cytosolic nuclear transcripts are polyadenylated with a poly(A) tail greater than 80 nucleotides [29, 54]. However, stored maternal nuclear transcripts possess a short poly(A) tail around 15–40 nucleotides that are elongated to over 80 nucleotides through cytoplasmic polyadenylation during activation [27, 29, 36, 40, 55]. In general, oligo(dT) capture approaches are not very efficient at capturing mRNAs with shorter poly(A) tails [35, 39, 40]. As far as we are aware, the capture efficiency of the Illumina® TruSeq® Stranded mRNA Sample Preparation Kit for short poly(A) tails such as those in stored maternal mRNA, has not been characterized. It is therefore likely that data collected following poly(A) retention represent primarily activated mRNAs with longer poly(A) tails and data following rRNA removal represent primarily the more abundant stored maternal transcripts with short poly(A) tails, which is consistent with a lower percentage of the transcript reads aligning to the nuclear transcriptome following poly(A) retention compared with rRNA removal.
Eggs of vertebrates possess large stockpiles of mitochondria that are active in providing energy during maturation and stockpiling ATP for energy to drive early embryonic events including cleavage [56,57,58,59]. Furthermore, there is extensive mitochondrial DNA (mtDNA) replication during maturation that then ceases until after ZGA. High proportions of mitochondrial transcripts as observed in our study have been reported in early cod and halibut embryos [48, 60]. The percentage of reads that were aligned as mitochondrial transcripts was reduced in the low-quality eggs following poly(A) retention as well as rRNA removal. The percentage of reads aligned to mitochondrial RNAs in the low-quality eggs following poly(A) retention was 23.6% compared with 47.2% in the high-quality eggs, and 2.4% compared with 3.9% respectively following rRNA removal (Table 2). A reduction in mitochondria and mtDNA has been identified as a cause of reduced fertility in several mammalian species and serves as a marker of oocyte quality in mammals [56, 57]. Furthermore, augmentation of mitochondrial number or mtDNA in maturing oocytes has been recognized as a method to improve oocyte quality in mammals including humans [56, 57, 61, 62]. We therefore used real-time quantitative PCR to measure abundance of mtDNA in our samples by measuring two mitochondrial genes, mt-atp6 and mt-cyb. Abundance of both genes were reduced 1.6 log2 fold change (log2FC) in the low-quality eggs (Fig. 2) suggesting a reduction in mitochondria or mtDNA may contribute to differences in the percentage of reads aligned to the mitochondrial transcriptome among egg quality groups.
Gene expression profiling revealed greater differences in polyadenylated mRNA than total mRNA abundance with fertility
We detected 44,330 and 39,133 transcripts with at least three normalized reads expressed in the libraries constructed by rRNA removal and poly(A) retention kits respectively. About 89% of the transcripts following rRNA removal and 77% of the transcripts following poly(A) retention were shared among the three treatment groups (Fig. 3) and many of those that were not shared had low transcript numbers. No DETs were identified comparing the 10 subfertile and 10 fertile or high-quality females from either the rRNA removal or poly(A) retention libraries by DESeq2 with the criteria of a false discovery rate (FDR) < 0.05. We therefore divided the subfertile group into a low- and medium-quality group (Table 1).
Only two DETs were identified among the three treatment groups; low-, medium-, and high-quality eggs, from our dataset sequenced from the libraries constructed with an rRNA removal kit by DESeq2 with the criteria of log2FC ≥ 1, FDR ≤ 0.05, and a total number of normalized reads in a given comparison being ≥500 (Additional file 1: Table S1). The transcript agfg1 (Arf-GAP domain and FG repeats-containing protein 1) was enriched in low-quality eggs compared with high-quality eggs (1.81 log2FC, FDR = 0.0172) and ahnak (neuroblast differentiation-associated protein ahnak-like) was enriched in low-quality eggs compared with medium-quality eggs (2.01 log2FC, FDR = 0.0026). AGFG1 has been shown to affect the accumulation of a small set of non-polyadenylated cellular mRNAs in mammalian cells , and mRNA of Drongo (Drosophila neural GTS1-like), the homolog of AGFG1 in the fruit fly, associates with Me31B (maternal expression of 31B) which is required for the translational repression of maternal mRNAs in the oocyte . AHNAK is considered to be involved in calcium flux regulation and has been proposed to interact with s100 proteins to regulate cellular Ca2+ homeostasis . Notably, an s100 protein, s100a1 (S100 calcium binding protein A1), is among the most reduced polyadenylated transcripts in low-quality eggs. Most important, the lack of DETs following rRNA removal supports the concept that low fertility was not a result of differences in the accumulation or degradation of transcripts throughout oogenesis or maturation.
A total of 1339 transcripts were differentially expressed among the egg quality groups when the libraries were constructed with a poly(A) retention kit, using the same criteria as with rRNA removal (Additional file 1: Table S1). There were only 2 DETs comparing medium- and high-quality eggs, esr2b (estrogen receptor beta 2) and pltp (phospholipid transfer protein), both of which were higher in medium quality eggs compared with both low- and high-quality eggs. There were 945 DETs comparing low- and high-quality eggs, and 1012 comparing low- and medium-quality eggs. There were 619 DETs shared between the low- versus high-quality egg comparison and the low- versus medium-quality comparison, and both DETs from the medium- versus high-quality comparison were also differentially expressed between the low- and medium-quality eggs. The shared DETs between the low- versus high- quality comparison and the low- versus medium-quality comparisons were consistent in direction of change for any given transcript. Taken together, there was little difference in the transcriptomes of medium-quality eggs with between 30 and 50% eyeing rates, and the high-quality eggs with greater than 80% eyeing rates; therefore, we focused our further analyses and discussions on the low- versus high-quality comparison. The 945 DETs included 732 unique gene descriptions, with many of the more differentially expressed transcripts sharing gene descriptions with multiple transcripts or locus tags.
Among those 945 DETs, 724 were decreased and 221 were increased in low-quality eggs (Additional file 1: Table S1). There were many nuclear genes associated with ribosome production and function among the most significantly reduced DETs. Over 100 of the 945 DETs were for 40s or 60s ribosomal proteins. In general, the functions of the transcripts were diverse even among the most differentially expressed transcripts (Fig. 4). The transcripts in Fig. 4 were associated with 62 different GO terms for Biological process and 30 for Molecular function, even with 14 of these transcripts having no associated GO terms. Among the 221 enriched transcripts in the low-quality eggs, the gene transcript with the largest increase, 2.07 log2FC, is tob1 (protein tob1-like) (Additional file 1: Table S1). This transcript, however, was highly variable among the low-quality groups (FDR = 0.0010). Whereas the high-quality eggs ranged from 46 to 280 normalized reads, one of the low-quality groups with no survival at eyeing and another with 1.3%, had over 1300 normalized reads each, whereas the other two had below 80 reads. TOB1 is characterized as an anti-proliferative protein whose activity is mediated through interactions with the Caf1a/Caf1b deadenylases leading to target mRNA deadenylation and decay [66,67,68]. Considerable differences in expression among the groups of low-quality eggs, including a transcript involved with mRNA deadenylation and decay, might suggest disparate causes leading to the reduction in quality or differences in the progression of viability among the low-quality families.
Because egg quality is a major issue in the aquaculture industry , many microarray and quantitative reverse transcription PCR studies have revealed sets of genes associated with egg quality . In rainbow trout, eggs with reduced egg quality in response to postovulatory aging had lower expression levels of tubb (tubulin b) and npm2 (nucleoplasmin-2) , and eggs with reduced quality in response to photoperiod treatment to accelerate the time of spawning had increased expression of pyc (pyruvate carboxylase) . Consistent with those findings, our data show transcript levels of tubb and npm2 are decreased, and pyc is increased, in low-quality eggs. On the other hand, we observed decreased expression of krt8 (keratin 8), krt18 (keratin 18), rpl24 (60s ribosomal protein l24), and apoc1 (apolipoprotein c1) with reduced egg quality, whereas krt8 and krt18 have been reported to increase with reduced egg quality due to postovulatory aging , and rpl24 and apoc1 to increase with reduced egg quality in response to hormone implantation or photoperiod manipulation to accelerate the time of spawning . The three studies used poly(A) retention and therefore likely measured changes in activated transcripts. In the present study total transcript abundance based on rRNA removal libraries did not significantly differ between high- and low-quality eggs for any of the seven genes mentioned, with the mean normalized reads differing by less than 10% for rpl24, pyc, npm2, and tubb, and 68, 80 and 146% for krt8, krt18 and apoc1 respectively. Moreover, the directions of the trends were similar with poly(A) retention and rRNA removal for each of these transcripts. The differences in the direction of the expression of specific genes among studies may be due to the treatments behind the reductions in egg quality or criteria for assessing developmental competence. Alternatively, differences may reflect where in the activation and deactivation cycle of the transcripts the eggs were sampled, and not just a measure of how many transcripts are activated overall. Further study of the kinetics of cytoplasmic adenylation and deadenylation in rainbow trout may help with interpretations of transcript levels and their associations with egg quality.
Even with poly(A) retention, only two transcripts were differentially expressed between medium- and high-quality eggs (esr2b and pltp) despite the medium-quality eggs ranging from 31 to 48% eyeing, and also representing only the lower 8% of clutches in terms of eyeing rate. The lack of DETs between medium- and high-quality eggs may suggest transcript variance did not contribute to the reduced egg quality, or transcript disparities were not detected due to the nature of the study or data analyses. Transcript signatures associated with less drastic differences in fertility have been identified in rainbow trout eggs in response to specific treatments [20, 22, 23]. Since the cause or causes of the reduced egg quality in the present study is undefined, there is likely a host of maladies among the egg clutches resulting in reduced fertility that do not share the same molecular pathways, making transcriptome signatures of any one difficult to discern. More robust analyses such as use of artificial neural networks and supervised machine learning which has been able to identify molecular signatures composed of many minor changes in transcript levels in oocytes of striped bass that associated with fertility  may be helpful. Regardless, the expressions of many genes have been linked to egg quality in the present study including genes not previously associated with egg quality, providing new insight into how an egg can become compromised.
Mitochondrial RNA expression
The rainbow trout mitochondrial genome encodes 13 polypeptides, 22 tRNAs, two rRNAs and has a noncoding D-loop region . Transcripts for the 13 mitochondrial protein coding genes and D-loop region were found to be significantly reduced in the low-quality eggs by RNA-seq with poly(A) retention libraries (Table 3). In addition, mt-tn (tRNA-Asn) was also decreased in low-quality eggs but few transcripts were detected for this or all other tRNA genes. There were also abundant reads for the rRNAs; mt-rnr1 (mitochondrion 12S) and mt-rnr2 (mitochondrion 16S), and despite the trends being similar for the rRNAs and mRNAs, the DESeq2 analysis did not detect a significant reduction in mt-rnr1 or mt-rnr2 in low-quality eggs (FDR > 0.1; Table 3). Although there are some reports of mitochondrial rRNAs being polyadenylated [72, 73], mitochondrial rRNAs are generally not polyadenylated or contain no more than 10 nucleotides of the tail in most vertebrate cells [52, 53, 74]. It is worth noting that whereas the rRNAs comprised about half of the mitochondrial RNA reads following rRNA removal (Additional file 1: Table S2), they comprised less than 5% of the mitochondrial RNA reads following poly(A) retention, supporting the idea that few of the rRNA transcripts were polyadenylated or possessed long poly(A) tails. We do not know if the antibodies used for rRNA removal recognize mt-rnr1 or mt-rnr2 to any extent. Mitochondrial rRNA transcripts are polyadenylated as part of the degradation process  and therefore the differences in polyadenylated rRNA transcripts among groups may represent differences in transcripts undergoing degradation.
The changes in expression of the 13 mitochondrial protein genes between high-quality and low-quality eggs are very consistent, ranging between − 1.34 and − 2.10 log2FC as calculated by DESeq2 (Table 3). This consistency may be due to coordinated regulation of the transcripts as they are all needed for oxidative phosphorylation [52, 57], or to a reduction in mitochondria in the low-quality eggs. However, DESeq2 uses shrinkage estimators  that resulted in reduced log2FC estimates for many DETs compared with calculating log2FC based only on normalized read values. If log2FC for the 13 protein mRNAs are calculated based simply on the normalized reads the average is about − 2.21 log2FC for the low-quality eggs compared with the high-quality eggs (Additional file 1: Table S2). This magnitude of change is greater than the − 1.6 log2FC difference in mitochondrial DNA (Fig. 2), supporting an additional reduction in expression beyond just a reduction in mitochondrial number although the measurements were conducted using different assays. It is also worth noting that using the same calculation approach, mt-rnr1 and mt-rnr2 are also about a 2.3 log2FC, suggesting similar circumstances behind the reduced levels of the rRNAs and the mRNAs in low-quality eggs.
Interestingly DESeq2 analysis of data from rRNA removal libraries revealed no significant differences in mitochondrial gene expression among egg quality groups (FDR > 0.999). Nevertheless, the numerical means for the protein and D-loop mRNAs, and rRNAs, were lower in the low-quality eggs compared to the high-quality eggs. Again, it is worth noting that whereas there was an average − 0.17 log2FC between high- and low-quality eggs for these transcripts with DESeq2 (Additional file 1: Table S1), the log2FC calculated based simply on the normalized reads averaged − 0.78 (Additional file 1: Table S2). This is in line with the differences in the proportion of aligned mitochondrial transcripts among treatments with rRNA removal libraries (Table 2). Differences in percentage of reads aligning with mitochondrial genes among groups are not consistent with a − 1.6 log2FC difference in mitochondrial DNA or mitochondrial numbers between low- and high-quality eggs. Among possible explanations for differences in mtDNA being greater than differences in total transcript reads for mitochondrial genes may lie in differences in the timing of transcription and mtDNA replication. During maturation there is an explosion in mtDNA replication, increasing 1000-fold in some species . It is not known when the transcripts being measured in ovulated eggs are transcribed or processed in relation to this increase in mtDNA replication. Perhaps a lower content of mtDNA in low-quality eggs has not yet been reflected in transcription rates. Alternatively, there could be a compensatory increase in mitochondrial transcription in those oocytes with reduced mtDNA. The greater reduction in polyadenylated mitochondrial transcripts than mtDNA in low-quality eggs may reflect a combination of less mtDNA and transcript processing within mitochondria.
The requirement for oxidative phosphorylation by mitochondria to provide ATP in the early embryo varies among species [56, 76]. In the zebrafish, the maternal ATP pool is insufficient to execute the ubiquitin proteasomal pathway required for protein degradation required to advance beyond the 32-cell stage . Although the 32-cell stage is well before ZGA in this species [45, 47], mitochondrial transcription has been shown to be active prior to mtDNA replication in zebrafish embryos . Mitochondria in zebrafish embryos are active and free fatty acids serve as substrate for oxidative phosphorylation to supply the required ATP . Less is known about how rainbow trout meet energy demands in the early embryo but oxidative metabolism is present even in the unfertilized eggs, and continues through early development . Although the timing of the mortality in zebrafish embryos deficient in the ability to produce ATP is similar to that of the present study, egg ATP levels were found to not correlate with fertility in rainbow trout eggs [10, 77]. The associations among mitochondrial mtDNA abundance, mitochondrial transcript levels, ATP levels, and egg quality in rainbow trout are unresolved. Mitochondria serve many functions in addition to ATP production that are essential to embryo survival such as sequestration and release of intracellular calcium [78, 79]. Furthermore, deficiencies in mitochondrial activity and mtDNA number have been shown to have separate although overlapping impacts on egg quality [56, 80].
Functional classification of differentially expressed genes
Gene ontology (GO) analysis of the 945 DETs from the comparison of polyadenylated transcript enriched libraries for low- versus high-quality eggs revealed one or more associated GO terms for 811 of the transcripts. The five most common GO terms in Biological process are ribosome biogenesis, translation, metabolic process, oxidation-reduction process, and DNA templated regulation of transcription (Fig. 5). Clustering of DETs based on associated GO terms, using the Database for Annotation, Visualization and Integrated Discovery (DAVID) gene functional classification algorithms under kappa value of 0.3, resulted in 31 functional modules [81, 82] (Table 4). The enrichment scores of the modules ranged from 1.52 to 13.40 with the number of genes in each module ranging from 4 to 164, which included 547 of the DETs in total.
The largest cluster included 164 transcripts involved in ribosome biogenesis of which only 7 were increased in low-quality eggs. As mentioned, over 100 DETs are for 40s and 60s ribosomal proteins. Moreover, the categories of the GO terms in gene clusters 2, 3, 4, 5, 7, 8, 9, and 31, are also associated with ribosome function. Clusters 2, 12, 14, 17, and 23 are involved in mitochondrial function. Together, many of the clusters are associated with translation and the production of energy by mitochondria to drive early cell division. The other gene clusters are mainly involved in regulation of transcription, cell division, apoptosis, and immune responses. Transcripts increased in low-quality eggs were distributed among 22 clusters.
There were 398 DETs that were not included in the clusters. This includes the 134 DETs without associated GO terms such as n1/n2 protein (histone-binding protein n1 n2-like), senp7 (sentrin-specific protease 7-like) and parp8 (kisutch poly [ADP-ribose] polymerase 8-like); and 264 that were orphan genes such as mettl3 (n6-adenosine-methyltransferase 70 kda subunit), bmp10 (bone morphogenetic protein 10-like), and haus3 (haus augmin-like complex subunit 3-like). The listed transcripts are among the top DETs based on FDR values and therefore may also serve as important indicators of egg quality (see Fig. 4).
The present study identifies differences in the transcriptome among ovulated eggs of different quality for which most of the mortality occurred between fertilization and the 32-cell stage, which is before ZGA. The identification of only two DETs by RNA-seq of libraries constructed by rRNA removal kits, compared with 1339 DETs derived from libraries following poly(A) retention kits, supports transcriptome differences with egg quality arose from differences in cytoplasmic polyadenylation or deadenylation of stored maternal transcripts as opposed to being the result of differences in the accumulation of maternal transcripts during oogenesis. Furthermore, few DETs were identified between medium- and high-quality eggs. Nine clusters of DETs which encompassed 375 DETs or about 40% of all DETs identified between high- and low-quality eggs, were associated with ribosome biogenesis and processing. The multitude of ribosome related DETs in low-quality eggs suggests inadequate ribosome production required for maternal mRNA translation, which could lead to a cascade of developmental dysfunction. This reduction in ribosomal gene expression was true of mitochondrial transcripts as well as nuclear transcripts. Moreover, mtDNA abundance was reduced by 1.6 log2FC in low-quality eggs compared with high-quality eggs further supporting the ability of the egg to provide energy following fertilization was compromised. In addition to genes associated with ribosome and mitochondrion biogenesis and function, GO analysis indicates levels of transcripts involved in the regulation of transcription, translation, cell division, apoptosis, and immune responses were altered in the low-quality eggs. Many of these genes have not previously been reported to contribute to egg quality in rainbow trout. The present study provides insights into how dysfunction of the egg transcriptome can affect developmental competence in fish eggs.
Eggs were collected from rainbow trout that were part of the selective breeding program at Troutlodge Inc. Sumner, WA, USA. Eggs from individual two-year-old broodstock rainbow trout were stripped into plastic bags. About 90 unfertilized eggs from each female were collected and immediately frozen in liquid nitrogen, and another 50 eggs were collected and placed into modified Davidson’s fixative  for examination to eliminate samples with overripe eggs or other abnormalities. The remaining eggs were fertilized with sperm harvested from neomales. The semen derived from each sire was used to fertilize eggs from two to three females. The fertilized eggs were incubated as individual families as part of the Troutlodge Inc. selective breeding program which evaluates eyeing rate at about 250 ATUs calculated as the sum of mean daily water temperature in degrees Celsius, which we refer to as eyeing in the present manuscript. The eggs were incubated at 10 °C for about the first 24 h post fertilization, after which time a sample of about 25–60 embryos were collected from each family and fixed in Stockard’s solution  to evaluate early embryonic survival and viability at about the 32-cell stage by enumerating the embryos reaching each stage of cell cleavage (Table 1). The frozen samples were kept in a − 80 °C freezer at Troutlodge Inc. until they were shipped on dry ice to the National Center for Cool and Cold Water Aquaculture (NCCCWA) after which they were again placed at − 80 °C until RNA isolation. The fixed samples were shipped at ambient temperature to NCCCWA for evaluation.
Selection of rainbow trout females for RNA-seq analysis of eggs
Selection of egg samples for RNA-seq analysis was based primarily on eyeing rate. A range of 130–218 individuals from each family were examined for survival and viability at eyeing. Dead and subviable eggs included those that were unfertilized, had precipitated yolk in response to shocking the eggs, or were considered to have poorly developed eyes. We considered an eyeing rate of 50% as the demarcation between fertile and subfertile families. Only 16 of the 192 families generated had eyeing rates that were less than 50% (Fig. 1). In addition, the sperm from the sires used for each of these matings also yielded at least one family with an eyeing rate greater than 78% confirming the sperm used for fertilization was not the cause of the poor eyeing rates. Ten of the subfertile families (0–47.5% eyeing) and 10 fertile or high-quality families with eyeing rates greater than 80% (80.6–97.7% eyeing) that shared sires with the ten subfertile families, were selected for RNA-seq analyses. Since there were no families with eyeing rates between 10 and 30%, we subsequently further divided the subfertile families as low- (0–4.2%) and medium- (30.6–47.5%) quality families or females. Visual examination of the fixed eggs from these 20 females revealed no obvious signs of poor egg quality before fertilization.
Assessment of early embryo development
All embryos collected at about 24 h post fertilization from each of the females selected for RNA-seq analysis were examined to determine viability. The fixed embryos were immersed in 0.5% methylene blue overnight. The cell number of each embryo was counted or confirmed to be greater than 32, using a stereo microscope (Nikon SMZ660). Those embryos with less than 32 cells were considered subviable.
RNA isolation and sequencing
RNAs were isolated from frozen eggs which were homogenized in Tri Reagent (Sigma, St. Louis, MO) with a Qiagen Retsch MM300 TissueLyser Shaker Mixer Grinder Agitator Mill (Retsch Inc., Haan, Germany). Total RNA was isolated following the manufacturer’s protocol with the modification of using Phase Lock Gel (5 PRIME, Inc., Gaithersburg, MD, USA) and Phase Separation Reagent (Molecular Research Center, Cincinnati, OH, USA) to separate the aqueous phase from the organic phase. The isolated RNAs were further purified by lithium chloride precipitation and treated with DNase. The RNA integrity was evaluated by gel electrophoresis (Additional file 2: Figure S1), a NanoDrop ND-1000 (Thermo SCIENTIFIC, Wilmington, DE, USA) with λ260/280 great than 1.98, and a 2100 Bioanalyzer (Santa Clara, CA, USA) with RNA Integrity Number between 7.2 and 9.1 (Additional file 1: Table S3). DNaseI treated RNA was used to construct libraries with the Illumina® TruSeq® Stranded Total RNA Library Prep Kit with Ribo-Zero Gold (rRNA removal) and TruSeq® Stranded mRNA Sample Prep Kit (poly(A) retention). The libraries were sequenced by HiSeq2500 with 100 nt paired-end reads. Raw reads were deposited in NCBI Sequence Read Archive database (SRA accession: SRP108797). We sequenced the libraries for rRNA removal in two batches. The first batch of four samples (99, 119, 108, and 129) had less than 4% of the reads aligned to rRNA whereas the second batch of 16 samples contained 14.4 to 55.7% rRNAs indicating the rRNAs were not effectively removed (Additional file 1: Table S4). Hence, the 16 samples were re-sequenced using a newly purchased kit which resulted in the effective removal of rRNAs to below an average of 6% of the reads. The data from the three rRNA removal sequencing runs were included in subsequent transcriptome analyses.
To classify reads as belonging to the nuclear transcriptome, the mitochondrial transcriptome, or as rRNA, after the adaptor was trimmed by bcl2fastq v220.127.116.11, the reads passing FastQC evaluation (http://www.bioinformatics.babraham.ac.uk/projects/fastqc), were respectively aligned to the nucleotide sequences corresponding to all CDS features annotated on the rainbow trout genome assembly (GCA_900005705.1)  with 72 additional genes selected from gene bank (Additional file 1: Table S5), the mitochondrial genome which includes 38 genes and the D-loop region [25, 86] and rRNA genes  by using Bowtie2 tool under default settings . To avoid redundancy between the transcripts assigned to the mitochondrial and nuclear transcriptomes, we used the mitochondrial transcripts as queries to blast the rainbow trout  transcriptome data mentioned above. The blast results showed that GSONMT00007417001 was aligned with mitochondrial mt-rnr2, mt-tl1, and mt-nd1, and GSONMT00007419001 aligned with mt-nd4l. Therefore, GSONMT00007417001 and GSONMT00007419001 were removed from the nuclear transcriptome reference list. In addition, although mt-rnr1 and mt-rnr2 transcripts were included in the mitochondrial transcriptome, some of the rRNAs in the rRNA reference list have high similarity to these genes and therefore some of the reads aligned under rRNA may have been mt-rnr1 and mt-rnr2. The raw reads aligned to nuclear transcriptome and mitochondrial mRNAs were merged and used as input to DESeq2 to identify DETs among groups at a FDR < 0.05, log2FC greater than 1, and total number of normalized reads within comparisons being greater than 500.
Gene ontology analysis of the identified DETs was conducted using Blast2GO PRO platform (BioBam Bioinformatics S.L., Spain) . An R script written according to the algorithms proposed by the Database for Annotation, Visualization and Integrated Discovery (DAVID) gene functional classification was used to cluster the identified DETs into functional gene modules based on the results of the GO analyses [81, 82]. Module enrichment scores were generated by calculating the geometric mean of the P-values which were derived from hypergeometric test of the input gene sets, followed by negative log transformation of the geometric mean . False discovery rate was calculated using the Benjamini-Hochberg procedure. Spearman correlation coefficient for eyeing rate and early embryo viability was estimated using R (R × 64 3.3.0).
Mitochondrial gene quantification by real-time quantitative PCR
The insoluble materials leftover following homogenization in Tri Reagent during RNA isolation was mixed with 180 μl of 1 × TE buffer and shaken at speed setting 30 for 2 min with a Qiagen TissueLyser and then incubated for 10 min at room temperature. The mixture was centrifuged at 12,000 rpm for 10 min, and the supernatant was transferred to a tube containing Phase Lock Gel. The tube was shaken and centrifuged as above. The supernatant was collected, and DNA was isolated using a Quick-DNA™ Universal Kit (Zymo Research, Irvine, CA, USA). The DNA isolate was digested with RNaseA and the DNA was purified using a ZR-Duet™ DNA/RNA MiniPrep kit (Zymo Research, Irvine, CA, USA). The DNA from about 90 eggs was eluted in 20 μl of elution buffer and then diluted by adding 20 μl of water. The relative quantity of the mt-atp6 and mt-cyb genes was measured on an ABI 7900HT sequence detection system (Applied Biosystems, Foster City, CA). Each reaction consisted of 1.5 μL of diluted DNA, 3 μl of each primer (5 μM) and 1× SYBR Green PCR Master Mix (Applied Biosystems, Foster City, CA). The thermal cycling profile was 50 °C for 2 min, 95 °C for 10 min and 40 cycles of 95 °C for 30 s, 60 °C for 20 s, and 72 °C for 30 s. A final dissociation step was performed to assess the specificity of the reaction. Relative quantification of the mitochondrial DNA was estimated by the standard curve method with three technical replications, and mean differences of the mitochondrial DNA were reported as relative change using the value for the high-quality eggs as a calibrator. An 18 s rRNA gene was used as the reference control . Primer sequences are shown in Additional file 1: Table S6.
Accumulated thermal units
Database for Annotation, Visualization and Integrated Discovery
Differentially expressed transcript
False discovery rate
Log2 fold change
National Center for Cool and Cold Water Aquaculture. Gene name abbreviations are presented in Additional file 1: Table S7
Zygotic genome activation
Brooks S, Tyler CR, Sumpter JP. Egg quality in fish: what makes a good egg? Rev Fish Biol Fish. 1997;7(4):387–416.
Vehvilainen H, Kause A, Koskinen H, Paananen T. Genetic architecture of rainbow trout survival from egg to adult. Genet Res. 2010;92(1):1–11.
Su GS, Liljedahl LE, Gall GAE. Genetic and environmental variation of female reproductive traits in rainbow trout (Oncorhynchus mykiss). Aquaculture. 1997;154(2):115–24.
Blom JH, Dabrowski K. Reproductive success of female rainbow-trout (Oncorhynchus-Mykiss) in response to graded dietary Ascorbyl monophosphate levels. Biol Reprod. 1995;52(5):1073–80.
Bromage NRaC, P.R.T. (ed.): Egg production in rainbow trout. Croom Helm., London; 1988.
Palace VP, Werner J. Vitamins a and E in the maternal diet influence egg quality and early life stage development in fish: a review. Sci Mar. 2006;70:41–57.
Contreras-Sanchez WM, Schreck CB, Fitzpatrick MS, Pereira CB. Effects of stress on the reproductive performance of rainbow trout (Oncorhynchus mykiss). Biol Reprod. 1998;58(2):439–47.
Campbell PM, Pottinger TG, Sumpter JP. Stress reduces the quality of gametes produced by rainbow-trout. Biol Reprod. 1992;47(6):1140–50.
Bonnet E, Fostier A, Bobe J. Characterization of rainbow trout egg quality: a case study using four different breeding protocols, with emphasis on the incidence of embryonic malformations. Theriogenology. 2007;67(4):786–94.
Aegerter S, Jalabert B. Effects of post-ovulatory oocyte ageing and temperature on egg quality and on the occurrence of triploid fry in rainbow trout, Oncorhynchus mykiss. Aquaculture. 2004;231(1–4):59–71.
Lahnsteiner F. Morphological, physiological and biochemical parameters characterizing the over-ripening of rainbow trout eggs. Fish Physiol Biochem. 2000;23(2):107–18.
Springate JRC, Bromage NR, Elliott JAK, Hudson DL. The timing of ovulation and stripping and their effects on the rates of fertilization and survival to eying, hatch and swim-up in the rainbow-trout (Salmo-Gairdneri R). Aquaculture. 1984;43(1–3):313–22.
Craik JCA, Harvey SM. Egg quality in rainbow-trout - the relation between egg viability, selected aspects of egg composition, and time of stripping. Aquaculture. 1984;40(2):115–34.
Tadros W, Lipshitz HD. The maternal-to-zygotic transition: a play in two acts. Development. 2009;136(18):3033–42.
Lyman-Gingerich J, Pelegri F. Maternal factors in fish oogenesis and embryonic development. In: Babin PJ, Cerda J, Labadie K, editors. The fish oocyte: from basic studies to biotechnological applications. Dordrecht, The Netherlands: Springer; 2007. p. 141–74.
Lubzens E, Bobe J, Young G, Sullivan CV. Maternal investment in fish oocytes and eggs: the molecular cargo and its contributions to fertility and early development. Aquaculture. 2017;472(1):37.
Sullivan CV, Chapman RW, Reading BJ, Anderson PE. Transcriptomics of mRNA and egg quality in farmed fish: some recent developments and future directions. Gen Comp Endocrinol. 2015;221:23–30.
Klangnurak W, Fukuyo T, Rezanujjaman MD, Seki M, Sugano S, Suzuki Y, Tokumoto T. Candidate gene identification of ovulation-inducing genes by RNA sequencing with an in vivo assay in zebrafish. PLoS One. 2018;13(5):e0196544.
Aegerter S, Jalabert B, Bobe J. mRNA stockpile and egg quality in rainbow trout (Oncorhynchus mykiss). Fish Physiol Biochem. 2003;28(1–4):317–8.
Aegerter S, Jalabert B, Bobe J. Large scale real-time PCR analysis of mRNA abundance in rainbow trout eggs in relationship with egg quality and post-ovulatory ageing. Mol Reprod Dev. 2005;72(3):377–85.
Aegerter S, Jalabert B, Bobe J. Messenger RNA stockpile of cyclin B, insulin-like growth factor I, insulin-like growth factor II, insulin-like growth factor receptor Ib, and p53 in the rainbow trout oocyte in relation with developmental competence. Mol Reprod Dev. 2004;67(2):127–35.
Bonnet E, Fostier A, Bobe J. Microarray-based analysis of fish egg quality after natural or controlled ovulation. BMC Genomics. 2007;8.
Bonnet E, Montfort J, Esquerre D, Hugot K, Fostier A, Bobe J. Effect of photoperiod manipulation on rainbow trout (Oncorhynchus mykiss) egg quality: a genomic study. Aquaculture. 2007;268(1–4):13–22.
Nagler JJ, Cavileer TD, Stoddard JW, Parsons JE. Maternal mRNA differences in unfertilized rainbow trout (Oncorhynchus mykiss) eggs from batches exhibiting variable embryonic survival. Cybium. 2008;32(2):233.
Ma H, Weber GM, Wei HR, Yao JB. Identification of mitochondrial genome-encoded small RNAs related to egg deterioration caused by postovulatory aging in rainbow trout. Mar Biotechnol. 2016;18(5):584–97.
Ma H, Weber GM, Hostuttler MA, Wei H, Wang L, Yao J. MicroRNA expression profiles from eggs of different qualities associated with post-ovulatory ageing in rainbow trout (Oncorhynchus mykiss). BMC Genomics. 2015;16:201.
Bachvarova RF. A maternal tail of poly(a) - the long and the short of it. Cell. 1992;69(6):895–7.
Pique M, Lopez JM, Foissac S, Guigo R, Mendez R. A combinatorial code for CPE-mediated translational control. Cell. 2008;132(3):434–48.
Villalba A, Coll O, Gebauer F. Cytoplasmic polyadenylation and translational control. Curr Opin Genet Dev. 2011;21(4):452–7.
Cui J, Sartain CV, Pleiss JA, Wolfner MF. Cytoplasmic polyadenylation is a major mRNA regulator during oogenesis and egg activation in Drosophila. Dev Biol. 2013;383(1):121–31.
Slater DW, Slater I, Gillespie D. Post-fertilization synthesis of Polyadenylic acid in sea-urchin embryos. Nature. 1972;240(5380):333.
Rosenthal ET, Tansey TR, Ruderman JV. Sequence-specific Adenylations and Deadenylations accompany changes in the translation of maternal messenger-Rna after fertilization of Spisula oocytes. J Mol Biol. 1983;166(3):309–27.
Paris J, Philippe M. Poly(a) metabolism and polysomal recruitment of maternal mRNAs during early Xenopus development. Dev Biol. 1990;140(1):221–4.
Brevini-Gandolfi TAL, Favetta LA, Mauri L, Luciano AM, Cillo F, Gandolfi F. Changes in poly(a) tail length of maternal transcripts during in vitro maturation of bovine oocytes and their relation with developmental competence. Mol Reprod Dev. 1999;52(4):427–33.
Blower MD, Jambhekar A, Schwarz DS, Toombs JA. Combining different mRNA capture methods to analyze the transcriptome: analysis of the Xenopus laevis transcriptome. PLoS One. 2013;8(10).
Gohin M, Fournier E, Dufort I, Sirard MA. Discovery, identification and sequence analysis of RNAs selected for very short or long poly a tail in immature bovine oocytes. Mol Hum Reprod. 2014;20(2):127–38.
Subtelny AO, Eichhorn SW, Chen GR, Sive H, Bartel DP. Poly(a)-tail profiling reveals an embryonic switch in translational control. Nature. 2014;508(7494):66.
Reyes JM, Chitwood JL, Ross PJ. RNA-Seq profiling of single bovine oocyte transcript abundance and its modulation by cytoplasmic polyadenylation. Mol Reprod Dev. 2015;82(2):103–14.
Meijer HA, Bushell M, Hill K, Gant TW, Willis AE, Jones P, de Moor CH. A novel method for poly(a) fractionation reveals a large population of mRNAs with a short poly(a) tail in mammalian cells. Nucleic Acids Res. 2007;35(19).
Cabada MO, Darnbrough C, Ford PJ, Turner PC. Differential accumulation of two size classes of poly(a) associated with messenger RNA during oogenesis in Xenopus laevis. Dev Biol. 1977;57(2):427–39.
Nagler JJ, Parsons JE, Cloud JG. Single pair mating indicates maternal effects on embryo survival in rainbow trout, Oncorhynchus mykiss. Aquaculture. 2000;184(1–2):177–83.
Jensen JOT. New mechanical shock sensitivity units in support of criteria for protection of salmonid eggs from blasting or seismic disturbance. Can Tech Rep Fish Aquat Sci. 2003(2452):27.
Jensen JOT, Alderdice DF. Comparison of mechanical shock sensitivity of eggs of five Pacific salmon (Oncorhynchus) species and steelhead trout (Salmo gairdneri). Aquaculture. 1989;78:19.
Marlow FL. Maternal Control of Development in Vertebrates. CA: My Mother Made Me Do It! San Rafael; 2010.
Kane DA, Kimmel CB. The zebrafish Midblastula transition. Development. 1993;119(2):447–56.
Zamir E, Kam Z, Yarden A. Transcription-dependent induction of G1 phase during the zebra fish midblastula transition. Mol Cell Biol. 1997;17(2):529–36.
Mathavan S, Lee SGP, Mak A, Miller LD, Murthy KRK, Govindarajan KR, Tong Y, Wu YL, Lam SH, Yang H, et al. Transcriptome analysis of zebrafish embryogenesis using microarrays. PLoS Genet. 2005;1(2):260–76.
Kleppe L, Edvardsen RB, Kuhl H, Malde K, Furmanek T, Drivenes O, Reinhardt R, Taranger GL, Wargelius A. Maternal 3'UTRs: from egg to onset of zygotic transcription in Atlantic cod. BMC Genomics. 2012;13:443.
Hall TE, Smith P, Johnston IA. Stages of embryonic development in the Atlantic cod Gadus morhua. J Morphol. 2004;259(3):255–70.
Kraeussling M, Wagner TU, Schartl M. Highly asynchronous and asymmetric cleavage divisions accompany early transcriptional activity in pre-blastula Medaka embryos. PLoS One. 2011;6(7).
Stoddard JW, Parsons JE, Nagler JJ. Early onset of embryonic mortality in sub-fertile families of rainbow trout (Oncorhynchus mykiss). Reprod Fertil Dev. 2005;17(8):785–90.
Fernandez-Silva P, Enriquez JA, Montoya J. Replication and transcription of mammalian mitochondrial DNA. Exp Physiol. 2003;88(1):41–56.
Temperley RJ, Wydro M, Lightowlers RN, Chrzanowska-Lightowlers ZM. Human mitochondrial mRNAs-like members of all families, similar but different. Bba-Bioenergetics. 2010;1797(6–7):1081–5.
Curanovic D, Cohen M, Singh I, Slagle CE, Leslie CS, Jaffrey SR. Global profiling of stimulus-induced polyadenylation in cells using a poly(a) trap. Nat Chem Biol. 2013;9(11):671.
Richter JD. Cytoplasmic polyadenylation in development and beyond. Microbiol Mol Biol R. 1999;63(2):446–+.
Wai T, Ao A, Zhang XY, Cyr D, Dufort D, Shoubridge EA. The role of mitochondrial DNA copy number in mammalian fertility. Biol Reprod. 2010;83(1):52–62.
Chappel S. The role of mitochondria from mature oocyte to viable blastocyst. Obstet Gynecol Int. 2013;2013:183024.
Artuso L, Romano A, Verri T, Domenichini A, Argenton F, Santorelli FM, Petruzzella V. Mitochondrial DNA metabolism in early development of zebrafish (Danio rerio). Bba-Bioenergetics. 2012;1817(7):1002–11.
St John J. The control of mtDNA replication during differentiation and development. Bba-Gen Subjects. 2014;1840(4):1345–54.
Bai J, Solberg C, Fernandes JMO, Johnston IA. Profiling of maternal and developmental-stage specific mRNA transcripts in Atlantic halibut Hippoglossus hippoglossus. Gene. 2007;386:202–10.
Schatten H, Sun QY, Prather R. The impact of mitochondrial function/dysfunction on IVF and new treatment possibilities for infertility. Reprod Biol Endocrin. 2014;12:1–11.
Shoubridge EA, Wai T. Mitochondrial DNA and the mammalian oocyte. Mitochondrion Germline Early Dev. 2007;77:87–111.
Sanchez-Velar N, Udofia EB, Yu Z, Zapp ML. hRIP, a cellular cofactor for rev function, promotes release of HIV RNAs from the perinuclear region. Genes Dev. 2004;18(1):23–34.
Catrina IE, Bayer LV, Yanez G, McLaughlin JM, Malaczek K, Bagaeva E, Marras SAE, Bratu DP. The temporally controlled expression of Drongo, the fruit fly homolog of AGFG1, is achieved in female germline cells via P-bodies and its localization requires functional Rab11. RNA Biol. 2016;13(11):1117–32.
Gentil BJ, Delphin C, Mbele GO, Deloulme JC, Ferro M, Garin J, Baudier J. The giant protein AHNAK is a specific target for the calcium- and zinc-binding S100B protein - potential implications for Ca2+ homeostasis regulation by S100B. J Biol Chem. 2001;276(26):23253–61.
Doidge R, Mittal S, Aslam A, Winkler GS. The anti-proliferative activity of BTG/TOB proteins is mediated via the Caf1a (CNOT7) and Caf1b (CNOT8) Deadenylase subunits of the Ccr4-not complex. PLoS One. 2012;7(12).
Matsuda S, KawamuraTsuzuku J, Ohsugi M, Yoshida M, Emi M, Nakamura Y, Onda M, Yoshida Y, Nishiyama A, Yamamoto T. Tob, a novel protein that interacts with p185(erbB2), is associated with antiproliferative activity. Oncogene. 1996;12(4):705–13.
Mauxion F, Chen CYA, Seraphin B, Shyu AB. BTG/TOB factors impact deadenylases. Trends Biochem Sci. 2009;34(12):640–7.
Migaud H, Bell G, Cabrita E, McAndrew B, Davie A, Bobe J, Herraez MP, Carrillo M. Gamete quality and broodstock management in temperate fish. Rev Aquacult. 2013;5:S194–223.
Chapman RW, Reading BJ, Sullivan CV. Ovary transcriptome profiling via artificial intelligence reveals a transcriptomic fingerprint predicting egg quality in striped bass, Morone saxatilis. PLoS One. 2014;9(5).
Zardoya R, GarridoPertierra A, Bautista JM. The complete nucleotide sequence of the mitochondrial DNA genome of the rainbow trout, Oncorhynchus mykiss. J Mol Evol. 1995;41(6):942–51.
Slomovic S, Laufer D, Geiger D, Schuster G. Polyadenylation and degradation of human mitochondrial RNA: the prokaryotic past leaves its mark. Mol Cell Biol. 2005;25(15):6427–35.
Baserga SJ, Linnenbach AJ, Malcolm S, Ghosh P, Malcolm ADB, Takeshita K, Forget BG, Benz EJ. Polyadenylation of a human mitochondrial ribosomal-Rna transcript detected by molecular-cloning. Gene. 1985;35(3):305–12.
Rorbach J, Minczuk M. The post-transcriptional life of mammalian mitochondrial RNA. Biochem J. 2012;444:357–73.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12).
Dutta A, Sinha DK. Zebrafish lipid droplets regulate embryonic ATP homeostasis to power early development. Open Biol. 2017;7(7).
Wendling NC, Bencic DC, Nagler JJ, Cloud JG, Ingermann RL. Adenosine triphosphate levels in steelhead (Oncorhynchus mykiss) eggs: an examination of turnover, localization and role. Comp Biochem Phys A. 2004;137(4):739–48.
Cagnone GLM, Tsai TS, Makanji Y, Matthews P, Gould J, Bonkowski MS, Elgass KD, Wong ASA, Wu LE, McKenzie M, et al. Restoration of normal embryogenesis by mitochondrial supplementation in pig oocytes exhibiting mitochondrial DNA deficiency. Sci Rep-Uk. 2016;6.
Dumollard R, Duchen M, Sardet C. Calcium signals and mitochondria at fertilisation. Semin Cell Dev Biol. 2006;17(2):314–23.
Ge HS, Tollner TL, Hu Z, Dai MM, Li XH, Guan HQ, Shan D, Zhang XJ, Lv JQ, Huang CJ, et al. The importance of mitochondrial metabolic activity and mitochondrial DNA replication during oocyte maturation in vitro on oocyte quality and subsequent embryo developmental competence. Mol Reprod Dev. 2012;79(6):392–401.
Huang DW, Sherman BT, Tan Q, Collins JR, Alvord WG, Roayaei J, Stephens R, Baseler MW, Lane HC, Lempicki RA. The DAVID gene functional classification tool: a novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol. 2007;8(9).
Ma H, Gao G, Weber GM. Use of DAVID algorithms for clustering custom annotated gene lists in a non-model organism, rainbow trout. BMC Res Notes. 2018;11(63):1–6.
Hershberger WK, Hostuttler MA. Variation in time to first cleavage in rainbow trout Oncorhynchus mykiss embryos: a major factor in induction of tetraploids. J World Aquacult Soc. 2005;36(1):96–102.
Velsen FPJ. Embryonic development in eggs of sockeye salmon, Oncorhynchus nerka. Can Spec Publ Fish Aquat Sci. 1980;49:1–19.
Berthelot C, Brunet F, Chalopin D, Juanchich A, Bernard M, Noel B, Bento P, Da Silva C, Labadie K, Alberti A, et al. The rainbow trout genome provides novel insights into evolution after whole-genome duplication in vertebrates. Nat Commun. 2014;5.
Kodama S, Yamada H, Annab L, Barrett JC. Elevated expression of mitochondrial cytochrome-B and Nadh dehydrogenase Subunit-4/4l genes in senescent human-cells. Exp Cell Res. 1995;219(1):82–6.
Abernathy J, Overturf K. Comparison of ribosomal RNA removal methods for transcriptome sequencing workflows in teleost fish. Anim Biotechnol. 2016;27(1):60–5.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–U354.
Conesa A, Gotz S. Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008;2008:619832.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Zhang HS, Maguire D, Swarts S, Sun WM, Yang SM, Wang W, Liu CM, Zhang M, Zhang D, Zhang L, et al. Replication of murine mitochondrial DNA following irradiation. Adv Exp Med Biol. 2009;645:43–8.
We thank Dr. Guangtu Gao for assistance in providing ribosomal reference sequences and help in data analysis and discussion; and Jill Birkett for assistance with embryo assessments. The use of trade, firm, or corporation names in this publication is for the information and convenience of the reader. Such use does not constitute an official endorsement or approval by the USDA or the ARS of any product or service to the exclusion of others that may be suitable. USDA is an equal opportunity provider and employer.
This work was supported by USDA-ARS CRIS Project 8082–31000-012. The funding bodies had no role in the design of the experiment, interpretation of the data, or writing of the manuscript.
Availability of data and materials
The sequences used in this study were deposited into NCBI Sequence Read Archive under accession number SPR108797 (https://www.ncbi.nlm.nih.gov/sra?term=SRP108797).
All animal experiments were conducted under a protocol approved by the USDA/ARS National Center for Cool and Cold Water Aquaculture Institutional Animal Care and Use Committee (protocol #50).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Differentially expressed transcripts among low-, medium- and high-quality eggs identified by DESeq2. Table S2. Mean and log2FC of normalized mitochondrial DET reads. Table S3. Quality analyses for RNAs treated with DNase and submitted for RNA sequencing. Table S4. Overview of RNA-seq read alignments from the libraries constructed with defective rRNA removal kit. Table S5. Additional gene transcripts used as reference. Table S6. Primers used in real-time quantitative PCR. Table S7. Gene name abbreviations. (XLSX 506 kb)
Figure S1. Electrophoresis of egg RNAs isolated from different families selected for RNA sequencing. The families labeled in red are from the low-quality group; the families labeled in green are from the medium quality group; the families labeled in black are from the high-quality group. About 400 ng/sample of RNA was loaded to each well. Family 10 was not used for RNA sequencing. (PPTX 122 kb)