Polymorphism in exons CpG rich regions of the cyp17-II gene affecting its mRNA expression and reproductive endocrine levels in female Japanese flounder (Paralichthys olivaceus)

Identification and Functional Analysis of Ras-Related Associated with Diabetes Gene (rrad) in Edwardsiella piscicida-Resistant Individuals of Japanese Flounder (Paralichthys olivaceus)

Ras-related associated with diabetes (RRAD) is a member of the Ras GTPase superfamily that plays a role in several cellular functions, such as cell proliferation and differentiation. In particular, the superfamily acts as an NF-κB signaling pathway inhibitor and calcium regulator to participate in the immune response pathway. A recent transcriptome study revealed that rrad was expressed in the spleen of disease-resistant Japanese flounder (Paralichthys olivaceus) individuals compared with disease-susceptible individuals, and the results were also verified by qPCR. Thus, the present study aimed to explore how rrad regulates antimicrobial immunity via the NF-κB pathway. First, the coding sequence of P. olivaceus rrad was identified. The sequence was 1092 bp in length, encoding 364 amino acids. Based on phylogenetic and structural relationship analyses, P. olivaceus rrad appeared to be more closely related to teleosts. Next, rrad expression differences between disease-resistant and disease-susceptible individuals in immune-related tissues were evaluated, and the results revealed that rrad was expressed preferentially in the spleen of disease-resistant individuals. In response to Edwardsiella piscicida infection, rrad expression in the spleen changed. In vitro, co-culture was carried out to assess the hypo-methylated levels of the rrad promoter in the disease-resistant spleen, which was consistent with the high mRNA expression. The siRNA-mediated knockdown of rrad performed with the gill cell line of P. olivaceus affected many rrad-network-related genes, i.e., dcp1b, amagt, rus1, rapgef1, ralbp1, plce1, rasal1, nckipsd, prkab2, cytbc-1, sh3, and others, as well as some inflammation-related genes, such as bal2 and Il-1β. In addition, flow cytometry analysis showed that rrad overexpression was more likely to induce cell apoptosis, with establishing a link between rrad's function and its potential roles in regulating the NF-κB pathway. Thus,. the current study provided some clarity in terms of understanding the immune response about rrad gene differences between disease-resistant and disease-susceptible P. olivaceus individuals. This study provides a molecular basis for fish rrad gene functional analysis and may serve as a reference for in-depth of bacterial disease resistance of teleost.

Gonadal Development and Differentiation of Hybrid F1 Line of Ctenopharyngodon idella (♀) × Squaliobarbus curriculus (♂)

The hybrid F1 offspring of Ctenopharyngodon idella (♂) and Squaliobarbus curriculus (♀) exhibit heterosis in disease resistance and also show abnormal sex differentiation. To understand the mechanism behind gonadal differentiation in the hybrid F1, we analyzed the transcriptomes of C. idella, S. curriculus, and the hybrid F1; screened for genes related to gonad development in these samples; and measured their expression levels. Our results revealed that compared to either C. idella or S. curriculus, the gene expressions in most sub-pathways of the SNARE interactions in the vesicular transport pathway in the hypothalamus, pituitary, and gonadal tissues of their hybrid F1 offspring were significantly up-regulated. Furthermore, insufficient transcription of genes involved in oocyte meiosis may be the main reason for the insufficient reproductive ability of the hybrid F1 offspring. Through transcriptome screening, we identified key molecules involved in gonad development, including HSD3B7, HSD17B1, HSD17B3, HSD20B2, CYP17A2, CYP1B1, CYP2AA12, UGT2A1, UGT1A1, and FSHR, which showed significant differences in expression levels in the hypothalamus, pituitary, and gonads of these fish. Notably, the expression levels of UGT1A1 in the gonads of the hybrid F1 were significantly higher than those in C. idella and S. curriculus. These results provide a scientific basis for further research on the gonadal differentiation mechanism of hybrid F1 offspring.

Exploring the influence of DNA methylation and single nucleotide polymorphisms of the Myostatin gene on growth traits in the hybrid grouper (Epinephelus fuscoguttatus (female) × Epinephelus polyphekadion (male))

