Functional Genomics and Proteomic Approaches for the Study of Gamete Formation and Viability in Farmed Finfish

Genome sequencing and analysis of black flounder (Paralichthys orbignyanus) reveals new insights into Pleuronectiformes genomic size and structure

Black flounder (Paralichthys orbignyanus, Pleuronectiformes) is a commercially significant marine fish with promising aquaculture potential in Argentina. Despite extensive studies on Black flounder aquaculture, its limited genetic information available hampers the crucial role genetics plays in the development of this activity. In this study, we first employed Illumina sequencing technology to sequence the entire genome of Black flounder. Utilizing two independent libraries-one from a female and another from a male-with 150 bp paired-end reads, a mean insert length of 350 bp, and over 35 X-fold coverage, we achieved assemblies resulting in a genome size of ~ 538 Mbp. Analysis of the assemblies revealed that more than 98% of the core genes were present, with more than 78% of them having more than 50% coverage. This indicates a somehow complete and accurate genome at the coding sequence level. This genome contains 25,231 protein-coding genes, 445 tRNAs, 3 rRNAs, and more than 1,500 non-coding RNAs of other types. Black flounder, along with pufferfishes, seahorses, pipefishes, and anabantid fish, displays a smaller genome compared to most other teleost groups. In vertebrates, the number of transposable elements (TEs) is often correlated with genome size. However, it remains unclear whether the sizes of introns and exons also play a role in determining genome size. Hence, to elucidate the potential factors contributing to this reduced genome size, we conducted a comparative genomic analysis between Black flounder and other teleost orders to determine if the small genomic size could be explained by repetitive elements or gene features, including the whole genome genes and introns sizes. We show that the smaller genome size of flounders can be attributed to several factors, including changes in the number of repetitive elements, and decreased gene size, particularly due to lower amount of very large and small introns. Thus, these components appear to be involved in the genome reduction in Black flounder. Despite these insights, the full implications and potential benefits of genome reduction in Black flounder for reproduction and aquaculture remain incompletely understood, necessitating further research.

Glucocorticoids in Fish Eggs: Variation, Interactions with the Environment, and the Potential to Shape Offspring Fitness

Wild and captive vertebrates face multiple stressors that all have the potential to induce chronic maternal stress (i.e., sustained, elevated plasma glucocorticoids), resulting in embryo exposure to elevated maternally derived glucocorticoids. In oviparous taxa such as fish, maternally derived glucocorticoids in eggs are known for their capacity to shape offspring phenotype. Using a variety of methodologies, scientists have quantified maternally derived levels of egg cortisol, the primary glucocorticoid in fishes, and examined the cascading effects of egg cortisol on progeny phenotype. Here we summarize and interpret the current state of knowledge on egg cortisol in fishes and the relationships linking maternal stress/state to egg cortisol and offspring phenotype/fitness. Considerable variation in levels of egg cortisol exists across species and among females within a species; this variation is hypothesized to be due to interspecific differences in reproductive life history and intraspecific differences in female condition. Outcomes of experimental studies manipulating egg cortisol vary both inter- and intraspecifically. Moreover, while exogenous elevation of egg cortisol (as a proxy for maternal stress) induces phenotypic changes commonly considered to be maladaptive (e.g., smaller offspring size), emerging work in other taxa suggests that there can be positive effects on fitness when the offspring's environment is taken into account. Investigations into (i) mechanisms by which egg cortisol elicits phenotypic change in offspring (e.g., epigenetics), (ii) maternal and offspring buffering capacity of cortisol, and (iii) factors driving natural variation in egg cortisol and how this variation affects offspring phenotype and fitness are all germane to discussions on egg glucocorticoids as signals of maternal stress.

Proteomic analysis of early-stage embryos: implications for egg quality in hapuku (Polyprion oxygeneios)

