Validation and Evaluation of Reference Genes for Quantitative Real-Time PCR in Macrobrachium Nipponense

Identification of the effects of alkalinity exposure on the gills of oriental river prawns, Macrobrachium nipponense

Macrobrachium nipponense is an important commercial freshwater species in China. However, the ability of alkali tolerance of M. nipponense is insufficient to culture in the major saline-alkali water source in China. Thus, it is urgently needed to perform the genetic improvement of alkali tolerance in this species. In the present study, we aimed to analyse the effects of alkali treatment on gills in this species after 96 h alkalinity exposure under the alkali concentrations of 0 mmol/L, 4 mmol/L, 8 mmol/L, and 12 mmol/L through performing the histological observations, measurement of antioxidant enzymes, metabolic profiling analysis, and transcriptome profiling analysis. The results of the present study revealed that alkali treatment stimulated the contents of malondialdehyde, glutathione, glutathione peroxidase in gills, indicating these antioxidant enzymes plays essential roles in the protection of body from the damage, caused by the alkali treatment. In addition, high concentration of alkali treatment (> 8 mmol/L) resulted in the damage of gill membrane and haemolymph vessel, affecting the normal respiratory function of gill. Metabolic profiling analysis revealed that Metabolic pathways, Biosynthesis of secondary metabolites, Biosynthesis of plant secondary metabolites, Microbial metabolism in diverse environments, Biosynthesis of amino acids were identified as the main enriched metabolic pathways of differentially expressed metabolites, which are consistent with the previous publications, treated by the various environmental factors. Transcriptome profiling analyses revealed that the alkali concentration of 12 mmol/L has more regulatory effects on the changes of gene expression than the other alkali concentrations. KEGG analysis revealed that Phagosome, Lysosome, Glycolysis/Gluconeogenesis, Purine Metabolism, Amino sugar and nucleotide sugar metabolism, and Endocytosis were identified as the main enriched metabolic pathways in the present study, predicting these metabolic pathways may be involved in the adaption of alkali treatment in M. nipponense. Phagosome, Lysosome, Purine Metabolism, and Endocytosis are immune-related metabolic pathways, while Glycolysis/Gluconeogenesis, and Amino sugar and nucleotide sugar metabolism are energy metabolism-related metabolic pathways. Quantitative PCR analyses of differentially expressed genes (DEGs) verified the accuracy of the RNA-Seq. Alkali treatment significantly stimulated the expressions of DEGs from the metabolic pathways of Phagosome and Lysosome, suggesting Phagosome and Lysosome play essential roles in the regulation of alkali tolerance in this species, as well as the genes from these metabolic pathways. The present study identified the effects of alkali treatment on gills, providing valuable evidences for the genetic improvement of alkali tolerance in M. nipponense.

Functional Study of the Role of the Methyl Farnesoate Epoxidase Gene in the Ovarian Development of Macrobrachium nipponense

Methyl farnesoate epoxidase (MFE) is a gene encoding an enzyme related to the last step of juvenile hormone biosynthesis. Mn-MFE cDNA has a total length of 1695 bp and an open reading frame (ORF) length of 1482 bp, encoding 493 amino acids. Sequence analysis showed that its amino acid sequence has a PPGP hinge, an FGCG structural domain, and other structural domains specific to the P450 family of enzymes. Mn-MFE was most highly expressed in the hepatopancreas, followed by the ovary and gill, weakly expressed in heart and muscle tissue, and barely expressed in the eyestalk and cranial ganglion. Mn-MFE expression remained stable during the larval period, during which it mainly played a critical role in gonadal differentiation. Expression in the ovary was positively correlated and expression in the hepatopancreas was negatively correlated with ovarian development. In situ hybridization (ISH) showed that the signal was expressed in the oocyte, nucleus, cell membrane and follicular cells, and the intensity of expression was strongest at stage O-IV. The knockdown of Mn-MFE resulted in a significantly lower gonadosomatic index and percentage of ovaries past stage O-III compared to the control group. However, no differences were found in the cumulative frequency of molting between the experimental and control groups. Moreover, the analysis of ovarian tissue sections at the end of the experiment showed differences between groups in development speed but not in subcellular structure. These results demonstrate that Mn-MFE promotes the ovarian development of Macrobrachium nipponense adults but has no effect on molting.

Reference genes for qPCR expression in black tiger shrimp, Penaeus monodon

BACKGROUND

Gene expression profiling via qPCR is an essential tool for unraveling the intricate molecular mechanisms underlying growth and development. Identifying and validating the most appropriate reference genes is essential for qPCR experiments. Nevertheless, there exists a deficiency in a thorough assessment of reference genes concerning the expression of the genes in the research in the context of the growth and development of the Black Tiger Shrimp, P. monodon. This popular marine crustacean is extensively raised for human consumption. In this study, we assessed the expression stability of seven reference genes (ACTB, 18S, EF-1α, AK, PK, cox1, and CLTC) in adult tissues (hepatopancreas, gills, and stomach) of small and large polymorphs of P. monodon.

METHODS AND RESULTS

The stability of gene expressions was assessed utilizing NormFinder, BestKeeper, and geNorm, and a comprehensive ranking of these genes was conducted through the online tool RefFinder. In the overall ranking, 18S and CLTC emerged as the most stable genes in the hepatopancreas and stomach, while CLTC and AK exhibited significant statistical reliability in the gills of adult P. monodon. The validation of these identified stable genes was carried out using a growth-associated gene, insr-1.

CONCLUSION

The results indicated that 18S and CLTC stand out as the most versatile reference genes for conducting qPCR analysis focused on the growth of P. monodon. This study represents the first comprehensive exploration that identifies and assesses reference genes for qPCR analysis in P. monodon, providing valuable tools for research involving similar crustaceans.

Identification of Candidate Male-Reproduction-Related Genes from the Testis and Androgenic Gland of Macrobrachium nipponense, Regulated by PDHE1, through Transcriptome Profiling Analysis

The previous publication identified that pyruvate dehydrogenase E1 (PDHE1) positively regulated the process of male reproduction in M. nipponense through affecting the expressions of insulin-like androgenic gland hormone. The present study aimed to identify the potential male-reproduction-related genes that were regulated by PDHE1 through performing the transcriptome profiling analysis in the testis and androgenic gland after the knockdown of the expressions of PDHE1 by the injection of dsPDHE1. Both RNA-Seq and qPCR analysis identified the significant decreases in PDHE1 expressions in the testis and androgenic gland in dsPDHE1-injected prawns compared to those in dsGFP-injected prawns, indicating the efficiency of dsPDHE1 in the present study. Transcriptome profiling analysis identified 56 and 127 differentially expressed genes (DEGs) in the testis and androgenic gland, respectively. KEGG analysis revealed that the energy-metabolism-related pathways represented the main enriched metabolic pathways of DEGs in both the testis and androgenic gland, including pyruvate metabolism, the Citrate cycle (TCA cycle), Glycolysis/Gluconeogenesis, and the Glucagon signaling pathway. Thus, it is predicted that these metabolic pathways and the DEGs from these metabolic pathways regulated by PDHE1 may be involved in the regulation of male reproduction in M. nipponense. Furthermore, four genes were found to be differentially expressed in both the testis and androgenic gland, of which ribosomal protein S3 was down-regulated and uncharacterized protein LOC113829596 was up-regulated in both the testis and androgenic gland in dsPDHE1-injected prawns. The present study provided valuable evidence for the establishment of an artificial technique to regulate the process of male reproduction in M. nipponense.

Evaluation of suitable reference genes for expression profile analyses of target genes in the coffee berry borer, Hypothenemus hampei (Ferrari) (Coleoptera: Curculionidae)

The coffee berry borer, Hypothenemus hampei (Ferrari) (Coleoptera: Curculionidae), is a major destructive insect pest of coffee, which impacts the coffee crops negatively. As a draft genome has been completed for this insect, most molecular studies on gene transcriptional levels under different experimental conditions will be conducted using real-time reverse-transcription quantitative polymerase chain reactions (RT-qPCR). However, the lack of suitable internal reference genes will affect the accuracy of RT-qPCR results. In this study, the expression stability of nine candidate reference genes was evaluated under different developmental stages, temperature stress, and Beauveria bassiana infection. Data analyses were completed by four commonly used programs, BestKeeper, NormFinder, geNorm, and RefFinder. The result showed that RPL3 and EF1α combination were recommended as the most stable reference genes for developmental stages. EF1α and RPS3a combination were the top two stable reference genes for B. bassiana infection. RPS3a and RPL3 combination performed as the optimal reference genes both in temperature stress and all samples. Our results should provide a good foundation for the expression profile analyses of target genes in the future, especially for molecular studies on insect genetic development, temperature adaptability, and immune mechanism to entomogenous fungi in H. hampei.

Effects of Alkalinity Exposure on Antioxidant Status, Metabolic Function, and Immune Response in the Hepatopancreas of Macrobrachium nipponense

