Molecular phenotype of zebrafish ovarian follicle by serial analysis of gene expression and proteomic profiling, and comparison with the transcriptomes of other animals

Comparative proteomic analysis of the ovarian fluid and eggs of Siberian sturgeon

BACKGROUND

Sturgeon species are living fossils that exhibit unique reproductive characteristics, and elucidation of the molecular processes governing the formation and quality of sturgeon eggs is crucial. However, comprehensive data on the protein composition of sturgeon ovarian fluid (OF) and eggs and their functional significance are lacking. To address this knowledge gap, the aim of the present study was to conduct a comprehensive comparative proteomic analysis of Siberian sturgeon OF and eggs using liquid chromatography-mass spectrometry (LC-MS/MS).

RESULTS

A total of 617 proteins were identified in OF, and 565 proteins were identified in eggs. A total of 772 proteins showed differential abundance. Among the differentially abundant proteins, 365 were more abundant in OFs, while 407 were more abundant in eggs. We identified 339 proteins unique to OFs and 287 proteins specific to eggs, and further investigated the top 10 most abundant proteins in each. The functional annotation of the OF proteins highlighted their predominant association with immune system processes, including the complement and coagulation cascade, neutrophil and leukocyte-mediated immunity, cholesterol metabolism, and regulation of the actin cytoskeleton. Analysis of egg proteins revealed enrichment in metabolic pathways, such as oxidative phosphorylation and fatty acid metabolism, and protein ubiquitination and translation. OF-specific proteins included extracellular matrix and secretory vesicles, and eggs were enriched in proteins localized to mitochondria and ribosome components.

CONCLUSIONS

This study presents the first comprehensive characterization of the protein composition of sturgeon OF and eggs and elucidates their distinct functional roles. These findings advance our understanding of sturgeon reproduction, OF-egg signaling and the origin of OF proteins. The mass spectrometry proteomics data have been deposited in the ProteomeXchange Consortium with the dataset identifier PXD044168 to ensure accessibility for further research.

Sex blind: bridging the gap between drug exposure and sex-related gene expression in Danio rerio using next-generation sequencing (NGS) data and a literature review to find the missing links in pharmaceutical and environmental toxicology studies

The sex of both humans and Danio rerio has previously been shown to affect the way individuals respond to drug exposure. Genes which allow identification of sex in juvenile zebrafish show potential to reveal these confounding variables between sex in toxicological and preclinical trials but the link between these is so far missing. These sex-specific, early expressed genes where expression is not altered by drug exposure must be carefully selected for this purpose. We aimed to discover genes which can be used in pharmaceutical trials and environmental toxicology studies to uncover sex-related variations in gene expression with drug application using the model organism Danio rerio. Previously published early sex determining genes from King et al. were evaluated as well as additional genes selected from our zebrafish Next-generation sequencing (NGS) data which are known from previously published works not to be susceptible to changes in expression with drug exposure. NGS revealed a further ten female-specific genes (vtg1, cyp17a1, cyp19a1a, igf3, ftz-f1, gdf9, foxl2a, Nr0b1, ipo4, lhcgr) and five male related candidate genes (FKBP5, apobb1, hbaa1, dmrt1, spata6) which are also expressed in juvenile zebrafish, 28 days post fertilisation (dpf). Following this, a literature review was performed to classify which of these early-expressed sex specific genes are already known to be affected by drug exposure in order to determine candidate genes to be used in pharmaceutical trials or environmental toxicology testing studies. Discovery of these early sex-determining genes in Danio rerio will allow identification of sex-related responses to drug testing to improve sex-specific healthcare and the medical treatment of human patients.

Revelation of candidate genes and molecular mechanism of reproductive seasonality in female rohu (Labeo rohita Ham.) by RNA sequencing

BACKGROUND

Carp fish, rohu (Labeo rohita Ham.) is important freshwater aquaculture species of South-East Asia having seasonal reproductive rhythm. There is no holistic study at transcriptome level revealing key candidate genes involved in such circannual rhythm regulated by biological clock genes (BCGs). Seasonality manifestation has two contrasting phases of reproduction, i.e., post-spawning resting and initiation of gonadal activity appropriate for revealing the associated candidate genes. It can be deciphered by RNA sequencing of tissues involved in BPGL (Brain-Pituitary-Gonad-Liver) axis controlling seasonality. How far such BCGs of this fish are evolutionarily conserved across different phyla is unknown. Such study can be of further use to enhance fish productivity as seasonality restricts seed production beyond monsoon season.

RESULT

A total of ~ 150 Gb of transcriptomic data of four tissues viz., BPGL were generated using Illumina TruSeq. De-novo assembled BPGL tissues revealed 75,554 differentially expressed transcripts, 115,534 SSRs, 65,584 SNPs, 514 pathways, 5379 transcription factors, 187 mature miRNA which regulates candidate genes represented by 1576 differentially expressed transcripts are available in the form of web-genomic resources. Findings were validated by qPCR. This is the first report in carp fish having 32 BCGs, found widely conserved in fish, amphibian, reptile, birds, prototheria, marsupials and placental mammals. This is due to universal mechanism of rhythmicity in response to environment and earth rotation having adaptive and reproductive significance.

CONCLUSION

This study elucidates evolutionary conserved mechanism of photo-periodism sensing, neuroendocrine secretion, metabolism and yolk synthesis in liver, gonadal maturation, muscular growth with sensory and auditory perception in this fish. Study reveals fish as a good model for research on biological clock besides its relevance in reproductive efficiency enhancement.

Neurodevelopment vs. the immune system: Complementary contributions of maternally-inherited gene transcripts and proteins to successful embryonic development in fish

The aim of this study was to investigate the respective contribution of maternally-inherited mRNAs and proteins to egg molecular cargo and to its developmental competence in fish using pikeperch as a model. Our study provides novel insights into the understanding of type-specific roles of maternally-inherited molecules in fish. Here we show, for the first time, that transcripts and proteins have distinct, yet complementary, functions in the egg of teleost fish. Maternally-inherited mRNAs would shape embryo neurodevelopment, while maternally-inherited proteins would rather be responsible for protecting the embryo against pathogens. Additionally, we observed that processes directly preceding ovulation may considerably affect the reproductive success by modifying expression level of genes crucial for proper embryonic development, being novel fish egg quality markers (e.g., smarca4 or h3f3a). These results are of major importance for understanding the influence of external factors on reproductive fitness in both captive and wild-type fish species.

Shedding new light on early sex determination in zebrafish

