Integrated analysis of miRNA and mRNA expression profiles in tilapia gonads at an early stage of sex differentiation

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.

Whole transcriptome analysis revealed the regulatory network and related pathways of non-coding RNA regulating ovarian atrophy in broody hens

Poultry broodiness can cause ovarian atresia, which has a detrimental impact on egg production. Non-coding RNAs (ncRNAs) have become one of the most talked-about topics in life sciences because of the increasing evidence of their novel biological roles in regulatory systems. However, the molecular mechanisms of ncRNAs functions and processes in chicken ovarian development remain largely unknown. Whole-transcriptome RNA sequencing of the ovaries of broodiness and laying chickens was thus performed to identify the ncRNA regulatory mechanisms associated with ovarian atresia in chickens. Subsequent analysis revealed that the ovaries of laying chickens and those with broodiness had 40 differentially expressed MicroRNA (miRNAs) (15 up-regulated and 25 down-regulated), 379 differentially expressed Long Noncoding RNA (lncRNAs) (213 up-regulated and 166 down-regulated), and 129 differentially expressed circular RNA (circRNAs) (63 up-regulated and 66 down-regulated). The competing endogenous RNAs (ceRNA) network analysis further revealed the involvement of ECM-receptor interaction, AGE-RAGE signaling pathway, focal adhesion, cytokine-cytokine receptor interaction, inflammatory mediator regulation of TRP channels, renin secretion, gap junction, insulin secretion, serotonergic synapse, and IL-17 signaling pathways in broodiness. Upon further analysis, it became evident that THBS1 and MYLK are significant candidate genes implicated in the regulation of broodiness. The expression of these genes is linked to miR-155-x, miR-211-z, miR-1682-z, gga-miR-155, and gga-miR-1682, as well as to the competitive binding of novel_circ_014674 and MSTRG.3306.4. The findings of this study reveal the existence of a regulatory link between non-coding RNAs and their competing mRNAs, which provide a better comprehension of the ncRNA function and processes in chicken ovarian development.

Tilapia, a good model for studying reproductive endocrinology

The Nile tilapia (Oreochromis niloticus), with a system of XX/XY sex determination, is a worldwide farmed fish with a shorter sexual maturation time than that of most cultured fish. Tilapia show a spawning cycle of approximately 14 days and can be artificially propagated in the laboratory all year round to obtain genetically all female (XX) and all male (XY) fry. Its genome sequence has been opened, and a perfect gene editing platform has been established. With a moderate body size, it is convenient for taking enough blood to measure hormone level. In recent years, using tilapia as animal model, we have confirmed that estrogen is crucial for female development because 1) mutation of star2, cyp17a1 or cyp19a1a (encoding aromatase, the key enzyme for estrogen synthesis) results in sex reversal (SR) due to estrogen deficiency in XX tilapia, while mutation of star1, cyp11a1, cyp17a2, cyp19a1b or cyp11c1 affects fertility due to abnormal androgen, cortisol and DHP levels in XY tilapia; 2) when the estrogen receptors (esr2a/esr2b) are mutated, the sex is reversed from female to male, while when the androgen receptors are mutated, the sex cannot be reversed; 3) the differentiated ovary can be transdifferentiated into functional testis by inhibition of estrogen synthesis, and the differentiated testis can be transdifferentiated into ovary by simultaneous addition of exogenous estrogen and androgen synthase inhibitor; 4) loss of male pathway genes amhy, dmrt1, gsdf causes SR with upregulation of cyp19a1a in XY tilapia. Disruption of estrogen synthesis rescues the male to female SR of amhy and gsdf but not dmrt1 mutants; 5) mutation of female pathway genes foxl2 and sf-1 causes SR with downregulation of cyp19a1a in XX tilapia; 6) the germ cell SR of foxl3 mutants fails to be rescued by estrogen treatment, indicating that estrogen determines female germ cell fate through foxl3. This review also summarized the effects of deficiency of other steroid hormones, such as androgen, DHP and cortisol, on fish reproduction. Overall, these studies demonstrate that tilapia is an excellent animal model for studying reproductive endocrinology of fish.

microRNA-mRNA Analysis Reveals Tissue-Specific Regulation of microRNA in Mangrove Clam (Geloina erosa)

Geloina erosa is an important benthic animal in the mangrove, serving as an indicator organism for coastal environmental pollution. This study aimed to investigate the tissue-specific expression of miRNAs and their regulatory roles in predicted targets in G. erosa. Through miRNA sequencing and co-expression network analysis, we extensively studied the miRNA expression in three tissues: gills, hepatopancreas, and muscle. The results revealed a total of 1412 miRNAs, comprising 1047 known miRNAs, and 365 newly predicted miRNAs. These miRNAs exhibited distinct tissue-specific expression patterns. In the miRNA target gene prediction, a total of 7404 potential predicted targets were identified, representing approximately 33% of all unique transcripts associated with miRNAs. Further co-expression network analysis revealed nine modules, each showing a positive correlation with specific tissues (gills, hepatopancreas, or muscle). The blue module showed a significant correlation with gills (r = 0.83, p-value = 0.006), the black module was significantly related to the hepatopancreas (r = 0.78, p-value = 0.01), and the purple module was significantly correlated with muscle (r = 0.83, p-value = 0.006). Within these modules, related miRNAs tended to cluster together, while their correlations with other modules were relatively weak. Functional enrichment analysis was performed on miRNAs and their predicted targets in each tissue. In the gills, miRNAs primarily regulate immune-related genes, substance transport, and cytoskeletal organization. In the hepatopancreas, miRNAs suppressed genes involved in shell formation and played a role in cellular motor activity and metabolism. In muscle, miRNAs participate in metabolism and photoreceptive processes, as well as immune regulation. In summary, this study provides valuable insights into the tissue-specific regulation of miRNAs in G. erosa, highlighting their potential roles in immune response, metabolism, and environmental adaptation. These findings offer important clues for understanding the molecular mechanisms and biological processes in G. erosa, laying the foundation for further validation and elucidation of these regulatory relationships.

