Genomic organization of ZP domain containing egg envelope genes in medaka (Oryzias latipes)

Low-temperature-induced disruption of reproductive axis and sperm vitality via stress axis in Monopterus albus

The ricefield eel (Monopterus albus) is inherently timid and highly sensitive to stress. Our previous studies have shown that low-temperature weather could significantly affect the sperm vitality of ricefield eels. This study aims to investigate the regulatory mechanism of low-temperature effects on testicular function and sperm vitality in ricefield eels. The ricefield eels were initially reared at low (10 °C) and normal (25 °C) temperatures for 24 h. Low temperatures were found to induce the expression of pituitary pro-opiomelanocortin (POMC) and testes insulin-like growth factor-binding protein 1 (IGFBP1) mRNA expression, suggesting that the reduction in sperm vitality could be attributed to the activation of the stress axis. Moreover, the results indicated a significant decrease in sperm occupancy and count in the testes, along with a reduced percentage of motile sperm. Subsequent transcriptome analysis showed substantial inhibition of reproductive hormone genes (gnrh1, lh, and fsh) in the brain and pituitary, and downregulation of meiosis-related genes (dmc1, rec8, and sycp3) in the testes. These findings suggest that low temperatures might disrupt testicular development and spermatogenesis by inhibiting the reproductive axis. Metabolomics analysis then demonstrated a significant reduction in the levels of metabolites related to glycolysis, fatty acid metabolism, and the tricarboxylic acid (TCA) cycle in the testes after low-temperature treatment. Interestingly, the expression of zona pellucida sperm-binding proteins 3 and 4 (ZP3 and ZP4), which may affect sperm vitality and spermatogenesis, was significantly induced by low temperatures in the testes. In conclusion, these findings suggested that low temperatures might affect testicular function and sperm vitality by simultaneously activating the stress axis and inhibiting the reproductive axis and energy metabolism in the testes.

Egg envelope formation of medaka Oryzias latipes requires ZP proteins originating from both the liver and ovary

Teleost oocytes are surrounded by a structure, called chorion or egg envelopes, which is composed of zona pellucida (ZP) proteins. As a result of the gene duplication in teleost, the expression site of the zp genes, coding the major component protein of egg envelopes, changed from the ovary to the maternal liver. In Euteleostei, there are three liver-expressed zp genes, named choriogenin (chg) h, chg hm, and chg l, and the composition of the egg envelope is mostly made up of these Chgs. In addition, ovary-expressed zp genes are also conserved in the medaka genomes, and their proteins have also been found to be minor components of the egg envelopes. However, the specific role of liver-expressed versus ovary-expressed zp genes was unclear. In the present study, we showed that ovary-synthesized ZP proteins first form the base layer of the egg envelope, and then Chgs polymerize inwardly to thicken the egg envelope. To analyze the effects of dysfunction of the chg gene, we generated some chg knockout medaka. All knockout females failed to produce normally fertilized eggs by the natural spawning. The egg envelopes lacking Chgs were significantly thinner, but layers formed by ZP proteins synthesized in the ovary were found in the thin egg envelope of knockout as well as wild-type eggs. These results suggest that the ovary-expressed zp gene is well conserved in all teleosts, including those species in which liver-derived ZP proteins are the major component, because it is essential for the initiation of egg envelope formation.

Comparative transcriptome sequencing analysis of female and male Decapterus macrosoma

Sexual growth dimorphism is a common phenomenon in teleost fish and has led to many reproductive strategies. Growth- and sex-related gene research in teleost fish would broaden our understanding of the process. In this study, transcriptome sequencing of shortfin scad Decapterus macrosoma was performed for the first time, and a high-quality reference transcriptome was constructed. After identification and assembly, a total of 58,475 nonredundant unigenes were obtained with an N50 length of 2,266 bp, and 28,174 unigenes were successfully annotated with multiple public databases. BUSCO analysis determined a level of 92.9% completeness for the assembled transcriptome. Gene expression analysis revealed 2,345 differentially expressed genes (DEGs) in the female and male D. macrosoma, 1,150 of which were female-biased DEGs, and 1,195 unigenes were male-biased DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the DEGs were mainly involved in biological processes including protein synthesis, growth, rhythmic processes, immune defense, and vitellogenesis. Then, we identified many growth- and sex-related genes, including Igf, Fabps, EF-hand family genes, Zp3, Zp4 and Vg. In addition, a total of 19,573 simple sequence repeats (SSRs) were screened and identified from the transcriptome sequences. The results of this study can provide valuable information on growth- and sex-related genes and facilitate further exploration of the molecular mechanism of sexual growth dimorphism.

Prediction of the Effect of Methylation in the Promoter Region of ZP2 Gene on Egg Production in Jinghai Yellow Chickens

Egg production in chickens is a quantitative trait. The aim of this study was to investigate the effect of promoter methylation of the Zona pellucida 2 (ZP2) gene on egg production. Real-time fluorescence quantification showed that the expression of the ZP2 gene in the ovaries of 300-day-old Jinghai yellow chickens in the high-laying group was significantly higher than that in the low-laying group (p < 0.01). A series of deletion fragments of the ZP2 gene promoter in Jinghai yellow chickens had different promoter activities in DF-1 cells, and the core region of the ZP2 gene promoter was found to be between −1552 and −1348. Four CpG islands in the promoter region of the ZP2 gene were detected by software prediction. The overall degree of methylation of the ZP2-1 amplified fragment was negatively correlated with mRNA expression to some extent (R = −0.197); the overall degree of methylation of the ZP2-2 amplified fragment was also negatively correlated with mRNA expression to some extent (R = −0.264), in which the methylation of methylcytosine (mC)-9, mC-20, and mC-21 sites was significantly negatively correlated with mRNA expression (p < 0.05). In addition, the mC-20 and mC-21 sites are located on the Sp1 transcription factor binding site, and it is speculated that these two sites may be the main sites for regulating transcription. In summary, the methylation sites mC-20 and mC-21 of the ZP2 gene may inhibit the binding of Sp1 and DNA, affect the transcription of the ZP2 gene, and then affect the number of eggs produced by the Jinghai yellow chickens.

Medaka, Oryzias latipes, egg envelopes are created by ovarian-expressed ZP proteins and liver-expressed choriogenins

The medaka (Oryzias latipes) egg envelope (chorion) is composed of three major glycoproteins, Zona Interna (ZI)-1, -2, and -3, that originate in the spawning female liver as the precursor proteins Choriogenin (Chg.)H, Chg.Hm, and Chg.L, respectively. These ZI and Chg. proteins contain a structural ZP protein domain that is conserved among the egg envelope proteins of all animals. While ovarian expression of ZP proteins (e.g., ZPCs and ZPB) has been reported in medakas, the functions of these proteins remain unknown. Thus, the present study aimed to determine whether the ovary-expressed medaka ZP protein, mZPC5, is involved in forming the chorion matrix.The mZPC5 gene (mzpc5) was expressed in the ovaries but not the livers of mature female medakas, as shown by reverse transcription-polymerase chain reaction assays with mzpc5-specific primers. In situ hybridization analysis revealed that ovarian mzpc5 expression was restricted to the ooplasm of early (stage I-III) previtellogenic oocytes, and its expression signal weakened with oocyte growth. Following sodium-dodecyl-sulfate polyacrylamide gel electrophoresis (SDS-PAGE) and Western blotting analysis with anti-mZPC5 antibodies, two immunoreactive proteins were detected in the ovary and chorion extracts. These proteins were approximately 50 and 74 kDa in size, like ZI-3 and ZI-2, respectively.Immunohistochemical assays using anti-mZPC5 and anti-Chg.H antibodies localized the mZPC5 protein in the ooplasm of early previtellogenic oocytes. With oocyte growth, mZPC5 tended to accumulate in the chorion, co-localizing with Chg.H.We previously showed that ovary-expressed ZP proteins could not compensate for Chg.L function loss in gene knock-out (chg.l -/-) medakas. As in our previous study, the chg.l-/- females produced oocytes with thin chorions, resulting in infertile soft eggs. However, in the present study, mZPC5 and Chg.H were co-localized in the chg.l-/- chorions. These results suggested that in the medaka previtellogenic oocyte, 1) mZPC5 is secreted from the ooplasm and deposited on the outer surface of its plasma membrane, creating the thin chorion layer; and 2) following the accumulation of liver-derived Chgs., the 3D structure of the chorion matrix is formed cooperatively with mZPC5 and Chgs. during oogenesis. More research is needed to confirm the functions of mZPC5 in chorion structure and physiology.

