Next Article in Journal
Construction of a Matrix Cancer-Associated Fibroblast Signature Gene-Based Risk Prognostic Signature for Directing Immunotherapy in Patients with Breast Cancer Using Single-Cell Analysis and Machine Learning
Next Article in Special Issue
Insulin-like Androgenic Gland Hormone Induced Sex Reversal and Molecular Pathways in Macrobrachium nipponense: Insights into Reproduction, Growth, and Sex Differentiation
Previous Article in Journal
Psoralen Alleviates Renal Fibrosis by Attenuating Inflammasome-Dependent NLRP3 Activation and Epithelial–Mesenchymal Transition in a Mouse Unilateral Ureteral Obstruction Model
Previous Article in Special Issue
The Trophoblast Compartment Helps Maintain Embryonic Pluripotency and Delays Differentiation towards Cardiomyocytes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Male Sex-Related Genes Regulated by SDHB in Macrobrachium nipponense Based on Transcriptome Analysis after an RNAi Knockdown

1
Wuxi Fisheries College, Nanjing Agricultural University, Wuxi 214081, China
2
Key Laboratory of Freshwater Fisheries and Germplasm Resources Utilization, Ministry of Agriculture and Rural Affairs, Freshwater Fisheries Research Center, Chinese Academy of Fishery Sciences, Wuxi 214081, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(17), 13176; https://doi.org/10.3390/ijms241713176
Submission received: 31 July 2023 / Revised: 18 August 2023 / Accepted: 22 August 2023 / Published: 24 August 2023

Abstract

:
The oriental river prawn (Macrobrachium nipponense) is a commercially important species in Asia. A previous study showed that the succinate dehydrogenase complex iron sulfur subunit B (SDHB) gene participates in testes development in this species through its effect on the expression of the insulin-like androgenic gland hormone gene. This study knocked-down the Mn-SDHB genes in M. nipponense using RNAi. A transcriptome analysis of the androgenic gland and testes was then performed to discover the male sex-related genes regulated by SDHB and investigate the mechanism of male sexual development in this species. More than 16,623 unigenes were discovered in each sample generated. In the androgenic gland, most of the differentially expressed genes were enriched in the hypertrophic cardiomyopathy pathway, while in the testes, they were enriched in the citrate cycle pathway. In addition, after Mn-SDHB knockdown, five genes were found to be downregulated in the androgenic gland in a series of biological processes associated with phosphorylated carbohydrate synthesis and transformations in the glycolysis/gluconeogenesis pathway. Moreover, a total of nine male sex-related genes were identified including Pro-resilin, insulin-like androgenic gland hormone, Protein mono-ADP-ribosyltransferase PAPR11, DNAJC2, C-type Lectin-1, Tyrosine-protein kinase Yes, Vigilin, and Sperm motility kinase Y-like, demonstrating the regulatory effects of Mn-SDHB, and providing a reference for the further study of the mechanisms of male development in M. nipponense.

1. Introduction

The oriental river prawn (Macrobrachium nipponense) is an economically important species in China, and is widely distributed in freshwater and low-salinity estuarine regions of China and other Asian countries [1,2,3]. The annual production of M. nipponense reached 225,321 tons in 2019 with a value of over 3 billion dollars, producing huge economic benefits [4]. A previous study showed that both the testes and ovaries of newly hatched M. nipponense can reach sexual maturity within 45 days after hatching [5]. This enables inbreeding between newborn shrimps, resulting in shrimps with a small market size and short life span, and therefore restricts the sustainable development of the M. nipponense industry. There is therefore an urgent need to understand the reproductive mechanisms of M. nipponense, in order to establish artificial techniques to regulate their gonad development processes.
Succinate dehydrogenase complex iron sulfur subunit B (SDHB) is a crucial protein subunit comprising one of the four components of succinate dehydrogenase, the only enzyme to be involved in the tricarboxylic acid cycle through its association with the inner mitochondrial membrane, where it catalyzes the oxidation of succinate [6,7]. Its particular iron–sulfur structure enables the succinate dehydrogenase complex to bind to two small integral membrane proteins of 13.5 and 15.5 kDa, which participate in the oxidative stress of oxidative phosphorylation, and play vital roles in the electron transport chain. The spermatogenic cells of Lymnaea stagnalis contain succinate dehydrogenase, which prevents metabolic acidosis in the testes through the anaerobic production of lactate and succinate by the Sertoli cells. They also play essential roles in developing rat testes [8]. SDHB also had multiple functions, including resistance to carboxin in Ustilago maydis [9], the prevention of superoxide generation and premature aging in Caenorhabditis elegans [10], and resistance to the fungicide boscalid in Sclerotinia sclerotiorum [11]. A previous study showed that knockdown of the expression of SDHB by RNAi in male M. nipponense can inhibit the development of the testes by affecting the expression of the insulin-like androgenic gland hormone (IAG) [12].
Many researchers have proposed that the eyestalk-androgenic gland-testis endocrine axis in male crustaceans plays an essential role in regulating male sexual differentiation and reproduction [13,14]. The androgenic gland and testes have therefore become the target tissues for studies of male differentiation and reproduction in crustacean species. The androgenic gland is a special organ in male crustaceans. The hormones it secretes play essential roles in testis formation and male secondary sexual characteristics [15,16,17]. IAG is specifically expressed in the androgenic gland, and has been shown to regulate the process of male sexual differentiation and reproduction in many crustaceans [18,19,20,21,22]. Sexual development and the regulatory mechanism of IAG in crustaceans are very complicated processes and rely on the expression of multiple genes. Recent studies have shown that several sex-related genes, including CFSH, Vtg, Wnt4, the Dmrt gene family, the Sox gene family, the cell cycle gene family, and Fem-1, participate in the regulation of IAG, either directly or indirectly. Nevertheless, further investigations are required to explore the interrelationships between sex-related genes and their association with IAG [23,24,25,26]. Knockdown of IAG using RNAi had a significant inhibitory effect on spermatogenesis in M. rosenbergii [27]. The testes play an essential regulatory role in reproduction, sexual maturity, and sex differentiation in M. nipponense. Some transcriptome analyses have been conducted on the testes of this species [28,29], and a number of genes selected from the testes have been shown to regulate its male reproductive processes [30,31,32].
In this study, we aimed to identify the genes regulated by SDHB in M. nipponense using transcriptome profiling analysis after SDHB knockdown using RNAi. The differentially expressed genes (DEGs) identified in this study may be involved in regulating male reproduction, providing valuable information regarding the mechanisms of male reproduction in crustacean species in general.

2. Results

2.1. Overview of the RNA-Seq of M. nipponense