Molecular biological, physiological, cytological, and epigenetic mechanisms of environmental sex differentiation in teleosts: A systematic review

Human activities have been exerting widespread stress and environmental risks in aquatic ecosystems. Environmental stress, including temperature rise, acidification, hypoxia, light pollution, and crowding, had a considerable negative impact on the life histology of aquatic animals, especially on sex differentiation (SDi) and the resulting sex ratios. Understanding how the sex of fish responds to stressful environments is of great importance for understanding the origin and maintenance of sex, the dynamics of the natural population in the changing world, and the precise application of sex control in aquaculture. This review conducted an exhaustive search of the available literature on the influence of environmental stress (ES) on SDi. Evidence has shown that all types of ES can affect SDi and universally result in an increase in males or masculinization, which has been reported in 100 fish species and 121 cases. Then, this comprehensive review aimed to summarize the molecular biology, physiology, cytology, and epigenetic mechanisms through which ES contributes to male development or masculinization. The relationship between ES and fish SDi from multiple aspects was analyzed, and it was found that environmental sex differentiation (ESDi) is the result of the combined effects of genetic and epigenetic factors, self-physiological regulation, and response to environmental signals, which involves a sophisticated network of various hormones and numerous genes at multiple levels and multiple gradations in bipotential gonads. In both normal male differentiation and ES-induced masculinization, the stress pathway and epigenetic regulation play important roles; however, how they co-regulate SDi is unclear. Evidence suggests that the universal emergence or increase in males in aquatic animals is an adaptation to moderate ES. ES-induced sex reversal should be fully investigated in more fish species and extensively in the wild. The potential aquaculture applications and difficulties associated with ESDi have also been addressed. Finally, the knowledge gaps in the ESDi are presented, which will guide the priorities of future research.

Exosome delivery to the testes for dmrt1 suppression: A powerful tool for sex-determining gene studies

Exosomes are endosome-derived extracellular vesicles about 100 nm in diameter. They are emerging as promising delivery platforms due to their advantages in biocompatibility and engineerability. However, research into and applications for engineered exosomes are still limited to a few areas of medicine in mammals. Here, we expanded the scope of their applications to sex-determining gene studies in early vertebrates. An integrated strategy for constructing the exosome-based delivery system was developed for efficient regulation of dmrt1, which is one of the most widely used sex-determining genes in metazoans. By combining classical methods in molecular biology and the latest technology in bioinformatics, isomiR-124a was identified as a dmrt1 inhibitor and was loaded into exosomes and a testis-targeting peptide was used to modify exosomal surface for efficient delivery. Results showed that isomiR-124a was efficiently delivered to the testes by engineered exosomes and revealed that dmrt1 played important roles in maintaining the regular structure and function of testis in juvenile fish. This is the first de novo development of an exosome-based delivery system applied in the study of sex-determining gene, which indicates an attractive prospect for the future applications of engineered exosomes in exploring more extensive biological conundrums.

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.

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.

Comprehensive Characterization of Circular RNAs in Ovary and Testis From Nile Tilapia

Circular RNA (circRNA) is an endogenous biomolecule in eukaryotes. It has tissue- and cell-specific expression patterns and can act as a microRNA sponge or competitive endogenous RNA. Although circRNA has been found in several species in recent years, the expression profiles in fish gonad are still not fully understood. We detected the expression of circRNA in the ovary, testis, and sex-changed gonad of tilapia by high-throughput deep sequencing, and circRNA-specific computing tools. A total of 20,607 circRNAs were obtained, of which 141 were differentially expressed in the testis and ovary. Among these circRNAs, 135 circRNAs were upregulated and 6 circRNAs were downregulated in female fish. In addition, GO annotation and KEGG pathway analysis of the host genes of circRNAs indicated that these host genes were mainly involved in adherens junction, androgen production, and reproductive development, such as ZP3, PLC, delta 4a, ARHGEF10, and HSD17b3. It is worth noting that we found that circRNAs in tilapia gonads have abundant miRNA-binding sites. Among them, 935 circRNAs have a regulatory effect on miR-212, 856 circRNAs have a regulatory effect on miR-200b-3p, and 529 circRNAs have a regulatory effect on miR-200b-5p. Thus, our findings provide a new evidence for circRNA–miRNA networks in the gonads in tilapia.