In contrast to established zebrafish gene annotations, the question of sex determination has still not been conclusively clarified for developing zebrafish, Danio rerio, larvae, 28 dpf or earlier. Recent studies indicate polygenic sex determination (PSD), with the genes being distributed throughout the genome. Early genetic markers of sex in zebrafish help unravel co-founding sex-related differences to apply to human health and environmental toxicity studies. A qPCR-based method was developed for six genes: cytochrome P450, family 17, subfamily A, polypeptide 1 (cyp17a1); cytochrome P450, family 19, subfamily A, polypeptide 1a (cyp19a1a); cytochrome P450, family 19, subfamily A, polypeptides 1b (cyp19a1b); vitellogenin 1 (vtg1); nuclear receptor subfamily 0, group B, member 1 (nr0b1), sry (sex-determining region Y)-box 9b (sox9b) and actin, beta 1 (actb1), the reference gene. Sry-box 9a (Sox9a), insulin-like growth factor 3 (igf3) and double sex and mab-3 related transcription factor 1 (dmrt1), which are also known to be associated with sex determination, were used in gene expression tests. Additionally, Next-Generation-Sequencing (NGS) sequenced the genome of two adult female and male and two juveniles. PCR analysis of adult zebrafish revealed sex-specific expression of cyp17a1, cyp19a1a, vtg1, igf3 and dmrt1, the first four strongly expressed in female zebrafish and the last one highly expressed in male conspecifics. From NGS, nine female and four male-fated genes were selected as novel for assessing zebrafish sex, 28 dpf. Differences in transcriptomes allowed allocation of sex-specific genes also expressed in juvenile zebrafish.

Transcriptome Analysis of Maternal Gene Transcripts in Unfertilized Eggs of Misgurnus anguillicaudatus and Identification of Immune-Related Maternal Genes

Maternal genes are important in directing early development and determining egg quality in fish. We here report the de novo transcriptome from four tissue libraries of the cyprinid loach, Misgurnus anguillicaudatus, and for the first time identified maternal gene transcripts in unfertilized eggs and suggest their immune system involvement. Expression profiles and functional enrichment revealed a total 24,116 transcripts were expressed as maternal transcripts in unfertilized eggs, which were involved in a wide range of biological functions and pathways. Comparison expression profiles and analysis of tissue specificity revealed that the large numbers of maternal transcripts were stored in unfertilized eggs near the late phase of ovarian maturation and before ovulation. Functional classification showed a total of 279 maternal immune-related transcripts classified with immune system process GO term and immune system KEGG pathway. qPCR analysis showed that transcript levels of identified maternal immune-related candidate genes were dynamically modulated during development and early ontogeny of M. anguillicaudatus. Taken together, this study could not only provide knowledge on the protective roles of maternal immune-related genes during early life stage of M. anguillicaudatus but could also be a valuable transcriptomic/genomic resource for further analysis of maternally provisioned genes in M. anguillicaudatus and other related teleost fishes.

Effects of ovarian fluid on sperm traits and its implications for cryptic female choice in zebrafish

In polyandrous mating systems, females maintain the opportunity to bias male fertilization success after mating in a process known as cryptic female choice. Mechanisms of cryptic female choice have been described both in internal and external fertilizers, and may affect fertilization processes at different stages before, during, and after fertilization. In internal fertilizers, females have substantial control over sperm storage and fertilization, whereas in external fertilizers, female control is limited. A key factor proposed to mediate cryptic female choice is the fluid surrounding the eggs, the ovarian fluid, as it may directly affect sperm performance. Here, we studied the role of ovarian fluid in post-mating sexual selection using the zebrafish, Danio rerio. Firstly, we assessed how ovarian fluid affects sperm swimming performance compared with freshwater. We focused on sperm motility, velocity, swimming trajectory, and longevity, all traits associated with competitive fertilization success in externally fertilizing fish. In a second step, we used a North Carolina II design to explore female, male, and female x male effects by testing sperm motility of 2 males in the ovarian fluid of 2 females in a total of 11 blocks. Our results indicate that the ovarian fluid affects sperm performance differently from freshwater. Specifically, sperm velocity, motility, and longevity were higher in the ovarian fluid than in freshwater, whereas sperm linearity and beat cross frequency showed the opposite pattern. Moreover, these effects varied according to male, female, and male x female identities, supporting the potential for cryptic female choice mediated by ovarian fluid in this species.

Genetic mechanism underlying sexual plasticity and its association with colour patterning in zebrafish (Danio rerio)

BACKGROUND

Elevated water temperature, as is expected through climate change, leads to masculinization in fish species with sexual plasticity, resulting in changes in population dynamics. These changes are one important ecological consequence, contributing to the risk of extinction in small and inbred fish populations under natural conditions, due to male-biased sex ratio. Here we investigated the effect of elevated water temperature during embryogenesis on sex ratio and sex-biased gene expression profiles between two different tissues, namely gonad and caudal fin of adult zebrafish males and females, to gain new insights into the molecular mechanisms underlying sex determination (SD) and colour patterning related to sexual attractiveness.

RESULTS

Our study demonstrated sex ratio imbalances with 25.5% more males under high-temperature condition, resulting from gonadal masculinization. The result of transcriptome analysis showed a significantly upregulated expression of male SD genes (e.g. dmrt1, amh, cyp11c1 and sept8b) and downregulation of female SD genes (e.g. zp2.1, vtg1, cyp19a1a and bmp15) in male gonads compared to female gonads. Contrary to expectations, we found highly differential expression of colour pattern (CP) genes in the gonads, suggesting the 'neofunctionalisation' of those genes in the zebrafish reproduction system. However, in the caudal fin, no differential expression of CP genes was identified, suggesting the observed differences in colouration between males and females in adult fish may be due to post-transcriptional regulation of key enzymes involved in pigment synthesis and distribution.

CONCLUSIONS

Our study demonstrates male-biased sex ratio under high temperature condition and support a polygenic SD (PSD) system in laboratory zebrafish. We identify a subset of pathways (tight junction, gap junction and apoptosis), enriched for SD and CP genes, which appear to be co-regulated in the same pathway, providing evidence for involvement of those genes in the regulation of phenotypic sexual dimorphism in zebrafish.

Transcriptome Profile Analysis on Ovarian Tissues of Autotetraploid Fish and Diploid Red Crucian Carp

Polyploidization can significantly alter the size of animal gametes. Autotetraploid fish (RRRR, 4nRR = 200) (4nRR) possessing four sets of chromosomes were derived from whole-genome duplication in red crucian carp (RR, 2n = 100) (RCC). The diploid eggs of the 4nRR fish were significantly larger than the eggs of RCC. To explore the differences between the ovaries of these two ploidies of fishes at the molecular level, we compared the ovary transcriptome profiles of 4nRR fish and RCC and identified differentially expressed genes (DEGs). A total of 19,015 unigenes were differentially expressed between 4nRR fish and RCC, including 12,591 upregulated and 6,424 downregulated unigenes in 4nRR fish. Functional analyses revealed that eight genes (CDKL1, AHCY, ARHGEF3, TGFβ, WNT11, CYP27A, GDF7, and CKB) were involved in the regulation of cell proliferation, cell division, gene transcription, ovary development and energy metabolism, suggesting that these eight genes were related to egg size in 4nRR fish and RCC. We validated the expression levels of these eight DEGs in 4nRR fish and RCC using quantitative PCR. The study results provided insights into the regulatory mechanisms underlying the differences in crucian carp egg sizes.