Investigations into the correlation between growth characteristics and DNA methylation levels, along with genetic variations, can provide fundamental insights to enhance growth performance in groupers. The Myostatin (mstn) gene plays a vital role in regulating skeletal muscle development and growth. This study scrutinized the DNA methylation levels of the mstn gene across hybrid groupers (E. fuscoguttatus (♀) × E. polyphekadion (♂)) and their parental species, to evaluate its impact on growth attributes in grouper fish. The nucleotide sequence of the mstn gene was directly sequenced in the hybrid grouper, exhibiting different growth performance to identify the single nucleotide polymorphisms (SNPs) of the mstn gene and explore their correlation with growth characteristics. The findings revealed no significant differences in global DNA methylation levels within muscle tissue among the hybrid grouper and parents. However, significant differences in DNA methylation sites were discovered between the hybrid grouper and E. polyphekadion at sites 824 and 1521 (located at exon 2 and intron 2, respectively), and between E. fuscoguttatus and E. polyphekadion at site 1521. These variations could potentially influence the mRNA expression of the mstn gene. The study also identified that SNP g.1003 T > C in exon 2 of the mstn gene was significantly associated with various growth traits including body weight, total length, body length, head length, caudal peduncle height, and body height (p < 0.01). Specimens with the TT genotype at site 1003 demonstrated superior growth performance compared to those with the TC genotype. Furthermore, microstructural analyses of muscle tissue showed that the average area and diameter of muscle fibers in TT genotype individuals were significantly greater than those in TC genotype individuals. Therefore, this research provides robust evidence linking the DNA methylation level and polymorphisms of the mstn gene with growth traits, which could be beneficial for grouper breeding programs.

Broodstock nutritional programming differentially affects the hepatic transcriptome and genome-wide DNA methylome of farmed gilthead sea bream (Sparus aurata) depending on genetic background

BACKGROUND

Broodstock nutritional programming improves the offspring utilization of plant-based diets in gilthead sea bream through changes in hepatic metabolism. Attention was initially focused on fatty acid desaturases, but it can involve a wide range of processes that remain largely unexplored. How all this can be driven by a different genetic background is hardly underlined, and the present study aimed to assess how broodstock nutrition affects differentially the transcriptome and genome-wide DNA methylome of reference and genetically selected fish within the PROGENSA® selection program.

RESULTS

After the stimulus phase with a low fish oil diet, two offspring subsets of each genetic background received a control or a FUTURE-based diet. This highlighted a different hepatic transcriptome (RNA-seq) and genome-wide DNA methylation (MBD-seq) pattern depending on the genetic background. The number of differentially expressed transcripts following the challenge phase varied from 323 in reference fish to 2,009 in genetically selected fish. The number of discriminant transcripts, and associated enriched functions, were also markedly higher in selected fish. Moreover, correlation analysis depicted a hyper-methylated and down-regulated gene expression state in selected fish with the FUTURE diet, whereas the opposite pattern appeared in reference fish. After filtering for highly represented functions in selected fish, 115 epigenetic markers were retrieved in this group. Among them, lipid metabolism genes (23) were the most reactive following ordering by fold-change in expression, rendering a final list of 10 top markers with a key role on hepatic lipogenesis and fatty acid metabolism (cd36, pitpna, cidea, fasn, g6pd, lipt1, scd1a, acsbg2, acsl14, acsbg2).

CONCLUSIONS

Gene expression profiles and methylation signatures were dependent on genetic background in our experimental model. Such assumption affected the magnitude, but also the type and direction of change. Thus, the resulting epigenetic clock of reference fish might depict an older phenotype with a lower methylation for the epigenetically responsive genes with a negative methylation-expression pattern. Therefore, epigenetic markers will be specific of each genetic lineage, serving the broodstock programming in our selected fish to prevent and mitigate later in life the risk of hepatic steatosis through changes in hepatic lipogenesis and fatty acid metabolism.

Dynamic epimarks in sex-related genes predict gonad phenotype in the European sea bass, a fish with mixed genetic and environmental sex determination

The integration of genomic and environmental influences into methylation patterns to bring about a phenotype is of central interest in developmental epigenetics, but many details are still unclear. The sex ratios of the species used here, the European sea bass, are determined by genetic and temperature influences. We created four families from parents known to produce offspring with different sex ratios, exposed larvae to masculinizing temperatures and examined, in juvenile gonads, the DNA methylation of seven genes related to sexual development by a targeted sequencing approach. The genes most affected by both genetics and environment were cyp19a1a and dmrt1, with contrasting sex-specific methylation and temperature responses. The relationship between cyp19a1a methylation and expression is relevant to the epigenetic regulation of vertebrate sex, and we report the evidence of such relationship only below a methylation threshold, ~ 80%, and that it was sex-specific: negatively correlated in females but positively correlated in males. From parents to offspring, the methylation in gonads was midway between oocytes and sperm, with bias towards oocytes for amh-r2, er-β2, fsh-r and cyp19a1a. In contrast, dmrt1 levels resembled those of sperm. The methylation of individual CpGs from foxl2, er-β2 and nr3c1 were conserved from parents to offspring, whereas those of cyp19a1a, dmrt1 and amh-r2 were affected by temperature. Utilizing a machine-learning procedure based on the methylation levels of a selected set of CpGs, we present the first, to our knowledge, system based on epigenetic marks capable of predicting sex in an animal with ~ 90% accuracy and discuss possible applications.