Lineage-specific evolution of zona pellucida genes in fish

The zona pellucida (ZP) protein constitutes the egg envelope, which surrounds the vertebrate embryo. We performed a comprehensive study on the molecular evolution of ZP genes in Teleostei by cloning and analyzing the expression of ZP genes in fish of Anguilliformes in Elopomorpha, Osteoglossiformes in Osteoglossomorpha, and Clupeiformes in Otocephala to cover unsurveyed fish groups in Teleostei. The present results confirmed findings from our previous reports that the principal organ of ZP gene expression changed from ovary to liver in the common ancestors of Clupeocephala. Even fish species that synthesize egg envelopes in the liver carry the ovary-expressed ZP proteins as minor egg envelope components that were produced by gene duplication during the early stage of Teleostei evolution. The amino acid repeat sequences located at the N-terminal region of ZP proteins are known to be the substrates of transglutaminase responsible for egg envelope hardening and hatching. A repeat sequence was found in zona pellucida Cs of phylogenetically early diverged fish. After changing the synthesis organ, its role is inherited by the N-terminal Pro-Gln-Xaa repeat sequence in liver-expressed zona pellucida B genes of Clupeocephala. These results suggest that teleost ZP genes have independently evolved to maintain fish-specific functions, such as egg envelope hardening and egg envelope digestion, at hatching.

Targeted deletion of liver-expressed Choriogenin L results in the production of soft eggs and infertility in medaka, Oryzias latipes

Egg envelopes (chorions) in medaka, Oryzias latipes, are composed of three major glycoproteins: ZI-1, - 2, and - 3. These gene-encoded chorion glycoproteins are expressed in the liver and/or ovarian oocytes of sexually mature female fish. In medaka, the glycoproteins produced in the female liver are induced by estrogen as Choriogenin (Chg.) H and Chg. H minor (m), which correspond to the zona pellucida (ZP) B (ZPB) protein in mammals, and Chg. L, which corresponds to ZPC in mammals. Chg. H, Chg. Hm, and Chg. L, are then converted to ZI-1, - 2, and - 3, respectively, during oogenesis in medaka ovaries.In the present study, we established a medaka line in which the chg.l gene was inactivated using the transcription activator-like effector nuclease (TALEN) technique. Neither intact chg.l transcripts nor Chg. L proteins were detected in livers of sexually mature female homozygotes for the mutation (homozygous chg.l knockout: chg.l-/-). The chg.l-/- females spawned string-like materials containing "smashed eggs." Closer examination revealed the oocytes in the ovaries of chg.l-/- females had thin chorions, particularly at the inner layer, despite a normal growth rate. In comparing chorions from normal (chg.l+/+) and chg.l-/- oocytes, the latter exhibited abnormal architecture in the chorion pore canals through which the oocyte microvilli pass. These microvilli mediate the nutritional exchange between the oocyte and surrounding spaces and promote sperm-egg interactions during fertilization. Thus, following in vitro fertilization, no embryos developed in the artificially inseminated oocytes isolated from chg.l-/- ovaries. These results demonstrated that medaka ZI-3 (Chg.L) is the major component of the inner layer of the chorion, as it supports and maintains the oocyte's structural shape, enabling it to withstand the pressures exerted against the chorion during spawning, and is essential for successful fertilization. Therefore, gene products of oocyte-specific ZP genes that may be expressed in medaka oocytes cannot compensate for the loss Chg. L function to produce offspring for this species.

Expression profiles of hepatic vitellogenin and gonadal zona pellucida subtypes in gray mullet (Mugil cephalus) with 17α-ethinylestradiol-induced gonadal abnormality

The subtypes of zona pellucida (zp), primarily expressed in female gonads, are considered novel molecular markers for testis-ova (or intersex), a type of gonadal abnormality caused by environmental estrogens (EEs) in Japanese medaka (Oryzias latipes). However, the association between testis-ova and the expression of gonadal zp subtypes is unclear in other teleost species, particularly in species studied in field surveys. In this study, 17α-ethinylestradiol (EE2) was orally administrated at 4-4000 ng/g body weight (BW)/day for 28 days to gray mullets (Mugil cephalus), and gonadal abnormalities were studied using histological analysis. The expression profiles of gonadal zp subtypes (zpb and zpc5) were analyzed to evaluate their suitability as gonadal abnormality markers by comparing with a hepatic vitellogenin (vtg) subtype (vtgAb). The oral administration of EE2 40 and 400 ng/g BW/day for 28 days induced significant gonadal zpb expression, and the gonads showed moderate abnormality (testis-ova). Conversely, the gonadal zpc5 levels decreased significantly in response to the oral administration of EE2 at 4000 ng/g BW/day for 28 days, and the gonads exhibited severe abnormalities. The hepatic vtgAb levels increased upon EE2 treatment regardless of gonadal abnormality. Therefore, the gonadal zpb levels and hepatic vtgAb levels served as appropriate markers for testis-ova and EE2 presence, respectively. However, the diagnosis of severe gonadal abnormality using gonadal zpc5 was moderately accurate. The findings suggest that the combination of vtgAb, zpb, and zpc5 is a potential marker for gonadal abnormality caused by EE contamination in gray mullet. That said, the potential of zpc5 should be reconsidered to determine if it shows greater accuracy in a larger or more diverse population.

Comparative transcriptomic analysis for identification of candidate sex-related genes and pathways in Crimson seabream (Parargyrops edita)

Teleost fishes display the largest array of sex-determining systems among animals, resulting in various reproductive strategies. Research on sex-related genes in teleosts will broaden our understanding of the process, and provide important insight into the plasticity of the sex determination process in vertebrates in general. Crimson seabream (Parargyrops edita Tanaka, 1916) is one of the most valuable and abundant fish resources throughout Asia. However, little genomic information on P. edita is available. In the present study, the transcriptomes of male and female P. edita were sequenced with RNA-seq technology. A total of 388,683,472 reads were generated from the libraries. After filtering and assembling, a total of 79,775 non redundant unigenes were obtained with an N50 of 2,921 bp. The unigenes were annotated with multiple public databases, including NT (53,556, 67.13%), NR (54,092, 67.81%), Swiss-Prot (45,265, 56.74%), KOG (41,274, 51.74%), KEGG (46,302, 58.04%), and GO (11,056, 13.86%) databases. Comparison of the unigenes of different sexes of P. edita revealed that 11,676 unigenes (9,335 in females, 2,341 in males) were differentially expressed between males and females. Of these, 5,463 were specifically expressed in females, and 1,134 were specifically expressed in males. In addition, the expression levels of ten unigenes were confirmed to validate the transcriptomic data by qRT-PCR. Moreover, 34,473 simple sequence repeats (SSRs) were identified in SSR-containing sequences, and 50 loci were randomly selected for primer development. Of these, 36 loci were successfully amplified, and 19 loci were polymorphic. Finally, our comparative analysis identified many sex-related genes (zps, amh, gsdf, sox4, cyp19a, etc.) and pathways (MAPK signaling pathway, p53 signaling pathway, etc.) of P. edita. This informative transcriptomic analysis provides valuable data to increase genomic resources of P. edita. The results will be useful for clarifying the molecular mechanism of sex determination and for future functional analyses of sex-associated genes.