Comparison of proteomic profiles in the ovary of Sterlet sturgeon (Acipenser ruthenus) during vitellogenic stages

One of the challenges of sturgeon aquaculture is that sturgeon takes an extended amount of time to reach sexual maturity. The pattern of the protein expression in relation to the late maturity of sturgeon can help to better understand changes in sexual maturity. 17β-estradiol (E2), testosterone (T) and vitellogenin (Vtg) levels were examined at all stages of sexual maturation in Sterlet sturgeon (Acipenser ruthenus). Two-dimensional gel electrophoresis and mass spectrometry analysis were used to show the pattern of the ovarian proteins. The T levels increased from the previtellogenic to the postvitellogenic stages (P < 0.05) and Vtg showed a decremental pattern in pre- and postvitellogenic, and atresia (not significantly). The analysis showed 900 protein spots, 19 of which were successfully identified and had significant differences between the previtellogenic and the vitellogenic groups (P < 0.05). Among the identified proteins, 40% involved in cell defense (heat shock protein, Glutathione peroxidase, natural killer enhancing factor, peroxiredoxin-2), 30% in transcription and translation (constitutive photomorphogenesis 9 and Ybx2), 20% in metabolism and energy production (triose-phosphate isomerase (TPI)) and 10% in transport (glycolipid transfer protein). In the vitellogenic stage, the proteins were related to metabolism and energy production (TPI, ES1, creatin kinase, enolase, nucleoside diphosphate kinase, 50%), cell defense (thioredoxin and dislophid isomerase, 20%) and transport (fatty acid binding protein, 10%). Our findings show changes in protein expression pattern from cell defense to metabolism during egg development.

Functional roles of eriocalyxin B in zebrafish revealed by transcriptome analysis

Eriocalyxin B (EriB) is a naturalent-kaurane diterpenoid obtained fromIsodon eriocalyxvar.laxiflora(family Lamiaceae), which exerted multiple biological activities (e.g.anti-tumor and anti-inflammatory)viathe alteration of gene expression and signaling transduction.

Transcriptomic Analysis for Differentially Expressed Genes in Ovarian Follicle Activation in the Zebrafish

In teleosts, the onset of puberty in females is marked by the appearance of the first wave of pre-vitellogenic (PV) follicles from the pool of primary growth (PG) follicles (follicle activation) in the ovary during sexual maturation. To understand the mechanisms underlying follicle activation and therefore puberty onset, we undertook this transcriptomic study to investigate gene expression profiles in the event. Our analysis revealed a total of 2,027 up-regulated and 859 down-regulated genes during the PG-PV transition. Gene Ontology (GO) analysis showed that in addition to basic cellular functions such as gene transcription, cell differentiation, and cell migration, other biological processes such as steroidogenesis, cell signaling and angiogenesis were also enriched in up-regulated genes; by comparison, some processes were down-regulated including piRNA metabolism, gene silencing and proteolysis. Further Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified a variety of signaling pathways that might play pivotal roles in PG-PV transition, including MAPK, TGF-β, Hedgehog, FoxO, VEGF, Jak-STAT, and phosphatidylinositol signaling pathways. Other pathways of particular interest included endocytosis and glycosaminoglycan biosynthesis. We also analyzed expression changes of genes expressed in different compartments viz. oocytes and follicle cells. Interestingly, most oocyte-specific genes remained unchanged in expression during follicle activation whereas a great number of genes specifically expressed in the follicle cells showed significant changes in expression. Overall, this study reported a comprehensive analysis for genes, biological processes and pathways involved in follicle activation, which also marks female puberty onset in the zebrafish when occurring for the first time in sexual maturation.

Proteomic responses to elevated ocean temperature in ovaries of the ascidian Ciona intestinalis

Ciona intestinalis, a common sea squirt, exhibits lower reproductive success at the upper extreme of the water temperatures it experiences in coastal New England. In order to understand the changes in protein expression associated with elevated temperatures, and possible response to global temperature change, we reared C. intestinalis from embryos to adults at 18°C (a temperature at which they reproduce normally at our collection site in Rhode Island) and 22°C (the upper end of the local temperature range). We then dissected ovaries from animals at each temperature, extracted protein, and measured proteomic levels using shotgun mass spectrometry (LC-MS/MS). 1532 proteins were detected at a 1% false discovery rate present in both temperature groups by our LC-MS/MS method. 62 of those proteins are considered up- or down-regulated according to our statistical criteria. Principal component analysis shows a clear distinction in protein expression pattern between the control (18°C) group and high temperature (22°C) group. Similar to previous studies, cytoskeletal and chaperone proteins are upregulated in the high temperature group. Unexpectedly, we find evidence that proteolysis is downregulated at the higher temperature. We propose a working model for the high temperature response in C. intestinalis ovaries whereby increased temperature induces upregulation of signal transduction pathways involving PTPN11 and CrkL, and activating coordinated changes in the proteome especially in large lipid transport proteins, cellular stress responses, cytoskeleton, and downregulation of energy metabolism.

Comparative proteomics analysis of spermary and ovary in Hyriopsis schlegelii

We provide the first large-scale quantitative proteomics analysis in Hyriopsis schlegelii. To investigate the proteins expressed in the gonads, a quantitative proteomics approach has been utilized to analyze differentially expressed proteins between the spermary and ovary. In this study, we identified and quantified 2416 proteins in the gonads of Hyriopsis schlegelii. Of these, 559 proteins showed significantly different expression between the spermary and ovary. Some specific proteins expressed in either the spermary or ovary were identified in Hyriopsis schlegelii. In addition, a series of proteins related to gametogenesis were also identified. Compared with previous reports, many proteins in Hyriopsis schlegelii identified here have different expression patterns between the spermary and ovary. The special hermaphroditism in Hyriopsis schlegelii may contribute to these inconsistent results. The provided proteomics data could be considered as a starting point for subsequent studies focusing on the proteins involved in sexual gland development and maturity.

Heat Shock Proteins and Maternal Contribution to Oogenesis and Early Embryogenesis

Early embryos develop from fertilized eggs using materials that are stored during oocyte growth and which can be defined as maternal contribution (molecules, factors, or determinants). Several heat shock proteins (HSPs) and the heat shock transcriptional factor (HSF) are part of the maternal contribution that is critical for successful embryogenesis and reproduction. A maternal role for heat shock-related genes was mainly demonstrated in genetic experimental organisms (e.g., fly, nematode, mouse). Nowadays, an increasing number of "omics" data are produced from a large panel of organisms implementing a catalog of maternal and/or embryonic HSPs and HSFs. However, for most of them, it remains to better understand their potential roles in this context. Existing and future genome-wide screens mainly set up to create loss-of-function are likely to improve this situation. This chapter will discuss available data from various experimental organisms following the developmental steps from egg production to early embryogenesis.

