Sex-biased miRNAs in gonad and their potential roles for testis development in yellow catfish

The potential regulatory role of the non-coding RNAs in regulating the exogenous estrogen-induced feminization in Takifugu rubripes gonad

Estrogen plays a pivotal role in the early stage of sex differentiation in teleost. However, the underlying mechanisms of estrogen-induced feminization process are still needed for further clarification. Here, the comparative analysis of whole-transcriptome RNA sequencing was conducted between 17beta-Estradiol induced feminized XY (E-XY) gonads and control gonads (C) in Takifugu rubripes. A total of 57 miRNAs, 65 lncRNAs, and 4 circRNAs were found to be expressed at lower levels in control-XY (C-XY) than that in control-XX (C-XX), and were up-regulated in XY during E2-induced feminization process. The expression levels of 24 miRNAs, and 55 lncRNAs were higher in C-XY than that in C-XX, and were down-regulated in E2-treated XY. Furthermore, a correlation analysis was performed between miRNA-seq and mRNA-seq data. In C-XX/C-XY, 114 differential expression (DE) miRNAs were predicted to target to 904 differential expression genes (DEGs), while in C-XY/E-XY, 226 DEmiRNAs were predicted to target to 2,048 DEGs. In C-XX/C-XY, and C-XY/E-XY, KEGG pathway enrichment analysis showed that those targeted genes were mainly enriched in MAPK signaling, calcium signaling, steroid hormone biosynthesis and ovarian steroidogenesis pathway. Additionally, the competitive endogenous RNA (ceRNA) regulatory network was constructed by 24 miRNAs, 21 lncRNAs, 4 circRNAs and 5 key sex-related genes. These findings suggested that the expression of critical genes in sex differentiation were altered in E2-treated XY T. rubripes may via the lncRNA-miRNA-mRNA regulation network to facilitate the differentiation and maintenance of ovaries. Our results provide a new insight into the comprehensive understanding of the effects of estrogen signaling pathways on sex differentiation in teleost gonads.

MicroRNAs are involved in ovarian physiology of greater amberjack (Seriola dumerili) under captivity

Gonad maturation is critical for the reproductive success of any organism, and in fish, captivity can significantly affect their reproductive performance, leading to maturation incompetence and spawning failure. The greater amberjack (Seriola dumerili), a fish species recently introduced to aquaculture fails to undergo oocyte maturation, ovulation, and spawning when reared in aquaculture facilities. Since confinement has been shown to influence gonad maturation and completion of the reproductive cycle, investigations into epigenetic mechanisms may shed light on the reasoning behind the reproductive dysfunctions of fish under captivity. Among the known important epigenetic regulators are small non-coding RNAs (sncRNAs), and in particular microRNAs (miRNAs). In this study, immature, maturing (late vitellogenesis), and spent ovaries of captive greater amberjack were collected, and the differential expression of miRNAs in the three different ovarian development stages was examined. Expression patterns of conserved and novel miRNAs were identified, and potential targets of highly differentially expressed miRNAs were detected. Additionally, read length distribution showed two prominent peaks in the three different ovarian maturation stages, corresponding to miRNAs and putative piwi-interacting RNAs (piRNAs), another type of ncRNAs with a germ-cell specific role. Furthermore, miRNA expression patterns and their putative target mRNAs are discussed, in relevance with the different ovarian maturation stages of captive greater amberjack. Overall, this study provides insights into the role of miRNAs in the reproductive dysfunctions observed in fish under captivity and highlights the importance of epigenetic mechanisms in understanding and managing the reproductive performance of economically important fish species.

Integrated analysis of lncRNA, miRNA and mRNA expression profiles reveals regulatory pathways associated with pig testis function

Long noncoding RNA (lncRNA) and microRNA (miRNA) are known to play pivotal roles in mammalian testicular function and spermatogenesis. However, their impact on porcine male reproduction has yet to be well unraveled. Here, we sequenced and identified lncRNA and miRNA expressed in the testes of Chinese indigenous Banna mini-pig inbred line (BMI) and introduced Western Duroc (DU) and Large White (LW) pigs. By pairwise comparison (BMI vs DU, BMI vs LW, and DU vs LW), we found the gene expression differences in the testes between Chinese local pigs and introduced Western commercial breeds were more striking than those between introduced commercial breeds. Furthermore, we found 1622 co-differentially expressed genes (co-DEGs), 122 co-differentially expressed lncRNAs (co-DELs), 39 co-differentially expressed miRNAs (co-DEMs) in BMI vs introduced commercial breeds (DU and LW). Functional analysis revealed that these co-DEGs and co-DELs/co-DEMs target genes were enriched in male sexual function pathways, including MAPK, AMPK, TGF-β/Smad, Hippo, NF-kappa B, and PI3K/Akt signaling pathways. Additionally, we established 10,536 lncRNA-mRNA, 11,248 miRNA-mRNA pairs, and 62 ceRNA (lncRNA-miRNA-mRNA) networks. The ssc-miR-1343 had the most interactive factors in the ceRNA network, including 20 mRNAs and 3 lncRNAs, consisting of 56 ceRNA pairs. These factors played extremely important roles in the regulation of testis function as key nodes in the interactive regulatory network. Our results provide insight into the functional roles of lncRNAs and miRNAs in porcine testis and offer a valuable resource for understanding the differences between Chinese indigenous and introduced Western pigs.

Circulating MicroRNAs Indicative of Sex and Stress in the European Seabass (Dicentrarchus labrax): Toward the Identification of New Biomarkers

MicroRNAs (miRNAs) constitute a new category of biomarkers. Studies on miRNAs in non-mammalian species have drastically increased in the last few years. Here, we explored the use of miRNAs as potential, poorly invasive markers, to identify sex and characterize acute stress in fish. The European seabass (Dicentrarchus labrax) was chosen as a model because of its rapid response to stress and its specific sex determination system, devoid of sexual chromosomes. We performed a small RNA-sequencing analysis in the blood plasma of male and female European seabass (mature and immature) as well as in the blood plasma of juveniles submitted to an acute stress and sampled throughout the recovery period (at 0 h, 0.5 h, 1.5 h and 6 h). In immature individuals, both miR-1388-3p and miR-7132a-5p were up-regulated in females, while miR-499a-5p was more abundant in males. However, no miRNAs were found to be differentially expressed between sexes in the blood plasma of mature individuals. For the acute stress analysis, five miRNAs (miR-155-5p, miR-200a-3p, miR-205-1-5p, miR-143-3p, and miR-223-3p) followed cortisol production over time. All miRNAs identified were tested and validated by RT-qPCR on sequenced samples. A complementary analysis on the 3'UTR sequences of the European seabass allowed to predict potential mRNA targets, some of them being particularly relevant regarding stress regulation, e.g., the glucocorticoid receptor 1 and the mineralocorticoid receptor. The present study provides new avenues and recommendations on the use of miRNAs as biomarkers of sex or stress of the European seabass, with potential application on other fish species.

High production of transfer RNAs identifies the presence of developing oocytes in ovaries and intersex testes of teleost fish

5S rRNA is highly transcribed in fish oocytes and this transcription levels can be used to identify the presence of oocytes in the intersex testes of fish exposed to xenoestrogens. Similar to 5S rRNA, tRNAs are transcribed by RNA polymerase III (Pol-III) in eukaryotes, so this study focuses in the analysis of the levels of expression of tRNAs in the gonads (ovaries and testes) of eight teleost species as a possible new oocyte molecular marker. Total RNA extracted from gonads of six commercial teleost species in the Biscay Bay, from the pollution sentinel species thicklip grey mullet (Chelon labrosus) known present intersex testes in response to xenoestrogens in Gernika estuary and from the laboratory model species Danio rerio were analysed through capillary electrophoresis. Bioanalyzer electropherograms were used to quantify the concentrations of tRNAs, 5S and 5.8S rRNA. All studied ovaries expressed significantly higher levels of tRNAs and 5S rRNA than testes. A tRNA to 5.8S rRNA index was calculated which differentiates ovaries from testes, and identifies some intersex testes in between testes and ovaries in mullets. The tRNA/5.8S ratio was highest in ovaries in previtellogenic stage, decreasing towards maturity. Thus, strong oocyte expression of tRNAs is an additional proof of high activity levels of Pol-III during early stages of oocyte development in teleost ovaries. Incidentally, we observed that miRNA concentrations were always higher in testes than ovaries. The indexing approach developed in the present study could have multiple applications in teleost reproduction research and in the development of early molecular markers of intersex condition.

The Male-Biased Expression of miR-2954 Is Involved in the Male Pathway of Chicken Sex Differentiation