The oriental river prawn Macrobrachium nipponense is an important freshwater economic species in China, producing huge economic benefits. However, M. nipponense shows lower alkali tolerance than fish species, thus genetic selection is urgently needed in order to improve alkali tolerance in this species. In the present study, the effects of alkalinity exposure on the hepatopancreas of M. nipponense were measured under the alkali concentrations of 0 (control), 4, 8, and 12 mmol/L with the exposure time of 96 h through histological observations, measurement of antioxidant enzymes, metabolic profiling analysis, and transcriptome profiling analysis. The present study identified that the low concentration of alkali treatment (<4 mmol/L) did not result in morphological changes in the hepatopancreas and activity changes in antioxidant enzymes, while high-alkali treatment (>8 mmol/L) damaged the normal structures of the lumen and vacuoles and significantly stimulated the levels of superoxide dismutase, catalase, and total antioxidant capacity, indicating these antioxidant enzymes play essential roles in the protection of the body from the damage caused by the alkali treatment. Metabolic profiling analysis revealed that the main enriched metabolic pathways of differentially expressed metabolites in the present study were consistent with the metabolic pathways caused by environmental stress in plants and other aquatic animals. Transcriptome profiling analysis revealed that the alkali concentration of <8 mmol/L did not lead to significant changes in gene expression. The main enriched metabolic pathways were selected from the comparison between 0 mmol/L vs. 12 mmol/L, and some significantly up-regulated genes were selected from these metabolic pathways, predicting these selected metabolic pathways and genes are involved in the adaptation to alkali treatment in M. nipponense. The expressions of Ras-like GTP-binding protein, Doublesex and mab-3 related transcription factor 1a, and Hypothetical protein JAY84 are sensitive to changes in alkali concentrations, suggesting these three genes participated in the process of alkali adaptation in M. nipponense. The present study identified the effects of alkalinity exposure on the hepatopancreas of M. nipponense, including the changes in antioxidant status and the expressions of metabolites and genes, contributing to further studies of alkali tolerance in this species.

Identification of genes regulated by 20-Hydroxyecdysone in Macrobrachium nipponense using comparative transcriptomic analysis

BACKGROUND

Macrobrachium nipponense is a freshwater prawn of economic importance in China. Its reproductive molt is crucial for seedling rearing and directly impacts the industry's economic efficiency. 20-hydroxyecdysone (20E) controls various physiological behaviors in crustaceans, among which is the initiation of molt. Previous studies have shown that 20E plays a vital role in regulating molt and oviposition in M. nipponense. However, research on the molecular mechanisms underlying the reproductive molt and role of 20E in M. nipponense is still limited.

RESULTS

A total of 240.24 Gb of data was obtained from 18 tissue samples by transcriptome sequencing, with > 6 Gb of clean reads per sample. The efficiency of comparison with the reference transcriptome ranged from 87.05 to 92.48%. A total of 2532 differentially expressed genes (DEGs) were identified. Eighty-seven DEGs associated with molt or 20E were screened in the transcriptomes of the different tissues sampled in both the experimental and control groups. The reliability of the RNA sequencing data was confirmed using Quantitative Real-Time PCR. The expression levels of the eight strong candidate genes showed significant variation at the different stages of molt.

CONCLUSION

This study established the first transcriptome library for the different tissues of M. nipponense in response to 20E and demonstrated the dominant role of 20E in the molting process of this species. The discovery of a large number of 20E-regulated strong candidate DEGs further confirms the extensive regulatory role of 20E and provides a foundation for the deeper understanding of its molecular regulatory mechanisms.

Evaluation and validation of reference genes for RT-qPCR gene expression in Naegleria gruberi

Naegleria gruberi is a free-living amoeboflagellate commonly found in freshwater and in soils around the world. It is a non-pathogenic relative of Naegleria fowleri, which is the etiologic agent of Primary Amoebic Meningoencephalitis (PAM). PAM occurs world-wide and it is considered a rare disease, but its fatality rate is high (96%) mainly because of delay in initiation of treatment due to misdiagnosis and lack of a specific treatment. The analysis of gene expression by quantitative real-time PCR in N. gruberi could be a highly efficient means to understand the pathogenicity of N. fowleri and also to find drug targets. Accurate RT-qPCR analysis requires correct normalization of gene expression data using reference genes (RG), whose expression should be constant under different experimental conditions. In this study, six genes, representing the most frequently used housekeeping genes, were selected for evaluation as reference genes in N. gruberi. The expression and stability of these genes was evaluated employing four algorithms (geNorm, NormFinder, BestKeeper and RefFinder). This work shows significant variations of the stability of RGs depending on the algorithms employed and on the experimental conditions (i.e. logarithmic, stationary, heat-shock and oxidative stress). The geNorm, NormFinder and RefFinder analysis of all the experimental conditions in combination revealed that ACT and G6PD were the most stable RGs. While BestKeeper analysis showed that 18S and TBP were the most stable RGs. Moreover, normalization of HSP90 gene expression with the most stable RGs resulted in an upregulation whereas when the normalization was done with the unstable RGs, the gene expression was not reliable. Hence, the implications of this study are relevant to gene expression studies in N. gruberi.

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

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.

Mn-XRN1 Has an Inhibitory Effect on Ovarian Reproduction in Macrobrachium nipponense

XRN1 is an exoribonuclease that degrades mRNA in the cytoplasm along the 5'-3' direction. A previous study indicated that it may be involved in the reproduction of Macrobrachium nipponense. Quantitative real-time PCR was used to detect the spatiotemporal expression pattern of Mn-XRN1. At the tissue level, Mn-XRN1 was significantly expressed in the ovary. During development, Mn-XRN1 was significantly expressed at the CS stage of the embryo, on the 10th day post-larval and in the O2 stage of ovarian reproduction. The in situ hybridization results showed the location of Mn-XRN1 in the ovary. The expression of Mn-VASA was significantly increased after in vivo injection of Mn-XRN1 dsRNA. This suggests that Mn-XRN1 negatively regulates the expression of Mn-VASA. Furthermore, we counted the number of M. nipponense at various stages of ovarian reproduction on different days after RNAi. The results showed that ovarian development was significantly accelerated. In general, the results of the present study indicate that Mn-XRN1 has an inhibitory effect on the ovarian maturation of M. nipponense. The inhibitory effect might be through negative regulation of Mn-VASA.

Deciphering Molecular Mechanisms Governing the Reproductive Molt of Macrobrachium nipponense: A Transcriptome Analysis of Ovaries across Various Molting Stages

The relationship between molting and reproduction has received more attention in economically important crustacean decapods. Molting and reproduction are synergistic events in Macrobrachium nipponense, but the molecular regulatory mechanisms behind them are unclear. In the current study, we performed Illumina sequencing for the ovaries of M. nipponense during the molt cycle (pre-molting, Prm; mid-molting, Mm; and post-molting, Pom). A total of 66.57 Gb of transcriptome data were generated through sequencing, resulting in the identification of 105,149 unigenes whose alignment ratio with the reference genome exceeded 87.57%. Differentially expressed genes (DEGs) were annotated through the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases for gene classification and pathway analysis. A total of twenty-six molt-related DEGs were found, and their expression patterns were examined across various molting stages. The KEGG enrichment analysis revealed that the key pathways involved in regulating the molting process of M. nipponense primarily include the mTOR, insect hormone biosynthesis, TGF-beta, and Wnt signaling pathways. Our transcriptomic data suggest that these pathways crosstalk with each other to regulate the synthesis and degradation of ecdysone throughout the molt cycle. The current study has deepened our understanding of the molecular mechanisms of crustacean molting and will serve as a basis for future studies of crustaceans and other molting animals.

RNAi Analysis of Potential Functions of Cyclin B3 in Reproduction of Male Oriental River Prawns (Macrobrachium nipponense)

Cyclin B3 (CycB3) is involved in the metabolic pathway of the cell cycle, playing essential roles in the regulation of cell proliferation and mitosis. CycB3 is also predicted to be involved in the reproduction of male oriental river prawns (Macrobrachium nipponense). In this study, the potential functions of CycB3 in M. nipponense were investigated using quantitative real-time PCR, RNA interference, and histological observations. The full-length DNA sequence of CycB3 in M. nipponense was 2147 base pairs (bp) long. An open reading frame of 1500 bp was found, encoding 499 amino acids. A highly conserved destruction box and two conserved cyclin motifs were found in the protein sequence of Mn-CycB3. Phylogenetic tree analysis revealed that this protein sequence was evolutionarily close to that of CycB3s of crustacean species. Quantitative real-time PCR analysis results suggested that CycB3 was involved in the process of spermiogenesis, oogenesis, and embryogenesis in M. nipponense. RNA interference analysis showed that CycB3 had a positive regulatory relationship with insulin-like androgenic gland hormone (IAG) in M. nipponense. In addition, sperm were rarely observed in the testis of double-stranded CycB3-injected prawns after 14 days of treatment, and sperm abundance was dramatically lower than that in the double-stranded GFP-injected prawns on the same day. This result indicated that CycB3 can regulate the testis reproduction in M. nipponense through inhibiting the IAG expressions. Overall, these results indicated that CycB3 plays essential roles in the regulation of male reproduction in M. nipponense, which may promote the studies of male reproduction in other crustacean species.

17β-Estradiol Induced Sex Reversal and Gonadal Transcriptome Analysis in the Oriental River Prawn (Macrobrachium nipponense): Mechanisms, Pathways, and Potential Harm

