Dynamics of miRNA transcriptome during gonadal development of zebrafish

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.

The intricacies of isomiRs: from classification to clinical relevance

MicroRNAs (miRNAs) and isoforms of their archetype, called isomiRs, regulate gene expression via complementary base-pair binding to messenger RNAs (mRNAs). The partially evolutionarily conserved isomiR sequence variations are differentially expressed among tissues, populations, and genders, and between healthy and diseased states. Aiming towards the clinical use of isomiRs as diagnostic biomarkers and for therapeutic purposes, several challenges need to be addressed, including (i) clarification of isomiR definition, (ii) improved annotation in databases with new standardization (such as the mirGFF3 format), and (iii) improved methods of isomiR detection, functional verification, and in silico analysis. In this review we discuss the respective challenges, and highlight the opportunities for clinical use of isomiRs, especially in the light of increasing amounts of next-generation sequencing (NGS) data.

The lncRNA lnc_AABR07044470.1 promotes the mitochondrial-damaged inflammatory response to neuronal injury via miR-214-3p/PERM1 axis in acute ischemic stroke

PURPOSE

We investigated the role of lnc_AABR07044470.1 on the occurrence and development of acute ischemic stroke (AIS) and neuronal injury by targeting the miR-214-3p/PERM1 axis to find a novel clinical drug target and prediction and treatment of AIS.

METHODS

The mouse AIS animal model was used in vivo experiments and hypoxia/reoxygenation cell model in vitro was established. Firstly, infarction volume and pathological changes of mouse hippocampal neurons were detected using HE staining. Secondly, rat primary neuron apoptosis was detected by flow cytometry assay. The numbers of neuron, microglia and astrocytes were detected using immunofluorescence (IF). Furthermore, binding detection was performed by bioinformatics database and double luciferase reporter assay. Lnc_AABR07044470.1 localization was performed using fluorescence in situ hybridization (FISH).Lnc_AABR07044470.1, miR-214-3pand PERM1mRNA expression was performed using RT-qPCR. NLRP3, ASC, Caspase-1 and PERM1 protein expression was performed using Western blotting. IL-1β was detected by ELISA assay.

RESULTS

Mouse four-vessel occlusion could easily establish the animal model, and AIS animal model had an obvious time-dependence. HE staining showed that, compared with the sham group, infarction volume and pathological changes of mouse hippocampal neurons were deteriorated in the model group. Furthermore, compared with the sham group, neurons were significantly reduced, while microglia and astrocytes were significantly activated. Moreover, the bioinformatics prediction and detection of double luciferase reporter confirmed the binding site of lnc_AABR07044470.1 to miR-214-3p and miR-214-3p to Perm1. lnc_AABR07044470.1 and PERM1 expression was significantly down-regulated and miR-214-3pexpression was significantly up-regulated in AIS animal model in vivo. At the same time, the expression of inflammasome NLRP3, ASC, Caspase-1 and pro-inflammatory factor IL-1β was significantly up-regulated in vivo and in vitro. The over-expression of lnc_AABR07044470.1 and miR-214-3p inhibitor could inhibit the neuron apoptosis and the expression of inflammasome NLRP3, ASC, Caspase-1 and pro-inflammatory factor IL-1β and up-regulate the expression of PERM1 in vitro. Finally, over-expression of lnc_AABR07044470.1 and miR-214-3p inhibitor transfected cell model was significant in relieving the AIS and neuronal injury.

CONCLUSION

Lnc_AABR07044470.1 promotes inflammatory response to neuronal injury via miR-214-3p/PERM1 axis in AIS.

Human Umbilical Cord Mesenchymal Stem Cell-Derived Exosomes Rescue Testicular Aging

BACKGROUND

Testicular aging is associated with diminished fertility and certain age-related ailments, and effective therapeutic interventions remain elusive. Here, we probed the therapeutic efficacy of exosomes derived from human umbilical cord mesenchymal stem cells (hUMSC-Exos) in counteracting testicular aging.

METHODS

We employed a model of 22-month-old mice and administered intratesticular injections of hUMSC-Exos. Comprehensive analyses encompassing immunohistological, transcriptomic, and physiological assessments were conducted to evaluate the effects on testicular aging. Concurrently, we monitored alterations in macrophage polarization and the oxidative stress landscape within the testes. Finally, we performed bioinformatic analysis for miRNAs in hUMSC-Exos.

RESULTS

Our data reveal that hUMSC-Exos administration leads to a marked reduction in aging-associated markers and cellular apoptosis while promoting cellular proliferation in aged testis. Importantly, hUMSC-Exos facilitated the restoration of spermatogenesis and elevated testosterone synthesis in aged mice. Furthermore, hUMSC-Exos could attenuate inflammation by driving the phenotypic shift of macrophages from M1 to M2 and suppress oxidative stress by reduced ROS production. Mechanistically, these efficacies against testicular aging may be mediated by hUMSC-Exos miRNAs.

CONCLUSIONS

Our findings suggest that hUMSC-Exos therapy presents a viable strategy to ameliorate testicular aging, underscoring its potential therapeutic significance in managing testicular aging.

The Discovery and Characterization of Conserved and Novel miRNAs in the Different Developmental Stages and Organs of Pikeperch (Sander lucioperca)

Micro RNAs (miRNAs) are short non-coding RNAs that act as post-transcriptional gene expression regulators. Genes regulated in vertebrates include those affecting growth and development or stress and immune response. Pikeperch (Sander lucioperca) is a species that is increasingly being considered for farming in recirculation aquaculture systems. We characterized the pikeperch miRNA repertoire to increase the knowledge of the genomic mechanisms affecting performance and health traits by applying small RNA sequencing to different developmental stages and organs. There were 234 conserved and 8 novel miRNA genes belonging to 104 families. A total of 375 unique mature miRNAs were processed from these genes. Many mature miRNAs showed high relative abundances or were significantly more expressed at early developmental stages, like the miR-10 and miR-430 family, let-7, the miRNA clusters 106-25-93, and 17-19-92. Several miRNAs associated with immune responses (e.g., slu-mir-731-5p, slu-mir-2188-5p, and slu-mir-8159-5p) were enriched in the spleen. The mature miRNAs slu-mir-203a-3p and slu-mir-205-5p were enriched in gills. These miRNAs are similarly abundant in many vertebrates, indicating that they have shared regulatory functions. There was also a significantly increased expression of the disease-associated miR-462/miR-731 cluster in response to hypoxia stress. This first pikeperch miRNAome reference resource paves the way for future functional studies to identify miRNA-associated variations that can be utilized in marker-assisted breeding programs.