iTRAQ proteomic analysis of salinity acclimation proteins in the gill of tropical marbled eel (Anguilla marmorata)

Osmoregulation plays an important role in the migration process of catadromous fish. The osmoregulatory mechanisms of tropical marbled eel (Anguilla marmorata), a typical catadromous fish, did not gain sufficient attention, especially at the molecular level. In order to enrich the protein database of A. marmorata, a proteomic analysis has been carried out by iTRAQ technique. Among 1937 identified proteins in gill of marbled eel, the expression of 1560 proteins (80 %) was quantified. Compared with the protein expression level in the gill of marbled eel in freshwater (salinity of 0 ‰), 336 proteins were up-regulated and 67 proteins were down-regulated in seawater (salinity of 25 ‰); 33 proteins were up-regulated and 32 proteins were down-regulated in brackish water (salinity of 10 ‰). These up-regulated proteins including Na(+)/K(+)-ATPase, V-type proton ATPase, sodium-potassium-chloride co-transporter and heat shock protein 90 were enriched in many KEGG-annotated pathways, which are related to different functions of the gill. The up-regulated oxidative phosphorylation and seleno-compound metabolism pathways involve the synthesis and consumption of ATP, which represents extra energy consumption. Another identified pathway is the ribosome pathway in which a large number of up-regulated proteins are involved. It is also more notable that tight junction and cardiac muscle contraction pathways may have correlation with ion transport in gill cells. This is the first report describing the proteome of A. marmorata for acclimating to the change of salinity. These results provide a functional database for migratory fish and point out some possible new interactions on osmoregulation in A. marmorata.

Comprehensive and quantitative proteomic analyses of zebrafish plasma reveals conserved protein profiles between genders and between zebrafish and human

Omic approaches have been increasingly used in the zebrafish model for holistic understanding of molecular events and mechanisms of tissue functions. However, plasma is rarely used for omic profiling because of the technical challenges in collecting sufficient blood. In this study, we employed two mass spectrometric (MS) approaches for a comprehensive characterization of zebrafish plasma proteome, i.e. conventional shotgun liquid chromatography-tandem mass spectrometry (LC-MS/MS) for an overview study and quantitative SWATH (Sequential Window Acquisition of all THeoretical fragment-ion spectra) for comparison between genders. 959 proteins were identified in the shotgun profiling with estimated concentrations spanning almost five orders of magnitudes. Other than the presence of a few highly abundant female egg yolk precursor proteins (vitellogenins), the proteomic profiles of male and female plasmas were very similar in both number and abundance and there were basically no other highly gender-biased proteins. The types of plasma proteins based on IPA (Ingenuity Pathway Analysis) classification and tissue sources of production were also very similar. Furthermore, the zebrafish plasma proteome shares significant similarities with human plasma proteome, in particular in top abundant proteins including apolipoproteins and complements. Thus, the current study provided a valuable dataset for future evaluation of plasma proteins in zebrafish.

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.

Transcriptome Profile Analysis of Ovarian Tissues from Diploid and Tetraploid Loaches Misgurnus anguillicaudatus

RNA sequencing and short-read assembly was utilized to produce a transcriptome of ovarian tissues from three-year-old diploid and tetraploid loaches (Misgurnus anguillicaudatus). A total of 28,369 unigenes were obtained, comprising 10,546 unigenes with length longer than 1000 bp. More than 73% of the unigenes were annotated through sequence comparison with databases. The RNA-seq data revealed that 2253 genes were differentially expressed between diploid and tetraploid loaches, including 1263 up-regulated and 990 down-regulated genes in tetraploid loach. Some differentially expressed genes, such as vitellogenin (Vtg), gonadotropin releasing hormone receptor type A (GnRHRA), steroidogenic acute regulatory protein (StAR), mitogen-activated protein kinase 14a (MAPK14a), ATP synthase subunit alpha (atp5a), and synaptonemal complex protein 1 (Scp1), were involved in regulation of cell proliferation, division, gene transcription, ovarian development and energy metabolism, suggesting that these genes were related to egg diameter of the loach. Results of transcriptome profiling here were validated using real time quantitative PCR in ten selected genes. The present study provided insights into the transcriptome profile of ovarian tissues from diploid and tetraploid loaches Misgurnus anguillicaudatus, which was made available to the research community for functional genomics, comparative genomics, polyploidy evolution and molecular breeding of this loach and other related species.

Whole genome data for omics-based research on the self-fertilizing fish Kryptolebias marmoratus

Genome resources have advantages for understanding diverse areas such as biological patterns and functioning of organisms. Omics platforms are useful approaches for the study of organs and organisms. These approaches can be powerful screening tools for whole genome, proteome, and metabolome profiling, and can be used to understand molecular changes in response to internal and external stimuli. This methodology has been applied successfully in freshwater model fish such as the zebrafish Danio rerio and the Japanese medaka Oryzias latipes in research areas such as basic physiology, developmental biology, genetics, and environmental biology. However, information is still scarce about model fish that inhabit brackish water or seawater. To develop the self-fertilizing killifish Kryptolebias marmoratus as a potential model species with unique characteristics and research merits, we obtained genomic information about K. marmoratus. We address ways to use these data for genome-based molecular mechanistic studies. We review the current state of genome information on K. marmoratus to initiate omics approaches. We evaluate the potential applications of integrated omics platforms for future studies in environmental science, developmental biology, and biomedical research. We conclude that information about the K. marmoratus genome will provide a better understanding of the molecular functions of genes, proteins, and metabolites that are involved in the biological functions of this species. Omics platforms, particularly combined technologies that make effective use of bioinformatics, will provide powerful tools for hypothesis-driven investigations and discovery-driven discussions on diverse aspects of this species and on fish and vertebrates in general.

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.

RNA-Seq technology and its application in fish transcriptomics

High-throughput sequencing technologies, also known as next-generation sequencing (NGS) technologies, have revolutionized the way that genomic research is advancing. In addition to the static genome, these state-of-art technologies have been recently exploited to analyze the dynamic transcriptome, and the resulting technology is termed RNA sequencing (RNA-seq). RNA-seq is free from many limitations of other transcriptomic approaches, such as microarray and tag-based sequencing method. Although RNA-seq has only been available for a short time, studies using this method have completely changed our perspective of the breadth and depth of eukaryotic transcriptomes. In terms of the transcriptomics of teleost fishes, both model and non-model species have benefited from the RNA-seq approach and have undergone tremendous advances in the past several years. RNA-seq has helped not only in mapping and annotating fish transcriptome but also in our understanding of many biological processes in fish, such as development, adaptive evolution, host immune response, and stress response. In this review, we first provide an overview of each step of RNA-seq from library construction to the bioinformatic analysis of the data. We then summarize and discuss the recent biological insights obtained from the RNA-seq studies in a variety of fish species.