Many expression data showed miRNAs have a potential function on regulating gonadal differentiation in animals, but their function is rarely studied in vivo, especially in chickens. Using the comprehensive expression profiles analysis, the specific male-biased miR-2954, which is significantly higher expressed in male embryos and gonads at all detected stages, was firstly screened during the early stages of chicken embryogenesis and gonadogenesis. In sex-reversed female gonads treated with aromatase inhibitors, the expression of miR-2954 was increased, which was consistent with the up-regulation of DMRT1 and SOX9. The injection of vivo-morpholino of miR-2954 significantly inhibited the expression of miR-2954 in chicken embryos, and the down-regulation of miR-2954 decreased the expression of testis-associated genes DMRT1 and SOX9, while the expression of ovary-associated genes and the gonadal morphology did not change obviously. These results confirm that miR-2954 coincides with testicular differentiation in chicken embryos, but whether it might be an upstream cell autonomous factor to sex development in birds still need to be further determined.

Comparative transcriptome profiles of four sexually size dimorphic fish

Sexual size dimorphism is widespread in fish species. Although sex growth differences in multiple species have been studied successively, the commonalities of regulatory mechanisms across sexually dimorphic species are unknown. In this study, we performed RNA-seq analysis of four representative fish (loach, half-smooth tongue sole, yellow catfish, and Nile tilapia) with significant growth differences between females and males. Clean reads were identified from four fish species, ranging from 45,718,052 to 57,733,120. Following comparison transcriptome analysis, there were 1,132 and 1,108, 1,290 and 1,102, 4,732 and 4,266, 748 and 192 differentially expressed genes (DEGs) in the brain and muscle of loach, half-smooth tongue sole, yellow catfish, and Nile tilapia, respectively. Furthermore, the expression levels were validated by quantitative real-time PCR (qRT-PCR). Comparative transcriptome profiles of four fish described here will provide fundamental information for further studies on the commonalities of sexually size dimorphic fish in regulating growth differences between females and males.

Identification and Characterization of Sex-Biased miRNAs in the Golden Pompano (Trachinotus blochii)

The golden pompano (Trachinotus blochii) is a marine fish of considerable commercial importance in China. It shows notable sexual size dimorphism; the growth rate of females is faster than that of males. Therefore, sex-biased research is of great importance in T. blochii breeding. However, there have been few studies on sex differentiation and mechanisms underlying sex determination in T. blochii. MicroRNAs (miRNAs) play crucial roles in sex differentiation and determination in animals. However, limited miRNA data are available on fish. In this study, two small RNA libraries prepared from the gonads of T. blochii were constructed and sequenced. The RNA-seq analysis yielded 1366 known and 69 novel miRNAs with 289 significantly differentially expressed miRNAs (p < 0.05). Gene ontology (GO) analysis confirmed that the TFIIA transcription factor complex (GO: 0005672) was the most significantly enriched GO term. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis showed that the differentially expressed miRNAs and target genes were mainly related to sex determination and gonadal developmental signaling pathways, specifically the Wnt signaling pathway, MAPK signaling pathway, and steroid biosynthetic pathway. MiRNA-mRNA co-expression network analysis strongly suggested a role for sex-biased miRNAs in sex determination/differentiation and gonadal development. For example, gata4, foxo3, wt1, and sf1 genes were found to be regulated by bta-miR-2898; esr2 and foxo3 by novel_176, and ar by oar-let-7b. Quantitative real-time polymerase chain reaction analysis of selected mRNAs and miRNAs validated the integrated analysis. This study established a set of sex-biased miRNAs that are potential regulatory factors in gonadal development in T. blochii. These results provide new insight into the function of miRNAs in sex differentiation and determination in T. blochii and highlight some key miRNAs for future studies.

Small RNA sequencing reveals sex-related miRNAs in Collichthys lucidus

Collichthys lucidus (C. lucidus) is an economically important fish species, exhibiting sexual dimorphism in its growth rate. However, there is a lack of research on its underlying sex-related mechanisms. Therefore, small RNA sequencing was performed to better comprehend these sex-related molecular mechanisms. In total, 171 differentially expressed miRNAs (DE-miRNAs) were identified between the ovaries and testes. Functional enrichment analysis revealed that the target genes of DE-miRNAs were considerably enriched in the p53 signaling, PI3K–Akt signaling, and TGF-beta signaling pathways. In addition, sex-related miRNAs were identified, and the expression of miR-430c-3p and miR-430f-3p was specifically observed in the gonads compared with other organs and their expression was markedly upregulated in the testes relative to the ovaries. Bmp15 was a target of miR-430c-3p and was greatly expressed in the ovaries compared with the testes. Importantly, miR-430c-3p and bmp15 co-expressed in the ovaries and testes. This research provides the first detailed miRNA profiles for C. lucidus concerning sex, likely laying the basis for further studies on sex differentiation in C. lucidus.

Tandem Mass Tag-Based Quantitative Proteomics Analysis of Gonads Reveals New Insight into Sexual Reversal Mechanism in Chinese Soft-Shelled Turtles

Chinese soft-shelled turtles display obvious sex dimorphism. The exogenous application of hormones (estradiol and methyltestosterone) can change the direction of gonadal differentiation of P. sinensis to produce sex reversed individuals. However, the molecular mechanism remains unclear. In this study, TMT-based quantitative proteomics analysis of four types of P. sinensis (female, male, pseudo-female, and pseudo-male) gonads were compared. Quantitative analysis of 6107 labeled proteins in the four types of P. sinensis gonads was performed. We identified 440 downregulated and 423 upregulated proteins between pseudo-females and males, as well as 394 downregulated and 959 upregulated proteins between pseudo-males and females. In the two comparisons, the differentially expressed proteins, including K7FKG1, K7GIQ2, COL4A6, K7F2U2, and K7FF80, were enriched in some important pathways, such as focal adhesion, endocytosis, apoptosis, extracellular matrix-receptor interaction, and the regulation of actin cytoskeleton, which were upregulated in pseudo-female vs. male and downregulated in pseudo-male vs. female. In pathways such as ribosome and spliceosome, the levels of RPL28, SRSF3, SNRNP40, and HNRNPK were increased from male to pseudo-female, while they decreased from female to pseudo-male. All differentially expressed proteins after sexual reversal were divided into six clusters, according to their altered levels in the four types of P. sinensis, and associated with cellular processes, such as embryonic development and catabolic process, that were closely related to sexual reversal. These data will provide clues for the sexual reversal mechanism in P. sinensis.

The miR-200 Family Targeting amh Affects the Gonadal Development of Japanese Flounder

Four members of the miR-200 family in Japanese flounder (Paralichthys olivaceus) have sex-biased expression patterns, but their target genes and how they work in the development of the gonads are rarely known. Anti-Müllerian hormone (AMH) can inhibit the development of Muller’s duct in female mammals and regulate the formation of gametes after sexual maturity. There is no Muller’s duct in teleosts, but the amh gene still exists. Knockout of amh results in sex reversal from male to female. Therefore, it is essential to explore the relationship between the miR-200 family and amh to clarify what role miR-200 plays in the development of the gonads. In Japanese flounder, the two binding sites for the miR-200 family in the 3′UTR of amh were found through bioinformatic prediction. Double luciferase and green fluorescent protein reporter experiments demonstrated amh to be directly targeted by miR-200a and miR-200b. Moreover, miR-200a and miR-200b reduced the expression of amh through site 1 rather than site 2. To explore the regulatory role of miR-200a in gonadal development, we further overexpressed miR-200a in the primary Sertoli cells of the testis. With the overexpression of miR-200a, the expression of amh decreased, while the expression of the other two male sex-related genes, dmrt1 (doublesex and mab-3 related transcription factor 1) and gsdf (diagonal soma driven factor), increased significantly. This result indicates that the miR-200 family regulates the gonadal differentiation and development by targeting amh in Japanese flounder.

Integrated analysis of the miRNA and mRNA expression profiles in Leiocassis longirostris at gonadal maturation

Leiocassis longirostris is a commercially important fish species that shows a sexually dimorphic growth pattern. A lack of molecular data from the gonads of this species has hindered research and selective breeding efforts. In this study, we conducted a comprehensive analysis of the expression profile of miRNA and mRNA to explore their regulatory roles in the gonadal maturation stage of L. longirostris. We identified 60 differentially expressed miRNAs and 20,752 differentially expressed genes by sequencing. A total of 90 miRNAs and 21 target genes involved in gonad development and sex determination were identified. Overall, the results of this study enhance our understanding of the molecular mechanisms underlying sex determination and differentiation and provide valuable genomic information for the selective breeding of L. longirostris.