Sex reversal induced by 17β-estradiol (E₂) has shown the potential possibility for monoculture technology development. The present study aimed to determine whether dietary supplementation with different concentrations of E₂ could induce sex reversal in M. nipponense, and select the sex-related genes by performing the gonadal transcriptome analysis of normal male (M), normal female (FM), sex-reversed male prawns (RM), and unreversed male prawns (NRM). Histology, transcriptome analysis, and qPCR were performed to compare differences in gonad development, key metabolic pathways, and genes. Compared with the control, after 40 days, feeding E₂ with 200 mg/kg at PL25 (PL: post-larvae developmental stage) resulted in the highest sex ratio (female: male) of 2.22:1. Histological observations demonstrated the co-existence of testis and ovaries in the same prawn. Male prawns from the NRM group exhibited slower testis development without mature sperm. RNA sequencing revealed 3702 differentially expressed genes (DEGs) between M vs. FM, 3111 between M vs. RM, and 4978 between FM vs. NRM. Retinol metabolism and nucleotide excision repair pathways were identified as the key pathways for sex reversal and sperm maturation, respectively. Sperm gelatinase (SG) was not screened in M vs. NRM, corroborating the results of the slice D. In M vs. RM, reproduction-related genes such as cathepsin C (CatC), heat shock protein cognate (HSP), double-sex (Dsx), and gonadotropin-releasing hormone receptor (GnRH) were expressed differently from the other two groups, indicating that these are involved in the process of sex reversal. Exogenous E₂ can induce sex reversal, providing valuable evidence for the establishment of monoculture in this species.

Sex Reversal Induced by Dietary Supplementation with 17α-Methyltestosterone during the Critical Period of Sex Differentiation in Oriental River Prawn (Macrobrachium nipponense)

The steroid 17α-methyltestosterone (MT) inhibits ovarian function and is often used to induce sex reversal artificially in vertebrates. In the present study, different concentrations of MT were added as dietary supplementation, and the effects on sex ratio, growth, and gonadal development were examined. After 40 days, the sex ratio (male:female) in each group increased at different degrees with 50 (1.36:1), 100 (1.57:1), and 200 (2.61:1) mg/kg MT, and neo-males with testis-ovary coexistence were observed in the 200 mg/kg MT group. Furthermore, 50 and 100 mg/kg MT could induce female reversion in neo-males. Histologically, the development of the testes in experimental groups was slower, but the ovaries of the experimental and control groups had similar developmental rates. The expression levels of DMRT11E, Foxl2, and SoxE1 in males at 200 mg/kg MT were 8.65-, 3.75-, and 3.45-fold greater than those of the control group. In crustaceans, sex reversal through vertebrate sex hormones can be observed. Neo-males (sex-reversed female prawns) were maintained by exogenous androgen, and over-reliance led to slow testis growth, small body size, and low growth rate, but sperm was still produced. In female prawns, MT inhibited ovary development and promoted growth.

Hepatopancreas Proteomic Analysis Reveals Key Proteins and Pathways in Regulatory of Ovary Maturation of Macrobrachium nipponense

A TMT-based (Tandem Mass Tag) liquid chromatography-tandem mass spectrometry (LC-MS/MS) proteomics approach was employed to explore differentially expressed proteins (DEPs) and KEGG pathways in hepatopancreas of 5 ovary stages. In total, 17,999 peptides were detected, among which 3395 proteins were identified. Further analysis revealed 26, 24, 37, and 308 DEPs in HE-I versus HE-II, HE-II versus HE-Ⅲ, HE-Ⅲ versus HE-Ⅳ, and HE-Ⅳ versus HE-Ⅴ, respectively (HE-I, HE-II, HE-III, HE-IV, and HE-V means hepatopancreas sampled from ovary stage I to V.). Gene ontology (GO) analysis indicated that DEPs were significantly enriched in "catalytic activity", "metabolic process", and "cell" of 4 comparison groups in turn. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results showed that in hepatopancreas, as the ovaries developed to maturation, carbohydrate metabolism, lipid metabolism, amino acid metabolism, and lysosome played important roles in turn. The mRNA expression of 15 selected DEPs were consistent with proteome results by qPCR analysis. Further mRNA expression investigation results suggested 4 proteins (fatty acid-binding protein, NPC intracellular cholesterol transporter 1, Serine hydroxymethyltransferase, and Crustapin) were involved in ovary maturation. These results enhance the understanding of the regulatory role of hepatopancreas in M. nipponense ovary maturation and provide new insights for understanding the crustacean regulation mechanisms.

CRISPR/Cas9 establishment-mediated targeted mutagenesis in Macrobrachium nipponense

Introduction: CRISPR/Cas9 is a gene-editing technology which could specifically cleave dsDNA and induce target gene mutation. CRISPR/Cas9 has been widely used in gene functional studies in many fields, such as medicine, biology, and agriculture due to its simple design, low cost, and high efficiency. Although it has been well developed in model fish and freshwater fish for gene function analysis, it is still novel in the studies dealing with economic crustacean species. Methods: In this study, we established a CRISPR/Cas9 system based on microinjection for M. nipponense, an important economic crustacean aquaculture species. The vitellogenin (Vg) gene and the eyeless (Ey) gene were selected as the targeted genes for mutation. Two sgRNAs were designed for Mn-Vg and Mn-Ey gene editing, respectively. Results and Discussion: For sg-Vg-1, the gastrula survival ratio was 8.69%, and the final hatching ratio was 4.83%. The blastula mutant ratio was 10%, and the hatching individual mutant ratio was 30%. For sg-Vg-2, the gastrula survival ratio was 5.85%, and the final hatching ratio was 3.89%. The blastula mutant ratio was 16.67%, and no mutant sequences were detected in hatching individuals. For sg-Ey-1, the gastrula survival ratio was 6.25%, and the final hatching ratio was 2.34%. The blastula mutant ratio was 10.00%, and the hatching individual mutant ratio was 66.67%. For sg-Ey-2, the gastrula survival ratio was 6.00%, and the final hatching ratio was 2.67%. No mutant sequence was detected in both blastula stage and hatching individuals. There were no significant morphological changes observed in the Mn-Vg group. Two deformed types were detected in sg-Ey-1-injected embryos. An evident developmental delay of the compound eye was detected in Ey-sg1-H1 in the zoea stage. The compound eyes of the Ey-sg1-H2 embryo could not form well-defined spheres, and the whole compound eye appeared to diffuse at the end of the late zoea stage. The establishment of a gene-editing platform based on CRISPR/Cas9 will not only provide an efficient and convenient method for gene function analysis but also provide a powerful tool for molecular-assisted breeding of Macrobrachium nipponense.

Transcriptome analysis provides novel insights into the immune mechanisms of Macrobrachium nipponense during molting

Molting is a basic physiological behavior of the Oriental river prawn (Macrobrachium nipponense), however, the gene expression patterns and immune mechanisms during the molting process of Oriental river prawn are unclear. In the current study, the gene expression levels of the hepatopancreas of the Oriental river prawn at different molting stages (pre-molting, Prm; mid-molting, Mm; and post-molting, Pom) were detected by mRNA sequencing. A total of 1721, 551, and 1054 differentially expressed genes (DEGs) were identified between the Prm hepatopancreas (PrmHe) and Mm hepatopancreas (MmHe), MmHe and Pom hepatopancreas (PomHe) and PrmHe and PomHe, respectively. The results showed that a total of 1151 DEGs were annotated into 316 signaling pathways, and the significantly enriched immune-related pathways were "Lysosome", "Hippo signaling pathway", "Apoptosis", "Autophagy-animal", and "Endocytosis". The qRT-PCR verification results of 30 randomly selected DEGs were consistent with RNA-seq. The expression patterns of eight immune related genes in different molting stages of the Oriental river prawn were analyzed by qRT-PCR. The function of Caspase-1 (CASP1) was further investigated by bioinformatics, qRT-PCR, and RNAi analysis. CASP1 has two identical conserved domains: histidine active site and pentapeptide motif, and the expression of CASP1 is the highest in ovary. The expression levels of triosephosphate isomerase (TPI), Cathepsin B (CTSB) and Hexokinase (HXK) were evaluated after knockdown of CASP1. This research provides a valuable basis to improve our understanding the immune mechanisms of Oriental river prawns at different molting stages. The identification of immune-related genes is of great significance for enhancing the immunity of the Oriental river prawn, or other crustaceans, by transgenic methods in the future.

RNA interference analysis of potential functions of cyclin A in the reproductive development of male oriental river prawns (Macrobrachium nipponense)

Cyclin A (CycA) plays essential roles in regulating multiple steps of the cell cycle, and it affects gonad development in mammals and invertebrates. Previous RNA interference (RNAi) analysis revealed that knocking-down the expression of CycA in female oriental river prawns (Macrobrachium nipponense) inhibited ovarian development. CycA was also predicted to have regulatory roles in reproductive development of male M. nipponense based on significant changes of Mn-CycA expression after eyestalk ablation. The goal of this study was to investigate the potential functions of CycA in the reproductive development of male M. nipponense using RNAi and histological observations. Quantitative real-time PCR analysis revealed that both single-side and double-side eyestalk ablation stimulated the expressions of Mn-CycA, and the expression was higher in prawns with double-side eyestalk ablation (p < 0.05). Mn-CycA expression was significantly higher in the testis and androgenic gland during the reproductive season than during the non-reproductive season (p < 0.05). In the RNAi analysis, Mn-CycA expression significantly decreased after prawns were injected with dsCycA, and the expression of insulin-like androgenic gland hormone (Mn-IAG) also decreased as Mn-CycA expression decreased. This result indicated that CycA positively regulated the expression of IAG in M. nipponense. Histological observations revealed that the number of sperm decreased dramatically to <5% of the total cells in the testis of the dsCycA-treated group compared to that of control group on day 14, indicating that knockdown of Mn-CycA expression inhibited testis development by affecting the expression of Mn-IAG in M. nipponense. These results highlighted the functions of CycA in male reproductive development of M. nipponense, which can be applied to future studies of male reproduction in other crustacean species.