MiR-202-5p Regulates Geese Follicular Selection by Targeting BTBD10 to Regulate Granulosa Cell Proliferation and Apoptosis

The regulation of granulosa cells (GCs) proliferation and apoptosis is the key step in follicular selection which determines the egg production performance of poultry. miR-202-5p has been reported to be involved in regulating the proliferation and apoptosis of mammalian ovarian GCs. However, its role in regulating the proliferation and apoptosis of goose GCs is still unknown. In the present study, the GCs of pre-hierarchical follicles (phGCs, 8-10 mm) and those of hierarchical follicles (hGCs, F2-F4) were used to investigate the role of miR-202-5p in cell proliferation and apoptosis during follicle selection. In phGCs and hGCs cultured in vitro, miR-202-5p was found to negatively regulate cell proliferation and positively regulate cell apoptosis. The results of RNA-seq showed that BTB Domain Containing 10 (BTBD10) is predicted to be a key target gene for miR-202-5p to regulate the proliferation and apoptosis of GCs. Furthermore, it is confirmed that miR-202-5p can inhibit BTBD10 expression by targeting its 3'UTR region, and BTBD10 was revealed to promote the proliferation and inhibit the apoptosis of phGCs and hGCs. Additionally, co-transfection with BTBD10 effectively prevented miR-202-5p mimic-induced cell apoptosis and the inhibition of cell proliferation. Meanwhile, miR-202-5p also remarkably inhibited the expression of Phosphatidylinositol-4,5-Bisphosphate 3-Kinase Catalytic Subunit Beta (PIK3CB) and AKT Serine/Threonine Kinase 1 (AKT1), while it was significantly restored by BTBD10. Overall, miR-202-5p suppresses the proliferation and promotes the apoptosis of GCs through the downregulation of PIK3CB/AKT1 signaling by targeting BTBD10 during follicular selection. Our study provides a theoretical reference for understanding the molecular mechanism of goose follicular selection, as well as a candidate gene for molecular marker-assisted breeding to improve the geese' egg production performance.

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.

Deciphering sex-specific miRNAs as heat-recorders in zebrafish

In the last decade, a plethora of microRNAs (miRNAs) has been reported in a wide variety of physiological processes, including reproduction, in many aquatic organisms. However, miRNAome alterations occurred by environmental cues due to water temperature increment have not yet been elucidated. With the aim to identify epigenetic regulations mediated by miRNAs in the gonads in a climate change scenario, the animal model zebrafish (Danio rerio) were subjected to high temperatures during sex differentiation, a treatment that results in male-skewed sex ratios in the adulthood. Once the fish reached adulthood, gonads were sequenced by high-throughput technologies and a total of 23 and 1 differentially expressed miRNAs in ovaries and testes, respectively, were identified two months after the heat treatment. Most of these heat-recorder miRNAs were involved in human sex-related cancer and about 400 predicted-target genes were obtained, some with reproduction-related functions. Their synteny in the zebrafish genome was, for more than half of the predicted target genes, in the chromosomes 7, 2, 4, 3 and 11 in the ovaries, chromosome 4 being the place where the sex-associated-region (sar) is localized in wild zebrafish. Further, spatial localization in the gonads of two selected heat-recorder miRNAs (miR-122-5p and miR-146-5p) showed exclusive expression in the ovarian germ cells. The present study expands the catalog of sex-specific miRNAs and deciphers, for the first time, thermosensitive miRNAs in the zebrafish gonads that might be used as potential epimarkers to predict environmental past events.

Integrated mRNA and miRNA Expression Profile Analysis of Female and Male Gonads in Acrossocheilus fasciatus

MicroRNAs (miRNAs) are regarded as key regulators in gonadal development and sex determination in diverse organisms. However, the functions of miRNAs in gonads of Acrossocheilus fasciatus, an economically important freshwater species in the south of China, are still unclear. Here, high-throughput sequencing was performed to investigate the mRNA and miRNAs on gonads of A. fasciatus. In total, 49,447 unigenes were obtained, including 11,635 differentially expressed genes (DEGs), among which 4147 upregulated genes and 7488 downregulated genes in the testis compared to the ovary, while 300 (237 known, and 63 novel) miRNAs with 36 differentially expressed miRNAs (DEMs) were identified, from which 17 upregulated and 19 downregulated DEMs. GO and KEGG enrichment analysis were performed to analyze the potential biological functions of DEGs and DEMs. Using qRT-PCR, 9 sex-related genes and 9 miRNAs were selected to verify the sequencing data. By dual-luciferase reporter assay, miR-22a-5p and miR-22b-5p interaction with piwil1, and miR-10d-5p interaction with piwil2 were identified. These findings could provide a reference for miRNA-regulated sex control of A. fasciatus and may reveal new insights into aquaculture and breeding concepts.

MicroRNA expression in male infertility

Male infertility is a multifactorial disorder that involves different physiopathological mechanisms and multiple genes. In this sense, we analyse the role of miRNAs in this pathology. Gene expression analysis can provide relevant information to detect biomarkers, signalling pathways, pathologic mechanisms, and potential therapeutic targets for the disease. In this review, we describe four miRNA microarrays related to patients who present infertility diseases, including azoospermia, asthenozoospermia, and oligoasthenozoospermic. We selected 13 miRNAs with altered expressions in testis tissue (hsa-miR-122-5p, hsa-miR-145-5p, hsa-miR-16-5p, hsa-miR-193a-3p, hsa-miR-19a-3p, hsa-miR-23a-3p, hsa-miR-30b-5p, hsa-miR-34b-5p, hsa-miR-34c-5p, hsa-miR-374b-5p, hsa-miR-449a, hsa-miR-574-3p and hsa-miR-92a-3p), and systematically examine the mechanisms of four relevant miRNAs (hsa-miR-16-5p, hsa-miR-19a-3p, hsa-miR-92a-3p and hsa-miR-30b-5p) which we found that regulated a large number of proteins. An interaction network was generated, and its connections allowed us to identify signalling pathways and interactions between proteins associated with male infertility. In this way, we confirm that the most affected and relevant pathway is the PI3K-Akt signalling.