Insights on the seasonal variations of reproductive features in the Eastern Atlantic Bluefin Tuna

The Atlantic Bluefin Tuna (ABFT, Thunnus thynnus) is one of the most intensely exploited fisheries resources in the world. In spite of the years of studies on ABFT, basic aspects of its reproductive biology remain uncertain. To gain insight regarding the seasonal changes of the reproductive characteristics of the eastern stock of ABFT, blood and tissue samples were collected from mature specimens caught in the Mediterranean basin during the reproductive (May-June) and non-reproductive season (Oct-Nov). Histological analysis of the gonads of May-June samples indicated that there were females which were actively spawning (contained post-ovulatory follicles) and females that were not actively spawning that had previtellogenic and fully vitellogenic oocytes. In males, testis were at early or late stage of spermatogenesis during the reproductive season. In Oct-Nov, ovaries contained mostly previtellogenic oocytes as well as β and α atretic follicles while the testis predominantly contained spermatogonia and few cysts with spermatocytes and spermatozoa. Gonadosomatic index (GSI) in females was highest among the actively spawning individuals while in males GSI was higher in early and late spermatogenic individuals compared to those that were spent. Plasma sex steroids levels varied with the reproductive season. In females, estradiol (E₂), was higher in May-June while testosterone (T) and progesterone (P) did not vary. In males, E₂ and T were higher in May-June while P levels were similar at the two sampling points. Circulating follicle stimulating hormone (FSH) was higher in Oct-Nov than in May-June both in males and females. Vitellogenin (VTG) was detected in plasma from both males and females during the reproductive season with levels in females significantly higher than in males. VTG was undetected in Oct-Nov samples. Since choriogenesis is an important event during follicle growth, the expression of three genes involved in vitelline envelope formation and hardening was measured and results showed significantly higher levels in ovaries in fish caught in May-June with respect to those sampled in Oct-Nov. In addition, a set of genes encoding for ion channels that are responsible for oocyte hydration and buoyancy, as well as sperm viability, were characterized at the two time points, and these were found to be more highly expressed in females during the reproductive season. Finally, the expression level of three mRNAs encoding for different lipid-binding proteins was analyzed with significantly higher levels detected in males, suggesting sex-specific expression. Our findings provide additional information on the reproductive biology of ABFT, particularly on biomarkers for the assessment of the state of maturation of the gonad, highlighting gender-specific signals and seasonal differences.

Expression of zona pellucida 3 gene is regulated by 17α-ethinylestradiol in adult topmouth culter Culter alburnus

Estrogen could lead to abnormal modulation or disruption of physical development, reproduction and sexual behavior in aquatic wildlife, especially in fish. Information on the toxicity of estrogens to native species in that can be used in site-specific risk assessments is scarce. In the present study, one zona pellucida 3 (ZP3) homologue termed CaZP3 was firstly identified from topmouth culter Culter alburnus, following its structural characteristics, tissue distribution and transcriptional modulation to 17α-ethinylestradiol (EE2) exposure were investigated. Meanwhile, vitellogenin (VTG) gene was employed to provide a comparison of the reactive ability to EE2 induction. The CaZP3 characterized with analogical functional domains such as ZP domain, SP, IHP, EHP, 12 cysteine residues, one N-linked glycosylation site and two conserved O-linked glycosylation sites and equal number of eight exons and seven introns with ZP3 counterparts of higher species. CaZP3 mRNA predominantly expressed in ovary, besides, highly expressed in female heart and male muscle and relatively high expressed in testis. CaZP3 has the lower reactive ability to EE2 induction in comparison with VTG, however, CaZP3 transcripts were significantly induced in gonads of both male and female culter by EE2 and could be used as an alternative biomarker to monitor EE2 activity. The present results supplement the database for toxicity of EE2, especially for fish species endemic to China and provide some useful information for the monitoring of EE2 activity in aquatic environment.

Gonadal Transcriptome Analysis of Pacific Abalone Haliotis discus discus: Identification of Genes Involved in Germ Cell Development

Little is known about the molecular mechanisms governing gonadal developmental processes in abalones. Here, we conducted transcriptome analysis of Pacific abalone Haliotis discus discus for gene discovery in the brain, ovary, testis, and unfertilized eggs. Among the annotated unigenes, 48.6% of unigenes were identified by Venn diagram analysis as having universal or tissue-specific expression. Twenty-three genes with gonad-biased gene ontology (GO) terms were first obtained. Secondly, 36 genes were found by screening known gene names related to germ cell development. Finally, 17 genes were obtained by querying the annotated unigene database for zygotically expressed gonadal genes (ovary and testis) and maternally expressed gonadal genes (ovary, testis, and unfertilized eggs) using keywords related to reproduction. To further verify tissue distribution pattern and subcellular localization of these genes, RT-PCR and in situ hybridization were performed using a unigene encoding a germ cell marker, vasa, as control. The results showed that vasa was expressed mainly in the early developmental stages of germ cells in both sexes. One of the candidate genes, vitelline envelope zona pellucida domain protein 12 (ZP12), was expressed in the primordial germ cells of immature gonad and early developmental stages of germ cells of the adult female. The results obtained from the present study suggest that vasa and ZP12 are involved in germ cell development of Pacific abalone and that ZP12 is an especially useful germ cell-specific marker in immature adults. The current gonadal transcriptome profile is an extensive resource for future reproductive molecular biology studies of this species.

Egg Coat Proteins Across Metazoan Evolution

All animal oocytes are surrounded by a glycoproteinaceous egg coat, a specialized extracellular matrix that serves both structural and species-specific roles during fertilization. Egg coat glycoproteins polymerize into the extracellular matrix of the egg coat using a conserved protein-protein interaction module-the zona pellucida (ZP) domain-common to both vertebrates and invertebrates, suggesting that the basic structural features of egg coats have been conserved across hundreds of millions of years of evolution. Egg coat proteins, as with other proteins involved in reproduction, are frequently found to be rapidly evolving. Given that gamete compatibility must be maintained for the fitness of sexually reproducing organisms, this finding is somewhat paradoxical and suggests a role for adaptive diversification in reproductive protein evolution. Here we review the structure and function of metazoan egg coat proteins, with an emphasis on the potential role their evolution has played in the creation and maintenance of species boundaries.

The transcriptomic signature of different sexes in two protogynous hermaphrodites: Insights into the molecular network underlying sex phenotype in fish

Sex differentiation is a puzzling problem in fish due to the variety of reproductive systems and the flexibility of their sex determination mechanisms. The Sparidae, a teleost family, reflects this remarkable diversity of sexual mechanisms found in fish. Our aim was to capture the transcriptomic signature of different sexes in two protogynous hermaphrodite sparids, the common pandora Pagellus erythrinus and the red porgy Pagrus pagrus in order to shed light on the molecular network contributing to either the female or the male phenotype in these organisms. Through RNA sequencing, we investigated sex-specific differences in gene expression in both species' brains and gonads. The analysis revealed common male and female specific genes/pathways between these protogynous fish. Whereas limited sex differences found in the brain indicate a sexually plastic tissue, in contrast, the great amount of sex-biased genes observed in gonads reflects the functional divergence of the transformed tissue to either its male or female character. Α common "crew" of well-known molecular players is acting to preserve either sex identity of the gonad in these fish. Lastly, this study lays the ground for a deeper understanding of the complex process of sex differentiation in two species with an evolutionary significant reproductive system.