Comparative proteomic profiling during ovarian development of the shrimp Metapenaeus ensis

Two-dimensional electrophoresis and mass spectrometry were used to identify proteins that are differentially expressed during ovarian maturation in Metapenaeus ensis. 87 spots with consistently significant quantitative differences (≥ 1.5-fold for vol%) among stage I, III and V ovaries were chosen for MS/MS analysis. 45 spots were significantly matched to known proteins in the database (Mascot score >40). Half of them were down-regulated, in contrast to 9 out of 45 proteins that were up-regulated as ovarian maturation proceeded. Functionally, these identified proteins could be classified into five major groups, including cytoskeleton (11 %), metabolism (18 %), signal transduction (32 %), gene expression (14 %) and immune response (7 %). Among the differentially expressed reproduction-related proteins, the mRNA expression level of cellular retinoic acid/retinol binding protein in M. ensis (MeCRABP) during ovarian maturation was further characterized by quantitative real-time PCR. It was down-regulated during ovarian maturation. In situ hybridization further revealed that MeCRABP transcript was localized in ooplasm of previtellogenic oocytes but not in vitellogenic oocytes. These results demonstrate the application of proteomic analysis for identification of proteins involved in shrimp ovarian maturation and they provide new insights into ovarian development.

Transcriptome of Atlantic cod (Gadus morhua L.) early embryos from farmed and wild broodstocks

Significant efforts have been made to elucidate factors affecting egg quality in fish. Recently, we have shown that eggs originating from wild broodstock (WB) of Atlantic cod (Gadus morhua L.) are of superior quality to those derived from farmed broodstock (FB), and this is associated with differences in the chemical composition of egg yolk. However, maternal transcripts, accumulated during oogenesis, have not been studied extensively in fish. The aim of the present study was to characterize putative maternal mRNA transcriptome in fertilized eggs of Atlantic cod and to compare transcript pools between WB and FB in order to investigate the relation between egg developmental potential and putative maternal mRNA deposits. We performed high-throughput 454 pyrosequencing. For each WB and FB group, five cDNA libraries were individually tagged and sequenced, resulting in 98,687 (WB) and 119,333 (FB) average reads per library. Sequencing reads were de novo assembled, annotated, and mapped. Out of 13,726 identified isotigs, 238 were differentially expressed between WB and FB, with 155 isotigs significantly upregulated in WB. The sequence reads were mapped to 11,340 different Atlantic cod transcripts and 158 sequences were differentially expressed between the 2 groups. Important transcripts involved in fructose metabolism, fatty acid metabolism, glycerophospholipid metabolism, and oxidative phosphorylation were differentially represented between the two broodstock groups, showing potential as biomarkers of egg quality in teleosts. Our findings contribute to the hypothesis that maternal mRNAs affect egg quality and, consequently, the early development of fish.

Analysis of protein expression in zebrafish during gonad differentiation by targeted proteomics

The molecular mechanisms governing sex determination and differentiation in the zebrafish (Danio rerio) are not fully understood. To gain more insights into the function of specific genes in these complex processes, the expression of multiple candidates needs to be assessed, preferably on the protein level. Here, we developed a targeted proteomics method based on selected reaction monitoring (SRM) to study the candidate sex-related proteins in zebrafish which were selected based on a global proteomics analysis of adult gonads and representational difference analysis of male and female DNA, as well as on published information on zebrafish and other vertebrates. We employed the developed SRM protocols to acquire time-resolved protein expression profiles during the gonad differentiation period in vas::EGFP transgenic zebrafish. Evidence on protein expression was obtained for the first time for several candidate genes previously studied only on the mRNA level or suggested by bioinformatic predictions. Tuba1b (tubulin alpha 1b), initially included in the study as one of the potential housekeeping proteins, was found to be preferentially expressed in the adult testis with nearly absent expression in the ovary. The revealed changes in protein expression patterns associated with gonad differentiation suggest that several of the examined proteins, especially Ilf2 and Ilf3 (interleukin enhancer-binding factors 2 and 3), Raldh3 (retinaldehyde dehydrogenase type 3), Zgc:195027 (low density lipoprotein-related receptor protein 3) and Sept5a (septin 5a), may play a specific role in the sexual differentiation in zebrafish.

Dynamics of the striped bass (Morone saxatilis) ovary proteome reveal a complex network of the translasome

We evaluated changes in the striped bass (Morone saxatilis) ovary proteome during the annual reproductive cycle using label-free quantitative mass spectrometry and a novel machine learning analysis based on K-means clustering and support vector machines. Modulated modularity clustering was used to group co-variable proteins into expression modules and Gene Ontology (GO) biological process and KEGG pathway enrichment analyses were conducted for proteins within those modules. We discovered that components of the ribosome along with translation initiation and elongation factors generally decrease as the annual ovarian cycle progresses toward ovulation, concomitant with a slight increase in components of the 26S-proteasome. Co-variation within more than one expression module of components from these two multi-protein complexes suggests that they are not only co-regulated, but that co-regulation occurs through more than one sub-network. These components also co-vary with subunits of the TCP-1 chaperonin system and enzymes of intermediary metabolic pathways, suggesting that protein folding and cellular bioenergetic state play important roles in protein synthesis and degradation. We provide further evidence to suggest that protein synthesis and degradation are intimately linked, and our results support function of a proteasome-ribosome supercomplex known as the translasome.

Transcriptomic analyses of sexual dimorphism of the zebrafish liver and the effect of sex hormones

The liver is one of the most sex-dimorphic organs in both oviparous and viviparous animals. In order to understand the molecular basis of the difference between male and female livers, high-throughput RNA-SAGE (serial analysis of gene expression) sequencing was performed for zebrafish livers of both sexes and their transcriptomes were compared. Both sexes had abundantly expressed genes involved in translation, coagulation and lipid metabolism, consistent with the general function of the liver. For sex-biased transcripts, from in addition to the high enrichment of vitellogenin transcripts in spawning female livers, which constituted nearly 80% of total mRNA, it is apparent that the female-biased genes were mostly involved in ribosome/translation, estrogen pathway, lipid transport, etc, while the male-biased genes were enriched for oxidation reduction, carbohydrate metabolism, coagulation, protein transport and localization, etc. Sexual dimorphism on xenobiotic metabolism and anti-oxidation was also noted and it is likely that retinol x receptor (RXR) and liver x receptor (LXR) play central roles in regulating the sexual differences of lipid and cholesterol metabolisms. Consistent with high ribosomal/translational activities in the female liver, female-biased genes were significantly regulated by two important transcription factors, Myc and Mycn. In contrast, Male livers showed activation of transcription factors Ppargc1b, Hnf4a, and Stat4, which regulate lipid and glucose metabolisms and various cellular activities. The transcriptomic responses to sex hormones, 17β-estradiol (E2) or 11-keto testosterone (KT11), were also investigated in both male and female livers and we found that female livers were relatively insensitive to sex hormone disturbance, while the male livers were readily affected. E2 feminized male liver by up-regulating female-biased transcripts and down-regulating male-biased transcripts. The information obtained in this study provides comprehensive insights into the sexual dimorphism of zebrafish liver transcriptome and will facilitate further development of the zebrafish as a human liver disease model.