RNA Interference Analysis of the Functions of Cyclin B in Male Reproductive Development of the Oriental River Prawn (Macrobrachium nipponense)

Cyclin B (CycB) plays essential roles in cell proliferation and promotes gonad development in many crustaceans. The goal of this study was to investigate the regulatory roles of this gene in the reproductive development of male oriental river prawns (Macrobrachium nipponense). A phylo-genetic tree analysis revealed that the protein sequence of Mn-CycB was most closely related to those of freshwater prawns, whereas the evolutionary distance from crabs was much longer. A quantitative PCR analysis showed that the expression of Mn-CycB was highest in the gonad of both male and female prawns compared to that in other tissues (p < 0.05), indicating that this gene may play essential roles in the regulation of both testis and ovary development in M. nipponense. In males, Mn-CycB expression in the testis and androgenic gland was higher during the reproductive season than during the non-reproductive season (p < 0.05), implying that CycB plays essential roles in the reproductive development of male M. nipponense. An RNA interference analysis revealed that the Mn-insulin-like androgenic gland hormone expression decreased as the Mn-CycB expression decreased, and that few sperm were detected 14 days after the dsCycB treatment, indicating that CycB positively affects testis development in M. nipponense. The results of this study highlight the functions of CycB in M. nipponense, and they can be applied to studies of male reproductive development in other crustacean species.

MnHR4 Functions during Molting of Macrobrachium nipponense by Regulating 20E Synthesis and Mediating 20E Signaling

HR4, a member of the nuclear receptor family, has been extensively studied in insect molting and development, but reports on crustaceans are still lacking. In the current study, the MnHR4 gene was identified in Macrobrachium nipponense. To further improve the molting molecular mechanism of M. nipponense, this study investigated whether MnHR4 functions during the molting process of M. nipponense. The domain, phylogenetic relationship and 3D structure of MnHR4 were analyzed by bioinformatics. Quantitative real-time PCR (qRT-PCR) analysis showed that MnHR4 was highly expressed in the ovary. In different embryo stages, the highest mRNA expression was observed in the cleavage stage (CS). At different individual stages, the mRNA expression of MnHR4 reached its peak on the fifteenth day after hatching (L15). The in vivo injection of 20-hydroxyecdysone (20E) can effectively promote the expression of the MnHR4 gene, and the silencing of the MnHR4 gene increased the content of 20E in M. nipponense. The regulatory role of MnHR4 in 20E synthesis and 20E signaling was further investigated by RNAi. Finally, the function of the MnHR4 gene in the molting process of M. nipponense was studied by counting the molting frequency. After knocking down MnHR4, the molting frequency of M. nipponense decreased significantly. It was proved that MnHR4 plays a pivotal role in the molting process of M. nipponense.

Selection and validation of reference genes for RT-qPCR-based analyses of Anastatus japonicus Ashmead (Hymenoptera: Helicopteridae)

RT-qPCR remains a vital approach for molecular biology studies aimed at quantifying gene expression in a range of physiological or pathological settings. However, the use of appropriate reference genes is essential to attain meaningful RT-qPCR results. Anastatus japonicus Ashmead (Hymenoptera: Helicopteridae) is an important egg parasitoid wasp and natural enemy of fruit bugs and forest caterpillars. While recent transcriptomic studies have analyzed gene expression profiles in A. japonicus specimens, offering a robust foundation for functional research focused on this parasitoid, no validated A. japonicus reference genes have yet been established, hampering further research efforts. Accordingly, this study sought to address this issue by screening for the most stable internal reference genes in A. japonicus samples to permit reliable RT-qPCR analyses. The utility of eight candidate reference genes (ACTIN, TATA, GAPDH, TUB, RPL13, RPS6, EF1α, RPS3a) was assessed under four different conditions by comparing developmental stages (larvae, pupae, adults), tissues (abdomen, chest, head), sex (male or female adults), or diapause states (diapause induction for 25, 35, 45, or 55 days, or diapause termination). RefFinder was used to calculate gene stability based on the integration of four algorithms (BestKeeper, Normfinder, geNorm, and ΔCt method) to determine the optimal RT-qPCR reference gene. Based on this approach, RPS6 and RPL13 were found to be the most reliable reference genes when assessing different stages of development, while ACTIN and EF1α were optimal when comparing adults of different sexes, RPL13 and EF1α were optimal when analyzing different tissues, and TATA and ACTIN were optimal for different diapause states. These results provide a valuable foundation for future RT-qPCR analyses of A. japonicus gene expression and function under a range of experimental conditions.

Hepatopancreas transcriptome analyses provide new insights into the molecular regulatory mechanism of fast ovary maturation in Macrobrachium nipponense

Background

Macrobrachium nipponense is an economically and ecologically important freshwater prawn that is widely farmed in China. In contrast to other species of marine shrimp, M. nipponense has a short sexual maturity period, resulting in not only high stocking densities, but also a reduced survival rate and increased risk of hypoxia. Therefore, there is an urgent need to study the molecular mechanisms underlying fast ovary maturation in this species.

Results

Comparative transcriptome analysis was performed using hepatopancreatic tissue from female M. nipponense across five ovarian maturation stages to explore differentially expressed genes and pathways involved in ovarian maturation. In total, 118.01 Gb of data were generated from 15 transcriptomes. Approximately 90.46% of clean reads were mapped from the M. nipponense reference genome. A comprehensive comparative analysis between successive ovarian maturation stages generated 230–5814 differentially expressed genes. Gene Ontology (GO) enrichment was highly concentrated in the “biological process” category in all four comparison groups, and mainly focused on energy synthesis and accumulation, energy decomposition and transport. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results showed that, among 20 significantly enriched KEGG pathways, nine were involved in the synthesis, degradation, and metabolism of carbohydrates, lipids, and other nutrient intermediates, suggesting that the hepatopancreas has an important role in energy supply during ovarian maturation. Furthermore, the “Insect hormone biosynthesis” pathway was found to have a dominant role in the development of the ovary from immaturity to maturity, supporting the hypothesis that ecdysteroid- and juvenile hormone-signaling pathways have an important role in hepatopancreas regulation of ovarian maturation.

Conclusion

Taken together, this study sheds light on the role of the hepatopancreas in the molecular regulation of ovary maturation in M. nipponense. The present study provided new insights for understanding the mechanisms of reproductive regulation in crustaceans.

Cathepsin D Plays a Vital Role in Macrobrachium nipponense of Ovary Maturation: Identification, Characterization, and Function Analysis

The oriental river prawn Macrobrachium nipponense is an economically important aquacultural species. However, its aquaculture is negatively impacted by the rapid sexual maturation of female M. nipponense. The fast sexual maturation produces a large number of offspring which leads to a reduction in resilience, a low survival rate, and an increased risk of hypoxia, this in turn, seriously affects the economic benefits of prawn farming. Cathepsin D is a lysosomal protease involved in the ovarian maturation of M. nipponense. In the current study, the cDNA of the gene encoding cathepsin D (Mn-CTSD) was cloned from M. nipponense. The total length was 2391 bp and consisted of an open reading frame (ORF) of 1158 bp encoding 385 amino acids. Sequence analysis confirmed the presence of conserved N-glycosylation sites and characteristic sequences of nondigestive cathepsin D. The qPCR analysis indicated that Mn-CTSD was highly expressed in all tissues tested, most significantly in the ovaries, whereas in situ hybridization showed that expression occurred mainly in oocyte nuclei. Analysis of its expression during development showed that Mn-CTSD peaked during the O-IV stage of ovarian maturation. For the RNAi interference experiment, female M. nipponense specimens in the ovary stage I were selected. Injection of Mn-CTSD double-stranded (ds)RNA into female M. nipponense decreased the expression of Mn-CTSD in the ovaries, such that the Gonad Somatic Index (GSI) of the experimental group was significantly lower than that of the control group (1.79% versus 4.57%; p < 0.05). Ovary development reached the O-III stage in 80% of the control group, compared with 0% in the experimental group. These results suggest that Mn-CTSD dsRNA inhibits ovarian maturation in M. nipponense, highlighting its important role in ovarian maturation in this species and suggesting an approach to controlling ovarian maturation during M. nipponense aquaculture.

Reference Genes Selection of Gymnosporangium yamadae during the Interaction with Apple Leaves