Estrogen therapy offsets thermal impairment of vitellogenesis, but not zonagenesis, in maiden spawning female Atlantic salmon (Salmo salar)

In female Atlantic salmon (Salmo salar), exposure to warm summer temperatures causes a reduction in plasma 17β-estradiol (E2), which impairs downstream vitellogenesis and zonagenesis, and reduces egg fertility and embryo survival. The aim of the present study was to determine whether E2-treatment could offset thermal impairment of endocrine function and maintain egg quality in maiden (first-time-spawning) S. salar reared at 22 °C. Treatment with E2 at 22 °C stimulated vitellogenin (vtg) gene expression and subsequent protein synthesis which promoted oocyte growth and increased egg size relative to untreated fish at 14 and 22 °C. However, E2-treatment at 22 °C was not associated with an increase in egg fertility and embryo survival relative to untreated fish at 22 °C, despite the positive effects of E2-treatment on vitellogenesis and oocyte growth. As there was no evidence to suggest that the estrogen receptor alpha expression was suppressed by high temperature, this could be due to the lack of stimulation on zonagenesis by E2-treatment observed at high temperature during oocyte development. Our results demonstrate that treatment with E2 is not able to maintain zonagenesis or egg quality in maiden S. salar at high temperature, even when vtg gene expression, protein synthesis and subsequent oocyte growth is promoted. This implies that the mechanisms regulating zonagenesis, but not vitellogenesis are impaired at elevated temperature in female S. salar broodstock, and highlights the remarkable complexity of thermally induced endocrine disruption in fish.

Estrogen alters gonadal soma-derived factor (Gsdf)/Foxl2 expression levels in the testes associated with testis-ova differentiation in adult medaka, Oryzias latipes

Testis-ova differentiation in sexually mature male medaka (Oryzias latipes) is easily induced by estrogenic chemicals, indicating that spermatogonia persist in sexual bipotentiality, even in mature testes in medaka. By contrast, the effects of estrogen on testicular somatic cells associated with testis-ova differentiation in medaka remain unclear. In this study, we focused on the dynamics of sex-related genes (Gsdf, Dmrt1, and Foxl2) expressed in Sertoli cells in the mature testes of adult medaka during estrogen-induced testis-ova differentiation. When mature male medaka were exposed to estradiol benzoate (EB; 800ng/L), testis-ova first appeared after EB treatment for 14days (observed as the first oocytes of the leptotene-zygotene stage). However, the testis remained structurally unchanged, even after EB treatment for 28days. Although Foxl2 is a female-specific sex gene, EB treatment for 7days induced Foxl2/FOXL2 expression in all Sertoli cell-enclosed spermatogonia before testis-ova first appeared; however, Foxl2 was not detected in somatic cells in control testes. Conversely, Sertoli-cell-specific Gsdf mRNA expression levels significantly decreased after EB treatment for 14days, and no changes were observed in DMRT1 localization following EB treatment, whereas Dmrt1 mRNA levels increased significantly. Furthermore, after EB exposure, FOXl2 and DMRT1 were co-localized in Sertoli cells during testis-ova differentiation, although FOXL2 localization was undetectable in Sertoli-cell-enclosed apoptotic testis-ova, whereas DMRT1 remained localized in Sertoli cells. These results indicated for the first time that based on the expression of female-specific sex genes, feminization of Sertoli cells precedes testis-ova differentiation induced by estrogen in mature testes in medaka; however, complete feminization of Sertoli cells was not induced in this study. Additionally, it is suggested strongly that Foxl2 and Gsdf expression constitute potential molecular markers for evaluating the effects of estrogenic chemicals on testicular somatic cells associated with estrogen-induced testis-ova differentiation in mature male medaka.

Neofunctionalization of zona pellucida proteins enhances freeze-prevention in the eggs of Antarctic notothenioids

The mechanisms by which the eggs of the Antarctic notothenioid fishes avoid freezing are not fully understood. Zona pellucida proteins (ZPs) are constituents of the chorion which forms a protective matrix surrounding the egg. Here we report occurrence of freezing temperature-related gene expansion and acquisition of unusual ice melting-promoting (IMP) activity in a family of Antarctic notothenioid ZPs (AnnotoZPs). Members of AnnotoZPs are shown to bind with ice and non-colligatively depress the melting point of a solution in a range of 0.26 to 0.65 °C at a moderate concentration. Eggs of zebrafishes expressing an AnnotoZP transgene show improved melting point depression and enhanced survival in freezing conditions. Mutational analyses in a representative AnnotoZP indicate the ZP domain and patches of acidic residues are essential structures for the IMP activity. AnnotoZPs, therefore, represent a group of macromolecules that prevent freezing by a unique ZP-ice interaction mechanism distinct from the known antifreeze proteins.

Gene expression profiles in testis of developing male Xenopus laevis damaged by chronic exposure of atrazine

As a widely used herbicide, atrazine (AZ) has been extensively studied for its adverse effects on the reproductive system, especially feminization in male animals. However, the relationship of gene expression changes and associated toxicological endpoints remains unclear. In this study, developing Xenopus laevis tadpoles were exposed to concentration of AZ at 0.1, 1, 10 or 100 μg/L continuously. Compared with froglets in the control group, there were no significant differences in body length, body weight, liver weight and hepatosomatic index (HSI) of males in groups treated with AZ for 90 d. At 100 μg/L AZ treatment caused a significant reduction of gonad weight and gonadosomatic index (GSI) of males (p < 0.01). In addition, AZ at all dose levels caused testicular degeneration, especially in froglets from the groups with 0.1 and 100 μg/L which exhibited U-shaped dose-response trend. We further investigated the gene expression changes associated with the testicular degeneration induced by AZ. We found that the expression of 1165 genes was significantly altered with 616 upregulated and 549 downregulated compared to the expression profile of the control animals. KEGG analysis showed that genes which were significantly affected by AZ are mainly involved in arginine and proline metabolism, cell cycle, riboflavin metabolism, spliceosome, base excision repair and progesterone-mediated oocyte maturation pathway. Our results show that AZ may affect reproductive and immune systems by interference with the related gene expression changes during the male X. laevis development. The findings may help to clarify the feminization mechanisms of AZ in male X. laevis.

Development of a promising fish model (Oryzias melastigma) for assessing multiple responses to stresses in the marine environment

With the increasing number of contaminants in the marine environment, various experimental organisms have been “taken into labs” by investigators to find the most suitable environmentally relevant models for toxicity testing. The marine medaka, Oryzias melastigma, has a number of advantages that make it a prime candidate for these tests. Recently, many studies have been conducted on marine medaka, especially in terms of their physiological, biochemical, and molecular responses after exposure to contaminants and other environmental stressors. This review provides a literature survey highlighting the steady increase of ecotoxicological research on marine medaka, summarizes the advantages of using O. melastigma as a tool for toxicological research, and promotes the utilization of this organism in future studies.

Molecular co-evolution of a protease and its substrate elucidated by analysis of the activity of predicted ancestral hatching enzyme

Background

Hatching enzyme is a protease that digests the egg envelope, enabling hatching of the embryo. We have comprehensively studied the molecular mechanisms of the enzyme action to its substrate egg envelope, and determined the gene/protein structure and phylogenetic relationships. Because the hatching enzyme must have evolved while maintaining its ability to digest the egg envelope, the hatching enzyme-egg envelope protein pair is a good model for studying molecular co-evolution of a protease and its substrate.