Testicular miRNAs and tsRNAs provide insight into gene regulation during overwintering and reproduction of Onychostoma macrolepis

The late overwintering period and breeding period are two important developmental stages of testis in Onychostoma macrolepis. Small non-coding RNAs (sncRNAs) are well-known regulators of biological processes associated with numerous biological processes. This study aimed to elucidate the roles of four sncRNA classes (microRNAs [miRNAs], Piwi-interacting RNAs [piRNAs], tRNA-derived small RNAs [tsRNAs], and rRNA-derived small RNAs [rsRNAs]) across testes in the late overwintering period (in March) and breeding period (in June) by high-throughput sequencing. The testis of O. macrolepis displayed the highest levels of piRNAs and lowest levels of rsRNAs. Compared with miRNAs and tsRNAs in June, tsRNAs in March had a higher abundance, while miRNAs in March had a much lower abundance. Bioinformatics analysis identified 1,362 and 1,340 differentially expressed miRNAs and tsRNAs, respectively. Further analysis showed that miR-200-1, miR-143-1, tRFi-Lys-CTT-1, and tRFi-Glu-CTC-1 could play critical roles during the overwintering and breeding periods. Our findings provided an unprecedented insight to reveal the epigenetic mechanism underlying the overwintering and reproduction process of male O. macrolepis.

A comprehensive overview of ovarian small non-coding RNAs in the late overwintering and breeding periods of Onychostoma macrolepis

The development of the ovary of Onychostoma macrolepis undergoes distinct annual cyclic changes in which small non-coding RNAs (sncRNAs) could play vital roles. In this study, four sncRNA classes in the ovary, including miRNA, piRNAs, tsRNA, and rsRNAs, were systematically profiled by high-throughput sequencing. In adult ovaries of O. macrolepis, 247 miRNAs and 235 tsRNAs were identified as differentially expressing in the late overwintering period (in March) and breeding period (in June). Some up-regulated sncRNAs in March, such as miR-125-1 and tRFi-Lys-CTT-1, could be involved in inhibiting biomolecule metabolism and enhancing stress tolerance during the overwintering period. Compared with the level expression of sncRNAs in March, some sncRNAs were up-regulated in June, such as miR-146-1 and tRFi-Gly-GCC-1, and could be involved in influencing molecular synthesis and metabolism, enhancing oocyte proliferation and maturation, accelerating ovarian development, and increasing fertilization of oocytes by regulating related target mRNAs. The results suggested that sncRNAs in the ovary of Onychostoma macrolepis not only reflect characteristics of the fish's physiology at different developmental periods, but also directly affect ovarian development and oocyte maturation during the breeding period. In conclusion, these results significantly advance our understanding of the roles of sncRNA during overwintering and reproduction periods, and provide a novel perspective for uncovering characteristics of the special overwintering ecology and reproductive physiology of an atypical cavefish.

Localisation and in silico based functional analysis of miR-202 in bull testis

Bull fertility is pivotal to the prosperity of the cattle industry worldwide. miR‐202 has been shown to be gonad specific and to have key roles in gonad function in different species. To further understand the involvement of miR‐202 in bull reproduction, this study aimed to establish its localization in bovine testicular tissue and to identify putative biological functions using bioinformatics approaches. We assessed the miR‐202 expression in paraffin‐embedded tissue samples collected form an abattoir using in situ hybridization. miR‐202 was present in Sertoli cells and in germ cells at different stages of development. Using available databases, a total of 466 predicted gene targets of miR‐202 were identified. Functional annotation revealed that miR‐202 target genes were mainly associated with protein modification and phosphorylation processes as well as longevity regulating pathway. Moreover, genes in the longevity regulating pathway mapped to PI3K/Akt/mTOR pathway which is involved in promoting proliferation of testicular cells and spermatogenesis. These findings suggest that miR‐202 plays important roles in regulating proliferation and viability of testicular cells including somatic and germ cells.

Effects of salinity acclimation on histological characteristics and miRNA expression profiles of scales in juvenile rainbow trout (Oncorhynchus mykiss)

Background

The scales serve as an ideal model for studying the regulatory mechanism of bone homeostasis in fish. To explore the effect of salinity acclimation on bone metabolism of juvenile rainbow trout (Oncorhynchus mykiss), three sampling time points during salinity acclimation (7D, 14D and 21D) were selected to detect variations in histological characteristics. In the histological analysis, osteoblast marker enzymes alkaline phosphatase (ALP), osteoclast marker tartrate-resistant acid phosphatase (TRAcP) and calcium salt deposit areas (Von Kossa′s) were detected. Changes in calcium (Ca), phosphorus (P) and the molar mass ratio of calcium to phosphorus (Ca/P) in the scales were also detected by Inductively Coupled Plasma Mass Spectrometry (ICP-MS). In addition, the global MicroRNA (miRNA) expression profiles during salinity acclimation were examined using Illumina sequencing platform because of their important regulatory roles in teleost biological processes.

Results

Twelve independent miRNA libraries were constructed, a total of 664 known and 92 putative novel miRNAs were identified. A total of 290 differentially expressed (DE) miRNAs were found in clusters with significant trends in the cluster analysis, and five types of clustering patterns were obtained; 22,374 DE predicted target genes of the aforementioned 290 DE miRNAs were obtained, 5957 of which clustered in six types of clustering patterns with a significant trend. To better understand the functions of the DE miRNAs, GO and KEGG analysis was performed on the 5957 target genes, as a result, they were significantly enriched in bone metabolism related signaling pathways such as MAPK signaling pathway, Calcium signaling pathway, Wnt signaling pathway, Mineral absorption and NF-kappa B signaling pathway. Six DE miRNAs were randomly selected and their expression were verified by quantitative real-time PCR (qRT-PCR), the expression trends were consistent with the results of transcriptome sequencing.

Conclusions

The DE miRNAs and DE target genes identified in this study might play an important role in regulation of bone metabolism during salinity acclimation, relative genes or pathways could serve as key candidates for further studies to elucidate molecular mechanism of teleost bone metabolism, and help performing salinity acclimation and developing marine culture of salmonid species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08531-7.

Sex-Biased miRNAs in the Gonads of Adult Chinese Alligator (Alligator sinensis) and Their Potential Roles in Sex Maintenance