Androgenic gland and testis samples were collected from both the control and the RNAi groups on day 7 after the injection of double-stranded RNA to form RNA-seq samples, namely SDHB_AG for the androgenic gland SDHB knockdown group; SDHB_T for the testis SDHB knockdown group; C_AG for the androgenic gland control group; and C_T for the testis control group. After quality control screening, about 84.34 Gb of the total raw reads (about 7.02 Gb for each sample) were generated. As the data show, the percentage of Q30 bases of all products exceeded 92.42%, indicating high quality sequences. Through alignment with the reference genome, the mapping rate exceeded 91.23% and the unigenes obtained exceeded 16,623 for each RNA-seq sample (Table 1).
All sequence reads were deposited in the Sequence Read Archive of the National Center for Biotechnology Information (NCBI) (accession SRR24335442-SRR24335459) under Bioproject PRJNA962267.

2.2. Identification and Functional Analysis of DEGs

Fragments per kilobase per million reads (FPKM) is a common measure for estimating gene expression levels in transcriptome sequencing data analysis. For the 12 samples in this study, the mean FPKM values for each gene ranged from 29.76 to 34.43, and the overall trends in gene expression were similar between the samples (Table S1, Figure S1).
The DEG analysis results revealed a total of 67 DEGs in the “SDHB_T vs. C_T” comparison. Of these, 49 genes were upregulated and 18 were downregulated. Similarly, the “SDHB_AG vs. C_AG” comparison identified 235 DEGs, consisting of 52 upregulated genes and 183 downregulated genes.

2.3. GO and KEGG Enrichment Analysis of DEGs

The gene products of the total unigenes and DEGs were clustered according to the GO database to describe their functional attributes. All of the unigenes and DEGs were clustered into cellular components, molecular functions, and biological processes (Figure 1). Most DEGs in the “SDHB_AG vs. C_AG” comparison involved terms including “cytoplasm” (32 DEGs), “extracellular region” (27 DEGs), “metal ion binding” (22 DEGs), “integral component of membrane” (22 DEGs), and “ATP binding” (18 DEGs). In the “SDHB_T vs. C_T” comparison, most of the DEGs involved terms including “cytoplasm” (12 DEGs), “ATP binding” (10 DEGs), “cytosol” (8 DEGs), and “protein serine/threonine kinase activity” (7 DEGs).
In the KEGG pathway enrichment analysis, DEGs in the “SDHB_AG vs. C_AG” comparison were enriched including those related to hypertrophic cardiomyopathy (HCM), dilated cardiomyopathy (DCM), and the PI3K-Akt signaling pathway. DEGs in the “SDHB_T vs. C_T” comparison were enriched including those related to the citrate cycle (TCA cycle), oxidative phosphorylation, and the MAPK signaling pathway. The top 20 pathways with the largest number of DEGs are shown in Figure 2.
Notably, some DEGs were enriched in energy-related biological pathways, including the TCA cycle and the glycolysis/gluconeogenesis pathway (Figure 2). In the TCA cycle pathway, the SDHB gene was downregulated in both the “SDHB_AG vs. C_AG” comparison and the “SDHB_T vs. C_T” comparison, as a direct result of the RNAi treatment. Additionally, in the “SDHB_AG vs. C_AG” comparison, five genes were found to be downregulated in a set of relatively continuous biological processes related to the phosphorylated carbohydrate synthesis and transformation in glycolysis/gluconeogenesis pathways, including hexokinase, fructose-1,6-bisphosphatase I (FBP), triosephosphate isomerase (TPI), phosphoenolpyruvate carboxykinase-GTP (PCK-GTP), and multiple inositol-polyphosphate phosphatase (MINPP1) (Figure 3).

2.4. Male Sexual Development-Related DEGs

In the “SDHB_AG vs. C_AG” and “SDHB_T vs. C_T” comparisons, a total of nine male development-related genes were screened (Table 2). Among them, Pro-resilin is a DEG expressed during the embryonic development of male and female crustaceans, while IAG and the serine protease inhibitor 1 (SERPINB1) are related to hormone regulating processes, and the other six genes are related to spermiogenesis and sperm motility (Table 2). As shown in Table 2, after the knockdown of SDHB in the androgenic gland of M. nipponense, Pro-resilin and IAG were downregulated, while Protein mono-ADP-ribosyltransferase PAPR11 (PARP11) was upregulated. In the testes of M. nipponense, DNAJC2 and C-type Lectin-1 (CTL-1) were upregulated, while Pro-resilin, Tyrosine-protein kinase Yes (YES), Vigilin, and sperm motility kinase Y-like (SMOK-Y-like) were upregulated.

2.5. Validation of DEGs by qRT-PCR

To validate the transcriptome results, six DEGs from each comparison that showed significantly different expression levels in the qRT-PCR analysis were selected at random. Positive numbers represent an upward trend, and negative numbers represent a downward trend. As shown in Figure 4, the expression patterns of the DEGs identified using qRT-PCR were generally similar to those obtained in the RNA-Seq analyses, although the relative expression levels were not completely consistent, thus proving that the transcriptome sequencing data were reliable.

3. Discussion

