Protein expression patterns and metal metabolites in a protogynous hermaphrodite fish, the ricefield eel (Monopterus albus).

BACKGROUND: The ricefield eel Monopterus albus undergoes a natural sex change from female to male during its life cycle, and previous studies have shown the potential mechanisms of this transition at the transcriptional and protein levels. However, the changes in protein levels have not been fully explored, especially in the intersexual stage. RESULTS: In the present study, the protein expression patterns in the gonadal tissues from five different periods, the ovary (OV), early intersexual stage gonad (IE), middle intersexual stage gonad (IM), late intersexual stage gonad (IL), and testis (TE), were determined by untargeted proteomics sequencing. A total of 5125 proteins and 394 differentially expressed proteins (DEPs) were detected in the gonadal tissues. Of the 394 DEPs, there were 136 between the OV and IE groups, 20 between the IM and IE groups, 179 between the IL and IM groups, and 59 between the TE and IL groups. Three candidate proteins, insulin-like growth factor 2 mRNA-binding protein 3 isoform X1 (Igf2bp3), triosephosphate isomerase (Tpi), and Cu-Zn superoxide dismutase isoform X1 [(Cu-Zn) Sod1], were validated by western blotting to verify the reliability of the data. Furthermore, metal metabolite-related proteins were enriched in the IL vs. IM groups and TE vs. IL groups, which had close relationships with sex change, including Cu2+-, Ca2+-, Zn2+- and Fe2+/Fe3+-related proteins. Analysis of the combined transcriptome data revealed consistent protein/mRNA expression trends for two metal metabolite-related proteins/genes [LOC109953912 and calcium Binding Protein 39 Like (cab39l)]. Notably, we detected significantly higher levels of Cu2+ during the sex change process, suggesting that Cu2+ is a male-related metal metabolite that may have an important function in male reproductive development. CONCLUSIONS: In summary, we analyzed the protein profiles of ricefield eel gonadal tissues in five sexual stages (OV, IE, IM, IL, and TE) and verified the plausibility of the data. After preforming the functional enrichment of metal metabolite-related DEPs, we detected the contents of the metal metabolites Zn2+, Cu2+, Ca2+, and Fe2+/Fe3+ at these five stages and screened for (Cu-Zn) Sod1 and Mmp-9 as possible key proteins in the sex reversal process.

Expression Patterns and Gonadotropin Regulation of the TGF-ß II Receptor (Bmpr2) during Ovarian Development in the Ricefield Eel Monopterus albus.

Bmpr2 plays a central role in the regulation of reproductive development in mammals, but its role during ovarian development in fish is still unclear. To ascertain the function of bmpr2 in ovarian development in the ricefield eel, we isolated and characterized the bmpr2 cDNA sequence; the localization of Bmpr2 protein was determined by immunohistochemical staining; and the expression patterns of bmpr2 in ovarian tissue incubated with FSH and hCG in vitro were analyzed. The full-length bmpr2 cDNA was 3311 bp, with 1061 amino acids encoded. Compared to other tissues, bmpr2 was abundantly expressed in the ovary and highly expressed in the early yolk accumulation (EV) stages of the ovary. In addition, a positive signal for Bmpr2 was detected in the cytoplasm of oocytes in primary growth (PG) and EV stages. In vitro, the expression level of gdf9, the ligand of bmpr2, in the 10 ng/mL FSH treatment group was significantly higher after incubation for 4 h than after incubation for different durations. However, bmpr2 expression in the 10 ng/mL FSH treatment group at 2 h, 4 h and 10 h was significantly lower. Importantly, the expression level of bmpr2 and gdf9 in the 100 IU/mL hCG group had similar changes that were significantly decreased at 4 h and 10 h. In summary, Bmpr2 might play a pivotal role in ovarian growth in the ricefield eel, and these results provide a better understanding of the function of bmpr2 in ovarian development and the basic data for further exploration of the regulatory mechanism of gdf9 in oocyte development.

Circular RNA expression profiles and CircSnd1-miR-135b/c-foxl2 axis analysis in gonadal differentiation of protogynous hermaphroditic ricefield eel Monopterus albus.

BACKGROUND: The expression and biological functions of circular RNAs (circRNAs) in reproductive organs have been extensively reported. However, it is still unclear whether circRNAs are involved in sex change. To this end, RNA sequencing (RNA-seq) was performed in gonads at 5 sexual stages (ovary, early intersexual stage gonad, middle intersexual stage gonad, late intersexual stage gonad, and testis) of ricefield eel, and the expression profiles and potential functions of circRNAs were studied. RESULTS: Seven hundred twenty-one circRNAs were identified, and the expression levels of 10 circRNAs were verified by quantitative real-time PCR (qRT-PCR) and found to be in accordance with the RNA-seq data, suggesting that the RNA-seq data were reliable. Then, the sequence length, category, sequence composition and the relationship between the parent genes of the circRNAs were explored. A total of 147 circRNAs were differentially expressed in the sex change process, and GO and KEGG analyses revealed that some differentially expressed (such as novel_circ_0000659, novel_circ_0004005 and novel_circ_0005865) circRNAs were closely involved in sex change. Furthermore, expression pattern analysis demonstrated that both circSnd1 and foxl2 were downregulated in the process of sex change, which was contrary to mal-miR-135b. Finally, dual-luciferase reporter assay and RNA immunoprecipitation showed that circSnd1 and foxl2 can combine with mal-miR-135b and mal-miR-135c. These data revealed that circSnd1 regulates foxl2 expression in the sex change of ricefield eel by acting as a sponge of mal-miR-135b/c. CONCLUSION: Our results are the first to demonstrate that circRNAs have potential effects on sex change in ricefield eel; and circSnd1 could regulate foxl2 expression in the sex change of ricefield eel by acting as a sponge of mal-miR-135b/c. These data will be useful for enhancing our understanding of sequential hermaphroditism and sex change in ricefield eel or other teleosts.

