Integrated analysis of methylome, transcriptome and miRNAome of three pig breeds

The Whole Genome DNA Methylation Signatures of Hindlimb Muscles in Chinese Alligators during Hibernation and Active Periods

Many ectotherms hibernate to increase their chances of survival during harsh winter conditions. The role of DNA methylation in regulating gene expression related to hibernation in ectotherms remains unclear. Here, we employed whole-genome bisulfite sequencing (WGBS) technology to construct a comprehensive genome-wide DNA methylation landscape of the hindlimb muscles in the Chinese alligator during hibernation and active periods. The results indicated that methylation modifications were most abundant at CG sites, identifying 9447 differentially methylated regions (DMRs) and 2329 differentially methylated genes (DMGs). KEGG pathway enrichment analysis of the DMGs revealed significant enrichment in major pathways such as the neurotrophin signaling pathway, the MAPK signaling pathway, the GnRH signaling pathway, the biosynthesis of amino acids, and the regulation of the actin cytoskeleton, which are closely related to lipid metabolism, energy metabolism, and amino acid metabolism. Among these, 412 differentially methylated genes were located in promoter regions, including genes related to energy metabolism such as ATP5F1C, ATP5MD, PDK3, ANGPTL1, and ANGPTL2, and genes related to ubiquitin-proteasome degradation such as FBXO28, FBXO43, KLHL40, and PSMD5. These findings suggest that methylation in promoter regions may play a significant role in regulating the adaptive hibernation mechanisms in the Chinese alligator. This study contributes to a further understanding of the epigenetic mechanisms behind the hibernation of the Chinese alligator.

DNA methylome and transcriptome identified Key genes and pathways involved in Speckled Eggshell formation in aged laying hens

BACKGROUND

The quality of poultry eggshells is closely related to the profitability of egg production. Eggshell speckles reflect an important quality trait that influences egg appearance and customer preference. However, the mechanism of speckle formation remains poorly understood. In this study, we systematically compared serum immune and antioxidant indices of hens laying speckled and normal eggs. Transcriptome and methylome analyses were used to elucidate the mechanism of eggshell speckle formation.

RESULTS

The results showed that seven differentially expressed genes (DEGs) were identified between the normal and speckle groups. Gene set enrichment analysis (GSEA) revealed that the expressed genes were mainly enriched in the calcium signaling pathway, focal adhesion, and MAPK signaling pathway. Additionally, 282 differentially methylated genes (DMGs) were detected, of which 15 genes were associated with aging, including ARNTL, CAV1, and GCLC. Pathway analysis showed that the DMGs were associated with T cell-mediated immunity, response to oxidative stress, and cellular response to DNA damage stimulus. Integrative analysis of transcriptome and DNA methylation data identified BFSP2 as the only overlapping gene, which was expressed at low levels and hypomethylated in the speckle group.

CONCLUSIONS

Overall, these results indicate that aging- and immune-related genes and pathways play a crucial role in the formation of speckled eggshells, providing useful information for improving eggshell quality.

Early fish domestication affects methylation of key genes involved in the rapid onset of the farmed phenotype

Animal domestication is a process of environmental modulation and artificial selection leading to permanent phenotypic modifications. Recent studies showed that phenotypic changes occur very early in domestication, i.e., within the first generation in captivity, which raises the hypothesis that epigenetic mechanisms may play a critical role on the early onset of the domestic phenotype. In this context, we applied reduced representation bisulphite sequencing to compare methylation profiles between wild Nile tilapia females and their offspring reared under farmed conditions. Approximately 700 differentially methylated CpG sites were found, many of them associated not only with genes involved in muscle growth, immunity, autophagy and diet response but also related to epigenetic mechanisms, such as RNA methylation and histone modifications. This bottom-up approach showed that the phenotypic traits often related to domestic animals (e.g., higher growth rate and different immune status) may be regulated epigenetically and prior to artificial selection on gene sequences. Moreover, it revealed the importance of diet in this process, as reflected by differential methylation patterns in genes critical to fat metabolism. Finally, our study highlighted that the TGF-β1 signalling pathway may regulate and be regulated by several differentially methylated CpG-associated genes. This could be an important and multifunctional component in promoting adaptation of fish to a domestic environment while modulating growth and immunity-related traits.

