Genome-wide identification, evolution of DNA methyltransferases and their expression under salinity stress in Larimichthys crocea.

DNA methyltransferases (Dnmts) are responsible for DNA methylation which influences patterns of gene expression and plays a crucial role in response to environmental changes. In this study, 7 LcDnmt genes were identified in the genome of large yellow croaker (Larimichthys crocea). The comprehensive analysis was conducted on gene structure, protein and location site of LcDnmts. LcDnmt proteins belonged to three groups (Dnmt1, Dnmt2, and Dnmt3) according to their conserved domains and phylogenetic analysis. Although Dnmt3 can be further divided into three sub groups (Dnmt3a, Dnmt3b, and Dnmt3l), there is no Dnmnt3l member in the large yellow croaker. Phylogenetic analysis revealed that the Dnmt family was highly conserved in teleosts. Expression patterns derived from the RNA-seq, qRT-PCR and Western blot analysis revealed that 2 LcDnmt genes (LcDnmt1 and LcDnmt3a2) significantly regulated under salinity stress in the liver, which was found to be dominantly expressed in the intestine and brain, respectively. These two genes may play an important role in the salinity stress of large yellow croaker and represent candidates for future functional analysis. Our results revealed the conservation of Dnmts during evolution and indicated a potential role of Dnmts in epigenetic regulation of response to salinity stress.

Effects of salinity stress on methylation of the liver genome and complement gene in large yellow croaker (Larimichthys crocea).

Salinity is an important environmental factor that affects the yield and quality of large yellow croaker (Larimichthys crocea) during aquaculture. Here, whole-genome bisulfite sequencing (WGBS), RNA-seq, bisulfite sequencing PCR (BSP), quantitative real-time PCR (qPCR), and dual luciferase reporter gene detection technologies were used to analyze the DNA methylation characteristics and patterns of the liver genome, the expression and methylation levels of important immune genes in large yellow croaker in response to salinity stress. The results of WGBS showed that the cytosine methylation of CG type was dominant, CpGIsland and repeat regions were important regions where DNA methylation occurred, and the DNA methylation in upstream 2k (2000bp upstream of the promoter) and repeat regions had different changes in the liver tissue of large yellow croaker in the response to the 12‰, 24‰, 36‰ salinity stress of 4 w (weeks). In the combined analysis of WGBS and transcriptome, the complement and coagulation cascade pathways were significantly enriched, in which the complement-related genes C7, C3, C5, C4, C1R, MASP1, and CD59 were mainly changed in response to salinity stress. In the studied area of MASP1 gene promoter, the methylation levels of many CpG sites as well as total cytosine were strongly negatively correlated with mRNA expression level. Methylation function analysis of MASP1 promoter further proved that DNA methylation could inhibit the activity of MASP1 promoter, indicating that salinity may affect the expressions of complement-related genes by DNA methylation of gene promoter region.

MiR-2014-5p and miR-1231-5p regulate muscle growth of Larimichthys crocea by targeting MSTN gene.

MicroRNAs (miRNAs) play an important role in regulating gene expression, and myostatin (MSTN) has been widely recognized as a key gene for muscle growth and development. Through high-throughput sequencing to study the effects of starvation on miRNA transcriptomes in Larimichthys crocea muscle tissue, we found that the expression of miR-2014, miR-1231 and miR-1470 were significantly different between fasting and normal feeding Larimichthys crocea. Bioinformatics analysis predicted that miR-2014, miR-1231 and miR-1470 target MSTN mRNA 3'UTR. To verify the accuracy of predictions, we constructed double luciferase plasmids containing MSTN 3'UTR and confirmed that miR-2014-5p and miR-1231-5p can inhibit MSTN expression by targeting MSTN 3'UTR using double luciferase experiments, while miR-1470 is not involved in regulation. Subsequent site-directed mutation experiments reflected the specificity of the target sequence. In addition, quantitative PCR experiments revealed that miR-2014-5p and miR-1231-5p may participate in the regulation of MSTN expression in fasting and refeeding period, respectively. These results implied that miRNA may take part in muscle growth regulation during starvation. It provides some insights into the molecular regulation mechanism of MSTN in response to starvation stress in fish.

Gill transcriptomes reveal expression changes of genes related with immune and ion transport under salinity stress in silvery pomfret (Pampus argenteus).

Salinity is a major ecological factor in the marine environment, and extremely important for the survival, development, and growth of fish. In this study, gill transcriptomes were examined by high-throughput sequencing at three different salinities (12 ppt as low salinity, 22 ppt as control salinity, and 32 ppt as high salinity) in an importantly economical fish silvery pomfret. A total of 187 genes were differentially expressed, including 111 up-regulated and 76 down-regulated transcripts in low-salinity treatment group and 107 genes differentially expressed, including 74 up-regulated and 33 down-regulated transcripts in high-salinity treatment group compared with the control group, respectively. Some pathways including NOD-like receptor signaling pathway, cytokine-cytokine receptor interaction, Toll-like receptor pathway, cardiac muscle contraction, and vascular smooth muscle contraction were significantly enriched. qPCR analysis further confirmed that mRNA expression levels of immune (HSP90A, IL-1ß, TNFα, TLR2, IP-10, MIG, CCL19, and IL-11) and ion transport-related genes (WNK2, NPY2R, CFTR, and SLC4A2) significantly changed under salinity stress. Low salinity stress caused more intensive expression changes of immune-related genes than high salinity. These results imply that salinity stress may affect immune function in addition to regulating osmotic pressure in silvery pomfret.

Dynamic alterations in methylation of global DNA and growth-related genes in large yellow croaker (Larimichthys crocea) in response to starvation stress.

DNA methylation is susceptible to various environmental factors such as salinity, temperature and nutritional conditions, and can affect gene function, organ metabolism, body growth and development. In order to explore the effect of starvation on growth-related genes in large yellow croaker (Larimichthys crocea), we studied methylation of the global DNA and growth-related genes (MSTN1,MSTN2,IGF1,IGF2) and the corresponding mRNA expressions, using ELISA-based technique, bisulfite sequencing PCR (BSP) technique and quantitative Real-time PCR (qRT-PCR) respectively. The results showed that the global DNA methylation levels were significantly different (p <0.05) between the experimental group and the control group at starvation 14d, 21d in muscle and at starvation 7d, 14d, 28d, and re-feeding 7d in liver. The CpG islands of MSTN1, MSTN2, IGF1 and IGF2 were enriched in exons rather than promoters. The proximal promoter of MSTN1 and IGF1 and the exon1 of MSTN2 had almost no methylation at all treatment stages. The methylation status in MSTN1 exon 1 and IGF2 exon 2 varied from different starvation time, and started to have significant differences on starvation 7d (p <0.05) both in liver and muscle. In the liver there was a strong positive correlation between IGF2 exon 2 methylation and global DNA methylation (r = 0.7558). The mRNA expression levels of these growth-related genes were significantly different at starvation 14d (p <0.05), but did not have significant correlation with the methylation of these exons. The results implied that exon methylation of these growth-related genes might affect post-transcriptional process.