In order to develop biomarkers that may help predict the egg quality of captive hapuku (Polyprion oxygeneios) and provide potential avenues for its manipulation, the present study (1) sequenced the proteome of early-stage embryos using isobaric tag for relative and absolute quantification analysis, and (2) aimed to establish the predictive value of the abundance of identified proteins with regard to egg quality through regression analysis. Egg quality was determined for eight different egg batches by blastomere symmetry scores. In total, 121 proteins were identified and assigned to one of nine major groups according to their function/pathway. A mixed-effects model analysis revealed a decrease in relative protein abundance that correlated with (decreasing) egg quality in one major group (heat-shock proteins). No differences were found in the other protein groups. Linear regression analysis, performed for each identified protein separately, revealed seven proteins that showed a significant decrease in relative abundance with reduced blastomere symmetry: two correlates that have been named in other studies (vitellogenin, heat-shock protein-70) and a further five new candidate proteins (78 kDa glucose-regulated protein, elongation factor-2, GTP-binding nuclear protein Ran, iduronate 2-sulfatase and 6-phosphogluconate dehydrogenase). Notwithstanding issues associated with multiple statistical testing, we conclude that these proteins, and especially iduronate 2-sulfatase and the generic heat-shock protein group, could serve as biomarkers of egg quality in hapuku.

Multiple vitellogenins and product yolk proteins in European sea bass (Dicentrarchus labrax): Molecular characterization, quantification in plasma, liver and ovary, and maturational proteolysis

Three complete vitellogenin (Vtg) polypeptides of European sea bass (Dicentrarchus labrax), an acanthomorph teleost spawning pelagic eggs in seawater, were deduced from cDNA and identified as VtgAa, VtgAb and VtgC based on current Vtg nomenclature and phylogeny. Label free quantitative mass spectrometry verified the presence of the three sea bass Vtgs or their product yolk proteins (YPs) in liver, plasma and ovary of postvitellogenic females. As evidenced by normalized spectral counts, VtgAb-derived protein was 2- to 5-fold more abundant, depending on sample type, than for VtgAa, while VtgC-derived protein was less abundant, albeit only 3-fold lower than for VtgAb in the ovary. Western blotting with Vtg type-specific antisera raised against corresponding gray mullet (Mugil cephalus) lipovitellins (Lvs) detected all three types of sea bass Vtg in the blood plasma of gravid females and/or estrogenized males and showed that all three forms of sea bass Lv undergo limited partial degradation during oocyte maturation. The comparatively high levels of VtgC-derived YPs in fully-grown oocytes and the maturational proteolysis of all three types of Lv differ from what has been reported for other teleosts spawning pelagic eggs in seawater but are similar to recent findings for two species of North American Moronidae, the striped bass (Morone saxatilis) and white perch (Morone americana), which spawn pelagic and demersal eggs, respectively in fresh water. Together with the high Vtg sequence homologies and virtually identical structural features of each type of Vtg between species, these findings indicate that the moronid multiple Vtg systems do not substantially vary with reproductive environment.

Gonadal Transcriptome Analysis in Sterile Double Haploid Japanese Flounder

Sterility is a serious problem that can affect all bionts. In teleosts, double haploids (DHs) induced by mitogynogenesis are often sterile. This sterility severely restricts the further application of DHs for production of clones, genetic analysis, and breeding. However, sterile DH individuals are good source materials for investigation of the molecular mechanisms of gonad development, especially for studies into the role of genes that are indispensable for fish reproduction. Here, we used the Illumina sequencing platform to analyze the transcriptome of sterile female DH Japanese flounder in order to identify major genes that cause sterility and to provide a molecular basis for an intensive study of gonadal development in teleosts. Through sequencing, assembly, and annotation, we obtained 52,474 contigs and found that 60.7% of these shared homologies with existing sequences. A total of 1225 differentially expressed unigenes were found, including 492 upregulated and 733 downregulated genes. Gene Ontology and KEGG analyses showed that genes showing significant upregulation, such as CYP11A1, CYP11B2, CYP17, CYP21, HSD3β, bcl2l1, and PRLR, principally correlated with sterol metabolic process, steroid biosynthetic process, and the Jak-stat signaling pathway. The significantly downregulated genes were primarily associated with immune response, antigen processing and presentation, cytokine–cytokine receptor interaction, and protein digestion and absorption. Using a co-expression network analysis, we conducted a comprehensive comparison of gene expression in the gonads of fertile and sterile female DH Japanese flounder. Identification of genes showing significantly different expression will provide further insights into DH reproductive dysfunction and oocyte maturation processes in teleosts.