Genome-Wide DNA Methylation Profile in Jejunum Reveals the Potential Genes Associated With Paratuberculosis in Dairy Cattle

Paratuberculosis in cattle causes substantial economic losses to the dairy industry. Exploring functional genes and corresponding regulatory pathways related to resistance or susceptibility to paratuberculosis is essential to the breeding of disease resistance in cattle. Co-analysis of genome-wide DNA methylation and transcriptome profiles is a critically important approach to understand potential regulatory mechanism underlying the development of diseases. In this study, we characterized the profiles of DNA methylation of jejunum from nine Holstein cows in clinical, subclinical, and healthy groups using whole-genome bisulfite sequencing (WGBS). The average methylation level in functional regions was 29.95% in the promoter, 29.65% in the 5’ untranslated region (UTR), 68.24% in exons, 71.55% in introns, and 72.81% in the 3’ UTR. A total of 3,911, 4,336, and 4,094 differentially methylated genes (DMGs) were detected in clinical vs. subclinical, clinical vs. healthy, and subclinical vs. healthy comparative group, respectively. Gene ontology (GO) and analysis based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) showed that these DMGs were significantly enriched in specific biological processes related to immune response, such as Th1 and Th2 cell differentiation, wnt, TNF, MAPK, ECM-receptor interaction, cellular senescence, calcium, and chemokine signaling pathways (q value <0.05). The integration of information about DMGs, differentially expressed genes (DEGs), and biological functions suggested nine genes CALCRL, TNC, GATA4, CD44, TGM3, CXCL9, CXCL10, PPARG, and NFATC1 as promising candidates related to resistance/susceptibility to Mycobacterium avium subspecies paratuberculosis (MAP). This study reports on the high-resolution DNA methylation landscapes of the jejunum methylome across three conditions (clinical, subclinical, and healthy) in dairy cows. Our investigations integrated different sources of information about DMGs, DEGs, and pathways, enabling us to find nine functional genes that might have potential application in resisting paratuberculosis in dairy cattle.

Improvement, identification, and target prediction for miRNAs in the porcine genome by using massive, public high-throughput sequencing data

Despite the broad variety of available microRNA (miRNA) research tools and methods, their application to the identification, annotation, and target prediction of miRNAs in nonmodel organisms is still limited. In this study, we collected nearly all public sRNA-seq data to improve the annotation for known miRNAs and identify novel miRNAs that have not been annotated in pigs (Sus scrofa). We newly annotated 210 mature sequences in known miRNAs and found that 43 of the known miRNA precursors were problematic due to redundant/missing annotations or incorrect sequences. We also predicted 811 novel miRNAs with high confidence, which was twice the current number of known miRNAs for pigs in miRBase. In addition, we proposed a correlation-based strategy to predict target genes for miRNAs by using a large amount of sRNA-seq and RNA-seq data. We found that the correlation-based strategy provided additional evidence of expression compared with traditional target prediction methods. The correlation-based strategy also identified the regulatory pairs that were controlled by nonbinding sites with a particular pattern, which provided abundant complementarity for studying the mechanism of miRNAs that regulate gene expression. In summary, our study improved the annotation of known miRNAs, identified a large number of novel miRNAs, and predicted target genes for all pig miRNAs by using massive public data. This large data-based strategy is also applicable for other nonmodel organisms with incomplete annotation information.

Integrated analysis of the methylome and transcriptome of chickens with fatty liver hemorrhagic syndrome

Background

