Transcriptomic signatures for ovulation in vertebrates.

The central roles of luteinizing hormone (LH), progestin and their receptors for initiating ovulation have been well established. However, signaling pathways and downstream targets such as proteases that are essential for the rupture of follicular cells are still unclear. Recently, we found anovulation in nuclear progestin receptor (Pgr) knockout (Pgr-KO) zebrafish, which offers a new model for examining genes and pathways that are important for ovulation and fertility. In this study, we examined expression of all transcripts using RNA-Seq in preovulatory follicular cells collected following the final oocyte maturation, but prior to ovulation, from wild-type (WT) or Pgr-KO fish. Differential expression analysis revealed 3567 genes significantly differentially expressed between WT and Pgr-KO fish (fold change⩾2, p<0.05). Among those, 1543 gene transcripts were significantly more expressed, while 2024 genes were significantly less expressed, in WT than those in Pgr-KO. We then retrieved and compared transcriptional data from online databases and further identified 661 conserved genes in fish, mice, and humans that showed similar levels of high (283 genes) or low (387) expression in animals that were ovulating compared to those with no ovulation. For the first time, ovulatory genes and their involved biological processes and pathways were also visualized using Enrichment Map and Cytoscape. Intriguingly, enrichment analysis indicated that the genes with higher expression were involved in multiple ovulatory pathways and processes such as inflammatory response, angiogenesis, cytokine production, cell migration, chemotaxis, MAPK, focal adhesion, and cytoskeleton reorganization. In contrast, the genes with lower expression were mainly involved in DNA replication, DNA repair, DNA methylation, RNA processing, telomere maintenance, spindle assembling, nuclear acid transport, catabolic processes, and nuclear and cell division. Our results indicate that a large set of genes (>3000) is differentially regulated in the follicular cells in zebrafish prior to ovulation, terminating programs such as growth and proliferation, and beginning processes including the inflammatory response and apoptosis. Further studies are required to establish relationships among these genes and an ovulatory circuit in the zebrafish model.

Identification, expression and antibacterial activities of an antimicrobial peptide NK-lysin from a marine fish Larimichthys crocea.

As fundamental immunologic mechanism, the innate immunity system is more important than the specific immunity system in teleost fishes during pathogens infection. Antimicrobial peptides are integral parts of the innate immune system, and play significant roles against pathogens infection. NK-lysin, the compounds of the natural killer cells and cytotoxic T cells, are potent and effective antimicrobial peptides widely distributed in animals. In this study, we reported the sequence characteristics, expression profiles and antibacterial activities of a NK-lysin gene (Lc-NK-lysin) from a commercially important marine fish, the large yellow croaker (Larimichthys crocea). The open reading frame of Lc-NK-lysin cDNA sequence was 447 bp in length, coding 148 amino acids. The genomic DNA of Lc-NK-lysin has the common features of NK-lysin family, consisting of five exons and four introns, and in its deduced mature peptide, there are six well-conserved cysteine residues and a Saposin B domain. Lc-NK-lysin was expressed in all tested tissues (skin, muscle, gill, brain, head kidney, heart, liver, spleen, stomach and intestine) with different expression patterns. In pathogens infection the expression profiles of Lc-NK-lysin varied significantly in gill, head kidney, spleen and liver, indicating its role in immune response. Two peptides (Lc-NK-lysin-1 and Lc-NK-lysin-2) divided from the core region of the Lc-NK-lysin mature polypeptide were chemically synthesized and their antibacterial activities were examined; the potential function on the inhibition of bacteria propagation was revealed. Our results suggested that Lc-NK-lysin is a typical member of the NK-lysin family and as an immune-related gene it involves in the immune response when pathogens invasion.

Dynamic methylation pattern of cyp19a1a core promoter during zebrafish ovarian folliculogenesis.

In vertebrates, the aromatase coded by the cyp19a1a gene can catalyze the conversion from androgens to estrogens. Thus, the regulatory mechanisms of cyp19a1a gene expression are a critical research field in reproductive endocrinology. In this study, we use zebrafish as a model to study the dynamic methylation levels of the cyp19a1a gene core promoter during zebrafish ovarian folliculogenesis. The results show that there is an apparent fluctuation of the methylation levels of zebrafish cyp19a1a core promoter. Moreover, the methylation levels are inversely correlated with the expression levels of cyp19a1a transcripts when the ovarian follicles develop from PV into the MV stage. Also, the CpG dinucleotides which are close to the transcriptional starting site may have provided a significant blocking effect on inhibiting the transcriptional function of RNA polymerase II. Taken together, the results from the present study strongly suggest that DNA methylation was one of mechanisms that are involved in the regulation of cyp19a1a gene expression during folliculogenesis. This methylation mechanism modifying transcriptional process accompanied with zebrafish ovarian folliculogenesis might also shed new light on the regulation of cyp19a1a expression during the ovarian developmental stage in other vertebrates.