Sex-biased dynamics of three-spined stickleback (Gasterosteus aculeatus) gene expression patterns

The study of the differences between sexes presents an excellent model to unravel how phenotypic variation is achieved from a similar genetic background. Sticklebacks are of particular interest since evidence of a heteromorphic chromosome pair has not always been detected. The present study investigated sex-biased mRNA and small non-coding RNA (sncRNA) expression patterns in the brain, adipose tissues, and gonads of the three-spined stickleback. The sncRNA analysis indicated that regulatory functions occurred mainly in the gonads. Alleged miRNA-mRNA interactions were established and a mapping bias of differential expressed transcripts towards chromosome 19 was observed. Key players previously shown to control sex determination and differentiation in other fish species but also genes like gapdh were among the transcripts identified. This is the first report in the three-spined stickleback demonstrating tissue-specific expression comprising both mRNA and sncRNA between sexes, emphasizing the importance of mRNA-miRNA interactions as well as new presumed genes not yet identified to have gender-specific roles.

MicroRNAs May Play an Important Role in Sexual Reversal Process of Chinese Soft-Shelled Turtle, Pelodiscus sinensis

The Chinese soft-shelled (Pelodiscus sinensis) turtle exhibits obvious sex dimorphism, which leads to the higher economic and nutritional value of male individuals. Exogenous hormones can cause the transformation from male to female phenotype during gonadal differentiation. However, the molecular mechanism related to the sexual reversal process is unclear. In this study, we compared the difference between the small RNAs of male, female, and pseudo-female turtles by small RNA-seq to understand the sexual reversal process of Chinese soft-shelled turtles. A certain dose of estrogen can cause the transformation of Chinese soft-shelled turtles from male to female, which are called pseudo-female individuals. The result of small RNA-seq has revealed that the characteristics of pseudo-females are very similar to females, but are strikingly different from males. The number of the microRNAs (miRNAs) of male individuals was significantly less than the number of female individuals or pseudo-female individuals, while the expression level of miRNAs of male individuals were significantly higher than the other two types. Furthermore, we found 533 differentially expressed miRNAs, including 173 up-regulated miRNAs and 360 down-regulated miRNAs, in the process of transformation from male to female phenotype. Cluster analysis of the total 602 differential miRNAs among females, males, and pseudo-females showed that miRNAs played a crucial role during the sexual differentiation. Among these differential miRNAs, we found 12 miRNAs related to gonadal development and verified their expression by qPCR. The TR-qPCR results confirmed the differential expression of 6 of the 12 miRNAs: miR-26a-5p, miR-212-5p, miR-202-5p, miR-301a, miR-181b-3p and miR-96-5p were involved in sexual reversal process, which was consistent with the results of omics. Using these six miRNAs and some of their target genes, we constructed a network diagram related to gonadal development. We suggest that these miRNAs may play an important role in the process of effective sex reversal, which would contribute to the breeding of all male strains of Chinese soft-shelled turtles.

Advances in Reproductive Endocrinology and Neuroendocrine Research Using Catfish Models

Catfishes, belonging to the order siluriformes, represent one of the largest groups of freshwater fishes with more than 4000 species and almost 12% of teleostean population. Due to their worldwide distribution and diversity, catfishes are interesting models for ecologists and evolutionary biologists. Incidentally, catfish emerged as an excellent animal model for aquaculture research because of economic importance, availability, disease resistance, adaptability to artificial spawning, handling, culture, high fecundity, hatchability, hypoxia tolerance and their ability to acclimate to laboratory conditions. Reproductive system in catfish is orchestrated by complex network of nervous, endocrine system and environmental factors during gonadal growth as well as recrudescence. Lot of new information on the molecular mechanism of gonadal development have been obtained over several decades which are evident from significant number of scientific publications pertaining to reproductive biology and neuroendocrine research in catfish. This review aims to synthesize key findings and compile highly relevant aspects on how catfish can offer insight into fundamental mechanisms of all the areas of reproduction and its neuroendocrine regulation, from gametogenesis to spawning including seasonal reproductive cycle. In addition, the state-of-knowledge surrounding gonadal development and neuroendocrine control of gonadal sex differentiation in catfish are comprehensively summarized in comparison with other fish models.

Integrated Analysis of miR-430 on Steroidogenesis-Related Gene Expression of Larval Rice Field Eel Monopterus albus

The present study aims to reveal the mechanism by which miR-430s regulate steroidogenesis in larval rice field eel Monopterus albus. To this end, M. albus embryos were respectively microinjected with miRNA-overexpressing mimics (agomir430a, agomir430b, and agomir430c) or miRNA-knockdown inhibitors (antagomir430a, antagomir430b, and antagomir430c). Transcriptome profiling of the larvae indicated that a total of more than 149 differentially expressed genes (DEGs) were identified among the eight treatments. Specifically, DEGs related to steroidogenesis, the GnRH signaling pathway, the erbB signaling pathway, the Wnt signaling pathway, and other pathways were characterized in the transcriptome. We found that steroidogenesis-related genes (hydroxysteroid 17-beta dehydrogenase 3 (17β-hsdb3), hydroxysteroid 17-beta dehydrogenase 7 (17β-hsdb7), hydroxysteroid 17-beta dehydrogenase 12 (17β-hsdb12), and cytochrome P450 family 19 subfamily a (cyp19a1b)) were significantly downregulated in miR-430 knockdown groups. The differential expressions of miR-430 in three gonads indicated different roles of three miR-430 (a, b, and c) isoforms in regulating steroidogenesis and sex differentiation. Mutation of the miR-430 sites reversed the downregulation of cytochrome P450 family 17 (cyp17), cyp19a1b, and forkhead box L2 (foxl2) reporter activities by miR-430, indicating that miR-430 directly interacted with cyp17, cyp19a1b, and foxl2 genes to inhibit their expressions. Combining these findings, we concluded that miR-430 regulated the steroidogenesis and the biosynthesis of steroid hormones by targeting cyp19a1b in larval M. albus. Our results provide a novel insight into steroidogenesis at the early stage of fish at the molecular level.

miR-133b targets tagln2 and functions in tilapia oogenesis

microRNAs (miRNAs) are important components of non-coding RNAs that participate in diverse life activities by regulating gene expression at the post transcriptional level through base complementary pairing with 3'UTRs of target mRNAs. miR-133b is a member of the miR-133 family, which play important roles in muscle differentiation and tumorigenesis. Recently, miR-133b was reported to affect estrogen synthesis by targeting foxl2 in mouse, while its role in fish reproduction remains to be elucidated. In the present study, we isolated the complete sequence of miR-133b, which was highly expressed in tilapia ovary at 30 and 90 dah (days after hatching) and subsequently decreased at 120 to 150 dah by qPCR. Interestingly, only a few oogonia were remained in the antagomir-133b treated tilapia ovary, while phase I and II oocytes were observed in the ovaries of the control group. Unexpectedly, the expression of foxl2 and cyp19a1a, as well as estradiol levels in serum were increased in the treated group. Furthermore, tagln2, an important factor for oogenesis, was predicted as the target gene of miR-133b, which was confirmed by dual luciferase reporter vector experiments. miR-133b and tagln2 were co-expressed in tilapia ovaries. Taken together, miR-133b may be involved in the early oogenesis of tilapia by regulating tagln2 expression. This study enriches the understanding of miR-133b function during oogenesis and lays a foundation for further study of the regulatory network during oogenesis.

Sex bias miRNAs in Cynoglossus semilaevis could play a role in transgenerational inheritance

Alterations of non-coding RNA profiling in spermatozoa are candidate mechanisms related to changes in paternal environment and progeny. Transgenerational inheritance of sex in pseudomales of Cynoglossus semilaevis, a fish with significant sex dimorphism, is a typical example of non-Mendelian inheritance. In the present study, miRNA profiles of spermatozoa were compared between male and pseudomale of C. semilaevis. Differential miRNAs in sperm from F0 and F1 generation also provides clues for revealing the possible role of non-coding RNA mediated transgenerational inheritance. Four sexual bias miRNAs, dre-miR-26a-5p, dre-miR-27b-3p, dre-miR-125b-5p,pol-199a-5p, were identified and verified in F0 and F1 generation of C. semilaevis. All of them were highly expressed in male sperm compared with pseudomale sperm. Function of target genes indicates that target genes of these differential RNAs are highly correlated with sex differentiation, gametogenesis and maintenance of secondary sexual characteristics. In a word, identification of epigenetic markers in gametes has great prospects in predicting susceptibility and properties in offsprings, and providing an indicator of parentalgenetic property.

Identification of sex differentiation-related microRNA and long non-coding RNA in Takifugu rubripes gonads