DNA methylation, a biochemical modification of cytosine, has an important role in lipid metabolism. Fatty liver hemorrhagic syndrome (FLHS) is a serious disease and is tightly linked to lipid homeostasis. Herein, we compared the methylome and transcriptome of chickens with and without FLHS.

Results

We found genome-wide dysregulated DNA methylation pattern in which regions up- and down-stream of gene body were hypo-methylated in chickens with FLHS. A total of 4155 differentially methylated genes and 1389 differentially expressed genes were identified. Genes were focused when a negative relationship between mRNA expression and DNA methylation in promoter and gene body were detected. Based on pathway enrichment analysis, we found expression of genes related to lipogenesis and oxygenolysis (e.g., PPAR signaling pathway, fatty acid biosynthesis, and fatty acid elongation) to be up-regulated with associated down-regulated DNA methylation. In contrast, genes related to cellular junction and communication pathways (e.g., vascular smooth muscle contraction, phosphatidylinositol signaling system, and gap junction) were inhibited and with associated up-regulation of DNA methylation.

Conclusions

In the current study, we provide a genome-wide scale landscape of DNA methylation and gene expression. The hepatic hypo-methylation feature has been identified with FLHS chickens. By integrated analysis, the results strongly suggest that increased lipid accumulation and hepatocyte rupture are central pathways that are regulated by DNA methylation in chickens with FLHS.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07305-3.

A gene prioritization method based on a swine multi-omics knowledgebase and a deep learning model