Proteome analysis of the fathead minnow (Pimephales promelas) reproductive testes

Proteomics is becoming more widely used as a tool in fish physiology and toxicology and can offer mechanistic insight into organism responses to environmental signals and stressors. Using a LTQ Orbitrap Velos MS/MS, we detected 1075 proteins in the reproductive testis of fathead minnow. Proteins localized to the testis included those with a role in spermatogenesis, DNA repair, gamete meiosis, and proteins that have methylation and phosporylation activity, which are important regulatory mechanisms required for sperm maturation. Enrichment analysis revealed that proteins involved in translation, excision DNA repair, and chromatin remodeling were significantly enriched in the testis (>25% protein coverage of the cellular pathways). Proteins involved in RNA metabolism, spliceosome assembly, metabolism, and DNA unwinding were localized to the testis, and the DEAD (Asp-Glu-Ala-Asp) box RNA-dependent helicase family was well represented in this reproductive tissue. Based upon common detected proteins and functional processes between FHMs and the more ancient sharks, other ray-finned fishes, and mammals, we hypothesize that biological processes involved in the testis (DNA unwinding, RNA processing, spliceosome assembly) have been conserved throughout vertebrate evolution. This study provides the foundation for more in depth proteomics studies investigating the effects of hormones and endocrine disruptors in the teleostean testes.

An ovary transcriptome for all maturational stages of the striped bass (Morone saxatilis), a highly advanced perciform fish

Background

The striped bass and its relatives (genus Morone) are important fisheries and aquaculture species native to estuaries and rivers of the Atlantic coast and Gulf of Mexico in North America. To open avenues of gene expression research on reproduction and breeding of striped bass, we generated a collection of expressed sequence tags (ESTs) from a complementary DNA (cDNA) library representative of their ovarian transcriptome.

Results

Sequences of a total of 230,151 ESTs (51,259,448 bp) were acquired by Roche 454 pyrosequencing of cDNA pooled from ovarian tissues obtained at all stages of oocyte growth, at ovulation (eggs), and during preovulatory atresia. Quality filtering of ESTs allowed assembly of 11,208 high-quality contigs ≥ 100 bp, including 2,984 contigs 500 bp or longer (average length 895 bp). Blastx comparisons revealed 5,482 gene orthologues (E-value < 10^-3), of which 4,120 (36.7% of total contigs) were annotated with Gene Ontology terms (E-value < 10^-6). There were 5,726 remaining unknown unique sequences (51.1% of total contigs). All of the high-quality EST sequences are available in the National Center for Biotechnology Information (NCBI) Short Read Archive (GenBank: SRX007394). Informative contigs were considered to be abundant if they were assembled from groups of ESTs comprising ≥ 0.15% of the total short read sequences (≥ 345 reads/contig). Approximately 52.5% of these abundant contigs were predicted to have predominant ovary expression through digital differential display in silico comparisons to zebrafish (Danio rerio) UniGene orthologues. Over 1,300 Gene Ontology terms from Biological Process classes of Reproduction, Reproductive process, and Developmental process were assigned to this collection of annotated contigs.

Conclusions

This first large reference sequence database available for the ecologically and economically important temperate basses (genus Morone) provides a foundation for gene expression studies in these species. The predicted predominance of ovary gene expression and assignment of directly relevant Gene Ontology classes suggests a powerful utility of this dataset for analysis of ovarian gene expression related to fundamental questions of oogenesis. Additionally, a high definition Agilent 60-mer oligo ovary 'UniClone' microarray with 8 × 15,000 probe format has been designed based on this striped bass transcriptome (eArray Group: Striper Group, Design ID: 029004).

DIGE and iTRAQ as biomarker discovery tools in aquatic toxicology

Molecular approaches in ecotoxicology have greatly enhanced mechanistic understanding of the impact of aquatic pollutants in organisms. These methods have included high throughput Omics technologies, including quantitative proteomics methods such as 2D differential in-gel electrophoresis (DIGE) and isobaric tagging for relative and absolute quantitation (iTRAQ). These methods are becoming more widely used in ecotoxicology studies to identify and characterize protein bioindicators of adverse effect. In teleost fish, iTRAQ has been used successfully in different fish species (e.g. fathead minnow, goldfish, largemouth bass) and tissues (e.g. hypothalamus and liver) to quantify relative protein abundance. Of interest for ecotoxicology is that many proteins commonly utilized as bioindicators of toxicity or stress are quantifiable using iTRAQ on a larger scale, providing a global baseline of biological effect from which to assess changes in the proteome. This review highlights the successes to date for high throughput quantitative proteomics using DIGE and iTRAQ in aquatic toxicology. Current challenges for the iTRAQ method for biomarker discovery in fish are the high cost and the lack of complete annotated genomes for teleosts. However, the use of protein homology from teleost fishes in protein databases and the introduction of hybrid LTQ-FT (Linear ion trap-Fourier transform) mass spectrometers with high resolution, increased sensitivity, and high mass accuracy are able to improve significantly the protein identification rates. Despite these challenges, initial studies utilizing iTRAQ for ecotoxicoproteomics have exceeded expectations and it is anticipated that the use of non-gel based quantitative proteomics will increase for protein biomarker discovery and for characterization of chemical mode of action.

Global proteomics analysis of testis and ovary in adult zebrafish (Danio rerio)

The molecular mechanisms controlling sex determination and differentiation in zebrafish (Danio rerio) are largely unknown. A genome-wide analysis may provide comprehensive insights into the processes involved. The mRNA expression in zebrafish gonads has been fairly well studied, but much less data on the corresponding protein expression are available, although the proteins are considered to be more relevant markers of gene function. Because mRNA and protein abundances rarely correlate well, mRNA profiles need to be complemented with the information on protein expression. The work presented here analyzed the proteomes of adult zebrafish gonads by a multidimensional protein identification technology, generating the to-date most populated lists of proteins expressed in mature zebrafish gonads. The acquired proteomics data partially confirmed existing transcriptomics information for several genes, including several novel transcripts. However, disagreements between mRNA and protein abundances were often observed, further stressing the necessity to assess the expression on different levels before drawing conclusions on a certain gene's expression and function. Several gene groups expressed in a sexually dimorphic way in zebrafish gonads were identified. Their potential importance for gonad development and function is discussed. The data gained in the current study provide a basis for further work on elucidating processes occurring during zebrafish development with use of high-throughput proteomics.