Although sex determination and differentiation are key developmental processes in animals, the involvement of non-coding RNA in the regulation of this process is still not clarified. The tiger pufferfish (Takifugu rubripes) is one of the most economically important marine cultured species in Asia, but analyses of miRNA and long non-coding RNA (lncRNA) at early sex differentiation stages have not been conducted yet. In our study, high-throughput sequencing technology was used to sequence transcriptome libraries from undifferentiated gonads of T. rubripes. In total, 231 (107 conserved, and 124 novel) miRNAs were obtained, while 2774 (523 conserved, and 2251 novel) lncRNAs were identified. Of these, several miRNAs and lncRNAs were predicted to be the regulators of the expression of sex-related genes (including fru-miR-15b/foxl2, novel-167, novel-318, and novel-538/dmrt1, novel-548/amh, lnc_000338, lnc_000690, lnc_000370, XLOC_021951, and XR_965485.1/gsdf). Analysis of differentially expressed miRNAs and lncRNAs showed that three mature miRNAs up-regulated and five mature miRNAs were down-regulated in male gonads compared to female gonads, while 79 lncRNAs were up-regulated and 51 were down-regulated. These findings could highlight a group of interesting miRNAs and lncRNAs for future studies and may reveal new insights into the function of miRNAs and lncRNAs in sex determination and differentiation.

Non-coding RNA Expression Patterns of Two Different Teleost Gonad Maturation Stages

Non-coding RNAs (ncRNAs) are involved in several different regulatory pathways including reproduction. In teleost fish, efficacious reproduction is heavily dependent on the completion of the reproductive cycle. The presence of ncRNA, however, and their expression dynamics and putative regulatory role in mature and immature gonads have not yet been extensively explored. Therefore, the abundance of ncRNAs in mature and immature female sharpsnout seabream (Diplodus puntazzo) was investigated. The sharpsnout seabream is a rudimentary hermaphrodite which, in captivity, displays dysfunctions in the gonad maturation process. Our analyses revealed a gonad specific read length distribution with two main peaks representing miRNAs (21-26 nt) and PIWI RNA (27-34 nt). Besides, distinct expression patterns for several ncRNA biotypes including microRNAs (miRNAs), PIWI RNAs (piRNAs), and ribosomal RNAs (rRNAs) were detected. Identified miRNA accounted to 938, corresponding to ~ 13% of obtained transcripts. Among the differential expressed ncRNAs, 10 (~ 7%) were annotated as miRNA, out of which 2 were found in higher abundance in immature gonads (miR-125c and miR-24) and 8 (miR-451, miR-7a, miR-122-1, miR190a, miR129, ENSGACT00000029608, ENSGACT00000029489, and ENSGACT00000029667) were found to be higher expressed in mature gonads. Putative miRNA targets, including long non-coding RNAs (lncRNAs) and genes, are proposed. Target genes are involved in several processes of fish oocyte development, such as steroidogenesis, proteolysis, and apoptosis, and may explain hormone regulation. This study demonstrates a gonad maturation biased ncRNA profile which in turn may support the role of ncRNAs in ovarian physiology and reproductive performance of fish, stressing the specific function of each RNA biotype in oocyte development.

A complement factor I (CFI) gene mediates innate immune responses in yellow catfish Pelteobagrus fulvidraco

This study isolated CFI gene from Pelteobagrus fulvidraco and named it PfCFI. The cDNA of PfCFI is 2374 bp long, including a 52 bp 5' untranslated sequence, a 222 bp 3' untranslated sequence, and an open reading frame (ORF) of 2100 bp encoding polypeptide consisting of 699 amino acids. Phylogenetic analysis revealed that the PfCFI was closely related to CFI of Ictalurus punctatus. Real-time quantitative reverse transcription-PCR (qRT-PCR) analysis indicate that there is the PfCFI gene which expressed in all the rest of tested tissues in varied levels, and mainly distributed in liver and least in heart. The reseachers induce the expressions level of PfCFI gene in liver, spleen, head kidney and blood at different points in time after challenged with lipopolysaccharide (LPS), and polyriboinosinic polyribocytidylic acid (poly I:C), respectively. Together these results suggested that CFI gene plays an important role in resistance to pathogens in yellow catfish immunity.

Gonad-Specific Transcriptomes Reveal Differential Expression of Gene and miRNA Between Male and Female of the Discus Fish (Symphysodon aequifasciatus)

The discus fish (Symphysodon aequifasciatus) is an ornamental fish that is well-known around the world. In artificial reproduction, they must be matched by one male and one female, whereas phenotype investigation indicated that there are no significant differences in appearance between males and females, which causes great difficulties in the mating during artificial reproduction. So, it is of great importance to establish artificial sex identification methods for the discus fish. The molecular mechanism of the sexual dimorphism of the discus fish was previously unknown. In this study, we constructed six cDNA libraries from three adult testes and three adult ovaries and performed RNA sequencing for identifying sex-biased candidate genes and microRNAs (miRNAs). A total of 50,082 non-redundant genes (unigenes) were identified, of which 18,570 unigenes were significantly overexpressed in testes, and 11,182 unigenes were significantly overexpressed in ovaries. A total of 551 miRNAs were identified, of which 47 miRNAs were differentially expressed between testes and ovaries. Eight differentially expressed unigenes, seven differentially expressed miRNAs and one non-differential miRNA were validated by quantitative real-time polymerase chain reaction. Twenty-four of these differentially expressed miRNAs and their 15 predicted target genes constituted 41 miRNA–mRNA interaction pairs, and some of vital sex-related metabolic pathways were also identified. These results revealed these differentially expressed genes and miRNAs between ovary and testis might be involved in regulating gonadal development, sex determination, gametogenesis, and physiological function maintenance, and there are complex regulatory networks between genes and miRNAs. It can help us understand the molecular mechanism of the sexual dimorphism and obtain a high-efficiency sex identification method in the artificial reproduction process of the discus fish.

Identification and Comparison of microRNAs in the Gonad of the Yellowfin Seabream (Acanthopagrus Latus)

Yellowfin seabream (Acanthopagrus latus) is a commercially important fish in Asian coastal waters. Although natural sex reversal has been described in yellowfin seabream, the mechanisms underlying sexual differentiation and gonadal development in this species remain unclear. MicroRNAs (miRNAs) have been shown to play crucial roles in gametogenesis and gonadal development. Here, two libraries of small RNAs, constructed from the testes and ovaries of yellowfin seabream, were sequenced. Across both gonads, we identified 324 conserved miRNAs and 92 novel miRNAs: 67 ovary-biased miRNAs, including the miR-200 families, the miR-29 families, miR-21, and miR-725; and 88 testis-biased miRNAs, including the let-7 families, the miR-10 families, miR-7, miR-9, and miR-202-3p. GO (Gene Ontology) annotations and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analyses of putative target genes indicated that many target genes were significantly enriched in the steroid biosynthesis pathway and in the reproductive process. Our integrated miRNA-mRNA analysis demonstrated a putative negatively correlated expression pattern in yellowfin seabream gonads. This study profiled the expression patterns of sex-biased miRNAs in yellowfin seabream gonads, and provided important molecular resources that will help to clarify the miRNA-mediated post-transcriptional regulation of sexual differentiation and gonadal development in this species.

Heterogenic Origin of Micro RNAs in Atlantic Salmon (Salmo salar) Seminal Plasma

The origin and contribution of seminal plasma RNAs into the whole semen RNA repertoire are poorly known, frequently being overlooked or neglected. In this study, we used high-throughput sequencing and RT-qPCR to profile microRNA (miRNA) constituents in the whole semen, as well as in fractionated spermatozoa and seminal plasma of Atlantic salmon (Salmo salar). We found 85 differentially accumulated miRNAs between spermatozoa and the seminal plasma. We identified a number of seminal plasma-enriched and spermatozoa-enriched miRNAs. We localized the expression of some miRNAs in juvenile and mature testes. Two abundant miRNAs, miR-92a-3p and miR-202-5p, localized to both spermatogonia and somatic supporting cells in immature testis, and they were also highly abundant in somatic cells in mature testis. miR-15c-5p, miR-30d-5p, miR-93a-5p, and miR-730-5p were detected only in mature testis. miRs 92a-3p, 202-5p, 15c-5p, and 30d-5p were also detected in a juvenile ovary. The RT-qPCR experiment demonstrated lack of correlation in miRNA transcript levels in seminal plasma versus blood plasma. Our results indicate that salmon semen is rich in miRNAs, which are present in both spermatozoa and seminal plasma. Testicular-supporting somatic cells are likely the source of seminal plasma enrichment, whereas blood plasma is unlikely to contribute to the seminal plasma miRNA repertoire.