Apple rust disease caused by Gymnosporangium yamadae is the one of the major threats to the development of the apple industry in China, but the pathogenic molecular mechanism of the disease remains unclear. It is imperative to screen out appropriate reference genes during the interaction between G. yamadae and apple leaves to analyze the gene expression patterns during the pathogenesis of G. yamadae. ACT, EF1, EF2, GAPDH, 40S, 60S, α-TUB, β-TUB and UBCE3 were selected as candidate reference genes based on the transcriptomic dataset of G. yamadae. The expression levels were tested by real-time quantitative PCR during time-course infection of apple leaves and the expression stabilities were evaluated by △Ct method as well as by three software (NormFinder, geNorm and BestKeeper) and one web-based analysis software (RefFinder). The expression stability of the candidate reference genes was further validated by using the effector candidate gene Cluster-3395.48660 as the target gene in RT-qPCR. According to the results by △Ct and BestKeeper, 40S, EF2 and EF1 were the most stable reference genes, while EF1, EF2 and GAPDH were the most stable reference genes based on the NormFinder analysis result. The geNorm recommended the most stable genes EF1, EF2 and α-TUB as reference genes. Comprehensive analysis results of the RefFinder indicated EF1, EF2 and α-TUB were the most suitable genes. Based on these results, EF1, EF2 and α-TUB were considered as reference genes for analyzing the gene expression profiles of Cluster-3395.48660 in different infection stages, and the results were consistent with the transcriptome data. All the results suggest that the combination of EF1, EF2 and α-TUB proved to be acceptable reference genes during the interaction between G. yamadae and apple leaves.

Selection of reference genes for quantitative real-time PCR normalization in the coffee white stem borer, Xylotrechus quadripes Chevrolat (Coleoptera: Cerambycidae)

The coffee white stem borer, Xylotrechus quadripes Chevrolat (Coleoptera: Cerambycidae), is a major destructive pest of Coffea arabica L. (Gentianales: Rubiaceae), widely planted in many Asian countries, including China. Quantitative real-time polymerase chain reaction (qRT-PCR) is a common method for quantitative analysis of gene transcription levels. To obtain accurate and reliable qRT-PCR results, it is necessary to select suitable reference genes to different experimental conditions for normalizing the target gene expression. However, the stability of the expression of reference genes in X. quadripes has rarely been studied. In this study, the expression stability of nine candidate reference genes were investigated under biotic and abiotic conditions for use in qRT-PCR's normalization. By integrating the results of four algorithms of NormFinder, BestKeeper, geNorm, and RefFinder, the optimal reference gene combinations in different experimental conditions were performed as follows: RPL10a and EIF3D were the optimal reference genes for developmental stage samples, EIF4E, RPL10a, and RPS27a for tissue samples, V-ATP and EF1α for the sex samples, EIF3D and V-ATP for temperature treatment, RPS27a and RPL10a for insecticide stress, and RPL10a, RPS27a, and EF1α for all the samples. This study will help to obtain the stable internal reference genes under biotic and abiotic conditions and lay the foundation for in-depth functional research of target genes or genomics on olfactory molecular mechanisms, temperature adaptability, and insecticide resistance in X. quadripes.

Selection and Validation of Reference Genes for Quantitative Real-Time PCR Normalization in Athetis dissimilis (Lepidoptera: Noctuidae) Under Different Conditions

Reference genes are the key to study gene expression patterns using quantitative real-time PCR (qRT-PCR). No studies on the reference genes of Athetis dissimilis, an important agricultural pest, have been reported. In order to determine the reference genes for qRT-PCR normalization in A. dissimilis under different conditions, 10 candidate genes [18S ribosomal protein (18S), 28S ribosomal protein (28S), arginine kinase (AK), elongation factor 1 alpha (EF1-α), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L32 (RPL32), ribosomal protein L40 (RPL40), alpha-tubulin (α-TUB), beta-actin (β-ACT), and beta-tubulin (β-TUB)] of A. dissimilis were selected to evaluate their stability as reference genes under different biotic and abiotic conditions by using five tools, geNorm, NormFinder, BestKeeper, ΔCt, and RefFinder. Furthermore, CSP1 and superoxide dismutase (SOD) were used as target genes to validate the candidate reference genes. The results showed that different reference genes were needed under different experimental conditions, among which, EF-1α, RPL40, and 18S are most suitable reference genes for studying genes related development stages of A. dissimilis, RPL40 and α-TUB for larval tissues, α-TUB and 28S for adult tissues, EF-1α and β-ACT for insecticidal treatments, β-ACT and 28S for temperature treatments, EF-1α and β-ACT for starvation treatments, RPL40 and 18S for dietary treatments, and 18S, 28S, and α-TUB for all the samples. These results provide suitable reference genes for studying gene expression in A. dissimilis under different experimental conditions, and also lay the foundation for further research into the function of related genes in A. dissimilis.

RNA Interference Analysis Reveals the Positive Regulatory Role of Ferritin in Testis Development in the Oriental River Prawn, Macrobrachium nipponense

Ferritin plays an essential role in organismic and cellular iron homeostasis in Macrobrachium nipponense. In this study, we aimed to investigate the role of ferritin in the sexual development of male M. nipponense. According to the qPCR analysis of different tissues and developmental stages, ferritin exhibited high expression levels in the testis and androgenic gland, from post-larval developmental stage 5 (PL5) to PL15, indicating that it may be involved in gonad differentiation and development, especially in male sexual development. In situ hybridization and qPCR analysis in various reproductive cycles of the testis indicated that ferritin may play an essential role in spermatogonia development in M. nipponense. RNAi analysis revealed that ferritin positively affected mRNA expression of the insulin-like androgenic gland (Mn-IAG) and the secretion of testosterone, and thus positively affected testis development in M. nipponense. This study highlighted the functions of ferritin in the sexual development of male M. nipponense and provided important information for the establishment of a technique to regulate the process of testis development in M. nipponense.

Identification of Important Genes Involved in the Sex-Differentiation Mechanism of Oriental River Prawn, Macrobrachium nipponense, During the Gonad Differentiation and Development Period

Identification of important genes, involved in the gonad differentiation and development, plays essential roles in the establishment of the artificial technique to regulate the process of testis development in M. nipponense. In this study, we aimed to determine the sensitive period of gonad differentiation and development through hematoxylin and eosin (HE) staining. The important genes involved in the gonad differentiation and development of M. nipponense were then identified through transcriptome profiling analysis during the sensitive period of gonad differentiation and development. HE staining analysis revealed that the sensitive period of gonad differentiation and development was from the post–larval developmental stages 5 (PL5) to PL25, which was dramatically faster than was for the other identified aquatic animals. The transcriptome profiling analysis predicted that phagosome, lysosome, oxidative phosphorylation, and glycolysis/gluconeogenesis play essential roles in the mechanism of gonad differentiation and development in M. nipponense. A total of 29 genes were further identified as the candidate genes, involved in the process of gonad differentiation and development in M. nipponense, based on the gene annotation and gene expression pattern. The qPCR analysis of Mn-JHEH, Mn-DHP, Mn-ALY, and Mn-SMA6 during the whole developmental process revealed that all of these four genes showed high expression levels during the sensitive period of gonad differentiation and development in M. nipponense. Mn-JHEH, Mn-DHP, and Mn-ALY showed higher expressions at PL25F than at PL25M, while Mn-SMA6 showed a higher expression at PL25M. The RNA interference (RNAi) analysis was further used to investigate the potential functions of SMA6 in male sexual development of M. nipponense. The RNAi analysis revealed that SMA6 positively regulated the testis development in M. nipponense by affecting the expression of Mn-IAG. This study provided valuable evidences for the establishment of the technique to regulate the process of gonad development in M. nipponense.

Evaluation of Optimal Reference Genes for qRT-PCR Analysis in Hyphantria cunea (Drury)

The relative quantification of gene expression is mainly achieved through reverse transcription-quantitative PCR (qRT-PCR); however, its reliability and precision rely on proper data normalization using one or more optimal reference genes. Hyphantria cunea (Drury) has been an invasive pest of forest trees, ornamental plants, and fruit trees in China for many years. Currently, the molecular physiological role of reference genes in H. cunea is unclear, which hinders functional gene study. Therefore, eight common reference genes, RPS26, RPL13, UBI, AK, RPS15, EIF4A, β-actin, α-tub, were selected to evaluate levels of gene expression stability when subjected to varied experimental conditions, including developmental stage and gender, different tissues, larvae reared on different hosts and different larval density. The geNorm, BestKeeper, ΔCt method, and NormFinder statistical algorithms were used to normalize gene transcription data. Furthermore, the stability/suitability of these candidates was ranked overall by RefFinder. This study provides a comprehensive evaluation of reference genes in H. cunea and could help select reference genes for other Lepidoptera species.

Identification of potential functions of polo-like kinase 1 in male reproductive development of the oriental river prawn (Macrobrachium nipponense) by RNA interference analysis