Results

Hatching enzymes from medaka (Oryzias latipes) and killifish (Fundulus heteroclitus) showed species-specific egg envelope digestion. We found that by introducing four medaka-type residue amino acid substitutions into recombinant killifish hatching enzyme, the mutant killifish hatching enzyme could digest medaka egg envelope. Further, we studied the participation of the cleavage site of the substrate in the species-specificity of hatching enzyme. A P2-site single amino acid substitution was responsible for the species-specificity. Estimation of the activity of the predicted ancestral enzymes towards various types of cleavage sites along with prediction of the evolutionary timing of substitutions allowed prediction of a possible evolutionary pathway, as follows: ancestral hatching enzyme, which had relatively strict substrate specificity, developed broader specificity as a result of four amino acid substitutions in the active site cleft of the enzyme. Subsequently, a single substitution occurred within the cleavage site of the substrate, and the recent feature of species-specificity was established in the hatching enzyme-egg envelope system.

Conclusions

The present study clearly provides an ideal model for protease-substrate co-evolution. The evolutionary process giving rise to species-specific egg envelope digestion of hatching enzyme was initiated by amino acid substitutions in the enzyme, resulting in altered substrate specificity, which later allowed an amino acid substitution in the substrate.

A specific genomic organization and a novel promoter sequence for both ZP2 and ZP3 gene expressions in the Pingxiang red transparent crucian carp, Carassius auratus var. pingxiangnensis

We cloned the full-length cDNA of ZP2 from Carassius auratus var. Pingxiangnensis (CaP_ZP2) and identified a cluster of three ZP genes through its DNA walker. These three genes, CaP_ZP3.1, CaP_ZP2 and CaP_ZP3.2 were located within a 10,855bp region and each comprised of eight exons spanning 1348bp, 1638bp and 1348bp, respectively. Protein bands of egg membrane between 40 and 70kDa were in concordance with the deduced amino acid of these three CaP_ZP genes cDNA and with their molecular mass. This is the first report that two CaP_ZP3 genes were separated by CaP_ZP2 gene. A novel sequence of 1097bp, located between CaP_ZP2 and CaP_ZP3.2, was inserted into the modified pAcGFP1-1 vector in the forward and reverse directions. Results showed that individual sequence served as promoters utilizing common regulatory elements in the forward and reverse directions for both CaP_ZP2 and CaP_ZP3.2 gene expressions. In situ hybridization against CaP_ZP2 confirmed that a strong positive signal was detected in the early development oocytes. Similarly, real-time PCR results also showed that CaP_ZP2 transcription increased mainly in a 4-8month ovary, but was decreased dramatically in a 9-12month ovary.

Expression of two zona pellucida genes is regulated by 17α-ethinylestradiol in adult rare minnow Gobiocypris rarus

Zona pellucida (ZP) proteins are glycoproteins synthesized in liver, ovary or in both tissues in fish. In the present study, we aimed to determine the responsiveness of ZP2 and ZP3 to 17α-ethinylestradiol (EE2) in adult rare minnow Gobiocypris rarus. The full length of ZP3 cDNA was firstly characterized and its tissue distribution revealed that ZP3 mRNA was predominantly expressed in ovary of G. rarus. The gene expression profiles of ZP2, ZP3 and vitellogenin (VTG) were analyzed in gonad and liver of adult G. rarus exposed to EE2 at 1, 5, 25, and 125 ng/L for 3 and 6 days. The results show that ZP2 is more sensitive than ZP3 in gonads of both genders, and VTG in liver is extremely sensitive to EE2 in male fish. However, at lower concentrations (1 and 5 ng/L), the ZP2 in testes shows higher responsiveness to EE2 compared with VTG in rare minnow. The 5' flanking regions of ZP2 and ZP3 were isolated and the comparison of transcription factors in the regions of ZP2 and ZP3 suggested that the disparity for the responsiveness of ZP2 and ZP3 to EE2 could partly be a result of differential cis-elements such as oocyte-specific protein (Osp1) binding sites or/and sex-determining region Y (SRY) binding site.

Use of gene expression data to determine effects on gonad phenotype in Japanese medaka after exposure to trenbolone or estradiol

Various aquatic bioassays using one of several fish species have been developed or are in the process of being developed by organizations like the US Environmental Protection Agency and the Office of Economic Cooperation and Development for testing potential endocrine-disrupting chemicals (EDCs). Often, these involve assessment of the gonad phenotype of individuals as a key endpoint that is inputted into a risk or hazard assessment. Typically, gonad phenotype is determined histologically, which involves specialized and time-consuming techniques. The methods detailed here utilize an entirely different methodology, reverse-transcription quantitative polymerase chain reaction, to determine the relative expression levels of 4 genes after exposure to either 17β-estradiol or 17β-trenbolone and, by extension, the effects of EDCs on the phenotypic status of the gonad. The 4 genes quantified, Sox9b, protamine, Fig1α, and ZPC1, are all involved in gonad development and maintenance in Japanese medaka (Oryzias latipes); these data were then inputted into a permutational multivariate analysis of variance to determine whether significant differences exist between treatment groups. This information in conjunction with the sexual genotype, which can be determined in medaka, can be used to determine adverse effects of exposure to EDCs in a similar fashion to the histologically determined gonad phenotype.

Proteomic characterization and evolutionary analyses of zona pellucida domain-containing proteins in the egg coat of the cephalochordate, Branchiostoma belcheri

Background

Zona pellucida domain-containing proteins (ZP proteins) have been identified as the principle constituents of the egg coat (EC) of diverse metazoan taxa, including jawed vertebrates, urochordates and molluscs that span hundreds of millions of years of evolutionary divergence. Although ZP proteins generally contain the zona pellucida (ZP) structural modules to fulfill sperm recognition and EC polymerization functions during fertilization, the primary sequences of the ZP proteins from the above-mentioned animal classes are drastically different, which makes it difficult to assess the evolutionary relationships of ZP proteins. To understand the origin of vertebrate ZP proteins, we characterized the egg coat components of Branchiostoma belcheri, an invertebrate species that belongs to the chordate subphylum Cephalochordata.

Results

Five ZP proteins (BbZP1-5) were identified by mass spectrometry analyses using the egg coat extracts from both unfertilized and fertilized eggs. In addition to the C-terminal ZP module in each of the BbZPs, the majority contain a low-density lipoprotein receptor domain and a von Willebrand factor type A (vWFA) domain, but none possess an EGF-like domain that is frequently observed in the ZP proteins of urochordates. Fluorescence in situ hybridization and immuno-histochemical analyses of B. belcheri ovaries showed that the five BbZPs are synthesized predominantly in developing eggs and deposited around the extracellular space of the egg, which indicates that they are bona fide egg coat ZP proteins. BbZP1, BbZP3 and BbZP4 are significantly more abundant than BbZP2 and BbZP5 in terms of gene expression levels and the amount of mature proteins present on the egg coats. The major ZP proteins showed high polymorphism because multiple variants are present with different molecular weights. Sequence comparison and phylogenetic analysis between the ZP proteins from cephalochordates, urochordates and vertebrates showed that BbZP1-5 form a monophyletic group and share no significant sequence similarities with the ZP proteins of urochordates and the ZP3 subtype of jawed vertebrates. By contrast, small regions of homology were identifiable between the BbZP and ZP proteins of the non-jawed vertebrate, the sea lamprey Petromyzon marinus. The lamprey ZP proteins were highly similar to the ZP1 and ZP2 subtypes of the jawed vertebrates, which suggests that the ZP proteins of basal chordates most likely shared a recent common ancestor with vertebrate ZP1/2 subtypes and lamprey ZP proteins.