Identification and Characterization of miRNAs and Their Predicted mRNAs in the Larval Development of Pearl Oyster Pinctada fucata

As an important economic shellfish, the pearl oyster, Pinctada fucata, and its larvae are an ideal model for studying molecular mechanisms of larval development in invertebrates. Larval development directly affects the quantity and quality of pearl oysters. MicroRNAs (miRNAs) may play important roles in development, but the effects of miRNA expression on P. fucata early development remain unknown. In this study, miRNA and mRNA transcriptomics of seven different P. fucata developmental stages were analyzed using Illumina RNA sequencing. A total of 329 miRNAs, including 87 known miRNAs and 242 novel miRNAs, and 33,550 unigenes, including 26,333 known genes and 7217 predicted new genes, were identified in these stages. A cluster analysis showed that the difference in the numbers of miRNAs was greatest between fertilized eggs and trochophores. In addition, the integrated mRNA transcriptome was used to predict target genes for differentially expressed miRNAs between adjacent developmental stages, and the target genes were subjected to a gene ontology enrichment analysis. Using the gene ontology annotation, 100 different expressed genes and 95 differentially expressed miRNAs were identified as part of larval development regulation. Real-time PCR was used to identify eight mRNAs and three miRNAs related to larval development. The present findings will be helpful for further clarifying the regulatory mechanisms of miRNA in invertebrate larval development.

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.

Highly Sensitive Detection of miR-21 through Target-Activated Catalytic Hairpin Assembly of X-Shaped DNA Nanostructures

MicroRNAs (miRNAs) are found in extremely low concentrations in cells, so highly sensitive quantitation is a great challenge. Herein, a simple dual-amplification strategy involving target-activated catalytic hairpin assembly (CHA) coupled with multiple fluorophores concentrated on one X-shaped DNA is reported. In this strategy, four hairpin probes (H1, H2, H3, and H4) are modified with FAM and BHQ1 at both sticky ends, while a circulating hairpin probe (H0) is used to activate CHA circuits once it binds to complementary sequences in the target miR-21 (T). The powerful dual-amplification cascades in Förster resonance energy transfer (FRET)-based nonenzymatic nucleic acid circuits are triggered by T-H0-activated formation of the X-shaped DNA nanostructure, freeing T-H0 for the next CHA reaction cycle. CHA circuits increase the fluorescence due to the wide distance between FAM and BHQ1 in the formed X-shaped DNA nanostructure, resulting in signal amplification and highly sensitive detection of miR-21, with a limit of detection (LOD, 3σ) of 0.025 nM, which is 25.6 or 57.6 times lower than that obtained through a single-amplification strategy without multiple fluorophores on one X-shaped DNA or CHA circuit. Furthermore, this cascade reaction was completed in 45 min, effectively avoiding target degradation. This new enzyme-free signal amplification strategy holds promising potential for sensitively detecting different DNA or RNA sequences by simply adapting the fragment of the H0 sequence complementary to the target.

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.

Integrated Analysis of mRNA- and miRNA-Seq in the Ovary of Rare Minnow Gobiocypris rarus in Response to 17α-Methyltestosterone

17α-Methyltestosterone (MT) is a synthetic androgen. The objective of this study was to explore the effects of exogenous MT on the growth and gonadal development of female rare minnow Gobiocypris rarus. Female G. rarus groups were exposed to 25–100 ng/L of MT for 7 days. After exposure for 7 days, the total weight and body length were significantly decreased in the 50-ng/L MT groups. The major oocytes in the ovaries of the control group were vitellogenic oocytes (Voc) and cortical alveolus stage oocytes (Coc). In the MT exposure groups, some fish had mature ovaries with a relatively lower proportion of mature oocytes, and the diameter of the perinucleolar oocytes (Poc) was decreased compared with those of the control group. Ovarian VTG, FSH, LH, 11-KT, E2, and T were significantly increased after exposure to 50 ng/L of MT for 7 days. Unigenes (73,449), 24 known mature microRNAs (miRNAs), and 897 novel miRNAs in the gonads of G. rarus were found using high-throughput sequencing. Six mature miRNAs (miR-19, miR-183, miR-203, miR-204, miR-205, and miR-96) as well as six differentially expressed genes (fabp3, mfap4, abca1, foxo3, tgfb1, and zfp36l1) that may be associated with ovarian development and innate immune response were assayed using qPCR. Furthermore, the miR-183 cluster and miR-203 were differentially expressed in MT-exposed ovaries of the different G. rarus groups. This study provides some information about the role of miRNA–mRNA pairs in the regulation of ovarian development and innate immune system, which will facilitate future studies of the miRNA–RNA-associated regulation of teleost reproduction.

Expression profile of miRNAs involved in the hepatoprotective effects of curcumin against oxidative stress in Nile tilapia