In this study, we attempted to identify male development-related genes and pathways regulated by SDHB in M. nipponense, based on the transcriptomes observed after Mn-SDHB knockdown using RNAi. This method can inhibit gene expression or translation through the short double-stranded RNA molecules in a cell’s cytoplasm [33], and has been widely used in gene function analysis in M. nipponense. Based on the RNAi treatment of Mn-SDHB dsRNA, a previous study found the inhibition of testis development and decrement in sperm in M. nipponense, indicating that SDHB may be involved in male sexual development in this species [12]. Because Mn-SDHB dsRNA was similarly employed in this study, based on both qRT-PCR and the study of the differential expression of the transcriptome, decreasing expression of Mn-SDHB was detected in the androgenic gland and testis in samples from day 7 following the Mn-SDHB dsRNA treatment, indicating the effectiveness of the dsRNA treatment. The qRT-PCR analysis of some representative DEGs suggested that the transcriptome sequencing and differential expression analysis were reliable.
In the androgenic gland of M. nipponense, five downregulated genes were detected in the glycolysis/gluconeogenesis pathway (which promotes the conversion of glucose into pyruvate, which releases free energy to form ATP, and pyruvate, which fuels the tricarboxylic acid cycle), as well as precursors for secondary metabolism, and amino acid and fatty acid biosynthesis. These DEGs are involved in a relatively continuous series of biological processes (Figure 3), indicating that the Mn-SDHB gene promotes the synthesis of phosphorylated carbohydrate and the transformations in these processes, and has an active role in the energy-related functions mentioned above. In addition, the glycolysis/gluconeogenesis pathway is involved in most biological processes, and has been shown to be an important pathway participating in male development in M. nipponense. In previous transcriptome studies in M. nipponense, the expression of many genes in the glycolysis/gluconeogenesis pathway was found to vary during the post-larval developmental stages 5 (PL5) to PL25, the sensitive periods for gonad development [34]. In addition, after ablation of the eyestalk, which is a key organ for secreting many hormones during male sexual development in crustaceans, the glycolysis/gluconeogenesis pathway is also one of the main sources of enrichment of the DEGs in the testes and androgenic gland of M. nipponense [35]. Moreover, the downregulated genes have also been shown to be associated with male sexual development in mammals. Phosphoenolpyruvate carboxykinase in the Leydig cells of prepubertal mouse testes have been shown to play an important role in steroidogenesis [36]. In a study of human disease, expression of the androgen receptor gene is abnormal in triosephosphate isomerase deficient patients, indicating that triosephosphate isomerase may affect the action of androgen [37]. There is also evidence that hexokinase, TPI, and MINPP1 are regulated by androgen [38,39,40]. Although the genes in the glycolysis/gluconeogenesis pathway have also been shown to play an important role in processes such as spermatogenesis or increasing sperm motility [41,42], this study detected no differences in the expression of these genes in the testes of M. nipponense, which might suggest that the different expression of Mn-SDHB does not affect these biological processes by regulating the glycolysis/gluconeogenesis pathway. According to previous in situ hybridization studies of the androgenic gland of M. nipponense, SDHB is expressed in the ejaculatory bulb surrounding the androgenic gland cells rather than in the androgenic gland cells themselves [43], indicating that SDHB might promote and support the formation of androgenic gland cells by influencing the glycolysis/gluconeogenesis process, and that it plays an essential role in maintaining the normal structure and function of the androgenic gland.
In this study, the differential expression of nine male sex-related genes was detected after Mn-SDHB knockdown. Among them is one gene that is differentially expressed during the embryonic development of male and female crustaceans. Pro-resilin is a type of resilin in a family of elastic proteins that includes elastin, as well as gluten, gliadin, abductin, and spider silks. Pro-resilin genes were found to be upregulated in males compared to females in the embryonic development of M. rosenbergii [44], a related species in the same genus as M. nipponense, indicating that it is important in the development of male prawn embryos. In this study, after Mn-SDHB knockdown, the pro-resilin gene was down-regulated in the androgenic gland and up-regulated in the testes of M. nipponense, showing that the effects of Mn-SDHB on pro-resilin vary in different male development-related tissues.
Two genes were related to the hormone regulatory process, IAG and SERPINB1. Coded by the IAG gene, insulin-like androgenic gland hormone is an insulin analogue of the insulin/insulin-like growth factor family, and has been shown to play essential roles in male sexual differentiation and development in crustacean species [27,45,46,47]. In the androgenic gland of M. nipponense, the expression of Mn-IAG was found to decrease after treatment with Mn-SDHB dsRNA, which is consistent with the results of previous qRT-PCR studies, and again points to the positive regulatory effect of Mn-SDHB on Mn-IAG in M. nipponense. SERPINB1 is an inhibitor of serine proteases, and resides mainly in the cytoplasm of neutrophils and monocytes, suppressing the enzymatic activities of serine proteases and preventing unwanted cellular damage during degranulation [48]. In rats, the SERPINB1 gene was apparently down-regulated after a Sertoli cell-selective knockout of the androgen receptor, suggesting that it participates in tubular restructuring and cell junction dynamics processes in Sertoli cells, controlled in part by androgens [49]. Mn-SERPINB1 was also downregulated in the androgenic gland of M. nipponense after Mn-SDHB knockdown.
Six genes were related to spermiogenesis and sperm motility. PAPR11 is a kind of ADP-ribosyltransferase poly polymerase, which plays a role in immune processes such as antiviral activity [50]. Deletion of the PARP11 gene affects the shape of the cell nucleus and causes the formation of abnormally shaped fertilization-incompetent sperm, resulting in teratozoospermia and male infertility in mice. Therefore, it was predicted to have functional relevance for nuclear envelope stability and nuclear reorganization during spermiogenesis [51]. C-type lectins characteristically require calcium-related proteins to mediate essential cell functions through binding to carbohydrates. In Ancylostoma ceylanicum, the CTL-1 protein is a male gender-specific C-type lectin identified from sperm and soluble protein extracts and is important in hookworm reproductive physiology [52]. In studies on humans and mice, the testes/sperm galactosyl receptor was identified as a C-type lectin with possible roles in cell–cell interaction during spermatogenesis and sperm-zona pellucida binding at fertilization [53]. DNAJ chaperone was initially identified as a heat shock protein in Escherichia coli, and a prokaryotic homologue of the eukaryotic heat shock protein 70. It assists other proteins in their folding, transport, and assembly into complexes, which are necessary for testis development and spermatogenesis [54]. DNAJC2 was described in bull sperm by Selvaraju in 2017 and is involved in stem cell differentiation and early embryonic development [55]. Tyrosine-protein kinase Yes Protein was shown to participate in the phosphorylation of tyrosine residues, whose truncated forms are related to spermiogenesis [55]. Vigilin, also known as high-density lipoprotein binding protein, is related to the processes of DNA damage repair [56], chromatin condensation, and gene silencing and according to a transcriptome study, its expression is the greatest in ovine sperm [55]. The SMOK-Y-like gene belongs to the sperm motility kinase family and is involved in sperm motility [57]. These male development-related genes were differentially expressed after the decreased expression of Mn-SDHB. However, only the DNAJC2 and CLT-1 found in the testes of M. nipponense were downregulated following Mn-SDHB treatment, indicating that the influence of Mn-SDHB on male development is complicated and requires further study.

4. Materials and Methods

4.1. Knockdown of the Expression of SDHB by RNAi in M. nipponense

The expression of SDHB of M. nipponense was knocked-down using RNAi. The specific RNAi primer for SDHB with a T7 promoter site was described in Table S2 in our previous study [12]. A Transcript Aid™ T7 High Yield Transcription kit (Fermentas Inc., Horsham, PA, USA) was used to synthesize Mn-SDHB dsRNA (dsSDHB) and GFP dsRNA (dsGFP), according to the manufacturer’s protocol. The dsGFP was used as a negative control [58]. A total of 300 mature male M. nipponense with body weights ranging from 3.26 to 4.76 g were collected. They were randomly divided into a dsSDHB treatment group (RNAi) and a dsGFP treatment group (control), and kept in separate 500 L tanks for 72 h prior to injection with either dsSDHB or dsGFP. The injected doses were 4 μg/g in both cases, consistent with the dose rates used in previous studies [59,60]. Androgenic gland and testis samples were collected from both the control and the RNAi groups on day 7 after dsSDHB and dsGFP injection, consistent with previous studies. qRT-PCR analysis revealed that the expression of Mn-SDHB decreased by over 90% in both the testes and androgenic glands of the dsSDHB-treated group, compared with the dsGFP group.