Hypothalamus-pituitary-gonad axis of rainbow trout (Oncorhynchus mykiss) during early ovarian development and under dense rearing condition

The objective of this study was to determine the hypothalamus-pituitary-gonad (HPG) axis of female rainbow trout (Oncorhynchus mykiss) during early ovarian development and under high rearing density. Trouts were sampled from 240 (ovarian stage II) to 540 (ovarian stage IV) days following hatching (DFH) as control group (Ctrl, 4.6-31.1kg/m(3)) to determine HPG axis during early ovarian development. Trouts from the same batch of fertilized eggs were reared in two higher densities during 240-540 DFH as stocking density 1 and 2 (SD1, 6.6-40.6kg/m(3); SD2, 8.6-49.3kg/m(3)) to elucidate effects of high density on reproductive parameters. Dopamine, E2 (estradiol), 17α,20β-P (17α,20β-dihydroxy4-pregnen-3-one) and P4 (progesterone) were evaluated by radioimmunoassay or ELISA. mRNA expression of hypothalamic gnrh-1, -2 (gonadotropin-releasing hormone-1, -2), pituitary gonadotropins (fsh/lh, follicle-stimulating hormone/luteinizing hormone) and their cognate receptors (fshr/lhr) in ovaries were examined by qRT-PCR. Our findings demonstrated mRNA expression of hypothalamic sgnrh-1, pituitary fsh and ovarian fshr increased in early ovarian development (360 DFH). Serum 17α,20β-P and pituitary lh mRNA expression first increased when trouts were in ovarian stage III (420 DFH). Ovaries were at different stages when reared in different densities. Long-term high density treatment (over 31.7kg/m(3)) resulted in decreased hypothalamic sgnrh-1, pituitary fsh, ovarian fshr, serum E2, and increased hypothalamus gnrh-2 and serum dopamine during vitellogenin synthesis, suggesting HPG of rainbow trout might be retarded under dense rearing condition.

Genetic polymorphisms and DNA methylation in exon 1 CpG-rich regions of PACAP gene and its effect on mRNA expression and growth traits in half smooth tongue sole (Cynoglossus semilaevis)

The pituitary adenylate cyclase activating polypeptide (PACAP) is a new type of hypophysiotropic hormone and plays an important role in regulating the synthesis and secretion of growth hormone and gonadotropin. The research on the relationship between PACAP and different growth traits would contribute to explain its function during the process of growth. Moreover, epigenetic modifications, especially DNA methylation at the CpG sites of the SNPs, play important roles in regulating gene expression. The results suggest that a SNP mutation (c.C151G) in the PACAP gene of male half smooth tongue sole (Cynoglossus semilaevis) is significantly associated with growth traits and serum physiological and biochemical parameters such as inorganic phosphorus (P < 0.05). The SNP is located in a CpG-rich region of exon 1. Intriguingly, the transition (C→G) added a new methylation site of PACAP gene. This SNP was also significantly related to the expression and methylation level of PACAP (P < 0.05). Individuals with GG genotype had faster growth rates than those of CG and CC genotypes. Moreover, GG genotype had significantly higher PACAP expression level and lower methylation level than CG and CC genotypes. In the serum indexes, only inorganic phosphorus content within GG genotypes was significantly higher than CC genotypes. This implied that the mutation and methylation status of PACAP gene could influence growth traits and this locus could be considered as a candidate genetic or epigenetic marker for Cynoglossus semilaevis molecular breeding.

Gonadal Transcriptome Analysis in Sterile Double Haploid Japanese Flounder

Sterility is a serious problem that can affect all bionts. In teleosts, double haploids (DHs) induced by mitogynogenesis are often sterile. This sterility severely restricts the further application of DHs for production of clones, genetic analysis, and breeding. However, sterile DH individuals are good source materials for investigation of the molecular mechanisms of gonad development, especially for studies into the role of genes that are indispensable for fish reproduction. Here, we used the Illumina sequencing platform to analyze the transcriptome of sterile female DH Japanese flounder in order to identify major genes that cause sterility and to provide a molecular basis for an intensive study of gonadal development in teleosts. Through sequencing, assembly, and annotation, we obtained 52,474 contigs and found that 60.7% of these shared homologies with existing sequences. A total of 1225 differentially expressed unigenes were found, including 492 upregulated and 733 downregulated genes. Gene Ontology and KEGG analyses showed that genes showing significant upregulation, such as CYP11A1, CYP11B2, CYP17, CYP21, HSD3β, bcl2l1, and PRLR, principally correlated with sterol metabolic process, steroid biosynthetic process, and the Jak-stat signaling pathway. The significantly downregulated genes were primarily associated with immune response, antigen processing and presentation, cytokine–cytokine receptor interaction, and protein digestion and absorption. Using a co-expression network analysis, we conducted a comprehensive comparison of gene expression in the gonads of fertile and sterile female DH Japanese flounder. Identification of genes showing significantly different expression will provide further insights into DH reproductive dysfunction and oocyte maturation processes in teleosts.