MicroRNA (miRNA) is a category of single-stranded non-coding small RNA (sRNA) that regulates gene expression by targeting mRNA. It plays a key role in the temperature-dependent sex determination of Chinese alligator (Alligator sinensis), a reptile whose sex is determined solely by the temperature during the incubation period and remains stable thereafter. However, the potential function of miRNAs in the gonads of adult Chinese alligators is still unclear. Here, we prepared and sequenced sRNA libraries of adult female and male alligator gonads, from breeding (in summer) and hibernating (in winter) animals. We obtained 130 conserved miRNAs and 683 novel miRNAs, which were assessed for sex bias in summer and winter; a total of 65 miRNAs that maintained sex bias in both seasons were identified. A regulatory network of sex-biased miRNAs and genes was constructed. Sex-biased miRNAs targeted multiple genes in the meiosis pathway of adult Chinese alligator oocytes and the antagonistic gonadal function maintenance pathway, such as MOS, MYT1, DMRT1, and GDF9. Our study emphasizes the function of miRNA in the epigenetic mechanisms of sex maintenance in crocodilians.

A cross-species comparative approach to assessing multi- and transgenerational effects of endocrine disrupting chemicals

A wide range of chemicals have been identified as endocrine disrupting chemicals (EDCs) in vertebrate species. Most studies of EDCs have focused on exposure of both male and female adults to these chemicals; however, there is clear evidence that EDCs have dramatic effects when mature or developing gametes are exposed, and consequently are associated with in multigenerational and transgenerational effects. Several publications have reviewed such actions of EDCs in subgroups of species, e.g., fish or rodents. In this review, we take a holistic approach synthesizing knowledge of the effects of EDCs across vertebrate species, including fish, anurans, birds, and mammals, and discuss the potential mechanism(s) mediating such multi- and transgenerational effects. We also propose a series of recommendations aimed at moving the field forward in a structured and coherent manner.

Screening and Validation of p38 MAPK Involved in Ovarian Development of Brachymystax lenok

Brachymystax lenok (lenok) is a rare cold-water fish native to China that is of high meat quality. Its wild population has declined sharply in recent years, and therefore, exploring the molecular mechanisms underlying the development and reproduction of lenoks for the purposes of artificial breeding and genetic improvement is necessary. The lenok comparative transcriptome was analyzed by combining single molecule, real-time, and next generation sequencing (NGS) technology. Differentially expressed genes (DEGs) were identified in five tissues (head kidney, spleen, liver, muscle, and gonad) between immature [300 days post-hatching (dph)] and mature [three years post-hatching (ph)] lenoks. In total, 234,124 and 229,008 full-length non-chimeric reads were obtained from the immature and mature sequencing data, respectively. After NGS correction, 61,405 and 59,372 non-redundant transcripts were obtained for the expression level and pathway enrichment analyses, respectively. Compared with the mature group, 719 genes with significantly increased expression and 1,727 genes with significantly decreased expression in all five tissues were found in the immature group. Furthermore, DEGs and pathways involved in the endocrine system and gonadal development were identified, and p38 mitogen-activated protein kinases (MAPKs) were identified as potentially regulating gonadal development in lenok. Inhibiting the activity of p38 MAPKs resulted in abnormal levels of gonadotropin-releasing hormone, follicle-stimulating hormone, and estradiol, and affected follicular development. The full-length transcriptome data obtained in this study may provide a valuable reference for the study of gene function, gene expression, and evolutionary relationships in B. lenok and may illustrate the basic regulatory mechanism of ovarian development in teleosts.

The miRNome function transitions from regulating developmental genes to transposable elements during pollen maturation

Animal and plant microRNAs (miRNAs) are essential for the spatio-temporal regulation of development. Together with this role, plant miRNAs have been proposed to target transposable elements (TEs) and stimulate the production of epigenetically active small interfering RNAs. This activity is evident in the plant male gamete containing structure, the male gametophyte or pollen grain. How the dual role of plant miRNAs, regulating both genes and TEs, is integrated during pollen development and which mRNAs are regulated by miRNAs in this cell type at a genome-wide scale are unknown. Here, we provide a detailed analysis of miRNA dynamics and activity during pollen development in Arabidopsis thaliana using small RNA and degradome parallel analysis of RNA end high-throughput sequencing. Furthermore, we uncover miRNAs loaded into the two main active Argonaute (AGO) proteins in the uninuclear and mature pollen grain, AGO1 and AGO5. Our results indicate that the developmental progression from microspore to mature pollen grain is characterized by a transition from miRNAs targeting developmental genes to miRNAs regulating TE activity.

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.

A Comparison of Reproductive Performances in Young and Old Females: A Case Study on the Atlantic Bluefin Tuna in the Mediterranean Sea

Simple Summary

The Atlantic bluefin tuna Thunnus thynnus is a species characterized by complex trans-oceanic migrations linked to size, which rely on the delicate trade-off between somatic growth and reproduction before and during the migratory movements to reach spawning grounds. Therefore, understanding the processes that drive reproduction and elucidating its age-related regulation is essential in the context of sustainable fishery management. In this study, carried out in the Mediterranean Sea, older bluefin tuna females were found to have greater reproductive performances than younger females according to a molecular biology approach (i.e., gene expression), a result that likely mirrors a better physical condition, different habitat usage or migratory behaviour. This result highlights the importance of preserving large females for their major reproductive contribution at a stock level. Furthermore, the gonad-specific mir-202, which belongs to a class of non-coding RNA, called miRNA, that regulate the post-transcription of protein-coding genes, was identified as a potential candidate to play a role in egg quality and quantity (i.e., fecundity) during ovarian maturation through age- or stage-dependent reproductive processes. Overall, the present study contributes to improve the sustainability of the Atlantic bluefin tuna fishery in the Mediterranean Sea.

Abstract