Molecular characterization and expression analysis of AMPK α subunit isoform genes from Scophthalmus maximus responding to salinity stress.

AMP-activated protein kinase (AMPK) is a highly conserved and multi-functional protein kinase that plays important roles in both intracellular energy balance and cellular stress response. In the present study, molecular characterization, tissue distribution and gene expression levels of the AMPK α1 and α2 genes from turbot (Scophthalmus maximus) under salinity stress are described. The complete coding regions of the AMPK α1 and α2 genes were isolated from turbot through degenerate primers in combination with RACE using muscle cDNA. The complete coding regions of AMPK α1 (1722 bp) and α2 (1674 bp) encoded 573 and 557 amino acids peptides, respectively. Multiple alignments, structural analysis and phylogenetic tree construction indicated that S. maximus AMPK α1 and α2 shared a high amino acid identity with other species, especially fish. AMPK α1 and α2 genes could be detected in all tested tissues, indicating that they are constitutively expressed. Salinity challenges significantly altered the gene expression levels of AMPK α1 and α2 mRNA in a salinity- and time-dependent manners in S. maximus gill tissues, suggesting that AMPK α1 and α2 played important roles in mediating the salinity stress in S. maximus. The expression levels of AMPK α1 and α2 mRNA were a positive correlation with gill Na+, K+-ATPase activities. These findings will aid our understanding of the molecular mechanism of juvenile turbot in response to environmental salinity changes.

Characterization of RAG1 and IgM (mu chain) marking development of the immune system in red-spotted grouper (Epinephelus akaara).

In vertebrates, lymphoid-specific recombinase protein encoded by recombination-activating genes (RAG1/2) plays a key role in V(D)J recombination of the T-cell receptor and B-cell receptor. In this study, both RAG1 and the immunoglobulin M (IgM) mu chain were cloned to characterize their potential role in the immune defense at developmental stages of red-spotted grouper, Epinephelus akaara. The open reading frame (ORF) of E. akaara RAG1 included 2778 nucleotide residues encoding a putative protein of 925 amino acids, while the ORF of the IgM mu chain had 1734 nucleotide residues encoding 578 amino acids including variable (VH) and constant (CH1-CH2-CH3-CH4) regions. E. akaara RAG1 was composed of a zinc-binding dimerization domain (ZDD) with a RING finger and zinc finger A (ZFA) in the non-core region and a nonamer-binding region (NBR), with a zinc finger B (ZFB), the central and C-terminal domains in the core region. Tridimensional models of the ZDD and NBR of E. akaara RAG1 were constructed for the first time in fishes, while a 3D model of the E. akaara IgM mu chain was also clarified. The RAG1 mRNA was only detected in the thymus and kidney of 4-month and 1.5-year old groupers using qPCR, and the RAG1 protein was confirmed using western blotting and immunohistochemistry. The IgM mu mRNA was examined in most tissues except the gonad. RAG1 and IgM mu gene expression were observed at 15 dph (days post-hatching) and 23 dph respectively, and increased to a higher level at 37 dph. In addition, this was the first time that the morphology of the E. akaara thymus was characterized. The oval-shaped thymus of 4-month old fish was clearly seen and there were amounts of T lymphocytes present. The results suggested that the immune action of E. akaara was likely to start to develop around 15 dph to 29 dph. The transcript level of the RAG1 gene and the number of lymphocytes in the thymus between 4-month and 1.5-year old groupers indicated that age-related thymic atrophy also occurs in fishes. The similar functional structures of RAG1 and IgM protein between fish and mammals indicated that teleost species share a similar mechanism of V(D)J recombination with higher vertebrates.

Impaired oocyte maturation and ovulation in membrane progestin receptor (mPR) knockouts in zebrafish.

Accumulating evidence suggest that membrane progestin receptor α (mPRα) is the membrane receptor mediating nongenomic progestin signaling that induces oocyte maturation in teleost. However, the involvement of other members of mPR family in oocyte maturation is still unclear. In this study, we found impaired oocyte maturation in zebrafish lacking mPRα1, mPRα2, mPRß, or mPRγ2. In contrast, no difference was observed in oocyte maturation in the single knockout of mPRγ1, mPRδ, or mPRε. To study possible redundant functions of different mPRs in oocyte maturation, we generated a zebrafish line lacking all seven kinds of mPRs (mprs-/-). We found oocyte maturation was further impaired in mprs-/-. In addition, oocyte ovulation delay was observed in mprs-/- females, which was associated with low levels of nuclear progestin receptor (Pgr), a key regulator for ovulation. We also found reduced fertility in mprs-/- female zebrafish. Furthermore, eggs spawned by mprs-/- females were of poor quality.