4.2. RNA Isolation, Library Construction, and Sequencing

The testes and androgenic glands of five M. nipponense from the dsSDHB-treated group (RNAi) and the dsGFP-treated group (control), were collected and individually pooled together to form a biological replicate, and three replicates were used for the transcriptome analysis.
Total RNA was extracted from each replicate using TRIzol reagent (Invitrogen, Waltham, MA, USA) according to the manufacturer’s protocol, followed by RNA purity and quantification evaluation using a Nano Drop 2000 spectrophotometer (Thermo Scientific, Waltham, MA, USA). RNA integrity was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). With the valid RNA samples (RIN ≥ 7, 28S:18S ≥ 0.7, >1 μg), libraries were then constructed using a VAHTS Universal V6 RNA-seq Library Prep Kit (Vazyme, Nanjing, China) according to the manufacturer’s instructions. The libraries were sequenced using an Illumina Novaseq 6000 platform (Illumina, San Diego, CA, USA) to generate the 150 bp paired-end reads.

4.3. Differential Gene Expression Analysis

Low-quality raw reads for each sample were removed from the data using fastp [61] with the default parameters. The clean reads obtained were mapped to the M. nipponense reference genome (Genbank access numbers: GCA_015110555.1 and GCA_015104395.1) using HISAT2 [62]. Gene expression was calculated using the FPKM method, where FPKM = cDNA fragments/mapped fragments (millions)/transcript length (kb), using HTSeq-count [63].
Differential expression analysis was performed using DESeq2 [64]. The false discovery rate (FDR) was calculated using the Benjamini–Hochberg correction method [65] with q-value < 0.05 and fold change >2 or <0.5 as the thresholds for significant DEGs.
Based on the hypergeometric distribution, with a threshold p-value < 0.05, the functional annotation and classification of DEGs were conducted according to the GO [66] and KEGG databases [67] to analyze the metabolic pathways of the enriched DEGs.

4.4. Quantitative Analysis

qRT-PCR was performed to evaluate the sequencing and data analysis and to validate the DEGs. Gonadal RNA was extracted (100 mg) using 1 mL TRIzol reagent (TaKaRa, Tokyo, Japan) and first-strand cDNA synthesis was performed using a Reverse Transcriptase M-MLV Kit (TaKaRa, Tokyo, Japan).
Genes were randomly selected from the DEGs for evaluation and their functional annotations, and the primers used are shown in Table S2. The qRT-PCR was performed using an iCycler iQ5 real-time PCR system (Bio-Rad, Hercules, CA, USA), with eukaryotic translation initiation factor 5 A (EIF) as the reference gene, because of its stable expression in a variety of situations [68]. The reaction was amplified with 35 cycles of 94 °C for 30 s, 50 °C for 30 s, and 72 °C for 1 min, followed by 10 min incubation at 72 °C as the extension step. Each sample had three replicates while each reaction had three controls: nuclease-free water; primer-free water; and template-free water. The system recorded fluorescence curves and data automatically, and dissociation curves of the amplified products were recorded at the end of each PCR. The mRNA expression levels were determined using the 2−ΔΔCT method [69].

4.5. Ethics Statement

The protocols of all experiments involving M. nipponense were approved by the Institutional Animal Care and Use Ethics Committee of the Freshwater Fisheries Research Center of the Chinese Academy of Fishery Sciences (Wuxi, China).

5. Conclusions

In this study, we knocked-down Mn-SDHB expression using RNAi and performed transcriptome studies of the androgenic gland and testes of M. nipponense. In the androgenic gland group treated with dsSDHB, five downregulated genes participating in phosphorylated carbohydrate synthesis and transformation were found in the glycolysis/gluconeogenesis pathway, which is important during the sexual development of male M. nipponense. Moreover, a total of nine male development-related genes were screened from the DEGs identified from the androgenic glands and testes of M. nipponense, demonstrating some of the regulatory effects of Mn-SDHB in the male development process of M. nipponense, and providing reference data for the further study of the mechanisms of male development in crustaceans in general.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms241713176/s1.

Author Contributions

Conceptualization, S.J. (Shubo Jin); Formal analysis, Z.G.; Investigation, Z.G.; Resources, W.Z. and S.J. (Sufei Jiang); Writing—original draft, Z.G.; Writing—review & editing, H.Y. and P.C.; Supervision, H.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by grants from the Natural Science Foundation of Jiangsu Province (BK20221207), the Central Public-Interest Scientific Institution Basal Research Fund CAFS (2020TD36), the Seed Industry Revitalization Project of Jiangsu Province (JBGS (2021) 118), the Jiangsu Agricultural Industry Technology System earmarked fund for CARS-48, and the New Cultivar Breeding Major Project of Jiangsu Province (PZCZ201745).

Institutional Review Board Statement

All experimental protocols and methods were approved in October 2019 (Authorization No. 20191015006) by the Animal Care and Use Ethics Committee in the Freshwater Fisheries Research Center (Wuxi, China).

Data Availability Statement

The data presented in this study are available on request from the corresponding author for scientific purposes.

Acknowledgments