Conclusions

The results document the spectra of zona pellucida domain-containing proteins of the egg coat of basal chordates. Particularly, the study provides solid evidence for an invertebrate origin of vertebrate ZP proteins and indicates that there are diverse domain architectures in ZP proteins of various metazoan groups.

Microarray gene expression profiles from mature gonad tissues of Atlantic bluefin tuna, Thunnus thynnus in the Gulf of Mexico

Background

Bluefin tunas are highly prized pelagic fish species representing a significant economic resource to fisheries throughout the world. Atlantic bluefin tuna (Thunnus thynnus) populations have significantly declined due to overexploitation. As a consequence of their value and population decline, T. thynnus has been the focus of considerable research effort concerning many aspects of their life history. However, in-depth understanding of T. thynnus reproductive biology is still lacking. Knowledge of reproductive physiology is a very important tool for determining effective fisheries and aquaculture management. Transcriptome techniques are proving powerful and provide novel insights into physiological processes. Construction of a microarray from T. thynnus ESTs sourced from reproductive tissues has provided an ideal platform to study the reproductive physiology of bluefin tunas. The aim of this investigation was to compare transcription profiles from the ovaries and testes of mature T. thynnus to establish sex specific variations underlying their reproductive physiology.

Results

Male and females T. thynnus gonad tissues were collected from the wild and histologically staged. Sub-samples of sexually mature tissues were also measured for their mRNA differential expression among the sexes using the custom microarray design BFT 4X44K. A total of 7068 ESTs were assessed for differential expression of which 1273 ESTs were significantly different (p<0.05) with >2 fold change in expression according to sex. Differential expression for 13 of these ESTs was validated with quantitative PCR. These include genes involved in egg envelope formation, hydration, and lipid transport/accumulation more highly expressed in ovaries compared with testis, while genes involved in meiosis, sperm motility and lipid metabolism were more highly expressed in testis compared with ovaries.

Conclusions

This investigation has furthered our knowledge of bluefin tunas reproductive biology by using a contemporary transcriptome approach. Gene expression profiles in T. thynnus sexually mature testes and ovaries were characterized with reference to gametogenesis and potential alternative functions. This report is the first application of microarray technology for bluefin tunas and demonstrates the efficacy by which this technique may be used for further characterization of specific biological aspects for this valuable teleost fish.

Variations in the gene expression of zona pellucida proteins, zpb and zpc, in female European eel (Anguilla anguilla) during induced sexual maturation

Vertebrate eggs are surrounded by an extracellular glycoprotein coat termed zona pellucida (ZP). Integrity of ZP is critical for a correct embryo development. Two zona pellucida protein genes (zpb and zpc) from European eel were characterized, specific qPCR assays developed and their expression in immature males and females carried out. An experimental group of silver-stage eel females was maintained at 18 °C and hormonally induced to sexual maturation by weekly injections of carp pituitary extract during 12 weeks. Changes in zpb and zpc expression during sexual maturation were studied in liver and ovary by qPCR. In liver, no changes were recorded during hormonal treatment, while in ovary expression of both genes decreased during sexual development. These results are a first step in the characterization of ZP in European eel and in the understanding of the mechanism underlying egg envelope formation.

Gene expression profiles in the testis associated with testis-ova in adult Japanese medaka (Oryziaslatipes) exposed to 17α-ethinylestradiol

The occurrence of oocytes in the testis (testis-ova) of several fish species is often associated with exposure of estrogenic chemicals. However, induction mechanisms of the testis-ova remain to be elucidated. To develop marker genes for detecting testis-ova in the testis, adult male medaka were exposed to nominal concentration of 100 ng L(-1) of 17α-ethinylestradiol (EE2) for 3-5 weeks, and 800 ng estradiol benzoate (EB) for 3 weeks (experiment I), and a measured concentration of 20 ng L(-1) EE2 for 1-6 weeks (experiment II). Histological analysis was performed for the testis, and microarray analyses were performed for the testis, liver and brain. Microarray analysis in the estrogen-exposed medaka liver showed vitellogenin and choriogenin as estrogen responsive genes. Testis-ova were induced in the testis after 4 weeks of exposure to 100 ng L(-1) EE2, 3 weeks of exposure to 800 ng EB, and 6 weeks of exposure to 20 ng L(-1) EE2. Microarray analysis of estrogen-exposed testes revealed up-regulation of genes related to zona pellucida (ZP) and the oocytes marker gene, 42Sp50. Using quantitative RT-PCR we confirmed that Zpc5 gene can be used as a marker for the detection of testis-ova in male medaka.

Characterization of rainbow trout gonad, brain and gill deep cDNA repertoires using a Roche 454-Titanium sequencing approach