Curcumin is a polyphenol with antioxidant activity that has been used to protect the health of fish livers. Our previous studies about comparative transcriptome have shown that curcumin can enhance the Nrf2-Keap1 signaling pathway and induce downstream anti-stress genes to maintain cell viability. However, the possible role of miRNAs in the protective mechanism of curcumin is not understood. In this study, the tilapia hepatocyte H₂O₂ stress model was used, and the miRNA expression profile for four groups (control group, curcumin group, H₂O₂ group, and protection group) were established by high-throughput sequencing. In our results, 278-333 types of Oreochromis niloticus miRNAs, 309-543 types of conserved miRNAs, and 535-746 types of novel miRNAs were identified in different samples. Differentially expressed miRNAs were identified by comparing miRNA expression profiles among the four groups. The expression levels were confirmed by q-PCR. The target genes of these differentially expressed miRNAs were predicted, and their functional annotations were enriched by GO and KEGG analysis, which revealed that many target genes were involved in "response to stimulus" and "antioxidant activity" in each pair of groups. Several miRNAs related to oxidative stress showed differential expression. For example, in the H₂O₂ group, the expression of miR-122 was decreased, and the expression of miR-21 and miR-489 increased significantly. In the curcumin group, the expression of miR-153b was decreased, and the expression of miR-200a and miR-29 was increased significantly. miR-153b, miR-200a, and miR-29 may be involved in the regulation of the Nrf2-Keap1 signaling pathway by curcumin. This work might provide insights into the molecular mechanisms of miRNA regulation of curcumin on the prevention and alleviation of liver diseases in fish.

Global expression response of genes in sex-undifferentiated Nile tilapia gonads after exposure to trace letrozole

The aromatase inhibitor letrozole can be found in rivers, effluents, and even drinking water. Studies have demonstrated that letrozole affects various metabolic pathways and may cause reproductive toxicity, especially in fish exposed during development. However, studies on the effect of a low concentration of letrozole at the whole-gonad transcriptomic level in the early stage of fish sexual development have not been investigated. The aim of our study was to explore the potential effects of a low concentration of letrozole on the gonad transcriptome of Nile tilapia at an early stage of sexual development. In this study, 9 dpf (days postfertilization) Nile tilapia were exposed to trace letrozole for 12 days. Letrozole exposure from 9 dpf to 21 dpf persistently altered phenotypic sex development and induced the male-biased sex ratio. The transcriptome results showed that 1173 differentially expressed genes (DEGs) were present in the female control vs 1.5 μg/L letrozole-treated female comparison group and that 1576 DEGs were present in the 1.5 μg/L letrozole-treated female vs male control comparison group. Differentially expressed gene enrichment analysis revealed several crucial pathways, including the drug metabolism-cytochrome P450 pathway, the ErbB-PI3K/Akt/mTOR pathway, and the calcium signalling pathway. Further analysis of these identified DEGs indicated that some key genes correlated with metabolism and epigenetic regulation were significantly affected by letrozole, such as UDP-glucuronosyltransferase (Ugt), glutathione S-transferase omega-1 (Gsto1), lysine-specific demethylase 6bb (Kdm6bb, original name is Kdm6a), jumonji and AT-rich interaction domain containing 2 (Jarid2b, original name is Jarid2), growth arrest and DNA damage inducible gamma (Gadd45g), and chromobox protein 7 (Cbx7). The qRT-PCR validation results for twelve DEGs showed that the Pearson's correlation of the log₁₀^{fold change} values between the qPCR and RNA-Seq results was 0.90, indicating the accuracy and reliability of the RNA-Seq results. Our study is the first to report the effect of letrozole on the transcriptome of gonads from fish during early-stage sexual development. These findings will be useful for understanding the toxic effects and molecular mechanisms of letrozole exposure at the early stage of gonad development on the sexual development of aquatic organisms.

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 and Characterization of Osmoregulation Related MicroRNAs in Gills of Hybrid Tilapia Under Three Types of Osmotic Stress

Researchers have increasingly suggested that microRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate gene expression and protein translation in organs and respond to abiotic and biotic stressors. To understand the function of miRNAs in osmotic stress regulation of the gills of hybrid tilapia (Oreochromis mossambicus ♀ × Oreochromis urolepis hornorum ♂), high-throughput Illumina deep sequencing technology was used to investigate the expression profiles of miRNAs under salinity stress (S, 25‰), alkalinity stress (A, 4‰) and salinity-alkalinity stress (SA, S: 15‰, A: 4‰) challenges. The results showed that 31, 41, and 27 upregulated and 33, 42, and 40 downregulated miRNAs (P < 0.05) were identified in the salt stress, alkali stress, and saline-alkali stress group, respectively, which were compared with those in the control group (C). Fourteen significantly differently expressed miRNAs were selected randomly and then validated by a quantitative polymerase chain reaction. On the basis of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis, genes related to osmoregulation and biosynthesis were enriched in the three types of osmotic stress. In addition, three miRNAs and three predicted target genes were chosen to conduct a quantitative polymerase chain reaction in the hybrid tilapia and its parents during 96-h osmotic stress. Differential expression patterns of miRNAs provided the basis for research data to further investigate the miRNA-modulating networks in osmoregulation of teleost.

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.

High-quality chromosome-level genomes of two tilapia species reveal their evolution of repeat sequences and sex chromosomes