An estradiol-17β/miRNA-26a/cyp19a1a regulatory feedback loop in the protogynous hermaphroditic fish, Epinephelus coioides

Cyp19a1a is a key gene responsible for the production of estradiol-17β (E₂), the main functional estrogen and a major downstream regulator of reproduction in teleost fish. It is widely known that CYP19 gene expression, aromatase activity, and E₂ production can influence gonadal differentiation and sex reversal in teleost fish, but the feedback mechanisms whereby E₂ regulates cyp19a1a remain poorly understood, especially regarding the potential roles of endogenous small RNA molecules (miRNAs). Here, we identified miR-26a-5p as a regulatory factor of its predicted target gene (cyp19a1a). In vitro and in vivo studies showed that miR-26a-5p can decrease cyp19a1a expression. Furthermore, high doses of E₂ act as a repressor of miR-26a-5p. This study proposes a regulatory feedback loop whereby E₂ regulates cyp19a1a through miR-26a-5p, and suggests that this positive feedback is an important aspect of the control of E₂ production.

De novo transcriptome assembly and sex-biased gene expression in the gonads of Amur catfish (Silurus asotus)

Amur catfish is extensively distributed and cultured in Asian countries. Despite of economic importance, the genomic information of this species remains limited. A reference transcriptome of Amur catfish was assembled and the sex-biased gene expression in the gonads was characterized using RNA-sequencing. The assembled transcriptome of Amur catfish consisted of 74,840 transcripts. The N50, mean length and max length of transcripts are 1970, 1235 and 16,748 bp. Putative sex-specific transcripts were identified and sex-specific expression of the representative genes was verified by RT-PCR. Differential expression analysis identified 5401 ovary-biased and 5618 testis-biased genes. The ovary-biased genes were mainly enriched in pathways such as RNA transport and ribosome biogenesis in eukaryotes. The testis-biased genes were enriched in calcium signaling and cytokine-cytokine receptor interaction, etc. Our data provide a valuable genomic resource for further investigating the genetic basis of sex determination, sex differentiation and sexual dimorphism of catfish.

Twenty years of transcriptomics, 17alpha-ethinylestradiol, and fish

In aquatic toxicology, perhaps no pharmaceutical has been investigated more intensely than 17alpha-ethinylestradiol (EE2), the active ingredient of the birth control pill. At the turn of the century, the fields of comparative endocrinology and endocrine disruption research witnessed the emergence of omics technologies, which were rapidly adapted to characterize potential hazards associated with exposures to environmental estrogens, such as EE2. Since then, significant advances have been made by the scientific community, and as a result, much has been learned about estrogen receptor signaling in fish from environmental xenoestrogens. Vitellogenin, the egg yolk precursor protein, was identified as a major estrogen-responsive gene, establishing itself as the premier biomarker for estrogenic exposures. Omics studies have identified a plethora of estrogen responsive genes, contributing to a wealth of knowledge on estrogen-mediated regulatory networks in teleosts. There have been ~40 studies that report on transcriptome responses to EE2 in a variety of fish species (e.g., zebrafish, fathead minnows, rainbow trout, pipefish, mummichog, stickleback, cod, and others). Data on the liver and testis transcriptomes dominate in the literature and have been the subject of many EE2 studies, yet there remain knowledge gaps for other tissues, such as the spleen, kidney, and pituitary. Inter-laboratory genomics studies have revealed transcriptional networks altered by EE2 treatment in the liver; networks related to amino acid activation and protein folding are increased by EE2 while those related to xenobiotic metabolism, immune system, circulation, and triglyceride storage are suppressed. EE2-responsive networks in other tissues are not as comprehensively defined which is a knowledge gap as regulated networks are expected to be tissue-specific. On the horizon, omics studies for estrogen-mediated effects in fish include: (1) Establishing conceptual frameworks for incorporating estrogen-responsive networks into environmental monitoring programs; (2) Leveraging in vitro and computational toxicology approaches to identify chemicals associated with estrogen receptor-mediated effects in fish (e.g., male vitellogenin production); (3) Discovering new tissue-specific estrogen receptor signaling pathways in fish; and (4) Developing quantitative adverse outcome pathway predictive models for estrogen signaling. As we look ahead, research into EE2 over the past several decades can serve as a template for the array of hormones and endocrine active substances yet to be fully characterized or discovered.

The integrative knowledge base for miRNA-mRNA expression in colorectal cancer