Effect of a Sodium Alginate Coating Infused with Tea Polyphenols on the Quality of Fresh Japanese Sea Bass (Lateolabrax japonicas) Fillets.

Sodium alginate (SA) and tea polyphenols (TP) are natural preservatives commonly used in the food industry, including the production of fish products. The effect of SA coating infused with TP on the quality of fresh Japanese sea bass (Lateolabrax japonicas) fillets was evaluated over a 20-day period at 4 °C. SA (1.5%, w/v) or TP (0.5%, w/v) treatment alone, and the SA coating infused with TP (SA-TP) all reduced microbial counts, with the SA-TP providing the greatest effect. Fish fillet samples treated with SA-TP had significantly lower levels of total volatile basic nitrogen, lipid oxidation, and protein decomposition during the storage period, relative to the remaining treatments. The samples treated with SA-TP had the highest sensory quality rating as well. Collectively, sodium alginate coating infused with tea polyphenols may represent a promising treatment for preservation of Japanese sea bass fillets during cold storage. PRACTICAL APPLICATION: The sodium alginate-tea polyphenols composite coating has strong potential to be used as a new biopreservative for maintaining fish fillet quality.

Complete mitochondrial genome and phylogenetic analysis of the barred mudskipper Periophthalmus argentilineatus (Perciformes, Gobiidae).

The complete mitochondrial genome sequence of the barred mudskipper Periophthalmus argentilineatus was first determined in this study. The circle genome was 16 509 bp long and consisted of 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and 1 control region. The nucleotide composition of the heavy strand of P. argentilineatus is 28.28% for A, 27.83% for C, 16.01% for T, and 27.88% for G, with a slight G + C bias of 55.71%. Only the NADH dehydrogenase subunit 6 (ND6) and eight tRNA genes were encoded on the light strand, other mitochondrial genes were all encoded on the heavy strand. The mitochondrial gene arrangement of P. argentilineatus is similar to those of most other gobies. The phylogenetic analysis using the neighbor-joining method showed that the kinship between Periophthalmus and Boleophthalmus is closer than those between Periophthalmus and other selected genera. Periophthalmus argentilineatus was clustered into one branch with other four species from the same genus Periophthalmus.

Complete mitochondrial genome and phylogenetic analysis of the serpent mudskipper Parapocryptes serperaster (Perciformes, Gobiidae): repetitive sequences in the control region.

The complete mitochondrial genome sequence of the serpent mudskipper Parapocryptes serperaster was first determined in this study. The circle genome was 17 243 bp long and consisted of 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes, and 1 control region. The nucleotide composition of the heavy strand of P. serperaster is 28.65% for A, 30.07% for C, 25.31% for T, and 15.97% for G, with a slight A + T bias of 53.96%. Only the NADH dehydrogenase subunit 6 (ND6) and eight tRNA genes were encoded on the light strand, all other mitochondrial genes were encoded on the heavy strand. A 178-bp tandem repeat, which started from the 827-bp site of the control region, was identified. Like most other gobies, P. serperaster also has the typical vertebrate mitochondrial gene arrangement. The phylogenetic analysis showed that the kinship between Parapocryptes and Boleophthalmus is closer than those between Parapocryptes and other selected genera.

Complete mitochondrial genome and phylogenic analysis of the mudskipper Scartelaos gigas (Perciformes, Gobiidae).

In this paper, the complete mitochondrial genome sequence of Scartelaos gigas was firstly determined. The circular genome (16 717 bp) comprises 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 1 control region. The overall base composition of S. gigas is 28.9% for C, 28.3% for A, 26.4% for T, 16.4% for G, with a slight A + T bias of 54.7%. In the control region, the termination-associated sequence and conserved sequence block domains were found, but the tandem repeat structure was not found. It has the typical vertebrate mitochondrial gene arrangement. The phylogenic analysis using the Neighbor-Joining method showed that the fishes belonging to Gobiidae, Odontoburidae, and Eleotridae formed three branches grouped with other fishes into one clade which separated from the mammals. We hope that the results from the present study will provide useful molecular information for the further studies on genetic structure and demographic history of S. gigas.