Estrogen-induced yolk precursors in European sea bass, Dicentrarchus labrax: Status and perspectives on multiplicity and functioning of vitellogenins

The estrogen-inducible egg yolk precursor, vitellogenin, of the European sea bass (Dicentrarchus labrax) has received considerable scientific attention by virtue of its central importance in determination of oocyte growth and egg quality in this important aquaculture species. However, the multiplicity of vitellogenins in the sea bass has only recently been examined. Recent cloning and homology analyses have revealed that the sea bass possesses the three forms of vitellogenin, VtgAa, VtgAb and VtgC, reported to occur in some other highly evolved teleosts. Progress has been made in assessing the relative abundance and special structural features of the three Vtgs and their likely roles in oocyte maturation and embryonic nutrition. This report discusses these findings in the context of our prior knowledge of vitellogenesis in this species and of the latest advances in our understanding of the evolution and function of multiple Vtgs in acanthomorph fishes.

Construction, De-Novo Assembly and Analysis of Transcriptome for Identification of Reproduction-Related Genes and Pathways from Rohu, Labeo rohita (Hamilton)

Rohu is a leading candidate species for freshwater aquaculture in South-East Asia. Unlike common carp the monsoon breeding habit of rohu restricts its seed production beyond season indicating strong genetic control over spawning. Genetic information is limited in this regard. The problem is exacerbated by the lack of genomic-resources. We identified 182 reproduction-related genes previously by Sanger-sequencing which were less to address the issue of seasonal spawning behaviour of this important carp. Therefore, the present work was taken up to generate transcriptome profile by mRNAseq. 16 GB, 72 bp paired end (PE) data was generated from the pooled-RNA of twelve-tissues from pre-spawning rohu using IlluminaGA-II-platform. There were 64.97 million high-quality reads producing 62,283 contigs and 88,612 numbers of transcripts using velvet and oases programs, respectively. Gene ontology annotation identified 940 reproduction-related genes consisting of 184 mainly associated with reproduction, 223 related to hormone-activity and receptor-binding, 178 receptor-activity and 355 embryonic-development related-proteins. The important reproduction-relevant pathways found in KEGG analysis were GnRH-signaling, oocyte-meiosis, steroid-biosynthesis, steroid-hormone biosynthesis, progesterone-mediated oocyte-maturation, retinol-metabolism, neuroactive-ligand-receptor interaction, neurotrophin-signaling and photo-transduction. Twenty nine simple sequence repeat containing sequences were also found out of which 12 repeat loci were polymorphic with mean expected-&-observed heterozygosity of 0.471 and 0.983 respectively. Quantitative RT-PCR analyses of 13-known and 6-unknown transcripts revealed differences in expression level between preparatory and post-spawning phase. These transcriptomic sequences have significantly increased the genetic-&-genomic resources for reproduction-research in Labeo rohita.

Expression profiles of Fsh-regulated ovarian genes during oogenesis in coho salmon

The function of follicle-stimulating hormone (Fsh) during oogenesis in fishes is poorly understood. Using coho salmon as a fish model, we recently identified a suite of genes regulated by Fsh in vitro and involved in ovarian processes mostly unexplored in fishes, like cell proliferation, differentiation, survival or extracellular matrix (ECM) remodeling. To better understand the role of these Fsh-regulated genes during oocyte growth in fishes, we characterized their mRNA levels at discrete stages of the ovarian development in coho salmon. While most of the transcripts were expressed at low levels during primary growth (perinucleolus stage), high expression of genes associated with cell proliferation (pim1, pcna, and mcm4) and survival (ddit4l) was found in follicles at this stage. The transition to secondary oocyte growth (cortical alveolus and lipid droplet stage ovarian follicles) was characterized by a marked increase in the expression of genes related to cell survival (clu1, clu2 and ivns1abpa). Expression of genes associated with cell differentiation and growth (wt2l and adh8l), growth factor signaling (inha), steroidogenesis (cyp19a1a) and the ECM (col1a1, col1a2 and dcn) peaked in vitellogenic follicles, showing a strong and positive correlation with transcripts for fshr. Other genes regulated by Fsh and associated with ECM function (ctgf, wapl and fn1) and growth factor signaling (bmp16 and smad5l) peaked in maturing follicles, along with increases in steroidogenesis-related gene transcripts. In conclusion, ovarian genes regulated by Fsh showed marked differences in their expression patterns during oogenesis in coho salmon. Our results suggest that Fsh regulates different ovarian processes at specific stages of development, likely through interaction with other intra- or extra-ovarian factors.

