Environmental sex determination (ESD) occurs in divergent, phylogenetically unrelated taxa, and in some species, co-occurs with genetic sex determination (GSD) mechanisms. Although epigenetic ...regulation in response to environmental effects has long been proposed to be associated with ESD, a systemic analysis on epigenetic regulation of ESD is still lacking. Using half-smooth tongue sole (Cynoglossus semilaevis) as a model-a marine fish that has both ZW chromosomal GSD and temperature-dependent ESD-we investigated the role of DNA methylation in transition from GSD to ESD. Comparative analysis of the gonadal DNA methylomes of pseudomale, female, and normal male fish revealed that genes in the sex determination pathways are the major targets of substantial methylation modification during sexual reversal. Methylation modification in pseudomales is globally inherited in their ZW offspring, which can naturally develop into pseudomales without temperature incubation. Transcriptome analysis revealed that dosage compensation occurs in a restricted, methylated cytosine enriched Z chromosomal region in pseudomale testes, achieving equal expression level in normal male testes. In contrast, female-specific W chromosomal genes are suppressed in pseudomales by methylation regulation. We conclude that epigenetic regulation plays multiple crucial roles in sexual reversal of tongue sole fish. We also offer the first clues on the mechanisms behind gene dosage balancing in an organism that undergoes sexual reversal. Finally, we suggest a causal link between the bias sex chromosome assortment in the offspring of a pseudomale family and the transgenerational epigenetic inheritance of sexual reversal in tongue sole fish.
Chinese tongue sole is a marine fish with ZW sex determination. Genome sequencing suggested that the Z-linked dmrt1 is a putative male determination gene, but direct genetic evidence is still ...lacking. Here we show that TALEN of dmrt1 efficiently induced mutations of this gene. The ZZ dmrt1 mutant fish developed ovary-like testis, and the spermatogenesis was disrupted. The female-related genes foxl2 and cyp19a1a were significantly increased in the gonad of the ZZ dmrt1 mutant. Conversely, the male-related genes Sox9a and Amh were significantly decreased. The dmrt1 deficient ZZ fish grew much faster than ZZ male control. Notably, we obtained an intersex ZW fish with a testis on one side and an ovary on the other side. This fish was chimeric for a dmrt1 mutation in the ovary, and wild-type dmrt1 in the testis. Our data provide the first functional evidence that dmrt1 is a male determining gene in tongue sole.
Blind-side hypermelanosis has emerged as a major concern in commercial rearing environments of the flatfish aquaculture industry. To date, the underlying molecular mechanisms are not well understood. ...To fill this gap, in this study, whole transcriptomic sequencing and analyses were performed using normal skins and hypermelanic skins of the blind side of Chinese tongue sole (
). Differentially expressed long non-coding RNAs (DElncRNAs), miRNAs (DEmiRNAs), and differentially expressed genes as well as their competing endogenous RNA (ceRNA) networks were identified. A total of 34 DElncRNAs, 226 DEmiRNAs, and 610 DEGs were identified. Finally, lncRNA-miRNA-mRNA regulatory networks (involving 29 DElncRNAs, 106 DEmiRNAs, and 162 DEGs) associated with blind-side hypermelanosis were constructed. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses of 162 DEGs in ceRNA networks identified DEGs (e.g., oca2, mc1r, and ihhb) in pigmentation-related biological processes and DEGs (e.g., ca4, glul, and fut9) in nitrogen metabolism, glycosphingolipid biosynthesis, and folate biosynthesis pathways, as well as their corresponding DElncRNAs and DEmiRNAs to potentially play key regulatory roles in blind-side hypermelanosis. In conclusion, this is the first study on the ceRNA regulatory network associated with blind-side hypermelanosis in flatfish. These new findings expand the spectrum of non-coding regulatory mechanisms underpinning blind-side hypermelanosis, which facilitates the further exploration of molecular regulatory mechanisms of malpigmentation in flatfish.
Background: To observe the activation strategies of the ankle muscles using surface electromyography (sEMG) during single-leg standing (SLS) and both-leg standing (BLS) on flat ground (FG), soft mat ...(SM), and BOSU ball (BB) surfaces. Methods: Thirty healthy young adults participated in the study. The muscle activities of the tibialis anterior (TA) and gastrocnemius medial (GM) were measured on the three surfaces during SLS and BLS. Electromyographic evaluations of the TA and GM were recorded during maximum voluntary isometric contractions (MVIC). Muscle activation was evaluated using MVIC%, and muscle co-contraction was evaluated using the co-contraction index (CI). Results: A statistically significant increase was observed in the MVIC% of the TA, GM, and CI on the three surfaces during SLS compared to BLS, except for the comparison of CI on BB between SLS and BLS (t = −1.35, p = 0.19). The MVIC% of the TA and GM during SLS and BLS on BB was significantly increased in comparison with FG and SM. The CI during BLS on BB increased compared to FG (t = 3.19, p < 0.01) and SM (t = 4.64, p < 0.01). The CI during BLS on SM (t = −1.46, p = 0.15) decreased when compared to FG but without statistical significance. Conclusions: SLS and unstable surfaces can induce greater muscle activation, and SLS can have a greater influence on ankle muscles.
In aquaculture, selective breeding for survival till harvest have become an alternative strategy for improving disease resistance and production. However, knowledge of genetic parameters of harvest ...survival, e.g., heritability and genetic correlations between survival and growth rate traits, is still scarce. The aims of this study were to estimate genetic parameters for harvest survival and growth rate traits under commercial farming conditions in olive flounder (Paralichthys olivaceus). Harvest survival was defined as a binary trait; growth traits were measured as average daily gain (ADG), specific growth rate (SGR), daily growth coefficient (DGC) and body weight (BW). Data from a population of 241 full-sib families (involving 39,904 individuals, four generations) were used. Heritabilities of survival were low but significant, which were 0.15 ± 0.04 and 0.22 ± 0.01 on observed and underlying scale, respectively. Heritability estimates for ADG, SGR and DGC were medium to high, which were 0.33 ± 0.06, 0.83 ± 0.07, 0.58 ± 0.07, respectively. While the heritability of BW is of low magnitude (0.17 ± 0.08). The genetic correlations between harvest survival and three growth rate traits (i.e., ADG, SGR and DGC) were very strong (ranging from 0.66 to 0.79), which is an exciting result. However, the genetic correlation between harvest survival and BW was much lower (0.17 ± 0.08). These results suggest that selection for harvest survival would consequentially result in concomitant increase of growth rate, and vice versa. Our findings revealed novel insights into the genetic improvement of growth rate and harvest survival through genetic selection in olive flounder.
Abstract
Background
Infectious diseases have caused huge economic loss and food security issues in fish aquaculture. Current management and breeding strategies heavily rely on the knowledge of ...regulative mechanisms underlying disease resistance. Though the intestinal microbial community was linked with disease infection, there is little knowledge about the roles of intestinal microbes in fish disease resistance.
Cynoglossus semilaevis
is an economically important and widely cultivated flatfish species in China. However, it suffers from outbreaks of vibriosis, which results in huge mortalities and economic loss.
Results
Here, we used
C. semilaevis
as a research model to investigate the host-microbiome interactions in regulating vibriosis resistance. The resistance to vibriosis was reflected in intestinal microbiome on both taxonomic and functional levels. Such differences also influenced the host gene expressions in the resistant family. Moreover, the intestinal microbiome might control the host immunological homeostasis and inflammation to enhance vibriosis resistance through the microbe-intestine-immunity axis. For example,
Phaeobacter
regulated its
hdhA
gene and host
cyp27a1
gene up-expressed in bile acid biosynthesis pathways, but regulated its
trxA
gene and host
akt
gene down-expressed in proinflammatory cytokines biosynthesis pathways, to reduce inflammation and resist disease infection in the resistant family. Furthermore, the combination of intestinal microbes and host genes as biomarkers could accurately differentiate resistant family from susceptible family.
Conclusion
Our study uncovered the regulatory patterns of the microbe-intestine-immunity axis that may contribute to vibriosis resistance in
C. semilaevis
. These findings could facilitate the disease control and selective breeding of superior germplasm with high disease resistance in fish aquaculture.
Blind-side hypermelanosis is an emerging concern across the flatfish aquaculture industry including Chinese tongue sole (
Cynoglossus semilaevis
). Circular RNAs (circRNAs) as endogenous non-coding ...RNAs have been acknowledged to play important roles in various biological processes. However, the underlying regulatory mechanisms of circRNAs involved in flatfish blind-side hypermelanosis remain unclear. In this study, to profile the circRNA expression pattern and circRNA-microRNA-messenger RNA (mRNA) network, high-throughput sequencing was performed by using blind-side normal and hypermelanotic skins of tongue sole. A total of 73 differentially expressed circRNAs were identified, and the competing endogenous RNA (ceRNA) network was constructed. Furthermore, circRNA host genes and mRNAs involved in ceRNA network were subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes enrichment analyses. Several GO terms and pathways of biological significance were identified and well addressed the generation of blind-side hypermelanosis to some extent. These interesting results extend the understanding of the functional profile of circRNAs and yield valuable insights into the molecular regulatory mechanisms of hypermelanosis in flatfish.
Edwardsiella tarda causes acute symptoms with ascites in Japanese flounder (Paralichthys olivaceus) and is a major problem for China's aquaculture sector. Genomic selection (GS) has been widely ...adopted in breeding industries because it shortens generation intervals and results in the selection of individuals that have great breeding potential with high accuracy. Based on an artificial challenge test and re-sequenced data of 1099 flounders, the aims of this study were to estimate the genetic parameters of resistance to E. tarda in Japanese flounder and to evaluate the accuracy of single-step GBLUP (ssGBLUP), weighted ssGBLUP (WssGBLUP), and BayesB for improving resistance to E. tarda by using three subsets of pre-selected single nucleotide polymorphisms (SNPs). In addition, SNPs that are associated with this trait were identified using a single-SNP genome-wide association study (GWAS) and WssGBLUP.
We estimated a heritability of 0.13 ± 0.02 for resistance to E. tarda in Japanese flounder. One million SNPs at fixed intervals were selected from 4,978,724 SNPs that passed quality controls. GWAS identified significant SNPs on chromosomes 14 and 24. WssGBLUP revealed that the putative quantitative trait loci on chromosomes 1 and 14 contained SNPs that explained more than 1% of the genetic variance. Three 50 k-SNP subsets were pre-selected based on different criteria. Compared with pedigree-based prediction (ABLUP), the three genomic methods evaluated resulted in at least 7.7% greater accuracy of predictions. The accuracy of these genomic prediction methods was almost unchanged when pre-selected trait-related SNPs were used for prediction.
Resistance to E. tarda in Japanese flounder has a low heritability. GWAS and WssGBLUP revealed that the genetic architecture of this trait is polygenic. Genomic prediction of breeding values performed better than ABLUP. It is feasible to implement genomic selection to increase resistance to E. tarda in Japanese flounder with 50 k SNPs. Based on the criteria used here, pre-selection of SNPs was not beneficial and other criteria for pre-selection should be considered.
Chinese tongue sole (
Cynoglossus semilaevis
) is an economically important marine fish in China. Generally, the eyeless side of the Chinese tongue sole is white and the side with eyes is brown after ...metamorphosis, hypermelanosis may still occur in the eyeless side in certain individuals after metamorphosis, which greatly decreases consumer acceptance and market price. In order to study the possibility of genetic improvement, we determined genomic markers in Chinese tongue sole using the genotyping-by-sequencing method and analyzed their association with hypermelanosis area. Genetic analysis showed that hypermelanosis was a complicated quantitative trait, and the estimated heritability for hypermelanosis incidence and area ratio were 0.16 and 0.21, respectively. Genomic selection analysis showed that selection based on hypermelanosis incidence and area ratio had similar reliabilities and prediction accuracies, indicating the feasibility of genetic improvement. Nine loci were significantly associated with hypermelanosis, few of which included genes or flanked genes potentially associated with skin disease, indicating the potential complicated genetic mechanisms underlying hypermelanosis in the Chinese tongue sole.
The tiger pufferfish (Takifugu rubripes), also known as fugu, has recently suffered from severe C. irritans infections under aquaculture environment, yet the underlying immune mechanisms against the ...parasite remain poorly understood. In this study, we conducted a comprehensive transcriptome analysis of the gill tissue from infected and uninfected fish using PacBio long-read (one pooled sample each for seriously infected and healthy individuals, respectively) and Illumina short-read (three pools for mildly infected, seriously infected, and healthy individuals, respectively) RNA sequencing technologies. After aligning sequence data to fugu’s reference genome, 47,307 and 34,413 known full-length transcripts were identified and profiled in healthy and infected fish, respectively. Similarly, we identified and profiled 1126 and 803 novel genes that were obtained from healthy and infected fish, respectively. Interestingly, we found a decrease in the number of alternative splicing (AS) events and long non-coding RNAs (lncRNAs) after infection with C. irritans, suggesting that they may be involved in the regulation of the immune response in fugu. There were 687 and 1535 differentially expressed genes (DEGs) in moderately and heavily infected fish, respectively, compared to uninfected fish. Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that immune-related DEGs in the two comparison groups were mainly enriched in cytokine-cytokine receptor interactions, ECM-receptor interactions, T-cell receptor signaling pathways, Th1 and Th2 cell differentiation, and Th17 cell differentiation pathways. Further analysis revealed that a large number of immune-related genes were downregulated in infected fish relative to uninfected ones, such as CCR7, IL7R, TNFRSF21, CD4, COL2A1, FOXP3B, and ITGA8. Our study suggests that C. irritans is potentially a highly efficient parasite that may disrupt the defense mechanisms of fugu against it. In addition, in combination of short-read RNA sequencing and previous genome-wide association analyses, we identified five key genes (NDUFB6, PRELID1, SMOX, SLC25A4, and DENND1B) that might be closely associated with C. irritans resistance. This study not only provides valuable resources of novel genic transcripts for further research, but also provides new insights into the immune mechanisms underlying C. irritans infection response in farmed fugu.