Tilapias are one of the most farmed fishes that are coined as "aquatic chicken" by the food industry. Nile tilapia and blue tilapia exhibit very recent transition of sex chromosome systems since their divergence approximately five million years ago, making them a great model for elucidating the molecular and evolutionary mechanisms of sex chromosome turnovers. Studies of their sex-determining pathways are also critical for developing genetic sex control in aquaculture. We report here the newly produced genomes of Nile tilapia and blue tilapia that integrate long-read sequencing and chromatin conformation data. The two nearly complete genomes have anchored over 97% of the sequences into linkage groups (LGs), and assembled majorities of complex repetitive regions including telomeres, centromeres and rDNA clusters. In particular, we inferred two episodes of repeat expansion at LG3 respectively in the ancestor of cichlids and that of tilapias. The consequential large heterochromatic region concentrated at one end of LG3 comprises tandem arrays of mRNA and small RNA genes, among which we have identified a candidate female determining gene Paics in blue tilapia. Paics shows female-specific patterns of single-nucleotide variants, copy numbers and expression patterns in gonads during early gonadogenesis. Our work provides a very important genomic resource for functional studies of cichlids, and suggested that unequal distribution of repeat content that impacts the local recombination rate might make some chromosomes more likely to become sex chromosomes.

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.

Genomics to accelerate genetic improvement in tilapia

Selective breeding of tilapia populations started in the early 1990s and over the past three decades tilapia has become one of the most important farmed freshwater species, being produced in more than 125 countries around the globe. Although genome assemblies have been available since 2011, most of the tilapia industry still depends on classical selection techniques using mass spawning or pedigree information to select for growth traits with reported genetic gains of up to 20% per generation. The involvement of international breeding companies and research institutions has resulted in the rapid development and application of genomic resources in the last few years. GWAS and genomic selection are expected to contribute to uncovering the genetic variants involved in economically relevant traits and increasing the genetic gain in selective breeding programs, respectively. Developments over the next few years will probably focus on achieving a deep understanding of genetic architecture of complex traits, as well as accelerating genetic progress in the selection for growth-, quality- and robustness-related traits. Novel phenotyping technologies (i.e. phenomics), lower-cost whole-genome sequencing approaches, functional genomics and gene editing tools will be crucial in future developments for the improvement of tilapia aquaculture.

miRTil: An Extensive Repository for Nile Tilapia microRNA Next Generation Sequencing Data

Nile tilapia is the third most cultivated fish worldwide and a novel model species for evolutionary studies. Aiming to improve productivity and contribute to the selection of traits of economic impact, biotechnological approaches have been intensively applied to species enhancement. In this sense, recent studies have focused on the multiple roles played by microRNAs (miRNAs) in the post-transcriptional regulation of protein-coding genes involved in the emergence of phenotypes with relevance for aquaculture. However, there is still a growing demand for a reference resource dedicated to integrating Nile Tilapia miRNA information, obtained from both experimental and in silico approaches, and facilitating the analysis and interpretation of RNA sequencing data. Here, we present an open repository dedicated to Nile Tilapia miRNAs: the “miRTil database”. The database stores data on 734 mature miRNAs identified in 11 distinct tissues and five key developmental stages. The database provides detailed information about miRNA structure, genomic context, predicted targets, expression profiles, and relative 5p/3p arm usage. Additionally, miRTil also includes a comprehensive pre-computed miRNA-target interaction network containing 4936 targets and 19,580 interactions.

Transcriptome Profiling and Analysis of Genes Associated with High Temperature-Induced Masculinization in Sex-Undifferentiated Nile Tilapia Gonad

Artificially high temperatures during critical thermosensitive periods (TSPs) can induce the sex reversal of Nile tilapia (Oreochromis niloticus) females into pseudomales; Nile tilapia is a GSD + TE (genotypic plus temperature effects) fish species. Previous studies have shown that water temperature affects the expression levels of many genes in the gonad or brain in various teleost species. However, few studies on the effect of temperature at the whole-gonad transcriptomic level in the early stage of sex differentiation have been reported in fish species exhibiting GSD + TE. In this study, RNA-Seq was performed to characterize the transcriptomic profile and identify genes exhibiting temperature- and sex-biased expressions in the Nile tilapia gonad at 21 dpf. A total of 42 genes were found to be associated with both high-temperature treatment and sex development, as the expression levels of these genes differed in both FC (female control) vs MC (male control) and FC vs FT (high temperature-treated females in the TSP). Among these genes, the transcriptional alterations of many male sex determination and differentiation genes, such as Dmrt1, Gsdf, and the DNA damage-inducible protein GADD45 alpha, suggested that the male pathway is initiated after high-temperature treatment and that its initiation may play a role in high temperature-induced masculinization in Nile tilapia. The qRT-PCR validation results for thirteen differentially expressed genes showed that the Pearson's correlation of the log₁₀ fold change values between the qPCR and RNA-Seq results was 0.70 (p < 0.001), indicating the accuracy and reliability of the RNA-Seq results. Our study provides insights into how high-temperature treatment induces the sex reversal of Nile tilapia females.

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.

Effects of Bisphenol A on Foxl2 gene expression and DNA damage in adult viviparous fish Goodeaatripinnis