Steroids in teleost fishes: A functional point of view

Steroid hormones are involved in the regulation of a variety of processes like embryonic development, sex differentiation, metabolism, immune responses, circadian rhythms, stress response, and reproduction in vertebrates. Teleost fishes and humans show a remarkable conservation in many developmental and physiological aspects, including the endocrine system in general and the steroid hormone related processes in particular. This review provides an overview of the current knowledge about steroid hormone biosynthesis and the steroid hormone receptors in teleost fishes and compares the findings to the human system. The impact of the duplicated genome in teleost fishes on steroid hormone biosynthesis and perception is addressed. Additionally, important processes in fish physiology regulated by steroid hormones, which are most dissimilar to humans, are described. We also give a short overview on the influence of anthropogenic endocrine disrupting compounds on steroid hormone signaling and the resulting adverse physiological effects for teleost fishes. By this approach, we show that the steroidogenesis, hormone receptors, and function of the steroid hormones are reasonably well understood when summarizing the available data of all teleost species analyzed to date. However, on the level of a single species or a certain fish-specific aspect of physiology, further research is needed.

Polymorphisms and DNA methylation level in the CpG site of the GHR1 gene associated with mRNA expression, growth traits and hormone level of half-smooth tongue sole (Cynoglossus semilaevis)

The objectives of the present study were to estimate the GHR1 gene mutations and methylation status of CpGs, and whether those mutations and methylation were involved in the regulation of GHR1 gene expression, hormone level and growth traits in half-smooth tongue sole (Cynoglossus semilaevis). Identification of single-nucleotide polymorphisms was performed on 43 male fish. Through polymerase chain reaction-single-strand conformation polymorphism and sequencing, two SNPs were found. SNP1 [c.G1357A (p.Val376Ile)] creating one CpG site located in exon 8 was named L1 locus, and SNP2 (c.G1479A) located in exon 9 was named L2 locus. Individuals were divided into three genotypes, AA, AG and GG according to L1 locus (GG genotype had one more CpG site because of the mutation), and into two genotypes, AA- and GG-based on L2 locus. The results showed that only L1 locus was significantly associated with body weight (P < 0.01), gonad weight (P ≤ 0.05), triiodothyronine (T3) level (P ≤ 0.05) and mRNA expression (P < 0.01). At L1 locus, newly created CpG site in GG genotype was highly methylated (93.3 %), while there was no difference of methylation level in the other two CpG sites among three genotypes. AA genotype and AG genotype having higher T3 level were significantly different (P ≤ 0.05) from GG genotype. There were significant differences among body weights of AA, AG and GG genotypes (P < 0.01). Gonad weights of AA genotype and AG genotype were significantly lower than GG genotype. The GHR1 mRNA expression of GG genotype was significantly lower than AA and AG genotypes (P < 0.01). These implied that mutations and methylation status of GHR1 gene might influence the hormone level, growth traits and gene expression in male half-smooth tongue sole and the L1 locus could be regarded as a potential candidate genetic and epigenetic marker in half-smooth tongue sole selection.

CpG methylation of dmrt1 and cyp19a promoters in relation to their sexual dimorphic expression in the Japanese flounder Paralichthys olivaceus

To better understand the effects of DNA methylation on the expression patterns of dmrt1 (Doublesex and Mab-3-related transcription factor 1) and cyp19a (Cytochrome P450 19a) in the Japanese flounder Paralichthys olivaceus, quantitative expressions, cellular distributions and cytosine-p-guanine (CpG) methylation patterns of these two genes in the gonads were analysed. The results showed that P. olivaceus dmrt1 expression was 70 times higher in the testis than in the ovary (P < 0·05). Its mRNA was detected clearly in spermatocytes and Sertoli cells of the testis, but weakly in the ovary. Paralichthys olivaceus cyp19a expression was 40 times higher in the ovary than in the testis (P < 0·01). Its mRNA was detected clearly in follicular cells of the ovary, but weakly in spermatocytes of the testis. The dmrt1 promoter CpGs were not methylated in the testis, whereas 57·69% were methylated in the ovary. For the cyp19a promoter CpGs, 97·5% were methylated in the testis and 73·33% were methylated in the ovary. These findings demonstrate that P. olivaceus dmrt1 and cyp19a are sex-related genes with sexual dimorphic expression, CpG methylation levels of the two genes are consistent with their expression quantities, and this epigenetic modification can influence the differential expression of genes in the gonads of P. olivaceus.