Global transcriptome analysis identifies regulated transcripts and pathways activated during oogenesis and early embryogenesis in Atlantic cod

The molecular mechanisms underlying oogenesis and maternally controlled embryogenesis in fish are not fully understood, especially in marine species. Our aim was to study the egg and embryo transcriptome during oogenesis and early embryogenesis in Atlantic cod. Follicles from oogenesis stages (pre-, early-, and late-vitellogenic), ovulated eggs, and two embryonic stages (blastula, gastrula) were collected from broodstock fish and fertilized eggs. Gene expression profiles were measured in a 44 K oligo microarray consisting of 23,000 cod genes. Hundreds of differentially expressed genes (DEGs) were identified in the follicle stages investigated, implicating a continuous accumulation and degradation of polyadenylated transcripts throughout oogenesis. Very few DEGs were identified from ovulated egg to blastula, showing a more stable maternal RNA pool in early embryonic stages. The highest induction of expression was observed between blastula and gastrula, signifying the onset of zygotic transcription. During early vitellogenesis, several of the most upregulated genes are linked to nervous system signaling, suggesting increasing requirements for ovarian synaptic signaling to stimulate the rapid growth of oocytes. Highly upregulated genes during late vitellogenesis are linked to protein processing, fat metabolism, osmoregulation, and arrested meiosis. One of the genes with the highest upregulation in the ovulated egg is involved in oxidative phosphorylation, reflecting increased energy requirements during fertilization and the first rapid cell divisions of early embryogenesis. In conclusion, this study provides a large-scale presentation of the Atlantic cod's maternally controlled transcriptome in ovarian follicles through oogenesis, ovulated eggs, and early embryos. Mol. Reprod. Dev. 81: 619–635, 2014. © 2014 Wiley Periodicals, Inc.

Genomic resources for flatfish research and their applications

Flatfishes are a group of teleosts of high commercial and environmental interest, whose biology is still poorly understood. The recent rapid development of different 'omic' technologies is, however, enhancing the knowledge of the complex genetic control underlying different physiological processes of flatfishes. This review describes the different functional genomic approaches and resources currently available for flatfish research and summarizes different areas where microarray-based gene expression analysis has been applied. The increase in genome sequencing data has also allowed the construction of genetic linkage maps in different flatfish species; these maps are invaluable for investigating genome organization and identifying genetic traits of commercial interest. Despite the significant progress in this field, the genomic resources currently available for flatfish are still scarce. Further intensive research should be carried out to develop larger genomic sequence databases, high-density microarrays and, more detailed, complete linkage maps, using second-generation sequencing platforms. These tools will be crucial for further expanding the knowledge of flatfish physiology, and it is predicted that they will have important implications for wild fish population management, improved fish welfare and increased productivity in aquaculture.

Oogenesis in teleosts: how eggs are formed

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

Molecular pathways during marine fish egg hydration: the role of aquaporins

The pre-ovulatory hydration of the oocyte of marine teleosts, a unique process among vertebrates that occurs concomitantly with meiosis resumption (oocyte maturation), is a critical process for the correct development and survival of the embryo. Increasing information is available on the molecular mechanisms that control oocyte maturation in fish, but the identification of the cellular processes involved in oocyte hydration has remained long ignored. During the past few years, a number of studies have identified the major inorganic and organic osmolytes that create a transient intra-oocytic osmotic potential for hydrating the oocytes, whereas water influx was believed to occur passively. Recent work, however, has uncovered the role of a novel molecular water channel (aquaporin), designated aquaporin-1b (Aqp1b), which facilitates water permeation and resultant swelling of the oocyte. The Aqp1b belongs to a teleost-specific subfamily of water-selective aquaporins, similar to mammalian aquaporin-1 (AQP1) that has possibly evolved by duplication of a common ancestor and further neofunctionalization in oocytes of marine teleosts for water uptake. Strikingly, Aqp1b shows specific regulatory domains at the cytoplasmic tail, which are key to the vesicular trafficking and temporal insertion of Aqp1b in the oocyte plasma membrane during the phase of maximal hydration. These findings are revealing that the mechanism of oocyte hydration in marine teleosts is a highly regulated process based on the interplay between the generation of inorganic and organic osmolytes and the controlled insertion of Aqp1b in the oocyte surface. The discovery of Aqp1b in teleosts provides an important insight into the molecular basis of the production of viable eggs in marine fish.

