Identification of suitable reference genes for quantitative real-time PCR normalization in blotched snakehead Channa maculata

Histological observations and transcriptome analyses reveal the dynamic changes in the gonads of the blotched snakehead (Channa maculata) during sex differentiation and gametogenesis

BACKGROUND

Blotched snakehead (Channa maculata) displays significant sexual dimorphism, with males exhibiting faster growth rates and larger body sizes compared to females. The cultivation of the all-male population of snakeheads holds substantial economic and ecological value. Nonetheless, the intricate processes governing the development of bipotential gonads into either testis or ovary in C. maculata remain inadequately elucidated. Therefore, it is necessary to determine the critical time window of sex differentiation in C. maculata, providing a theoretical basis for sex control in production practices.

METHODS

The body length and weight of male and female C. maculata were measured at different developmental stages to reveal when sexual dimorphism in growth initially appears. Histological observations and spatiotemporal comparative transcriptome analyses were performed on ovaries and testes across various developmental stages to determine the crucial time windows for sex differentiation in each sex and the sex-related genes. Additionally, qPCR and MG2C were utilized to validate and locate sex-related genes, and levels of E2 and T were quantified to understand sex steroid synthesis.

RESULTS

Sexual dimorphism in growth became evident starting from 90 dpf. Histological observations revealed that morphological sex differentiation in females and males occurred between 20 and 25 dpf or earlier and 30-35 dpf or earlier, respectively, corresponding to the appearance of the ovarian cavity or efferent duct anlage. Transcriptome analyses revealed divergent gene expression patterns in testes and ovaries after 30 dpf. The periods of 40-60 dpf and 60-90 dpf marked the initiation of molecular sex differentiation in females and males, respectively. Male-biased genes (Sox11a, Dmrt1, Amh, Amhr2, Gsdf, Ar, Cyp17a2) likely play crucial roles in male sex differentiation and spermatogenesis, while female-biased genes (Foxl2, Cyp19a1a, Bmp15, Figla, Er) could be pivotal in ovarian differentiation and development. Numerous biological pathways linked to sex differentiation and gametogenesis were also identified. Additionally, E2 and T exhibited sexual dimorphism during sex differentiation and gonadal development. Based on these results, it is hypothesized that in C. maculata, the potential male sex differentiation pathway, Sox11a-Dmrt1-Sox9b, activates downstream sex-related genes (Amh, Amhr2, Gsdf, Ar, Cyp17a2) for testicular development, while the antagonistic pathway, Foxl2/Cyp19a1a, activates downstream sex-related genes (Bmp15, Figla, Er) for ovarian development.

CONCLUSIONS

This study provides a comprehensive overview of gonadal dynamic changes during sex differentiation and gametogenesis in C. maculata, establishing a scientific foundation for sex control in this species.

Evaluation of housekeeping genes as references for quantitative real-time PCR analysis of European eel, Anguilla anguilla

Eels are important aquaculture species for which an increasing number of reference genes are being identified and applied. In this study, five housekeeping genes [RPL7 (ribosomal protein L7), 18 S (18 S ribosomal RNA), EF1A (elongation factor 1α), ACTB (β-actin) and GAPDH (glyceraldehyde-3-phosphate dehydrogenase)] were chosen to evaluate their reliability as reference genes for quantitative real-time PCR (qPCR) for the study of Anguilla anguilla. The expression of the selected genes in different eel tissues was determined using qPCR at different growth stages or upon challenge by Anguillid herpesvirus (AngHV), and the expression levels of these genes were then compared and evaluated using the geNorm and NormFinder algorithms. Then, RefFinder was used to comprehensively rank the examined housekeeping genes. Interestingly, the expression of the evaluated housekeeping genes exhibited tissue-dependent and treatment-dependent variations. In different growth periods A. anguilla tissues, the most stable genes were the following: ACTB in mucus; 18 S in skin and kidney; RPL7 in muscle, gill, intestine and brain; EF1A in heart and liver; and GAPDH in spleen. In contrast, in AngHV-challenged A. anguilla tissues, the most stable genes were the following: 18 S in mucus; RPL7 in skin, gill, heart, spleen, kidney and intestine; EF1A in muscle and liver; and ACTB in brain. Further comparison analysis indicated that the expression of RPL7 and EF1A was stable in multiple A. anguilla tissues in different growth periods and in eels challenged by AngHV. Nonetheless, the expression level of GAPDH in eel tissues was lower, and it was unstable in several tissues. These results indicated that the selection of reference genes for qPCR analysis in A. anguilla should be made in accordance with experimental parameters, and both RPL7 and EF1A could be used as reference genes for qPCR study of A. anguilla at different growth stages or upon challenge by AngHV. The reference genes identified in this study could improve the accuracy of qPCR data and facilitate further studies aimed at understanding the biology of eels.

Characterization, expression and CpG methylation analysis of Dmrt1 and its response to steroid hormone in blotched snakehead (Channa maculata)