Polo-like kinase 1 (Plk1) has multiple functions in the cell cycle, including in the maturation of centrosomes during the G2/M transition, the separation of centrosomes, and the activation of cyclin-dependent kinase 1 expression and spindle assembly. In this study, we investigated the potential regulatory roles of Plk1 in the reproductive development of the male oriental river prawn (Machrobrachium nipponense). The full cDNA sequence of Mn-Plk1 was 2360 base pairs long, with an open reading frame of 1836 base pairs encoding 611 amino acids. Protein sequence alignment identified a conserved serine/threonine kinase domain and two Polo-boxes. Phylogenetic tree analysis revealed that Mn-Plk1 had the closest evolutionary distance with Plk1s of freshwater prawns and then with those of crustacean species, whereas the evolutionary distance with mollusks was much more distant. Quantitative PCR analysis predicted that Mn-Plk1 plays essential roles in the regulation of gonad development. RNA interference analysis and histological observations showed that expression of insulin-like androgenic gland hormone decreased as the expression of Mn-Plk1 decreased, and fewer than 5% of cells were sperm cells at day 14 in the dsPlk1 injected prawns. This result indicated that Plk1 positively regulated testis development in M. nipponense by affecting the expression of this hormone. Our results highlight the functions of Plk1 in M. nipponense and provide valuable information that can be applied to establish artificial techniques to regulate testis development in this species.

A Novel Legumain-Like Protease in Macrobrachium nipponense: Identification, Characterization, and Function Analysis in Ovary Maturation

Legumain, also called aspartic endopeptidase (AEP), is a member of the cysteine protease family and is involved in various physiological processes. In this study, we analyzed the characteristics of a novel legumain-like (named Mn-Lel) in the female oriental river prawn, Macrobrachium nipponense, which is involved in ovary maturation. The Mn-Lel is 1,454 bp in length, including a 1,290-bp open reading frame that encodes 430 amino acids. qPCR analysis indicated that Mn-Lel is specifically highly expressed in the hepatopancreas and ovaries of female prawns. It is rarely expressed in embryogenesis, weakly expressed in early larval development stages, and then significantly increased after metamorphosis, which indicated that Mn-Lel is not a maternal gene and mainly plays a role in adults. During the different ovarian stages, Mn-Lel expression in the hepatopancreas had no obvious rules, while its expression in the ovaries had a significant peak in stage III. In situ hybridization studies revealed that Mn-Lel is localized in the oocyte of the ovary. Changes in the gonadosomatic index confirmed the inhibitory effects of Mn-Lel dsRNA on ovary maturation. These results suggest that Mn-Lel has a key role in promoting ovary maturation.

Integrated Metabolomics and Transcriptomic Analysis of Hepatopancreas in Different Living Status Macrobrachium nipponense in Response to Hypoxia

As the basic element of aerobic animal life, oxygen participates in most physiological activities of animals. Hypoxia stress is often the subject of aquatic animal research. Macrobrachium nipponense, an economically important aquatic animal in southern China, has been affected by hypoxia for many years and this has resulted in a large amount of economic loss due to its sensitivity to hypoxia; Metabolism and transcriptome data were combined in the analysis of the hepatopancreas of M. nipponense in different physiological states under hypoxia; A total of 108, 86, and 48 differentially expressed metabolites (DEMs) were found in three different comparisons (survived, moribund, and dead shrimps), respectively. Thirty-two common DEMs were found by comparing the different physiological states of M. nipponense with the control group in response to hypoxia. Twelve hypoxia-related genes were identified by screening and analyzing common DEMs. GTP phosphoenolpyruvate carboxykinase (PEPCK) was the only differentially expressed gene that ranked highly in transcriptome analysis combined with metabolome analysis. PEPCK ranked highly both in transcriptome analysis and in combination with metabolism analysis; therefore, it was considered to have an important role in hypoxic response. This manuscript fills the one-sidedness of the gap in hypoxia transcriptome analysis and reversely deduces several new genes related to hypoxia from metabolites. This study contributes to the clarification of the molecular process associated with M. nipponense under hypoxic stress.

Public transcriptome database-based selection and validation of reliable reference genes for breast cancer research

Background

Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) is the most sensitive technique for evaluating gene expression levels. Choosing appropriate reference genes (RGs) is critical for normalizing and evaluating changes in the expression of target genes. However, uniform and reliable RGs for breast cancer research have not been identified, limiting the value of target gene expression studies. Here, we aimed to identify reliable and accurate RGs for breast cancer tissues and cell lines using the RNA-seq dataset.

Methods

First, we compiled the transcriptome profiling data from the TCGA database involving 1217 samples to identify novel RGs. Next, ten genes with relatively stable expression levels were chosen as novel candidate RGs, together with six conventional RGs. To determine and validate the optimal RGs we performed qRT-PCR experiments on 87 samples from 11 types of surgically excised breast tumor specimens (n = 66) and seven breast cancer cell lines (n = 21). Five publicly available algorithms (geNorm, NormFinder, ΔCt method, BestKeeper, and ComprFinder) were used to assess the expression stability of each RG across all breast cancer tissues and cell lines.

Results

Our results show that RG combinations SF1 + TRA2B + THRAP3 and THRAP3 + RHOA + QRICH1 showed stable expression in breast cancer tissues and cell lines, respectively, and that they displayed good interchangeability. We propose that these combinations are optimal triplet RGs for breast cancer research.

Conclusions

In summary, we identified novel and reliable RG combinations for breast cancer research based on a public RNA-seq dataset. Our results lay a solid foundation for the accurate normalization of qRT-PCR results across different breast cancer tissues and cells.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12938-021-00963-8.

Identification of candidate genes from androgenic gland in Macrobrachium nipponense regulated by eyestalk ablation

The eyestalk of crustaceans, such as Macrobrachium nipponense, contains many neurosecretory hormones affecting the process of reproduction, molting, metabolism of glucose, and other functions. In this study, important metabolic pathways and candidate genes involved in male sexual development were selected from M. nipponense. The methodology involved performing long-read and next generation transcriptome sequencing of genes from the androgenic gland after eyestalk ablation. qPCR analysis revealed that the mRNA expression of Mn-IAG was significantly increased after ablation of both the single-side (SS) and double-side (DS) eyestalk, compared with the control group (CG). The long-read transcriptome generated 49,840 non-redundant transcripts. A total of 1319, 2092 and 4351 differentially expressed genes (DEGs) were identified between CG versus SS, SS versus DS and CG versus DS, respectively. These data indicated that ablation of the double-sided eyestalk played stronger regulatory roles than the single-side ablation on male sexual development in M. nipponense. This was consistent with the qPCR analysis. Cell Cycle, Cellular Senescence, Oxidative Phosphorylation, Glycolysis/Gluconeogenesis and Steroid Hormone Biosynthesis were the primary enriched metabolic pathways in all three comparisons, and the important genes from these metabolic pathways were also selected. qPCR permitted secondary confirmation of ten DEGs identified through RNA-seq. RNAi-mediated silencing analyses of Hydroxysteroid dehydrogenase like 1 (HSDL1) revealed that HSDL1 has a positive regulatory effect on testes development. This study provides valuable insight into male sexual development in M. nipponense, including metabolic pathways and genes, paving the way for advanced studies on male sexual development in this species and in other crustaceans.

Transcriptome analysis of hepatopancreas from different living states oriental river prawn (Macrobrachium nipponense) in response to hypoxia

As an important economical freshwater prawn, Macrobrachium nipponense has difficulty with adapting to hypoxia. In this study, comparative transcriptome analysis was used for the first time to explore the differences between different living states of Macrobrachium nipponense under hypoxia. A total of 94.22 Gb clean reads were obtained and assembled into 54,688 unigenes. A total of 224, 266, and 750 differently expressed genes were found in the comparison of the control and death groups, the control and moribund groups, and the control and survived groups, respectively. Three signal pathways closely related to hypoxia were found by enriching of the signal pathways in three comparison groups. In addition, much attention was focused on the differential genes in these pathways. Oxidative stress related genes, such as 70 kDa heat shock protein, phosphoenolpyruvate carboxykinase and cyclooxygenase were differentially expressed in different comparisons. After comparing with previous studies, cyclooxygenase was found to be an important hypoxia-related gene that is fully involved in the hypoxic response. Interestingly, two new genes with no Nr annotation were found in this manuscript. This manuscript will enrich our understanding of oxidative stress response to hypoxia and provide a theoretical basis for the subsequent solution of apoptosis caused by hypoxia.

Identification and Characterization of the Succinate Dehydrogenase Complex Iron Sulfur Subunit B Gene in the Oriental River Prawn, Macrobrachium nipponense

Previous studies have revealed that SDHB has potential functions in the male sexual differentiation and development in M. nipponense through providing ATP. In this study, the functions of Mn-SDHB were further analyzed in depth using quantitative polymerase chain reaction (qPCR), in situ hybridization, western-blot, and RNA interference (RNAi), combined with the histological observations. The full-genome sequence of Mn-SDHB was 54,608 bp at Chromosome 34, including 7 introns and 6 exons. The full-length cDNA sequence of Mn-SDHB was 1,268 base pairs (bp) long with an open reading frame of 807 bp, encoding for 268 amino acids. The highest expression level of Mn-SDHB in different tissues was observed in the testis, and male prawns at post-larval developmental stage 25 during different developmental stages, indicating that SDHB was potentially involved in the male sexual development in M. nipponense. In situ hybridization and western-blot analysis indicated that SDHB plays essential roles in the testis development. The in situ hybridization analysis also implies the potential roles of Mn-SDHB in ovarian development. The expressions of Mn-IAG were decreased after Mn-SDHB dsRNA injection, indicating SDHB has the positive regulatory effects on IAG in M. nipponese. Thus, SDHB was involved in the mechanism of the male sexual development. The testis development was inhibited, and sperms were rarely observed after 10 days of Mn-SDHB dsRNA injection, indicating SDHB has positive effects on the male sexual development in M. nipponense. This study highlights the functions of SDHB in M. nipponense, which provide new insights for the future studies of the male sexual development in other crustacean species.