In the Mediterranean Sea, a demographic substructure of the Atlantic bluefin tuna Thunnus thynnus has emerged over the last decade, with old and young individuals exhibiting different horizontal movements and spatial–temporal patterns of gonad maturation. In the present study, histology and molecular reproductive markers were integrated with the gonad-specific mir-202 gene expression and ovarian localization to provide a comprehensive picture of the reproductive performances in young and old females and investigate the role played by the mir-202 during gonadal maturation. During the reproductive period, old females (>100 kg; 194.6 ± 33.9 cm straight fork length; 11.3 ± 2.7 years old) were found to have greater reproductive performances than younger females (<80 kg; 139.3 ± 18.8 cm straight fork length; 8.4 ± 1.1 years old) according to gene expression results, suggesting a prolonged spawning season, earlier arrival on spawning grounds and/or better condition in older females. The mir-202-5p showed no global changes; it was abundantly expressed in granulosa cells and faintly present in the ooplasm. On the other hand, the mir-202-3p expression profile reflected levels of oocyte maturation molecular markers (star, lhr) and both histological and molecular (casp3) levels of follicular atresia. Overall, old females exhibited greater reproductive performances than younger females, likely reflecting different reproductive dynamics linked to the physical condition, habitat usage and migratory behaviour. These results highlight the importance of preserving large and old females in the context of fishery management. Finally, the mir-202 appears to be a good candidate to regulate the reproductive output of this species in an autocrine/paracrine manner through either stage- or age-dependent processes.

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.

GH Overexpression Alters Spermatic Cells MicroRNAome Profile in Transgenic Zebrafish

Overexpression of growth hormone (GH) in gh-transgenic zebrafish of a highly studied lineage F0104 has earlier been reported to cause increased muscle growth. In addition to this, GH affects a broad range of cellular processes in transgenic fish, such as morphology, physiology, and behavior. Reports show changes such as decreased sperm quality and reduced reproductive performance in transgenic males. It is hypothesized that microRNAs are directly involved in the regulation of fertility potential during spermatogenesis. The primary aim of our study was to verify whether gh overexpression disturbs the sperm miRNA profile and influences the sperm quality in transgenic zebrafish. We report a significant increase in body weight of gh-transgenic males along with associated reduced sperm motility and other kinetic parameters in comparison to the non-transgenic group. MicroRNA transcriptome sequencing of gh-transgenic zebrafish sperms revealed expressions of 186 miRNAs, among which six miRNA were up-regulated (miR-146b, miR-200a-5p, miR-146a, miR-726, miR-184, and miR-738) and sixteen were down-regulated (miR-19d-3p, miR-126a-5p, miR-126b-5p, miR-22a-5p, miR-16c-5p, miR-20a-5p, miR-126b-3p, miR-107a-3p, miR-93, miR-2189, miR-202–5p, miR-221–3p, miR-125a, miR-125b-5p, miR-126a-3p, and miR-30c-5p) in comparison to non-transgenic zebrafish. Some of the dysregulated miRNAs were previously reported to be related to abnormalities in sperm quality and reduced reproduction ability in other species. In this study, an average of 134 differentially expressed miRNAs-targeted genes were predicted using the in silico approach. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis demonstrated that the genes of affected pathways were primarily related to spermatogenesis, sperm motility, and cell apoptosis. Our results suggested that excess GH caused a detrimental effect on sperm microRNAome, consequently reducing the sperm quality and reproductive potential of zebrafish males.

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

BACKGROUND

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

RESULT

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

CONCLUSION

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

The brominated flame retardant, TBCO, impairs oocyte maturation in zebrafish (Danio rerio)

The brominated flame retardant, 1,2,5,6-tetrabromocyclooctane (TBCO), has been shown to decrease fecundity in Japanese medaka (Oryzias latipes) and there is indirect evidence from analysis of the transcriptome and proteome that this effect might be due to impaired oogenesis. An assay for disruption of oocyte maturation by chemical stressors has not been developed in Japanese medaka. Thus, using zebrafish (Danio rerio) as a model, objectives of the present study were to determine whether exposure to TBCO has effects on maturation of oocytes and to investigate potential mechanisms. Sexually mature female zebrafish were given a diet of 35.3 or 628.8 μg TBCO / g food for 14 days after which, stage IV oocytes were isolated to assess maturation in response to maturation inducing hormone. To explore potential molecular mechanisms, abundances of mRNAs of a suite of genes that regulate oocyte maturation were quantified by use of quantitative real-time PCR, and abundances of microRNAs were determined by use of miRNAseq. Ex vivo maturation of oocytes from fish exposed to TBCO was significantly less than maturation of oocytes from control fish. The percentage of oocytes which matured from control fish and those exposed to low and high TBCO were 89, 71, and 67%, respectively. Among the suite of genes known to regulate oocyte maturation, mRNA abundance of insulin like growth factor-3 was decreased by 1.64- and 3.44-fold in stage IV oocytes from females given the low and high concentrations of TBCO, respectively, compared to the control group. Abundances of microRNAs regulating the expression of proteins that regulate oocyte maturation, including processes related to insulin-like growth factor, were significantly different in stage IV oocytes from fish exposed to TBCO. Overall, results of this study indicated that impaired oocyte maturation might be a mechanism of reduced reproductive performance in TBCO-exposed fish. Results also suggested that effects of TBCO on oocyte maturation might be due to molecular perturbations on insulin-like growth factor signaling and expression of microRNAs.

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.

Integrated mRNA and miRNA expression profile analysis of female and male gonads in Hyriopsis cumingii

Hyriopsis cumingii is an important species for freshwater pearl cultivation in China. In terms of pearl production, males have larger pearls and better glossiness than females, but there are few reports focusing on the sex of H. cumingii. In this study, six mRNA and six microRNA (miRNA) libraries were prepared from ovaries and testes. Additionally, 28,502 differentially expressed genes (DEGs) and 32 differentially expressed miRNAs (DEMs) were identified. Compared with testis, 14,360 mRNAs and 20 miRNAs were up-regulated in ovary, 14,142 mRNAs and 12 miRNAs were down-regulated. In DEGs, the known genes related to sex determinism and/or differentiation were also identified, such as DMRT1, SOX9, SF1 for males, FOXL2 for females, and other potentially significant candidate genes. Three sex-related pathways have also been identified, which are Wnt, Notch, and TGF-beta. In 32 DEMs, the three miRNAs (miR-9-5p, miR-92, miR-184) were paid more attention, they predicted 28 target genes, which may also be candidates for sex-related miRNAs and genes. Differential miRNAs target genes analysis reveals the pathway associated with oocyte meiosis and spermatogenesis. Overall, the findings of the study provide significant insights to enhance our understanding of sex differentiation and/or sex determination mechanisms for H. cumingii.

The shift from early to late types of ribosomes in zebrafish development involves changes at a subset of rRNA 2'-O-Me sites