The analyses of multi-omics data have revealed candidate genes for objective traits. However, they are integrated poorly, especially in non-model organisms, and they pose a great challenge for prioritizing candidate genes for follow-up experimental verification. Here, we present a general convolutional neural network model that integrates multi-omics information to prioritize the candidate genes of objective traits. By applying this model to Sus scrofa, which is a non-model organism, but one of the most important livestock animals, the model precision was 72.9%, recall 73.5%, and F1-Measure 73.4%, demonstrating a good prediction performance compared with previous studies in Arabidopsis thaliana and Oryza sativa. Additionally, to facilitate the use of the model, we present ISwine (http://iswine.iomics.pro/), which is an online comprehensive knowledgebase in which we incorporated almost all the published swine multi-omics data. Overall, the results suggest that the deep learning strategy will greatly facilitate analyses of multi-omics integration in the future.

Yuhua Fu et al. develop a CNN model that integrates multi-omics information to prioritize candidate genes of objective traits. Their model performs well when applied to important livestock non-model animals like Sus scrofa. Finally, the authors present ISwine, an online comprehensive knowledgebase which includes all published swine omics data to facilitate the integration of heterogeneous data.

An integration analysis based on genomic, transcriptomic and QTX information reveals credible candidate genes for fat-related traits in pigs

Meat quality improvement is of great interest to researchers in pig breeding and many researchers have identified abundant associated quantitative trait loci, genes and polymorphisms (QTXs) for fat-related traits. However, it is challenging to determine credible candidate genes from a mass of associations. The efficiency of identification of credible candidate genes in these QTXs is restricted by limited integration analyses of data from multiple omics. In this study, we constructed a 'candidate gene map' of fat-related traits in pigs based on published literature and the latest genome. In total, 6,861 QTXs, which covered 9,323 genes on the pig genome, were used. Combining the QTX hotspots and pathway analysis, we identified 180 candidate genes that may regulate the fat-related traits, and choose PNPLA2, PPARG, SREBF1, ACACA, PPARD and PPARA as credible candidate genes. In addition, we discussed the importance of incorporating transcriptome data and genomic data in causal gene identification, and the multi-omics information can effectively improve the credibility of identified candidate genes.

A Transcriptome Analysis Reveals that Hepatic Glycolysis and Lipid Synthesis Are Negatively Associated with Feed Efficiency in DLY Pigs

Feed efficiency (FE) is an important trait in the porcine industry. Therefore, understanding the molecular mechanisms of FE is vital for the improvement of this trait. In this study, 6 extreme high-FE and 6 low-FE pigs were selected from 225 Duroc × (Landrace × Yorkshire) (DLY) pigs for transcriptomic analysis. RNA-seq analysis was performed to determine differentially expressed genes (DEGs) in the liver tissues of the 12 individuals, and 507 DEGs were identified between high-FE pigs (HE- group) and low-FE pigs (LE- group). A gene ontology (GO) enrichment and pathway enrichment analysis were performed and revealed that glycolytic metabolism and lipid synthesis-related pathways were significantly enriched within DEGs; all of these DEGs were downregulated in the HE- group. Moreover, Weighted gene co-expression analysis (WGCNA) revealed that oxidative phosphorylation, thermogenesis, and energy metabolism-related pathways were negatively related to HE- group, which might result in lower energy consumption in higher efficiency pigs. These results implied that the higher FE in the HE- group may be attributed to a lower glycolytic, energy consumption and lipid synthesizing potential in the liver. Furthermore, our findings suggested that the inhibition of lipid synthesis and glucose metabolic activity in the liver may be strategies for improving the FE of DLY pigs.

Photoperiodic effect on the testicular transcriptome in broiler roosters

Information about the effects of photoperiod on the testicular transcriptome of broiler roosters is limited. The aim of the present study was to explore the effect of different photoperiodic regimes on gene expression in the testes of broiler breeder roosters. One hundred and twenty Arbor Acres broiler breeder roosters aged 20 weeks were assigned to one of three groups (n = 40) and subjected to different photoperiodic regimes: control (CTR; 12.5 L:11.5 D), short day (SD; 8 L:16 D) and long day (LD; 16 L:8 D). After 4 weeks, the testes of 10 randomly selected birds from each group were dissected, sliced and haematoxylin-eosin stained. The testicular transcriptome of roosters from the SD and LD groups was determined by RNA sequencing (RNA-Seq), and the results were confirmed using quantitative real-time PCR. The seminiferous tubule area and sperm count increased significantly with the prolongation of photoperiod (p < .01). Additionally, the RNA-Seq results indicated that 387 genes were upregulated and 1,052 genes were downregulated in the LD group compared with those in the SD group. Several crucial genes involved in rooster testicular development and reproduction were also screened, including heat shock proteins 90, extracellular regulated protein kinases 1, phosphatidylinositol 3-kinase, adenosine 5'-monophosphate -activated protein kinase, BCL-6 and Smad3. The differentially expressed genes were enriched in the mammalian targets of rapamycin (mTOR), forkhead box (FoxO), transforming growth factor beta (TGF-β) and insulin signalling pathway. In conclusion, a 16 hr photoperiod for 4 weeks increased the seminiferous tubule duct area and promoted spermatogenesis in the rooster's testicles, and the mTOR, FoxO, TGF-β and insulin signalling pathways may be involved.

Environmental epigenetics and epigenetic inheritance in domestic farm animals

Epigenetics refers to molecular factors and processes around DNA that can affect genome activity and gene expression independent of DNA sequence. Epigenetic mechanisms drive developmental processes and have also been shown to be tied to disease development. Many epigenetic studies have been done using plants, rodent, and human models, but fewer have focused on domestic livestock species. The goal of this review is to present current epigenetic findings in livestock species (cattle, pigs, sheep and poultry). Much of this research examined epigenetic effects following exposure to toxicants, nutritional changes or infectious disease in those animals directly exposed, or in the offspring they produced. A limited number of studies in domestic animals have examined epigenetic transgenerational inheritance in the absence of continued exposures. One example used a porcine model to investigate the effect that feeding males a diet supplemented with micronutrients had on liver DNA methylation and muscle mass in grand-offspring (the transgenerational F2 generation). Further research into how epigenetic mechanisms affect the health and production traits of domestic livestock and their offspring is important to elucidate.