Transcriptome analysis of five ovarian stages reveals gonad maturation in female Macrobrachium nipponense

Background

Macrobrachium nipponense is an economically important species of freshwater shrimp in China. Unlike other marine shrimps, the ovaries in adult female M. nipponense can mature rapidly and periodically during the reproductive period, but the resulting high stocking densities and environmental deterioration can negatively impact the harvest yield and economic benefits. To better understand ovary development in female M. nipponense, we performed systematic transcriptome sequencing of five different stages of ovarian maturation.

Results

We obtained 255,966 Gb of high quality transcriptome data from 15 samples. Of the 105,082 unigenes that were selected, 30,878 were successfully annotated. From these unigenes, we identified 17 differentially expressed genes and identified three distinct gene expression patterns related to different biological processes. We found that cathepins, legumains, and cystatin were enriched in the lysosome pathway, and they are related to vitellogenin hydrolysis. Additionally, we found that myosin heavy chain 67 participated in oocyte excretion.

Conclusions

We provide the first detailed transcriptome data relating to the ovarian maturation cycle in M. nipponense. Our results provide important reference information about the genomics, molecular biology, physiology, and population genetics of M. nipponense and other crustaceans. It is conducive to further solve the problem of M. nipponense rapid ovarian maturation from the aspects of energy supply and cell division.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07737-5.

Selection and Validation of Reference Genes for Gene Expression Analysis in Tuta absoluta Meyrick (Lepidoptera: Gelechiidae)

Simple Summary

Reference genes are critical for standardizing expression data of RT-qPCR across samples of organisms under different experimental conditions. However, most commonly used reference genes may not be stably expressed leading to a risk of misinterpretation of results. In our study, nine reference genes were evaluated in Tuta absoluta (a destructive pest of tomato) at different developmental stages, tissues, 20-hydroxyecdysone (20E) and insecticide treatments. Finally, the expression profile of indicator gene EcR after 20E treatment was evaluated to verify the accuracy of the results. This study is essential for improving accuracy and reliability to normalize gene expression data in T. absoluta and provides a useful strategy for other insects.

Abstract

The tomato leaf miner, Tuta absoluta is a destructive pest of tomato. The leaf-mining activities of its larvae can cause significant yield losses. Real-time quantitative polymerase chain reaction (RT-qPCR) is commonly used to measure gene expression, and the selection of stable reference genes for calibration and standardization is critical for accurate use of RT-qPCR. We studied the stable expression of nine common housekeeping genes in T. absoluta. These were examined at different developmental stages, in larval tissues, as well as those induced by exposure to 20E and insecticides. Four dedicated algorithms (geNorm, BestKeeper, NormFinder, and ΔCt method) and online tool (RefFinder) were used to analyze and rank the tested reference genes. Based on the standardized gene expression data of target gene ecdysone receptor (EcR), the applicability of specific reference genes was verified. The results clarify that the optimal internal reference genes vary greatly under different experimental conditions. GAPDH and RPS11 were the best reference genes for developmental stages; RPL28 and RPL10 for different tissues; EF1α and RPL28 for 20E treatment; EF1α and RPL7A for insecticide treatments. The most suitable reference genes in all experimental conditions are EF1α and RPL28.

Transcriptome Profiling Analysis of the Testis After Eyestalk Ablation for Selection of the Candidate Genes Involved in the Male Sexual Development in Macrobrachium nipponense

The eyestalk of crustacean species secretes many hormones, affecting the process of reproduction, molting, metabolism of glucose, and other functions in crustaceans. In this study, important metabolic pathways and candidate genes involved in the male sexual development were identified through performing the transcriptome profiling analysis of the testis after the ablation of eyestalk from Macrobrachium nipponense. The histological observations revealed that the testis development became vigorous after eyestalk ablation, indicating that the hormones secreted by the eyestalk have negative effects on the testis development in M. nipponense. Transcriptome profiling analysis revealed that 1,039, 1,226, and 3,682 differentially expressed genes (DEGs) were identified between normal prawns (CG) vs single-side eyestalk ablation prawns (SS), SS vs double-side eyestalk ablation prawns (DS), and CG vs DS, respectively, indicating that the ablation of double-side eyestalk has more significant regulatory roles on male sexual development than that of single-side ablation, which was consistent with the histological observations. Lysosome, Apoptosis, Glycolysis/Gluconeogenesis, and Insulin signaling pathway were the main enriched metabolic pathways in all of these three comparisons, and the important genes from these metabolic pathways were also selected. The qPCR verifications of 10 DEGs from these metabolic pathways were the same as those of RNA-seq. The qPCR, in situ hybridization, and RNA interference analysis of Mn-NF_kBα revealed that NF_kBα has a positive regulatory effect on testis development. This study provided new insights on male sexual development in M. nipponense, promoting the studies on male sexual development in other crustaceans as well.

Identification and Evaluation of Suitable Reference Genes for RT-qPCR Analysis in Hippodamia variegata (Coleoptera: Coccinellidae) Under Different Biotic and Abiotic Conditions

Reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR) is an accurate and convenient technique for quantifying expression levels of the target genes. Selection of the appropriate reference gene is of the vital importance for RT-qPCR analysis. Hippodamia variegata is one of the most important predatory natural enemies of aphids. Recently, transcriptome and genome sequencings of H. variegata facilitate the gene functional studies. However, there has been rare investigation on the detection of stably expressed reference genes in H. variegata. In the current study, by using five analytical tools (Delta Ct, geNorm, NormFinder, BestKeeper, and RefFinder), eight candidate reference genes, namely, Actin, EF1α, RPL7, RPL18, RPS23, Tubulin-α, Tubulin-β, and TufA, were evaluated under four experimental conditions including developmental stages, tissues, temperatures, and diets. As a result, a specific set of reference genes were recommended for each experimental condition. These findings will help to improve the accuracy and reliability of RT-qPCR data, and lay a foundation for further exploration on the gene function of H. variegata.

Comparative transcriptome analysis of lethality in response to RNA interference of the oriental river prawn (Macrobrachium nipponense)

A previous study identified slow-tonic S2 tropomyosin and slow tropomyosin isoform as sex-related genes in Macrobrachium nipponense. Their functions were analyzed using RNA interference. However, more than half of the specimens died approximately 8-12 h after injection of the respective double-stranded RNAs (dsRNAs), and HE staining indicated that the heart and gills were the most likely tissues responsible for the resultant deaths. In the current study, we conducted a comparative transcriptomic study of the gills and hearts of M. nipponense to identify potential target genes associated with acute death after dsRNA injection. A total of 68,772 annotated unigenes were generated. In the heart, differentially expressed genes (DEGs) were mainly enriched in glycolysis/gluconeogenesis and oxidative phosphorylation, while the most relevant pathways in the gills were lysosome, phagosome, and peroxisome. Ten DEGs were screened out and analyzed under lethal hypoxic stress. Among these, fructose 1, 6-biphosphate-aldolase (FBA), glyceraldehyde 3-phosphate dehydrogenase (GDPDH), alcohol dehydrogenase class-3 (ADC3), ATP-synthase subunit 9 (ATPS9), and acid ceramidase-like (ACL) were all differentially expressed under hypoxic conditions. This study shed light on the lethal mechanism caused by interference with tropomyosin genes in M. nipponense, and identifies the related pathways and key genes that could help to improve stress resistance and tolerance in M. nipponense.

MnFtz-f1 Is Required for Molting and Ovulation of the Oriental River Prawn Macrobrachium nipponense

Molting and ovulation are the basic processes responsible for the growth and reproduction of Macrobrachium nipponense; however, the molecular mechanisms of molting and ovulation in M. nipponense are poorly understood. The present study aimed to use MnFtz-f1 as the starting point to study the molting and ovulation phenomena in M. nipponense at the molecular level. The full-length MnFtz-f1 cDNA sequence was 2,198 base pairs (bp) in length with an open reading frame of 1,899 bp encoding 632 amino acids. Quantitative real-time PCR analysis showed that MnFtz-f1 was highly expressed in the ovary at the cleavage stage and on the fifth day after hatching. In vivo administration of 20-hydroxyecdysone (20E) showed that 20E effectively inhibited the expression of the MnFtz-f1 gene, and the silencing of the MnFtz-f1 gene reduced the content of 20E in the ovary. In situ hybridization (ISH) analysis revealed the localization of MnFtz-f1 in the ovary. Silencing of MnFtz-f1 by RNA interference (RNAi) resulted in significant inhibition of the expression of the vitellogenin (Vg), Spook, and Phantom genes, thus confirming that MnFtz-f1 had a mutual regulatory relationship with Vg, Spook, and Phantom. After RNAi, the molting frequency and ovulation number of M. nipponense decreased significantly, which demonstrated that MnFtz-f1 played a pivotal role in the process of molting and ovulation.

Molecular Characterization of a Novel Cathepsin L in Macrobrachium nipponense and Its Function in Ovary Maturation