"miRNA colorectal cancer" (https://mirna-coadread.omics.si/) is a freely available web application for studying microRNA and mRNA expression and their correlation in colorectal cancer. To the best of our knowledge, "miRNA colorectal cancer" has the largest knowledge base of miRNA-target gene expressions and correlations in colorectal cancer, based on the largest available sample size from the same source of data. Data from high-throughput molecular profiling of 295 colon and rectum adenocarcinoma samples from The Cancer Genome Atlas was analyzed and integrated into our knowledge base. The objective of developing this web application was to help researchers to discover the behavior and role of miRNA-target gene interactions in colorectal cancer. For this purpose, results of differential expression and correlation analyses of miRNA and mRNA data collected in our knowledge base are available through web forms. To validate our knowledge base experimentally, we selected genes FN1, TGFB2, RND3, ZEB1 and ZEB2 and miRNAs hsa-miR-200a/b/c-3p, hsa-miR-141-3p and hsa-miR-429. Both approaches revealed a negative correlation between miRNA hsa-miR-200b/c-3p and its target gene FN1 and between hsa-miR-200a-3p and its target TGFB2, thus supporting the usefulness of the developed knowledge base.

A non-mammalian Toll-like receptor 26 (TLR26)gene mediates innate immune responses in yellow catfish Pelteobagrus fulvidraco

In this study, we identified a fish-specific Toll-like receptor (TLR) in Pelteobagrus fulvidraco, an economically important freshwater fish in China. This TLR, PfTLR26, was shown to be encoded by a 3084 bp open reading frame (ORF), producing a polypeptide 1027 amino acids in length. The PfTLR26 protein contains a signal peptide, eight leucine-rich repeat (LRR) domains, two LRR_TYP domains in the extracellular region, and a Toll/interleukin (IL)-1 receptor (TIR) domain in the cytoplasmic region, consistent with the characteristic TLR domain architecture. This predicted 117.1 kDa protein was highly homologous to those of other fish, with phylogenetic analysis revealing the closest relation to TLR26 of Ictalurus punctatus. Real-time quantitative reverse transcription-PCR (qRT-PCR) analysis showed that the PfTLR26 gene was expressed in all tissues tested, with the highest expression levels seen in the head kidney and blood, and the lowest seen in muscle. PfTLR26 exhibited significant upregulation in liver, spleen, head kidney, and blood at different time points following challenge with the common TLR agonists lipopolysaccharide (LPS) and polyriboinosinic polyribocytidylic acid (Poly I:C). Taken together, these results suggest that PfTLR26 may be an important component of the P. fulvidraco innate immune system, participating in the transduction of TLR signaling under pathogen stimulation.

Gonadal microRNA Expression Profiles and Their Potential Role in Sex Differentiation and Gonadal Maturation of Mud Crab Scylla paramamosain

Although the sexual dimorphism in terms of gonadal development and gametogenesis of mud crab has been described, the internal regulating mechanism and sex differentiation process remain unclear. A comparative gonadal miRNA transcriptomic study was conducted to identify miRNAs that are differentially expressed between testes and ovaries, and potentially uncover miRNAs that might be involved in sex differentiation and gonadal maturation mechanisms of mud crabs (Scylla paramamosain). A total of 10 known miRNAs and 130 novel miRNAs were identified, among which 54 were differentially expressed. Target gene prediction revealed a significant enrichment in 30 KEGG pathways, including some reproduction-related pathways, e.g. phosphatidylinositol signalling system and inositol phosphate metabolism pathways. Further analysis on six differentially expressed known miRNAs, six differentially expressed novel miRNAs and their reproduction-related putative target genes shows that both miRNAs and putative target genes showed stage-specific expression during gonadal maturation, suggesting their potential regulatory roles in sex differentiation and reproductive development. This study reveals the sex-biased miRNA profile and establishes a solid foundation for understanding the sex differentiation and gonadal maturation mechanisms of S. paramamosain.

Seminal Plasma Exosomes: Promising Biomarkers for Identification of Male and Pseudo-Males in Cynoglossus semilaevis

In mammals, small RNAs enclosed in exosomes have been identified as appropriate signatures for disease diagnosis. However, there is limited information on exosomes derived from seminal plasma, and few studies have reported analyzed the composition of exosomes and enclosed small RNAs in fish. The half-smooth tongue sole (Cynoglossus semilaevis) is an economically important fish for aquaculture, and it exhibits sexual dimorphism: the female gender show higher growth rates and larger body sizes than males. Standard karyotype analysis and artificial gynogenesis tests have revealed that this species uses heterogametic sex determination (ZW/ZZ), and so-called sex-reversed pseudo-males exist. In this study, we successfully identified exosomes in the seminal plasma of C. semilaevis; to the best of our knowledge, this is the first report of exosomes in fish seminal plasma. Analysis of the nucleotide composition showed that miRNAs were dominant in the exosomes, and the miRNAs were sequenced and compared to identify signature miRNAs as sexual biomarkers. Moreover, target genes of the signature miRNAs were predicted by sequence matching and annotation. Finally, four miRNAs (dre-miR-141-3P, dre-miR-10d-5p, ssa-miR-27b-3p, and ssa-miR-23a-3p) with significant differential expression in the males and pseudo-males were selected from the signature candidate miRNAs as markers for sex identification, and their expression profiles were verified using real-time quantitative PCR. Our findings could provide an effective detection method for sex differentiation in fish.

Characterization of the salmon louse Lepeophtheirus salmonis miRNome: Sex-biased differences related to the coding and non-coding RNA interplay

The salmon louse Lepeophtheirus salmonis is a marine ectoparasite that has a detrimental impact on salmon farms. Genomic knowledge of adult stages is critical to understand the reproductive success and lifecycle completion of this species. Here, we report a comprehensive characterization of the L. salmonis miRNome with emphasis on the sex-differences of the parasite. Small-RNA sequencing was conducted on males and females, and mRNA-sequencing was also conducted to identify miRNA-targets at these stages. Based on bioinformatics analyses, 3101 putative miRNAs were found in L. salmonis, including precursors and variants. The most abundant and over-expressed miRNAs belonged to the bantam, mir-100, mir-1, mir-263a and mir-276 families, while the most differentially expressed mRNAs corresponded to genes related to reproduction and other biological processes involved in cell-differentiation. Target analyses revealed that the most up-regulated miRNAs in males can act by inhibiting the expression of genes related to female differentiation such as vitellogenin genes. Target prediction and expression patterns suggested a pivotal role of miRNAs in the reproductive development of L. salmonis.

Whole-Genome Sequencing of Chinese Yellow Catfish Provides a Valuable Genetic Resource for High-Throughput Identification of Toxin Genes

Naturally derived toxins from animals are good raw materials for drug development. As a representative venomous teleost, Chinese yellow catfish (Pelteobagrus fulvidraco) can provide valuable resources for studies on toxin genes. Its venom glands are located in the pectoral and dorsal fins. Although with such interesting biologic traits and great value in economy, Chinese yellow catfish is still lacking a sequenced genome. Here, we report a high-quality genome assembly of Chinese yellow catfish using a combination of next-generation Illumina and third-generation PacBio sequencing platforms. The final assembly reached 714 Mb, with a contig N50 of 970 kb and a scaffold N50 of 3.65 Mb, respectively. We also annotated 21,562 protein-coding genes, in which 97.59% were assigned at least one functional annotation. Based on the genome sequence, we analyzed toxin genes in Chinese yellow catfish. Finally, we identified 207 toxin genes and classified them into three major groups. Interestingly, we also expanded a previously reported sex-related region (to ≈6 Mb) in the achieved genome assembly, and localized two important toxin genes within this region. In summary, we assembled a high-quality genome of Chinese yellow catfish and performed high-throughput identification of toxin genes from a genomic view. Therefore, the limited number of toxin sequences in public databases will be remarkably improved once we integrate multi-omics data from more and more sequenced species.

Molecular characterization and expression analysis of tumor necrosis factor receptor-associated factors 3 and 6 in large yellow croaker (Larimichthys crocea)

The large yellow croaker (Larimichthys crocea) has a well-developed innate immune system. To gain a better understanding of the defense mechanisms involved in this system, we studied tumor necrosis factor receptor-associated factors (TRAFs), which play important roles in the Toll-like receptor (TLR) pathway. We characterized the full-length open reading frames and protein structures of TRAF3 and TRAF6 to determine their identities, and conducted phylogenetic analysis to determine their evolutionary relationships. To assess the roles of TRAFs in innate immune responses in the large yellow croaker, we performed quantitative reverse-transcription PCR (qRT-PCR) to characterize expression profiles in a range of tissues at different stages after challenge with polyinosinic polycytidylic acid (poly I:C) and Vibrio anguillarum. Following poly I:C challenge, the expression levels of TRAF3 and TRAF6 were highest in the kidneys and lowest in the spleen, whereas after infection with V. anguillarum, TRAF6 expression was the highest in the kidneys and lowest in the liver.

MiR-205 Mediated Cu-Induced Lipid Accumulation in Yellow Catfish Pelteobagrus fulvidraco

The present working hypothesis is that the Cu-induced changes in lipid metabolism may be mediated by miRNAs. Here, we describe the miRNA profile of the liver tissues of yellow catfish exposed to waterborne Cu, based on larger-scale sequencing of small RNA libraries. We identified a total of 172 distinct miRNAs. Among these miRNAs, compared to the control, mRNA expression levels of 16 miRNAs (miR-203a, 205, 1788-3p, 375, 31, 196a, 203b-3p, 2187-5p, 196d, 459-3p, 153a and miR-725, and two novel-miRNAs: chr4-1432, chr-7684) were down-regulated, and mRNA levels of miR-212 and chr20-5274 were up-regulated in Cu-exposed group. The functions of their target genes mainly involved ether lipid metabolism, glycerophospholipid metabolism, linoleic acid metabolism and α-linolenic acid metabolism. Cu exposure inhibited the expression of miR-205, whose predicted target genes were enriched in the pathway of lipid metabolism, including fas, lxrα, ddit3, lamp2, casp3a and baxa. These potential target genes were further verified by Dual-luciferase reporter gene assay. Using primary hepatocytes of yellow catfish, Cu incubation down-regulated miR-205 expression, and increased TG contents and FAS activity. LXR antagonist effectively ameliorate the Cu-induced change of TG content and FAS activity. These data suggest that down-regulation of the miRNA-205 may be an important step in Cu-induced changes in lipid metabolism in yellow catfish.

MiR-202 controls female fecundity by regulating medaka oogenesis

Female gamete production relies on coordinated molecular and cellular processes that occur in the ovary throughout oogenesis. In fish, as in other vertebrates, these processes have been extensively studied both in terms of endocrine/paracrine regulation and protein expression and activity. The role of small non-coding RNAs in the regulation of animal reproduction remains however largely unknown and poorly investigated, despite a growing interest for the importance of miRNAs in a wide variety of biological processes. Here, we analyzed the role of miR-202, a miRNA predominantly expressed in male and female gonads in several vertebrate species. We studied its expression in the medaka ovary and generated a mutant line (using CRISPR/Cas9 genome editing) to determine its importance for reproductive success with special interest for egg production. Our results show that miR-202-5p is the most abundant mature form of the miRNA and that it is expressed in granulosa cells and in the unfertilized egg. The knock out (KO) of mir-202 gene resulted in a strong phenotype both in terms of number and quality of eggs produced. Mutant females exhibited either no egg production or produced a dramatically reduced number of eggs that could not be fertilized, ultimately leading to no reproductive success. We quantified the size distribution of the oocytes in the ovary of KO females and performed a large-scale transcriptomic analysis approach to identified dysregulated molecular pathways. Together, cellular and molecular analyses indicate that the lack of miR-202 impairs the early steps of oogenesis/folliculogenesis and decreases the number of large (i.e. vitellogenic) follicles, ultimately leading to dramatically reduced female fecundity. This study sheds new light on the regulatory mechanisms that control the early steps of follicular development, including possible targets of miR-202-5p, and provides the first in vivo functional evidence that a gonad-predominant microRNA may have a major role in female reproduction.

The role of small non-coding RNAs in the regulation of animal reproduction remains poorly investigated, despite a growing interest for the importance of miRNAs in a wide variety of biological processes. Here, we analyzed the role of miR-202, a miRNA predominantly expressed in gonads in vertebrate. We studied its expression in the medaka ovary and knocked out the mir-202 gene to study its importance for reproductive success. We showed that the lack of miR-202 results in the sterility of both females and males. In particular, it led to a drastic reduction of both the number and the quality of eggs produced by females. Mutant females exhibited either no egg production or produced a drastically reduced number of eggs that could not be fertilized, ultimately leading to no reproductive success. Quantitative histological and molecular analyses indicated that mir-202 KO impairs oocyte development and is also associated with the dysregulation of many genes that are critical for reproduction. This study sheds new light on the regulatory mechanisms that control oogenesis, including possible targets of miR-202-5p, and provides the first in vivo functional evidence that a gonad-predominant microRNA may have a major role in female reproduction.

A myeloid differentiation factor 88 gene from yellow catfish Pelteobagrus fulvidraco and its molecular characterization in response to polyriboinosinic polyribocytidylic acid and lipopolysaccharide challenge

Myeloid differentiation factor 88 (MyD88) is an adaptor protein of Toll-like receptor (TLR) signalling pathways that activates the innate immune system. Herein, MyD88 was identified in the economically important freshwater fish Pelteobagrus fulvidraco. The complete 2156 bp PfMyD88 cDNA includes a 147 bp 5'-untranslated region (UTR), a 1133 bp 3'-UTR, and an open reading frame (ORF) of 876 bp encoding a 291 residue protein containing Death and Toll/interleukin-1 receptor (TIR) domains. The deduced protein sequence shares 88.8%, 73.8% and 59.3% identity with orthologs in Ictalurus punctatus, Danio rerio and Homo sapiens, respectively. qRT-PCR revealed expression in all tested tissues, highest in trunk kidney, followed by spleen, and lowest in muscle. After challenge with lipopolysaccharide (LPS) or polyriboinosinic polyribocytidylic acid (Poly I:C), PfMyD88 expression was up-regulated in blood, liver, head kidney and spleen. Thus, PfMyD88 acts in innate immunity in P. fulvidraco.

Discovery and functional characterization of microRNAs and their potential roles for gonadal development in spotted knifejaw, Oplegnathus punctatus

The spotted knifejaw (Oplegnathus punctatus) is a newly emerging economical fishery species in China. Studies focused on the regulation of gonadal development and gametogenesis of spotted knifejaw are still insufficient. As a key post-transcriptional regulator, miRNAs have been shown to play important roles in development and reproduction systems. In this study, small RNA deep sequencing in ovary and testis of spotted knifejaw were performed to screen miRNA expression patterns. After sequencing and bioinformatics analysis, a total of 247 conserved known miRNAs and 41 novel miRNAs were identified in spotted knifejaw gonads for the first time. In addition, 36 miRNAs were differentially expressed between testis and ovary. The putative target genes of differentially expressed (DE) miRNAs were significantly enriched in several pathways related to sexual differentiation and gonadal development, such as steroid hormone biosynthesis. Sequencing data was validated through qRT-PCR analysis of selected DE miRNAs. Dual-luciferase reporter analyses of filtered miRNA-target gene pairs confirmed that opu-miR-27b-3p targeted in piwi2 and mov10l1 3' UTRs and down-regulated their expressions in spotted knifejaw. The notion that mov10l1 and piwi2 enhance germ cells proliferation and regulate gonadal development and gametogenesis suggests that opu-miR-27b-3p may attenuated this process in the gonads of spotted knifejaw. These findings provided insights into regulatory roles of gonadal miRNAs and supplied fundamental resources for further studies on miRNA-mediated post-transcriptional regulation in reproductive system of spotted knifejaw.

An miR-200 Cluster on Chromosome 23 Regulates Sperm Motility in Zebrafish

Besides its well-documented roles in cell proliferation, apoptosis, and carcinogenesis, the function of the p53-microRNA axis has been recently revealed in the reproductive system. Recent studies indicated that miR-200 family members are dysregulated in nonobstructive azoospermia patients, whereas their functions remain poorly documented. The aim of this study was to investigate the function of the miR-200 family on zebrafish testis development and sperm activity. There was no substantial difference in testis morphology and histology between wild-type (WT) and knockout zebrafish with deletion of miR-200 cluster on chromosome 6 (chr6-miR-200-KO) or on chromosome 23 (chr23-miR-200-KO). Interestingly, compared with WT zebrafish, the chr6-miR-200-KO zebrafish had no difference on sperm motility, whereas chr23-miR-200-KO zebrafish showed significantly improved sperm motility. Consistently, ectopic expression of miR-429a, miR-200a, and miR-200b, which are located in the miR-200 cluster on chromosome 23, significantly reduced motility traits of sperm. Several sperm motility-related genes, such as amh, wt1a, and srd5a2b have been confirmed as direct targets of miR-200s on chr23. 17α-ethynylestradiol (EE2) exposure resulted in upregulated expression of p53 and miR-429a in testis and impairment of sperm motility. Strikingly, in p53 mutant zebrafish testis, the expression levels of miR-200s on chr23 were significantly reduced and accompanied by a stimulation of sperm motility. Moreover, the upregulation of miR-429a associated with EE2 treatment was abolished in testis with p53 mutation. And the impairment of sperm activity by EE2 treatment was also eliminated when p53 was mutated. Together, our results reveal that miR-200 cluster on chromosome 23 controls sperm motility in a p53-dependent manner.

miRNAome expression profiles in the gonads of adult Melopsittacus undulatus

The budgerigar (Melopsittacus undulatus) is one of the most widely studied parrot species, serving as an excellent animal model for behavior and neuroscience research. Until recently, it was unknown how sexual differences in the behavior, physiology, and development of organisms are regulated by differential gene expression. MicroRNAs (miRNAs) are endogenous short non-coding RNA molecules that can post-transcriptionally regulate gene expression and play a critical role in gonadal differentiation as well as early development of animals. However, very little is known about the role gonadal miRNAs play in the early development of birds. Research on the sex-biased expression of miRNAs in avian gonads are limited, and little is known about M. undulatus. In the current study, we sequenced two small non-coding RNA libraries made from the gonads of adult male and female budgerigars using Illumina paired-end sequencing technology. We obtained 254 known and 141 novel miRNAs, and randomly validated five miRNAs. Of these, three miRNAs were differentially expressed miRNAs and 18 miRNAs involved in sexual differentiation as determined by functional analysis with GO annotation and KEGG pathway analysis. In conclusion, this work is the first report of sex-biased miRNAs expression in the budgerigar, and provides additional sequences to the avian miRNAome database which will foster further functional genomic research.

Identification and expression analysis of microRNAs in medaka gonads

Gonad development is a highly regulated, coordinated biological process and increasing evidences have indicated that microRNA (miRNA) may be involved in this dynamic program. Medaka (Oryzias latipes) is a good model for reproductive research as it has distinct sex determining genes, however, research in gonadal miRNAs is lacked. In this study, two small RNA libraries from the ovaries and testes were constructed and sequenced. A total of 285 conserved and 388 novel miRNAs were obtained, among which 142 mature miRNAs were significantly (> two-fold change) up or down regulated in the testis compared to the ovary. Quantitative Reverse Transcription-Polymerase Chain Reaction (qRT-PCR) analysis showed that miR-430c, miR-26a and miR-202-5p were expressed in a gonad-specific or sex-biased pattern. Fluorescence in situ hybridization (FISH) indicated that miR-202-5p was present throughout spermatogenesis and was only detected at the early stages of oogenesis, this sex biased expression pattern suggested that miR-202-5p might be a crucial candidate in male differentiation and development. Our study provides the repertoire, a comprehensive annotation of miRNAs from gonads and a reference for functional studies of miRNAs in medaka.

Leucine mediates autophagosome-lysosome fusion and improves sperm motility by activating the PI3K/Akt pathway

Amino acid supplementation is an efficient and effective strategy to increase sperm quality. In our research, a comparative study was conducted to screen free amino acids to improve sperm motility, and we found that leucine was the most efficient one. Leucine treatment increases sperm motility depending on the activation of PI3K/Akt signaling pathway, while the chemical inhibitor of PI3K/Akt signal could reduce the amount of pAkt activated by leucine treatment. Moreover, leucine treatment improved the expression of P62 and LC3-II, substantially suppressed the autophagy process in zebrafish testis. In vitro studies showed that leucine could reduce the fusion of autophagosome and lysosome that was indicated by the co-localization of EGFP-LC3 and lysosome marker. Two chemical modulators of autophagy, such as LY294002 (the inhibitor of PI3K/Akt signal) and chloroquine were administered to investigate the process of autophagy on zebrafish sperm motility. LY294002 inhibited autophagosome formation to reduced sperm motility, while chloroquine inhibited the fusion of autophagosome and lysosome to improve sperm motility. Our data suggest that short-term treatment with leucine could increase zebrafish sperm motility by affecting the autophagy and inhibiting the fusion of autophagosome and lysosomes, depending on the activation of PI3K/Akt signaling pathway.

Integrative testis transcriptome analysis reveals differentially expressed miRNAs and their mRNA targets during early puberty in Atlantic salmon

BACKGROUND

Our understanding of the molecular mechanisms implementing pubertal maturation of the testis in vertebrates is incomplete. This topic is relevant in Atlantic salmon aquaculture, since precocious male puberty negatively impacts animal welfare and growth. We hypothesize that certain miRNAs modulate mRNAs relevant for the initiation of puberty. To explore which miRNAs regulate mRNAs during initiation of puberty in salmon, we performed an integrated transcriptome analysis (miRNA and mRNA-seq) of salmon testis at three stages of development: an immature, long-term quiescent stage, a prepubertal stage just before, and a pubertal stage just after the onset of single cell proliferation activity in the testis.

RESULTS

Differentially expressed miRNAs clustered into 5 distinct expression profiles related to the immature, prepubertal and pubertal salmon testis. Potential mRNA targets of these miRNAs were predicted with miRmap and filtered for mRNAs displaying negatively correlated expression patterns. In summary, this analysis revealed miRNAs previously known to be regulated in immature vertebrate testis (miR-101, miR-137, miR-92b, miR-18a, miR-20a), but also miRNAs first reported here as regulated in the testis (miR-new289, miR-30c, miR-724, miR-26b, miR-new271, miR-217, miR-216a, miR-135a, miR-new194 and the novel predicted n268). By KEGG enrichment analysis, progesterone signaling and cell cycle pathway genes were found regulated by these differentially expressed miRNAs. During the transition into puberty we found differential expression of miRNAs previously associated (let7a/b/c), or newly associated (miR-15c, miR-2184, miR-145 and the novel predicted n7a and b) with this stage. KEGG enrichment analysis revealed that mRNAs of the Wnt, Hedgehog and Apelin signaling pathways were potential regulated targets during the transition into puberty. Likewise, several regulated miRNAs in the pubertal stage had earlier been associated (miR-20a, miR-25, miR-181a, miR-202, let7c/d/a, miR-125b, miR-222a/b, miR-190a) or have now been found connected (miR-2188, miR-144, miR-731, miR-8157 and the novel n2) to the initiation of puberty.

CONCLUSIONS

This study has - for the first time - linked testis maturation to specific miRNAs and their inversely correlated expressed targets in Atlantic salmon. The study indicates a broad functional conservation of already known miRNAs and associated pathways involved in the transition into puberty in vertebrates. The analysis also reveals miRNAs not previously associated with testis tissue or its maturation, which calls for further functional studies in the testis.

Transcriptome analysis of yellow catfish (Pelteobagrus fulvidraco) liver challenged with polyriboinosinic polyribocytidylic acid (poly I:C)

Yellow catfish (Pelteobagrus fulvidraco) is one of the most important economic freshwater species in China. However, infection by bacterial pathogenic diseases has caused high mortality and great economic loss in aquaculture. It is necessary for disease control to know more about the P. fulvidraco immune system and its related genes in response to bacterial or viral infections. In this study, the transcriptomic profiles of liver from P. fulvidraco stimulated by polyriboinosinic polyribocytidylic acid (poly I:C) was analyzed using high-throughput sequencing method. After assembly and annotation, total 67,447 unigenes were acquired, with an average length of 1091 bp. Under the infection of poly I:C, 522 differentially expressed genes (DEGs) were identified, including 307 up-regulated genes and 215 down-regulated genes. To further investigate the immune-related DEGs, Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment were performed. The result of GO enrichment indicated gene response to external stimulus, regulation of response to stimulus, cellular response to stimulus, immune response and immune system progress. Significant KEGG enrichment analysis identified major immune related pathways. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis revealed that 13 immune response genes were identified to be up-regulated after 12 h of poly I:C stimulation compared to controls. Taken together, the results of our study are beneficial for better understanding of the immune system and defense mechanisms of yellow catfish in response to poly I:C infection.

Integrated analysis of mRNA-seq and miRNA-seq reveals the potential roles of sex-biased miRNA-mRNA pairs in gonad tissue of dark sleeper (Odontobutis potamophila)

Background

The dark sleeper (Odontobutis potamophila) is an important commercial fish species which shows a sexually dimorphic growth pattern. However, the lack of sex transcriptomic data is hindering further research and genetically selective breeding of the dark sleeper. In this study, integrated analysis of mRNA and miRNA was performed on gonad tissue to elucidate the molecular mechanisms of sex determination and differentiation in the dark sleeper.

Results

A total of 143 differentially expressed miRNAs and 16,540 differentially expressed genes were identified. Of these, 8103 mRNAs and 75 miRNAs were upregulated in testes, and 8437 mRNAs and 68 miRNAs were upregulated in ovaries. Integrated analysis of miRNA and mRNA expression profiles predicted more than 50,000 miRNA-mRNA interaction sites, and among them 27,583 negative miRNA-mRNA interactions. A number of sex related genes were targeted by sex-biased miRNAs. The relationship between 15 sex-biased genes and 15 sex-biased miRNAs verified by using qRT-PCR were described. Additionally, a number of SNPs were revealed through the transcriptome data.

Conclusions

The overall results of this study facilitate our understanding of the molecular mechanism underlying sex determination and differentiation and provide valuable genomic information for selective breeding of the dark sleeper.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3995-9) contains supplementary material, which is available to authorized users.