The authors are thankful for the Jiangsu Province Platform for the Conservation and Utilization of Agricultural Germplasm.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cai, Y.; Shokita, S. Report on a collection of freshwater shrimps (Crustacea: Decapoda: Caridea) from the Philippines, with descriptions of four new species. Raffles Bull. Zool. 2006, 54, 245–270. [Google Scholar]
  2. Salman, S.D.; Page, T.J.; Naser, M.D.; Yasser, A.A.G. The invasion of Macrobrachium nipponense (De Haan, 1849)(Caridea: Palaemonidae) into the southern Iraqi marshes. Aquat. Invasions 2006, 1, 109–115. [Google Scholar] [CrossRef]
  3. Ma, K.; Feng, J.; Lin, J.; Li, J. The complete mitochondrial genome of Macrobrachium nipponense. Gene 2011, 487, 160–165. [Google Scholar] [CrossRef]
  4. Zhang, X.L.; Cui, L.F.; Li, S.M.; Liu, X.Z.; Han, X.; Jiang, K.Y. Bureau Of Fisheries, Ministry of Agriculture, Prc Fisheries Economic Statistics; China Fishery Yearbook; China Agricultural Press: Beijing, China, 2020; p. 24. [Google Scholar]
  5. Jin, S.B.; Zhang, Y.; Guan, H.H.; Fu, H.T.; Jiang, S.F.; Xiong, Y.W.; Qiao, H.; Zhang, W.; Gong, Y.; Wu, Y. Histological observation of gonadal development during post-larva in oriental river prawn, Macrobrachium nipponense. Chin. J. Fish 2016, 29, 11–16. [Google Scholar]
  6. Au, H.C.; Ream-Robinson, D.; Bellew, L.A.; Broomfield, P.L.E.; Saghbini, M.; Scheffler, I.E. Structural organization of the gene encoding the human iron-sulfur subunit of succinate dehydrogenase. Gene 1995, 159, 249–253. [Google Scholar] [CrossRef]
  7. Kita, K.; Oya, H.; Gennis, R.B.; Ackrell, B.A.C.; Kasahara, M. Human complex II (succinate-ubiquinone oxidoreductase): cDNA cloning of iron sulfur (Ip) subunit of liver mitochondria. Biochem. Biophys. Res. Commun. 1990, 166, 101–108. [Google Scholar] [CrossRef] [PubMed]
  8. Gaździk, T.; Kamiński, M.; Plewka, D.; Plewka, A. Evaluation of localization and intensity of the reactions of succinic dehydrogenase, NADH tetrazole reductase and lactate dehydrogenase in developing rat testis: I. Physiological conditions. Acta Histochem. 1986, 78, 1–6. [Google Scholar] [CrossRef]
  9. Broomfield, P.L.E.; Hargreaves, J.A. A single amino-acid change in the iron-sulphur protein subunit of succinate dehydrogenase confers resistance to carboxin in Ustilago maydis. Curr. Genet. 1992, 22, 117–121. [Google Scholar] [CrossRef] [PubMed]
  10. Huang, J.; Lemire, B.D. Mutations in the C. elegans succinate dehydrogenase iron–sulfur subunit promote superoxide generation and premature aging. J. Mol. Biol. 2009, 387, 559–569. [Google Scholar] [CrossRef]
  11. Wang, Y.; Duan, Y.; Wang, J.; Zhou, M. A new point mutation in the iron–sulfur subunit of succinate dehydrogenase confers resistance to boscalid in S clerotinia sclerotiorum. Mol. Plant Pathol. 2015, 16, 653–661. [Google Scholar]
  12. Jin, S.; Hu, Y.; Fu, H.; Jiang, S.; Xiong, Y.; Qiao, H.; Zhang, W.; Gong, Y.; Wu, Y. Identification and characterization of the succinate dehydrogenase complex iron sulfur subunit B gene in the Oriental River Prawn, Macrobrachium nipponense. Front. Genet. 2021, 12, 698318. [Google Scholar] [CrossRef]
  13. Khalaila, I.; Manor, R.; Weil, S.; Granot, Y.; Keller, R.; Sagi, A. The eyestalk–androgenic gland–testis endocrine axis in the crayfish Cherax quadricarinatus. Gen. Comp. Endocrinol. 2002, 127, 147–156. [Google Scholar] [CrossRef]
  14. Song, C.; Liu, L.; Hui, M.; Liu, Y.; Liu, H.; Cui, Z. Primary molecular basis of androgenic gland endocrine sex regulation revealed by transcriptome analysis in Eriocheir sinensis. J. Oceanol. Limnol. 2019, 37, 223–234. [Google Scholar]
  15. Okuno, A.; Hasegawa, Y.; Nishiyama, M.; Ohira, T.; Ko, R.; Kurihara, M.; Matsumoto, S.; Nagasawa, H. Preparation of an active recombinant peptide of crustacean androgenic gland hormone. Peptides 2002, 23, 567–572. [Google Scholar] [PubMed]
  16. Sroyraya, M.; Chotwiwatthanakun, C.; Stewart, M.J.; Soonklang, N.; Kornthong, N.; Phoungpetchara, I.; Hanna, P.J.; Sobhon, P. Bilateral eyestalk ablation of the blue swimmer crab, Portunus pelagicus, produces hypertrophy of the androgenic gland and an increase of cells producing insulin-like androgenic gland hormone. Tissue Cell 2010, 42, 293–300. [Google Scholar] [CrossRef]
  17. Sagi, A.; Cohen, D.A.N.; Milner, Y. Effect of androgenic gland ablation on morphotypic differentiation and sexual characteristics of male freshwater prawns, Macrobrachium rosenbergii. Gen. Comp. Endocrinol. 1990, 77, 15–22. [Google Scholar] [CrossRef]
  18. Huang, X.; Ye, H.; Huang, H.; Yang, Y.; Gong, J. An insulin-like androgenic gland hormone gene in the mud crab, Scylla paramamosain, extensively expressed and involved in the processes of growth and female reproduction. Gen. Comp. Endocrinol. 2014, 204, 229–238. [Google Scholar] [PubMed]
  19. Li, S.; Li, F.; Sun, Z.; Xiang, J. Two spliced variants of insulin-like androgenic gland hormone gene in the Chinese shrimp, Fenneropenaeus chinensis. Gen. Comp. Endocrinol. 2012, 177, 246–255. [Google Scholar] [CrossRef]
  20. Liu, F.; Shi, W.; Ye, H.; Liu, A.; Zhu, Z. RNAi reveals role of insulin-like androgenic gland hormone 2 (IAG2) in sexual differentiation and growth in hermaphrodite shrimp. Front. Mar. Sci. 2021, 8, 666763. [Google Scholar] [CrossRef]
  21. Ma, K.-Y.; Li, J.-L.; Qiu, G.-F. Identification of putative regulatory region of insulin-like androgenic gland hormone gene (IAG) in the prawn Macrobrachium nipponense and proteins that interact with IAG by using yeast two-hybrid system. Gen. Comp. Endocrinol. 2016, 229, 112–118. [Google Scholar]
  22. Zhou, T.; Wang, W.; Wang, C.; Sun, C.; Shi, L.; Chan, S.F. Insulin-like androgenic gland hormone from the shrimp Fenneropenaeus merguiensis: Expression, gene organization and transcript variants. Gene 2021, 782, 145529. [Google Scholar]
  23. Farhadi, A.; Cui, W.; Zheng, H.; Li, S.; Zhang, Y.; Ikhwanuddin, M.; Ma, H. The regulatory mechanism of sexual development in decapod crustaceans. Front. Mar. Sci. 2021, 8, 679687. [Google Scholar]
  24. Farhadi, A.; Fang, S.; Zhang, Y.; Cui, W.; Fang, H.; Ikhwanuddin, M.; Ma, H. The significant sex-biased expression pattern of Sp-Wnt4 provides novel insights into the ovarian development of mud crab (Scylla Paramamosain). Int. J. Biol. Macromol. 2021, 183, 490–501. [Google Scholar]
  25. Sun, R.; Li, Y. A sex-reversing factor: Insulin-like androgenic gland hormone in decapods. Rev. Aquac. 2021, 13, 1352–1366. [Google Scholar]
  26. Farhadi, A. Sex determination and developmental mechanism of crustaceans and shellfish, volume II. Front. Endocrinol. 2023, 14, 1155209. [Google Scholar] [CrossRef] [PubMed]
  27. Ventura, T.; Manor, R.; Aflalo, E.D.; Weil, S.; Rosen, O.; Sagi, A. Timing sexual differentiation: Full functional sex reversal achieved through silencing of a single insulin-like gene in the prawn, Macrobrachium rosenbergii. Biol. Reprod. 2012, 86, 90–91. [Google Scholar] [CrossRef] [PubMed]
  28. Jin, S.; Hu, Y.; Fu, H.; Sun, S.; Jiang, S.; Xiong, Y.; Qiao, H.; Zhang, W.; Gong, Y.; Wu, Y. Analysis of testis metabolome and transcriptome from the oriental river prawn (Macrobrachium nipponense) in response to different temperatures and illumination times. Comp. Biochem. Physiol. Part D Genom. Proteom. 2020, 34, 100662. [Google Scholar]
  29. Qiao, H.; Fu, H.; Jin, S.; Wu, Y.; Jiang, S.; Gong, Y.; Xiong, Y. Constructing and random sequencing analysis of normalized cDNA library of testis tissue from oriental river prawn (Macrobrachium nipponense). Comp. Biochem. Physiol. Part D Genom. Proteom. 2012, 7, 268–276. [Google Scholar]
  30. Zhang, Y.P.; Qiao, H.; Zhang, W.Y.; Sun, S.M.; Jiang, S.F.; Gong, Y.S.; Xiong, Y.W.; Jin, S.B.; Fu, H.T. Molecular cloning and expression analysis of two sex-lethal homolog genes during development in the oriental river prawn, Macrobrachium nipponense. Genet. Mol. Res 2013, 12, 4698–4711. [Google Scholar] [CrossRef]
  31. Zhang, Y.; Jiang, S.; Xiong, Y.; Shengming, S.U.N.; Hui, Q.; Shubo, J.; Yongsheng, G.; Hongtuo, F.U. Molecular cloning and expression analysis of extra sex combs gene during development in Macrobrachium nipponense. Turk. J. Fish. Aquat. Sci. 2013, 13, 2. [Google Scholar] [CrossRef]
  32. Zhang, Y.P.; Fu, H.T.; Qiao, H.; Jin, S.B.; Gong, Y.S.; Jiang, S.F.; Xiong, Y.W. cDNA cloning, characterization and expression analysis of a transformer-2 gene in the oriental river prawn, Macrobrachium nipponense. J. World Aquacult. Soc 2013, 44, 338–349. [Google Scholar]
  33. Kusaba, M. RNA interference in crop plants. Curr. Opin. Biotechnol. 2004, 15, 139–143. [Google Scholar] [CrossRef] [PubMed]
  34. Jin, S.; Zhang, W.; Xiong, Y.; Jiang, S.; Qiao, H.; Gong, Y.; Wu, Y.; Fu, H. Identification of Important Genes Involved in the Sex-Differentiation Mechanism of Oriental River Prawn, Macrobrachium nipponense, During the Gonad Differentiation and Development Period. Front. Genet. 2022, 13, 797796. [Google Scholar] [CrossRef] [PubMed]
  35. Jin, S.; Fu, Y.; Hu, Y.; Fu, H.; Jiang, S.; Xiong, Y.; Qiao, H.; Zhang, W.; Gong, Y.; Wu, Y. Transcriptome profiling analysis of the testis after eyestalk ablation for selection of the candidate genes involved in the male sexual development in Macrobrachium nipponense. Front. Genet. 2021, 12, 675928. [Google Scholar] [PubMed]
  36. Ahn, S.W.; Gang, G.-T.; Tadi, S.; Nedumaran, B.; Kim, Y.D.; Park, J.H.; Kweon, G.R.; Koo, S.-H.; Lee, K.; Ahn, R.-S. Phosphoenolpyruvate carboxykinase and glucose-6-phosphatase are required for steroidogenesis in testicular Leydig cells. J. Biol. Chem. 2012, 287, 41875–41887. [Google Scholar]
  37. Oláh, J.; Orosz, F.; Puskás, L.G.; Hackler Jr, L.; Horányi, M.; Polgár, L.; Hollán, S.; Ovádi, J. Triosephosphate isomerase deficiency: Consequences of an inherited mutation at mRNA, protein and metabolic levels. Biochem. J. 2005, 392, 675–683. [Google Scholar] [CrossRef]
  38. Moon, J.-S.; Jin, W.-J.; Kwak, J.-H.; Kim, H.-J.; Yun, M.-J.; Kim, J.-W.; Park, S.W.; Kim, K.-S. Androgen stimulates glycolysis for de novo lipid synthesis by increasing the activities of hexokinase 2 and 6-phosphofructo-2-kinase/fructose-2, 6-bisphosphatase 2 in prostate cancer cells. Biochem. J. 2011, 433, 225–233. [Google Scholar]
  39. Brooks, D.E. Activity and androgenic control of glycolytic enzymes in the epididymis and epididymal spermatozoa of the rat. Biochem. J. 1976, 156, 527–537. [Google Scholar] [CrossRef]
  40. Singhal, R.L.; Ling, G.M. Metabolic control mechanisms in mammalian systems. IV. Androgenic induction of hexokinase and glucose-6-phosphate dehydrogenase in rat seminal vesicles. Can. J. Physiol. Pharmacol. 1969, 47, 233–239. [Google Scholar] [CrossRef]
  41. du Plessis, S.S.; Agarwal, A.; Mohanty, G.; Van der Linde, M. Oxidative phosphorylation versus glycolysis: What fuel do spermatozoa use? Asian J. Androl. 2015, 17, 230. [Google Scholar]
  42. Rato, L.; Alves, M.G.; Socorro, S.; Duarte, A.I.; Cavaco, J.E.; Oliveira, P.F. Metabolic regulation is important for spermatogenesis. Nat. Rev. Urol. 2012, 9, 330–338. [Google Scholar] [CrossRef]
  43. Li, F.; Qiao, H.; Fu, H.; Sun, S.; Zhang, W.; Jin, S.; Jiang, S.; Gong, Y.; Xiong, Y.; Wu, Y. Identification and characterization of opsin gene and its role in ovarian maturation in the oriental river prawn Macrobrachium nipponense. Comp. Biochem. Physiol. Part B Biochem. Mol. Biol. 2018, 218, 1–12. [Google Scholar]
  44. Grinshpan, N.; Abayed, F.A.A.; Wahl, M.; Ner-Gaon, H.; Manor, R.; Sagi, A.; Shay, T. The transcriptional landscape of the giant freshwater prawn: Embryonic development and early sexual differentiation mechanisms. Front. Endocrinol. 2022, 13, 1059936. [Google Scholar]
  45. Ventura, T.; Manor, R.; Aflalo, E.D.; Weil, S.; Khalaila, I.; Rosen, O.; Sagi, A. Expression of an androgenic gland-specific insulin-like peptide during the course of prawn sexual and morphotypic differentiation. Int. Sch. Res. Not. 2011, 2011, 476283. [Google Scholar]
  46. Ventura, T.; Manor, R.; Aflalo, E.D.; Weil, S.; Raviv, S.; Glazer, L.; Sagi, A. Temporal silencing of an androgenic gland-specific insulin-like gene affecting phenotypical gender differences and spermatogenesis. Endocrinology 2009, 150, 1278–1286. [Google Scholar] [PubMed]
  47. Rosen, O.; Manor, R.; Weil, S.; Gafni, O.; Linial, A.; Aflalo, E.D.; Ventura, T.; Sagi, A. A sexual shift induced by silencing of a single insulin-like gene in crayfish: Ovarian upregulation and testicular degeneration. PLoS ONE 2010, 5, e15281. [Google Scholar]
  48. Choi, Y.J.; Kim, S.; Choi, Y.; Nielsen, T.B.; Yan, J.; Lu, A.; Ruan, J.; Lee, H.-R.; Wu, H.; Spellberg, B. SERPINB1-mediated checkpoint of inflammatory caspase activation. Nat. Immunol. 2019, 20, 276–287. [Google Scholar]
  49. Denolet, E.; De Gendt, K.; Allemeersch, J.; Engelen, K.; Marchal, K.; Van Hummelen, P.; Tan, K.A.L.; Sharpe, R.M.; Saunders, P.T.K.; Swinnen, J.V. The effect of a Sertoli cell-selective knockout of the androgen receptor on testicular gene expression in prepubertal mice. Mol. Endocrinol. 2006, 20, 321–334. [Google Scholar] [CrossRef]
  50. Guo, T.; Zuo, Y.; Qian, L.; Liu, J.; Yuan, Y.; Xu, K.; Miao, Y.; Feng, Q.; Chen, X.; Jin, L. ADP-ribosyltransferase PARP11 modulates the interferon antiviral response by mono-ADP-ribosylating the ubiquitin E3 ligase β-TrCP. Nat. Microbiol. 2019, 4, 1872–1884. [Google Scholar]
  51. Meyer-Ficca, M.L.; Ihara, M.; Bader, J.J.; Leu, N.A.; Beneke, S.; Meyer, R.G. Spermatid head elongation with normal nuclear shaping requires ADP-ribosyltransferase PARP11 (ARTD11) in mice. Biol. Reprod. 2015, 92, 80–81. [Google Scholar] [CrossRef]
  52. Brown, A.C.; Harrison, L.M.; Kapulkin, W.; Jones, B.F.; Sinha, A.; Savage, A.; Villalon, N.; Cappello, M. Molecular cloning and characterization of a C-type lectin from Ancylostoma ceylanicum: Evidence for a role in hookworm reproductive physiology. Mol. Biochem. Parasitol. 2007, 151, 141–147. [Google Scholar] [CrossRef] [PubMed]
  53. Goluboff, E.T.; Mertz, J.R.; Tres, L.L.; Kierszenbaum, A.L. Galactosyl receptor in human testis and sperm is antigenically related to the minor C-type (Ca2+-dependent) lectin variant of human and rat liver. Mol. Reprod. Dev. 1995, 40, 460–466. [Google Scholar] [CrossRef] [PubMed]
  54. Qiu, X.B.; Shao, Y.M.; Miao, S.; Wang, L. The diversity of the DnaJ/Hsp40 family, the crucial partners for Hsp70 chaperones. Cell. Mol. Life Sci. CMLS 2006, 63, 2560–2570. [Google Scholar] [CrossRef] [PubMed]
  55. Selvaraju, S.; Parthipan, S.; Somashekar, L.; Kolte, A.P.; Krishnan Binsila, B.; Arangasamy, A.; Ravindra, J.P. Occurrence and functional significance of the transcriptome in bovine (Bos taurus) spermatozoa. Sci. Rep. 2017, 7, 42392. [Google Scholar] [CrossRef]
  56. Banday, S.; Pandita, R.K.; Mushtaq, A.; Bacolla, A.; Mir, U.S.; Singh, D.K.; Jan, S.; Bhat, K.P.; Hunt, C.R.; Rao, G. Autism-associated vigilin depletion impairs DNA damage repair. Mol. Cell. Biol. 2021, 41, e00082-00021. [Google Scholar] [CrossRef]
  57. Lyon, M.F. Transmission ratio distortion in mice. Annu. Rev. Genet. 2003, 37, 393–408. [Google Scholar] [CrossRef]
  58. Zhang, S.B.; Jiang, P.; Wang, Z.Q.; Long, S.R.; Liu, R.D.; Zhang, X.; Yang, W.; Ren, H.J.; Cui, J. DsRNA-mediated silencing of Nudix hydrolase in Trichinella spiralis inhibits the larval invasion and survival in mice. Exp. Parasitol. 2016, 162, 35–42. [Google Scholar] [CrossRef]
  59. Jiang, F.W.; Fu, H.T.; Qiao, H.; Zhang, W.Y.; Jiang, S.F.; Xiong, Y.W.; Sun, S.M.; Gong, Y.S.; Jin, S.B. The RNA interference regularity of transformer-2 gene of oriental river prawn Macrobrachium nipponense. Chin. Agric. Sci. Bull 2014, 30, 32–37. [Google Scholar]
  60. Li, F.; Bai, H.; Xiong, Y.; Fu, H.; Jiang, S.; Jiang, F.; Jin, S.; Sun, S.; Qiao, H.; Zhang, W. Molecular characterization of insulin-like androgenic gland hormone-binding protein gene from the oriental river prawn Macrobrachium nipponense and investigation of its transcriptional relationship with the insulin-like androgenic gland hormone gene. Gen. Comp. Endocrinol. 2015, 216, 152–160. [Google Scholar] [CrossRef]
  61. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  62. Kim, D.; Paggi, J.M.; Park, C.; Bennett, C.; Salzberg, S.L. Graph-based genome alignment and genotyping with HISAT2 and HISAT-genotype. Nat. Biotechnol. 2019, 37, 907–915. [Google Scholar] [CrossRef] [PubMed]
  63. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  64. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed]
  65. Ferreira, J.A.; Zwinderman, A.H. On the benjamini–hochberg method. Ann. Statist. 2006, 34, 1827–1849. [Google Scholar] [CrossRef]
  66. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T. Gene ontology: Tool for the unification of biology. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [PubMed]
  67. Kanehisa, M.; Goto, S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  68. Hu, Y.; Fu, H.; Qiao, H.; Sun, S.; Zhang, W.; Jin, S.; Jiang, S.; Gong, Y.; Xiong, Y.; Wu, Y. Validation and evaluation of reference genes for quantitative real-time PCR in Macrobrachium Nipponense. Int. J. Mol. Sci. 2018, 19, 2258. [Google Scholar] [CrossRef]
  69. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