During zebrafish development, an early type of rRNA is gradually replaced by a late type that is substantially different in sequence. We applied RiboMeth-seq to rRNA from developmental stages for profiling of 2'-O-Me, to learn if changes in methylation pattern were a component of the shift. We compiled a catalog of 2'-O-Me sites and cognate box C/D guide RNAs comprising 98 high-confidence sites, including 10 sites that were not known from other vertebrates, one of which was specific to late-type rRNA. We identified a subset of sites that changed in methylation status during development and found that some of these could be explained by availability of their cognate SNORDs. Sites that changed during development were enriched in the novel sites revealed in zebrafish. We propose that the early type of rRNA is a specialized form and that its structure and ribose methylation pattern may be an adaptation to features of development, including translation of specific maternal mRNAs.

Kisspeptin Influences the Reproductive Axis and Circulating Levels of microRNAs in Senegalese Sole

Kisspeptin regulates puberty and reproduction onset, acting upstream of the brain-pituitary-gonad (HPG) axis. This study aimed to test a kisspeptin-based hormonal therapy on cultured Senegalese sole (G1) breeders, known to have reproductive dysfunctions. A single intramuscular injection of KISS2-10 decapeptide (250 µg/kg) was tested in females and males during the reproductive season, and gonad maturation, sperm motility, plasma levels of gonadotropins (Fsh and Lh) and sex steroids (11-ketotestosterone, testosterone and estradiol), as well as changes in small non-coding RNAs (sncRNAs) in plasma, were investigated. Fsh, Lh, and testosterone levels increased after kisspeptin injection in both sexes, while sperm analysis did not show differences between groups. Let7e, miR-199a-3p and miR-100-5p were differentially expressed in females, while miR-1-3p miRNA was up-regulated in kisspeptin-treated males. In silico prediction of mRNAs targeted by miRNAs revealed that kisspeptin treatment might affect paracellular transporters, regulate structural and functional polarity of cells, neural networks and intracellular trafficking in Senegalese sole females; also, DNA methylation and sphingolipid metabolism might be altered in kisspeptin-treated males. Results demonstrated that kisspeptin stimulated gonadotropin and testosterone secretion in both sexes and induced an unanticipated alteration of plasma miRNAs, opening new research venues to understand how this neuropeptide impacts in fish HPG axis.

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.

Male reproductive dysfunction in Solea senegalensis: new insights into an unsolved question

Senegalese sole (Solea senegalensis) is a species with a high commercial value that exhibits a reproductive dysfunction in males born and raised in captivity (F1) that hinders their sustainable culture. The present study evaluates the sperm quality and dopaminergic pathway of males born in the wild environment and of F1 males. Traditional sperm analyses were performed, finding only significant differences in curvilinear velocity (VCL) and no significant differences in viability and total motility. No differences in global sperm methylation were observed either in spermatozoa or brain between the two groups (F1 and wild-born males). However, our results point to a different sperm molecular signature between wild fish and fish born in captivity, specifically the differential expression in miR-let7-d and miR-200a-5p between these two groups. miR-let7-d has been correlated with spermatogenesis and sex preferences, whereas the miR-200 family is implied in target innervation of dopaminergic neurons in zebrafish. When we analysed the dopaminergic pathway, no differences were found in terms of different mRNA expression of dopaminergic markers. However, some differences were detected in terms of tyrosine hydroxylase protein expression by western blot analysis, thus suggesting an altered post-transcriptional regulation in F1 males. The results of this study suggest that an altered sperm miRNA signature in F1 males could be one possible mode of transmission of reproductive dysfunction to the progeny.

Distinct expression patterns of seven crucial microRNAs during early embryonic development in medaka (Oryzias latipes)

MicroRNAs (i.e. miRNAs) are small non-coding RNAs that play essential modulation roles in embryonic development in vertebrates. Paternal and maternal miRNAs contribute to the development of post-fertilization embryo and zygotic genome activation. The pattern of expression and their roles in embryonic development of medaka are not clearly understood. The present study, therefore, examined a temporal expression of seven miRNAs, ola-let-7a, ola-miR-202-3p, ola-miR-126-3p, ola-miR-122, ola-miR-92a, ola-miR-125a-3p and ola-miR-430a in sperm, oocytes, and embryos during early developmental stages. Three unique expression patterns of miRNAs were observed. ola-let7a, ola-miR-202-3p and ola-miR-126-3p showed both paternal and maternal expression, and ola-miR-122, ola-miR-92a, ola-miR-125a-3p showed maternal expression only. The expression of six out of seven miRNAs significantly decreased after maternal-zygotic transition (MZT), whereas ola-miR-430a expression initiated only after MZT. The temporal dynamic expression of these miRNAs suggests their potential roles in early embryogenesis and genome-zygotic activation in medaka.

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.

High-throughput sequencing analysis identified microRNAs associated with egg production in ducks ovaries

BACKGROUND

MicroRNAs (miRNAs) exist widely and are involved in multiple biological processes in ducks, whereas the regulatory mechanism of miRNAs in egg laying of ducks has remained unclear. This study aims to reveal key miRNAs involved in the regulation of egg production in duck ovaries.

METHODS

High-throughput sequencing was performed on four egg-type duck ovaries and four egg-meat-type duck ovaries at the start of the egg-laying stage. Quantitative reverse transcription PCR (qRT-PCR) validation was performed on differentially expressed miRNAs (DE miRNAs). Gene network of DEmiRNA-mRNA-pathway was constructed by Cytoscape.

RESULTS

A total of 251 know miRNAs and 1,972 novel miRNAs were obtained from whole clean reads. Among the known miRNAs, we identified 21 DEmiRNAs, including eight down-regulated and 13 up-regulated miRNAs in egg-type ducks compared with egg-meat-type ducks. Among the novel miRNAs, we identified 70 DEmiRNAs, including 58 down-regulated and 12 up-regulated in egg-type ducks compared with egg-meat-type ducks. The expression patterns of four miRNAs were verified by qRT-PCR. The DEmiRNAs were involved in the function of response to folic acid and the pathway of valine, leucine and isoleucine degradation. Specific target genes of DEmiRNAs enrichment was found in some egg-laying regulation pathways, such as dopaminergic synapse, ovarian steroidogenesis and oocyte meiosis. The DEmiRNA-mRNA-pathway network including three DEmiRNAs, nine mRNAs and 11 pathways. apl-miR-194-5p and apl-miR-215-5p may be potential key miRNAs in regulating egg laying.