Bisphenol A (BPA) is an emerging pollutant of global concern. Viviparous fish Goodea atripinnis is endemic to the Central Mexican Plateau where BPA was detected; however, few studies examined the influence of this chemical on native viviparous fish. The effects of BPA (sublethal dose) were determined on DNA integrity and Foxl2 expression in G. atripinnis gonads, and interactions of BPA with FOXL2 protein. Genotoxicity analysis revealed that % comets, at 14 and 28 days and comet tail length (at 14 days) were significantly higher in exposed compared to controls. In general, the % DNA tail was not markedly higher in BPA-treated fish; however, tail moment related to tail length exhibited significant increases in DNA damage. RT-qPCR assays showed Foxl2 overexpression after 14 and 28 days of exposure in females; while in males, Foxl2 was overexpressed after 28 days. In silico analysis demonstrated that BPA interacted with seven residues located in FOXL2 homeodomain. In summary, sublethal BPA doses induced DNA damage and changes in Foxl2 expression in gonadal cells of G. atripinnis, which may adversely affect reproduction in BPA-exposed wild populations. Foxl2 overexpression and BPA-FOXL2 interaction suggested alterations in processes involving Foxl2. Viviparous fish may thus serve as potential non-conventional models for assessing pollutants effects.

Do G-protein coupled estrogen receptor and bisphenol A analogs influence on Leydig cell epigenetic regulation in immature boar testis ex vivo?

Organotypic culture of testicular fragments from 7-day-old male pigs (Polish White Large) was used. Tissues were treated with an antagonist of G-protein coupled estrogen receptor (GPER) (G-15; 10 nM), and bisphenol A (BPA), and its analogs (TBBPA, TCBPA; 10 nM) alone or in combination and analyzed using electron and light (stainings for collagen fibers, lipid droplet and autophagy markers) microscopes. In addition, mRNA and protein abundances and localization of molecules required for miRNA biogenesis and function (Drosha, Exportin 5; EXPO5, Dicer, and Argonaute 2; AGO2) were assessed together with calcium ion (Ca²⁺) and estradiol concentrations. Regardless of GPER blockade and/or treatment with BPA, TBBPA and TCBPA, there were no changes in Leydig cell morphology. Also, there were no changes in lipid droplet content and distribution but there were changes in lipid and autophagy protein abundance. In the interstitial tissue, there was an increase of collagen content, especially after treatment with BPA analogs and G-15 + BPA. Independent of the treatment, there was downregulation of EXPO5 and Dicer genes but the Drosha and AGO2 genes were markedly upregulated as a result of treatment with G-15 + BPA and TCBPA, respectively. There was always a lesser abundance of EXPO5 and AGO2 proteins regardless of treatment. There was markedly greater abundances of Drosha after G-15 + BPA treatment, and this also occurred for Dicer after treatment with G-15 + TCBPA. Immunolocalization of miRNA proteins indicated there was a cytoplasmic-nuclear pattern in control and treated cells. There was an increase of Ca²⁺ concentrations after treatment with G-15 and BPA analogs. Estradiol secretion decreased after antagonist and chemical treatments when these were administered alone, however, there was an increase in estradiol secretion after treatment with combinations of these compounds.

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.

Comparative expression analysis of let-7 microRNAs during ovary development in Megalobrama amblycephala

As a critical regulator of gene expression, let-7 family miRNAs have been reported to be involved in multiple physiological processes. In this study, in order to elucidate the putative regulatory effect of let-7 miRNAs during fish gonadal development and to identify which member is crucial for this regulation, the expression of ten members including let-7a/b/c/d/e/f/g/h/i/j were quantified in ovary, pituitary, and brain tissues during the different ovarian developmental stages of blunt snout bream Megalobrama amblycephala. According to the data from analysis of expression patterns, let-7a showed the highest expression value in almost all the tested ovaries, pituitaries, and brains, with let-7b and let-7d moderately expressed, following by other let-7 miRNAs. In terms of the differential expression levels of ten let-7 miRNAs at each developmental stage, the results showed that let-7a/b/d/f/h expression gradually increased during the ovary development from stage I to V and dropped significantly at the phase VI in ovary tissues. However, the expression of let-7a/b/e/f in pituitary increased during the ovary development from stage I to IV and declined at stage V. Among all the let-7 miRNAs, let-7a/b/d had the highest expression and their expression patterns were consistent with the gonad development of M. amblycephala. Furthermore, the mostly predicted target genes of let-7 miRNAs are significantly enriched in signaling pathways closely related to gonadal development through KEGG enrichment analysis. These results indicate that let-7 miRNA members, especially let-7a/b/d, may play important roles in the regulation of ovary development in M. amblycephala through negatively regulating expression of their target genes.

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.

Profiling expression changes of genes associated with temperature and sex during high temperature-induced masculinization in the Nile tilapia brain