Progestin increases the expression of gonadotropins in pituitaries of male zebrafish.

Our previous study showed that the in vivo positive effects of 17α,20ß-dihydroxy-4-pregnen-3-one (DHP), the major progestin in zebrafish, on early spermatogenesis was much stronger than the ex vivo ones, which may suggest an effect of DHP on the expression of gonadotropins. In our present study, we first observed that fshb and lhb mRNA levels in the pituitary of male adult zebrafish were greatly inhibited by 3 weeks exposure to 10nM estradiol (E2). However, an additional 24h 100nM DHP exposure not only reversed the E2-induced inhibition, but also significantly increased the expression of fshb and lhb mRNA. These stimulatory effects were also observed in male adult fish without E2 pretreatment, and a time course experiment showed that it took 24h for fshb and 12h for lhb to respond significantly. Because these stimulatory activities were partially antagonized by a nuclear progesterone receptor (Pgr) antagonist mifepristone, we generated a Pgr-knockout (pgr(-/-)) model using the TALEN technique. With and without DHP in vivo treatment, fshb and lhb mRNA levels of pgr(-/-) were significantly lower than those of pgr(+/+) Furthermore, ex vivo treatment of pituitary fragments of pgr(-/-) with DHP stimulated lhb, but not fshb mRNA expression. Results from double-colored fluorescent in situ hybridization showed that pgr mRNA was expressed only in fshb-expressing cells. Taken together, our results indicated that DHP participated in the regulation of neuroendocrine control of reproduction in male zebrafish, and exerted a Pgr-mediated direct stimulatory effect on fshb mRNA at pituitary level.

Molecular cloning, characterization and expression of Lc-Sox11a in large yellow croaker Larimichthys crocea.

Sox genes play important roles in various developmental processes such as sex determination, embryogenesis, oogenesis, neurogenesis, and larval development. In order to clarify the roles of Sox genes in the developmental process of large yellow croaker, the full-length cDNA of the Sox11a gene (Lc-Sox11a) was cloned for the first time. Bioinformatics analysis indicated that Lc-Sox11a contains a protein of 366 amino acids with a Ser-rich region, a C-terminal conserved region, and a high mobility group box. The expression of Lc-Sox11a in different tissues of both sexes and in different developmental embryonic stages revealed that Lc-Sox11a were expressed with tissue and gender specificity, of which the expression level in female was ovary>brain>eye>gill; in male was brain>testis>gill. The gender differences occurred in the brain and eye with the male brain>female brain, female eye>male eye. Moreover, the expression of Lc-Sox11a in the gonad and brain at different growth stages was detected. Significant up-regulated expression of Lc-Sox11a was found in the ovary and the male brain at 1000dph (days post hatching) compared with 270dph and 635dph. However, significant down-regulated expression of Lc-Sox11a occurred in the testis with growth. Besides, the expression of Lc-Sox11a in the female brain showed a trend of first rising then falling, with the highest peak in 635dph. The results of in situ hybridization displayed that Lc-Sox11a was widely distributed only in cytoplasm of oocytes at each stage in oogenesis. In early stage of oocytes, Lc-Sox11a was expressed weakly and evenly. As the appearance of vacuoles and synthesis of yolks, positive signals of Lc-Sox11a distributed intensively in the residual cytoplasm. In spermatogenesis, Lc-Sox11a was distributed in cytoplasm of all male germ cells except spermatozoon with spermatogonium>spermatocyte>spermatid. During embryogenesis, Lc-Sox11a was expressed in most embryonic stages, the highest expression occurred in the formation-of-eye-lens stage, closely followed by the closure-of-blastopore stage, then the beginning-of-heart-pulsation stage. The results of whole mount in situ hybridization showed that the expression of Lc-Sox11a began to increase beginning with the multiple-cell stage, with the major distribution of Lc-Sox11a in the brain and eye areas in the pre-hatching stage.

Nuclear progestin receptor (pgr) knockouts in zebrafish demonstrate role for pgr in ovulation but not in rapid non-genomic steroid mediated meiosis resumption.