The Functions of MicroRNA-200 Family in Ovarian Cancer: Beyond Epithelial-Mesenchymal Transition

The majority of studies on microRNA-200 family members (miR-200s) in human cancers are based on the premise that miR-200s maintain epithelial cell integrity by suppressing epithelial-mesenchymal transition (EMT) through direct inhibition of mesenchymal transcription factors zinc finger E-box-binding homeobox 1/2 (ZEB1/ZEB2) and transforming growth factor-β (TGF-β), a potent inducer of EMT. Hence, downregulation of miR-200 in cancer cells promotes EMT and cancer metastasis. Yet, miR-200s are highly expressed in ovarian cancer, and ovarian cancer metastasizes primarily by dissemination within the pelvic cavity. In this review, we will refocus the epithelial property of ovarian cancer cells and the role of miR-200s in safeguarding this property, as well as the diverse roles of miR-200s in inclusion cyst formation, cancer cell growth, collective movement, angiogenesis, exosome-mediated cell communication, and chemoresponse. Taken together, miR-200s play a significant role in the initiation, progression and metastasis of ovarian cancer and may serve as diagnostic biomarkers and a target in therapeutic development.

The Functions of MicroRNA-200 Family in Ovarian Cancer: Beyond Epithelial-Mesenchymal Transition