Fish sex-determining mechanisms can be classified as genotypic (GSD), temperature (TSD), or genotypic plus temperature effects (GSD+TE). Previous studies have shown that culturing water temperature during thermosensitive periods (TSP) could affect the expression of many genes in the gonad in some fish. However, few studies have focused on gene expression changes in the brain after temperature treatment during TSP in fish species. In this study, three families were developed by crossing XX neomales with XX females and one of them was used for transcriptome analysis. The results showed that a total of 105, 3164 and 4666 DEGs were respectively obtained in FC (female control) vs. FT (high temperature-treated females at TSP), FC vs. MC (male control), and MC vs. FT comparison groups. By profiling analysis, we show that the mRNA expression levels of 16 differentially expressed genes (DEGs) exhibited significant downregulation or upregulation after high temperature treatment and reached a similar level as that in MC. Among the 16 DEGs, LOC100699848 (lysine specific demethylase 6A) and Jarid2 contained JmjC domain, showing the possible important role of JmjC domain in response to temperature treatment in Nile tilapia. Kdm6b (lysine demethylase 6B) and Jarid2 have been shown to play important roles in reptile TSD, showing the relative conservation of underlying regulation mechanisms between TSD in reptile and TSD or GSD+TE in fish species. Finally, the transcriptome profiling was validated by quantitative real-time PCR in nine selected genes. These results provide a direction for investigating the GSD+TE molecular mechanism in fish species.

Comparative transcriptomic analysis of hypothalamus-pituitary-liver axis in bighead carp (Hypophthalmichthys nobilis) with differential growth rate

Background

Growth rate is one of the most important features for aquaculture species and deciphering its regulation mechanism has great significance both in genetics and in economics. Hypothalamus-pituitary growth axis (HP growth axis) or neuro-endocrine axis plays a vital role in growth regulation in different aquaculture animals.

Results

In this study, the HP and liver transcriptomes of two female groups (H and L) with phenotypically extreme growth rate were sequenced using RNA-Seq. A total of 30,524 and 22,341 genes were found expressed in the two tissues, respectively. The average expression levels for the two tissues were almost the same, but the median differed significantly. A differential expression analysis between H and L groups identified 173 and 204 differentially expressed genes (DEGs) in HP and liver tissue, respectively. Pathway analysis revealed that DEGs in HP tissue were enriched in regulation of cell proliferation and angiogenesis while in liver tissue these genes were overrepresented in sterol biosynthesis and transportation. Genomic overlapping analyses found that 4 and 5 DEGs were within growth-related QTL in HP and liver tissue respectively. A deeper analysis of these 9 genes revealed 3 genes were functionally linked to the trait of interest. The expression of 2075 lncRNAs in HP tissue and 1490 in liver tissue were also detected, and some of lncRNAs were highly expressed in the two tissues.

Conclusions

Above all, the results of the present study greatly contributed to the knowledge of the regulation of growth and then assisted the design of new selection strategies for bighead carp with improved growth-related traits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5691-4) contains supplementary material, which is available to authorized users.

Cluster expansion of apolipoprotein D (ApoD) genes in teleost fishes

BACKGROUND

Gene and genome duplication play important roles in the evolution of gene function. Compared to individual duplicated genes, gene clusters attract particular attention considering their frequent associations with innovation and adaptation. Here, we report for the first time the expansion of the apolipoprotein D (ApoD) ligand-transporter genes in a cluster manner specific to teleost fishes.

RESULTS

Based on comparative genomic and transcriptomic analyses, protein 3D structure comparison, positive selection detection and breakpoints detection, the single ApoD gene in the ancestor expanded into two clusters following a dynamic evolutionary pattern in teleost fishes. Orthologous genes show conserved expression patterns, whereas lineage-specific duplicated genes show tissue-specific expression patterns and even evolve new gene expression profiles. Positive selection occurred in branches before and after gene duplication, especially for lineage-specific duplicated genes. Cluster analyses based on protein 3D structure comparisons, especially comparisons of the four loops at the opening side, show gene duplication-segregating patterns. Duplicated ApoD genes are predicted to be associated with forkhead transcription factors and MAPK genes. ApoD clusters are located next to the breakpoints of genome rearrangements.

CONCLUSIONS

Here, we report the expansion of ApoD genes specific to teleost fishes in a cluster manner for the first time. Neofunctionalization and subfunctionalization were observed at both the protein and expression levels after duplication. Evidence from different aspects-i.e., abnormal expression-induced disease in humans, fish-specific expansion, predicted associations with forkhead transcription factors and MAPK genes, specific expression patterns in tissues related to sexual selection and adaptation, duplicated genes under positive selection and their location next to the breakpoints of genome rearrangements-suggests the potentially advantageous roles of ApoD genes in teleost fishes. The cluster expansion of ApoD genes specific to teleost fishes provides thus an ideal evo-devo model for studying gene duplication, cluster maintenance and new gene function emergence.

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.

Genome-wide microRNA screening in Nile tilapia reveals pervasive isomiRs' transcription, sex-biased arm switching and increasing complexity of expression throughout development