Dmrt1 is an important transcriptional regulator that plays critical role in male gonadogenesis, testicular differentiation and development. In this study, Dmrt1 was cloned from blotched snakehead (Channa maculata), which is designated as CmDmrt1. CmDmrt1 encoded a putative protein with 293 amino acids and presented an extremely conserved DM domain. It was nearly expressed in the gonads, and the expression was more than 15 times higher in the testis than in the ovary. 1851 bp promoter sequence of CmDmrt1 was characterized and the methylation levels of the CpG sites were analyzed to detect sex-related differences. A significant negative correlation between CmDmrt1 expression and CpG methylation level of its promoter was found in the testis and ovary. During gonadal development, CmDmrt1 transcription displayed strong male-biased expression patterns, increased with the maturation of testis and reached the peak at 195 days after hatching (dah), which indicates a significant role of Dmrt1 in spermatogenesis. Steroid treatment could influence CmDmrt1 expression, and long-term 17β-estradiol (E2) treatment could induce the male-to-female secondary sex reversal (SSR), which resulted in the differentiated testis transformed to ovary or ovotestis. Meanwhile, CmDmrt1 expression was down-regulated to fairly low level in the ovary of the SSR XY fish, which was similar to that in normal XX females ovary. Our research illustrates that Dmrt1 is linked to testis differentiation and spermatogenesis in blotched snakehead, providing information for functional studies on sex differentiation and gonadal development of C. maculata, and scientific basis for the production practice of all-male snakehead breeding.

Comparative transcriptome analysis on four types of gonadal tissues of blotched snakehead (Channa maculata)

Blotched snakehead (Channa maculata) is an economically important freshwater fish in China, of which males grow much faster than females. To illuminate the molecular mechanism of sex differentiation and gonad development, RNA-Sequencing was performed to identify sex-related genes and pathway in gonads of 6-month-old normal XX females (XX-F), normal XY males (XY-M), XY sex reversal females (XY-F) and YY super-males (YY-M). The analysis showed that many differentially expressed genes (DEGs) had similar expression patterns in XY-F and XX-F, which were different from XY-M and YY-M. qRT-PCR indicated that Amh, Dmrt1, and Sox9 had relatively high expression in testes of XY-M and YY-M. Taking Amh as an example, there was a relative fold change of 1.0 in XX-F, 2.1 fold change in XY-F, 36.1 fold change in XY-M, and 26.0 fold change in YY-M. Cyp19a1a, Figla, and Foxl2 were highly expressive in ovaries of XX-F and XY-F. Taking Figla as an example, there was a relative fold change of 557 in XX-F, 304.5 fold change in XY-F, 5.6 fold change in XY-M, and 4.4 fold change in YY-M. KEGG analysis revealed many DEGs distributed in pathways related to sex differentiation, steroid hormone synthesis and growth, etc. Significant variation and trends in relative expression levels tested by qRT-PCR were consistent with those recorded by RNA-Sequencing. This is the first time that transcriptome of snakehead has been investigated systematically and in an integrated way. Large quantities of candidate genes involved in sex differentiation, gonad development and growth dimorphism were identified. The study provides useful resources for understanding sex differentiation and growth dimorphism, potentially assisting mono-sex production of snakehead in aquaculture.

Selection of reference genes for quantitative real-time PCR normalisation in largemouth bass Micropterus salmoides fed on alternative diets

The partial cDNA sequences of eight reference genes (actb, tuba1, gapdh58, gapdh59, eef1a1, RNA 18 s, pabpc1, ube2I) were cloned from largemouth bass Micropterus salmoides. The expression levels of these eight genes were compared in the various tissues (eye, spleen, kidney, gill, muscle, brain, liver, heart, gut and gonad) of M. salmoides fed on forage fish. The results showed that the candidate genes exhibited tissue-specific expression to various degrees and the stability ranking order was eef1a1 > tuba1 > RNA 18 s > pabpc1 > ube2I > actb > gapdh58 > gapdh59 among tissue types. Four candidate genes eef1a1, tuba1, RNA 18 s and actb were used to analyse the stability in liver tissues of largemouth bass between the forage-fish group and the formulated-feed group. The candidate genes also showed some changes in expression levels in the livers, while eef1a1 and tuba1 had the most stable expression in livers of fish fed on alternative diets within 10 candidates. So eef1a1 and tuba1 were recommended as optimal reference gene in quantitative real-time PCR analysis to normalise the expression levels of target genes in tissues and lives of the M. salmoides fed on alternative diets. In livers, the expression levels of gck normalised by eef1a1 and tuba1 showed the significant up-regulation in formulated feed group (P < 0.05) than those in forage-fish group. While sex difference has no significant effects on the expression levels of gck in both groups.

The DNA methylation level is associated with the superior growth of the hybrid fry in snakehead fish (Channa argus × Channa maculata)

Hybrid vigour, or heterosis, refers to the increased productivity and growth rate of hybrid offsprings relative to the parents. Various heterosis have been well exploited in fish for fisheries. However, the molecular mechanisms underlying heterosis are largely unknown in fish. In this study, two inbred and hybrid lines between the northern snakehead (NS, Channa argus) and blotched snakehead (BS, Channa maculata) were generated. The analysis on various growth traits, including body length, head length, and body height, showed that hybrid fry obviously exhibited a spontaneous growth heterosis over the inbred. Moreover, the methylation-sensitive amplification polymorphism (MSAP) analysis revealed that the DNA methylation levels were negatively related to the body growth in all fry. Especially, the DNA methylation levels in the hybrid fry were significantly lower than those in the inbred. Additionally, qRT-PCR showed that the snakehead fish Dnmt3a mRNA was initially detectable in embryos at 12 hpf and gradually increased as developing. Intriguingly, the level of Dnmt3a mRNA expression was found to be closely correlated to the DNA methylation level in embryos/fry. The results of this study firstly demonstrated the correlations between growth heterosis, DNA methylation level and Dnmt3a mRNA expression in fish fry. The findings of this study implied that the hybrids' heterosis formation is probably accompanied by DNA methylation alterations and modulated by Dnmt3a gene in fish. This study would provide new clues for further investigations on mechanisms behind heterosis formation in fish hybrid.