The majority of studies on microRNA-200 family members (miR-200s) in human cancers are based on the premise that miR-200s maintain epithelial cell integrity by suppressing epithelial-mesenchymal transition (EMT) through direct inhibition of mesenchymal transcription factors zinc finger E-box-binding homeobox 1/2 (ZEB1/ZEB2) and transforming growth factor-β (TGF-β), a potent inducer of EMT. Hence, downregulation of miR-200 in cancer cells promotes EMT and cancer metastasis. Yet, miR-200s are highly expressed in ovarian cancer, and ovarian cancer metastasizes primarily by dissemination within the pelvic cavity. In this review, we will refocus the epithelial property of ovarian cancer cells and the role of miR-200s in safeguarding this property, as well as the diverse roles of miR-200s in inclusion cyst formation, cancer cell growth, collective movement, angiogenesis, exosome-mediated cell communication, and chemoresponse. Taken together, miR-200s play a significant role in the initiation, progression and metastasis of ovarian cancer and may serve as diagnostic biomarkers and a target in therapeutic development.

The Functions of MicroRNA-200 Family in Ovarian Cancer: Beyond Epithelial-Mesenchymal Transition

The majority of studies on microRNA-200 family members (miR-200s) in human cancers are based on the premise that miR-200s maintain epithelial cell integrity by suppressing epithelial-mesenchymal transition (EMT) through direct inhibition of mesenchymal transcription factors zinc finger E-box-binding homeobox 1/2 (ZEB1/ZEB2) and transforming growth factor-β (TGF-β), a potent inducer of EMT. Hence, downregulation of miR-200 in cancer cells promotes EMT and cancer metastasis. Yet, miR-200s are highly expressed in ovarian cancer, and ovarian cancer metastasizes primarily by dissemination within the pelvic cavity. In this review, we will refocus the epithelial property of ovarian cancer cells and the role of miR-200s in safeguarding this property, as well as the diverse roles of miR-200s in inclusion cyst formation, cancer cell growth, collective movement, angiogenesis, exosome-mediated cell communication, and chemoresponse. Taken together, miR-200s play a significant role in the initiation, progression and metastasis of ovarian cancer and may serve as diagnostic biomarkers and a target in therapeutic development.