Rainbow trout, Oncorhynchus mykiss, is an important aquaculture species worldwide and, in addition to being of commercial interest, it is also a research model organism of considerable scientific importance. Because of the lack of a whole genome sequence in that species, transcriptomic analyses of this species have often been hindered. Using next-generation sequencing (NGS) technologies, we sought to fill these informational gaps. Here, using Roche 454-Titanium technology, we provide new tissue-specific cDNA repertoires from several rainbow trout tissues. Non-normalized cDNA libraries were constructed from testis, ovary, brain and gill rainbow trout tissue samples, and these different libraries were sequenced in 10 separate half-runs of 454-Titanium. Overall, we produced a total of 3million quality sequences with an average size of 328bp, representing more than 1Gb of expressed sequence information. These sequences have been combined with all publicly available rainbow trout sequences, resulting in a total of 242,187 clusters of putative transcript groups and 22,373 singletons. To identify the predominantly expressed genes in different tissues of interest, we developed a Digital Differential Display (DDD) approach. This approach allowed us to characterize the genes that are predominantly expressed within each tissue of interest. Of these genes, some were already known to be tissue-specific, thereby validating our approach. Many others, however, were novel candidates, demonstrating the usefulness of our strategy and of such tissue-specific resources. This new sequence information, acquired using NGS 454-Titanium technology, deeply enriched our current knowledge of the expressed genes in rainbow trout through the identification of an increased number of tissue-specific sequences. This identification allowed a precise cDNA tissue repertoire to be characterized in several important rainbow trout tissues. The rainbow trout contig browser can be accessed at the following publicly available web site (http://www.sigenae.org/).

Expression of zona pellucida B proteins in juvenile rare minnow (Gobiocypris rarus) exposed to 17α-ethinylestradiol, 4-nonylphenol and bisphenol A

Zona pellucida (ZP) containing proteins are glycoproteins in teleost chorion and are encoded by several gene subfamilies, mainly including ZPA, ZPB, ZPC and ZPX genes. In teleost species, ZP genes are expressed either in liver under regulation of estrogen or in ovary. In the present study, five ZP gene isoforms were isolated and characterized in Gobiocypris rarus. The putative amino acid sequences of these ZP gene isoforms contain the typical trefoil motif and a ZP domain. These five G. rarus ZP gene isoforms were named as grZPB.1, grZPB.2, grZPB.3, grZPB.4 and grZPB.5. Real-time quantitative reverse transcription PCR (RT-qPCR) analysis indicated that all these ZP mRNA isoforms were exclusively expressed in ovary. G. rarus juveniles at the age of 21 days postfertilization were exposed to 17α-ethinylestradiol (EE2; 0.01, 0.1 and 1 nM), 4-nonylphenol (4-NP; 10, 100 and 1000 nM) or bisphenol A (BPA; 0.1, 1 and 10nM) for 3 days. mRNA expressions of ZPB isoforms following the exposure to xenoestrogen were detected by RT-qPCR. Data were analyzed by the 2(-△△Cq) method. The results indicate that induction by 0.1-1nM EE2 on mRNA expression of the grZPB isoforms is weaker than for vitellogenin. 4-NP exposures at three concentrations had differential effects on the grZPBs. BPA at three concentrations weakly induced mRNA expression of the grZPB isoforms.

Molecular characterization and expression pattern of three zona pellucida 3 genes in the Chinese sturgeon, Acipenser sinensis

Chinese sturgeon (Acipenser sinensis) is a rare and endangered species and also an important resource for the sturgeon aquaculture industry. SMART cDNAs were synthesized from the ovary of A. sinensis, and the full-length cDNAs of three zona pellucida glycoprotein 3 genes (the new gene named AsZP3) were cloned and sequenced. AsZP3.1, AsZP3.2, and AsZP3.3 were 1,388 base pairs (bp), 1,288, and 1,290 bp in length, respectively, and they could be translated into proteins with 440, 394, and 398 amino acids, respectively. High level of amino acids sequence identity was seen between AsZP3.2 and AsZP3.3 (about 82%), but they share low identity with AsZP3.1 (26 and 23%, respectively). The AsZP3.1 has 42-50% amino acids sequence identity values with other fish and lower values with higher vertebrates (38%); AsZP3.2 and AsZP3.3 shared about 30-44% sequence identity with higher vertebrates and other fish. RT-PCR analysis indicated that AsZP3.1 displayed a wide tissue distribution at the mRNA levels including liver, kidney, spleen, heart, and ovary, but AsZP3.2 and AsZP3.3 mRNAs were expressed exclusively in the gonad. All three AsZP3 mRNAs were not detected during embryogenesis and early larval development; furthermore, they were not detected in the gonads of 1- and 2-year-old Chinese sturgeons. All three AsZP3 mRNAs were detected in the testes of 3-year-old males and in the ovaries of 4- and 5-year-old female Chinese sturgeons.

Conservation of the egg envelope digestion mechanism of hatching enzyme in euteleostean fishes

We purified two hatching enzymes, namely high choriolytic enzyme (HCE; EC 3.4.24.67) and low choriolytic enzyme (LCE; EC 3.4.24.66), from the hatching liquid of Fundulus heteroclitus, which were named Fundulus HCE (FHCE) and Fundulus LCE (FLCE). FHCE swelled the inner layer of egg envelope, and FLCE completely digested the FHCE-swollen envelope. In addition, we cloned three Fundulus cDNAs orthologous to cDNAs for the medaka precursors of egg envelope subunit proteins (i.e. choriogenins H, H minor and L) from the female liver. Cleavage sites of FHCE and FLCE on egg envelope subunit proteins were determined by comparing the N-terminal amino acid sequences of digests with the sequences deduced from the cDNAs for egg envelope subunit proteins. FHCE and FLCE cleaved different sites of the subunit proteins. FHCE efficiently cleaved the Pro-X-Y repeat regions into tripeptides to dodecapeptides to swell the envelope, whereas FLCE cleaved the inside of the zona pellucida domain, the core structure of egg envelope subunit protein, to completely digest the FHCE-swollen envelope. A comparison showed that the positions of hatching enzyme cleavage sites on egg envelope subunit proteins were strictly conserved between Fundulus and medaka. Finally, we extended such a comparison to three other euteleosts (i.e. three-spined stickleback, spotted halibut and rainbow trout) and found that the egg envelope digestion mechanism was well conserved among them. During evolution, the egg envelope digestion by HCE and LCE orthologs was established in the lineage of euteleosts, and the mechanism is suggested to be conserved.

Evolution of the teleostean zona pellucida gene inferred from the egg envelope protein genes of the Japanese eel, Anguilla japonica

A fish egg envelope is composed of several glycoproteins, called zona pellucida (ZP) proteins, which are conserved among vertebrate species. Euteleost fishes synthesize ZP proteins in the liver, while otocephalans synthesize them in the growing oocyte. We investigated ZP proteins of the Japanese eel, Anguilla japonica, belonging to Elopomorpha, which diverged earlier than Euteleostei and Otocephala. Five major components of the egg envelope were purified and their partial amino acid sequences were determined by sequencing. cDNA cloning revealed that the eel egg envelope was composed of four ZPC homologues and one ZPB homologue. Four of the five eel ZP (eZP) proteins possessed a transmembrane domain, which is not found in the ZP proteins of Euteleostei and Otocephala that diverged later, but is found in most other vertebrate ZP proteins. This result suggests that fish ZP proteins originally possessed a transmembrane domain and lost it during evolution. Northern blotting and RT-PCR revealed that all of the eZP transcripts were present in the ovary, but not in the liver. Phylogenetic analyses of fish zp genes showed that ezps formed a group with other fish zp genes that are expressed in the ovary, and which are distinct from the group of genes expressed in the liver. Our results support the hypothesis that fish ZP proteins were originally synthesized in the ovary, and then the site of synthesis was switched to the liver during the evolutionary pathway to Euteleostei.

Comparative analysis of the testis and ovary transcriptomes in zebrafish by combining experimental and computational tools

Studies on the zebrafish model have contributed to our understanding of several important developmental processes, especially those that can be easily studied in the embryo. However, our knowledge on late events such as gonad differentiation in the zebrafish is still limited. Here we provide an analysis on the gene sets expressed in the adult zebrafish testis and ovary in an attempt to identify genes with potential role in (zebra)fish gonad development and function. We produced 10,533 expressed sequence tags (ESTs) from zebrafish testis or ovary and downloaded an additional 23,642 gonad-derived sequences from the zebrafish EST database. We clustered these sequences together with over 13,000 kidney-derived zebrafish ESTs to study partial transcriptomes for these three organs. We searched for genes with gonad-specific expression by screening macroarrays containing at least 2600 unique cDNA inserts with testis-, ovary- and kidney-derived cDNA probes. Clones hybridizing to only one of the two gonad probes were selected, and subsequently screened with computational tools to identify 72 genes with potentially testis-specific and 97 genes with potentially ovary-specific expression, respectively. PCR-amplification confirmed gonad-specificity for 21 of the 45 clones tested (all without known function). Our study, which involves over 47,000 EST sequences and specialized cDNA arrays, is the first analysis of adult organ transcriptomes of zebrafish at such a scale. The study of genes expressed in adult zebrafish testis and ovary will provide useful information on regulation of gene expression in teleost gonads and might also contribute to our understanding of the development and differentiation of reproductive organs in vertebrates.

Real-time PCR-based prediction of gonad phenotype in medaka

An important endpoint in aquatic bioassays for potential endocrine disrupting chemicals (EDCs) is the gonadal phenotype of exposed fish, with special interest in intersex and sex-reversed individuals. Traditionally, the assessment of gonad phenotype is done via histology, which involves specialized and time-consuming techniques. The method detailed here increases the efficiency of the analysis by first determining the relative expression of four genes involved in gonad development/maintenance in Japanese medaka (Oryzias latipes), and then by using principal component analysis, assigning a phenotype to each gonad based upon the gene expression data. The gonad phenotype and the sexual genotype, which can be determined in medaka, can then be compared to assess potential adverse effects of exposure to endocrine disrupting chemicals.

Endocrine disruptors modulate expression of hepatic choriogenin genes in the hermaphroditic fish, Kryptolebias marmoratus

Choriogenins (Chgs) are precursors of inner layer of egg envelope that are synthesized in fish liver in response to estrogens. Therefore, study of their expression serves as biomarker of exposure to endocrine disrupting chemicals (EDCs). The self-fertilizing fish, Kryptolebias marmoratus has been established as a model species for testing the action of EDCs. To use this fish as a model for assessing estrogenic activity of EDCs on Chg expression, two K. marmoratus choreogenin genes, Km-ChgH and Km-ChgL were cloned and their expression was analyzed in different tissues and in developmental stages by real-time RT-PCR. Expression levels of liver mRNA were compared between hermaphrodites and secondary males after exposure to EDCs. Km-ChgH and Km-ChgL genes that were predominantly expressed in liver contained zona pellucida (ZP) domains. During embryonic development, low expression of mRNA was observed at stage 1 (2 dpf) that reached highest level at stage 4 (12 dpf) or stage 5 (5 h post hatching). The expression of Km-Chg mRNAs was highly increased in liver exposed to natural estrogen, 17alpha-estradiol (E2) as well as EDCs such as bisphenol A and 4-n-nonylphenol in both the gender types. Another EDC, 4-tert-octylphenol, showed modulatory effect only on Km-ChgH in hermaphrodites. Tamoxifen, an antagonist of the estrogen receptor showed no effect on expression of Chg genes in either of the gender types of K. marmoratus. These findings indicate that Km-Chg genes would be associated with estrogen and measurement of their expression would serve as a surrogate biomarker of exposure to environmental EDCs.

Choriogenin mRNA as a sensitive molecular biomarker for estrogenic chemicals in developing brackish medaka (Oryzias melastigma)

Teleost choriogenins, precursors of the inner layer subunits of the egg envelope, are regarded as sensitive biomarkers for estrogenic pollutants. In this study, two full-length cDNAs, omChgH and omChgL, which encode the choriogenin H and L forms, respectively, were isolated from a brackish medaka, Oryzias melastigma. 17beta-Estradiol (E2; 10 microg/L)-dependent expression of omChgH and omChgL was observed starting at embryonic stage 34 and restricted exclusively to the liver. In hatchlings, E2 induction of omChgH was stronger than that of omChgL. Static exposure of adult fish to E2 (0, 1, 10, 100, and 500 ng/L), 17alpha-ethinylestradiol (EE2; 0, 1, 10, 100 and 500 ng/L), 4-nonylphenol (NP; 0, 1, 10, 100, and 200 microg/L), and bisphenol A (BPA; 0, 1, 10, 100, and 200 microg/L) in artificial seawater for 7 days resulted in dose-dependent induction of both genes in the liver. In the male livers, the sensitivity of omChgH to these estrogenic compounds was higher than that of omChgL; the lowest-observed-effect concentrations (LOECs) of E2, EE2, NP, and BPA on omChgH were 10 ng/L, 10 ng/L, 100 microg/L and 100 microg/L, respectively, and on omChgL were 100 ng/L, 100 ng/L, 100 microg/L, and 200 microg/L, respectively. All these suggest that omChgH can be used as a highly sensitive biomarker for monitoring estrogenic chemicals in the marine environment.

Transcript expression profiles of Takifugu rubripes spermatozoa and eggs by expressed sequence tag analysis

Two cDNA libraries from Takifugu rubripes spermatozoa and eggs were constructed and a total of 620 expressed sequence tag (EST) clones were generated from the two libraries: 300 clones are from the spermatozoa library and 320 clones are from the eggs library. The most abundant cDNA clones in the two libraries were identified. A total of 207 'contigs' (or single) EST clones were found to share significant sequence identity with known sequences in the GenBank database, representing at least 51 different genes. In order to understand the two types of germ cells further, the expression profiles of the identified clones in these cDNA libraries were analyzed. Furthermore, the presence of specific messenger RNAs in the spermatozoa and eggs has been demonstrated with BLAST analysis; the spermatozoa and egg library can supply unique and novel cDNA sequences in the Takifugu rubripes EST project. Another aim of this study is to identify cDNA clones that can be used as molecular markers for the analysis of the spermatogenesis and oogenesis in Takifugu rubripes. Six potential clones (S1-3 from spermatozoa and E1-3 from eggs) were selected to analyze their expression patterns by reverse transcription (RT)-PCR analyses. Half of these showed a specific expression in the expected tissue. Two of the clones were found by RT-PCR and in situ hybridization to be expressed specifically in the testis or ovary, and they maybe suitable molecular markers for the analysis of spermatogenesis and oogenesis.

Comparative genomics approach to the expression of figalpha, one of the earliest marker genes of oocyte differentiation in medaka (Oryzias latipes)

We analyzed molecular cascades of sex differentiation in medaka gonads by examining the transcriptional regulation of the oocyte-expressed gene, figalpha. We first confirmed that figalpha is one of the earliest marker genes of oocyte differentiation by quantitative RT-PCR and in situ hybridization. Expression of putative figalpha target genes, zpc4 and zpb, followed that of figalpha. A meiosis-specific gene, scp3, showed expression temporally and spatially similar to figalpha. To characterize the cis-regulatory sequences of figalpha, we compared genomic organizations of vertebrate figalpha genes. Both number and sequence homology of the C-terminal exons showed divergence, suggesting their less important roles. In the frog, Xenopus tropicalis, and in many teleosts, figalpha is located between hexokinase 2 and beta-adducin. We compared this genomic region for potential cis-regulatory elements and found no DNA stretches with high homology. In spite of this lack of sequence similarities, fluorescent protein transgenes surrounded with figalpha flanking sequences from the compact genomes of fugu or Tetraodon faithfully reproduced the endogenous expression of figalpha in the medaka oocytes, indicating conserved regulatory mechanisms.

ZP genes in avian species illustrate the dynamic evolution of the vertebrate egg envelope

The vertebrate egg envelope is composed of a set of related proteins, encoded by the ZP genes. The apparent simplicity of the egg envelope is in contrast to the number of ZP genes identified by conventional cloning and data mining of genome sequences from a number of vertebrates. The vertebrate ZP genes fall into five classes, ZP1, ZP2, ZP3, ZPD and ZPAX. Analysis of chicken genome and EST sequence data has revealed the presence of seven distinct ZP genes, falling into these classes that are expressed in the female reproductive system. Comparison with the repertoire of ZP genes in other vertebrates suggests a major source of diversity in the composition of the egg envelope is a continual process of amplification, diversification and attrition of ZP gene sequences.

Zona Pellucida Domain of ZPB1 controls specific binding of ZPB1 and ZPC in Japanese quail (Coturnix japonica)

The extracellular matrix surrounding avian oocytes, referred to as the perivitelline membrane (PL), exhibits a three-dimensional network of fibrils between granulosa cells and the oocyte. We previously reported that one of its components, ZPC, is synthesized in granulosa cells that are specifically incorporated into the PL; this incorporation might be mediated by a specific interaction with ZPB1, another PL constituent, which is synthesized in the liver. In order to extend our previous findings, we established an expression system for quail ZPB1 using a mammalian cell line, and several ZPB1 mutants lacking the zona pellucida (ZP) domain or the glutamine-rich repeat region were produced. Western blot analysis of the immunoprecipitated materials with anti-ZPC antiserum indicated that ZPB1 was coimmunoprecipitated with the antiserum in the presence of ZPC. Ligand blotting also revealed the specific binding of ZPC and ZPB1 and indicated that the binding of these two components might be mediated via an ionic interaction. An analysis using recombinant ZPB1 demonstrated that the ZPB1 lacking the ZP domain did not bind to ZPC, whereas the mutant missing the glutamine-rich repeat region retained its capacity for binding. Furthermore, although the ZPB1 lacking the N-terminal half of the ZP domain was able to bind to ZPC, the deletion of the C-terminal half completely abolished ZPB1 binding to ZPC. These results suggested that the C-terminal half of the ZP domain of ZPB1 contains a binding site for ZPC, and that it appears to be involved in insoluble PL fibril formation in the quail ovary.