- 1Key Laboratory of Integrated Rice-fish Farming, Ministry of Agriculture and Rural Affairs, Shanghai Ocean University, Shanghai, China
- 2Key Laboratory of Exploration and Utilization of Aquatic Genetic Resources, Ministry of Education, Shanghai Ocean University, Shanghai, China
- 3Institute of Eco-Environmental Protection, Shanghai Academy of Agricultural Sciences, Shanghai, China
- 4College of Life Sciences, Huzhou University, Zhejiang, China
- 5Suqian Institute of Agricultural Sciences, Jiangsu Academy of Agricultural Sciences, Suqian, China
- 6Huzhou Fengshengwan Aquatic Seed Industry Co. Ltd., Zhejiang, China
Macrobrachium rosenbergii (M. rosenbergii), as a species of common prawn, is a delicacy that is consumed all over the world. By interacting with the target gene 3′-untranslated region (3'-UTR), microRNAs (miRNAs) regulate its expression and ultimately participate in the regulation of reproductive development. However, research focusing on miRNA regulation during gonadal development in M. rosenbergii received very little attention. To explore the association between miRNA and reproduction, we performed RNA sequencing (RNA-seq) on brain and gonad organs in male and female M. rosenbergii. A total of 494 miRNAs were obtained in RNA-seq, including 31 and 59 differentially expressed (DE) miRNAs in the brain and gonads, respectively. Furthermore, 9 DE miRNAs were randomly selected from the brain and gonads, and qRT-PCR was conducted to validate the results of RNA-seq. Interestingly, dpu-miR-133 was found to be substantially expressed in the male brain and testis but poorly expressed in the female brain, ovary, and other organs. Analysis of dpu-miR-133 by Targetscan and MiRanda predicted to target 5-HT1. Furthermore, the dual-luciferase reporter assay manifested that dpu-miR-133 can combine with 5-HT1. Overall, our research work provides basic data for further study on the miRNA-mediated regulation of brain, gonad, and reproductive development of study M. rosenbergii.
1 Highlights
RNA sequencing was performed and analyzed on the adult brain and gonads of M. rosenbergii.
In total, 31 and 59 differentially expressed miRNAs were identified from the brain and gonad of M. rosenbergii, respectively.
Dual-luciferase reporter assay proved that a testis biased expressed miRNA-133 targeting 5-HT1.
2 Introduction
MicroRNAs (miRNAs) are small, non-coding, single-stranded RNA molecules that are produced by transcription from intragenic or intergenic genomic sequences. miRNAs are first transcribed into primary miRNA, which is then cut by Drosha enzyme to hairpin precursor miRNA, and finally in the cytoplasm, mature miRNA is cleaved with the help of Dicer enzyme (O'Brien et al., 2018). miRNAs as considered key post-transcriptional regulators, specifically target the 3'-UTR of mRNA to cleave it or to restrict translation, which in turn regulates several physiological and developmental events (Bartel, 2004; Brennecke et al., 2005). Interestingly, each mRNA has recognition sequences for many miRNAs, and a single miRNA can target thousands of them (Miranda et al., 2006; Zhu et al., 2015). miRNAs are involved to regulate several aspects of life, such as cell communication, metabolism regulation, and cell differentiation (Betel et al., 2008). Furthermore, recent studies indicate that miRNAs are essential for the germ cell differentiation and the development of the reproductive system (He et al., 2012; Hossain et al., 2012; Meng et al., 2018; Zou et al., 2020).
Gamete growth as well as production in the gonads is a prerequisite for reproduction, which entails an intricate and exactly harmonized process, consisting of changes in the morphology and function of different kinds of follicles or spermatogenesis cells (Piprek, 2010). In order to control mechanisms, the expression of many DEGs must be strictly regulated at the transcriptional and post-transcriptional stages. In Decapod Crustaceans, numerous research has been carried out on transcriptional control of gonadal development as well as gametophyte occurrence (Uawisetwathana et al., 2011; He et al., 2012; Qiao et al., 2012; Meng et al., 2015; Peng et al., 2015). Being the most significant post-transcriptional regulator, there has widespread attention on miRNAs in the field of reproduction in recent years. There is growing evidence that miRNAs are vital for all stages of reproductive development, consisting of follicular genesis, spermatogenesis, and steroid production (Macias et al., 2009; Hossain et al., 2012). Biogenesis of miRNAs was inhibited during the initial phases of proliferation and differentiation, which led to cause problems with the maturation of oocytes and disruption of spermatogenesis (Murchison et al., 2007; Hayashi et al., 2008). Currently, miRNAs related to gonads are found in Eriocheir sinensis and Japanese penaeus (Song et al., 2014; He et al., 2015).
Apart from the gonads, the brain which is a part of the central nervous system is also concerned with the development of the gonads. Discrepancies in the expression of brain genes may lead to sex development differences (Dennis, 2004; Davies and Wilkinson, 2006). So, brain events are likely to regulate the development of the gonads (Francis, 1992). Recent studies have shown that nearly half of the known miRNAs are expressed in mammalian brains (Krichevsky et al., 2003; Landgraf et al., 2007). In Paralichthys olivaceus, transcriptome analyses find that there are differential gender bias genes in female and male brains, and nearly half of the genes are targeted by the DEMs during the gonadal differentiation (Zou et al., 2020). So far, information on the role of miRNAs is limited in gonadal differentiation. In short, it is very important to study the miRNAs and target genes in the brain.
The main mechanisms of sex determination in crustaceans are genotype, environmental, and genotype sex determination influenced by the environment (Vogt, 2018), Although there is no direct evidence to confirm the existence of sex chromosome in M. rosenbergii, it can be inferred from previous studies that M. rosenbergii belongs to WZ/ZZ chromosome sex determination type (Staelens et al., 2008). It is noteworthy that the Androgenic gland hormone (AGH) secreted by the Androgenic gland in M. rosenbergii is not only a sex-determining factor but also a regulatory factor of male sex differentiation. Sex-related genes of M. rosenbergii have been screened by cDNA library construction and transcriptome sequencings, such as MR-IAG(Ventura et al., 2009), MR-IR (Sharabi et al., 2016), Mr-Mrr (Phoungpetchara et al., 2012) and MRPINK(Cao et al., 2007). Researchers have conducted transcriptome sequencing analysis of ovary and testis from18 M. rosenbergii based on 454 sequencing platforms and obtained more than 750,000 high-quality clean reads and 44,407 Contigs. Among them, 112 and 270 genes were differentially expressed in ovary and testis of M. rosenbergii, respectively, suggesting that these genes may be related to sex differentiation (Jung et al., 2016).
Macrobrachium rosenbergii (M. rosenbergii) is one of the largest freshwater prawn species that live in both fresh and saltwater. Breeding of this species originated in Malaysia in the early 1960s and is now widely cultivated in South and Western Pacific and Southeast Asia (Bonami and Widada, 2011; Jiang and Qiu, 2013). The growth rate of female shrimps decreases due to the need to provide a large number of nutrients for ovarian development, while the male shrimps can keep growing faster all the time. Therefore, under the same cultural conditions, the average body weight of male shrimps is about twice that of female shrimps. While improvements in productivity and technology have led to huge advances, many development problems have also emerged due to inbreeding and bacterial contamination in breeding and farming techniques. Therefore, it is necessary to conduct in-depth research on the gonadal development mechanism of M. rosenbergii and provide a basis for breeding and related technologies optimization.
In this study, for RNA-seq, we collected the brain and gonadal tissues of male and female M. rosenbergii individuals, and the results were verified by qRT-PCR. Results obtained from Targetscan and MiRanda predicted the targeting relationship between dpu-miR-133 and 5-HT1. Furthermore, the dual-luciferase reporter gene test confirmed that dpu-miR-133 could bind to 5-HT1. These results are helpful for further research into the miRNA-mediated regulation of M. rosenbergii brain, gonad, and reproduction.
3 Materials and methods
3.1 Animals and sample collection
M. rosenbergii in good health were collected from Huzhou Fengsheng Bay Aquatic Products Co., LTD. and shipped back to our lab. These prawns were reared in a water circulation system at a temperature of 26°C with 14 h of light and 10 h of dark light on the campus of Shanghai Ocean University (Lin’gang, Shanghai, China). At the time of sampling, M. rosenbergii were four months old and had reached the stage of sexual differentiation. After dissection, samples of brain, gonads, muscles, liver, intestines, heart, and gills from both female and male shrimp are frozen in liquid nitrogen at once and kept at −80°C. Research in M. rosenbergii was conducted strictly in accordance with the regulations of the Ethics Committee for Experimental Animal Research of Shanghai Ocean University.
3.3 Construction of small RNA library
By using TRizol® (Invitrogen, Carlsbad, CA), total RNA was isolated from the brain and gonad of M. rosenbergii. Verification of RNA integrity by Agilent 2,100 Bioanalyzer (Agilent Technologies, United States). Small RNA libraries were prepared according to the TruSeq Small RNA Sample Prep Kit (Illumina, San Diego, United States) using 1 μg total RNA. The four libraries were then sequenced using Illumina Solexa technology at OE Biotech. Co., Ltd. (Shanghai, China).
3.4 Analysis of sequencing data and miRNAs identification
The data processing is carried out following the described procedure (Yue et al., 2016). Screen the original readings by wiping out low-mass sequences, splicing common RNA families (rRNA, tRNA, snRNA, snoRNA), and repeating sequences. Then we applied strict principles to identify conservative and new miRNAs, consisting of sequences between 18 and 26 nt in length with >10 reads. In addition, conserved miRNAs did not allow seed sequence mismatch, with any more than two mismatches at other locations. All new miRNAs had a similar stem-like structure to that of typical miRNA precursors. The RNAfold software (http://rna.tbi.univie.ac.at//cgi-bin/RNAWebSuite/RNAfold.cgi) was used for the prediction of the secondary structure of new sequences. miRNAs readings were then standardized for better comparisons, assessing differences in MiRNAs between males and females. Target genes were determined by TargetScan and Miranda software to better understand the differentially expressed miRNAs. Then we use GO and KEGG enrichment analyses.
3.5 miRNA Extraction and qRT-PCR
The PrimeScript™ First Strand cDNA Synthesis Kit (Takara, Shiga, Japan) was used to reverse transcribe total RNA into cDNA. Using Primer Premier 6.0 to design the divergence primer (Supplementary Table S1). The miRNA was subjected to a reverse transcription PCR (RT-PCR) chain reaction. The reaction mixture is made up of a 1 μL cDNA template, 0.5 μL reverse primers 10 mm, and 18 μL premixed Taq (Takara Taq version2.0 + dye). A total of 35 cycles, each containing 95°C denatured 15 s, 60°C annealed 15 s, and 72°C extended 15 s, were used to obtain PCR product agarose gel DNA extraction kit (Takara, Shiga, Japan) from a 1.5% agarose gel using MiniBEST, and further verified by Sanger sequencing. To study the expression of miRNAs in M. rosenbergii, 9 DE miRNAs in the brain and gonad were randomly selected. Using ABI7500 real-time fluorescence quantitative polymerase chain reaction system and TB Green®PreMix Ex Taq™II (Japan), we verified the tissue distribution of nine miRNAs in the muscle, liver, gill, heart, intestine, both male and female-brain, ovary, and testis of the adult of M. rosenbergii. The qRT-PCR amplification procedure is consistent with previous descriptions (Sun et al., 2020). Following amplification, the products were evaluated using dissociation curves. With U6 as the endogenous reference gene, the expression level was normalized by the 2−ΔΔCt method. IBM SPSS Statistics 25.0 software was used to test the association between qRT-PCR and RNA-seq.
3.6 Prediction of miR-133 targeting gene for 5-HT1
Depending on seed region pairing and minimal free energy, miR-133 probable targets were predicted using TargetScan (Lewis et al., 2003) and bioinformatics analysis of RNA22 v2 (http://www.mybiosoftware.com/).
3.7 Dual-luciferase reporter assay
miR-133 has a binding site to 5-HT1 3'-UTR. To validate miR-133 targets, a luciferase reporter plasmid was constructed. Clone a miR-133 binding sequence containing 5-HT1 3'-UTR and insert the pmirGLO vector. By mutating the predicted miR-133 seed region, a 5-HT1 3'-UTR mutant reporter gene was cloned by PCR, and the primers were shown in Supplementary Table S1. For 5-HT1 mutant, the binding sequence (CCGGTTGAGAAGGACCAAATGCAGTTA) was replaced by CCGCATGTGAAGCTGCAAATGCAGTTA using the primers of 5-HT1-MT-FW and 5-HT1-MT-RV, which was conducted into pmirGLO vector through XbaI and SalI sites. The plasmids were sequence-verified. The culture and transfection methods of human embryonic kidney 293 (HEK293) cells used in the experiment were as previously described (Sun et al., 2020). HEK293 cells were co-transfected with wild-type (WT) or mutant-type (MT) plasmids and miR-133 mimics (GenePharma, China) using FuGENE®HD (Promega, United States). The Dual-Glo® Luciferase Assay System (Promega, United States) was used to quantify the luciferase activity 36 h after transfection. At least three times each experiment was conducted.
4 Result
4.1 Characterization of gonadal small RNA
In this research, using a sample of ovary, testis, and brain of both gender, four libraries were created from small RNA. Initially, the four libraries generated raw reads, 32.84, 47.58, 28.63, and 29.78 M. After removing low-quality sequences, it had 22.18 M, 27.25 M, 24,49 M, and 24.65 M clean reads were included in the relevant libraries. In the library, non-coding RNAs are categorized based on their sequence similarity to short RNA sequences in Rfam: rRNA, snRNA, tRNA, known miRNA, and un-annotated. (Table 1). To identify miRNAs in the brain and gonads of M. rosenbergii, we compared the sequences of the brain and gonads of M. rosenbergii with that of mature miRNAs in miRBase21.0, resulting in 271 conserved miRNAs and 223 new miRNAs (Table 2). In this study, we found dpu-miR-100 with 3,416,104 reads in the female brain and dpu-miR-1 with 775,001 reads in the testis of males were the most abundant and conservative miRNAs. In addition, tcf-miR-184-3p, dpu-miR-71, and dpu-miR-276 were captured over a hundred thousand times. Furthermore, the box beard plot was used to understand the distribution of miRNA in the four libraries (Supplementary Figure S1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990677/fgene-13-990677-HTML/image_m/fgene-13-990677-t001.jpg)
TABLE 1. Reads distribution of four small RNA libraries. Note: FB, female brain; MB, male brain; O, ovary; T, testis.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990677/fgene-13-990677-HTML/image_m/fgene-13-990677-t002.jpg)
TABLE 2. Statistics of known and novel miRNAs in four small RNA libraries. Note: FB, female brain; MB, male brain; O, ovary; T, testis.
4.2 Distribution of miRNA length
The distribution of miRNA length presents a distinctly different bimodal pattern in terms. According to an analysis of the length distribution, 22 nt was the main size, and 21 and 23 nt is followed, which is the same as the typical characteristics of Disher enzymes. Another peak for 32 nt mainly represents endogenous small non-coding RNA molecules, called longer piwi-interacting RNAs (piRNAs). As shown, the peak represented by the 21–23 nt size class was much higher than the peaks of the longest and shortest size class, indicating the abundance of small RNAs in the brain and gonad of M. rosenbergii (Figure 1A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990677/fgene-13-990677-HTML/image_m/fgene-13-990677-g001.jpg)
FIGURE 1. Identification of miRNAs. (A) The distribution of miRNA length from the M. rosenbergii. (B) Statistic of differently expressed miRNA. (C) Volcanic map of M. rosenbergii differences in the brain. (D) Volcanic map of M. rosenbergii differences in the gonad. MB: male brain; FB: female brain; T: testis; O: ovary.
4.3 Target gene prediction and bioinformatic analysis
By comparing the expression of miRNAs in various libraries, we were able to obtain DE miRNAs. On the other hand, comparing the DE miRNAs between the female and male brain, 31 mature miRNAs with remarkable up-or down-regulation in the female brain as compared to the male brain, including 19 female brain-biased and 12 male brain-biased miRNAs. Furthermore, comparing the DE miRNAs between ovaries and testis, 59 miRNAs were found in total, of which 29 miRNAs in the ovary were higher than that in testis, while 30 miRNAs expression was lower in the ovary (Figure 1B). The volcanic map shows significant DEMs in the brains and gonads of males and females (Figures 1C,D).
Overall, 116 and 1,255 genes were considered potential targets for expression levels of miRNAs in the brain and gonad, respectively. Most DE miRNAs have multiple target genes, and most of them are regulated by multiple miRNAs. We conducted an enrichment analysis of GO as well as KEGG, to identify the biological role of the DE mRNAs. In KEGG pathway analysis, we set the p-value at 0.05 as the threshold. As compared male brain VS. female brain, several signaling pathways, such as purine metabolism, glycosylphosphatidylinositol, and various types of N-glycan biosynthesis were enrichment (Figure 2A). While in the gonads, Ubiquitin mediated proteolysis and protein processing in the endoplasmic reticulum, Glycine, Serine, and threonine metabolism were also involved in its development (Figure 2B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990677/fgene-13-990677-HTML/image_m/fgene-13-990677-g002.jpg)
FIGURE 2. Top 30 significantly enriched KEGG pathways of putative target genes between the female and male macrobrachium rosenbergii in the brain (A) and gonad (B). Significant enrichment of the KEGG pathway (p < 0.05). The color represents the size of the p-value and the point represents pathway. The redder the color, the smaller the p-value, the larger the number of points, and the more enriched genes. MB: male brain; FB: female brain; T: testis; O: ovary.
At level 2, GO categories are divided into three groups of potential. In the three categories of terms of male brain VS. female brain, the most abundant GO terms are “regulation of transcription by RNA polymerase II”, “nucleus”, “ATP binding”, and “metal ion binding” (Figure 3A). The GO terms such as “positive regulation of transcription”, “DNA-templated nucleus”, “metal ion binding”, and “DNA binding” were rich in the ovary as compared to the testis (Figure 3B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990677/fgene-13-990677-HTML/image_m/fgene-13-990677-g003.jpg)
FIGURE 3. Top 30 GO terms of unique genes targeted by DEMs from the female and male macrobrachium rosenbergii in the brain (A) and gonad (B). GO analysis falls into three categories: biological processes, molecular functions, and cellular components. MB: male brain; FB: female brain; T: testis; O: ovary.
4.4 qRT-PCR verifies DEMs
In order to verify the expression of DE miRNAs in M. rosenbergii, 9 DE miRNAs were randomly selected from the brain and gonad and performed qRT-PCR to verify their distribution in the eye, brain, muscle, liver, gut, ovary, and testis for further investigation (Figures 4,5). The relative expression of U6 single-stranded RNA was measured as the endogenous reference gene, and the relative expression of DE miRNAs screened out was consistent with the sequencing results. In general, DE miRNAs expression was significantly different in the brain and gonad of both genders of M. rosenbergii.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990677/fgene-13-990677-HTML/image_m/fgene-13-990677-g004.jpg)
FIGURE 4. The expression of DEMs in the brain of M. rosenbergii. The endogenous reference gene was U6, and its expression was determined by comparative CT (ΔΔCT) methods. T: testis; O: ovary; M: muscle; L: liver; G: gill; MB: male brain; FB: female brain. *: 0.01 < p < 0.05; **: 0.001 < p < 0.01; ***: p < 0.001.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990677/fgene-13-990677-HTML/image_m/fgene-13-990677-g005.jpg)
FIGURE 5. The expression of DEMs in the gonad of M. rosenbergii. The endogenous reference gene was U6, and its expression was determined by comparative CT (ΔΔCT) methods. M: muscle; B: brain; L: liver; G gill; I: intestines; O: ovary; T: testis. *: 0.01 < p < 0.05; **: 0.001 < p < 0.01; ***: p < 0.001.
4.5 5-HT1 is Regulated by miR-133
Through the above software analysis, we found that miR-133 targets the 5-HT1-3'UTR conservatively (Figure 6A). To further validate the targeting relationship between 5-HT1 and miR-133 in vitro, the pmirGLO -5-HT1-3'UTR and miR-133 mimics were co-transfected into the HEK293 cells. The results confirmed that expression of 5-HT1 was significantly inhibited under the regulation of miR-33, while there were no differences between mutant-type (MT) and negative controls. (Figure 6B). The findings strongly imply that 5-HT1 potentially regulates the development of gonads and gamete generation in M. rosenbergii.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/990677/fgene-13-990677-HTML/image_m/fgene-13-990677-g006.jpg)
FIGURE 6. Double luciferase reporter gene analysis of miR-133 and 5-HT1. (A). The binding site of miR-133 to 5-HT1 3’-UTR. (B). Relative luciferase activity in HEK293 cells. GraphPad Prism 7 draws a histogram. The means ± standard deviations of the three independent experiments are shown in the data. *p < 0.05.
5 Discussion
Even though miRNAs are crucial for controlling gene expression throughout gonadal development and sexual differentiation (Qiu et al., 2018; Sun et al., 2020), Up till now, little research has reported the involvement of miRNAs in crustacean reproductive development (He et al., 2012; Song et al., 2014). Therefore, this research aimed to find DEMs in M. rosenbergii and to assess differential expression patterns in the brain and gonad.
In current research work, 494 miRNAs were detected totally in the brain and gonads of M. rosenbergii, including 271 known miRNAs and 223 novel miRNAs, contributing to the limited database of miRNAs in crustaceans. In general, the abundance of novel miRNAs is lower than that of conserved miRNAs, suggesting that novel miRNAs may play a role or be of lower functional importance in specific gonad stages (Ro et al., 2007). The distribution analysis of miRNA length populations suggested that the size distribution of all four libraries was similar, with a peak at 22 nt, while only the testis library had a peak at 27–29 nt. Peaks in 21–22 nt length represent representative sizes of Dicer derivatives, while peaks at 27–29 nt may be from piRNAs. Usually, piRNAs are about 25–30 nt in length and have high expression in mice, insects, and fish in the gonads (Grimson et al., 2008; Klattenhoff and Theurkauf, 2008). The sperm of M. rosenbergii contained small RNA of similar size (27–29 nt RNA peak), suggesting that piRNA may regulate the development of testis.
Among the 494 miRNAs, we focused on DEMs and identified 49 female-biased miRNAs and 31 male-biased miRNAs. These DEMs might be involved in gender differentiation and gonadal development. For example, miR-7, miR-71, let-7-5p, and miR-184 are necessary for the development of the brain (Diotel et al., 2011). The effect of estrogens, its metabolites, and testosterone have different effects on the development of males’ and females’ brains (Negri-Cesi et al., 2004; He et al., 2012). Previously, researchers found high expression of miR-7 in the pituitary gland of mice and pigs (Li et al., 2011; He et al., 2020). In sows, miR-7 inhibits the synthesis and secretion of Follicle-Stimulating hormones (He et al., 2018). Mice knocked out of MiR-7 present with little gonads production and infertility (Ahmed et al., 2017). miR-71 plays a role in the neurological system and mediates the impact of the germ line on lifespan (Boulias and Horvitz, 2012). Our findings demonstrate that the expressions of miR-7 and miR-71 in the brain of the male in M. rosenbergii were higher than those in the female, suggesting that they might involve in male reproductive development by regulating brain hormones. Previous studies shows, the wide expression of let-7 family in the gonads of insects, fish, and birds (Bizuayehu et al., 2012; Gu et al., 2014; Quah et al., 2015). In our study, let-7-5p, which is a member of the let-7 family, was highly expressed in the testis as compared to the ovaries, suggesting that let-7-5p may participate in the male reproduction of M. rosenbergii. Apart from it, miR-184-3p is expressed highly in the ovary compared to the testis, suggesting that miR-184-3p might participate in the female reproduction of M. rosenbergii. In agreement with our findings, recent research has demonstrated that miR-184 plays a role in the female reproductive system development, and the loss of miR-184 may lead to serious defects in Drosophila eggs (Iovino et al., 2009). At the same time, miR-184 also plays a role in male reproductive development by contributing to the formation of germ cells (Wu et al., 2011).
miR-133b belongs to the miR-133 family, and the research by predecessors has shown that miR-133b is expressed in both mouse and human oocytes and follicles (Dai et al., 2013; Xiao et al., 2014). After gonadotropin-releasing hormone treatment, miR-133b is up-regulated in porcine pituitary cells, suggesting that it may affect gonadotropin release by targeting gonadotropin β and C-Jun (Ye et al., 2013). In tilapia, miR-133b may modulate tagln2 expression to achieve inhibition of primitive follicle formation (Ma et al., 2021). In crustaceans, similarly, miR-133 is differentially expressed in crab oocytes during meiosis and maturation and regulates the 3'-UTR of the crab cyclin B (Song et al., 2014). In this research, the expression of miR-133 in M. rosenbergii was shown to be highly expressed in the brain and ovaries of females, and relatively low or no expression in the brain and testis of males. The analysis of dual-luciferase reporter had shown that dpu-miR-133 may down-regulate the expression of 5-HT1 by predicting targets in the 3'-UTR.
As an important biogenic amine, 5-hydroxytryptamine (5-HT) is involved in the important physiological processes in decapod crustaceans. The serotonin receptor family is classified into 7 subfamilies (5-HT1-5-HT7) (Cortes-Altamirano et al., 2018). 5-HT can induce ovarian maturation in Penaeus vanamei, (Vaca and Alfaro, 2000), Penaeus monodon (Wongprasert et al., 2006), Platycephalus indicus (Tomy et al., 2016), Procambarus clarkii (Kulkarni et al., 1992) and M. rosenbergii (Soonthornsumrith et al,, 2018). 5-HT can promote the secretion of estradiol-17β (E2) and progesterone (P4) in the ovary (Benmansour et al., 2014). We predict that 5-HT may participate in ovarian maturation by localizing specific brain regions and ovaries of shrimp, inducing the secretion of estrogen, which may promote the development of yolk and thus lead to the growth and maturation of oocytes (Soonthornsumrith et al., 2018).
In conclusion, we reported known and novel miRNAs in the brain and gonads of female and male M. rosenbergii. Based on the analysis of miRNAs, 80 DEMs were identified. The analysis of GO and KEGG showed that some DE miRNAs were involved in gonadal development. Through software prediction, miR-133 was found to target 5-HT1, and its targeting relationship was verified by a dual-luciferin report. The data presented in this research can provide the basis for additional information on the genetic mechanism of M. rosenbergii and valuable information for the reproductive control technology of M. rosenbergii.
6 Institutional review board statement
The current research work was carried out in accordance with the Declaration of Helsinki and the Shanghai Ocean University Animal Care and Use Committee granted consent for the study with the approval code SHOU-2021-118.
Data availability statement
The datasets (ANALYZED) for this study can be found in the Figshare https://doi.org/10.6084/m9.figshare.20388261.v1e.
Ethics statement
The animal study was reviewed and approved by Shanghai Ocean University Animal Care and Use Committee with approval number SHOU-2021-118.
Author contributions
Methodology, investigation data curation, analysis, Writing: JX and DL; Visualization: SY and XW; Resources: BL and LG; Review and editing: MJ and HX; Funding, writing—review and editing: ML and WZ; All authors have read and agreed to the published version of the manuscript.
Funding
This research was funded by the Agricultural Science and Technology Independent Innovation Fund project of Jiangsu Province: Technological Innovation and Demonstration of Rice-Loach -Eel Green Co-cultivation in Northern Jiangsu province [CX (20) 2031] and the National Natural Science Foundation of China (31672700).
Conflict of interest
BL was empolyed by the company Huzhou Fengshengwan Aquatic Seed Industry Co. Ltd.
The remaining authors declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.990677/full#supplementary-material
Abbreviations
DEM, differentially expressed miRNA; FB, female brain; GO, gene ontology; KEGG, kyoto encyclopedia of genes and genomes; MB, male brain; qRT-PCR, quantitative real-time PCR; nt, nucleotide.
References
Ahmed, K., LaPierre, M. P., Gasser, E., Denzler, R., Yang, Y., Rulicke, T., et al. (2017). Loss of microRNA-7a2 induces hypogonadotropic hypogonadism and infertility. J. Clin. Investig. 127 (3), 1061–1074. doi:10.1172/JCI90031
Bartel, D. P. (2004). MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 116 (2), 281–297. doi:10.1016/s0092-8674(04)00045-5
Benmansour, S., Privratsky, A. A., Adeniji, O. S., and Frazer, A. (2014). Signaling mechanisms involved in the acute effects of estradiol on 5-HT clearance. Int. J. Neuropsychopharmacol. 17 (5), 765–777. doi:10.1017/S146114571300165X
Betel, D., Wilson, M., Gabow, A., Marks, D. S., and Sander, C. (2008). The microRNA.org resource: targets and expression. Nucleic Acids Res. 36, D149–D153. doi:10.1093/nar/gkm995
Bizuayehu, T. T., Babiak, J., Norberg, B., Fernandes, J. M., Johansen, S. D., and Babiak, I. (2012). Sex-biased miRNA expression in Atlantic halibut (Hippoglossus hippoglossus) brain and gonads. Sex. Dev. 6 (5), 257–266. doi:10.1159/000341378
Bonami, J. R., and Widada, J. S. (2011). Viral diseases of the giant fresh water prawn Macrobrachium rosenbergii: A review. J. Invertebr. Pathol. 106 (1), 131–142. doi:10.1016/j.jip.2010.09.007
Boulias, K., and Horvitz, H. R. (2012). The C. elegans MicroRNA mir-71 acts in neurons to promote germline-mediated longevity through regulation of DAF-16/FOXO. Cell. Metab. 15 (4), 439–450. doi:10.1016/j.cmet.2012.02.014
Brennecke, J., Stark, A., Russell, R. B., and Cohen, S. M. (2005). Principles of microRNA-target recognition. PLoS Biol. 3 (3), e85. doi:10.1371/journal.pbio.0030085
Cao, J. X., Dai, J. Q., Dai, Z. M., Yin, G. L., and Yang, W. J. (2007). A male reproduction-related kazal-type peptidase inhibitor gene in the prawn, Macrobrachium rosenbergii: Molecular characterization and expression patterns. Mar. Biotechnol. 9 (1), 45–55. doi:10.1007/s10126-006-6026-4
Cortes-Altamirano, J. L., Olmos-Hernandez, A., Jaime, H. B., Carrillo-Mora, P., Bandala, C., Reyes-Long, S., et al. (2018). Review: 5-HT1, 5-HT2, 5-HT3 and 5-HT7 receptors and their role in the modulation of pain response in the central nervous system. Curr. Neuropharmacol. 16 (2), 210–221. doi:10.2174/1570159X15666170911121027
Dai, A. Y., Sun, H. X., Fang, T., Zhang, Q., Wu, S. G., Jiang, Y., et al. (2013). MicroRNA-133b stimulates ovarian estradiol synthesis by targeting Foxl2. FEBS Lett. 587 (15), 2474–2482. doi:10.1016/j.febslet.2013.06.023
Davies, W., and Wilkinson, L. S. (2006). It is not all hormones: Alternative explanations for sexual differentiation of the brain. Brain Res. 1126 (1), 36–45. doi:10.1016/j.brainres.2006.09.105
DeFalco, T., and Capel, B. (2009). Gonad morphogenesis in vertebrates: divergent means to a convergent end. Annu. Rev. Cell. Dev. Biol. 25, 457–482. doi:10.1146/annurev.cellbio.042308.13350
Dennis, C. (2004). Brain development: The most important sexual organ. Nature 427 (6973), 390–392. doi:10.1038/427390a
Diotel, N., Do Rego, J. L., Anglade, I., Vaillant, C., Pellegrini, E., Gueguen, M. M., et al. (2011). Activity and expression of steroidogenic enzymes in the brain of adult zebrafish. Eur. J. Neurosci. 34 (1), 45–56. doi:10.1111/j.1460-9568.2011.07731.x
Francis, R. C. (1992). Sexual lability in teleosts Developmental factors. Q. Rev. Biol. 67 (1), 1–18.
Grimson, A., Srivastava, M., Fahey, B., Woodcroft, B. J., Chiang, H. R., King, N., et al. (2008). Early origins and evolution of microRNAs and Piwi-interacting RNAs in animals. Nature 455 (7217), 1193–1197. doi:10.1038/nature07415
Gu, Y., Zhang, L., and Chen, X. (2014). Differential expression analysis of Paralichthys olivaceus microRNAs in adult ovary and testis by deep sequencing. Gen. Comp. Endocrinol. 204, 181–184. doi:10.1016/j.ygcen.2014.05.019
Hayashi, K., Chuva de Sousa Lopes, S. M., Kaneda, M., Tang, F., Hajkova, P., Lao, K., et al. (2008). MicroRNA biogenesis is required for mouse primordial germ cell development and spermatogenesis. PLoS One 3 (3), e1738. doi:10.1371/journal.pone.0001738
He, J., Xu, S., Ji, Z., Sun, Y., Cai, B., Zhang, S., et al. (2020). The role of miR-7 as a potential switch in the mouse hypothalamus-pituitary-ovary axis through regulation of gonadotropins. Mol. Cell. Endocrinol. 518, 110969. doi:10.1016/j.mce.2020.110969
He, J., Zhang, J., Wang, Y., Liu, W., Gou, K., Liu, Z., et al. (2018). MiR-7 mediates the zearalenone signaling pathway regulating FSH synthesis and secretion by targeting FOS in female pigs. Endocrinology 159 (8), 2993–3006. doi:10.1210/en.2018-00097
He, L., Wang, Q., Jin, X., Wang, Y., Chen, L., Liu, L., et al. (2012). Transcriptome profiling of testis during sexual maturation stages in Eriocheir sinensis using Illumina sequencing. PLoS One 7 (3), e33735. doi:10.1371/journal.pone.0033735
He, L., Wang, Y. L., Li, Q., Yang, H. D., Duan, Z. L., and Wang, Q. (2015). Profiling microRNAs in the testis during sexual maturation stages in Eriocheir sinensis. Anim. Reprod. Sci. 162, 52–61. doi:10.1016/j.anireprosci.2015.09.008
Hossain, M. M., Sohel, M. M., Schellander, K., and Tesfaye, D. (2012). Characterization and importance of microRNAs in mammalian gonadal functions. Cell. Tissue Res. 349 (3), 679–690. doi:10.1007/s00441-012-1469-6
Huang, J., Ju, Z., Li, Q., Hou, Q., Wang, C., Li, J., et al. (2011). Solexa sequencing of novel and differentially expressed microRNAs in testicular and ovarian tissues in Holstein cattle. Int. J. Biol. Sci. 7 (7), 1016–1026. doi:10.7150/ijbs.7.1016
Iovino, N., Pane, A., and Gaul, U. (2009). miR-184 has multiple roles in Drosophila female germline development. Dev. Cell. 17 (1), 123–133. doi:10.1016/j.devcel.2009.06.008
Jiang, X. H., and Qiu, G. F. (2013). Female-only sex-linked amplified fragment length polymorphism markers support ZW/ZZ sex determination in the giant freshwater prawn Macrobrachium rosenbergii. Anim. Genet. 44 (6), 782–785. doi:10.1111/age.12067
Jung, H., Yoon, B. H., Kim, W. J., Kim, D. W., Hurwood, D. A., Lyons, R. E., et al. (2016). Optimizing hybrid de Novo transcriptome assembly and extending genomic Resources for giant freshwater prawns (Macrobrachium rosenbergii): The identification of genes and markers associated with reproduction. Int. J. Mol. Sci. 17 (5), E690. doi:10.3390/ijms17050690
Klattenhoff, C., and Theurkauf, W. (2008). Biogenesis and germline functions of piRNAs. Development 135 (1), 3–9. doi:10.1242/dev.006486
Krichevsky, A. M., King, K. S., Donahue, C. P., Khrapko, K., and Kosik, K. S. (2003). A microRNA array reveals extensive regulation of microRNAs during brain development. RNA 9 (10), 1274–1281. doi:10.1261/rna.5980303
Kulkarni, G. K., Nagabhushanam, R., Amaldoss, G., Jaiswal, R. G., and Fingerman, M. (1992). In vivostimulation of ovarian development in the red swamp crayfish, Procambarus clarkii(Girard), by 5-hydroxytryptamine. Invertebr. Reproduction Dev. 21 (3), 231–239. doi:10.1080/07924259.1992.9672242
Landgraf, P., Rusu, M., Sheridan, R., Sewer, A., Iovino, N., Aravin, A., et al. (2007). A mammalian microRNA expression atlas based on small RNA library sequencing. Cell. 129 (7), 1401–1414. doi:10.1016/j.cell.2007.04.040
Lewis, B. P., Shih, I. H., Jones-Rhoades, M. W., Bartel, D. P., and Burge, C. B. (2003). Prediction of mammalian microRNA targets. Cell. 115 (7), 787–798. doi:10.1016/s0092-8674(03)01018-3
Li, M., Liu, Y., Wang, T., Guan, J., Luo, Z., Chen, H., et al. (2011). Repertoire of porcine microRNAs in adult ovary and testis by deep sequencing. Int. J. Biol. Sci. 7 (7), 1045–1055. doi:10.7150/ijbs.7.1045
Ma, Z. S., Yang, J., Zhang, Q. Q., Xu, C. M., Wei, J., Sun, L. N., et al. (2021). miR-133b targets tagln2 and functions in tilapia oogenesis. Comp. Biochem. Physiol. B Biochem. Mol. Biol. 256, 110637. doi:10.1016/j.cbpb.2021.110637
Macias, S., Michlewski, G., and Caceres, J. F. (2009). Hormonal regulation of microRNA biogenesis. Mol. Cell. 36 (2), 172–173. doi:10.1016/j.molcel.2009.10.006
Meng, X. L., Liu, P., Jia, F. L., Li, J., and Gao, B. Q. (2015). De novo transcriptome analysis of Portunus trituberculatus ovary and testis by RNA-seq: Identification of genes involved in gonadal development. PLoS One 10 (6), e0128659. doi:10.1371/journal.pone.0128659
Meng, X., Zhang, X., Li, J., and Liu, P. (2018). Identification and comparative profiling of ovarian and testicular microRNAs in the swimming crab Portunus trituberculatus. Gene 640, 6–13. doi:10.1016/j.gene.2017.10.026
Miranda, K. C., Huynh, T., Tay, Y., Ang, Y. S., Tam, W. L., Thomson, A. M., et al. (2006). A pattern-based method for the identification of MicroRNA binding sites and their corresponding heteroduplexes. Cell. 126 (6), 1203–1217. doi:10.1016/j.cell.2006.07.031
Murchison, E. P., Stein, P., Xuan, Z., Pan, H., Zhang, M. Q., Schultz, R. M., et al. (2007). Critical roles for Dicer in the female germline. Genes. Dev. 21 (6), 682–693. doi:10.1101/gad.1521307
Negri-Cesi, P., Colciago, A., Celotti, F., and Motta, M. (2004). Sexual differentiation of the brain: role of testosterone and its active metabolites. J. Endocrinol. Investig. 27, 120–127.
O'Brien, J., Hayder, H., Zayed, Y., and Peng, C. (2018). Overview of MicroRNA biogenesis, mechanisms of actions, and circulation. Front. Endocrinol. 9, 402. doi:10.3389/fendo.2018.00402
Peng, J., Wei, P., Zhang, B., Zhao, Y., Zeng, D., Chen, X., et al. (2015). Gonadal transcriptomic analysis and differentially expressed genes in the testis and ovary of the Pacific white shrimp (Litopenaeus vannamei). BMC Genomics 16, 1006. doi:10.1186/s12864-015-2219-4
Phoungpetchara, I., Tinikul, Y., Poljaroen, J., Changklungmoa, N., Siangcham, T., Sroyraya, M., et al. (2012). Expression of the male reproduction-related gene (Mar-Mrr) in the spermatic duct of the giant freshwater prawn, Macrobrachium rosenbergii. Cell. Tissue Res. 348 (3), 609–623. doi:10.1007/s00441-012-1380-1
Piprek, R. P. (2010). Molecular and cellular machinery of gonadal differentiation in mammals. Int. J. Dev. Biol. 54 (5), 779–786. doi:10.1387/ijdb.092939rp
Qiao, H., Fu, H., Jin, S., Wu, Y., Jiang, S., Gong, Y., et al. (2012). Constructing and random sequencing analysis of normalized cDNA library of testis tissue from oriental river prawn (Macrobrachium nipponense). Comp. Biochem. Physiol. Part D. Genomics Proteomics 7 (3), 268–276. doi:10.1016/j.cbd.2012.04.003
Qiu, W., Zhu, Y., Wu, Y., Yuan, C., Chen, K., and Li, M. (2018). Identification and expression analysis of microRNAs in medaka gonads. Gene 646, 210–216. doi:10.1016/j.gene.2017.12.062
Quah, S., Breuker, C. J., and Holland, P. W. (2015). A diversity of conserved and novel ovarian MicroRNAs in the speckled wood (Pararge aegeria). PLoS One 10 (11), e0142243. doi:10.1371/journal.pone.0142243
Ro, S., Song, R., Park, C., Zheng, H., Sanders, K. M., and Yan, W. (2007). Cloning and expression profiling of small RNAs expressed in the mouse ovary. RNA 13 (12), 2366–2380. doi:10.1261/rna.754207
Soonthornsumrith, B., Saetan, J., Kruangkum, T., Thongbuakaew, T., Senarai, T., Palasoon, R., et al. (2018). Three-dimensional organization of the brain and distribution of serotonin in the brain and ovary, and its effects on ovarian steroidogenesis in the giant freshwater prawn, Macrobrachium rosenbergii. Invert Neurosci. 18 (2), 5. doi:10.1007/s10158-018-0209-3
Sharabi, O., Manor, R., Weil, S., Aflalo, E. D., Lezer, Y., Levy, T., et al. (2016). Identification and characterization of an insulin-like receptor involved in Crustacean reproduction. Endocrinology 157 (2), 928–941. doi:10.1210/en.2015-1391
Song, Y. N., Shi, L. L., Liu, Z. Q., and Qiu, G. F. (2014). Global analysis of the ovarian microRNA transcriptome: Implication for miR-2 and miR-133 regulation of oocyte meiosis in the Chinese mitten crab, Eriocheir sinensis (Crustacea:Decapoda). Bmc Genomics 15, 547. doi:10.1186/1471-2164-15-547
Soonthornsumrith, B., Saetan, J., Kruangkum, T., Thongbuakaew, T., Senarai, T., Palasoon, R., et al. (2018). Three-dimensional organization of the brain and distribution of serotonin in the brain and ovary, and its effects on ovarian steroidogenesis in the giant freshwater prawn, Macrobrachium rosenbergii. Invert. Neurosci. 18 (2), 5. doi:10.1007/s10158-018-0209-3
Staelens, J., Rombaut, D., Vercauteren, I., Argue, B., Benzie, J., and Vuylsteke, M. (2008). High-density linkage maps and sex-linked markers for the black tiger shrimp (Penaeus monodon). Genetics 179 (2), 917–925. doi:10.1534/genetics.107.080150
Sun, L. L., Zhong, Y., Qiu, W. W., Guo, J., Gui, L., and Li, M. Y. (2020). MiR-26 regulates ddx3x expression in medaka (Oryzias latipes) gonads. Comp. Biochem. Physiology B-Biochemistry Mol. Biol. 246–247, 110456. doi:10.1016/j.cbpb.2020.110456
Tomy, S., Saikrithi, P., James, N., Balasubramanian, C. P., Panigrahi, A., Otta, S. K., et al. (2016). Serotonin induced changes in the expression of ovarian gene network in the Indian white shrimp, Penaeus indicus. Aquaculture 452, 239–246. doi:10.1016/j.aquaculture.2015.11.003
Uawisetwathana, U., Leelatanawit, R., Klanchui, A., Prommoon, J., Klinbunga, S., and Karoonuthaisiri, N. (2011). Insights into eyestalk ablation mechanism to induce ovarian maturation in the black tiger shrimp. PLoS One 6 (9), e24427. doi:10.1371/journal.pone.0024427
Vaca, A. A., and Alfaro, J. (2000). Ovarian maturation and spawning in the white shrimp, Penaeus vannamei, by serotonin injection. Aquaculture 182 (3-4), 373–385. doi:10.1016/s0044-8486(99)00267-7
Ventura, T., Manor, R., Aflalo, E. D., Weil, S., Raviv, S., Glazer, L., et al. (2009). Temporal silencing of an androgenic gland-specific insulin-like gene affecting phenotypical gender differences and spermatogenesis. Endocrinology 150 (3), 1278–1286. doi:10.1210/en.2008-0906
Vogt, G. (2018). Investigating the genetic and epigenetic basis of big biological questions with the parthenogenetic marbled crayfish: A review and perspectives. J. Biosci. 43 (1), 189–223. doi:10.1007/s12038-018-9741-x
Wongprasert, K., Asuvapongpatana, S., Poltana, P., Tiensuwan, M., and Withyachumnarnkul, B. (2006). Serotonin stimulates ovarian maturation and spawning in the black tiger shrimp Penaeus monodon. Aquaculture 261 (4), 1447–1454. doi:10.1016/j.aquaculture.2006.08.044
Wu, J., Bao, J., Wang, L., Hu, Y., and Xu, C. (2011). MicroRNA-184 downregulates nuclear receptor corepressor 2 in mouse spermatogenesis. BMC Dev. Biol. 11, 64. doi:10.1186/1471-213X-11-64
Xiao, G. H., Xia, C. L., Yang, J., Liu, J. Q., Du, H. Z., Kang, X. J., et al. (2014). MiR-133b regulates the expression of the actin protein TAGLN2 during oocyte growth and maturation: A potential target for infertility therapy. Plos One 9 (6), e100751. doi:10.1371/journal.pone.0100751
Ye, R. S., Xi, Q. Y., Qi, Q., Cheng, X., Chen, T., Li, H., et al. (2013). Differentially expressed miRNAs after GnRH treatment and their potential roles in FSH regulation in porcine anterior pituitary cell. PLoS One 8 (2), e57156. doi:10.1371/journal.pone.0057156
Yue, Y., Meng, Y. J., Ma, H. R., Hou, S. Q., Cao, G. Z., Hong, W. L., et al. (2016). A large family of Dscam genes with tandemly arrayed 5 ' cassettes in Chelicerata. Nat. Commun. 7, 11252. doi:10.1038/ncomms11252
Zhu, L., Gao, J., Huang, K., Luo, Y., Zhang, B., and Xu, W. (2015). miR-34a screened by miRNA profiling negatively regulates Wnt/β-catenin signaling pathway in Aflatoxin B1 induced hepatotoxicity. Sci. Rep. 5, 16732. doi:10.1038/srep16732
Zou, Y., Wu, Z., Fan, Z., Liang, D., Wang, L., Song, Z., et al. (2020). Analyses of mRNA-seq and miRNA-seq of the brain reveal the sex differences of gene expression and regulation before and during gonadal differentiation in 17beta-estradiol or 17alpha-methyltestosterone-induced olive flounder (Paralichthys olivaceus). Mol. Reprod. Dev. 87 (1), 78–90. doi:10.1002/mrd.23303
Keywords: macrobrachium rosenbergii, brain, gonad, RNA-seq, miRNA, reproduction, MiR-133, 5-HT1
Citation: Xia J, Liu D, Zhou W, Yi S, Wang X, Li B, Jawad M, Xu H, Gui L and Li M (2022) Comparative transcriptome analysis of brain and gonad reveals reproduction-related miRNAs in the giant prawn, Macrobrachium rosenbergii. Front. Genet. 13:990677. doi: 10.3389/fgene.2022.990677
Received: 10 July 2022; Accepted: 04 August 2022;
Published: 26 August 2022.
Edited by:
Feng He, Ocean University of China, ChinaReviewed by:
Hongyan Xu, College of Fisheries, Southwest University, Chongqing, ChinaQiaomu Hu, Yangtze River Fisheries Research Institute (CAFS), China
Copyright © 2022 Xia, Liu, Zhou, Yi, Wang, Li, Jawad, Xu, Gui and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wenzong Zhou, zhouwz001@163.com; Mingyou Li, myli@shou.edu.cn
†These authors have contributed equally to this work and share the first authorship