CONCLUSIONS

This study provided miRNAs profiles in ducks about egg laying and establish a theoretical basis for subsequent selection or modification of duck phenotypes at the molecular level.

Identification and comparative profiling of gonadal microRNAs in the adult pigeon (Columba livia)

1. MicroRNAs are small noncoding RNA molecules that play crucial roles in gene expression. However, the comparative profiling of testicular and ovarian microRNAs in birds are rarely reported, particularly in pigeon.2. In this study, Illumina next-generation sequencing technology was used to sequence miRNA libraries of the gonads from six healthy adult utility pigeons. A total of 344 conserved known miRNAs and 32 novel putative miRNAs candidates were detected. Compared with those of ovaries, 130 differentially expressed (DE) miRNAs were identified in the testes. Among them, 70 miRNAs showed down-regulation in the ovaries, while another 60 miRNAs were up-regulated.3. Combining the results of the expression of target gene measurements and pathway enrichment analyses, it was revealed that some DEmiRNAs from the gonad samples involved in sexual differentiation and development (such as cli-miR-210-3p and cli-miR-214-3p) could down-regulate AR (androgen receptor). Cli-miR-181b-5p, cli-miR-9622-3p and cli-miR-145-5p were highly expressed in both the ovaries and testes, which could co-target HOXC9, and were related to regulation of primary metabolic processes. KEGG enrichment analysis showed that DEmiRNAs may play biological and sex-related roles in pigeon gonads.4. The expression profiles of testicular and ovarian miRNA in adult pigeon gonads are presented for the first time, and the findings may contribute to a better understanding of gonadal expression in poultry.

Developmental fluoxetine exposure in zebrafish reduces offspring basal cortisol concentration via life stage-dependent maternal transmission

Fluoxetine (FLX) is a pharmaceutical used to treat affective disorders in humans, but as environmental contaminant also affects inadvertently exposed fish in urban watersheds. In humans and fish, acute FLX treatment and exposure are linked to endocrine disruption, including effects on the reproductive and stress axes. Using the zebrafish model, we build on the recent finding that developmental FLX exposure reduced cortisol production across generations, to determine possible parental and/or life-stage-dependent (age and/or breeding experience) contributions to this phenotype. Specifically, we combined control and developmentally FLX-exposed animals of both sexes (F0) into four distinct breeding groups mated at 5 and 9 months, and measured offspring (F1) basal cortisol at 12 dpf. Basal cortisol was lower in F1 descended from developmentally FLX-exposed F0 females bred at 5, but not 9 months, revealing a maternal, life-stage dependent effect. To investigate potential molecular contributions to this phenotype, we profiled maternally deposited transcripts involved in endocrine stress axis development and regulation, epigenetic (de novo DNA methyltransferases) and post-transcriptional (miRNA pathway components and specific miRNAs) regulation of gene expression in unfertilized eggs. Maternal FLX exposure resulted in decreased transcript abundance of glucocorticoid receptor, dnmt3 paralogues and miRNA pathway components in eggs collected at 5 months, and increased transcript abundance of miRNA pathway components at 9 months. Specific miRNAs predicted to target stress axis transcripts decreased (miR-740) or increased (miR-26, miR-30d, miR-92a, miR-103) in eggs collected from FLX females at 5 months. Increased abundance of miRNA-30d and miRNA-92a persisted in eggs collected from FLX females at 9 months. Clustering and principal component analyses of egg transcript profiles separated eggs collected from FLX-females at 5 months from other groups, suggesting that oocyte molecular signatures, and miRNAs in particular, may serve as predictive tools for the offspring phenotype of reduced basal cortisol in response to maternal FLX exposure.

Ginsenoside Rg3 protects mouse leydig cells against triptolide by downregulation of miR-26a

BACKGROUND

Ginsenoside Rg3 has been reported to exert protection function on germ cells. However, the mechanisms by which Rg3 regulates apoptosis in mouse Leydig cells remain unclear. In addition, triptolide (TP) has been reported to induce infertility in male rats. Thus, this study aimed to investigate the protective effect of Rg3 against TP-induced toxicity in MLTC-1 cells.

METHODS

CCK-8, immunofluorescence assay, Western blotting and flow cytometry were used to detect cell proliferation and cell apoptosis, respectively. In addition, the dual luciferase reporter system assay was used to detect the interaction between miR-26a and GSK3β in MLTC-1 cells.

RESULTS

TP significantly inhibited the proliferation of MLTC-1 cells, while the inhibitory effect of TP was reversed by Rg3. In addition, TP markedly induced apoptosis in MLTC-1 cells via increasing the expressions of Bax, active caspase 3, Cyto c and active caspase 9, and decreasing the level of Bcl-2. However, Rg3 alleviated TP-induced apoptosis of MLTC-1 cells. Moreover, the level of miR-26a was obviously downregulated by Rg3 treatment. The protective effect of Rg3 against TP-induced toxicity in MLTC-1 cells was abolished by miR-26a upregulation. Meanwhile, dual-luciferase assay showed GSK3β was the direct target of miR-26a in MLTC-1 cells. Overexpression of miR-26a markedly decreased the level of GSK3β. As expected, upregulation of miR-26a could abrogate the protective effects of Rg3 against TP-induced cytotoxicity via inhibiting the expression of GSK3β.

CONCLUSION

These results indicated that Rg3 could protect MLTC-1 against TP by downregulation of miR-26a. Therefore, Rg3 might serve as a potential agent for the treatment of male hypogonadism.

Identification of Novel MicroRNAs and Characterization of MicroRNA Expression Profiles in Zebrafish Ovarian Follicular Cells