Knock-down of cathepsin D affects the retinal pigment epithelium, impairs swim-bladder ontogenesis and causes premature death in zebrafish

The lysosomal aspartic protease Cathepsin D (CD) is ubiquitously expressed in eukaryotic organisms. CD activity is essential to accomplish the acid-dependent extensive or partial proteolysis of protein substrates within endosomal and lysosomal compartments therein delivered via endocytosis, phagocytosis or autophagocytosis. CD may also act at physiological pH on small-size substrates in the cytosol and in the extracellular milieu. Mouse and fruit fly CD knock-out models have highlighted the multi-pathophysiological roles of CD in tissue homeostasis and organ development. Here we report the first phenotypic description of the lack of CD expression during zebrafish (Danio rerio) development obtained by morpholino-mediated knock-down of CD mRNA. Since the un-fertilized eggs were shown to be supplied with maternal CD mRNA, only a morpholino targeting a sequence containing the starting ATG codon was effective. The main phenotypic alterations produced by CD knock-down in zebrafish were: 1. abnormal development of the eye and of retinal pigment epithelium; 2. absence of the swim-bladder; 3. skin hyper-pigmentation; 4. reduced growth and premature death. Rescue experiments confirmed the involvement of CD in the developmental processes leading to these phenotypic alterations. Our findings add to the list of CD functions in organ development and patho-physiology in vertebrates.

The zebrafish transcriptome during early development

Background

The transition from fertilized egg to embryo is accompanied by a multitude of changes in gene expression, and the transcriptional events that underlie these processes have not yet been fully characterized. In this study RNA-Seq is used to compare the transcription profiles of four early developmental stages in zebrafish (Danio rerio) on a global scale.

Results

An average of 79 M total reads were detected from the different stages. Out of the total number of reads 65% - 73% reads were successfully mapped and 36% - 44% out of those were uniquely mapped. The total number of detected unique gene transcripts was 11187, of which 10096 were present at 1-cell stage. The largest number of common transcripts was observed between 1-cell stage and 16-cell stage. An enrichment of gene transcripts with molecular functions of DNA binding, protein folding and processing as well as metal ion binding was observed with progression of development. The sequence data (accession number ERP000635) is available at the European Nucleotide Archive.

Conclusion

Clustering of expression profiles shows that a majority of the detected gene transcripts are present at steady levels, and thus a minority of the gene transcripts clusters as increasing or decreasing in expression over the four investigated developmental stages. The three earliest developmental stages were similar when comparing highly expressed genes, whereas the 50% epiboly stage differed from the other three stages in the identity of highly expressed genes, number of uniquely expressed genes and enrichment of GO molecular functions. Taken together, these observations indicate a major transition in gene regulation and transcriptional activity taking place between the 512-cell and 50% epiboly stages, in accordance with previous studies.

Proteins in ecotoxicology - how, why and why not?

The growing interest in the application of proteomic technologies to solve toxicology issues and its relevance in ecotoxicology research has resulted in the emergence of "ecotoxicoproteomics". There is a general consensus that ecotoxicoproteomics is a powerful tool to spot early molecular events involved in toxicant responses, which are responsible for the adverse effects observed at higher levels of biological organization, thus contributing to elucidate the mode of action of stressors and to identify specific biomarkers. Ultimately, early-warning indicators can then be developed and deployed in "in situ" bioassays and in environmental risk assessment. The number of field experiments or laboratory trials using ecologically relevant test-species and involving proteomics has been, until recently, insufficient to allow a critical analysis of the real benefits of the application of this approach to ecotoxicology. This article intends to present an overview on the applications of proteomics in the context of ecotoxicology, focusing mainly on the prospective research to be done in invertebrates. Although these represent around 95% of all animal species and in spite of the key structural and functional roles they play in ecosystems, proteomic research in invertebrates is still in an incipient stage. We will review applications of ecotoxicoproteomics by evaluating the technical methods employed, the organisms and the contexts studied, the advances achieved until now and lastly the limitations yet to overcome will be discussed.

Egg and sperm quality in fish

Fish egg quality can be defined as the ability of the egg to be fertilized and subsequently develop into a normal embryo. Similarly, sperm quality can be defined as its ability to successfully fertilize an egg and subsequently allow the development of a normal embryo. In the wild or under aquaculture conditions, the quality of fish gametes can be highly variable and is under the influence of a significant number of external factors or broodstock management practices. For these reasons, the topic of gamete quality has received increasing attention. Despite the significant efforts made towards a better understanding of the factors involved in the control of gamete quality, the picture is far from being complete and the control of gamete quality remains an issue in the aquaculture industry. Some of the factors responsible for the observed variability of gamete quality remain largely unknown or poorly understood. In addition very little is known about the cellular and molecular mechanisms involved in the control of egg and sperm quality. In the present review, the molecular and cellular characteristics of fish gametes are presented with a special interest for the mechanisms that could participate in the regulation of gamete quality. Then, after defining egg and sperm quality, and how can it can be accurately estimated or predicted, we provide an overview of the main factors that can impact gamete quality in teleosts.

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.

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.

Comparative proteomic and regulatory network analyses of the elongating pig conceptus

Embryo loss during peri-implantation can approach 20% in swine following artificial insemination or natural mating and coincides with rapid conceptus elongation. The objective of the present study was to establish a comprehensive profile of the abundant proteins of the pig conceptus at the time prior to implantation and identify stage-specific changes during elongation. The abundant proteins of a homogenous population of gestational day-11 ovoid (0.7-1 cm) and gestational day-12 filamentous (15-20 cm) porcine concepti were compared by extracting proteins from three independent conceptus pools and separating the proteins by 2-DE. Proteins in 305 spots were analyzed by MALDI-TOF or additionally by LC-MS/MS and 275 were positively identified representing 174 distinct proteins. The proteins could be classified into the following functional categories: cell proliferation/differentiation, cytoskeleton, metabolism, and stress response. Based on spot density, 35 proteins associated with cell proliferation, differentiation, apoptosis, and embryo/maternal signaling, were found to be differentially expressed between ovoid and filamentous concepti. A comparison of the protein expression profile with transcriptomic data from pig concepti of the same developmental stages identified similarities and dissimilarities between protein and mRNA expression profiles. This proteomic study helps to elucidate the biological mechanisms underlying the early embryonic development of the pig.

Identifying genes related to choriogenesis in insect panoistic ovaries by Suppression Subtractive Hybridization