Cathepsin L genes, which belonged to cysteine proteases, were a series of multifunctional protease and played important roles in a lot of pathological and physiological processes. In this study, we analyzed the characteristics a cathepsin L (named Mn-CL2) in the female oriental river prawn, Macrobrachium nipponense which was involved in ovary maturation. The Mn-CL2 was1,582 bp in length, including a 978 bp open reading frame that encoded 326 amino acids. The Mn-CL2 was classified into the cathepsin L group by phylogenetic analysis. Real-time PCR (qPCR) analysis indicated that Mn-CL2 was highly expressed in the hepatopancreas and ovaries of female prawns. During the different ovarian stages, Mn-CL2 expression in the hepatopancreas and ovaries peaked before ovarian maturation. In situ hybridization studies revealed that Mn-CL2 was localized in the oocyte of the ovary. Injection of Mn-CL2 dsRNA significantly reduced the expression of vitellogenin. Changes in the gonad somatic index also confirmed the inhibitory effects of Mn-CL2 dsRNA on ovary maturation. These results suggest that Mn-CL2 has a key role in promoting ovary maturation.

Three neuroparsin genes from oriental river prawn, Macrobrachium nipponense, involved in ovary maturation

In this study, we identified three neuroparsin (NP) genes in Macrobrachium nipponense: Mn-NP1, Mn-NP2, and Mn-NP3, encoding 99, 100, and 101 amino acid proteins, respectively. Multiple sequence alignments showed that these genes contained 12 cysteine residues, of which 11 were at conserved positions. The total sequence identity between the genes was 47.5%, and they showed a high degree of sequence identity (> 54% similarity) with other crustacean genes. Phylogenetic tree analysis showed that Mn-NPs were clustered at different branches, indicating that Mn-NPs may have different functions. Tissue distribution data revealed that the three genes were present in males and females during the breeding and nonbreeding season, but their expression patterns differed. Mn-NP1 was highly expressed in the breeding season, in the male testis, and highly expressed in the nonbreeding season, in the female ovary. Mn-NP3 exhibited biased female expression in the breeding and nonbreeding season, with dominant expression in the ovary. All Mn-NPs were detected during embryo development, but with different expression patterns. These data indicated that Mn-NP1 may function during embryonic development, and that Mn-NP2 may be expressed during early embryo cell division, and late larval development. Mn-NP3 expression patterns reflected maternal inheritance, and may be associated with ovarian maturation. These expression data suggested that Mn-NP1 and Mn-NP2 are negatively correlated with ovarian development, with inhibition roles during this development. Mn-NP3 may be involved in vitellogenesis.

Reference gene and tropomyosin expression in mud crab Scylla olivacea, Scylla paramamosain and Scylla tranquebarica

Tropomyosin, a muscle tissue protein is a major allergen in most of shellfish including mud crab. Quantitative real time-PCR (qRT-PCR) using a stable reference gene is the most sensitive approach to produce accurate relative gene expression that has yet to be demonstrated for allergenic tropomyosin in mud crab species. This study was conducted to identify the suitable reference gene and tropomyosin expression in different body parts of local mud crabs, Scylla olivacea, Scylla paramamosain and Scylla tranquebarica. Myosin, 18S rRNA, GADPH and EF1α were selected as candidate reference genes and their expression was measured in the abdomen, walking leg and cheliped tissues of local Scylla spp. The expression stability was analyzed using the comparative delta-Ct method, BestKeeper, NormFinder and geNorm then comprehensively ranked by RefFinder algorithm. Findings showed that EF1α was the most suitable reference gene across three mud crab species. Meanwhile, the abdomen, walking leg and cheliped selected their own suitable reference gene either Myosin, 18S rRNA, EF1α or GADPH. Overall, tropomyosin was the highest in S. tranquebarica, whereas the least was in S. paramamosain. Interestingly, tropomyosin was the highest in the abdomen of all mud crab species. This is the first analysis on reference genes selection for qRT-PCR data normalization of tropomyosin expression in mud crab. These results will provide more accurate findings for further gene expression and allergen analysis in Scylla spp.

Molecular cloning, characterization and functional analysis of two neuropeptide F genes from the oriental river prawn (Macrobrachium nipponense)

In invertebrates, neuropeptide F (NPF) has many regulatory functions, similar to NPY, its homologous peptide. In this study, two neuropeptide F genes were identified in Macrobrachium nipponense: Mn-NPF1 and Mn-NPF2. Mn-NPF2 shared the same amino acid sequence with Mn-NPF1, except for a 37 amino acid insert in the middle of the NPF region. The quantitative-PCR (qPCR) results indicated that Mn-NPF1 expression was positively correlated with ovarian maturation, whereas Mn-NPF2 had opposing expression patterns. Both Mn-NPFs were poorly expressed at early embryonic stages, but enhanced expression levels were observed up to day 10 after hatching, when the gonads began to differentiate. Ovary in situ hybridization (ISH) analyses showed that both Mn-NPFs were present at all stages, but were differentially localized to distinct compartments. Temperature gradient studies showed that both Mn-NPFs were implicated in the seasonal regulation of reproduction. A double-stranded (ds) RNA-Mn-NPF2 injection led to a significant 38.5% increase in the vitellogenin (VG) transcript (P < 0.05). These results demonstrated that Mn-NPF2 plays an important role in inhibiting ovarian maturation.

Validation of Real-time PCR Reference Genes of Muscle Metabolism in Harvested Spiny-Cheek Crayfish (Faxonius limosus) Exposed to Seasonal Variation

Real-time quantitative reverse transcription PCR (RT-qPCR) is a sensitive and broadly used technique of assessing gene activity. To obtain a reliable result, stably expressed reference genes are essential for normalization of transcripts in various samples. To our knowledge, this is the first systematic analysis of reference genes for normalization of RT-qPCR data in spiny-cheek crayfish (Faxonius limosus). In this study, expression of five candidate reference genes (actb, β-actin; gapdh, glyceraldehyde-3-phosphate dehydrogenase; eif, eukaryotic translation initiation factor 5a; ef-1α, elongation factor-1α; and tub, α-tubulin) in muscle samples from male and female F. limosus in spring and autumn was analyzed. Additionally, the most stable reference genes were used for accurate normalization of five target genes, i.e., tnnc, troponin c; ak, arginine kinase; fr, ferritin; ccbp-23, crustacean calcium-binding protein 23; and actinsk8, skeletal muscle actin 8. Results obtained using the geNorm and NormFinder algorithms showed high consistency, and differences in the activity of the selected actb with eif genes were successfully identified. The spring and autumn activities of the target genes (except ak) in the muscle tissue of males and females differed significantly, showing that both sexes are immensely involved in an array of breeding behaviors in spring, and females intensively recover in the autumn season. Characterization of first reference genes in spiny-cheek crayfish will facilitate more accurate and reliable expression studies in this key species.

Molecular cloning, expression, and in situ hybridization analysis of MnGPx-3 and MnGPx-4 from oriental river prawn, Macrobrachium nipponense, in response to hypoxia and reoxygenation

Glutathione peroxidase (GPx) has been the focus of increased research because of its important role as an antioxidant and in reactive oxygen species (ROS) induced damage repair. Studies on GPxs have relevance with Macrobrachium nipponense because it has poor tolerance to hypoxia in Macrobrachium nipponense. The two subunits named as MnGPx-3 and MnGPx-4 according to the glutathione peroxidase nomenclature system. Both full-length cDNAs were cloned from the hepatopancreas. In this study, we analyzed the expression of two GPxs in Macrobrachium nipponense in response to changes in environmental oxygen. Expression levels of MnGPx-3 and MnGPx-4 indicated that both have strong responses to hypoxia. In situ hybridization showed that MnGPx-3 and MnGPx-4 were located in secretory and storage cells in hepatopancreas. These results suggest that GPx gene is expressed and released by secretory cells and released response to hypoxia. In the gill tissue, however, GPxs are located in blood cells, suggesting that they perform different functions in different tissues or organs. The results of in situ hybridization were consistent with those of quantitative Real-time PCR. This study provides a basis for understanding the oxidative stress response in M. nipponense under hypoxia.

Identification and evaluation of suitable reference genes for gene expression analysis in rubber tree leaf

Gene expression profiles are increasingly applied to investigate molecular mechanism for which, normalization with suitable reference genes is critical. Previously we have reported several suitable reference genes for laticifer samples from rubber tree, however, little is known in leaf. The main objective of this current study was to identify some stable expression reference genes at various developmental stages of leaf, as well as during abiotic (high and low temperature extremes) and biotic stresses (pathogen stress). Gene expression profilings identified the ubiquitin-proteasome system as excellent potential as reference genes for rubber tree leaf. Among a total of 30 tested genes investigated, 24 new candidate (including 11 genes involved in the ubiquitin-proteasome system), 4 previously identified and 2 specific genes, were further evaluated using quantitative real-time PCR. Our results indicated that the new candidate genes had better expression stability comparing with others. For instance, an ubiquitin conjugating enzyme (RG0099) and three ubiquitin-protein ligases (RG0928, RG2190 and RG0118) expressed stably in all samples, and were confirmed to be suitable reference genes for rubber tree leaf under four different conditions. Finally, we suggest that using more than one reference gene may be appropriate in gene expression studies when employing different software to normalize gene expression data. Our findings have significant implications for the reliability of data obtained from genomics studies in rubber tree and perhaps in other species.