Identification, characterization and functional analysis of gonadal long noncoding RNAs in a protogynous hermaphroditic teleost fish, the ricefield eel (Monopterus albus).

BACKGROUND: An increasing number of long noncoding RNAs (lncRNAs) have been found to play important roles in sex differentiation and gonad development by regulating gene expression at the epigenetic, transcriptional and posttranscriptional levels. The ricefield eel, Monopterus albus, is a protogynous hermaphroditic fish that undergoes a sequential sex change from female to male. However, the roles of lncRNA in the sex change is unclear. RESULTS: Herein, we performed RNA sequencing to analyse lncRNA expression patterns in five different stages of M. albus development to investigate the roles of lncRNAs in the sex change process. A total of 12,746 lncRNAs (1503 known lncRNAs and 11,243 new lncRNAs) and 2901 differentially expressed lncRNAs (DE-lncRNAs) were identified in the gonads. The target genes of the DE-lncRNAs included foxo1, foxm1, smad3, foxr1, camk4, ar and tgfb3, which were mainly enriched in signalling pathways related to gonadal development, such as the insulin signalling pathway, MAPK signalling pathway, and calcium signalling pathway. We selected 5 highly expressed DE-lncRNAs (LOC109952131, LOC109953466, LOC109954337, LOC109954360 and LOC109958454) for full length amplification and expression pattern verification. They were all expressed at higher levels in ovaries and intersex gonads than in testes, and exhibited specific time-dependent expression in ovarian tissue incubated with follicle-stimulating hormone (FSH) and human chorionic gonadotropin (hCG). The results of quantitative real-time PCR (qRT-PCR) analysis and a dual-luciferase assay showed that znf207, as the gene targeted by LOC109958454, was expressed in multiple tissues and gonadal developmental stages of M. albus, and its expression was also inhibited by the hormones FSH and hCG. CONCLUSIONS: These results provide new insights into the role of lncRNAs in gonad development, especially regarding natural sex changes in fish, which will be useful for enhancing our understanding of sequential hermaphroditism and sex changes in the ricefield eel (M. albus) and other teleosts.

Expression and regulation of Smad2 by gonadotropins in the protogynous hermaphroditic ricefield eel (Monopterus albus).

Smad2, a receptor-activated Smad, plays a critical role in regulating gametogenesis. In this study, a smad2 homologue was identified and sequenced from ricefield eel ovary cDNA, and its mRNA and protein expression levels were analysed during oocyte development. The cDNA sequence of ricefield eel smad2 consisted of 1863 bp encoding a 467-amino acid protein that had high sequence homology with Smad proteins in other teleosts, especially in Poeciliopsis prolifica. The results of real-time quantitative PCR (RT-qPCR) analysis revealed that smad2 is expressed in the ovary during gonad development, increased continuously until the early vitellogenic stage in the ovaries, and then decreased with ovary maturation. Smad2 protein immunoreactivity was localized in the cytoplasm of follicular cells, oogonia, and primary growth stage oocytes. In vitro experiments revealed that follicle-stimulating hormone (FSH) and human chorionic gonadotropin (hCG) promoted smad2 expression in ovary tissue in a time- and dose-dependent manner, respectively. In summary, Smad2 plays a potentially vital role in ricefield eel ovary development.

Genetic analysis of wild Ancherythroculter nigrocauda in tributaries and the main stream of the upper Yangtze River basin of China.

Biodiversity is threatened by several factors that are often associated with overfishing, water pollution and hydroelectric dams, among other environmental impacts. The present study aimed to evaluate the genetic aspects of wild groups of Ancherythroculter nigrocauda using the mitochondrial cytochrome c oxidase subunit I (coI) and cytochrome b (cytb) genes and the d-loop region. We collected 89 representative individuals from three geographically distinct ranges of the Upper Yangtze River, including the Longxi River (LOR), Laixi River (LAR), and Hejiang range of the Yangtze River (HJ). The genetic analysis results showed that the three populations of A. nigrocauda had high levels of haplotype diversity (0.3434-0.951) and low levels of nucleotide diversity (0.00074-0.00412) based on the single gene sequences and the combination of gene sequences. Haplotype genealogy showed that only one haplotype (Hap-2) was shared by these three geographic groups, and 2-3 were shared by two groups; the other haplotypes were group-specific. The genetic distance within and between the populations was low; however, most of the molecular variance came from within the populations. Furthermore, high gene flow (>1.0) was found in HJ vs LOR and HJ vs LAR based on the d-loop region sequence and combination. These results suggested that there was a decrease in the degree of A. nigrocauda genetic diversity in the upper Yangtze River, and the genetic protection of the populations should be highlighted in the future.

The complete mitochondrial DNA of Gymnocypris potanini firmispinatus and comparative mitogenomic analyses of the Gymnocypris species.

The complete mitochondrial genome of Gymnocypris potanini firmispinatus was sequenced and compared with others Gymnocypris species. The mitochondrial genome, consisting of 16,680 base pairs (bp), encoded 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs and a non-coding control region, as those found in other Gymnocypris species. These results can provide useful data for further studies on taxonomic status, molecular systematics and stock evaluation.

The complete mitochondrial DNA of Gymnocypris potanini Herzensten and comparative mitogenomic analyses in Gymnocypris species.

In the study, the complete mitochondrial genome of Gymnocypris potanini Herzensten was sequenced and compared with other Gymnocypris species. The mitochondrial genome, consisting of 16,749 base pairs (bp), encoded 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a noncoding control region, similar as that found in other Gymnocypris species. These results can provide useful information for further studies on taxonomic status, molecular systematics, and stock evaluation.