Background

Insect ovarioles are classified into two categories: panoistic and meroistic, the later having apparently evolved from an ancestral panoistic type. Molecular data on oogenesis is practically restricted to meroistic ovaries. If we aim at studying the evolutionary transition from panoistic to meroistic, data on panoistic ovaries should be gathered. To this end, we planned the construction of a Suppression Subtractive Hybridization (SSH) library to identify genes involved in panoistic choriogenesis, using the cockroach Blattella germanica as model.

Results

We constructed a post-vitellogenic ovary library by SSH to isolate genes involved in choriogenesis in B. germanica. The tester library was prepared with an ovary pool from 6- to 7-day-old females, whereas the driver library was prepared with an ovary pool from 3- to 4-day-old females. From the SSH library, we obtained 258 high quality sequences which clustered into 34 unique sequences grouped in 19 contigs and 15 singlets. The sequences were compared against non-redundant NCBI databases using BLAST. We found that 44% of the unique sequences had homologous sequences in known genes of other organisms, whereas 56% had no significant similarity to any of the databases entries. A Gene Ontology analysis was carried out, classifying the 34 sequences into different functional categories. Seven of these gene sequences, representative of different categories and processes, were chosen to perform expression studies during the first gonadotrophic cycle by real-time PCR. Results showed that they were mainly expressed during post-vitellogenesis, which validates the SSH technique. In two of them corresponding to novel genes, we demonstrated that they are specifically expressed in the cytoplasm of follicular cells in basal oocytes at the time of choriogenesis.

Conclusion

The SSH approach has proven to be useful in identifying ovarian genes expressed after vitellogenesis in B. germanica. For most of the genes, functions related to choriogenesis are postulated. The relatively high percentage of novel genes obtained and the practical absence of chorion genes typical of meroistic ovaries suggest that mechanisms regulating chorion formation in panoistic ovaries are significantly different from those of meroistic ones.

Upregulated mRNA expression of desaturase and elongase, two enzymes involved in highly unsaturated fatty acids biosynthesis pathways during follicle maturation in zebrafish

BACKGROUND

Although unsaturated fatty acids such as eicosapentaenoic acid (EPA, C20:5n-3), docosahexaenoic acid (DHA, C22:6n-3) and arachidonic acid (ARA, C20:4n-6), collectively known as the highly unsaturated fatty acids (HUFA), play pivotal roles in vertebrate reproduction, very little is known about their synthesis in the ovary. The zebrafish (Danio rerio) display capability to synthesize all three HUFA via pathways involving desaturation and elongation of two precursors, the linoleic acid (LA, C18:2n-6) and linolenic acid (LNA, C18:3n-3). As a prerequisite to gain full understanding on the importance and regulation of ovarian HUFA synthesis, we described here the mRNA expression pattern of two enzymes; desaturase (fadsd6) and elongase (elovl5), involved in HUFA biosynthesis pathway, in different zebrafish ovarian follicle stages. Concurrently, the fatty acid profile of each follicle stage was also analyzed.

METHODS

mRNA levels of fadsd6 and elovl5 in different ovarian follicle stages were determined by semi-quantitative reverse transcription-polymerase chain reaction (RT-PCR) assays. For analysis of the ovarian follicular fatty acid composition, gas chromatography was used.

RESULTS

Our results have shown that desaturase displayed significant upregulation in expression during the oocyte maturation stage. Expression of elongase was significantly highest in pre-vitellogenic follicles, followed by maturation stage. Fatty acid composition analysis of different ovarian follicle stages also showed that ARA level was significantly highest in pre-vitellogenic and matured follicles. DHA level was highest in both late vitellogenic and maturation stage.

CONCLUSION

Collectively, our findings seem to suggest the existence of a HUFA synthesis system, which could be responsible for the synthesis of HUFA to promote oocyte maturation and possibly ovulation processes. The many advantages of zebrafish as model system to understand folliculogenesis will be useful platform to further elucidate the regulatory and mechanism aspects of ovarian HUFA synthesis.

Genomic resources for a commercial flatfish, the Senegalese sole (Solea senegalensis): EST sequencing, oligo microarray design, and development of the Soleamold bioinformatic platform

Background

The Senegalese sole, Solea senegalensis, is a highly prized flatfish of growing commercial interest for aquaculture in Southern Europe. However, despite the industrial production of Senegalese sole being hampered primarily by lack of information on the physiological mechanisms involved in reproduction, growth and immunity, very limited genomic information is available on this species.

Results

Sequencing of a S. senegalensis multi-tissue normalized cDNA library, from adult tissues (brain, stomach, intestine, liver, ovary, and testis), larval stages (pre-metamorphosis, metamorphosis), juvenile stages (post-metamorphosis, abnormal fish), and undifferentiated gonads, generated 10,185 expressed sequence tags (ESTs). Clones were sequenced from the 3'-end to identify isoform specific sequences. Assembly of the entire EST collection into contigs gave 5,208 unique sequences of which 1,769 (34%) had matches in GenBank, thus showing a low level of redundancy. The sequence of the 5,208 unigenes was used to design and validate an oligonucleotide microarray representing 5,087 unique Senegalese sole transcripts. Finally, a novel interactive bioinformatic platform, Soleamold, was developed for the Senegalese sole EST collection as well as microarray and ISH data.

Conclusion

New genomic resources have been developed for S. senegalensis, an economically important fish in aquaculture, which include a collection of expressed genes, an oligonucleotide microarray, and a publicly available bioinformatic platform that can be used to study gene expression in this species. These resources will help elucidate transcriptional regulation in wild and captive Senegalese sole for optimization of its production under intensive culture conditions.

Over-represented sequences located on 3' UTRs are potentially involved in regulatory functions

Eukaryotic gene expression must be coordinated for the proper functioning of biological processes. This coordination can be achieved both at the transcriptional and post-transcriptional levels. In both cases, regulatory sequences placed at either promoter regions or on UTRs function as markers recognized by regulators that can then activate or repress different groups of genes according to necessity. While regulatory sequences involved in transcription are quite well documented, there is a lack of information on sequence elements involved in post-transcriptional regulation. We used a statistical over-representation method to identify novel regulatory elements located on UTRs. An exhaustive search approach was used to calculate the frequency of all possible n-mers (short nucleotide sequences) in 16,160 human genes of NCBI RefSeq sequences and to identify any peculiar usage of n-mers on UTRs. After a stringent filtering process, we identified 2,772 highly over-represented n-mers on 3' UTRs. We provide evidence that these n-mers are potentially involved in regulatory functions. Identified n-mers overlap with previously identified binding sites for HuR and TIA-1 and, ARE and GRE sequences. We determine also that n-mers overlap with predicted miRNA target sites. Finally, a method to cluster n-mer groups allowed the identification of putative gene networks.