MicroRNAs (miRNAs) are small non-coding RNAs that regulate gene expression primarily at the post-transcriptional levels and thereby play important roles in regulating many physiological and developmental processes. Oocyte maturation in fish is induced by hormones produced from the hypothalamus, pituitary, and ovary. Gonadotropin-releasing hormone (GnRH) stimulates the secretion of luteinizing hormone (LH), which in turn, induces the secretion of maturation-inducing hormone (MIH) from the ovary. It is documented that small early vitellogenic (or stage IIIa) follicles are unable to undergo oocyte maturation whereas oocytes in mid- to late vitellogenic (stage IIIb) follicles can be induced by LH and MIH to become mature. To determine whether miRNAs may be involved in the growth and acquisition of maturational competency of ovarian follicles, we determined the miRNA expression profiles in follicular cells collected from stage IIIa and IIIb follicles using next-generation sequencing. It was found that miRNAs are abundantly expressed in the follicular cells from both stages IIIa and IIIb follicles. Furthermore, bioinformatics analysis revealed the presence of 214 known, 31 conserved novel and 44 novel miRNAs in zebrafish vitellogenic ovarian follicular cells. Most mature miRNAs in follicular cells were found to be in the length of 22 nucleotides. Differential expression analysis revealed that 11 miRNAs were significantly up-regulated, and 13 miRNAs were significantly down-regulated in the stage IIIb follicular cells as compared with stage IIIa follicular cells. The expression of four of the significantly regulated miRNAs, dre-miR-22a-3p, dre-miR-16a, dre-miR-181a-3p, and dre-miR-29a, was validated by real-time PCR. Finally, gene enrichment and pathway analyses of the predicted targets of the significantly regulated miRNAs supported the involvement of several key signaling pathways in regulating ovarian function, including oocyte maturation. Taken together, this study identifies novel zebrafish miRNAs and characterizes miRNA expression profiles in somatic cells within the zebrafish ovarian follicles. The differential expression of miRNAs between stage IIIa and IIIb follicular cells suggests that these miRNAs are important regulators of zebrafish ovarian follicle development and/or oocyte maturation.

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.

Changes in the morphology and gene expression of developing zebrafish gonads

Zebrafish gonadal sexual differentiation is an important but poorly understood subject. The difficulty in investigating zebrafish sexual development lies in its sex determination plasticity, the lack of morphological tools to distinguish juvenile females from males, and the lack of sex chromosomes in laboratory strains. Zebrafish sexual differentiation starts at around 8 days post-fertilization when germ cells start to proliferate. The number of germ cells determines the future sex of the gonad. Gonads with more germ cells differentiate into ovaries, whereas a reduced germ cell number leads to male-biased sexual differentiation. Genes controlling sexual differentiation in pre-meiotic gonads encode proteins such as transcription factors, the transforming growth factor (TGF)-β family of signaling proteins, and RNA-binding proteins. These proteins coordinately control germ cell proliferation/meiosis/maintenance and gonadal somatic cell differentiation, leading to stepwise differentiation of gonads. Morphological changes in differentiating gonads are characterized by the appearance of oocytes containing condensed chromatin, followed by incorporation of vitellogenin and oocyte maturation. Marker genes and morphological characteristics help distinguish the steps in zebrafish gonadal differentiation during this important sex-determining stage.

The molecular pathways underlying early gonadal development

The body of knowledge surrounding reproductive development spans the fields of genetics, anatomy, physiology and biomedicine, to build a comprehensive understanding of the later stages of reproductive development in humans and animal models. Despite this, there remains much to learn about the bi-potential progenitor structure that the ovary and testis arise from, known as the genital ridge (GR). This tissue forms relatively late in embryonic development and has the potential to form either the ovary or testis, which in turn produce hormones required for development of the rest of the reproductive tract. It is imperative that we understand the genetic networks underpinning GR development if we are to begin to understand abnormalities in the adult. This is particularly relevant in the contexts of disorders of sex development (DSDs) and infertility, two conditions that many individuals struggle with worldwide, with often no answers as to their aetiology. Here, we review what is known about the genetics of GR development. Investigating the genetic networks required for GR formation will not only contribute to our understanding of the genetic regulation of reproductive development, it may in turn open new avenues of investigation into reproductive abnormalities and later fertility issues in the adult.

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.

Epigenetics in teleost fish: From molecular mechanisms to physiological phenotypes

While the field of epigenetics is increasingly recognized to contribute to the emergence of phenotypes in mammalian research models across different developmental and generational timescales, the comparative biology of epigenetics in the large and physiologically diverse vertebrate infraclass of teleost fish remains comparatively understudied. The cypriniform zebrafish and the salmoniform rainbow trout and Atlantic salmon represent two especially important teleost orders, because they offer the unique possibility to comparatively investigate the role of epigenetic regulation in 3R and 4R duplicated genomes. In addition to their sequenced genomes, these teleost species are well-characterized model species for development and physiology, and therefore allow for an investigation of the role of epigenetic modifications in the emergence of physiological phenotypes during an organism's lifespan and in subsequent generations. This review aims firstly to describe the evolution of the repertoire of genes involved in key molecular epigenetic pathways including histone modifications, DNA methylation and microRNAs in zebrafish, rainbow trout, and Atlantic salmon, and secondly, to discuss recent advances in research highlighting a role for molecular epigenetics in shaping physiological phenotypes in these and other teleost models. Finally, by discussing themes and current limitations of the emerging field of teleost epigenetics from both theoretical and technical points of view, we will highlight future research needs and discuss how epigenetics will not only help address basic research questions in comparative teleost physiology, but also inform translational research including aquaculture, aquatic toxicology, and human disease.

miR-34a Regulates Sperm Motility in Zebrafish

Increasing attention has been focused on the role of microRNAs in post-transcription regulation during spermatogenesis. Recently, the miR-34 family has been shown to be involved in the spermatogenesis, but the clear function of the miR-34 family in spermatogenesis is still obscure. Here we analyzed the function of miR-34a, a member of the miR-34 family, during spermatogenesis using miR-34a knockout zebrafish generated by the clustered regularly interspaced short palindromic repeats/associated protein 9 (CRISPR/Cas9) system. miR-34a knockout zebrafish showed no obvious defects on testis morphology and sperm quantity. However, we found a significant increase in progressive sperm motility that is one of the pivotal factors influencing in vitro fertilization rates, in the knockout zebrafish. Moreover, breeding experiments showed that, when miR-34a-knockout male zebrafish mated with the wide-type females, they had a higher fertilization rate than did the wide-type males. Glycogen synthase kinase-3a (gsk3a), a potential sperm motility regulatory gene was predicted to be targeted by miR-34a, which was further supported by luciferase reporter assays, since a significant decrease of luciferase activity was detected upon ectopic overexpression of miR-34a. Our findings suggest that miR-34a downregulates gsk3a by targeting its 3' untranslated region, and miR-34a/gsk3a interaction modulates sperm motility in zebrafish. This study will help in understanding in the role of the miR-34 family during spermatogenesis and will set paths for further studies.