Progestins, progesterone derivatives, are the most critical signaling steroid for initiating final oocyte maturation (FOM) and ovulation, in order to advance fully-grown immature oocytes to become fertilizable eggs in basal vertebrates. It is well-established that progestin induces FOM at least partly through a membrane receptor and a non-genomic steroid signaling process, which precedes progestin triggered ovulation that is mediated through a nuclear progestin receptor (Pgr) and genomic signaling pathway. To determine whether Pgr plays a role in a non-genomic signaling mechanism during FOM, we knocked out Pgr in zebrafish using transcription activator-like effector nucleases (TALENs) and studied the oocyte maturation phenotypes of Pgr knockouts (Pgr-KOs). Three TALENs-induced mutant lines with different frame shift mutations were generated. Homozygous Pgr-KO female fish were all infertile while no fertility effects were evident in homozygous Pgr-KO males. Oocytes developed and underwent FOM normally in vivo in homozygous Pgr-KO female compared to the wild-type controls, but these mature oocytes were trapped within the follicular cells and failed to ovulate from the ovaries. These oocytes also underwent normal germinal vesicle breakdown (GVBD) and FOM in vitro, but failed to ovulate even after treatment with human chronic gonadotropin (HCG) or progestin (17α,20ß-dihydroxyprogesterone or DHP), which typically induce FOM and ovulation in wild-type oocytes. The results indicate that anovulation and infertility in homozygous Pgr-KO female fish was, at least in part, due to a lack of functional Pgr-mediated genomic progestin signaling in the follicular cells adjacent to the oocytes. Our study of Pgr-KO supports previous results that demonstrate a role for Pgr in steroid-dependent genomic signaling pathways leading to ovulation, and the first convincing evidence that Pgr is not essential for initiating non-genomic progestin signaling and triggering of meiosis resumption.

Mudskipper genomes provide insights into the terrestrial adaptation of amphibious fishes.

Mudskippers are amphibious fishes that have developed morphological and physiological adaptations to match their unique lifestyles. Here we perform whole-genome sequencing of four representative mudskippers to elucidate the molecular mechanisms underlying these adaptations. We discover an expansion of innate immune system genes in the mudskippers that may provide defence against terrestrial pathogens. Several genes of the ammonia excretion pathway in the gills have experienced positive selection, suggesting their important roles in mudskippers' tolerance to environmental ammonia. Some vision-related genes are differentially lost or mutated, illustrating genomic changes associated with aerial vision. Transcriptomic analyses of mudskippers exposed to air highlight regulatory pathways that are up- or down-regulated in response to hypoxia. The present study provides a valuable resource for understanding the molecular mechanisms underlying water-to-land transition of vertebrates.

Variation with semilunar periodicity of plasma steroid hormone production in the mudskipper Boleophthalmus pectinirostris.

Variation in the production of the plasma steroid hormones E(2), 17alpha-OHP and T in females and T and 11-KT in males, was investigated in the mudskipper Boleophthalmus pectinirostris during the spawning season. Females with oocytes at the vitellogenic stage (GSI 5.97-6.86%) and mature males with GSI of 0.255-0.288% were collected at intervals of 3-4 days within the two complete semilunar cycles from May 31 to June 30, 2006. The results showed that variations in the levels of plasma steroid hormones were synchronized obviously with semilunar periodicity in both females and males. Each steroid hormone level exhibited two cycles, each cycle with a peak. In females, the first peaks in plasma E(2), 17alpha-OHP and T levels were observed 3 days after the first lunar quarter, and the second ones, 4 days after the last lunar quarter. In males, the first peaks of plasma T and 11-KT levels occurred 3 days after the first lunar quarter, and the second ones, at the last lunar quarter. The fact that, in the present study, changes in the levels of plasma steroid hormones were synchronized with semilunar periodicity, although the fish were at the same stages of gonadal development, suggests that variation of plasma steroid hormones is basically regulated by biological rhythms (Zeitgebers), and that tidal movement (with its semilunar periodicity) is the major environmental factor stimulating steroid hormone production in B. pectinirostris.

[Effects of BaP exposure on ultrastructures of hepatic cells of Boleophthalmus pectinirostris].

The changes of ultrastructures of hepatic cells of Boleophthalmus pectinirostris were investigated after the fish were exposed under benzo(a) pyrene in different concentrations under experimental condition. The results showed that the organelles in hepatic cells of B. pectinirostris were damaged to different extents after the fish was exposed under lower concentration of BaP (0.5 mg.L-1) for up to 7 d, in which, mitochondria and endoplasmic reticulum were the chief organelles affected by BaP exposure. While the fish was exposed under higher concentration of BaP (5 mg.L-1) for 2 h, almost all of the organelles including mitochondria and endoplasmic reticulum in hepatic cells of B. pectinirostr were affected by BaP exposure. The structures of liver cells were seriously damaged. It was demonstrated that BaP could produce multiorganalle lesions in hepatic cells of B. pectinirostris, and the severity extent of such lesions was dependent on the concentration level of BaP.