Maternal transfer and transcriptional onset of immune genes during ontogenesis in Atlantic cod

The immune system in teleosts is not completely developed during embryonic and larval stages and immune competence is assumed to be restricted. This study is the first to address whether immune transcripts are maternally transferred to offspring and when immune genes are transcriptionally active in Atlantic cod (Gadus morhua). In unfertilised eggs, transcripts encoding lysozyme and cathelicidin were found indicating maternal transfer of antibacterial transcripts. Lysozyme activity was also present at this stage suggesting the presence of a functional protein. Transcripts of two other putative antibacterial genes (hepcidin and pentraxin) and antiviral genes (ISG15 and LGP2) were absent in unfertilised eggs. The transcriptional onset of these genes occurred during the gastrula period. Transcripts of the heavy chain constant regions of the immunoglobulin (Ig) D, membrane-associated and secreted form of IgM were absent in unfertilised eggs. Transcription of the heavy chain locus commenced at low levels during the segmentation period indicating the onset of B-cell development. Most innate immune genes showed an increase in transcription around hatch and first feeding, indicating a preparation for increased pathogen exposure at this time. Prior to and during metamorphosis all genes showed a pronounced elevation in transcript levels indicating a further maturation of the immune system during this period.

New insights into molecular pathways associated with flatfish ovarian development and atresia revealed by transcriptional analysis

Background

The Senegalese sole (Solea senegalensis) is a marine flatfish of increasing commercial interest. However, the reproduction of this species in captivity is not yet controlled mainly because of the poor knowledge on its reproductive physiology, as it occurs for other non-salmonid marine teleosts that exhibit group-synchronous ovarian follicle development. In order to investigate intra-ovarian molecular mechanisms in Senegalese sole, the aim of the present study was to identify differentially expressed genes in the ovary during oocyte growth (vitellogenesis), maturation and ovarian follicle atresia using a recently developed oligonucleotide microarray.

Results

Microarray analysis led to the identification of 118 differentially expressed transcripts, of which 20 and 8 were monitored by real-time PCR and in situ hybridization, respectively. During vitellogenesis, many up-regulated ovarian transcripts had putative mitochondrial function/location suggesting high energy production (NADH dehydrogenase subunits, cytochromes) and increased antioxidant protection (selenoprotein W2a), whereas other regulated transcripts were related to cytoskeleton and zona radiata organization (zona glycoprotein 3, alpha and beta actin, keratin 8), intracellular signalling pathways (heat shock protein 90, Ras homolog member G), cell-to-cell and cell-to-matrix interactions (beta 1 integrin, thrombospondin 4b), and the maternal RNA pool (transducer of ERBB2 1a, neurexin 1a). Transcripts up-regulated in the ovary during oocyte maturation included ion transporters (Na⁺-K⁺-ATPase subunits), probably required for oocyte hydration, as well as a proteinase inhibitor (alpha-2-macroglobulin) and a vesicle calcium sensor protein (extended synaptotagmin-2-A). During follicular atresia, few transcripts were found to be up-regulated, but remarkably most of them were localized in follicular cells of atretic follicles, and they had inferred roles in lipid transport (apolipoprotein C-I), chemotaxis (leukocyte cell-derived chemotaxin 2,), angiogenesis (thrombospondin), and prevention of apoptosis (S100a10 calcium binding protein).

Conclusion

This study has identified a number of differentially expressed genes in the ovary that were not previously found to be regulated during ovarian development in marine fish. Specifically, we found evidence, for the first time in teleosts, of the activation of chemoattractant, angiogenic and antiapoptotic pathways in hypertrophied follicular cells at the onset of ovarian atresia.