Figure 1. GO classification of the unigenes and DEGs discovered: (A) “SDHB_AG vs. C_AG”; (B) “SDHB_T vs. C_T”. The abscissa shows the second level terms in the three GO categories. The ordinates show the number of genes annotated to each term and the percentage of all genes.
Figure 1. GO classification of the unigenes and DEGs discovered: (A) “SDHB_AG vs. C_AG”; (B) “SDHB_T vs. C_T”. The abscissa shows the second level terms in the three GO categories. The ordinates show the number of genes annotated to each term and the percentage of all genes.
Ijms 24 13176 g001
Figure 2. KEGG enrichment of DEGs: (A) “SDHB_AG vs. C_AG”; (B) “SDHB_T vs. C_T”. The size of the dots represents the number of genes. The color of each dot represents the p-value. The abscissa shows the enrichment score. The ordinate shows the number of genes annotated to each term and the percentage of all genes.
Figure 2. KEGG enrichment of DEGs: (A) “SDHB_AG vs. C_AG”; (B) “SDHB_T vs. C_T”. The size of the dots represents the number of genes. The color of each dot represents the p-value. The abscissa shows the enrichment score. The ordinate shows the number of genes annotated to each term and the percentage of all genes.
Ijms 24 13176 g002
Figure 3. Schematic diagram of a partial glycolysis/gluconeogenesis pathway. The green boxes represent genes that are downregulated in the androgenic gland of Macrobrachium nipponense.
Figure 3. Schematic diagram of a partial glycolysis/gluconeogenesis pathway. The green boxes represent genes that are downregulated in the androgenic gland of Macrobrachium nipponense.
Ijms 24 13176 g003
Figure 4. qRT-PCR validation of RNA-Seq data: (A) “SDHB_AG vs. C_AG”; (B) “SDHB_T vs. C_T”. r is the Pearson correlation coefficient between RNA-Seq group and qRT-PCR group.
Figure 4. qRT-PCR validation of RNA-Seq data: (A) “SDHB_AG vs. C_AG”; (B) “SDHB_T vs. C_T”. r is the Pearson correlation coefficient between RNA-Seq group and qRT-PCR group.
Ijms 24 13176 g004
Table 1. Statistics of the RNA-seq data in each of the three replicates.
Table 1. Statistics of the RNA-seq data in each of the three replicates.
SampleClean Bases (Gb)Q30 (%)GC (%)Mapped Reads (%)Unigene Number
C_AG17.1592.4945.8991.5816,623
C_AG26.8992.4745.3491.2316,731
C_AG36.9492.4845.1491.7617,035
C_T17.0392.6843.1391.9818,951
C_T26.7692.7843.1191.2819,421
C_T36.9692.7243.3791.3220,053
SDHB_AG16.9892.5845.792.1916,623
SDHB_AG26.8692.4245.3392.0416,731
SDHB_AG37.1492.8745.3291.3117,035
SDHB_T17.392.6143.5991.6218,951
SDHB_T27.1392.6543.5991.4819,421
SDHB_T37.292.8443.4892.2220,053
Table 2. Male development-related genes screened from the differentially expressed genes.
Table 2. Male development-related genes screened from the differentially expressed genes.
Gene NameAccession NumberDifferential Expression
SDHB_AG vs. C_AGSDHB_T vs. C_T
Pro-resilinXM_043027950.1downup
IAGXM_045267991.1down
SERPINB1XM_045251131.1down
PARP11XM_047641508.1up
DNAJC2XM_042365393.1 down
CTL-1XM_042349745.1 down
YESXM_043000203.1 up
VigilinXM_027372514.1 up
SMOK-Y-likeXM_027369036.1 up
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gao, Z.; Zhang, W.; Jiang, S.; Yuan, H.; Cai, P.; Jin, S.; Fu, H. Identification of Male Sex-Related Genes Regulated by SDHB in Macrobrachium nipponense Based on Transcriptome Analysis after an RNAi Knockdown. Int. J. Mol. Sci. 2023, 24, 13176. https://doi.org/10.3390/ijms241713176

AMA Style

Gao Z, Zhang W, Jiang S, Yuan H, Cai P, Jin S, Fu H. Identification of Male Sex-Related Genes Regulated by SDHB in Macrobrachium nipponense Based on Transcriptome Analysis after an RNAi Knockdown. International Journal of Molecular Sciences. 2023; 24(17):13176. https://doi.org/10.3390/ijms241713176

Chicago/Turabian Style

Gao, Zijian, Wenyi Zhang, Sufei Jiang, Huwei Yuan, Pengfei Cai, Shubo Jin, and Hongtuo Fu. 2023. "Identification of Male Sex-Related Genes Regulated by SDHB in Macrobrachium nipponense Based on Transcriptome Analysis after an RNAi Knockdown" International Journal of Molecular Sciences 24, no. 17: 13176. https://doi.org/10.3390/ijms241713176

APA Style

Gao, Z., Zhang, W., Jiang, S., Yuan, H., Cai, P., Jin, S., & Fu, H. (2023). Identification of Male Sex-Related Genes Regulated by SDHB in Macrobrachium nipponense Based on Transcriptome Analysis after an RNAi Knockdown. International Journal of Molecular Sciences, 24(17), 13176. https://doi.org/10.3390/ijms241713176

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop