Identification of differentially expressed ovarian genes during primary and early secondary oocyte growth in coho salmon, Oncorhynchus kisutch
Reproductive Biology and Endocrinologyvolume 6, Article number: 2 (2008)
The aim of this study was to identify differentially expressed ovarian genes during primary and early secondary oocyte growth in coho salmon, a semelparous teleost that exhibits synchronous follicle development.
Reciprocal suppression subtractive hybridization (SSH) libraries were generated from ovaries with perinucleolus (P) or cortical alveolus (CA) stage follicles and selected genes were assessed with quantitative PCR (qPCR). An assessment of changes in RNA composition during oocyte growth and its relationship to transcript levels was also conducted.
SSH revealed several differentially expressed genes during early oogenesis, some which will not likely be utilized until 1–3 years later in salmon. Zona pellucida glycoprotein (zp) genes, vitellogenin receptor (vldlr) isoforms, cathepsin B (ctsba), cyclin E (ccne), a DnaJ transcript (dnaja2), and a ferritin subunit (fth3) were significantly elevated at the P stage, while a C-type lectin, retinol dehydrogenase (rdh1), and a coatomer protein subunit (cope) were upregulated at the CA stage. Putative follicle cell transcripts such as anti-Müllerian hormone (amh), lipoprotein lipase (lpl), apolipoprotein E (apoe), gonadal soma-derived growth factor (gsdf) and follicle-stimulating hormone receptor (fshr) also increased significantly at the CA stage. The analysis of RNA composition during oocyte growth showed that the total RNA yield and proportion of messenger RNA relative to non-polyadenylated RNAs declined as oogenesis progressed. This influenced apparent transcript levels depending on the type of RNA template used and normalization method.
In coho salmon, which exhibit a dramatic change in oocyte size and RNA composition during oogenesis, use of messenger RNA as template and normalization of qPCR data to a housekeeping gene, ef1a, yielded results that best reflected transcript abundance within the ovarian follicle. Synthesis of zp transcripts and proteins involved in yolk incorporation and processing occurred during primary growth, while increased expression of a CA component and genes related to lipid incorporation occurred concomitant with the appearance of CA, but prior to lipid accumulation. Significant increases in transcripts for fshr, gsdf, and amh at the CA stage suggest a role of FSH and TGFβ peptides in previtellogenic oocyte growth and puberty onset in female salmon.
Oocyte growth is a period of intense RNA synthesis, replication and redistribution of cytoplasmic organelles, and nutrient incorporation in oviparous vertebrates. In teleost fish, this period may encompass a significant portion of the lifespan, lasting well over a decade in some species. Despite this, research on fish oogenesis has primarily focused on vitellogenesis, final maturation and ovulation, while stages of primary and early secondary oocyte growth remain largely unexplored [1, 2]. For example, it remains unclear what endocrine and/or intraovarian factors regulate oocyte growth and how this period may influence timing of puberty, fecundity, egg quality, and early embryogenesis.
Similar to primordial follicle development in mammals, primary oocyte growth in fish begins with the onset of meiosis and subsequent meiotic arrest in the diplotene stage of the first prophase. The oocytes are then completely enveloped by a monolayer of presumptive granulosa cells and a thin theca cell layer and epithelial sheath are added to the surface, forming the basic follicle structure [2, 3]. As the follicle develops, the nucleus of the oocyte increases in size and numerous ribosome-producing nucleoli appear around its periphery ("perinucleolus" stage). Intense RNA synthesis occurs over this period and much of the RNA present in the fully grown oocyte is thought to be synthesized at this time [4–6]. During primary growth alone in fish, the oocyte volume may increase as much as 1,000- to 5,000-fold .
Initiation of secondary growth is signified by the appearance and accumulation of cortical alveoli (formerly yolk vesicles). These endogenously synthesized secretory vesicles, analogous to cortical granules in invertebrates and other vertebrates, are derived from Golgi bodies and play important roles in the fertilization response and early embryogenesis . Upon fertilization, cortical alveoli fuse with the oocyte membrane and discharge their glycoprotein contents into the perivitelline space to prevent polyspermy and entry of microbes or pathogens. Cortical alveoli increase in number during early secondary growth, initially forming a ring around the periphery of the oocyte and then accumulating inward to the nucleus. In most fishes, a brief period of oocyte lipid deposition (lipid droplet stage) occurs late in the cortical alveolus stage and prior to significant yolk incorporation. Vitellogenesis (yolk incorporation) marks the final phase of secondary growth, during which dramatic follicle growth occurs as the oocyte sequesters vitellogenin, a hepatically derived yolk protein precursor, from the bloodstream .
Through recent large-scale genomic studies mainly conducted on zebrafish, salmonids, and Fugu pufferfish [8–13], a number of ovarian genes have been sequenced and catalogued in databases making it possible to identify many fish mRNAs, profile their expression, and determine their function(s). Genes involved in sex differentiation and early gametogenesis [14–16], and final oocyte maturation  have received considerable attention, while other studies have focused on specific gene families such as TGFβ superfamily members , zona pellucida glycoproteins , and vitellogenin receptor (very low density lipoprotein receptor, vldlr) . Through these studies highly expressed ovarian genes, such as zona pellucida glycoprotein (zp) genes and egg lectins have been revealed. However, relatively few ovarian genes have been profiled in fish and little is known about temporal gene expression during oocyte growth.
The objective of this study was to identify differentially expressed genes during primary and early secondary oocyte growth with the ultimate goal of establishing what regulates these stages of oogenesis and identifying the intraovarian factors that drive (or block) puberty onset. Coho salmon, Oncorhynchus kisutch, was selected as a model because it is a semelparous species (spawns only once in its life and then dies) that exhibits synchronous follicle development. This unique reproductive life history allows for stage specific transcript analysis of a homogeneous clutch of follicles, which is not possible in iteroparous species like rainbow trout (O. mykiss) and Atlantic salmon (Salmo salar). Salmonids are also one of the best studied groups of fish and the large available repository of trout and Atlantic salmon ESTs facilitates identification and functional annotation of coho salmon transcripts. In this study, suppression subtractive hybridization (SSH) was conducted with ovaries containing follicles in primary or early secondary growth and putatively regulated genes identified by SSH were screened and then validated with real-time quantitative PCR (qPCR). Quantitative PCRs were developed for 17 genes identified by SSH and additional assays were developed for candidate genes of interest including known follicle cell transcripts, such as follicle-stimulating hormone receptor (fshr) and anti-Müllerian hormone (amh). Lastly, a thorough assessment of changes in the mRNA/total RNA ratio during oocyte growth was conducted to determine how gene expression results are influenced by rapid growth of follicles during oogenesis. This is germane to ovarian transcript analyses in other species and potentially other rapidly growing and/or differentiating tissues.
Animals and sampling
Coho salmon (2004 and 2003 brood) were reared at the Northwest Fisheries Science Center (Seattle, WA) in recirculated fresh water under a simulated natural photoperiod and fed a standard ration of a commercial diet. These salmon typically spawn in December at 3 years of age and the same cohorts (2004 brood = cohort 1, 2003 brood = cohort 2) were used for all experiments. In October 2005, female salmon (N = 43 fish of cohort 1; N = 20 fish of cohort 2) were euthanized and their ovaries removed and weighed. At this time, cohort 1 fish were 0+ age (10 months old), 72–103 mm fork length (FL) and 4.5–13.3 g body mass with an ovary mass of 0.029 ± 0.001 (mean ± SEM) and gonadosomatic index (GSI) of 0.32 ± 0.01. Cohort 2 fish were 1+ age (22 months old), 194–235 mm FL and 93.5–157.8 g body mass with an ovary mass of 0.400 ± 0.021 g and GSI of 0.33 ± 0.01. A piece of ovary from each fish was fixed in Bouins for paraffin histology , and the remaining tissue was snap frozen for later RNA isolation. Histology revealed that in October 2005, ovaries of cohort 1 possessed primary growth follicles at the perinucleolus (P) stage containing minimal Balbiani material and no cortical alveoli (Fig. 1). Ovaries of cohort 2 contained follicles that were early to mid-cortical alveolus (CA) stage with cortical alveoli filling greater than 50% of the ooplasm. Samples were selected for SSH that represented primary growth (P stage, cohort 1) and early secondary oocyte growth (mid-CA stage, cohort 2). Salmon used for qPCR validation of the SSH results (N = 10 fish/stage) possessed ovaries that were in the same stages as those used for SSH. For the template RNA assessment (see below), fish were sampled as described above at a later time point in August 2006 (N = 5/cohort). At this time, fish from cohort 1 were 1+ age (20 months old) and follicles were generally late P stage with most oocytes showing very few cortical alveoli in the periphery of the ooplasm, while cohort 2 fish were 2+ age (30 months old) and oocytes were in the yolk granule stage.
Fish were reared and handled according to the policies and guidelines of the University of Washington Institutional Animal Care and Use Committee (IACUC Protocol #2313-09).
RNA isolation and cDNA synthesis
Ovarian total RNA was isolated with Tri-Reagent (Molecular Research Center, Cincinnati, OH). For SSH, total RNA pools were generated from 18 P stage and 3 CA stage ovaries (whole ovaries minus the part for fixation) and mRNA was further isolated via the PolyATract mRNA Isolation System (Promega, Madison, WI). For individual samples for qPCR, the MicroPoly(A)Purist kit (Ambion, Austin, TX) was used. Samples for qPCR were reverse transcribed with SuperScript II (Invitrogen, Carlsbad, CA) and either 100 ng mRNA or 1000 ng total RNA in 20 μl reactions.
Subtractive hybridization and analysis
SSH was performed with 1 μg of mRNA per stage using the PCR-Select cDNA Subtraction and Advantage cDNA PCR kits (Clontech, Palo Alto, CA). Subtracted cDNA samples were cloned in TOPO pCR 2.1 (Invitrogen) and plated on LB agar containing 100 μg/ml kanamycin. A total of 288 clones per library were randomly selected for sequencing with Big Dye Terminator (Applied Biosystems (ABI), Foster City, CA) on an ABI 3730 sequencer using the M13 reverse primer. Sequence chromatogram files were trimmed for quality (phred), vector-screened (cross match), and analyzed in the following manner to determine transcript identities: sequences ≤ 30 bp were removed and the remaining sequences were analyzed for redundancy with CAP3 , and analyzed locally using blastx against the NCBI nonredundant (nr) protein database, blastn against the NCBI nucleotide (nt) database, blastn against the EST database (est_others), and blastx against the Gene Ontology (GO) database . The e-value cutoff was 10-5 for blast searches. Nomenclature of the Zebrafish Information Network (ZFIN) was used where possible, except for the zp genes for which we followed the recommended nomenclature of Spargo and Hope .
Approximately 60 PCR primer sets were designed from sequences obtained by SSH with MacVector software (Accelrys, San Diego, CA) and undigested SSH cDNA was used as template to screen for differentially expressed genes with semi-quantitative PCR as described by Goetz et al. . Genes that appeared to be differentially expressed in the initial screen were measured with real-time qPCR using primers listed in Table 1. Assays were run on an ABI 7700 Sequence Detector in 96-well plates using the standard cycling conditions: 50°C for 2 min, 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min. Reactions consisted of 1× Power SYBR Green PCR master mix (ABI), 150 nM of each gene-specific forward and reverse primer, and 0.01 ng cDNA template (based on the mRNA loaded into the RT reactions). For total RNA samples used in the template RNA assessment, 1 ng cDNA was loaded per reaction.
Triplicate standard curve samples generated from a serial dilution of pooled RNA from ovaries of 6 previtellogenic coho salmon were included in each plate. The standard curve (log input cDNA template versus cycle threshold) was linear for cDNA template ranging from 0.001–1 ng. Results were analyzed using the relative standard curve method (for details see the ABI Prism 7700 Sequence Detection System User Bulletin #2, P/N 4303859). Experimental samples were run in duplicate for the housekeeping gene (see below) and not replicated for the target genes. No amplification controls (NAC) that lacked reverse transcriptase and a no template control (NTC) that lacked template altogether were evaluated in each plate to confirm the absence of genomic DNA in the RNA preparations and the absence of PCR carryover contamination, respectively. Negative control samples showed either no detectability or negligible values (>10 Ct separation from samples). Melt curve analysis was included for each target gene to ensure that a single product was amplified. In cases where more than one product was generated the PCR was redesigned. In addition, a qPCR product from each plate was directly sequenced to verify that the target was successfully amplified. Target gene results were normalized to elongation factor-1 alpha (ef1a), which served as a housekeeping gene and has been used in previous gonadal studies [e.g., [15, 26]]. Ef1a transcript levels were not different across stages when mRNA was used as template (P = 0.5). The lowest normalized value for each gene was arbitrarily set to 1 to enhance data presentation.
Template RNA assessment
The RNA composition and size of ovarian follicles change dramatically during oocyte growth, both of which can influence transcript levels depending on the method of normalization. Therefore, we compared transcript levels when measured from total RNA versus mRNA as template, correcting for follicle number and RNA recovery. Follicles were counted prior to RNA isolation and the expression of several genes was measured in total RNA and mRNA preparations from the same individual ovary samples. Briefly, fish were euthanized and body mass, FL, and ovary mass recorded. Fragments weighing ~40 mg taken from the middle region of the ovary were snap frozen for RNA analysis and an additional piece was fixed in Bouins for histology. Triplicate fragments of fresh ovary were weighed and the follicles counted under a dissecting microscope to determine the average follicle mass for each fish. Fecundity and RNA yield per follicle estimates were based on average follicle mass data. Gene expression was assessed in total RNA and mRNA preparations for a housekeeping gene, ef1a; a downregulated gene, ferritin H-3 (fth3); and an upregulated gene, fshr.
Gene expression results were compared by unpaired t-tests and Welch's correction was applied in cases where variances were unequal across stages (Prism 4, GraphPad Software, San Diego, CA). The minimum level of statistical significance was set at P < 0.05.
General library statistics
The SSH library putatively enriched in cDNAs/transcripts that were more abundant at the P stage yielded 275 sequences (>30 bp) with an average size of 420 bp that clustered to form 55 contigs (of two or greater sequences with >90% identity) and 155 singletons. Approximately 80% of the sequences were identified as similar to annotated genes, while 15% had similarity to unannotated genes (entries of unknown function, e.g. hypothetical proteins) and 5% were novel (no significant identity to GenBank sequences). The SSH library putatively enriched in genes upregulated at the CA stage yielded 267 sequences with an average size of 401 bp, which clustered to form 47 contigs and 103 singletons. Approximately 71% of the sequences were similar to annotated genes, 16% were similar to unannotated genes, and 13% were novel. The two predominant putative genes in the P and CA libraries with a frequency of ~10% were fth3 and serum lectin isoform 2 (lcal), respectively.
Expressed sequence tags (ESTs) were submitted to the NCBI database [GenBank: EX152144–EX152685]. Genes further assessed with qPCR are listed in Table 2.
Gene expression analysis
Zona pellucida genes, including an EST similar to factor in the germline alpha (figla), were more abundant in P stage follicles compared to CA stage follicles (Fig. 2). All genes falling into this cascade (figla, zpx, zpc, eeg) showed a similar pattern, but the fold-difference across stages was most pronounced for figla and zpx where mRNA levels were elevated ~2-fold at the P stage. Several genes associated with lipoprotein uptake and yolk processing exhibited this pattern, including vldlr, somatic lipoprotein receptor (vldlro), and cathepsin B (ctsba) (Fig. 3). Cathepsin D (ctsd) and cathepsin Z (ctsz), however, were not differentially expressed across stages when assessed with qPCR (Fig. 3). Lipoprotein lipase (lpl) and apolipoprotein E (apoe) mRNA levels were significantly higher during the CA stage.
Genes encoding cortical alveoli components had variable expression patterns (Fig. 4). Lcal was significantly elevated in CA stage samples, whereas rhamnose binding lectin STL3 (lrham) and alveolin (alv) were not statistically different across stages.
A number of other genes not falling into the above cascades were differentially expressed across stages (Fig. 5). Fth3, DnaJ subfamily A member 2 (dnaja2), and cyclin E (ccne) transcripts were significantly elevated at the P stage, while retinol dehydrogenase 1 (rdh1) and coatomer protein epsilon subunit (cope), showed the opposite pattern and were higher at the CA stage. Anterior pharynx defective 1b (aph1b) was not differentially expressed based on qPCR results. Finally, three transcripts likely originating from the follicle cells, fshr, amh, and gonadal soma-derived growth factor (gsdf) were dramatically upregulated at the CA stage.
Template RNA assessment
The mean quantity of total RNA isolated per gram of ovary for cohorts 1 and 2 were 4750 and 1509 μg, respectively, and thus the total RNA yield was about 3 times greater per g ovary for cohort 1 (Table 3). However, because there were 40 times more follicles per tissue mass for cohort 1, the total RNA yielded per follicle was lower relative to cohort 2. The average mRNA yielded per μg total RNA was also greater for cohort 1 relative to cohort 2. Again, due to the difference in follicle number the mRNA yielded per follicle was actually lower for cohort 1. Based on mRNA yields, total RNA from cohort 1 and 2 consisted of ~12% and 2% mRNA, respectively.
Ef1a transcript levels were similar across stages based on template mRNA (0 to 1.2-fold change), but decreased (1.5 to 4.5-fold) in more advanced ovaries based on total RNA template (Table 3). This decrease was most pronounced when total RNA based data were expressed on a per follicle basis. Other housekeeping genes, such as acidic ribosomal protein, showed results similar to ef1a (data not shown). Fth3 declined in more advanced ovaries with all means of analysis, but this decline was most pronounced with total RNA template. Fshr showed a 6.9 to 7.9-fold increase with mRNA template, while total RNA based samples showed either no increase or a diminished increase (0 to 4.3-fold increase). Results were most skewed in total RNA preparations when presented on a per follicle basis and most closely resembled mRNA preparations when normalized to ef1a.
The goal of this study was to identify differentially expressed ovarian genes during the little studied periods of primary and early secondary oocyte growth in coho salmon. A number of differentially expressed genes were successfully identified through SSH and qPCR that are known to play important roles during oogenesis, such as zona pellucida development, sequestration and processing of lipoproteins, cell cycle control, and the fertilization response. Interestingly, some genes involved in vitellogenesis, fertilization, and embryogenesis were more highly expressed during primary growth than early secondary growth. This pattern is intriguing because proteins encoded by many of these genes are not thought to be utilized until 1–3 years later in this species. Through work mainly conducted in Xenopus it has been well documented that maternal RNAs are deposited in the oocyte early in oogenesis and then stored as messenger ribonucleoprotein particles (mRNPs), which contain proteins such as Y-box proteins and DEAD-box RNA helicases that mask RNA from the translational apparatus [5, 27, 28]. In the present study, the ccne transcript was upregulated during primary growth and would likely represent a maternal RNA based on its role in early embryogenesis (see below). However, other transcripts, such as vldlr and cathepsins, associated with later aspects of oogenesis itself were also highly expressed during primary growth. Although protein data are necessary to resolve this issue, one possible explanation for the early transcription of such genes is that some oogenesis related mRNAs are subject to masking during early oogenesis as described for traditional maternal RNAs.
Other transcripts such as lpl, apoe, lcal, rdh1, amh, gsdf, and fshr increased significantly during early secondary growth when CA are abundant in the oocytes, but lipid droplets are not yet evident. Some of these transcripts were likely derived from the follicle cells and showed a dramatic increase across stages (2 to 6 fold). These data together with results of previous studies discussed in detail below, suggest an increased role of the follicle cells during early secondary growth in coho salmon and potential regulation of this transition by FSH and TGFβ family peptides.
Based on GO annotation and examination of the literature, the majority of genes revealed by SSH were of oocyte origin. The SSH method employed enriches for rare transcripts through PCR steps; however it appears that copious mRNA from the oocytes overshadowed differentially expressed genes originating from the follicle/interstitial cells. This idea is supported by the observed differences in abundance of follicle cell transcripts, such as amh and gsdf, which were candidate genes obtained from a follicle/interstitial cell enriched library . As first noted by Goetz and colleagues  these findings suggest that it is necessary to enrich for these cell layers prior to SSH to increase the likelihood of revealing rare transcripts of the granulosa or theca/interstitial cells.
Analyzing transcript levels across stages of oogenesis
One problem with studying the ontogeny of ovarian gene expression, especially in oviparous vertebrates with large eggs, is the dramatic change in oocyte size and RNA composition that occurs during oogenesis. Changes in the size and number of follicles per tissue mass and instability of housekeeping genes complicate across stage comparisons making it difficult to interpret transcript abundance data in a biologically relevant way [e.g., [21, 31–34]]. To gain insight into this problem, we quantified RNA yields per tissue mass and per follicle, and transcript abundance in total RNA and mRNA preparations. As oogenesis progressed, total RNA yielded per tissue mass and the proportion of mRNA relative to non-polyadenylated RNA declined (Table 3). Furthermore, the yield of RNA per follicle was higher in more advanced follicles but the amount of mRNA did not increase proportional to follicle size. This change in RNA composition had a significant effect on apparent transcript levels for a variety of genes expressed in the oocyte and follicle cells depending on whether mRNA or total RNA was used as template for cDNA synthesis. Housekeeping genes, such as ef1a, were similarly expressed across stages when mRNA was used as the template, irrespective of whether data were calculated per unit template RNA or per follicle. Our results indicate that if total RNA is used as a template without normalization to a housekeeping gene one can get misleading results, such as an apparent decline in transcripts that are actually stably expressed, an accentuated decline in downregulated genes, or no change in transcripts for upregulated genes. Results were most skewed in total RNA preparations when presented on a per follicle basis and most closely resembled mRNA based results when normalized to ef1a. In summary, use of mRNA as template and normalization of qPCR data to ef1a generated results that best reflected transcript abundance within the follicle.
Zona pellucida glycoprotein genes
The zona pellucida is the acellular membrane that not only encloses the oocyte, but is also critical to optimal oocyte growth and preventing polyspermy after fertilization. The zona pellucida in vertebrates consists of highly sulfated zona pellucida glycoproteins (ZPs). Four ZP subfamilies, ZPA, ZPX, ZPB and ZPC, have been described and are each found in fish except ZPA [19, 24]. Several coho salmon zp genes were revealed by SSH and all exhibited significantly elevated mRNA levels during primary growth relative to early secondary growth (Fig. 2). Consistent with this, other studies have shown that zp transcripts are abundant in oocytes during early oogenesis and make up a significant portion of the ovarian transcriptome [7–9, 13, 31].
Few studies in fish have focused on the zp-related transcription factor, figla. In mammals, however, figla is required for primordial follicle development and transcription of the zp genes [6, 35]. Like the zp genes, figla is highly expressed in primary oocytes and has been localized to the ooplasm in medaka fish, Oryzias latipes . In the present study, mRNA levels for both the zp genes and figla EST were higher during primary growth. The similar profile of figla and the zp genes in this study along with work in other species [35, 36] suggests that transcription of this family of genes is highly coordinated.
Lipoprotein uptake and processing
Several transcripts associated with vitellogenesis, such as vldlr and vldlro, were highly expressed during primary growth and declined significantly by early secondary growth (Fig. 3). Vldlr has been widely studied and is responsible for uptake of hepatically-derived vitellogenin by the oocyte . On the other hand, very little is known about vldlro, which contains an O-linked sugar domain, is expressed in the ovary and somatic tissues of some fishes [33, 37, 38], and is thought to mediate uptake of lipoproteins other than vitellogenin . Findings of the present study are consistent with work in rainbow trout which showed that transcription of vldlr began shortly after female differentiation  and declined by early vitellogenesis, becoming nearly undetectable by mid-vitellogenesis . Given that vldlr mRNA levels are low during much of the period of active uptake of vitellogenin, it is widely believed that the Vldlr protein is recycled during vitellogenesis [20, 37, 40].
Two genes, lpl and apoe, associated with lipid or lipoprotein uptake by oocytes increased 2–3 fold during early secondary growth. Lpl cleaves fatty acids from plasma lipoproteins and thus facilitates lipid transport across biological membranes, while Apoe is a lipid binding protein that mediates recognition and internalization of plasma lipoproteins by cell surface receptors. The timing of this increase is interesting considering that the CA stage just precedes the appearance of lipid inclusions within the salmon oocyte. Few studies have focused on Lpl and Apoe in the teleost ovary. However, studies in rainbow trout [41, 42] demonstrated that lpl mRNA levels and Lpl activity increased steadily during vitellogenesis and peaked at late vitellogenesis. In sea bass (D. labrax), lpl mRNA was localized to the follicle cells and high lpl mRNA levels and Lpl activity coincided with the appearance of oocyte lipid inclusions . Thus mounting evidence suggests that Lpl may be involved in lipid uptake associated with secondary oocyte growth in fish. To our knowledge, apoe mRNA has been measured in the ovary in only one other study in fish , which demonstrated expression during very early oogenesis. Together these data indicate there is increased expression of genes associated with lipid transport prior to significant lipid incorporation.
Cathepsins are lysosomal enzymes that in the oocyte are responsible for proteolytically cleaving vitellogenin into its constituent yolk proteins and play a role in oocyte reabsorption during atresia [2, 44, 45]. Coho salmon cathepsins revealed by SSH showed different patterns of expression during oocyte growth. The ctsb transcript was more abundant during primary growth relative to early secondary growth. Both ctsd and ctsz mRNA levels were not different across stages and would be considered false positives of SSH, which are commonly encountered with this technique . Ctsd is considered the major protease responsible for cleaving vitellogenin into yolk proteins . However, in the barfin flounder, Ctsb was implicated in this role . In the teleost, Fundulus heteroclitus, ctsz expression was relatively stable during oocyte growth and maturation , and in rainbow trout ctsz was upregulated during maturation  and correlated with egg quality . Data generally suggest that ovarian cathepsins may be regulated post-transcriptionally and thus transcript abundance may not correlate well with cathepsin enzymatic activity [42, 44, 50].
Cortical alveoli components
Studies in Xenopus have shown that 70% of the proteins contained in cortical granules are lectins . In fish oocytes, lectins have a number of biological functions including block to polyspermy at fertilization and defense against pathogens . In the present study, lcal, a Ca2+-dependent or C-type lectin, was more highly expressed at the CA stage and was the predominant transcript from the CA stage SSH library, which is not surprising given the abundance of CA evident in the ovarian histology. Indeed, C-type lectins are often highly represented in ovary transcriptomes in fish [9, 16, 52]. In contrast, however, two other CA components, lrham and alv, were not differentially expressed across stages. Lrham is an oocyte lectin implicated in block of polyspermy that was localized to CA of rainbow trout . Alv is a metalloproteinase first identified in medaka that is released by CA after fertilization to induce zona pellucida (chorion) hardening . Interestingly, when examined as a group, the CA component genes identified in the present study were not transcribed coordinately during oocyte growth as shown in rainbow trout .
Other ovarian transcripts
A subunit of ferritin heavy polypeptide, fth3, was the predominant transcript in the P stage library and qPCR verified that mRNA levels were elevated during primary growth. Ferritins are important because they store and transport iron atoms . Since iron is a critical constituent of metalloproteins, such as enzymes and oxygen carriers, it is essential to all organisms. Free iron, however, can be highly toxic to cells and thus ferritin protects cells from the damaging effects of iron, but makes it readily available. Studies have documented high levels of ferritin subunit mRNAs in the ovary [13, 49] but little is known about their specific role.
Like ferritin, dnaja2 and ccne showed significantly higher transcript levels during primary growth. Dnaja2 encodes a chaperone protein associated with unfolded protein and heat shock protein binding, while cyclins are positive cell cycle regulators that appear to be profuse in the fish ovary [30, 57]. Ccne is transcribed and stored during oocyte growth in goldfish and is thought to be important to the first embryonic cell cycles . The high expression of ccne during primary growth in this study is consistent with its early transcription in goldfish and suggests it is a classic maternal mRNA, as shown for some other cyclins .
Rdh1 and cope showed an opposite profile with transcript levels higher during early secondary growth. Retinol dehydrogenases are involved in the synthesis of retinoic acid, the active form of vitamin A, which regulates cellular growth and differentiation, embryogenesis, and reproduction in vertebrates . Because of the diverse functions of retinol dehydrogenases, it is unclear what role Rdh1 may play during previtellogenic growth. However, based on the recent detection of other retinoid-related transcripts, such as retinol dehydrogenase type II and retinol binding protein in the trout ovary , several players in this cascade are present during oogenesis and likely play an important role during oogenesis and/or early embryogenesis.
The cope transcript encodes the epsilon subunit of a coatomer complex protein. One well characterized coatomer protein is clathrin, which is involved in receptor-mediated endocytosis and the transport of proteins from the Golgi network . At this point it is not possible to determine what process the coatomer complex gene identified here is associated with, but perhaps it is involved in uptake of vitellogenin by receptor-mediated endocytosis, protein trafficking, and/or the immune response.
Transcripts for fshr, amh, gsdf, and lpl, which are all likely produced in follicle cells, exhibited the largest change in abundance across the P and CA stages with levels increasing 4–5 fold. The increase in fshr during this transition from primary to secondary growth together with previous data showing a progressive increase in plasma FSH from the P to lipid droplet stage and effects of FSH on ovarian steroidogenesis [21, 61–63], suggest FSH plays an important role in the endocrine control of this phase of oogenesis in coho salmon. It is not known, however, whether the expression of follicle cell transcripts such as those identified in this study is regulated by FSH and/or other endocrine or paracrine factors.
Transcript levels for two members of the TGFβ superfamily, amh and gsdf, also increased 4–5 fold from the P to CA stage. In female mammals, AMH is produced in the granulosa cells, increases at puberty onset reaching peak levels in small antral follicles, diminishes in later stages, and is no longer detectable during the FSH-dependent final stages of follicle growth or in atretic follicles . Recently, the structure of Amh has been characterized in fish [65, 66] and studies have primarily focused on its expression during sex differentiation [39, 65]. During oogenesis in zebrafish, amh has been localized to granulosa cells where transcript levels peak at the CA stage and progressively decline at onset of yolk incorporation, reaching non-detectable levels by late vitellogenesis. Gsdf is a recently identified gonad-specific cytokine that appears to exist only in teleosts . In rainbow trout, gsdf mRNA was localized to somatic cells of the genital ridge during embryogenesis and Sertoli and granulosa cells during gametogenesis . Gsdf plays a role in primordial germ cell and spermatogonial proliferation, but its role in the ovary is unclear, as it did not induce oocyte proliferation in trout. Based on increased gsdf transcript levels during secondary growth in coho salmon and the ability of Gsdf to stimulate germ cell proliferation in trout, it is possible this factor plays a role in granulosa cell proliferation that occurs during this period.
This study sheds light on differentially expressed ovarian genes during previtellogenic oocyte growth in coho salmon and provides a platform for future studies on the regulation of this process. Major gene families represented in the SSH libraries included zp genes, lipoprotein receptors, yolk proteases, and CA components, most of which appear to be derived from the oocyte. Transcript abundance was measured for selected genes identified by SSH and candidate genes likely expressed in follicle cells, providing for some genes the first transcriptional profile during early oogenesis in fish. Interestingly, a number of oogenesis related transcripts that will not be utilized until 1–3 years later in coho salmon were highly expressed during primary growth. Clearly, further studies are necessary to determine when these mRNAs are translated and if these transcripts may be subject to masking during early oogenesis. Finally, we observed increased expression of genes encoding FSH receptor, TGFβ family peptides, proteins involved in lipid uptake, and a CA component during the transition from primary to secondary oocyte growth. This period in coho salmon coincides with increased FSH signaling. However, the degree to which FSH regulates differentially expressed genes identified in this study is not known and will be the subject of our future investigations.
Tyler CR, Sumpter JP: Oocyte growth and development in teleosts. Rev Fish Biol Fish. 1996, 6: 287-318. 10.1007/BF00122584.
Patiño R, Sullivan CV: Ovarian follicle growth, maturation, and ovulation in teleost fish. Fish Physiol Biochem. 2002, 26: 57-70. 10.1023/A:1023311613987.
Wallace RA, Selman K: Cellular and dynamic aspects of oocyte growth in teleosts. Am Zool. 1981, 21: 325-343.
Wallace RA, Selman K: Ultrastructural aspects of oogenesis and oocyte growth in fish and amphibians. J Electron Microsc Tech. 1990, 16: 175-201. 10.1002/jemt.1060160302.
Pelegri F: Maternal factors in zebrafish development. Dev Dyn. 2003, 228: 535-554. 10.1002/dvdy.10390.
Song JL, Wessel GM: How to make an egg: transcriptional regulation in oocytes. Differentiation. 2005, 73: 1-17. 10.1111/j.1432-0436.2005.07301005.x.
Wessel GM, Brooks JM, Green E, Haley S, Voronina E, Wong J, Zaydfudim V, Conner S: The biology of cortical granules. Int Rev Cytol. 2001, 209: 117-206.
Zeng S, Gong Z: Expressed sequence tag analysis of expression profiles of zebrafish testis and ovary. Gene. 2002, 294: 45-53. 10.1016/S0378-1119(02)00791-6.
Li Y, Chia JM, Bartfai R, Christoffels A, Yue GH, Ding K, Ho MY, Hill JA, Stupka E, Orban L: Comparative analysis of the testis and ovary transcriptomes in zebrafish by combining experimental and computational tools. Comp Funct Genom. 2004, 5: 403-418. 10.1002/cfg.418.
Rise ML, von Schalburg KR, Brown GD, Mawer MA, Devlin RH, Kuipers N, Busby M, Beetz-Sargent M, Alberto R, Gibbs AR, Hunt P, Shukin R, Zeznik JA, Nelson C, Jones SRM, Smailus DE, Jones SJM, Schein JE, Marra MA, Butterfield YSN, Stott JM, Ng SHS, Davidson WS, Koop BF: Development and application of a salmonid EST database and cDNA microarray: data mining and interspecific hybridization characteristics. Genome Res. 2004, 14: 478-490. 10.1101/gr.1687304.
von Schalburg KR, Rise ML, Brown GD, Davidson WS, Koop BF: A comprehensive survey of the genes involved in maturation and development of the rainbow trout ovary. Biol Reprod. 2005, 72: 687-699. 10.1095/biolreprod.104.034967.
Govoroun M, Le Gac F, Guiguen Y: Generation of a large scale repertoire of expressed sequence tags (ESTs) from normalized rainbow trout cDNA libraries. BMC Genom. 2006, 7: 196-10.1186/1471-2164-7-196.
Knoll-Gellida A, André M, Gattegno T, Forgue J, Admon A, Babin PJ: Molecular phenotype of zebrafish ovarian follicle by serial analysis of gene expression and proteomic profiling, and comparison with the transcriptomes of other animals. BMC Genom. 2006, 7: 46-10.1186/1471-2164-7-46.
Kanamori A: Systematic identification of genes expressed during early oogenesis in medaka. Mol Reprod Dev. 2000, 55: 31-36. 10.1002/(SICI)1098-2795(200001)55:1<31::AID-MRD5>3.0.CO;2-7.
Baron D, Houlgatte R, Fostier A, Guiguen Y: Large-scale temporal gene expression profiling during gonadal differentiation and early gametogenesis in rainbow trout. Biol Reprod. 2005, 73: 959-966. 10.1095/biolreprod.105.041830.
Wen C, Zhang Z, Ma W, Xu M, Wen Z, Peng J: Genome-wide identification of female-enriched genes in zebrafish. Dev Dynam. 2005, 232: 171-179. 10.1002/dvdy.20210.
Bobe J, Nguyen T, Jalabert B: Targeted gene expression profiling in the rainbow trout (Oncorhynchus mykiss) ovary during maturational competence acquisition and oocyte maturation. Biol Reprod. 2004, 71: 73-82. 10.1095/biolreprod.103.025205.
Ge W: Intrafollicular paracrine communication in the zebrafish ovary: the state of the art of an emerging model for the study of vertebrate folliculogenesis. Mol Cell Endocrinol. 2005, 237: 1-10. 10.1016/j.mce.2005.03.012.
Modig C, Modesto T, Canario A, Cerdà J, von Hofsten J, Olsson P-E: Molecular characterization and expression pattern of zona pellucida proteins in gilthead seabream (Sparus aurata). Biol Reprod. 2006, 75: 717-725. 10.1095/biolreprod.106.050757.
Perazzolo LM, Coward K, Davail B, Normand E, Tyler CR, Pakdel F, Schneider WJ, Le Menn F: Expression and localization of messenger ribonucleic acid for the vitellogenin receptor in ovarian follicles throughout oogenesis in the rainbow trout, Oncorhynchus mykiss. Biol Reprod. 1999, 60: 1057-1068. 10.1095/biolreprod60.5.1057.
Campbell B, Dickey J, Beckman B, Young G, Pierce A, Fukada H, Swanson P: Previtellogenic oocyte growth in salmon: relationships among body growth, plasma insulin-like growth factor-1, estradiol-17beta, follicle-stimulating hormone and expression of ovarian genes for insulin-like growth factors, steroidogenic-acute regulatory protein and receptors for gonadotropins, growth hormone, and somatolactin. Biol Reprod. 2006, 75: 34-44. 10.1095/biolreprod.105.049494.
Huang X, Madan A: CAP3: a DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.
The Gene Ontology Consortium: Gene ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.
Spargo SC, Hope RM: Evolution and nomenclature of the zona pellucida gene family. Biol Reprod. 2003, 68: 358-362. 10.1095/biolreprod.102.008086.
Goetz FW, Iliev DB, McCauley LAR, Liarte CQ, Tort LB, Planas JV, MacKenzie S: Analysis of genes isolated from lipopolysaccharide-stimulated rainbow trout (Oncorhynchus mykiss) macrophages. Mol Immunol. 2004, 41: 1199-1210. 10.1016/j.molimm.2004.06.005.
Luckenbach JA, Early LW, Rowe AH, Borski RJ, Daniels HV, Godwin J: Aromatase cytochrome P450: cloning, intron variation, and ontogeny of gene expression in southern flounder (Paralichthys lethostigma). J Exp Zoolog A Comp Exp Biol. 2005, 303 (8): 643-656. 10.1002/jez.a.198.
Hake LE, Richter JD: Translational regulation of maternal mRNA. Biochem Biophys Acta. 1997, 1332: M31-M38.
Duranthon V, Renard J-P: Storage and functional recovery of molecules in oocytes. Biology and Pathology of the Oocyte, Its Role in Fertility and Reproductive Medicine. Edited by: Trounson AO, Gosden RG. 2003, Cambridge, UK: Cambridge University Press, 81-112.
Luckenbach JA, Goetz FW, Swanson P: Expressed sequence tags (ESTs) of follicle/interstitial cell enriched ovarian tissue from previtellogenic coho salmon. Cybium.
Goetz FW, McCauley L, Goetz GW, Norberg B: Using global genome approaches to address problems in cod mariculture. ICES J Mar Sci. 2006, 63: 393-399. 10.1016/j.icesjms.2005.10.006.
Chang YS, Hsu CC, Wang SC, Tsao CC, Huang FL: Molecular cloning, structural analysis, and expression of carp ZP2 gene. Mol Reprod Dev. 1997, 46: 258-267. 10.1002/(SICI)1098-2795(199703)46:3<258::AID-MRD4>3.0.CO;2-O.
Kusakabe M, Nakamura I, Young G: 11β-hydroxysteroid dehydrogenase complementary deoxyribonucleic acid in rainbow trout: cloning, sites of expression, and seasonal changes in gonads. Endocrinology. 2003, 144: 2534-2545. 10.1210/en.2002-220446.
Agulleiro MJ, André M, Morais S, Cerdà J, Babin PJ: High transcript level of fatty acid-binding protein 11 but not of very low-density lipoprotein receptor is correlated to ovarian follicle atresia in a teleost fish (Solea senegalensis). Biol Reprod. 2007, 77: 504-516. 10.1095/biolreprod.107.061598.
Mittelholzer C, Andersson E, Consten D, Hirai T, Nagahama Y, Norberg B: 20β-hydroxysteroid dehydrogenase and CYP19A1 are differentially expressed during maturation in Atlantic cod (Gadus morhua). J Mol Endocrinol. 2007, 39: 319-328. 10.1677/JME-07-0070.
Liang L, Soyal SM, Dean J: FIGα, a germ cell specific transcription factor involved in the coordinate expression of the zona pellucida genes. Development. 1997, 124: 4939-4947.
Onichtchouk D, Aduroja K, Belting H-G, Gnügge L, Driever W: Transgene driving GFP expression from the promoter of the zona pellucida gene zpc is expressed in oocytes and provides an early marker for gonad differentiation in zebrafish. Dev Dyn. 2003, 228: 393-404. 10.1002/dvdy.10392.
Prat F, Coward K, Sumpter JP, Tyler CR: Molecular characterization and expression of two ovarian lipoprotein receptors in the rainbow trout, Oncorhynchus mykiss. Biol Reprod. 1998, 58: 1146-1153. 10.1095/biolreprod58.5.1146.
Li A, Sadasivam M, Ding JL: Receptor-ligand interaction between vitellogenin receptor (VtgR) and vitellogenin (Vtg), implications on low density lipoprotein receptor and apolipoprotein B/E. J Biol Chem. 2003, 278: 2799-2806. 10.1074/jbc.M205067200.
Baron D, Guiguen Y: Gene expression during gonadal sex differentiation in rainbow trout (Oncorhynchus mykiss): from candidate gene studies to high throughout genomic approach. Fish Physiol Biochem. 2003, 28: 119-123. 10.1023/B:FISH.0000030496.21741.6c.
Hiramatsu N, Chapman RW, Lindzey JK, Haynes MR, Sullivan CV: Molecular characterization and expression of vitellogenin receptor from white perch (Morone americana). Biol Reprod. 2004, 70: 1720-1730. 10.1095/biolreprod.103.023655.
Black D, Skinner ER: Changes in plasma lipoproteins and tissue lipoprotein lipase and salt-resistant lipase activities during spawning in the rainbow trout (Salmo gairdnerii R.). Comp Biochem Physiol B. 1987, 88 (1): 261-267. 10.1016/0305-0491(87)90111-8.
Kwon JY, Prat F, Randall C, Tyler CR: Molecular characterization of putative yolk processing enzymes and their expression during oogenesis and embryogenesis in rainbow trout (Oncorhynchus mykiss). Biol Reprod. 2001, 65: 1701-1709. 10.1095/biolreprod65.6.1701.
Ibáñez AJ, Peinado-Onsurbe J, Sánchez E, Prat F: The role of lipoprotein lipase (LPL) in the incorporation of neutral lipids into the oocytes of the European sea bass (Dicentrarchus labrax L.) during gonadal development. Fish Physiol Biochem. 2003, 28: 291-293. 10.1023/B:FISH.0000030558.82463.99.
Carnevali O, Cionna C, Tosti L, Lubzens E, Maradonna F: Role of cathepsins in ovarian follicle growth and maturation. Gen Comp Endocrinol. 2006, 146: 195-203. 10.1016/j.ygcen.2005.12.007.
Wood AW, Van Der Kraak G: Yolk proteolysis in rainbow trout oocytes after serum-free culture: evidence for a novel biochemical mechanism of atresia in oviparous vertebrates. Mol Reprod Dev. 2003, 65: 219-227. 10.1002/mrd.10272.
Goetz FW: The "ups" and "downs" in using subtractive cloning techniques to isolate regulated genes in fish. Integr Comp Biol. 2003, 43: 786-793. 10.1093/icb/43.6.786.
Matsubara T, Nagae M, Ohkubo N, Andoh T, Sawaguchi S, Hiramatsu N, Sullivan CV, Hara A: Multiple vitellogenins and their unique roles in marine teleosts. Fish Physiol Biochem. 2003, 28: 295-299. 10.1023/B:FISH.0000030559.71954.37.
Fabra M, Cerdà J: Ovarian cysteine proteinases in the teleost Fundulus heteroclitus: molecular cloning and gene expression during vitellogenesis and oocyte maturation. Mol Reprod Dev. 2004, 67: 282-294. 10.1002/mrd.20018.
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: 377-385. 10.1002/mrd.20361.
Brooks S, Tyler CR, Carnevali O, Coward K, Sumpter JP: Molecular characterization of ovarian cathepsin D in the rainbow trout, Oncorhynchus mykiss. Gene. 1997, 201: 45-54. 10.1016/S0378-1119(97)00426-5.
Tateno H, Saneyoshi A, Ogawa T, Muramoto K, Kamiya H, Saneyoshi M: Isolation and characterization of rhamnose-binding lectins from eggs of steelhead trout (Oncorhynchus mykiss) homologous to low density lipoprotein receptor superfamily. J Biol Chem. 1998, 273: 19190-19197. 10.1074/jbc.273.30.19190.
Dong C-H, Yang S-T, Yang Z-A, Zhang L, Gui J-F: A C-type lectin associated and translocated with cortical granules during oocyte maturation and egg fertilization in fish. Dev Biol. 2004, 265: 341-354. 10.1016/j.ydbio.2003.08.028.
Tateno H, Yamaguchi T, Ogawa T, Muramoto K, Watanabe T, Kamiya H, Saneyoshi M: Immunohistochemical localization of rhamnose-binding lectins in the steelhead trout (Oncorhynchus mykiss). Dev Comp Immunol. 2002, 26: 543-550. 10.1016/S0145-305X(02)00007-1.
Shibata Y, Iwamatsu T, Oba Y, Kobayashi D, Tanaka M, Nagahama Y, Suzuki N, Yoshikuni M: Identification and cDNA cloning of alveolin, an extracellular metalloproteinase, which induces chorion hardening of medaka (Oryzias latipes) eggs upon fertilization. J Biol Chem. 2000, 275: 8349-8354. 10.1074/jbc.275.12.8349.
Asahina S, Sato C, Kitajima K: Developmental expression of a sialyltransferase responsible for sialylation of cortical alveolus glycoprotein during oogenesis in rainbow trout (Oncorhynchus mykiss). J Biochem. 2004, 136: 189-198. 10.1093/jb/mvh106.
Crichton RR: Proteins of iron storage and transport. Adv Protein Chem. 1990, 40: 281-363.
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: 127-135. 10.1002/mrd.10384.
Yoshida N, Yamashita M: Non-dependence of cyclin E/Cdk2 kinase activity on the initiation of oocyte maturation in goldfish. Dev Growth Differ. 2000, 42: 285-294. 10.1046/j.1440-169x.2000.00506.x.
Sammar M, Levi L, Hurvitz A, Lubzens E: Studies on retinol-binding protein during vitellogenesis in the rainbow trout (Oncorhynchus mykiss). Gen Comp Endocrinol. 2005, 141: 141-151. 10.1016/j.ygcen.2004.12.014.
Schekman R, Orci L: Coat proteins and vesicle budding. Science. 1996, 271: 1526-1533. 10.1126/science.271.5255.1526.
Swanson P, Bernard M, Nozaki M, Suzuki K, Kawauchi H, Dickhoff WW: Gonadotropins I and II in juvenile coho salmon. Fish Physiol Biochem. 1989, 7: 169-176. 10.1007/BF00004704.
Swanson P, Suzuki K, Kawauchi H, Dickhoff WW: Isolation and characterization of coho salmon gonadotropins, GTH I and GTH II. Biol Reprod. 1991, 44: 29-38. 10.1095/biolreprod44.1.29.
Swanson P: Salmon gonadotropins: reconciling old and new ideas. Proceedings of the Fourth International Symposium on the Reproductive Physiology of Fish. Edited by: Scott AP, Sumpter JP, Kime DE, Rolfe MS. 1991, University of East Anglia, Norwich, England, 2-7.
Visser JA, de Jong FH, Laven JSE, Themmen APN: Anti-Müllerian hormone: a new marker for ovarian function. Reproduction. 2006, 131: 1-9. 10.1530/rep.1.00529.
Rodríguez-Marí A, Yan Y-L, BreMiller RA, Wilson C, Cañestro C, Postlethwait JH: Characterization and expression pattern of zebrafish anti-Müllerian hormone (amh) relative to sox9a, sox 9b, and cyp19a1a, during gonad development. Gene Expression Patterns. 2005, 5: 655-667. 10.1016/j.modgep.2005.02.008.
Halm S, Rocha A, Miura T, Prat F, Zanuy S: Anti-müllerian hormone (AMH/AMH) in the European sea bass: its gene structure, regulatory elements, and the expression of alternatively-spliced isoforms. Gene. 2007, 388: 148-158. 10.1016/j.gene.2006.10.018.
Sawatari E, Shikina S, Takeuchi T, Yoshizaki G: A novel transforming growth factor-β superfamily member expressed in gonadal somatic cells enhances primordial germ cell and spermatogonial proliferation in rainbow trout (Oncorhynchus mykiss). Dev Biol. 2007, 301: 266-275. 10.1016/j.ydbio.2006.10.001.
We thank G.W. Goetz for bioinformatics assistance, J.T. Dickey and T.J. Stringer for lab assistance, and B. Gadberry and P. Parkins for fish husbandry. This research was supported by a grant from the National Oceanic and Atmospheric Administration's Oceans and Human Health Initiative to PS.
JAL participated in the study design, execution of experiments, data analysis and interpretation, and drafted the manuscript. DBI performed the SSH and participated in data analysis. FWG participated in the study design and assisted with data analysis and interpretation. PS participated in the study design, acquired funding, and assisted with data analysis and interpretation. All the authors read, edited, and approved the final manuscript.