MicroRNAs (miRNAs) are key regulators of gene expression in multicellular organisms. The elucidation of miRNA function and evolution depends on the identification and characterization of miRNA repertoire of strategic organisms, as the fast-evolving cichlid fishes. Using RNA-seq and comparative genomics we carried out an in-depth report of miRNAs in Nile tilapia (Oreochromis niloticus), an emergent model organism to investigate evo-devo mechanisms. Five hundred known miRNAs and almost one hundred putative novel vertebrate miRNAs have been identified, many of which seem to be teleost-specific, cichlid-specific or tilapia-specific. Abundant miRNA isoforms (isomiRs) were identified with modifications in both 5p and 3p miRNA transcripts. Changes in arm usage (arm switching) of nine miRNAs were detected in early development, adult stage and even between male and female samples. We found an increasing complexity of miRNA expression during ontogenetic development, revealing a remarkable synchronism between the rate of new miRNAs recruitment and morphological changes. Overall, our results enlarge vertebrate miRNA collection and reveal a notable differential ratio of miRNA arms and isoforms influenced by sex and developmental life stage, providing a better picture of the evolutionary and spatiotemporal dynamics of miRNAs.

Transcriptome-seq provides insights into sex-preference pattern of gene expression between testis and ovary of the crucifix crab (Charybdis feriatus)

The crucifix crab, Charybdis feriatus, which mainly inhabits Indo-Pacific region, is regarded as one of the most high-potential species for domestication and incorporation into the aquaculture sector. However, the regulatory mechanisms of sex determination and differentiation of this species remain unclear. To identify candidate genes involved in sex determination and differentiation, high throughput sequencing of transcriptome from the testis and ovary of C. feriatus was performed by the Illumina platform. After removing adaptor primers, low-quality sequences and very short (<50 nt) reads, we obtained 80.9 million and 66.2 million clean reads from testis and ovary, respectively. A total of 86,433 unigenes were assembled, and ~43% (37,500 unigenes) were successfully annotated to the NR, NT, Swiss-Prot, KEGG, COG, GO databases. By comparing the testis and ovary libraries, we obtained 27,636 differentially expressed genes. Some candidate genes involved in the sex determination and differentiation of C. feriatus were identified, such as vasa, pgds, vgr, hsp90, dsx-f, fem-1, and gpr. In addition, 88,608 simple sequence repeats were obtained, and 61,929 and 77,473 single nucleotide polymorphisms from testis and ovary were detected, respectively. The transcriptome profiling was validated by quantitative real-time PCR in 30 selected genes, which showed a good consistency. The present study is the first high-throughput transcriptome sequencing of C. feriatus. These findings will be useful for future functional analysis of sex-associated genes and molecular marker-assisted selections in C. feriatus.

Comparative microRNA-seq Analysis Depicts Candidate miRNAs Involved in Skin Color Differentiation in Red Tilapia

Differentiation and variation in body color has been a growing limitation to the commercial value of red tilapia. Limited microRNA (miRNA) information is available on skin color differentiation and variation in fish so far. In this study, a high-throughput Illumina sequencing of sRNAs was conducted on three color varieties of red tilapia and 81,394,491 raw reads were generated. A total of 158 differentially expressed miRNAs (|log₂(fold change)| ≥ 1 and q-value ≤ 0.001) were identified. Target prediction and functional analysis of color-related miRNAs showed that a variety of putative target genes—including slc7a11, mc1r and asip—played potential roles in pigmentation. Moreover; the miRNA-mRNA regulatory network was illustrated to elucidate the pigmentation differentiation, in which miR-138-5p and miR-722 were predicted to play important roles in regulating the pigmentation process. These results advance our understanding of the molecular mechanisms of skin pigmentation differentiation in red tilapia.

Identification and differential expression of microRNAs in testis and ovary of Amur sturgeon (Acipenser schrenckii)

BACKGROUND

MicroRNAs (miRNAs) cooperate with sex-related genes in post-transcriptional regulation and play extremely important roles in the establishment of sexually dimorphic traits in animals. However, the gonad miRNAs and expression patterns of miRNAs in sturgeon have not been investigated.

METHODS

In the present study, we used high-throughput small RNA sequencing (RNA-Seq) to discover gonad miRNAs from the ovaries and testes of Amur sturgeons (Acipenser schrenckii). Further, microarray and real-time PCR assays were performed to identify the expression patterns of gonad miRNAs.

RESULTS

As a result, a total of 679 conserved and 51 novel miRNAs were successfully discovered in the gonads of A. schrenckii. Moreover, we found wide sequence variations (isomiRs) in gonad miRNAs, including 5' and 3' isomiRs. Our microarray analysis further characterized the 730 miRNAs expression profiles, which indicated that 117 differentially expressed miRNAs were detected with sex-biased patterns: 71 testis-biased and 46 ovary-biased miRNAs. Based on bioinformatics prediction, we found that there were functional differences between the testis-biased and ovary-biased miRNA targets involved in reproductive-related GO and KEGG pathways. Further, the association of the differentially expressed miRNAs and sex-related target mRNAs was uncovered. Finally, the expression patterns of 11 sex-biased miRNAs and 7 sex-related targets were validated in testes and ovaries using real-time PCR. Putative, negatively expressed miRNA-mRNA relationships were confirmed, such as Dmrt1 and asc-miR-2779, AR and asc-miR-203b-3p, foxl2 and asc-miR-30d.

CONCLUSION

This study provides information regarding the gonad miRNAs in sturgeon. The differential expression miRNAs in the gonads will help us to further understand the role of miRNA-mediated post-transcriptional regulation in the ovary and testis of sturgeon.

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.